#include "cs_defs.h"#include "cs_math.h"#include <math.h>#include "bft_mem.h"#include "bft_printf.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_turbulence_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:
 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[], const cs_real_3_t beta[], 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[], const cs_real_3_t beta[], 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... | |
| 
 | static | 
Add a resulspension event.
TODO add additional info to events.
| [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 | 
| 
 | 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
| [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 | 
| 
 | static | 
Integration of SDEs by 1st order time scheme.
| [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) | 
| [in] | beta | proportional to the gradient of T_lag | 
| [out] | terbru | Diffusion coefficient accounting for Brownian (molecular) effect | 
| 
 | 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.
| [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 | 
| 
 | 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
| [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 | 
| 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[], | ||
| const cs_real_3_t | beta[], | ||
| cs_lnum_t * | nresnew | ||
| ) | 
Integration of particle equations of motion:
| [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 | Diffusion coefficient accounting for Brownian (molecular) effect | 
| [in] | vislen | nu/u* = y/y+ | 
| [in] | beta | proportional to the gradient of T_lag | 
| 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} \]](form_496.png) 
When there is interaction with a boundary face, the integration degenerates to order 1 (even if the 2nd order scheme is active).
| [in] | attr | attribute/variable | 
| [in] | tcarac | variable characteristic time | 
| [in] | pip | right-hand side associated with SDE |