9.1
general documentation
Functions
cs_turbulence_bc.cpp File Reference
#include "base/cs_defs.h"
#include <assert.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include "bft/bft_error.h"
#include "bft/bft_printf.h"
#include "base/cs_log.h"
#include "base/cs_field.h"
#include "base/cs_field_pointer.h"
#include "base/cs_map.h"
#include "base/cs_math.h"
#include "base/cs_mem.h"
#include "base/cs_parall.h"
#include "base/cs_param_types.h"
#include "mesh/cs_mesh_location.h"
#include "turb/cs_turbulence_model.h"
#include "turb/cs_turbulence_bc.h"
+ Include dependency graph for cs_turbulence_bc.cpp:

Functions

void cs_turbulence_bc_init_pointers (void)
 Initialize turbulence model boundary condition pointers. More...
 
void cs_turbulence_bc_free_pointers (void)
 Free memory allocations for turbulence boundary condition pointers. More...
 
void cs_turbulence_bc_ke_hyd_diam (double uref2, double dh, double rho, double mu, double *ustar2, double *k, double *eps)
 Calculation of $ u^\star $, $ k $ and $\varepsilon $ from a diameter $ D_H $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall (use for inlet boundary conditions). More...
 
void cs_turbulence_bc_ke_turb_intensity (double uref2, double t_intensity, double dh, double *k, double *eps)
 Calculation of $ k $ and $\varepsilon$ from a diameter $ D_H $, a turbulent intensity $ I $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall (for inlet boundary conditions). More...
 
void cs_turbulence_bc_inlet_hyd_diam (cs_lnum_t face_id, double uref2, double dh, double rho, double mu)
 Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall. More...
 
void cs_turbulence_bc_inlet_turb_intensity (cs_lnum_t face_id, double uref2, double t_intensity, double dh)
 Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $, a turbulent intensity $ I $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall. More...
 
void cs_turbulence_bc_inlet_k_eps (cs_lnum_t face_id, double k, double eps)
 Set inlet boundary condition values for turbulence variables based on given k and epsilon values. More...
 
void cs_turbulence_bc_set_hmg_neumann (cs_lnum_t face_id)
 Assign homogeneous Neumann turbulent boundary conditions to a given face. More...
 
void cs_turbulence_bc_set_uninit_inlet_hyd_diam (cs_lnum_t face_id, double vel_dir[], double shear_dir[], double uref2, double dh, double rho, double mu)
 Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall, only if not already set. More...
 
void cs_turbulence_bc_set_uninit_inlet_turb_intensity (cs_lnum_t face_id, double uref2, double t_intensity, double dh)
 Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $, a turbulent intensity $ I $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall, only if not already set. More...
 
void cs_turbulence_bc_set_uninit_inlet_k_eps (cs_lnum_t face_id, double k, double eps, double vel_dir[], double shear_dir[])
 Set inlet boundary condition values for turbulence variables based on given k and epsilon values only if not already initialized. More...
 

Detailed Description

Base turbulence boundary conditions.

Function Documentation

◆ cs_turbulence_bc_free_pointers()

void cs_turbulence_bc_free_pointers ( void  )

Free memory allocations for turbulence boundary condition pointers.

◆ cs_turbulence_bc_init_pointers()

void cs_turbulence_bc_init_pointers ( void  )

Initialize turbulence model boundary condition pointers.

◆ cs_turbulence_bc_inlet_hyd_diam()

void cs_turbulence_bc_inlet_hyd_diam ( cs_lnum_t  face_id,
double  uref2,
double  dh,
double  rho,
double  mu 
)

Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall.

We use the laws from Idel'Cik, i.e. the head loss coefficient $ \lambda $ is defined by:

\[ |\dfrac{\Delta P}{\Delta x}| = \dfrac{\lambda}{D_H} \frac{1}{2} \rho U_{ref}^2 \]

then the relation reads $u^\star = U_{ref} \sqrt{\dfrac{\lambda}{8}}$. $\lambda $ depends on the hydraulic Reynolds number $ Re = \dfrac{U_{ref} D_H}{ \nu} $ and is given by:

  • for $ Re < 2000 $

    \[ \lambda = \dfrac{64}{Re} \]

  • for $ Re > 4000 $

    \[ \lambda = \dfrac{1}{( 1.8 \log_{10}(Re)-1.64 )^2} \]

  • for $ 2000 < Re < 4000 $, we complete by a straight line

    \[ \lambda = 0.021377 + 5.3115. 10^{-6} Re \]

From $ u^\star $, we can estimate $ k $ and $ \varepsilon$ from the well known formulae of developped turbulence

Parameters
[in]face_idboundary face id
[in]uref2square of the reference flow velocity
[in]dhhydraulic diameter $ D_H $
[in]rhomass density $ \rho $
[in]mudynamic viscosity $ \nu $

◆ cs_turbulence_bc_inlet_k_eps()

void cs_turbulence_bc_inlet_k_eps ( cs_lnum_t  face_id,
double  k,
double  eps 
)

Set inlet boundary condition values for turbulence variables based on given k and epsilon values.

Parameters
[in]face_idboundary face id
[in]kturbulent kinetic energy
[in]epsturbulent dissipation

◆ cs_turbulence_bc_inlet_turb_intensity()

void cs_turbulence_bc_inlet_turb_intensity ( cs_lnum_t  face_id,
double  uref2,
double  t_intensity,
double  dh 
)

Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $, a turbulent intensity $ I $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall.

