cs_lagr_sde.cpp File Reference

#include "base/cs_defs.h"
#include "base/cs_math.h"
#include <math.h>
#include "bft/bft_mem.h"
#include "bft/bft_printf.h"
#include "mesh/cs_mesh.h"
#include "mesh/cs_mesh_quantities.h"
#include "base/cs_physical_constants.h"
#include "pprt/cs_physical_model.h"
#include "base/cs_prototypes.h"
#include "base/cs_random.h"
#include "base/cs_thermal_model.h"
#include "turb/cs_turbulence_model.h"
#include "lagr/cs_lagr.h"
#include "lagr/cs_lagr_adh.h"
#include "lagr/cs_lagr_deposition_model.h"
#include "lagr/cs_lagr_event.h"
#include "lagr/cs_lagr_roughness.h"
#include "lagr/cs_lagr_tracking.h"
#include "lagr/cs_lagr_prototypes.h"
#include "lagr/cs_lagr_sde.h"

Include dependency graph for cs_lagr_sde.cpp:

Functions
static void	_add_resuspension_event (cs_lagr_event_set_t *events, cs_lagr_particle_set_t *particles, cs_lnum_t p_id, cs_lnum_t face_id, const cs_real_t particle_velocity[3])
	Add a resulspension event.
static void	_lages1 (cs_real_t dtp, int n_phases, const cs_real_t *const *taup, const cs_real_3_t *const *tlag, const cs_real_3_t *const *piil, const cs_real_33_t *const *bx, const cs_real_3_t *const *vagaus, const cs_real_t brgaus[], const cs_real_3_t force_p[], const cs_real_3_t *const *beta, cs_real_t *terbru)
	Integration of SDEs by 1st order time scheme.
static void	_lages2 (cs_real_t dtp, const cs_real_t *const *taup_i, const cs_real_3_t *const *tlag_i, const cs_real_3_t *const *piil_i, const cs_real_33_t *const *bx_i, cs_real_t tsfext[], const cs_real_3_t *const *vagaus, const cs_real_t brgaus[], const cs_real_3_t force_p[], cs_real_t *terbru)
	Integration of SDEs by 2nd order scheme.
static void	_lagesd (cs_real_t dtp, cs_lnum_t p_id, const cs_real_t taup[], const cs_real_3_t piil[], const cs_real_3_t *const *vagaus, const cs_real_t romp[], const cs_real_3_t force_p[], cs_real_t tempf, const cs_real_t vislen[], cs_lagr_event_set_t *events, cs_real_t *depint, cs_lnum_t *nresnew)
	Deposition submodel.
static void	_lagdep (cs_real_t dtp, const cs_real_t *const *taup_i, const cs_real_3_t *const *tlag_i, const cs_real_3_t *const *piil_i, const cs_real_33_t *const *bx_i, const cs_real_3_t *const *vagaus, const cs_real_t romp[], const cs_real_3_t force_p[], const cs_real_t vislen[], cs_lnum_t *nresnew)
	Deposition submodel.
void	cs_lagr_sde (cs_real_t dt_p, const cs_real_t **taup, const cs_real_3_t **tlag, const cs_real_3_t **piil, const cs_real_33_t **bx, cs_real_t tsfext[], cs_real_3_t *force_p, cs_real_t terbru[], const cs_real_t vislen[], const cs_real_3_t **beta, cs_lnum_t *nresnew)
	Integration of particle equations of motion:
void	cs_lagr_sde_attr (cs_lagr_attribute_t attr, cs_real_t *tcarac, cs_real_t *pip)
	Integration of a stochastic differential equation (SDE) for a user particle variable (attribute).

Function Documentation

_add_resuspension_event()

void _add_resuspension_event ( cs_lagr_event_set_t * events, cs_lagr_particle_set_t * particles, cs_lnum_t p_id, cs_lnum_t face_id, const cs_real_t particle_velocity[3] )

static

Add a resulspension event.

TODO add additional info to events.

Parameters

[in]	events	pointer to events set
[in]	particles	pointer to particle set
[in]	p_id	particle id
[in]	face_id	associated face id
[in]	particle_velocity	velocity after event

_lagdep()

void _lagdep ( cs_real_t dtp, const cs_real_t *const * taup_i, const cs_real_3_t *const * tlag_i, const cs_real_3_t *const * piil_i, const cs_real_33_t *const * bx_i, const cs_real_3_t *const * vagaus, const cs_real_t romp[], const cs_real_3_t force_p[], const cs_real_t vislen[], cs_lnum_t * nresnew )

static

Deposition submodel.

Main subroutine of the submodel 1/ Calculation of the normalized wall-normal distance of the boundary-cell particles 2/ Sorting of the particles with respect to their normalized wall-normal distance

• If y^+ > depint : the standard Langevin model is applied
• If y^+ < depint : the deposition submodel is applied

Parameters

[in]	dtp	time step
[in]	taup	dynamic characteristic time
[in]	tlag	fluid characteristic time
[in]	piil	term in integration of UP SDEs
[in]	bx	turbulence characteristics
[in]	vagaus	gaussian random variables
[in]	romp	particles associated density
[in]	force_p	forces per mass unit on particles (m/s^2)
[in]	vislen	FIXME

