derive equilibrium from macro
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | iField |
current field |
||
real(kind=rk), | intent(in) | :: | mass_dens(nFields) |
mass density of all species |
||
real(kind=rk), | intent(in) | :: | moleFraction(nFields) |
molefraction of all species |
||
real(kind=rk), | intent(in) | :: | velocity(3,nFields) |
velocity of all species |
||
real(kind=rk), | intent(in) | :: | inv_thermodyn_fac(nFields,nFields) |
inverse of thermodynamic factor |
||
type(mus_scheme_layout_type), | intent(in) | :: | layout |
scheme layout contains stencil definition and lattice weight |
||
integer, | intent(in) | :: | nFields |
number of species |
||
real(kind=rk), | intent(in) | :: | phi(nFields) |
molecular weight ratio of iField |
||
real(kind=rk), | intent(in) | :: | paramBInv |
free parameter B |
||
real(kind=rk), | intent(in) | :: | resi_coeff(nFields,nFields) |
resistivity coefficients |
||
real(kind=rk), | intent(in) | :: | theta_eq |
parameter to tune mixture velocity in equilibrium quadratic term |
return equilibrium