#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_field.h"
#include "cs_field_pointer.h"
#include "cs_log.h"
#include "cs_map.h"
#include "cs_parall.h"
#include "cs_parameters.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
cs_turb_model_t *	cs_get_glob_turb_model (void)
	Provide access to global turbulence model structure cs_glob_turb_model. More...

void	cs_turb_compute_constants (void)
	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 acces to cs_glob_turb_ref_values. More...

cs_turb_rans_model_t *	cs_get_glob_turb_rans_model (void)
	Provide acces to cs_glob_turb_rans_model. More...

cs_turb_les_model_t *	cs_get_glob_turb_les_model (void)
	Provide acces to cs_glob_turb_les_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...

Variables
const double	cs_turb_xkappa = 0.42

const double	cs_turb_vdriest = 25.6

const double	cs_turb_cstlog = 5.2

const double	cs_turb_cstlog_rough = 8.5

double	cs_turb_cstlog_alpha

const double	cs_turb_apow = 8.3

const double	cs_turb_bpow = 1.0/7.0

double	cs_turb_dpow = -1.

double	cs_turb_cmu = 0.09

double	cs_turb_cmu025

const double	cs_turb_ce1 = 1.44

const double	cs_turb_ce2 = 1.92

const double	cs_turb_ce4 = 1.20

const double	cs_turb_sigmak = 1.0

double	cs_turb_sigmae = -1.

double	cs_turb_crij1 = 1.80

double	cs_turb_crij2 = 0.60

double	cs_turb_crij3 = 0.55

const double	cs_turb_crijp1 = 0.50

const double	cs_turb_crijp2 = 0.30

const double	cs_turb_cssge2 = 1.83

const double	cs_turb_cssgs1 = 1.70

const double	cs_turb_cssgs2 = -1.05

const double	cs_turb_cssgr1 = 0.90

const double	cs_turb_cssgr2 = 0.80

const double	cs_turb_cssgr3 = 0.65

const double	cs_turb_cssgr4 = 0.625

const double	cs_turb_cssgr5 = 0.20

const double	cs_turb_cebms1 = 1.70

const double	cs_turb_cebms2 = 0.

const double	cs_turb_cebmr1 = 0.90

const double	cs_turb_cebmr2 = 0.80

const double	cs_turb_cebmr3 = 0.65

const double	cs_turb_cebmr4 = 0.625

const double	cs_turb_cebmr5 = 0.20

double	cs_turb_csrij

const double	cs_turb_cebme2 = 1.83

const double	cs_turb_cebmmu = 0.22

const double	cs_turb_xcl = 0.122

const double	cs_turb_xa1 = 0.1

const double	cs_turb_xct = 6.0

const double	cs_turb_xceta = 80.0

const double	cs_turb_cpale1 = 1.44

const double	cs_turb_cpale2 = 1.83

const double	cs_turb_cpale3 = 2.3

const double	cs_turb_cpale4 = 0.4

const double	cs_turb_cpalse = 1.5

const double	cs_turb_cpalmu = 0.22

const double	cs_turb_cpalc1 = 1.7

const double	cs_turb_cpalc2 = 0.9

const double	cs_turb_cpalct = 4.0

const double	cs_turb_cpalcl = 0.164

const double	cs_turb_cpalet = 75.0

const double	cs_turb_ckwsk1 = 1.0/0.85

const double	cs_turb_ckwsk2 = 1.0

const double	cs_turb_ckwsw1 = 2.0

const double	cs_turb_ckwsw2 = 1.0/0.856

const double	cs_turb_ckwbt1 = 0.075

const double	cs_turb_ckwbt2 = 0.0828

double	cs_turb_ckwgm1 = -1.

double	cs_turb_ckwgm2 = -1.

const double	cs_turb_ckwa1 = 0.31

const double	cs_turb_ckwc1 = 10.0

const double	cs_turb_csab1 = 0.1355

const double	cs_turb_csab2 = 0.622

const double	cs_turb_csasig = 2.0/3.0

const double	cs_turb_csav1 = 7.1

double	cs_turb_csaw1 = -1.

const double	cs_turb_csaw2 = 0.3

const double	cs_turb_csaw3 = 2.0

const double	cs_turb_cssr1 = 1.0

const double	cs_turb_cssr2 = 2.0

const double	cs_turb_cssr3 = 1.0

const double	cs_turb_ccaze2 = 1.83

const double	cs_turb_ccazsc = 0.119

const double	cs_turb_ccaza = 4.3

const double	cs_turb_ccazb = 5.130

const double	cs_turb_ccazc = 0.453

const double	cs_turb_ccazd = 0.682

