Advection relaxation routine based on the L-Model for the D3Q19 Lattice Boltzmann model for the generic advection and anisotropic-dispersion equation.
Irina Ginzburg (2005), "Equilibrium-type and link-type lattice Boltzmann models for generic advection and anisotropic-dispersion equation", Advances in Water Resources, Volume 28, Issue 11
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.
compute the link-wise omega the directions of omegas are 1 !< west x- 2 !< south y- 3 !< bottom z- 4 !< east x+ 5 !< north y+ 6 !< top z+ 7 !< z-,y- 8 !< z+,y- 9 !< z-,y+ 10 !< z+,y+ 11 !< x-,z- 12 !< x+,z- 13 !< x-,z+ 14 !< x+,z+ 15 !< y-,x- 16 !< y+,x- 17 !< y-,x+ 18 !< y+,x+ 19 !< rest density is last
| 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 |