Moment based open boundary condition from Sam Bennent PhD thesis "A Lattice Boltzmann Model for Diffusion of Binary Gas Mixtures" for incompressible LBM model
Usage
boundary_condition = {
{ label = 'inlet',
kind = 'pressure_momentsbased',
pressure = 'p0'}
}
variable = {
name = 'p0',
ncomponents = 1,
vartype = 'st_fun',
st_fun = 1.0
}
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.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| 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 |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| real(kind=rk), | private | :: | pressB(globBC%nElems(iLevel)) | ||||
| real(kind=rk), | private | :: | rho | ||||
| real(kind=rk), | private | :: | vel(3) | ||||
| real(kind=rk), | private | :: | inv_press | ||||
| integer, | private | :: | iELem | ||||
| integer, | private | :: | iDir | ||||
| integer, | private | :: | nLinks | ||||
| integer, | private | :: | iLink | ||||
| integer, | private | :: | |||||
| 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 | :: | first_moments(3) | ||||
| real(kind=rk), | private | :: | second_moments(6) | ||||
| real(kind=rk), | private | :: | third_moments(6) | ||||
| real(kind=rk), | private | :: | fourth_moments(3) | ||||
| integer, | private | :: | nMom(4) | ||||
| integer, | private | :: | posInBuffer | ||||
| integer, | private | :: | bcPress_pos |