#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_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	_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_3_t gradpr[], const cs_real_t romp[], const cs_real_t brgaus[], cs_real_t terbru, cs_real_3_t fextla)

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_3_t gradpr[], const cs_real_t romp[], const cs_real_t brgaus[], cs_real_t terbru, cs_real_3_t fextla)

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_3_t gradpr[], cs_real_t romp, cs_real_t tempf, const cs_real_t vislen[], 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_3_t gradpr[], const cs_real_t romp[], cs_real_3_t fextla, 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

◆ _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_3_t	gradpr[],
		const cs_real_t	romp[],
		cs_real_3_t *	fextla,
		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]	gradpr	pressure gradient
[in]	romp	particles associated density
[in]	fextla	external user forces (m/s2)
[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_3_t	gradpr[],
		const cs_real_t	romp[],
		const cs_real_t	brgaus[],
		cs_real_t *	terbru,
		cs_real_3_t *	fextla
	)

static

◆ _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_3_t	gradpr[],
		const cs_real_t	romp[],
		const cs_real_t	brgaus[],
		cs_real_t *	terbru,
		cs_real_3_t *	fextla
	)

static

◆ _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_3_t	gradpr[],
		cs_real_t	romp,
		cs_real_t	tempf,
		const cs_real_t	vislen[],
		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]	gradpr	pressure gradient
[in]	romp	particles associated density
[in]	tempf	temperature of the fluid (K)
[in]	vislen	FIXME
[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).

$\frac{dV}{dt} = \frac{V - PIP}{TCARAC}$

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

Functions

Function Documentation

◆ _lagdep()

◆ _lages1()

◆ _lages2()

◆ _lagesd()

◆ cs_lagr_sde()

◆ cs_lagr_sde_attr()