mus_turbVisc_Smagorinsky_fromPreColPDF_incomp_2D Subroutine

public subroutine mus_turbVisc_Smagorinsky_fromPreColPDF_incomp_2D(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

Type IntentOptional Attributes Name
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_2d~~CallsGraph proc~mus_turbvisc_smagorinsky_fromprecolpdf_incomp_2d mus_turbVisc_Smagorinsky_fromPreColPDF_incomp_2D proc~getequilibriumincomp getEquilibriumIncomp proc~mus_turbvisc_smagorinsky_fromprecolpdf_incomp_2d->proc~getequilibriumincomp proc~secondmom_2d secondMom_2D proc~mus_turbvisc_smagorinsky_fromprecolpdf_incomp_2d->proc~secondmom_2d cxdirrk cxdirrk proc~getequilibriumincomp->cxdirrk