This module contains function to compute eddy viscosity using Vreman LES turbulence model. Vreman, A. W. (2004). An eddy-viscosity subgrid-scale model for turbulent shear flow: Algebraic theory and applications. Physics of Fluids, 16(10), 3670–3681. model. author: Kannan Masilamani
Calculate eddy viscosity with Vreman model for 3D stencil Fortran implementation of this model: http://www.vremanresearch.nl/Vreman_Subgridmodel_Fortran.txt
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
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 |
||
type(mus_Grad_type), | intent(in) | :: | Grad | Object that contains pointers to calculate gradients |
Calculate eddy viscosity with Vreman model for 2D stencil model \todo add reference and formula
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
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 |
||
type(mus_Grad_type), | intent(in) | :: | Grad | Object that contains pointers to calculate gradients |