equilFromMacroWTDF Function

private pure function equilFromMacroWTDF(iField, mass_dens, moleFraction, velocity, inv_thermodyn_fac, layout, nFields, phi, paramBInv, resi_coeff, theta_eq) result(fEq)

derive equilibrium from macro

Arguments

Type IntentOptional 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 Value real(kind=rk), (layout%fStencil%QQ)

return equilibrium