#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"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, 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. | |
| 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. | |
| 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. | |
| 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. | |
| 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: | |
| 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). | |
|
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_504.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 |