getShearStressTensorIncomp_forPdfSubset Function

private pure function getShearStressTensorIncomp_forPdfSubset(subset, omega, layout, rho0) result(tau)

Calculate the viscous shear stress (exclude pressure) This function returns a one-dimensional array with 6 entries: tau(1:6) = [ Sxx, Syy, Szz, Sxy, Syz, Sxz ]

The formula is: \n where is the stress in the -direction on a face normal to the -axis,\n is the non-equilibirium density.\n For more information, please refer to:\n Krueger T, Varnik F, Raabe D. Shear stress in lattice Boltzmann simulations. Physical Review E. 2009;79(4):1-14.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: subset(:)
real(kind=rk), intent(in) :: omega
type(mus_scheme_layout_type), intent(in) :: layout
real(kind=rk), intent(in) :: rho0

Return Value real(kind=rk)(6)


Calls

proc~~getshearstresstensorincomp_forpdfsubset~~CallsGraph proc~getshearstresstensorincomp_forpdfsubset getShearStressTensorIncomp_forPdfSubset proc~getequilibriumincomp getEquilibriumIncomp proc~getshearstresstensorincomp_forpdfsubset->proc~getequilibriumincomp proc~convprepost convPrePost proc~getshearstresstensorincomp_forpdfsubset->proc~convprepost proc~secondmom secondMom proc~getshearstresstensorincomp_forpdfsubset->proc~secondmom

Called by

proc~~getshearstresstensorincomp_forpdfsubset~~CalledByGraph proc~getshearstresstensorincomp_forpdfsubset getShearStressTensorIncomp_forPdfSubset interface~getshearstresstensor getShearStressTensor interface~getshearstresstensor->proc~getshearstresstensorincomp_forpdfsubset proc~getwss2d getWSS2D proc~getwss2d->interface~getshearstresstensor proc~getwss getWSS proc~getwss->interface~getshearstresstensor

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: fEq(layout%fStencil%QQ)
real(kind=rk), private :: density
real(kind=rk), private :: vel(3)
real(kind=rk), private :: fNeq(layout%fStencil%QQ)
integer, private :: iDir