Calculate eddy viscosity with WALE (Wall-Adapting Local Eddy-viscosity) model \todo add reference and formula
gradient of velocity
| 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 |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | private | :: | iElem | ||||
| real(kind=rk), | private | :: | gradU(2,2,vlen) | ||||
| real(kind=rk), | private | :: | gradU_sqr(2,2,vlen) | ||||
| real(kind=rk), | private | :: | SR(3) | ||||
| real(kind=rk), | private | :: | Sd(3) | ||||
| real(kind=rk), | private | :: | onehalf_trSd | ||||
| real(kind=rk), | private | :: | Sd_sqr | ||||
| real(kind=rk), | private | :: | SR_sqr | ||||
| real(kind=rk), | private | :: | OP1 | ||||
| real(kind=rk), | private | :: | OP2 | ||||
| integer, | private | :: | ndims | ||||
| integer, | private | :: | nChunks | ||||
| integer, | private | :: | iChunks | ||||
| integer, | private | :: | nChunkElems | ||||
| integer, | private | :: | low_bound | ||||
| integer, | private | :: | elemPos | ||||
| real(kind=rk), | private | :: | visc_coeff |