getNEq_diffusive Function

public function getNEq_diffusive(layout, omega, Sxx) result(nEq)

Calculate the non-equilibrium part of pdf from strain rate tensor based on the diffusive scaling

According to \cite Junk:2005cr \n The non-equilibrium part of pdf is set by \n where is the viscosity, is the strain rate tensor and and is the number of dimension. \n Notice here that strain rate tensor above has to be a traceless tensor, i.e. . In current implementation, the above equation is slightly modified so that the strain rate tensor is not required to be traceless anymore. In this way, calculated by this routine can recover the input strain rate tensor no matter it is traceless or not.\n Specificly the in above equation is modified slightly, i.e. This routine has a unit test program utest/mus_fNeq_diffusive_test


type(mus_scheme_layout_type), intent(in) :: layout
real(kind=rk), intent(in) :: omega
real(kind=rk), intent(in) :: Sxx(3,3)

Strain rate tensor. It is a symmetric 3x3 matrix

Return Value real(kind=rk)(layout%fStencil%QQ)


proc~~getneq_diffusive~~CallsGraph proc~getneq_diffusive getNEq_diffusive proc~convprepost convPrePost proc~getneq_diffusive->proc~convprepost

Called by

proc~~getneq_diffusive~~CalledByGraph proc~getneq_diffusive getNEq_diffusive proc~mus_init_pdf mus_init_pdf proc~mus_init_pdf->proc~getneq_diffusive proc~mus_init_byic mus_init_byIC proc~mus_init_byic->proc~mus_init_pdf proc~mus_init_flow mus_init_flow proc~mus_init_flow->proc~mus_init_byic proc~mus_initialize mus_initialize proc~mus_initialize->proc~mus_init_flow



integer, private :: iVal
integer, private :: jVal
integer, private :: iDir
real(kind=rk), private :: strain(3,3)