8.1
general documentation
Functions | Variables
cs_turbulence_model.c File Reference
#include "cs_defs.h"
#include <assert.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include "bft_mem.h"
#include "bft_error.h"
#include "bft_printf.h"
#include "cs_assert.h"
#include "cs_field.h"
#include "cs_field_pointer.h"
#include "cs_function.h"
#include "cs_gradient.h"
#include "cs_log.h"
#include "cs_log_iteration.h"
#include "cs_math.h"
#include "cs_mesh.h"
#include "cs_parall.h"
#include "cs_mesh_location.h"
#include "cs_time_step.h"
#include "cs_wall_functions.h"
#include "cs_turbulence_model.h"
+ Include dependency graph for cs_turbulence_model.c:

Functions

void cs_turb_model_init (void)
 Initialize turbulence model structures. More...
 
void cs_set_type_order_turbulence_model (void)
 Initialize type and order members of turbulence model structure. More...
 
cs_turb_model_t * cs_get_glob_turb_model (void)
 Provide write access to turbulence model structure. More...
 
void cs_set_glob_turb_model (void)
 Set global pointer to turbulence model structure. More...
 
void cs_turb_compute_constants (int phase_id)
 Compute turbulence model constants, some of which may depend on the model choice. More...
 
cs_turb_ref_values_t * cs_get_glob_turb_ref_values (void)
 Provide access to cs_glob_turb_ref_values. More...
 
cs_turb_rans_model_t * cs_get_glob_turb_rans_model (void)
 Provide access to cs_glob_turb_rans_model. More...
 
cs_turb_les_model_t * cs_get_glob_turb_les_model (void)
 Provide access to cs_glob_turb_les_model. More...
 
cs_turb_hybrid_model_t * cs_get_glob_turb_hybrid_model (void)
 Provide access to cs_glob_turb_hybrid_model. More...
 
void cs_turb_model_log_setup (void)
 Print the turbulence model parameters to setup.log. More...
 
void cs_turb_constants_log_setup (void)
 Print the turbulent constants to setup.log. More...
 
void cs_turb_init_ref_quantities (void)
 Compute characteristic length for turbulence if not already done. More...
 
void cs_clip_turbulent_fluxes (int flux_id, int variance_id)
 Clipping for the turbulence flux vector. More...
 
void cs_turbulence_function_k (int location_id, cs_lnum_t n_elts, const cs_lnum_t *elt_ids, void *input, void *vals)
 Return or estimate the value of the turbulent kinetic energy over specified elements. More...
 
void cs_turbulence_function_eps (int location_id, cs_lnum_t n_elts, const cs_lnum_t *elt_ids, void *input, void *vals)
 Return or estimate the value of the turbulent dissipation over specified elements. More...
 
void cs_turbulence_function_rij (int location_id, cs_lnum_t n_elts, const cs_lnum_t *elt_ids, void *input, void *vals)
 Return or estimate the value of the Reynolds stresses over specified elements. More...
 

Variables

double cs_turb_xkappa = 0.42
 
double cs_turb_vdriest = 25.6
 
double cs_turb_cstlog = 5.2
 
double cs_turb_cstlog_rough = 8.5
 
double cs_turb_cstlog_alpha
 
double cs_turb_apow = 8.3
 
double cs_turb_bpow = 1.0/7.0
 
double cs_turb_dpow = -1.
 
double cs_turb_cmu = 0.09
 
double cs_turb_cmu025 = 0.547722557
 
double cs_turb_ce1 = 1.44
 
double cs_turb_ce2 = 1.92
 
double cs_turb_ce4 = 1.20
 
double cs_turb_crij1 = 1.80
 
double cs_turb_crij2 = 0.60
 
double cs_turb_crij_c0 = 3.5
 
double cs_turb_crij3 = 0.55
 
double cs_turb_crijp1 = 0.50
 
double cs_turb_crijp2 = 0.30
 
double cs_turb_cssge2 = 1.83
 
