mus_turbVisc_Smagorinsky_fromPreColPDF_incomp Subroutine

  1. mus_Smagorinsky_module.f90
  2. mus_Smagorinsky_module
  3. mus_turbVisc_Smagorinsky_fromPreColPDF_incomp

public subroutine mus_turbVisc_Smagorinsky_fromPreColPDF_incomp(turbVisc, turbConfig, state, neigh, auxField, densPos, velPos, nSize, nSolve, nScalars, nAuxScalars, layout, dxL, dtL, viscKine)

Calculate eddy viscosity with smagorinsky model for incompressible model using pre-collision PDF

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(out) :: turbVisc(:)

output: turbulent viscosity
type(mus_turbulence_config_type), intent(in) :: turbConfig

Contains turbulenct coefficients
real(kind=rk), intent(in) :: state(:)

state array
integer, intent(in) :: neigh(:)

neigh array to obtain precollision pdf
real(kind=rk), intent(in) :: auxField(:)

Auxiliary field variable array
integer, intent(in) :: densPos

position of density in auxField
integer, intent(in) :: velPos(3)

position of velocity components in auxField
integer, intent(in) :: nSize

number of elements in state array
integer, intent(in) :: nSolve

Number of element to solve in this level
integer, intent(in) :: nScalars

number of scalars in state array
integer, intent(in) :: nAuxScalars

number of scalars in auxField array
type(mus_scheme_layout_type), intent(in) :: layout

scheme layout
real(kind=rk), intent(in) :: dxL

current level lattice element size
real(kind=rk), intent(in) :: dtL

current level lattice time step size
real(kind=rk), intent(in) :: viscKine(:)

Background kinematic viscosity in lattice divided by dtL

Calls

proc~~mus_turbvisc_smagorinsky_fromprecolpdf_incomp~~CallsGraph proc~mus_turbvisc_smagorinsky_fromprecolpdf_incomp mus_turbVisc_Smagorinsky_fromPreColPDF_incomp proc~getequilibriumincomp getEquilibriumIncomp proc~mus_turbvisc_smagorinsky_fromprecolpdf_incomp->proc~getequilibriumincomp proc~secondmom secondMom proc~mus_turbvisc_smagorinsky_fromprecolpdf_incomp->proc~secondmom

Contents

Variables


Variables

TypeVisibilityAttributesNameInitial
integer, private :: iElem
integer, private :: iDir
integer, private :: QQ
integer, private :: elemOff
real(kind=rk), private :: visc_coeff
real(kind=rk), private :: rho
real(kind=rk), private :: vel(3)
real(kind=rk), private :: f_preCol(layout%fStencil%QQ)

precollision PDF
real(kind=rk), private :: fEq(layout%fStencil%QQ)
real(kind=rk), private :: nEq(layout%fStencil%QQ)
real(kind=rk), private :: nEqTens(6)
real(kind=rk), private :: nEqTensMag
real(kind=rk), private :: viscKineTerm
real(kind=rk), private :: nEqTensTerm