Calculate kinematic viscosity from nonNewtonian Carreau-Yasuda model for incompressible model. $\mu = \mu_\inf + (\mu_0-\mu_\inf)(1+(\lambdashearRate)a)^((n-1)/a)$
Shear rate is computed from strain rate which is computed from nonEquilibrium PDF which in turn computed from pre-collision PDF
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(mus_nNwtn_type), | intent(in) | :: | nNwtn | contains non-Newtonian model parameters loaded from config file |
||
| real(kind=rk), | intent(inout) | :: | viscKine(:) | non-Netonian model |
||
| real(kind=rk), | intent(in) | :: | omega(:) | Kinematic viscosity omega from last timestep |
||
| real(kind=rk), | intent(in) | :: | state(:) | state array |
||
| integer, | intent(in) | :: | neigh(:) | neigh array to obtain precollision pdf |
||
| real(kind=rk), | intent(in) | :: | auxField(:) | Auxiliary field variable array |
||
| integer, | intent(in) | :: | densPos | position of density in auxField |
||
| integer, | intent(in) | :: | velPos(3) | position of velocity components in auxField |
||
| integer, | intent(in) | :: | nSize | number of elements in state array |
||
| integer, | intent(in) | :: | nSolve | Number of element to solve in this level |
||
| integer, | intent(in) | :: | nScalars | number of scalars in state array |
||
| integer, | intent(in) | :: | nAuxScalars | number of scalars in auxField array |
||
| type(mus_scheme_layout_type), | intent(in) | :: | layout | scheme layout |
||
| type(mus_convertFac_type), | intent(in) | :: | convFac | conversion factor to convert lattice to physical units |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | private | :: | iElem | ||||
| integer, | private | :: | iDir | ||||
| integer, | private | :: | |||||
| integer, | private | :: | elemOff | ||||
| real(kind=rk), | private | :: | rho | ||||
| real(kind=rk), | private | :: | vel(3) | ||||
| real(kind=rk), | private | :: | f_preCol(layout%fStencil%QQ) | precollision PDF |
|||
| real(kind=rk), | private | :: | fEq(layout%fStencil%QQ) | ||||
| real(kind=rk), | private | :: | nEq(layout%fStencil%QQ) | ||||
| real(kind=rk), | private | :: | nEqTens(6) | ||||
| real(kind=rk), | private | :: | nEqTensMag | ||||
| real(kind=rk), | private | :: | shearRate | ||||
| real(kind=rk), | private | :: | strainRate | ||||
| real(kind=rk), | private | :: | v0_vInf | ||||
| real(kind=rk), | private | :: | coeffSR | ||||
| real(kind=rk), | private | :: | viscDynaPhy | ||||
| real(kind=rk), | private | :: | viscTerm |