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, smagor)

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

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(ncepdp,6)	ckupdc,
		double precision, dimension(ncesmp,nvar)	smacel,
		double precision, dimension(ncelet)	smagor
	)

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
[in,out]	smagor	smagorinski's constant in the cas of a dynamic model