moments_outflow Subroutine

  1. mus_bc_fluid_experimental_module.f90
  2. mus_bc_fluid_experimental_module
  3. moments_outflow

public subroutine moments_outflow(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.

do iElem = 1, globBC%nElems( iLevel )
  !ux = 0.0_rk
  do iDir = 1,QQ
    fTmp_2(iDir) = state(                                                  &

&nelems (( idir-1) neigh+ globbc%elemlvl(ilevel)%elem%val( ielem ))+( ifield-1) qq+ nscalars0 & & ) end do rhoTmp = sum(fTmp_2) call derVarPos%velFromState( state = fTmp_2, & & iField = iField, & & nElems = 1, & & varSys = varSys, & & layout = layout, & & res = uxTmp ) write(,) 'iElem ', iElem write(,) 'density ', rhoTmp, cs2rhoTmp ! write(,) 'velocity ', uxTmp ! write(,) 'velocity ', uxTmp*rhoTmp end do

Arguments

TypeIntentOptionalAttributesName
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

Contents

Variables


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: rho
real(kind=rk), private :: press
real(kind=rk), private :: vel(3)
integer, private :: iELem
integer, private :: iDir
integer, private :: QQ
integer, private :: nLinks
integer, private :: iLink
integer, private :: elemPos
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 :: uxB_1(3*globBC%nElems(iLevel))
real(kind=rk), private :: rho_1(globBC%nElems(iLevel))