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 |