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 | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|

integer, | private | :: | iElem | ||||

integer, | private | :: | ndims | ||||

real(kind=rk), | private | :: | b_11 | ||||

real(kind=rk), | private | :: | b_12 | ||||

real(kind=rk), | private | :: | b_22 | ||||

real(kind=rk), | private | :: | gradU_sqr | ||||

real(kind=rk), | private | :: | visc_coeff | ||||

real(kind=rk), | private | :: | bbeta | ||||

real(kind=rk), | private | :: | gradU(2,2,vlen) | gradient of velocity |
|||

integer, | private | :: | nChunks | ||||

integer, | private | :: | iChunks | ||||

integer, | private | :: | nChunkElems | ||||

integer, | private | :: | low_bound | ||||

integer, | private | :: | elempos |