8.0
general documentation
visecv.f90 File Reference

Computes the secondary viscosity contribution \(\kappa -\dfrac{2}{3} \mu\) in order to compute: More...

Functions/Subroutines

subroutine visecv (secvif, secvib)
 

Detailed Description

Computes the secondary viscosity contribution \(\kappa -\dfrac{2}{3} \mu\) in order to compute:

\[ \grad\left( (\kappa -\dfrac{2}{3} \mu) \trace( \gradt(\vect{u})) \right) \]

with:

  • \( \mu = \mu_{laminar} + \mu_{turbulent} \)
  • \( \kappa \) is the volume viscosity (generally zero)
Remarks
In LES, the tensor \(\overline{\left(\vect{u}-\overline{\vect{u}}\right)\otimes\left(\vect{u} -\overline{\vect{u}}\right)}\) is modeled by \(\mu_t \overline{\tens{S}}\) and not by \(\mu_t\overline{\tens{S}}-\dfrac{2}{3}\mu_t \trace\left(\overline{\tens{S}}\right)\tens{1}+\dfrac{2}{3}k\tens{1}\) so that no term \(\mu_t \dive \left(\overline{\vect{u}}\right)\) is needed.

Please refer to the visecv section of the theory guide for more informations.

Function/Subroutine Documentation

◆ visecv()

subroutine visecv ( double precision, dimension(nfac)  secvif,
double precision, dimension(nfabor)  secvib 
)
Parameters
[in,out]secviflambda*surface at interior faces
[in,out]secviblambda*surface at boundary faces