getShearStressTensor_forElemFromState Function

private pure function getShearStressTensor_forElemFromState(state, neigh, elem, omega, layout, iField, varPos, nScalars) 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) :: state(:)
integer, intent(in) :: neigh(:)
integer, intent(in) :: elem
real(kind=rk), intent(in) :: omega
type(mus_scheme_layout_type), intent(in) :: layout
integer, intent(in) :: iField
integer, intent(in) :: varPos(:)
integer, intent(in) :: nScalars

Return Value real(kind=rk)(6)


Calls

proc~~getshearstresstensor_forelemfromstate~~CallsGraph proc~getshearstresstensor_forelemfromstate getShearStressTensor_forElemFromState proc~convprepost convPrePost proc~getshearstresstensor_forelemfromstate->proc~convprepost interface~getequilibrium getEquilibrium proc~getshearstresstensor_forelemfromstate->interface~getequilibrium proc~secondmom secondMom proc~getshearstresstensor_forelemfromstate->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_forelemfromstate~~CalledByGraph proc~getshearstresstensor_forelemfromstate getShearStressTensor_forElemFromState interface~getshearstresstensor getShearStressTensor interface~getshearstresstensor->proc~getshearstresstensor_forelemfromstate 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)
integer, private :: iDir
integer, private :: nElems
integer, private :: QQ