const double	cs_turb_xlesfl = 2.0

const double	cs_turb_ales = 1.0

const double	cs_turb_bles = 1.0/3.0

double	cs_turb_csmago = 0.065

const double	cs_turb_xlesfd = 1.5

double	cs_turb_smagmx = -1.

double	cs_turb_smagmn = 0.

const double	cs_turb_cdries = 26.0

const double	cs_turb_cv2fa1 = 0.05

const double	cs_turb_cv2fe2 = 1.85

const double	cs_turb_cv2fmu = 0.22

const double	cs_turb_cv2fc1 = 1.4

const double	cs_turb_cv2fc2 = 0.3

const double	cs_turb_cv2fct = 6.0

const double	cs_turb_cv2fcl = 0.25

const double	cs_turb_cv2fet = 110.0

double	cs_turb_cwale = 0.25

const double	cs_turb_xiafm = 0.7

const double	cs_turb_etaafm = 0.4

const double	cs_turb_c1trit = 4.15

const double	cs_turb_c2trit = 0.55

const double	cs_turb_c3trit = 0.5

const double	cs_turb_c4trit = 0.

const double	cs_turb_cthafm = 0.236

const double	cs_turb_cthdfm = 0.31

Detailed Description

Base turbulence model data.

Function Documentation

◆ cs_get_glob_turb_les_model()

cs_turb_les_model_t* cs_get_glob_turb_les_model ( void )

Provide acces 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 access to global turbulence model structure cs_glob_turb_model.

It is needed to initialize structure with GUI.

◆ cs_get_glob_turb_rans_model()

cs_turb_rans_model_t* cs_get_glob_turb_rans_model ( void )

Provide acces 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 acces to cs_glob_turb_ref_values.

needed to initialize structure with GUI

◆ cs_turb_compute_constants()

void cs_turb_compute_constants ( void )

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

◆ cs_turb_constants_log_setup()

void cs_turb_constants_log_setup ( void )

Print the turbulent constants to setup.log.

◆ cs_turb_model_log_setup()

void cs_turb_model_log_setup ( void )

Print the turbulence model parameters to setup.log.

Variable Documentation

◆ cs_turb_ales

const 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

const double cs_turb_apow = 8.3

Werner and Wengle coefficient

◆ cs_turb_bles

const 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

const double cs_turb_bpow = 1.0/7.0

Werner and Wengle coefficient

◆ cs_turb_c1trit

const double cs_turb_c1trit = 4.15

Coefficient of turbulent DFM flow model.

◆ cs_turb_c2trit

const double cs_turb_c2trit = 0.55

Coefficient of turbulent DFM flow model.

◆ cs_turb_c3trit

const double cs_turb_c3trit = 0.5

Coefficient of turbulent DFM flow model.

◆ cs_turb_c4trit

const double cs_turb_c4trit = 0.

Coefficient of turbulent DFM flow model.

◆ cs_turb_ccaza

const double cs_turb_ccaza = 4.3

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazb

const double cs_turb_ccazb = 5.130

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazc

const double cs_turb_ccazc = 0.453

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazd

const double cs_turb_ccazd = 0.682

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccaze2

const double cs_turb_ccaze2 = 1.83

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_ccazsc

const double cs_turb_ccazsc = 0.119

Constants of the Cazalbou rotation/curvature correction.

◆ cs_turb_cdries

const 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

const 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

const 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

const 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

const double cs_turb_cebme2 = 1.83

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_cebmmu

const double cs_turb_cebmmu = 0.22

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_cebmr1

const double cs_turb_cebmr1 = 0.90

◆ cs_turb_cebmr2

const double cs_turb_cebmr2 = 0.80

◆ cs_turb_cebmr3

const double cs_turb_cebmr3 = 0.65

◆ cs_turb_cebmr4

const double cs_turb_cebmr4 = 0.625

◆ cs_turb_cebmr5

const double cs_turb_cebmr5 = 0.20

◆ cs_turb_cebms1

const double cs_turb_cebms1 = 1.70

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_cebms2

const double cs_turb_cebms2 = 0.

Constant of the Rij-epsilon EBRSM.

◆ cs_turb_ckwa1

const 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

const 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

const 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

const 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

const 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

const 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

const 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

const 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 except for the v2f model (see cs_turb_cv2fmu for the value of $C_\mu$ in case of v2f modelling). Useful if and only if iturb = 20, 21, 30, 31 or 60 ( $k-\varepsilon$ , $R_{ij}-\varepsilon$ or $k-\omega$ ).

◆ cs_turb_cmu025

double cs_turb_cmu025

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