mrt_advRel_d3q19f3_MSLiquid_WTDF Subroutine

  1. mus_compute_MSLiquid_module.f90
  2. mus_MSLiquid_module
  3. mrt_advRel_d3q19f3_MSLiquid_WTDF

public subroutine mrt_advRel_d3q19f3_MSLiquid_WTDF(fieldProp, inState, outState, auxField, neigh, nElems, nSolve, level, layout, params, varSys, derVarPos)

Optimized Advection relaxation routine for the multispecies mrt model for d3q19 with 3 species with thermodynamic factor

This routine contains the implementation based on the paper "A Lattice Boltzmann Scheme for liquid mixtures - Part II: Discretization and Numerics, Jens Zudrop, Sabine Roller, Pietro Asinari. "\n MRT paper "Lattice Boltzmann liquid mixture modeling for electrodialytic engineering applications - J. Zudrop, K. Masilamani, S. Roller and P. Asinari"\n

This subroutine interface must match the abstract interface definition kernel in scheme/mus_scheme_type_module.f90 in order to be callable via compute function pointer.

Arguments

TypeIntentOptionalAttributesName
type(mus_field_prop_type), intent(in) :: fieldProp(:)

Array of field properties (fluid or species)
real(kind=rk), intent(in) :: inState(nElems*varSys%nScalars)

input pdf vector
real(kind=rk), intent(out) :: outState(nElems*varSys%nScalars)

output pdf vector
real(kind=rk), intent(inout) :: auxField(nElems*varSys%nAuxScalars)

Auxiliary field computed from pre-collision state Is updated with correct velocity field for multicomponent models
integer, intent(in) :: neigh(nElems*layout%fStencil%QQ)

connectivity vector
integer, intent(in) :: nElems

number of elements in state Array
integer, intent(in) :: nSolve

number of elements solved in kernel
integer, intent(in) :: level

current level
type(mus_scheme_layout_type), intent(in) :: layout

current layout
type(mus_param_type), intent(in) :: params

global parameters
type(tem_varSys_type), intent(in) :: varSys

variable system definition
type(mus_derVarPos_type), intent(in) :: derVarPos(:)

position of derived quantities in varsys for all fields

Calls

proc~~mrt_advrel_d3q19f3_msliquid_wtdf~~CallsGraph proc~mrt_advrel_d3q19f3_msliquid_wtdf mrt_advRel_d3q19f3_MSLiquid_WTDF interface~mus_calc_thermfactor mus_calc_thermFactor proc~mrt_advrel_d3q19f3_msliquid_wtdf->interface~mus_calc_thermfactor interface~mus_calc_ms_diffmatrix mus_calc_MS_DiffMatrix proc~mrt_advrel_d3q19f3_msliquid_wtdf->interface~mus_calc_ms_diffmatrix proc~invert_matrix invert_matrix proc~mrt_advrel_d3q19f3_msliquid_wtdf->proc~invert_matrix proc~mus_calc_thermfactor_single mus_calc_thermFactor_single interface~mus_calc_thermfactor->proc~mus_calc_thermfactor_single proc~mus_calc_ms_diffmatrix_single mus_calc_MS_DiffMatrix_single interface~mus_calc_ms_diffmatrix->proc~mus_calc_ms_diffmatrix_single interface~calc_ms_diff_matrix_from_moledens calc_ms_diff_matrix_from_moledens proc~mus_calc_ms_diffmatrix_single->interface~calc_ms_diff_matrix_from_moledens interface~calc_therm_factor_loc calc_therm_factor_loc proc~mus_calc_thermfactor_single->interface~calc_therm_factor_loc

Contents

