spc_inlet_eq Subroutine

  1. mus_bc_species_module.f90
  2. mus_bc_species_module
  3. spc_inlet_eq

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

Inlet species velocity equilibrium boundary with specified mixture averaged mass velocity and its molefraction mixture kinematic pressure is extrapolated here. Density and velocity of all fields are used to compute equilibrium Usage

boundary_condition = {
 { label = 'outlet',
    kind = 'spc_inlet_eq',
    velocity = 'inlet_vel',
    mole_fraction = 'inlet_mole'
 }
}
variable = {
 {
   label = 'inlet_vel',
   ncomponents = 3,
   vartype = 'st_fun',
   st_fun = {0.1,0.0,0.0}
 },
 {
   label = 'inlet_mole',
   ncomponents = 1,
   vartype = 'st_fun',
   st_fun = 0.1
 }
}

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.

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 :: fTmp(layout%fStencil%QQ*globBC%nElems(iLevel)*varSys%nStateVars)
real(kind=rk), private :: fEq(layout%fStencil%QQ*globBC%nElems(iLevel))
real(kind=rk), private :: mass_dens(globBC%nElems(iLevel)*varSys%nStateVars)
real(kind=rk), private :: uxB(globBC%nElems(iLevel)*varSys%nStateVars*3)

Velocity on boundary element
real(kind=rk), private :: velocity(3,globBC%nElems(iLevel)*varSys%nStateVars)
real(kind=rk), private :: spc_vel(globBC%nElems(iLevel)*3)
real(kind=rk), private :: inv_vel
real(kind=rk), private :: moleFrac(globBC%nElems(iLevel))
integer, private :: iElem
integer, private :: iDir
integer, private :: iFieldLoc
integer, private :: nFields
integer, private :: pos
integer, private :: QQ
integer, private :: offset
integer, private :: bcVel_pos
integer, private :: bcMoleFrac_pos
integer, private :: elemPos
integer, private :: posInBuffer