getEquilibrium_forElemfromState Function

private pure function getEquilibrium_forElemfromState(state, elem, layout, varPos, nScalars, neigh) result(equil)

Calculate the equilibrium distribution function in all directions

The equilibrim distribution function is:\n \n where is the weight in each direction,\n is the macroscopic value of density,\n is the speed of sound,\n is the lattice unit velocity in each direction,\n is the macroscopic value of velocity.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: state(:)
integer, intent(in) :: elem
type(mus_scheme_layout_type), intent(in) :: layout
integer, intent(in) :: varPos(:)
integer, intent(in) :: nScalars
integer, intent(in) :: neigh(:)

Return Value real(kind=rk)(layout%fStencil%QQ)


Calls

proc~~getequilibrium_forelemfromstate~~CallsGraph proc~getequilibrium_forelemfromstate getEquilibrium_forElemfromState interface~getvelocity getVelocity proc~getequilibrium_forelemfromstate->interface~getvelocity interface~getdensity getDensity proc~getequilibrium_forelemfromstate->interface~getdensity proc~getvelocity_forpdfsubset getVelocity_forPdfSubset interface~getvelocity->proc~getvelocity_forpdfsubset proc~getvelocity_forelemfromstate_noforce getVelocity_forElemFromState_noForce interface~getvelocity->proc~getvelocity_forelemfromstate_noforce proc~getdensity_forelemfromstate getDensity_forElemFromState interface~getdensity->proc~getdensity_forelemfromstate proc~getdensity_forpdfsubset getDensity_forPdfSubset interface~getdensity->proc~getdensity_forpdfsubset

Called by

proc~~getequilibrium_forelemfromstate~~CalledByGraph proc~getequilibrium_forelemfromstate getEquilibrium_forElemfromState interface~getequilibrium getEquilibrium interface~getequilibrium->proc~getequilibrium_forelemfromstate proc~set_pdfacoustic set_pdfAcoustic proc~set_pdfacoustic->interface~getequilibrium proc~getshearstresstensor_forelemfromstate getShearStressTensor_forElemFromState proc~getshearstresstensor_forelemfromstate->interface~getequilibrium proc~getshearratetensor_acoustic_lbm getShearRateTensor_acoustic_lbm proc~getshearratetensor_acoustic_lbm->interface~getequilibrium proc~mus_turbvisc_smagorinsky_fromprecolpdf mus_turbVisc_Smagorinsky_fromPreColPDF proc~mus_turbvisc_smagorinsky_fromprecolpdf->interface~getequilibrium proc~getshearstresstensor_forpdfsubset getShearStressTensor_forPdfSubset proc~getshearstresstensor_forpdfsubset->interface~getequilibrium interface~getshearratetensor_acoustic getShearRateTensor_acoustic interface~getshearratetensor_acoustic->proc~getshearratetensor_acoustic_lbm interface~getshearstresstensor getShearStressTensor interface~getshearstresstensor->proc~getshearstresstensor_forelemfromstate interface~getshearstresstensor->proc~getshearstresstensor_forpdfsubset proc~getwss2d getWSS2D proc~getwss2d->interface~getshearstresstensor proc~getwss getWSS proc~getwss->interface~getshearstresstensor

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: rho
real(kind=rk), private :: vel(3)
real(kind=rk), private :: ucx
real(kind=rk), private :: usq
integer, private :: iDir