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)$

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:

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

subroutine visecv	(	double precision, dimension(nfac)	secvif,
		double precision, dimension(nfabor)	secvib
	)

Parameters

[in,out]	secvif	lambda*surface at interior faces
[in,out]	secvib	lambda*surface at boundary faces