Variables


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: pdfTmp(QQ,3)
integer, private :: iElem
integer, private :: iField
integer, private :: nScalars
integer, private :: iField_2
integer, private :: iField_3
type(mus_varSys_data_type), private, pointer:: fPtr
type(mus_scheme_type), private, pointer:: scheme
real(kind=rk), private, dimension(3):: rsigma
real(kind=rk), private, dimension(3):: uxsigma
real(kind=rk), private, dimension(3):: uysigma
real(kind=rk), private, dimension(3):: uzsigma
real(kind=rk), private, dimension(3):: qxsigma
real(kind=rk), private, dimension(3):: qysigma
real(kind=rk), private, dimension(3):: qzsigma
real(kind=rk), private, dimension(3):: gqxsigma
real(kind=rk), private, dimension(3):: gqysigma
real(kind=rk), private, dimension(3):: gqzsigma
real(kind=rk), private, dimension(3):: uxstar
real(kind=rk), private, dimension(3):: uystar
real(kind=rk), private, dimension(3):: uzstar
real(kind=rk), private, dimension(3):: molWeight
real(kind=rk), private, dimension(3):: molWeight_inv
real(kind=rk), private, dimension(3):: phi
real(kind=rk), private, dimension(3):: num_dens
real(kind=rk), private, dimension(3):: moleFrac
real(kind=rk), private :: usqr
real(kind=rk), private :: totNum_dens_inv
real(kind=rk), private :: totMass_dens
real(kind=rk), private :: velAvg_x
real(kind=rk), private :: velAvg_y
real(kind=rk), private :: velAvg_z
real(kind=rk), private :: velQuadTerm_x
real(kind=rk), private :: velQuadTerm_y
real(kind=rk), private :: velQuadTerm_z
real(kind=rk), private :: omega
real(kind=rk), private :: omegadiv2
real(kind=rk), private :: omega_fac
real(kind=rk), private :: omega_o
real(kind=rk), private :: theta_eq
real(kind=rk), private :: theta_eq_spc
real(kind=rk), private :: paramB_inv
real(kind=rk), private :: B_12
real(kind=rk), private :: B_13
real(kind=rk), private :: B_23
real(kind=rk), private :: Bratio_12
real(kind=rk), private :: Bratio_13
real(kind=rk), private :: Bratio_23
real(kind=rk), private :: mbbEq_12
real(kind=rk), private :: mbbEq_13
real(kind=rk), private :: mbbEq_21
real(kind=rk), private :: mbbEq_23
real(kind=rk), private :: mbbEq_31
real(kind=rk), private :: mbbEq_32
real(kind=rk), private :: chi12
real(kind=rk), private :: chi13
real(kind=rk), private :: chi21
real(kind=rk), private :: chi23
real(kind=rk), private :: chi31
real(kind=rk), private :: chi32
real(kind=rk), private :: rho_d
real(kind=rk), private :: rho_dm
real(kind=rk), private :: rho_o
real(kind=rk), private :: rho_om
real(kind=rk), private :: sum1_1
real(kind=rk), private :: sum2_1
real(kind=rk), private :: sum3_1
real(kind=rk), private :: sum4_1
real(kind=rk), private :: sum5_1
real(kind=rk), private :: sum6_1
real(kind=rk), private :: sum7_1
real(kind=rk), private :: sum8_1
real(kind=rk), private :: sum9_1
real(kind=rk), private :: sum1_2
real(kind=rk), private :: sum2_2
real(kind=rk), private :: sum3_2
real(kind=rk), private :: sum4_2
real(kind=rk), private :: sum5_2
real(kind=rk), private :: sum6_2
real(kind=rk), private :: sum7_2
real(kind=rk), private :: sum8_2
real(kind=rk), private :: sum9_2
real(kind=rk), private :: omegaMoments(QQ,QQ,3)
real(kind=rk), private :: fneq(QQ)
real(kind=rk), private :: fneq_om(QQ)
integer, private :: dens_pos(3)
integer, private :: mom_pos(3,3)
integer, private :: elemOff
real(kind=rk), private, dimension(3, 3):: matA
real(kind=rk), private, dimension(3, 3):: invA
real(kind=rk), private, dimension(3, 3):: resi_coeff
real(kind=rk), private, dimension(3, 3):: diff_coeff
real(kind=rk), private, dimension(3, 3):: thermodynamic_fac
real(kind=rk), private, dimension(3, 3):: inv_thermodyn_fac
real(kind=rk), private :: temp
real(kind=rk), private :: press
real(kind=rk), private :: phy_moleDens_fac