bgk_HybridRecursiveRegularized_d3q19 Subroutine

  1. mus_compute_d3q19_module.f90
  2. mus_d3q19_module
  3. bgk_HybridRecursiveRegularized_d3q19

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

Projected Recursive Regularized relaxation routine for the D3Q19 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. works for both d3q19 and d3q27

TODO:tr_SR = div2_3 * (gradU(1, 1, 1) + gradU(2, 2, 1) + gradU(3, 3, 1))

Arguments

TypeIntentOptionalAttributesName
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

Calls

proc~~bgk_hybridrecursiveregularized_d3q19~~CallsGraph proc~bgk_hybridrecursiveregularized_d3q19 bgk_HybridRecursiveRegularized_d3q19 proc~second_order_moments_3d second_order_moments_3D proc~bgk_hybridrecursiveregularized_d3q19->proc~second_order_moments_3d proc~f_f_eq_regularized_4th_ord_d3q19_corr f_f_eq_regularized_4th_ord_d3q19_corr proc~bgk_hybridrecursiveregularized_d3q19->proc~f_f_eq_regularized_4th_ord_d3q19_corr proc~f_f_eq_regularized_2nd_ord_d3q19 f_f_eq_regularized_2nd_ord_d3q19 proc~f_f_eq_regularized_4th_ord_d3q19_corr->proc~f_f_eq_regularized_2nd_ord_d3q19

Contents

Variables


Variables

TypeVisibilityAttributesNameInitial
integer, private :: iElem
integer, private :: iDir
type(mus_gradData_type), private, pointer:: gradData

gradient data
type(mus_varSys_data_type), private, pointer:: fPtr
type(mus_scheme_type), private, pointer:: scheme
real(kind=rk), private :: f(QQ)

TODO:, tr_SR
real(kind=rk), private :: SR(6)

TODO:, tr_SR
real(kind=rk), private :: gradU(3,3,1)

TODO:, tr_SR
real(kind=rk), private :: SOM(6)

TODO:, tr_SR
real(kind=rk), private :: SOM_neq(6)

TODO:, tr_SR
real(kind=rk), private :: rho
real(kind=rk), private :: u_x
real(kind=rk), private :: u_y
real(kind=rk), private :: u_z
real(kind=rk), private :: a12xx
real(kind=rk), private :: a12xy
real(kind=rk), private :: a12yy
real(kind=rk), private :: a12zz
real(kind=rk), private :: a12xz
real(kind=rk), private :: a12yz
real(kind=rk), private :: omega
real(kind=rk), private :: taup
real(kind=rk), private :: cmpl_o
real(kind=rk), private :: feq(QQ)
real(kind=rk), private :: f1(QQ)
real(kind=rk), private :: sigma
integer, private :: denspos
integer, private :: velpos(3)
integer, private :: elemOff
integer, private :: nScalars