double cs_turb_cssgs1 = 1.70
 
double cs_turb_cssgs2 = -1.05
 
double cs_turb_cssgr1 = 0.90
 
double cs_turb_cssgr2 = 0.80
 
double cs_turb_cssgr3 = 0.65
 
double cs_turb_cssgr4 = 0.625
 
double cs_turb_cssgr5 = 0.20
 
double cs_turb_cebms1 = 1.70
 
double cs_turb_cebms2 = 0.
 
double cs_turb_cebmr1 = 0.90
 
double cs_turb_cebmr2 = 0.80
 
double cs_turb_cebmr3 = 0.65
 
double cs_turb_cebmr4 = 0.625
 
double cs_turb_cebmr5 = 0.20
 
double cs_turb_csrij
 
double cs_turb_cebme2 = 1.83
 
double cs_turb_cebmmu = 0.22
 
double cs_turb_xcl = 0.122
 
double cs_turb_xa1 = 0.1
 
double cs_turb_xct = 6.0
 
double cs_turb_xceta = 80.0
 
double cs_turb_cpale1 = 1.44
 
double cs_turb_cpale2 = 1.83
 
double cs_turb_cpale3 = 2.3
 
double cs_turb_cpale4 = 0.4
 
double cs_turb_cpalc1 = 1.7
 
double cs_turb_cpalc2 = 0.9
 
double cs_turb_cpalct = 4.0
 
double cs_turb_cpalcl = 0.164
 
double cs_turb_cpalet = 75.0
 
double cs_turb_ckwsk1 = 1.0/0.85
 
double cs_turb_ckwsk2 = 1.0
 
double cs_turb_ckwsw1 = 2.0
 
double cs_turb_ckwsw2 = 1.0/0.856
 
double cs_turb_ckwbt1 = 0.075
 
double cs_turb_ckwbt2 = 0.0828
 
double cs_turb_ckwgm1 = -1.
 
double cs_turb_ckwgm2 = -1.
 
double cs_turb_ckwa1 = 0.31
 
double cs_turb_ckwc1 = 10.0
 
double cs_turb_cddes = 0.65
 
double cs_turb_csas = 0.11
 
double cs_turb_csas_eta2 = 3.51
 
double cs_turb_chtles_bt0 = 0.48
 
double cs_turb_csab1 = 0.1355
 
double cs_turb_csab2 = 0.622
 
double cs_turb_csasig = 2.0/3.0
 
double cs_turb_csav1 = 7.1
 
double cs_turb_csaw1 = -1.
 
double cs_turb_csaw2 = 0.3
 
double cs_turb_csaw3 = 2.0
 
double cs_turb_cssr1 = 1.0
 
double cs_turb_cssr2 = 2.0
 
double cs_turb_cssr3 = 1.0
 
double cs_turb_ccaze2 = 1.83
 
double cs_turb_ccazsc = 0.119
 
double cs_turb_ccaza = 4.3
 
double cs_turb_ccazb = 5.130
 
double cs_turb_ccazc = 0.453
 
double cs_turb_ccazd = 0.682
 
double cs_turb_xlesfl = 2.0
 
double cs_turb_ales = 1.0
 
double cs_turb_bles = 1.0/3.0
 
double cs_turb_csmago = 0.065
 
double cs_turb_xlesfd = 1.5
 
double cs_turb_csmago_max = -1.
 
double cs_turb_csmago_min = 0.
 
double cs_turb_cdries = 26.0
 
double cs_turb_cv2fa1 = 0.05
 
double cs_turb_cv2fe2 = 1.85
 
double cs_turb_cv2fc1 = 1.4
 
double cs_turb_cv2fc2 = 0.3
 
double cs_turb_cv2fct = 6.0
 
double cs_turb_cv2fcl = 0.25
 
double cs_turb_cv2fet = 110.0
 
double cs_turb_cnl1 = 0.8
 
double cs_turb_cnl2 = 11.
 
double cs_turb_cnl3 = 4.5
 
