Calculate kinematic viscosity from nonNewtonian Carreau-Yasuda 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(:) |
output is physical kinematic viscosity will be overwritten by 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 |