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 |