double cs_turb_cnl4 = 1e3
 
double cs_turb_cnl5 = 1.
 
double cs_turb_cwale = 0.25
 
double cs_turb_xiafm = 0.7
 
double cs_turb_etaafm = 0.4
 
double cs_turb_c1trit = 4.15
 
double cs_turb_c2trit = 0.55
 
double cs_turb_c3trit = 0.5
 
double cs_turb_c4trit = 0.
 
double cs_turb_cthafm = 0.236
 
double cs_turb_cthdfm = 0.31
 
double cs_turb_cthebdfm = 0.22
 
double cs_turb_xclt = 0.305
 

Detailed Description

Base turbulence model data.

Function Documentation

◆ cs_clip_turbulent_fluxes()

void cs_clip_turbulent_fluxes ( int  flux_id,
int  variance_id 
)

Clipping for the turbulence flux vector.

Parameters
[in]flux_idturbulent flux index
[in]variance_idscalar variance index

◆ cs_get_glob_turb_hybrid_model()

cs_turb_hybrid_model_t* cs_get_glob_turb_hybrid_model ( void  )

Provide access to cs_glob_turb_hybrid_model.

needed to initialize structure with GUI

◆ cs_get_glob_turb_les_model()

cs_turb_les_model_t* cs_get_glob_turb_les_model ( void  )

Provide access to cs_glob_turb_les_model.

needed to initialize structure with GUI

◆ cs_get_glob_turb_model()

cs_turb_model_t* cs_get_glob_turb_model ( void  )

Provide write access to turbulence model structure.

◆ cs_get_glob_turb_rans_model()

cs_turb_rans_model_t* cs_get_glob_turb_rans_model ( void  )

Provide access to cs_glob_turb_rans_model.

needed to initialize structure with GUI

◆ cs_get_glob_turb_ref_values()

cs_turb_ref_values_t* cs_get_glob_turb_ref_values ( void  )

Provide access to cs_glob_turb_ref_values.

needed to initialize structure with GUI

◆ cs_set_glob_turb_model()

void cs_set_glob_turb_model ( void  )

Set global pointer to turbulence model structure.

This global pointer provides a read-only access to the structure.

◆ cs_set_type_order_turbulence_model()

void cs_set_type_order_turbulence_model ( void  )

Initialize type and order members of turbulence model structure.

◆ cs_turb_compute_constants()

void cs_turb_compute_constants ( int  phase_id)

Compute turbulence model constants, some of which may depend on the model choice.

Parameters
[in]phase_idturbulent phase id (-1 for single phase flow)

◆ cs_turb_constants_log_setup()

void cs_turb_constants_log_setup ( void  )

Print the turbulent constants to setup.log.

◆ cs_turb_init_ref_quantities()

void cs_turb_init_ref_quantities ( void  )

Compute characteristic length for turbulence if not already done.

◆ cs_turb_model_init()

void cs_turb_model_init ( void  )

Initialize turbulence model structures.

◆ cs_turb_model_log_setup()

void cs_turb_model_log_setup ( void  )

Print the turbulence model parameters to setup.log.

◆ cs_turbulence_function_eps()

void cs_turbulence_function_eps ( int  location_id,
cs_lnum_t  n_elts,
const cs_lnum_t *  elt_ids,
void *  input,
void *  vals 
)

Return or estimate the value of the turbulent dissipation over specified elements.

Returned values are zero for turbulence models other than RANS.

This function matches the cs_eval_at_location_t function profile.

Parameters
[in]location_idbase associated mesh location id
[in]n_eltsnumber of associated elements
[in]elt_idsids of associated elements, or NULL if no filtering is required
[in,out]inputignored
[in,out]valspointer to output values (size: n_elts*dimension)

◆ cs_turbulence_function_k()

void cs_turbulence_function_k ( int  location_id,
cs_lnum_t  n_elts,
const cs_lnum_t *  elt_ids,
void *  input,
void *  vals 
)

