pressure_nonEqExpol Subroutine

  1. mus_bc_fluid_nonEqExpol_module.f90
  2. mus_bc_fluid_nonEqExpol_module
  3. pressure_nonEqExpol

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

Element-wise Dirichlet pressure non-equilibrium boundary condition for straight boundary and updates all directions. For straight wall, values are extrapolated along boundary normal instead of along the link and qVal=0.

Notation: b (fictious boundary) w (physical bolundary or surface) f (local element) ff (overnext fluidneighbor)

Usage

boundary_condition = {
  { label = 'outlet',
    kind = 'pressure_nonEqExpol',
    pressure = 'press,
    curved = true
  }
}
variable = {
  name = 'press',
  ncomponents = 1,
  vartype = 'st_fun',
  st_fun = 0.0
}

It is based on the paper: Guo, Z.; Zheng, C. & Shi B. (2002). An extrapolation method for boundary conditions in lattice Boltzmann method. Physics of Fluids 14, 2007 (2002); https://doi.org/10.1063/1.1471914

More informations concerning the nonEqExpol can be found in the corresponding subroutine.

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
integer, private :: iDir
real(kind=rk), private :: feq_b
real(kind=rk), private :: rho_b
real(kind=rk), private :: rho_w(globBC%nElems(iLevel))
real(kind=rk), private :: feq_ff
real(kind=rk), private :: rho_ff
real(kind=rk), private :: rhoInv
real(kind=rk), private :: vel_ff(3)
real(kind=rk), private :: pdf_ff(layout%fStencil%QQ)
integer, private :: bcPress_pos
integer, private :: QQ
integer, private :: iElem
integer, private :: elemPos
real(kind=rk), private :: inv_rho_phy