getShearStressTensor_forPdfSubset Function

private pure function getShearStressTensor_forPdfSubset(subset, omega, layout, varPos) 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
integer, intent(in) :: varPos(:)

Return Value real(kind=rk)(6)


Calls

proc~~getshearstresstensor_forpdfsubset~~CallsGraph proc~getshearstresstensor_forpdfsubset getShearStressTensor_forPdfSubset proc~convprepost convPrePost proc~getshearstresstensor_forpdfsubset->proc~convprepost interface~getequilibrium getEquilibrium proc~getshearstresstensor_forpdfsubset->interface~getequilibrium proc~secondmom secondMom proc~getshearstresstensor_forpdfsubset->proc~secondmom proc~getequilibrium_forpdfsubset getEquilibrium_forPdfSubset interface~getequilibrium->proc~getequilibrium_forpdfsubset proc~getequilibrium_forelemfromstate getEquilibrium_forElemfromState interface~getequilibrium->proc~getequilibrium_forelemfromstate proc~geteqbydensvel getEqByDensVel interface~getequilibrium->proc~geteqbydensvel proc~getvelocity_forpdfsubset getVelocity_forPdfSubset proc~getequilibrium_forpdfsubset->proc~getvelocity_forpdfsubset proc~getdensity_forpdfsubset getDensity_forPdfSubset proc~getequilibrium_forpdfsubset->proc~getdensity_forpdfsubset interface~getvelocity getVelocity proc~getequilibrium_forelemfromstate->interface~getvelocity interface~getdensity getDensity proc~getequilibrium_forelemfromstate->interface~getdensity interface~getvelocity->proc~getvelocity_forpdfsubset proc~getvelocity_forelemfromstate_noforce getVelocity_forElemFromState_noForce interface~getvelocity->proc~getvelocity_forelemfromstate_noforce interface~getdensity->proc~getdensity_forpdfsubset proc~getdensity_forelemfromstate getDensity_forElemFromState interface~getdensity->proc~getdensity_forelemfromstate

Called by

proc~~getshearstresstensor_forpdfsubset~~CalledByGraph proc~getshearstresstensor_forpdfsubset getShearStressTensor_forPdfSubset interface~getshearstresstensor getShearStressTensor 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 :: fEq(layout%fStencil%QQ)
real(kind=rk), private :: fNeq(layout%fStencil%QQ)