Return or estimate the value of the turbulent kinetic energy over specified elements.

Returned values are zero for turbulence models other than RANS.

This function matches the cs_eval_at_location_t function profile.

Parameters
[in]location_idbase associated mesh location id
[in]n_eltsnumber of associated elements
[in]elt_idsids of associated elements, or NULL if no filtering is required
[in,out]inputignored
[in,out]valspointer to output values (size: n_elts*dimension)

◆ cs_turbulence_function_rij()

void cs_turbulence_function_rij ( int  location_id,
cs_lnum_t  n_elts,
const cs_lnum_t *  elt_ids,
void *  input,
void *  vals 
)

Return or estimate the value of the Reynolds stresses over specified elements.

Returned values are zero for turbulence models other than RANS.

This function matches the cs_eval_at_location_t function profile.

Parameters
[in]location_idbase associated mesh location id
[in]n_eltsnumber of associated elements
[in]elt_idsids of associated elements, or NULL if no filtering is required
[in,out]inputignored
[in,out]valspointer to output values (size: n_elts*dimension)

Variable Documentation

◆ cs_turb_ales

double cs_turb_ales = 1.0

Constant used to define, for each cell $\Omega_i$, the width of the (implicit) filter:

  • $\overline{\Delta}=xlesfl(ales*|\Omega_i|)^{bles}$

Useful if and only if iturb = 40 or 41.

◆ cs_turb_apow

double cs_turb_apow = 8.3

Werner and Wengle coefficient

◆ cs_turb_bles

double cs_turb_bles = 1.0/3.0

Constant used to define, for each cell $Omega_i$, the width of the (implicit) filter:

  • $\overline{\Delta}=xlesfl(ales*|\Omega_i|)^{bles}$

Useful if and only if iturb = 40 or 41.

◆ cs_turb_bpow

double cs_turb_bpow = 1.0/7.0

Werner and Wengle coefficient

◆ cs_turb_c1trit

double cs_turb_c1trit = 4.15

Coefficient of turbulent DFM flow model.

◆ cs_turb_c2trit

double cs_turb_c2trit = 0.55

Coefficient of turbulent DFM flow model.

◆ cs_turb_c3trit

double cs_turb_c3trit = 0.5

Coefficient of turbulent DFM flow model.

◆ cs_turb_c4trit

double cs_turb_c4trit = 0.

Coefficient of turbulent DFM flow model.

◆ cs_turb_ccaza

double cs_turb_ccaza = 4.3

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazb

double cs_turb_ccazb = 5.130

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazc

double cs_turb_ccazc = 0.453

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazd

double cs_turb_ccazd = 0.682

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccaze2

double cs_turb_ccaze2 = 1.83

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazsc

double cs_turb_ccazsc = 0.119

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_cddes

double cs_turb_cddes = 0.65

Constant $ C_{DDES} $ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST) and hybrid_turb=1.

◆ cs_turb_cdries

double cs_turb_cdries = 26.0

