9.1
general documentation
Functions | Variables
cs_boundary_conditions_set_coeffs_turb.cpp File Reference
#include "base/cs_defs.h"
#include <stdarg.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>
#include <math.h>
#include <float.h>
#include "bft/bft_printf.h"
#include "atmo/cs_air_props.h"
#include "base/cs_ale.h"
#include "atmo/cs_atmo.h"
#include "base/cs_boundary_conditions_set_coeffs.h"
#include "base/cs_dispatch.h"
#include "base/cs_field_default.h"
#include "base/cs_field_pointer.h"
#include "base/cs_internal_coupling.h"
#include "base/cs_mem.h"
#include "mesh/cs_mesh.h"
#include "base/cs_parall.h"
#include "base/cs_parameters.h"
#include "base/cs_profiling.h"
#include "base/cs_physical_constants.h"
#include "pprt/cs_physical_model.h"
#include "rayt/cs_rad_transfer.h"
#include "base/cs_thermal_model.h"
#include "base/cs_turbomachinery.h"
#include "turb/cs_turbulence_bc.h"
#include "turb/cs_turbulence_model.h"
#include "base/cs_wall_functions.h"
#include "base/cs_boundary_conditions_set_coeffs_turb.h"
+ Include dependency graph for cs_boundary_conditions_set_coeffs_turb.cpp:

Functions

static void _cs_boundary_conditions_set_coeffs_turb_scalar (cs_field_t *f_sc, int isvhb, cs_real_t byplus[], cs_real_t bdplus[], cs_real_t buk[], cs_real_t buet[], cs_real_t bcfnns[], cs_real_t bdlmo[], cs_real_t hbord[], cs_real_t theipb[], cs_real_t *tetmax, cs_real_t *tetmin, cs_real_t *tplumx, cs_real_t *tplumn)
 Compute boundary coefficients for smooth/rough walls for scalar. More...
 
static void _cs_boundary_conditions_set_coeffs_turb_vector (cs_field_t *f_v, cs_real_t byplus[], cs_real_t bdplus[], cs_real_t buk[])
 Compute boundary coefficients for smooth walls for vector. More...
 
void cs_boundary_conditions_set_coeffs_turb (int isvhb, cs_real_t velipb[][3], cs_real_t rijipb[][6], cs_real_t visvdr[], cs_real_t hbord[], cs_real_t theipb[])
 Boundary conditions for smooth walls (icodcl = CS_BC_WALL_MODELLED). More...
 

Variables

static cs_lnum_t _ntlast = -1
 
static cs_lnum_t _iaff = 0
 
static int _kbfid = -1
 
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_boundary_conditions_set_coeffs_turb_scalar()

static void _cs_boundary_conditions_set_coeffs_turb_scalar ( cs_field_t *  f_sc,
int  isvhb,
cs_real_t  byplus[],
cs_real_t  bdplus[],
cs_real_t  buk[],
cs_real_t  buet[],
cs_real_t  bcfnns[],
cs_real_t  bdlmo[],
cs_real_t  hbord[],
cs_real_t  theipb[],
cs_real_t *  tetmax,
cs_real_t *  tetmin,
cs_real_t *  tplumx,
cs_real_t *  tplumn 
)
static

Compute boundary coefficients for smooth/rough walls for scalar.

Parameters
[in]f_scscalar field
[in]isvhbid of field whose exchange coeffient should be saved at the walls, or -1.
[in]byplusdimensionless distance to the wall
[in]bdplusdimensionless shift to the wall for scalable wall functions
[in]bukdimensionless velocity
[in]buetboundary ustar value
[in]bcfnnsboundary correction factor
[in]bdlmoboundary Monin Obukhov length inverse, or nullptr
[in,out]hbordexchange coefficient at boundary
[in]theipbvalue of thermal scalar at $ \centip $ of boundary cells
[out]tetmaxmaximum local ustar value
[out]tetminminimum local ustar value
[out]tplumxmaximum local tplus value
[out]tplumnminimum local tplus value

◆ _cs_boundary_conditions_set_coeffs_turb_vector()

static void _cs_boundary_conditions_set_coeffs_turb_vector ( cs_field_t *  f_v,
cs_real_t  byplus[],
cs_real_t  bdplus[],
cs_real_t  buk[] 
)
static

Compute boundary coefficients for smooth walls for vector.

Parameters
[in]f_vvector field
[in]byplusdimensionless distance to the wall
[in]bdplusdimensionless shift to the wall for scalable wall functions
[in]bukdimensionless velocity

◆ cs_boundary_conditions_set_coeffs_turb()

void cs_boundary_conditions_set_coeffs_turb ( int  isvhb,
cs_real_t  velipb[][3],
cs_real_t  rijipb[][6],
cs_real_t  visvdr[],
cs_real_t  hbord[],
cs_real_t  theipb[] 
)

Boundary conditions for smooth walls (icodcl = CS_BC_WALL_MODELLED).

The wall functions may change the value of the diffusive flux.

The values at a boundary face $ \fib $ stored in the face center $ \centf $ of the variable $ P $ and its diffusive flux $ Q $ are written as:

\[ P_{\face} = A_P^g + B_P^g P_{\centi} \]

and

\[ Q_{\face} = A_P^f + B_P^f P_{\centi} \]

where $ P_\centi $ is the value of the variable $ P $ at the neighboring cell.

Warning:

  • For a vector field such as the velocity $ \vect{u} $ the boundary conditions may read:

    \[ \vect{u}_{\face} = \vect{A}_u^g + \tens{B}_u^g \vect{u}_{\centi} \]

    and

    \[ \vect{Q}_{\face} = \vect{A}_u^f + \tens{B}_u^f \vect{u}_{\centi} \]

    where $ \tens{B}_u^g $ and $ \tens{B}_u^f $ are 3x3 tensor matrix which coupled velocity components next to a boundary.

Please refer to the wall boundary conditions section of the theory guide for more informations, as well as the clptur section.

Parameters
[in]isvhbid of field whose exchange coeffient should be saved at the walls, or -1.
[in]velipbvalue of the velocity at $ \centip $ of boundary cells
[in]rijipbvalue of $ R_{ij} $ at $ \centip $ of boundary cells
[out]visvdrdynamic viscosity after V. Driest damping in boundary cells
[out]hbordexchange coefficient at boundary
[in]theipbvalue of thermal scalar at $ \centip $ of boundary cells

Variable Documentation

◆ _iaff

cs_lnum_t _iaff = 0
static

◆ _iv2t

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

◆ _jv2t

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

◆ _kbfid

int _kbfid = -1
static

◆ _ntlast

cs_lnum_t _ntlast = -1
static