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


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


proc~~equilfrommacrowtdf~~CallsGraph proc~equilfrommacrowtdf equilFromMacroWTDF proc~equilvelfrommacrowtdf equilVelFromMacroWTDF proc~equilfrommacrowtdf->proc~equilvelfrommacrowtdf

Called by

proc~~equilfrommacrowtdf~~CalledByGraph proc~equilfrommacrowtdf equilFromMacroWTDF proc~deriveequilwtdf_msliquid deriveEquilWTDF_MSLiquid proc~deriveequilwtdf_msliquid->proc~equilfrommacrowtdf



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