Van Driest constant appearing in the van Driest damping function applied to the Smagorinsky constant:

  • $ (1-\exp^{(-y^+/cdries}) $.

Useful if and only if iturb = 40 or 41.

◆ cs_turb_ce1

double cs_turb_ce1 = 1.44

Constant $C_{\varepsilon 1}$ for all the RANS turbulence models except for the v2f and the $k-\omega$ models. Useful if and only if iturb= 20, 21, 30 or 31 ( $k-\varepsilon$ or $R_{ij}-\varepsilon$).

◆ cs_turb_ce2

double cs_turb_ce2 = 1.92

Constant $C_{\varepsilon 2}$ for the $k-\varepsilon$ and $R_{ij}-\varepsilon$ LRR models. Useful if and only if iturb = 20, 21 or 30 ( $k-\varepsilon$ or $R_{ij}-\varepsilon$ LRR).

◆ cs_turb_ce4

double cs_turb_ce4 = 1.20

Coefficient of interfacial coefficient in k-eps, used in Lagrange treatment.

Constant $C_{\varepsilon 4}$ for the interfacial term (Lagrangian module) in case of two-way coupling. Useful in case of Lagrangian modelling, in $k-\varepsilon$ and $R_{ij}-\varepsilon$ with two-way coupling.

◆ cs_turb_cebme2

double cs_turb_cebme2 = 1.83

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_cebmmu

double cs_turb_cebmmu = 0.22

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_cebmr1

double cs_turb_cebmr1 = 0.90

◆ cs_turb_cebmr2

double cs_turb_cebmr2 = 0.80

◆ cs_turb_cebmr3

double cs_turb_cebmr3 = 0.65

◆ cs_turb_cebmr4

double cs_turb_cebmr4 = 0.625

◆ cs_turb_cebmr5

double cs_turb_cebmr5 = 0.20

◆ cs_turb_cebms1

double cs_turb_cebms1 = 1.70

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_cebms2

double cs_turb_cebms2 = 0.

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_chtles_bt0

double cs_turb_chtles_bt0 = 0.48

Constant $ \beta_0 $ for the HTLES model. Useful if and only if iturb=60 ( $k-\omega$ SST) or if iturb=51 ( $BL-v^2-k$) and hybrid_turb=4.

◆ cs_turb_ckwa1

double cs_turb_ckwa1 = 0.31

Specific constant of k-omega SST. Constant $a_1$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

◆ cs_turb_ckwbt1

double cs_turb_ckwbt1 = 0.075

Constant $\beta_1$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

◆ cs_turb_ckwbt2

double cs_turb_ckwbt2 = 0.0828

Constant $\beta_2$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

◆ cs_turb_ckwc1

double cs_turb_ckwc1 = 10.0

Constant $ c_1 $ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST). Specific constant of k-omega SST.

◆ cs_turb_ckwgm1

double cs_turb_ckwgm1 = -1.

$\frac{\beta_1}{C_\mu}-\frac{\kappa^2}{\sqrt{C_\mu}\sigma_{\omega 1}}$. Constant $\gamma_1$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

Warning
: $\gamma_1$ is calculated before the call to usipsu. Hence, if $\beta_1$, $C_\mu$, $\kappa$ or $\sigma_{\omega 1}$ is modified in usipsu, cs_turb_ckwgm1 must also be modified in accordance.

◆ cs_turb_ckwgm2

double cs_turb_ckwgm2 = -1.

$\frac{\beta_2}{C_\mu}-\frac{\kappa^2}{\sqrt{C_\mu}\sigma_{\omega 2}}$. Constant $\gamma_2$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

Warning
: $\gamma_2$ is calculated before the call to usipsu. Hence, if $\beta_2$, $C_\mu$, $\kappa$ or $\sigma_{\omega 2}$ is modified in usipsu, cs_turb_ckwgm2 must also be modified in accordance.

◆ cs_turb_ckwsk1

double cs_turb_ckwsk1 = 1.0/0.85

Constant $\sigma_{k1}$ for the $k-\omega$ SST model. Useful if and only if iturb=60.

◆ cs_turb_ckwsk2

double cs_turb_ckwsk2 = 1.0

Constant $\sigma_{k2}$ for the $k-\omega$ SST model. Useful if and only if iturb=60.

◆ cs_turb_ckwsw1

double cs_turb_ckwsw1 = 2.0

Constant $\sigma_{\omega 1}$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

◆ cs_turb_ckwsw2

double cs_turb_ckwsw2 = 1.0/0.856

Constant $\sigma_{\omega 2}$ for the $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST).

◆ cs_turb_cmu

double cs_turb_cmu = 0.09

Constant $C_\mu$ for all the RANS turbulence models. Warning: different value for v2f models. Useful only for RANS models ( $k-\varepsilon$, $R_{ij}-\varepsilon$ or $k-\omega$).

◆ cs_turb_cmu025

double cs_turb_cmu025 = 0.547722557

$ C_\mu^\frac{1}{4} $

