equilFromMacro Function

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

derive equilibrium from macro

Arguments

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

return equilibrium


Calls

proc~~equilfrommacro~~CallsGraph proc~equilfrommacro equilFromMacro proc~equilvelfrommacro equilVelFromMacro proc~equilfrommacro->proc~equilvelfrommacro

Called by

proc~~equilfrommacro~~CalledByGraph proc~equilfrommacro equilFromMacro proc~deriveeqmsliquid_fromstate deriveEqMSLiquid_FromState proc~deriveeqmsliquid_fromstate->proc~equilfrommacro proc~deriveequilmsliquid deriveEquilMSLiquid proc~deriveequilmsliquid->proc~equilfrommacro proc~deriveequilmsliquid_frommacro deriveEquilMSLiquid_FromMacro proc~deriveequilmsliquid_frommacro->proc~equilfrommacro proc~deriveequilmsliquid_fromaux deriveEquilMSLiquid_fromAux proc~deriveequilmsliquid_fromaux->proc~equilfrommacro

Contents


Variables

TypeVisibilityAttributesNameInitial
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