#include "cs_config.h"
#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_wall_functions.h"
Functions | |
void | wallfunctions (const int *const iwallf, const cs_lnum_t *const ifac, const cs_real_t *const l_visc, const cs_real_t *const t_visc, const cs_real_t *const vel, const cs_real_t *const y, const cs_real_t *const rough_d, const cs_real_t *const rnnb, const cs_real_t *const kinetic_en, int *iuntur, cs_lnum_t *nsubla, cs_lnum_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) |
void | hturbp (const int *const iwalfs, const cs_real_t *const l_visc, const cs_real_t *const prl, const cs_real_t *const prt, const cs_real_t *const rough_t, const cs_real_t *const uk, const cs_real_t *const yplus, const cs_real_t *const dplus, cs_real_t *htur, cs_real_t *yplim) |
cs_wall_functions_t * | cs_get_glob_wall_functions (void) |
void | cs_wall_functions_velocity (cs_wall_f_type_t iwallf, cs_lnum_t ifac, 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_lnum_t *nsubla, cs_lnum_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 ![]() ![]() | |
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... | |
Wall functions descriptor and computation.
cs_wall_functions_t* cs_get_glob_wall_functions | ( | void | ) |
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:
[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_lnum_t | ifac, | ||
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_lnum_t * | nsubla, | ||
cs_lnum_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] | ifac | face number |
[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 | ![]() |
[out] | dplus | dimensionless shift to the wall for scalable wall functions |
void hturbp | ( | const int *const | iwalfs, |
const cs_real_t *const | l_visc, | ||
const cs_real_t *const | prl, | ||
const cs_real_t *const | prt, | ||
const cs_real_t *const | rough_t, | ||
const cs_real_t *const | uk, | ||
const cs_real_t *const | yplus, | ||
const cs_real_t *const | dplus, | ||
cs_real_t * | htur, | ||
cs_real_t * | yplim | ||
) |
void wallfunctions | ( | const int *const | iwallf, |
const cs_lnum_t *const | ifac, | ||
const cs_real_t *const | l_visc, | ||
const cs_real_t *const | t_visc, | ||
const cs_real_t *const | vel, | ||
const cs_real_t *const | y, | ||
const cs_real_t *const | rough_d, | ||
const cs_real_t *const | rnnb, | ||
const cs_real_t *const | kinetic_en, | ||
int * | iuntur, | ||
cs_lnum_t * | nsubla, | ||
cs_lnum_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 | ||
) |