Routines and parameter definitions for the isothermal acoustic Eq D3Q19 model
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | private, | parameter | :: | = | 19 |
Definition of the discrete velocity set |
|
| integer, | private, | parameter | :: | qN00 | = | 1 | |
| integer, | private, | parameter | :: | q0N0 | = | 2 | |
| integer, | private, | parameter | :: | q00N | = | 3 | |
| integer, | private, | parameter | :: | q100 | = | 4 | |
| integer, | private, | parameter | :: | q010 | = | 5 | |
| integer, | private, | parameter | :: | q001 | = | 6 | |
| integer, | private, | parameter | :: | q0NN | = | 7 | |
| integer, | private, | parameter | :: | q0N1 | = | 8 | |
| integer, | private, | parameter | :: | q01N | = | 9 | |
| integer, | private, | parameter | :: | q011 | = | 10 | |
| integer, | private, | parameter | :: | qN0N | = | 11 | |
| integer, | private, | parameter | :: | q10N | = | 12 | |
| integer, | private, | parameter | :: | qN01 | = | 13 | |
| integer, | private, | parameter | :: | q101 | = | 14 | |
| integer, | private, | parameter | :: | qNN0 | = | 15 | |
| integer, | private, | parameter | :: | qN10 | = | 16 | |
| integer, | private, | parameter | :: | q1N0 | = | 17 | |
| integer, | private, | parameter | :: | q110 | = | 18 | |
| integer, | private, | parameter | :: | q000 | = | 19 | |
| real(kind=rk), | private, | parameter | :: | f1 | = | 2.0_rk/5.0_rk | |
| real(kind=rk), | private, | parameter | :: | f2 | = | 1.0_rk/30.0_rk | |
| real(kind=rk), | private, | parameter | :: | f8 | = | 1.0_rk/30.0_rk |
Advection relaxation routine for the D3Q19 model with BGK for the isothermal acoustic equation.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| 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 |