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
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | private | :: | iDir | ||||
| integer, | private | :: | |||||
| 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 |