Calculation of turbulent viscosity for a dynamic Smagorinsky LES model. More...

Functions/Subroutines
subroutine	visdyn (nvar, nscal, ncepdp, ncesmp, icepdc, icetsm, itypsm, dt, ckupdc, smacel)

Detailed Description

Calculation of turbulent viscosity for a dynamic Smagorinsky LES model.

$smago = \dfrac{L_{ij}M_{ij}}{M_{ij}M_{ij}}$

$\mu_T = \rho smago L^2 \sqrt{2 S_{ij}S_{ij}}$

$S_{ij} = \dfrac{\der{u_i}{x_j} + \der{u_j}{x_i}}{2}$

We have at edge faces types at previous time step (except at first time step, when tables itypfb and itrifb have not been filled).

Please refer to the dynamic Smagorinsky model section of the theory guide for more informations.

Function/Subroutine Documentation

subroutine visdyn	(	integer	nvar,
		integer	nscal,
		integer	ncepdp,
		integer	ncesmp,
		integer, dimension(ncepdp)	icepdc,
		integer, dimension(ncesmp)	icetsm,
		integer, dimension(ncesmp,nvar)	itypsm,
		double precision, dimension(ncelet)	dt,
		double precision, dimension(6,ncepdp)	ckupdc,
		double precision, dimension(ncesmp,nvar)	smacel
	)

Parameters

[in]	nvar	total number of variables
[in]	nscal	total number of scalars
[in]	ncepdp	number of cells with head loss
[in]	ncesmp	number of cells with mass source term
[in]	icepdc	number of ncepdp cells with losses
[in]	icetsm	number of cells with mass source
[in]	itypsm	type of mass source for the variable (cf. cs_user_mass_source_terms)
[in]	dt	time step (per cell)
[in]	ckupdc	work array for head losses
[in]	smacel	value of variables associated to the mass source for ivar = ipr, smacel = mass flux