getShearRateTensor_diffusive_forPdfSubset Function

private pure function getShearRateTensor_diffusive_forPdfSubset(f, omega, layout) result(S)

Calculate the strain rate tensor through 2nd moment. This function returns a two-dimensional 3 x 3 symmetirc array:

 S = | Sxx, Sxy, Sxz,|
     | Syx, Syy, Syz |
     | Szx, Szy, Szz |

The formula to compute this can be found in paper of Junk 2005 after equation A.15 where is the second velocity moment, is the identity matrix and .\n Note that is not shear rate yet.

\todo: this section can become a function that be called by several places.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: f(layout%fStencil%QQ)

pdf array ( post-collision value )

real(kind=rk), intent(in) :: omega

relaxation parameter

type(mus_scheme_layout_type), intent(in) :: layout

stencil layout

Return Value real(kind=rk)(3,3)

output array: strain rate tensor


Calls

proc~~getshearratetensor_diffusive_forpdfsubset~~CallsGraph proc~getshearratetensor_diffusive_forpdfsubset getShearRateTensor_diffusive_forPdfSubset proc~get_moment get_moment proc~getshearratetensor_diffusive_forpdfsubset->proc~get_moment proc~convprepost convPrePost proc~getshearratetensor_diffusive_forpdfsubset->proc~convprepost proc~get_momentvector get_momentVector proc~get_moment->proc~get_momentvector proc~mus_imomvector mus_iMomVector proc~get_momentvector->proc~mus_imomvector

Called by

proc~~getshearratetensor_diffusive_forpdfsubset~~CalledByGraph proc~getshearratetensor_diffusive_forpdfsubset getShearRateTensor_diffusive_forPdfSubset interface~getshearratetensor_diffusive getShearRateTensor_diffusive interface~getshearratetensor_diffusive->proc~getshearratetensor_diffusive_forpdfsubset

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: vel(3)
real(kind=rk), private :: diagVal(3,3)
real(kind=rk), private :: rho
integer, private :: iVal
integer, private :: jVal
integer, private :: order(3)