spc_moleFrac Subroutine

public subroutine spc_moleFrac(me, state, bcBuffer, globBC, levelDesc, tree, nSize, iLevel, sim_time, neigh, layout, fieldProp, varPos, nScalars, varSys, derVarPos, physics, iField, mixture)

Mole fraction boundary condition Usage

boundary_condition = {
 { label = 'outlet',
    kind = 'spc_moleFrac',
    moleFraction = 0.0

Post collision pdf of incoming link is updated with equilibrium functions which is similar to initial condition. \f$ \bar{f^c_k} = f^{eq}_k(\rho_k, u_k) + \frac{\lambda}{2} (f^{eq}_k(\rho_k, u_k) + f^{eq}_k(\rho_k, u^{eq}_k})) \f$ Here, $u_k$ - velocity of species k from original pdf computed from LSE $u^{eq}_k$ - equilibrium velocity of species k

This subroutine's interface must match the abstract interface definition boundaryRoutine in bc/mus_bc_header_module.f90 in order to be callable via [[boundary_type::fnct]] function pointer.


class(boundary_type) :: me

global boundary type

real(kind=rk), intent(inout) :: state(:)

Current state vector of iLevel

real(kind=rk), intent(in) :: bcBuffer(:)

state values of boundary elements of all fields of iLevel

type(glob_boundary_type), intent(in) :: globBC

scheme global boundary type

type(tem_levelDesc_type), intent(in) :: levelDesc

iLevel descriptor

type(treelmesh_type), intent(in) :: tree

Treelm Mesh

integer, intent(in) :: nSize

size of state array ( in terms of elements )

integer, intent(in) :: iLevel

the level On which this boundary was invoked

type(tem_time_type), intent(in) :: sim_time

global time information

integer, intent(in) :: neigh(:)

connectivity array corresponding to state vector

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

stencil layout information

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

fluid parameters and properties

integer, intent(in) :: varPos(:)

pointer to field variable in the state vector

integer, intent(in) :: nScalars

number of Scalars in the scheme var system

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

scheme variable system

type(mus_derVarPos_type), intent(in) :: derVarPos

position of derived quantities in varsys

type(mus_physics_type), intent(in) :: physics

scheme global boundary type

integer, intent(in) :: iField

current field

type(mus_mixture_type), intent(in) :: mixture

mixture info



real(kind=rk), private :: moleFrac(globBC%nElems(iLevel))
real(kind=rk), private :: rho
real(kind=rk), private :: press
integer, private :: iElem
integer, private :: iDir
integer, private :: iFieldLoc
integer, private :: nFields
integer, private :: pos
integer, private :: QQ
integer, private :: bcMoleFrac_pos
integer, private :: elemPos
integer, private :: posInBuffer
real(kind=rk), private :: fTmp_all(layout%fStencil%QQ*varSys%nStateVars)
real(kind=rk), private :: fEq
real(kind=rk), private :: fEqStar
real(kind=rk), private :: ucx
real(kind=rk), private :: ucxStar
real(kind=rk), private :: ucxQuad
real(kind=rk), private :: ucxQuadStar
real(kind=rk), private :: usq
real(kind=rk), private :: usqStar
real(kind=rk), private :: velAvg(3)
real(kind=rk), private :: velQuad(3)
real(kind=rk), private :: velQuadStar(3)
real(kind=rk), private :: eqVel(3)
real(kind=rk), private :: vel(3,varSys%nStateVars)
real(kind=rk), private :: mom(3,varSys%nStateVars)
real(kind=rk), private :: mass_dens(varSys%nStateVars)
real(kind=rk), private :: num_dens(varSys%nStateVars)
real(kind=rk), private :: totMassDens
real(kind=rk), private :: moleFrac_loc(varSys%nStateVars)
real(kind=rk), private :: molWeightInv(varSys%nStateVars)
real(kind=rk), private :: phi(varSys%nStateVars)
real(kind=rk), private :: resi_coeff(varSys%nStateVars,varSys%nStateVars)
real(kind=rk), private :: weight0_inv
real(kind=rk), private :: paramBInv
type(mus_varSys_data_type), private, pointer:: fPtr