Parameters
[in]face_idboundary face id
[in]uref2square of the reference flow velocity
[in]t_intensityturbulent intensity $ I $
[in]dhhydraulic diameter $ D_H $

◆ cs_turbulence_bc_ke_hyd_diam()

void cs_turbulence_bc_ke_hyd_diam ( double  uref2,
double  dh,
double  rho,
double  mu,
double *  ustar2,
double *  k,
double *  eps 
)

Calculation of $ u^\star $, $ k $ and $\varepsilon $ from a diameter $ D_H $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall (use for inlet boundary conditions).

Both $ u^\star $ and $ (k,\varepsilon )$ are returned, so that the user may compute other values of $ k $ and $ \varepsilon $ with $ u^\star $.

We use the laws from Idel'Cik, i.e. the head loss coefficient $ \lambda $ is defined by:

\[ |\dfrac{\Delta P}{\Delta x}| = \dfrac{\lambda}{D_H} \frac{1}{2} \rho U_{ref}^2 \]

then the relation reads $u^\star = U_{ref} \sqrt{\dfrac{\lambda}{8}}$. $\lambda $ depends on the hydraulic Reynolds number $ Re = \dfrac{U_{ref} D_H}{ \nu} $ and is given by:

  • for $ Re < 2000 $

    \[ \lambda = \dfrac{64}{Re} \]

  • for $ Re > 4000 $

    \[ \lambda = \dfrac{1}{( 1.8 \log_{10}(Re)-1.64 )^2} \]

  • for $ 2000 < Re < 4000 $, we complete by a straight line

    \[ \lambda = 0.021377 + 5.3115. 10^{-6} Re \]

From $ u^\star $, we can estimate $ k $ and $ \varepsilon$ from the well known formulae of developped turbulence

\[ k = \dfrac{u^{\star 2}}{\sqrt{C_\mu}} \]

\[ \varepsilon = \dfrac{ u^{\star 3}}{(\kappa D_H /10)} \]

Parameters
[in]uref2square of the reference flow velocity
[in]dhhydraulic diameter $ D_H $
[in]rhomass density $ \rho $
[in]mudynamic viscosity $ \nu $
[out]ustar2square of friction speed
[out]kcalculated turbulent intensity $ k $
[out]epscalculated turbulent dissipation $ \varepsilon $

◆ cs_turbulence_bc_ke_turb_intensity()

void cs_turbulence_bc_ke_turb_intensity ( double  uref2,
double  t_intensity,
double  dh,
double *  k,
double *  eps 
)

Calculation of $ k $ and $\varepsilon$ from a diameter $ D_H $, a turbulent intensity $ I $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall (for inlet boundary conditions).

\[ k = 1.5 I {U_{ref}}^2 \]

\[ \varepsilon = 10 \dfrac{{C_\mu}^{0.75} k^{1.5}}{ \kappa D_H} \]

Parameters
[in]uref2square of the reference flow velocity
[in]t_intensityturbulent intensity $ I $
[in]dhhydraulic diameter $ D_H $
[out]kcalculated turbulent intensity $ k $
[out]epscalculated turbulent dissipation $ \varepsilon $

◆ cs_turbulence_bc_set_hmg_neumann()

void cs_turbulence_bc_set_hmg_neumann ( cs_lnum_t  face_id)

Assign homogeneous Neumann turbulent boundary conditions to a given face.

This is useful for outgoing flow.

◆ cs_turbulence_bc_set_uninit_inlet_hyd_diam()

void cs_turbulence_bc_set_uninit_inlet_hyd_diam ( cs_lnum_t  face_id,
double  vel_dir[],
double  shear_dir[],
double  uref2,
double  dh,
double  rho,
double  mu 
)

Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall, only if not already set.

Apart from assigning values where not already initialized, this function is similar to cs_turbulence_bc_inlet_hyd_diam.

Parameters
[in]face_idboundary face id
[in]vel_dirvelocity direction (not normalized, or nullptr)
[in]shear_dirshear direction (or nullptr), it also contains the level of anisotropy (Rnt = a_nt k)
[in]uref2square of the reference flow velocity
[in]dhhydraulic diameter $ D_H $
[in]rhomass density $ \rho $
[in]mudynamic viscosity $ \nu $

◆ cs_turbulence_bc_set_uninit_inlet_k_eps()

void cs_turbulence_bc_set_uninit_inlet_k_eps ( cs_lnum_t  face_id,
double  k,
double  eps,
double  vel_dir[],
double  shear_dir[] 
)

Set inlet boundary condition values for turbulence variables based on given k and epsilon values only if not already initialized.

Parameters
[in]face_idboundary face id
[in]kturbulent kinetic energy
[in]epsturbulent dissipation
[in]vel_dirvelocity direction
[in]shear_dirshear direction, it also contains the level of anisotropy (Rnt = a_nt k)

◆ cs_turbulence_bc_set_uninit_inlet_turb_intensity()

void cs_turbulence_bc_set_uninit_inlet_turb_intensity ( cs_lnum_t  face_id,
double  uref2,
double  t_intensity,
double  dh 
)

Set inlet boundary condition values for turbulence variables based on a diameter $ D_H $, a turbulent intensity $ I $ and the reference velocity $ U_{ref} $ for a circular duct flow with smooth wall, only if not already set.

Apart from assigning values where not already initialized, this function is similar to cs_turbulence_bc_inlet_turb_intensity.

Parameters
[in]face_idboundary face id
[in]uref2square of the reference flow velocity
[in]t_intensityturbulent intensity $ I $
[in]dhhydraulic diameter $ D_H $