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 |
||

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 | molecular weight ratio of iField |
||

real(kind=rk), | intent(in) | :: | paramBInv | free parameter B |
||

real(kind=rk), | intent(in) | :: | resi_coeff(nFields) | resistivity coefficients |
||

real(kind=rk), | intent(in) | :: | theta_eq | parameter to tune mixture velocity in equilibrium quadratic term |

return equilibrium

Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|

integer, | private | :: | iDir | ||||

integer, | private | :: | QQ | ||||

real(kind=rk), | private | :: | totMass_densInv | ||||

real(kind=rk), | private | :: | ucx | ||||

real(kind=rk), | private | :: | usq | ||||

real(kind=rk), | private | :: | ucxQuadTerm | ||||

real(kind=rk), | private | :: | velAvg(3) | ||||

real(kind=rk), | private | :: | velQuadTerm(3) | ||||

real(kind=rk), | private | :: | eqVel(3) | ||||

real(kind=rk), | private | :: | weight0Inv | Inverse of lattice weight ar restPosition |