◆ cs_turb_cnl1

double cs_turb_cnl1 = 0.8

Constants for the Baglietto et al. quadratic k-epsilon model. Useful if and only if iturb = CS_TURB_K_EPSILON_QUAD

◆ cs_turb_cnl2

double cs_turb_cnl2 = 11.

◆ cs_turb_cnl3

double cs_turb_cnl3 = 4.5

◆ cs_turb_cnl4

double cs_turb_cnl4 = 1e3

◆ cs_turb_cnl5

double cs_turb_cnl5 = 1.

◆ cs_turb_cpalc1

double cs_turb_cpalc1 = 1.7

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpalc2

double cs_turb_cpalc2 = 0.9

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpalcl

double cs_turb_cpalcl = 0.164

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpalct

double cs_turb_cpalct = 4.0

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpale1

double cs_turb_cpale1 = 1.44

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpale2

double cs_turb_cpale2 = 1.83

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpale3

double cs_turb_cpale3 = 2.3

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpale4

double cs_turb_cpale4 = 0.4

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_cpalet

double cs_turb_cpalet = 75.0

Specific constant of v2f "BL-v2k" (or phi-alpha).

◆ cs_turb_crij1

double cs_turb_crij1 = 1.80

Constant $C_1$ for the $R_{ij}-\varepsilon$ LRR model. Useful if and only if iturb=30 ( $R_{ij}-\varepsilon$ LRR).

◆ cs_turb_crij2

double cs_turb_crij2 = 0.60

Constant $C_2$ for the $R_{ij}-\varepsilon$ LRR model. Useful if and only if iturb=30 ( $R_{ij}-\varepsilon$ LRR).

◆ cs_turb_crij3

double cs_turb_crij3 = 0.55

Constant $C_3$ for the $R_{ij}-\varepsilon$ models. Value is 0.55 for SSG and LRR, 0.6 for EBRSM.

◆ cs_turb_crij_c0

double cs_turb_crij_c0 = 3.5

Rotta constant $C_0$ for the $R_{ij}-\varepsilon$ model. Useful for the Lagrangian model. The value is set from $C_1$ if and only if iturb=CS_TURB_RIJ_EPSILON_LRR ( $R_{ij}-\varepsilon$ LRR) and $C_2=0$.

◆ cs_turb_crijp1

double cs_turb_crijp1 = 0.50

Constant $C_1^\prime$ for the $R_{ij}-\varepsilon$ LRR model, corresponding to the wall echo terms. Useful if and only if iturb=30 and cs_turb_rans_model_t::irijec=1 ( $R_{ij}-\varepsilon$ LRR).

◆ cs_turb_crijp2

double cs_turb_crijp2 = 0.30

Constant $C_2^\prime$ for the $R_{ij}-\varepsilon$ LRR model, corresponding to the wall echo terms. Useful if and only if iturb=30 and cs_turb_rans_model_t::irijec=1 ( $R_{ij}-\varepsilon$ LRR).

◆ cs_turb_csab1

double cs_turb_csab1 = 0.1355

Specific constant of Spalart-Allmaras.

◆ cs_turb_csab2

double cs_turb_csab2 = 0.622

Specific constant of Spalart-Allmaras.

◆ cs_turb_csas

double cs_turb_csas = 0.11

Constant $ C_{SAS}$ for the hybrid $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST) and hybrid_turb=3.

◆ cs_turb_csas_eta2

double cs_turb_csas_eta2 = 3.51

constant $ C_{DDES}$ for the hybrid $k-\omega$ SST model. Useful if and only if iturb=60 ( $k-\omega$ SST) and hybrid_turb=3.

◆ cs_turb_csasig

double cs_turb_csasig = 2.0/3.0

Specific constant of Spalart-Allmaras.

◆ cs_turb_csav1

double cs_turb_csav1 = 7.1

Specific constant of Spalart-Allmaras.

◆ cs_turb_csaw1

double cs_turb_csaw1 = -1.

