cs_lagr_sde.c File Reference

#include "cs_defs.h"
#include "cs_math.h"
#include <math.h>
#include "bft_mem.h"
#include "cs_mesh.h"
#include "cs_mesh_quantities.h"
#include "cs_physical_constants.h"
#include "cs_physical_model.h"
#include "cs_prototypes.h"
#include "cs_random.h"
#include "cs_thermal_model.h"
#include "cs_lagr.h"
#include "cs_lagr_adh.h"
#include "cs_lagr_deposition_model.h"
#include "cs_lagr_event.h"
#include "cs_lagr_roughness.h"
#include "cs_lagr_tracking.h"
#include "cs_lagr_prototypes.h"
#include "cs_lagr_sde.h"

Include dependency graph for cs_lagr_sde.c:

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. More...
static void	_lages1 (cs_real_t dtp, const cs_real_t taup[], const cs_real_3_t tlag[], const cs_real_3_t piil[], const cs_real_33_t bx[], const cs_real_33_t vagaus[], const cs_real_t brgaus[], const cs_real_3_t force_p[], cs_real_t *terbru)
	Integration of SDEs by 1st order time scheme. More...
static void	_lages2 (cs_real_t dtp, 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[], const cs_real_33_t vagaus[], const cs_real_t brgaus[], const cs_real_3_t force_p[], cs_real_t *terbru)
	Integration of SDEs by 2nd order scheme. More...
static void	_lagesd (cs_real_t dtp, cs_lnum_t ip, const cs_real_t taup[], const cs_real_3_t piil[], const cs_real_33_t 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. More...
static void	_lagdep (cs_real_t dtp, const cs_real_t taup[], const cs_real_3_t tlag[], const cs_real_3_t piil[], const cs_real_33_t bx[], const cs_real_33_t vagaus[], const cs_real_t romp[], const cs_real_3_t force_p[], const cs_real_t vislen[], cs_lnum_t *nresnew)
	Deposition submodel. More...
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[], const cs_real_3_t gradpr[], const cs_real_33_t gradvf[], cs_real_t terbru[], const cs_real_t vislen[], cs_lnum_t *nresnew)
	Integration of particle equations of motion: More...
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). More...

Function Documentation

_add_resuspension_event()

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]  )

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()

static void _lagdep ( cs_real_t  dtp, const cs_real_t  taup[], const cs_real_3_t  tlag[], const cs_real_3_t  piil[], const cs_real_33_t  bx[], const cs_real_33_t  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	taup times forces on particles (m/s)
[in]	vislen	FIXME

_lages1()

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

static

Integration of SDEs by 1st order time scheme.

Parameters

[in]	dtp	time step
[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	taup times forces on particles (m/s)
[out]	terbru

_lages2()

static void _lages2 ( cs_real_t  dtp, 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[], const cs_real_33_t  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	taup times forces on particles (m/s)
[out]	terbru

_lagesd()

static void _lagesd ( cs_real_t  dtp, cs_lnum_t  ip, const cs_real_t  taup[], const cs_real_3_t  piil[], const cs_real_33_t  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]	ip	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	taup times forces on particles (m/s)
[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[],
		const cs_real_3_t	gradpr[],
		const cs_real_33_t	gradvf[],
		cs_real_t	terbru[],
		const cs_real_t	vislen[],
		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
[in]	gradpr	pressure gradient
[in]	gradvf	fluid velocity gradient
[out]	terbru	FIXME
[in]	vislen	FIXME

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