bgk_advRel_generic_incomp Subroutine

public subroutine bgk_advRel_generic_incomp(fieldProp, inState, outState, auxField, neigh, nElems, nSolve, level, layout, params, varSys, derVarPos)

Advection relaxation routine for the BGK model with an explicit calculation of all equilibrium quantities. Slow and simple. This routine should only be used for testing purposes

This subroutine interface must match the abstract interface definition kernel in scheme/mus_scheme_type_module.f90 in order to be callable via compute function pointer.

Generic fetching step: Streaming for pull Local copy for push relaxation parameter

Pre-calculate velocitiy terms

Calculate equilibrium distribution functions fEq



type(mus_field_prop_type), intent(in) :: fieldProp(:)

Array of field properties (fluid or species)

real(kind=rk), intent(in) :: inState(nElems*varSys%nScalars)

input pdf vector

real(kind=rk), intent(out) :: outState(nElems*varSys%nScalars)

output pdf vector

real(kind=rk), intent(inout) :: auxField(nElems*varSys%nAuxScalars)

Auxiliary field computed from pre-collision state Is updated with correct velocity field for multicomponent models

integer, intent(in) :: neigh(nElems*layout%fStencil%QQ)

connectivity vector

integer, intent(in) :: nElems

number of elements in state Array

integer, intent(in) :: nSolve

number of elements solved in kernel

integer, intent(in) :: level

current level

type(mus_scheme_layout_type), intent(in) :: layout

current layout

type(mus_param_type), intent(in) :: params

global parameters

type(tem_varSys_type), intent(in) :: varSys

variable system definition

type(mus_derVarPos_type), intent(in) :: derVarPos(:)

position of derived quantities in varsys for all fields



integer, private :: iElem
integer, private :: iDir
integer, private :: QQ
integer, private :: nScalars
real(kind=rk), private :: pdfTmp(layout%fStencil%QQ)
real(kind=rk), private :: ux(3)
real(kind=rk), private :: rho
real(kind=rk), private :: usq
real(kind=rk), private :: ucx
real(kind=rk), private :: eqState(layout%fStencil%QQ)
real(kind=rk), private :: omega
integer, private :: dens_pos
integer, private :: vel_pos(3)
integer, private :: elemOff