8.1
general documentation
Functions/Subroutines
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 ( real(c_double), dimension(*)  secvif,
real(c_double), dimension(*)  secvib 
)
Parameters
[in,out]secviflambda*surface at interior faces
[in,out]secviblambda*surface at boundary faces