_lages1()

void _lages1 ( cs_real_t dtp, int n_phases, const cs_real_t *const * taup, const cs_real_3_t *const * tlag, const cs_real_3_t *const * piil, const cs_real_33_t *const * bx, const cs_real_3_t *const * vagaus, const cs_real_t brgaus[], const cs_real_3_t force_p[], const cs_real_3_t *const * beta, cs_real_t * terbru )

static

Integration of SDEs by 1st order time scheme.

Parameters

[in]	dtp	time step
[in]	n_phases	number of carrier phases
[in]	taup	dynamic characteristic time
[in]	tlag	lagrangian fluid characteristic time
[in]	piil	term in integration of UP SDEs
[in]	bx	caracteristiques de la turbulence
[in]	vagaus	gaussian random variables
[in]	brgaus	gaussian random variables
[in]	force_p	forces per mass unit on particles (m/s^2)
[in]	beta	proportional to the gradient of T_lag
[out]	terbru	Diffusion coefficient accounting for Brownian (molecular) effect

_lages2()

void _lages2 ( cs_real_t dtp, const cs_real_t *const * taup_i, const cs_real_3_t *const * tlag_i, const cs_real_3_t *const * piil_i, const cs_real_33_t *const * bx_i, cs_real_t tsfext[], const cs_real_3_t *const * vagaus, const cs_real_t brgaus[], const cs_real_3_t force_p[], cs_real_t * terbru )

static

Integration of SDEs by 2nd order scheme.

When there has beed interaction with a boundary face, the velocity and velocity seen computations are forced to 1st order.

Parameters

[in]	taup	temps caracteristique dynamique
[in]	tlag	temps caracteristique fluide
[in]	piil	terme dans l'integration des eds up
[in]	bx	caracteristiques de la turbulence
[in]	tsfext	infos pour couplage retour dynamique
[in]	vagaus	variables aleatoires gaussiennes
[in]	brgaus	gaussian variable for brownian movement
[in]	romp	masse volumique des particules
[in]	force_p	forces per mass unit on particles (m/s^2)
[out]	terbru

_lagesd()

void _lagesd ( cs_real_t dtp, cs_lnum_t p_id, const cs_real_t taup[], const cs_real_3_t piil[], const cs_real_3_t *const * vagaus, const cs_real_t romp[], const cs_real_3_t force_p[], cs_real_t tempf, const cs_real_t vislen[], cs_lagr_event_set_t * events, cs_real_t * depint, cs_lnum_t * nresnew )

static

Deposition submodel.

1/ Modification of the coordinate system (global ->local) 2/ Call of subroutine lagcli 3/ Integration of the stochastic differential equations in the 2 directions different from the normal to the boundary face 4/ Modification of the coordinate system (local ->global) 5/ Update of the particle position

Parameters

[in]	dtp	time step
[in]	p_id	particle id
[in]	taup	dynamic characteristic time
[in]	tlag	fluid characteristic time
[in]	piil	term in integration of UP SDEs
[in]	vagaus	gaussian random variables
[in]	romp	particles associated density
[in]	force_p	forces per mass unit on particles (m/s^2)
[in]	tempf	temperature of the fluid (K)
[in]	vislen	FIXME
[in]	events	associated events set
[in]	depint	interface location near-wall/core-flow

cs_lagr_sde()

void cs_lagr_sde	(	cs_real_t	dt_p,
		const cs_real_t **	taup,
		const cs_real_3_t **	tlag,
		const cs_real_3_t **	piil,
		const cs_real_33_t **	bx,
		cs_real_t	tsfext[],
		cs_real_3_t *	force_p,
		cs_real_t	terbru[],
		const cs_real_t	vislen[],
		const cs_real_3_t **	beta,
		cs_lnum_t *	nresnew )

Integration of particle equations of motion:

• Standard Model : First or second order
• Deposition submodel (Guingo & Minier, 2008) if needed

Parameters

[in]	dt_p	lagrangian time step
[in]	taup	dynamic characteristic time
[in]	tlag	fluid characteristic time
[in]	piil	term in integration of U-P SDEs
[in]	bx	turbulence characteristics
[out]	tsfext	info for return coupling source terms
[out]	force_p	forces per mass unit on particles (m/s^2)
[out]	terbru	Diffusion coefficient accounting for Brownian (molecular) effect
[in]	vislen	nu/u* = y/y+
[in]	beta	proportional to the gradient of T_lag

cs_lagr_sde_attr()

void cs_lagr_sde_attr	(	cs_lagr_attribute_t	attr,
		cs_real_t *	tcarac,
		cs_real_t *	pip )

Integration of a stochastic differential equation (SDE) for a user particle variable (attribute).

When there is interaction with a boundary face, the integration degenerates to order 1 (even if the 2nd order scheme is active).

Parameters

[in]	attr	attribute/variable
[in]	tcarac	variable characteristic time
[in]	pip	right-hand side associated with SDE