#include "cs_defs.h"#include <assert.h>#include <errno.h>#include <stdio.h>#include <stdarg.h>#include <string.h>#include <math.h>#include <stdlib.h>#include <float.h>#include <mpi.h>#include "bft_error.h"#include "bft_mem.h"#include "bft_printf.h"#include "cs_log.h"#include "cs_mesh.h"#include "cs_mesh_quantities.h"#include "cs_turbulence_model.h"#include "cs_domain.h"#include "cs_field.h"#include "cs_field_operator.h"#include "cs_field_pointer.h"#include "cs_field_default.h"#include "cs_wall_functions.h" Include dependency graph for cs_wall_functions.c:
 Include dependency graph for cs_wall_functions.c:| Functions | |
| cs_wall_functions_t * | cs_get_glob_wall_functions (void) | 
| void | cs_wall_functions_velocity (cs_wall_f_type_t iwallf, cs_real_t l_visc, cs_real_t t_visc, cs_real_t vel, cs_real_t y, cs_real_t rough_d, cs_real_t rnnb, cs_real_t kinetic_en, int *iuntur, cs_gnum_t *nsubla, cs_gnum_t *nlogla, cs_real_t *ustar, cs_real_t *uk, cs_real_t *yplus, cs_real_t *ypup, cs_real_t *cofimp, cs_real_t *dplus) | 
| Compute the friction velocity and  /  .  More... | |
| void | cs_wall_functions_scalar (cs_wall_f_s_type_t iwalfs, cs_real_t l_visc, cs_real_t prl, cs_real_t prt, cs_real_t rough_t, cs_real_t uk, cs_real_t yplus, cs_real_t dplus, cs_real_t *htur, cs_real_t *yplim) | 
| Compute the correction of the exchange coefficient between the fluid and the wall for a turbulent flow.  More... | |
| void | cs_immersed_boundary_wall_functions (int f_id, cs_real_t *st_exp, cs_real_t *st_imp) | 
| Compute boundary contributions for all immersed boundaries.  More... | |
Wall functions descriptor and computation.
| cs_wall_functions_t* cs_get_glob_wall_functions | ( | void | ) | 
Compute boundary contributions for all immersed boundaries.
This is similar to clptur.
| [in] | f_id | field id of the variable | 
| [out] | st_exp | explicit source term | 
| [out] | st_imp | implicit part of the source term | 
| void cs_wall_functions_scalar | ( | cs_wall_f_s_type_t | iwalfs, | 
| cs_real_t | l_visc, | ||
| cs_real_t | prl, | ||
| cs_real_t | prt, | ||
| cs_real_t | rough_t, | ||
| cs_real_t | uk, | ||
| cs_real_t | yplus, | ||
| cs_real_t | dplus, | ||
| cs_real_t * | htur, | ||
| cs_real_t * | yplim | ||
| ) | 
Compute the correction of the exchange coefficient between the fluid and the wall for a turbulent flow.
This is function of the dimensionless distance to the wall  .
.
Then the return coefficient reads:
![\[ h_{tur} = Pr \dfrac{y^+}{T^+} \]](form_300.png) 
| [in] | iwalfs | type of wall functions for scalar | 
| [in] | l_visc | kinematic viscosity | 
| [in] | prl | laminar Prandtl number | 
| [in] | prt | turbulent Prandtl number | 
| [in] | rough_t | scalar roughness lenghth scale | 
| [in] | uk | velocity scale based on TKE | 
| [in] | yplus | dimensionless distance to the wall | 
| [in] | dplus | dimensionless distance for scalable wall functions | 
| [out] | htur | corrected exchange coefficient | 
| [out] | yplim | value of the limit for   | 
| void cs_wall_functions_velocity | ( | cs_wall_f_type_t | iwallf, | 
| cs_real_t | l_visc, | ||
| cs_real_t | t_visc, | ||
| cs_real_t | vel, | ||
| cs_real_t | y, | ||
| cs_real_t | rough_d, | ||
| cs_real_t | rnnb, | ||
| cs_real_t | kinetic_en, | ||
| int * | iuntur, | ||
| cs_gnum_t * | nsubla, | ||
| cs_gnum_t * | nlogla, | ||
| cs_real_t * | ustar, | ||
| cs_real_t * | uk, | ||
| cs_real_t * | yplus, | ||
| cs_real_t * | ypup, | ||
| cs_real_t * | cofimp, | ||
| cs_real_t * | dplus | ||
| ) | 
Compute the friction velocity and  /
 /  .
. 
| [in] | iwallf | wall function type | 
| [in] | l_visc | kinematic viscosity | 
| [in] | t_visc | turbulent kinematic viscosity | 
| [in] | vel | wall projected cell center velocity | 
| [in] | y | wall distance | 
| [in] | rough_d | roughness length scale (not sand grain roughness) | 
| [in] | rnnb |   | 
| [in] | kinetic_en | turbulent kinetic energy (cell center) | 
| [in,out] | iuntur | indicator: 0 in the viscous sublayer | 
| [in,out] | nsubla | counter of cell in the viscous sublayer | 
| [in,out] | nlogla | counter of cell in the log-layer | 
| [out] | ustar | friction velocity | 
| [out] | uk | friction velocity | 
| [out] | yplus | dimensionless distance to the wall | 
| [out] | ypup | yplus projected vel ratio | 
| [out] | cofimp |  to ensure a good turbulence production | 
| [out] | dplus | dimensionless shift to the wall for scalable wall functions |