This routine checks whether omega is within the stability limit. If not it will terminate the simulation with error message. Using limits given in Tölke, J., Freudiger, S., & Krafczyk, M. (2006). An adaptive scheme using hierarchical grids for lattice Boltzmann multi-phase flow simulations. Computers & Fluids, 35(8–9), 820–830. For BGK: 2/3 < omega < 1.976 For MRT: 2/3 < omega < 1.999
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| character(len=*), | intent(in) | :: | schemeRelax |
scheme relaxation type |
||
| type(mus_relaxationParam_type), | intent(in) | :: | omLvlKine(minLevel:maxLevel) |
array of kinematic relaxation parameter on all levels |
||
| integer, | intent(in) | :: | nSolve(minLevel:maxLevel) |
Number of elements to solve in compute kernel |
||
| integer, | intent(in) | :: | minLevel |
minlevel and maxlevel |
||
| integer, | intent(in) | :: | maxLevel |
minlevel and maxlevel |
||
| type(tem_general_type), | intent(inout) | :: | general |
Contains proc, simControl, solveHead |