programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
Functions
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_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

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]dtptime step
[in]taupdynamic characteristic time
[in]tlagfluid characteristic time
[in]piilterm in integration of UP SDEs
[in]bxturbulence characteristics
[in]vagausgaussian random variables
[in]gradprpressure gradient
[in]rompparticles associated density
[in]fextlaexternal user forces (m/s2)
[in]vislenFIXME
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
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
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]dtptime step
[in]ipparticle id
[in]taupdynamic characteristic time
[in]tlagfluid characteristic time
[in]piilterm in integration of UP SDEs
[in]vagausgaussian random variables
[in]gradprpressure gradient
[in]rompparticles associated density
[in]tempftemperature of the fluid (K)
[in]vislenFIXME
[in]depintinterface location near-wall/core-flow
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_plagrangian time step
[in]taupdynamic characteristic time
[in]tlagfluid characteristic time
[in]piilterm in integration of U-P SDEs
[in]bxturbulence characteristics
[out]tsfextinfo for return coupling source terms
[in]gradprpressure gradient
[in]gradvffluid velocity gradient
[out]terbruFIXME
[in]vislenFIXME
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]attrattribute/variable
[in]tcaracvariable characteristic time
[in]pipright-hand side associated with SDE