proc_calc_nNwtn_visc_fromPreColPDF Abstract Interface

abstract interface

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

private subroutine proc_calc_nNwtn_visc_fromPreColPDF(nNwtn, viscKine, omega, state, neigh, auxField, densPos, velPos, nSize, nSolve, nScalars, nAuxScalars, layout, convFac)


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