mus_turbVisc_Vreman_3D Subroutine

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

Calculate eddy viscosity with Vreman model for 3D stencil Fortran implementation of this model: http://www.vremanresearch.nl/Vreman_Subgridmodel_Fortran.txt

$$\nu_{turb}=c_v (\Delta_x)^2 \cdot \left( \sqrt{\frac{B_\beta}{\alpha_{ij}\alpha{ij}}} \right)$$ with $$c_v = 2.5*C_s^2$$ , $$C_s$$ - Smagorinsky constant,

$$B_\beta = \beta_{11}\beta_{22} - \beta^2_{12} + \beta_{11}\beta_{33} - \beta^2_{13} + \beta_{22}\beta_{33} - \beta^2_{23}$$ ,

$$\beta_{ij} = \alpha_{mi}\alpha_{mj}$$ ,

$$\alpha_{ij} = \frac{\partial \Bar{u}_j}{\partial x_i}$$ . $$\alpha_{ij}$$ - Resolved velocity gradient.

Arguments

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