#include "base/cs_defs.h"
#include <assert.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include "bft/bft_printf.h"
#include "base/cs_array.h"
#include "base/cs_base.h"
#include "base/cs_field.h"
#include "base/cs_field_default.h"
#include "base/cs_field_pointer.h"
#include "base/cs_field_operator.h"
#include "alge/cs_gradient.h"
#include "base/cs_halo.h"
#include "base/cs_halo_perio.h"
#include "turb/cs_les_filter.h"
#include "base/cs_log.h"
#include "base/cs_math.h"
#include "base/cs_mem.h"
#include "mesh/cs_mesh.h"
#include "mesh/cs_mesh_location.h"
#include "mesh/cs_mesh_quantities.h"
#include "base/cs_parall.h"
#include "base/cs_physical_constants.h"
#include "pprt/cs_physical_model.h"
#include "base/cs_prototypes.h"
#include "turb/cs_turbulence_model.h"
#include "turb/cs_les_mu_t.h"

Include dependency graph for cs_les_mu_t.cpp:

Functions
void	cs_les_mu_t_smago_dyn_prepare (cs_real_t s_n[], cs_real_t sf_n[], cs_real_t f_vel[][3], cs_real_t mijmij[], cs_real_t mijlij[])
	Calculation of Mij:Mij and Mij:Lij for dynamic Smagorinsky model.
void	cs_les_mu_t_smago_dyn (void)
	Calculation of the turbulent viscosity for a dynamic Smagorinsky LES model.
void	cs_les_mu_t_smago_const (void)
	Calculation of turbulent viscosity for a Smagorinsky LES model.
void	cs_les_mu_t_wale (void)
	Compute the turbulent viscosity for the WALE LES model.
void	cs_les_mu_t_ksgs (void)
	Calculation of turbulent viscosity for a k-SGS LES model !
void	cs_les_mu_t_tausgs (void)
	Calculation of turbulent viscosity for a tau-SGS LES model !

Variables
static const cs_lnum_t	_iv2t [6] = {0, 1, 2, 0, 1, 0}
static const cs_lnum_t	_jv2t [6] = {0, 1, 2, 1, 2, 2}

Function Documentation

◆ cs_les_mu_t_ksgs()

void cs_les_mu_t_ksgs ( void )

Calculation of turbulent viscosity for a k-SGS LES model !

◆ cs_les_mu_t_smago_const()

void cs_les_mu_t_smago_const ( void )

Calculation of turbulent viscosity for a Smagorinsky LES model.

Calculation of the turbulent viscosity for a Smagorinsky LES model.

$\mu_T = \rho (C_{S} l)^2 \sqrt{2 S_{ij}S_{ij}}$

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

!

◆ cs_les_mu_t_smago_dyn()

void cs_les_mu_t_smago_dyn ( void )

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

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

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

!

◆ cs_les_mu_t_smago_dyn_prepare()

void cs_les_mu_t_smago_dyn_prepare	(	cs_real_t	s_n[],
		cs_real_t	sf_n[],
		cs_real_t	f_vel[][3],
		cs_real_t	mijmij[],
		cs_real_t	mijlij[] )

Calculation of Mij:Mij and Mij:Lij for dynamic Smagorinsky model.

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

Parameters

[out]	s_n	strain rate (sqrt(2SijSij))
[out]	sf_n	filtered strain rate
[out]	f_vel	filtered velocity
[out]	mijmij	Mij:Mij
[out]	mijlij	Mij:Lij

!

◆ cs_les_mu_t_tausgs()

void cs_les_mu_t_tausgs ( void )

Calculation of turbulent viscosity for a tau-SGS LES model !

◆ cs_les_mu_t_wale()

void cs_les_mu_t_wale ( void )

Compute the turbulent viscosity for the WALE LES model.

The turbulent viscosity is: $\mu_T = \rho (C_{wale} L)^2 * \dfrac{(\tens{S}:\tens{Sd})^{3/2}} {(\tens{S} :\tens{S})^(5/2) +(\tens{Sd}:\tens{Sd})^(5/4)}$ with $\tens{S} = \frac{1}{2}(\gradt \vect{u} + \transpose{\gradt \vect{u}})$ and $\tens{Sd} = \deviator{(\symmetric{(\tens{S}^2)})}$