Specific constant of Spalart-Allmaras.

◆ cs_turb_csaw2

double cs_turb_csaw2 = 0.3

Specific constant of Spalart-Allmaras.

◆ cs_turb_csaw3

double cs_turb_csaw3 = 2.0

Specific constant of Spalart-Allmaras.

◆ cs_turb_csmago

double cs_turb_csmago = 0.065

Smagorinsky constant used in the Smagorinsky model for LES. The sub-grid scale viscosity is calculated by $\displaystyle\mu_{sg}= \rho C_{smago}^2\bar{\Delta}^2\sqrt{2\bar{S}_{ij}\bar{S}_{ij}}$ where $\bar{\Delta}$ is the width of the filter and $\bar{S}_{ij}$ the filtered strain rate.

Useful if and only if iturb = 40.

Note
In theory Smagorinsky constant is 0.18. For a channel, 0.065 value is rather taken.

◆ cs_turb_csmago_max

double cs_turb_csmago_max = -1.

Maximum allowed value for the variable $C$ appearing in the LES dynamic model. Any larger value yielded by the calculation procedure of the dynamic model will be clipped to $ smagmx$.

Useful if and only if iturb = 41.

◆ cs_turb_csmago_min

double cs_turb_csmago_min = 0.

Minimum allowed value for the variable $C$ appearing in the LES dynamic model. Any smaller value yielded by the calculation procedure of the dynamic model will be clipped to $ smagmn$.

Useful if and only if iturb = 41.

◆ cs_turb_csrij

double cs_turb_csrij

Constant $C_s$ for the $R_{ij}-\varepsilon$ LRR model. Useful if and only if iturb=30 ( $R_{ij}-\varepsilon$ LRR).

◆ cs_turb_cssge2

double cs_turb_cssge2 = 1.83

Constant $C_{\varepsilon 2}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgr1

double cs_turb_cssgr1 = 0.90

Constant $C_{r1}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgr2

double cs_turb_cssgr2 = 0.80

Constant $C_{r2}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgr3

double cs_turb_cssgr3 = 0.65

Constant $C_{r3}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgr4

double cs_turb_cssgr4 = 0.625

constant $C_{r4}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgr5

double cs_turb_cssgr5 = 0.20

Constant $C_{r1}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgs1

double cs_turb_cssgs1 = 1.70

Constant $C_{s1}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssgs2

double cs_turb_cssgs2 = -1.05

Constant $C_{s2}$ for the $R_{ij}-\varepsilon$ SSG model. Useful if and only if iturb=31 ( $R_{ij}-\varepsilon$ SSG).

◆ cs_turb_cssr1

double cs_turb_cssr1 = 1.0

Constant of the Spalart-Shur rotation/curvature correction.

◆ cs_turb_cssr2

double cs_turb_cssr2 = 2.0

Constant of the Spalart-Shur rotation/curvature correction.

◆ cs_turb_cssr3

double cs_turb_cssr3 = 1.0

Constant of the Spalart-Shur rotation/curvature correction.

◆ cs_turb_cstlog

double cs_turb_cstlog = 5.2

Constant of logarithmic smooth law function: $ \dfrac{1}{\kappa} \ln(y^+) + cstlog $ ( $ cstlog = 5.2 $).

Constant of the logarithmic wall function. Useful if and only if iturb >= 10 (mixing length, $k-\varepsilon$, $R_{ij}-\varepsilon$, LES, v2f or $k-\omega$).

◆ cs_turb_cstlog_alpha

double cs_turb_cstlog_alpha

Constant $ \alpha $ for logarithmic law function switching from rough to smooth: $ \dfrac{1}{\kappa} \ln(y u_k/(\nu + \alpha \xi u_k)) + cstlog $ ( $ \alpha = \exp \left( -\kappa (8.5 - 5.2) \right) $).

Useful if and only if iturb >= 10 (mixing length, $k-\varepsilon$, $R_{ij}-\varepsilon$, LES, v2f or $k-\omega$).

