mus_turbVisc_Vreman_3D Subroutine

public subroutine mus_turbVisc_Vreman_3D(turbVisc, turbConfig, gradData, auxField, velPos, nSolve, nAuxScalars, dxL, dtL)

Calculate eddy viscosity with Vreman model for 3D stencil Fortran implementation of this model: \todo add reference and formula


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

output: turbulent viscosity

type(mus_turbulence_config_type), intent(in) :: turbConfig

Contains turbulenct coefficients

type(mus_gradData_type), intent(in) :: gradData

gradient data

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

Auxiliary field variable array

integer, intent(in) :: velPos(3)

position of velocity components in auxField

integer, intent(in) :: nSolve

Number of element to solve in this level

integer, intent(in) :: nAuxScalars

number of scalars in auxField array

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

current level lattice element size

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

current level lattice time step size


proc~~mus_turbvisc_vreman_3d~~CallsGraph proc~mus_turbvisc_vreman_3d mus_turbVisc_Vreman_3D proc~getgradu getGradU proc~mus_turbvisc_vreman_3d->proc~getgradu



real(kind=rk), private :: b_11
real(kind=rk), private :: b_12
real(kind=rk), private :: b_13
real(kind=rk), private :: b_22
real(kind=rk), private :: b_23
real(kind=rk), private :: b_33
real(kind=rk), private :: visc_coeff
real(kind=rk), private :: bbeta
integer, private :: iElem
integer, private :: ndims
real(kind=rk), private :: gradU(3,3,vlen)

gradient of velocity

real(kind=rk), private :: gradU_sqr
integer, private :: nChunks
integer, private :: iChunks
integer, private :: nChunkElems
integer, private :: low_bound
integer, private :: elempos