Set parameters for multispecies
Relaxation parameter for each level kine_viscosity = dtc^2(1/omega - 0.5) => omega = dt / (viscosity/c^2 + dt/2) omega for each level is stored at fluid%omLvl
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(mus_field_type), | intent(inout) | :: | field(nFields) | |||
| integer, | intent(in) | :: | nFields | |||
| type(mus_mixture_type), | intent(inout) | :: | mixture | |||
| type(mus_scheme_header_type), | intent(in) | :: | header | |||
| type(mus_scheme_layout_type), | intent(in) | :: | layout | |||
| integer, | intent(in) | :: | iLevel | |||
| type(tem_time_type), | intent(in) | :: | tNow | solver general info |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| real(kind=rk), | private | :: | omega_diff | ||||
| real(kind=rk), | private | :: | omega_kine | ||||
| real(kind=rk), | private | :: | omega_bulk | ||||
| real(kind=rk), | private | :: | fac | ||||
| integer, | private | :: | iField | ||||
| integer, | private | :: | iDir | ||||
| real(kind=rk), | private, | dimension( layout%fStencil%QQ, layout%fStencil%QQ ) | :: | identity | |||
| real(kind=rk), | private, | dimension( layout%fStencil%QQ, layout%fStencil%QQ ) | :: | tmpMatrix | |||
| real(kind=rk), | private, | dimension( layout%fStencil%QQ, layout%fStencil%QQ ) | :: | tmpMatrixInv |