◆ cs_turb_cstlog_rough

double cs_turb_cstlog_rough = 8.5

Constant of logarithmic rough law function: $ \dfrac{1}{\kappa} \ln(y/\xi) + cstlog_{rough} $ ( $ cstlog_{rough} = 8.5 $).

Constant of the logarithmic wall function. Useful if and only if iturb >= 10 (mixing length, $k-\varepsilon$, $R_{ij}-\varepsilon$, LES, v2f or $k-\omega$).

◆ cs_turb_cthafm

double cs_turb_cthafm = 0.236

Constant of GGDH and AFM on the thermal scalar.

◆ cs_turb_cthdfm

double cs_turb_cthdfm = 0.31

Constant of GGDH and AFM on the thermal scalar.

◆ cs_turb_cthebdfm

double cs_turb_cthebdfm = 0.22

◆ cs_turb_cv2fa1

double cs_turb_cv2fa1 = 0.05

Constant $a_1$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cv2fc1

double cs_turb_cv2fc1 = 1.4

Constant $C_1$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cv2fc2

double cs_turb_cv2fc2 = 0.3

Constant $C_2$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cv2fcl

double cs_turb_cv2fcl = 0.25

Constant $C_L$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cv2fct

double cs_turb_cv2fct = 6.0

Constant $C_T$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cv2fe2

double cs_turb_cv2fe2 = 1.85

Constant $C_{\varepsilon 2}$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cv2fet

double cs_turb_cv2fet = 110.0

Constant $C_\eta$ for the v2f $\varphi$-model. Useful if and only if iturb=50 (v2f $\varphi$-model).

◆ cs_turb_cwale

double cs_turb_cwale = 0.25

Constant of the WALE LES method.

◆ cs_turb_dpow

double cs_turb_dpow = -1.

Werner and Wengle coefficient

◆ cs_turb_etaafm

double cs_turb_etaafm = 0.4

Coefficient of turbulent AFM flow model.

◆ cs_turb_vdriest

double cs_turb_vdriest = 25.6

Van Driest constant. (= 25.6)

Useful if and only if cs_glob_wall_functions::iwallf = CS_WALL_F_2SCALES_VDRIEST. (Two scales log law at the wall using Van Driest mixing length expression).

◆ cs_turb_xa1

double cs_turb_xa1 = 0.1

Constant in the expression of Ce1' for the Rij-epsilon EBRSM.

◆ cs_turb_xceta

double cs_turb_xceta = 80.0

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_xcl

double cs_turb_xcl = 0.122

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_xclt

double cs_turb_xclt = 0.305

constant of EB-AFM and EB-DFM (0.122*2.5, See F. Dehoux thesis)

◆ cs_turb_xct

double cs_turb_xct = 6.0

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_xiafm

double cs_turb_xiafm = 0.7

Coefficient of turbulent AFM flow model.

◆ cs_turb_xkappa

double cs_turb_xkappa = 0.42

Karman constant. (= 0.42)

Useful if and only if iturb >= 10. (mixing length, $k-\varepsilon$, $R_{ij}-\varepsilon$, LES, v2f or $k-\omega$).

◆ cs_turb_xlesfd

double cs_turb_xlesfd = 1.5

Ratio between explicit and explicit filter width for a dynamic model. Constant used to define, for each cell $\Omega_i$, the width of the explicit filter used in the framework of the LES dynamic model: $\widetilde{\overline{\Delta}}=xlesfd\overline{\Delta}$.

Useful if and only if iturb = 41.

◆ cs_turb_xlesfl

double cs_turb_xlesfl = 2.0

Constant used in the definition of LES filtering diameter: $ \delta = \text{xlesfl} . (\text{ales} . volume)^{\text{bles}}$ cs_turb_xlesfl is a constant used to define, for each cell $\Omega_i$, the width of the (implicit) filter: $\overline{\Delta}=xlesfl(ales*|\Omega_i|)^{bles}$
Useful if and only if iturb = 40 or 41