moments_inflow Subroutine

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

Moment based velocity boundary condition from Sam Bennent PhD thesis "A Lattice Boltzmann Model for Diffusion of Binary Gas Mixtures"

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

! add small correction term to density which comes from momentum in the ! flow direction while substituing unknowns pdfs terms obtained ! from moments BC derivation to compute density. ! With this correction term density is recovered correctly at this node dens_correction = dot_product(vel, globBC%elemLvl(iLevel)%normal%val(:, iElem)) ! correction term for Multispecies ! rho_new = rho_old/(1.0 - dens_correction) rho = rho + rho0*dens_correction


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 :: rho
real(kind=rk), private :: vel(3)
real(kind=rk), private :: dens_correction
real(kind=rk), private :: vel_b(globBC%nElems(iLevel)*3)
real(kind=rk), private :: pressB(globBC%nElems(iLevel))
integer, private :: iELem
integer, private :: iDir
integer, private :: tmpDir
integer, private :: nLinks
integer, private :: iLink
integer, private :: elemPos
integer, private :: QQ
integer, private, allocatable:: missing_links(:)
real(kind=rk), private, allocatable:: rhs(:)
real(kind=rk), private, allocatable:: unKnown_fTmp(:)
real(kind=rk), private :: fTmp(layout%fStencil%QQ)
real(kind=rk), private :: moments(layout%fStencil%QQ)
real(kind=rk), private :: fTmp_2(layout%fStencil%QQ)
real(kind=rk), private :: rhoTmp
real(kind=rk), private :: uxTmp(3)
real(kind=rk), private :: inv_vel
real(kind=rk), private :: inv_press
integer, private :: bcVel_pos
integer, private :: bcPress_pos