Interface to calculate kinematic viscosity for non-Newtonian model. Viscosity is computed from shear rate which is computed from strain rate which is computed from nonEquilibrium PDF which in turn is 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: kinematic viscosity from 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 |