cs_wall_functions.cpp File Reference

#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.cpp:

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...

Function Documentation

cs_get_glob_wall_functions()

cs_wall_functions_t * cs_get_glob_wall_functions

(

void

)

cs_immersed_boundary_wall_functions()

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.

This is similar to clptur.

Parameters

[in]	f_id	field id of the variable
[out]	st_exp	explicit source term
[out]	st_imp	implicit part of the source term

cs_wall_functions_scalar()

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:

Parameters

[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

cs_wall_functions_velocity()

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 / .

Parameters

[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