programmer's documentation
cs_lagr.h
Go to the documentation of this file.
1 #ifndef __CS_LAGR_H__
2 #define __CS_LAGR_H__
3 
4 /*============================================================================
5  * Functions and types for the Lagrangian module
6  *============================================================================*/
7 
8 /*
9  This file is part of Code_Saturne, a general-purpose CFD tool.
10 
11  Copyright (C) 1998-2016 EDF S.A.
12 
13  This program is free software; you can redistribute it and/or modify it under
14  the terms of the GNU General Public License as published by the Free Software
15  Foundation; either version 2 of the License, or (at your option) any later
16  version.
17 
18  This program is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
20  FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
21  details.
22 
23  You should have received a copy of the GNU General Public License along with
24  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
25  Street, Fifth Floor, Boston, MA 02110-1301, USA.
26 */
27 
28 /*----------------------------------------------------------------------------*/
29 
30 #include "cs_defs.h"
31 
32 #include "assert.h"
33 #include "cs_base.h"
34 #include "cs_field.h"
35 
36 /*----------------------------------------------------------------------------*/
37 
39 
40 /*============================================================================
41  * Type definitions
42  *============================================================================*/
43 
46 typedef enum {
47 
61 
63 
66 typedef enum {
73 
75 /*---------------------*/
76 
77 typedef struct {
78 
79  int nusbrd;
82  int nflagm;
83  int ndlaim;
85  int nvgaus;
86  int nbrgau;
87 
88  int ncharm2;
89  int nlayer;
92 
94 /*---------------------*/
95 
96 typedef struct {
97 
98  int ntersl;
99  int nvisbr;
101 } cs_lagr_dim_t;
102 
104 /*-----------------------------------------------------*/
105 
106 typedef struct {
107 
120  int iilagr;
121 
122  /* indicates the steady (=1) or unsteady (=0) state of the
123  continuous phase flow
124  in particular, \ref isttio = 1 is needed in order to:
125  calculate steady statistics in the volume or at the boundaries
126  (starting respectively from the iterations \ref nstist)
127  and calculate time-averaged two-way coupling source terms (from the
128  time step \ref nstits).
129  Useful if \ref iilagr=1 or \ref iilagr=2 (if \ref iilagr=3,
130  then \ref isttio=1 automatically) */
131  int isttio;
132 
138  int isuila;
139 
141  int t_order;
142 
149  int modcpl;
150 
154  int idirla;
155 
160  int idistu;
161 
166  int idiffl;
167 
173  int ilapoi;
174 
191 
194 
196 
198 /*-----------------------------------------------------------*/
199 
200 typedef struct {
201 
212  int physical_model; /* FIXME: => enum: CS_LAGR_PHYS_STD,
213  CS_LAGR_PHYS_COAL,
214  CS_LAGR_PHYS_HEAT... */
216 
218  int dlvo;
219 
223 
227 
228  /* - 0: no clogging model
229  - 1: clogging model */
230  int clogging;
231 
233  int fouling;
234 
237 
239 
240 /* ========================================================================== */
241 
242 typedef struct {
243 
247 
251 
254 
257 
260 
263 
266 
269 
272 
275 
278 
281 
284 
287 
290 
292 
293 /* ========================================================================== */
294 
295 typedef struct {
296 
297  /* activation (=1) or not (=0) of an evolution equation on the particle
298  temperature (in degrees Celsius).
299  Useful if \ref physical_model=1 and if there is a thermal scalar associated with
300  the continuous phase
301  */
302  int itpvar;
303 
304  /* activation (=1) or not (=0) of an evolution equation on the particle
305  diameter. Useful if \ref physical_model = 1
306  */
307  int idpvar;
308 
309  /* activation (=1) or not (=0) of an evolution equation on the particle mass
310  Useful if \ref physical_model = 1
311  */
312  int impvar;
313 
314  /* initialization temperature (in degree Celsius) for the particles already
315  present in the calculation domain when an evolution equation on
316  the particle temperature is activated during a calculation (\ref physical_model =
317  1 and \ref itpvar = 1).
318  Useful if \ref isuila = 1 and \ref itpvar = 0 in the previous calculation
319  */
321 
322  /* initialization value for the specific heat (\f$ J.kg^{-1}.K^{-1} \f$)
323  of the particles already present
324  in the calculation domain when an evolution equation
325  on the particle temperature is activated during a calculation
326  (\ref physical_model = 1 and \ref itpvar = 1).
327  Useful if \ref isuila = 1 and \ref itpvar = 0 in the previous calculation
328  */
330 
332 
333 /* ========================================================================== */
334 
335 typedef struct {
336 
337  /* - 0: no resuspension model
338  - 1: resuspension model */
339  int ireent;
340 
341  /* - 0: no head losses calculation for influence of the deposit on the flow
342  - 1: head losses calculation for influence of the deposit on the flow */
343  int iflow;
344 
345  /* Parameters of the particle resuspension model*/
351 
353 
354 /* ========================================================================== */
355 
356 typedef struct {
357 
358  /* number of particle classes*/
359  int nbrclas;
360  /* diameter of particles formed by precipitation*/
362  /* density of particles formed by precipitation*/
364  /* number of precipitated particles */
365  int *nbprec;
366  /* */
368  /* number of precipitated particles */
370 
372 
373 /* ========================================================================== */
374 
375 typedef struct {
376 
377  /* Parameter of the particle clogging model */
381 
383 
384 /* ========================================================================== */
385 
386 typedef struct {
387 
388  /* current step id (for 2nd order scheme) */
389  int nor;
390 
391  /* duration of a Lagrangian iteration */
393 
394  /* physical time of the Lagrangian simulation */
396 
398 
399 /* ========================================================================== */
400 
401 typedef struct {
402 
405 
409 
416 
421 
426 
431 
436 
439 
441  int cluster;
442 
445 
447  cs_real_t velocity[3];
448 
451 
454 
457 
460 
463 
466 
469 
472 
475 
478 
481 
484 
487 
490 
493 
495 
496 /* ========================================================================== */
497 
498 /* ========================================================================== */
499 
500 typedef struct {
501 
505  int ltsdyn;
506 
509  int ltsmas;
510 
511  /* if \ref physical_model = 1 and \ref itpvar = 1, \ref ltsthe
512  activates (=1) or not (=0) the two-way coupling on temperature.
513  if \ref physical_model = 2, \ref ltsthe activates (=1) or not (=0) the
514  two-way coupling on the eulerian variables related to pulverised
515  coal combustion.
516  Useful if \ref iilagr = 2 */
517  int ltsthe;
518 
520  int itsvx;
521 
523  int itsvy;
524 
526  int itsvz;
527 
530  int itsli;
531 
532  /* explicit source term for the turbulent dissipation and the
533  turbulent energy if the \f$k-\varepsilon\f$ turbulence model is used
534  for the continuous phase */
535  int itske;
536 
540  int itsr11;
541 
545  int itsr12;
546 
550  int itsr13;
551 
555  int itsr22;
556 
560  int itsr23;
561 
565  int itsr33;
566 
569  int itste;
570 
573  int itsti;
574 
576  int itsmas;
577 
578  /* source term for the light volatile matters */
579 //TODO
580  int *itsmv1;//ncharm2
581 
582  /* source term for the heavy volatile matters */
583 //TODO
584  int *itsmv2;//ncharm2
585 
587  int itsco;
588 
590  int itsfp4;
591 
604  int nstits;
605 
607  int npts;
608 
611  int ntxerr;
612 
615 
618 
621 
623 
624 /* ========================================================================== */
625 
626 /* Structures useful to deal with boundary conditions
627  For USLABO => _boundary_track_treatment */
628 
629 typedef struct {
630 
631  int n_b_zones; /* NFRLAG */
633 
634  cs_lnum_t *b_zone_id; /* ILFLAG */
635  int *b_zone_classes; /* IUSNCL */
636  int *b_zone_natures; /* IUSCLB */
637 
638  int *b_face_zone_id; /* IFRLAG */
639 
641 
642  cs_real_t *particle_flow_rate; /* DEBLAG -> post-processing use */
643 
645 
646 /* Structures useful to deal with iternal conditions */
647 
648 typedef struct {
649 
651 
653 
654 /* ========================================================================== */
655 
656 typedef struct {
657 
658  /* activates (=1) or not (=0) the option of coal particle fouling.
659  It then is necessary to specify the domain boundaries
660  on which fouling may take place. Useful if \ref physical_model = 2*/
661  int iencra;
662 
663  /* encrustation data*/
664  int npencr;
665  // TODO cf particles->n_part_fou in cs_lagr_tracking.c
666 
667  /* encrustation data*/
668 //TODO
669  cs_real_t *enc1;//ncharm2
670  /* encrustation data*/
671 //TODO
672  cs_real_t *enc2;//ncharm2
673 
674  /* limit temperature (in degree Celsius) below which the coal particles do
675  not cause any fouling (if the fouling model is activated).
676  Useful if \ref physical_model = 2 and \ref iencra = 1*/
677 //TODO
678  cs_real_t *tprenc;//ncharm2
679 
680  /* ash critical viscosity in \f$ kg.m^{-1}.s^{-1} \f$, in the fouling model
681  cf J.D. Watt et T. Fereday (J.Inst.Fuel, Vol.42-p99).
682  Useful if \ref physical_model = 2 and \ref iencra = 1*/
683 //TODO
684  cs_real_t *visref;//ncharm2
685 
686  /* encrustation data */
688 
690 
691 /* ========================================================================== */
692 
693 typedef struct {
694 
697 
701 
704 
707 
710 
713 
716 
718 
719 /* ========================================================================== */
720 
721 typedef struct {
722 
723  /* brownnian motion activation */
724  int lamvbr;
725 
727 
728 /* ========================================================================== */
729 
730 typedef struct {
731 
735  int nusbor;
736 
743  int npstf;
744 
752  int npstft;
753 
760  int inbrbd;
761 
767  int iflmbd;
768 
773  int iangbd;
774 
779  int ivitbd;
780 
785  int iclgst;
786 
788  int iencnbbd;
789 
791  int iencmabd;
792 
794  int iencdibd;
795 
797  int iencckbd;
798 
800  int inbr;
801 
803  int iflm;
804 
806  int iang;
807 
809  int ivit;
810 
812  int ires;
813 
815  int iflres;
816 
818  int iencnb;
819 
821  int iencma;
822 
824  int iencdi;
825 
827  int iencck;
828 
830  int *iusb;
831 
832  /* the recordings in \ref bound_stat at every particle/boundary interaction are
833  cumulated values (possibly reset to zero at every iteration in the
834  unsteady case). They must therefore be divided by a quantity to
835  get boundary statistics. The user can choose between two average types:
836  - = 0: no average is applied to the recorded cumulated values.
837  - = 1: a time-average is calculated. The cumulated value
838  is divided by the physical duration in the case of steady
839  averages (\ref isttio=1). The cumulated value is divided by the
840  value of the last time step in the case of unsteady averages
841  (\ref isttio=0), and also in the case of steady averages while the
842  absolute iteration number is inferior to \ref nstist.
843  - = 2: a particulate average is calculated. The cumulated value is divided
844  by the number of particle/boundary interactions (in terms of
845  statistical weight) recorded in \ref bound_stat "bound_stat"(nfabor,inbr).
846  This average can only be calculated when \ref inbrbd=1.
847  Only the cumulated value is recorded in the restart file.
848  Useful if \ref iensi3=1 */
849  int *imoybr;
850 
852  int inclg;
853 
855  int inclgt;
856 
858  int iclogt;
859 
861  int iclogh;
862 
864  int iscovc;
865 
866  /* id for mean of particle deposition height */
867  int ihdepm;
868 
869  /* id for variance of particle deposition height */
870  int ihdepv;
871 
872  /* id for mean diameter of deposited particles */
873  int ihdiam;
874 
875  /* id for sum of deposited particle diameters */
876  int ihsum;
877 
886 
894  char **nombrd;
895 
897 
898 /* ========================================================================== */
899 
900 typedef struct {
901 
902  /* \anchor iensi3
903  activation (=1) or not (=0) of the recording of the particle/boundary
904  interactions in \ref bound_stat, and of the calculation of the
905  statistics at the corresponding boundaries.
906  By default, the statistics are unsteady (reset to zero at every
907  time step). They may be steady if \ref isttio=1 (i.e.
908  calculation of a cumulated value over time, and then calculation of an
909  average over time or over the number of interactions with the boundary).*/
910  int iensi3;
911 
912  /* associates (=1) or not (=0) the variable "velocity of the locally
913  undisturbed fluid flow field" with the output of particles or
914  trajectories. */
915  int ivisv1;
916 
917  /* associates (=1) or not (=0) the variable "particle velocity"
918  with the output of particles or trajectories. */
919  int ivisv2;
920 
921  /* associates (=1) or not (=0) the variable "residence time"
922  with the output of particles or trajectories. */
923  int ivistp;
924 
925  /* associates (=1) or not (=0) the variable "particle diameter"
926  with the output of particles or trajectories. */
927  int ivisdm;
928 
929  /* associates (=1) or not (=0) the variable "particle temperature"
930  with the output of particles or trajectories. */
931  int iviste;
932 
933  /* associates (=1) or not (=0) the variable "particle mass"
934  with the output of particles or trajectories. */
935  int ivismp;
936 
937  /* associates (=1) or not (=0) the variable "shrinking core diameter of
938  the coal particles" with the output of particles or trajectories.
939  useful only if \ref physical_model = 2 */
940  int ivisdk;
941 
942  /* associates (=1) or not (=0) the variable "mass of reactive coal of the
943  coal particles" with the output of particles or trajectories.
944  useful only if \ref physical_model = 2 */
945  int ivisch;
946 
947  /* associates (=1) or not (=0) the variable "mass of coal of the
948  coal particles" with the output of particles or trajectories.
949  useful only if \ref physical_model = 2 */
950  int ivisck;
951 
953  int iviswat;
954 
956 
957 /* ========================================================================== */
958 
959 typedef struct {
960 
961  /* Turbulence model */
962  int iturb;
963  int itytur;
964 
965  /* cpincl */
966  int ncharb;
967 
968  /* ppppar */
969  int ncharm;
970  int nozppm;
971 
972  /* radiation */
973  int iirayo;
974 
975  /* icp */
976  int icp;
977 
978  /* diftl0 */
980 
981  /* cmu */
983 
984  /* visls0 */
986 
987  /*****************
988  * Useful fields *
989  *****************/
990 
991  /* wall ustar */
993 
994  /* Fluid density */
996 
997  /* Fluid pressure */
999 
1000  /* Fluid temparature */
1004 
1005  /* Fluid velocity */
1007 
1008  /* Fluid viscosity */
1010 
1011  /* Fluid viscosity */
1013 
1014  /* Fluid specific heat capacity */
1016 
1017  /* Radiat. */
1019 
1020  /* Combustion */
1024 
1025  /* Turbulence */
1026  /* Turbulent intensity */
1028 
1029  /* Turbulent dissipation */
1031 
1032  /* Omega from k-omega SST model*/
1034 
1035  /* Reynolds stress component Rxx */
1037  /* Reynolds stress component Ryy */
1039  /* Reynolds stress component Rzz */
1041 
1043 
1044 /* external data relative to coal combustion */
1045 
1046 typedef struct {
1047 
1048  int ih2o; // cpincl
1049  int io2; // cpincl
1050  int ico; // cpincl
1051 
1052  int iatc; // ppthch
1053  cs_real_t prefth; // ppthch
1054  cs_real_t trefth; // ppthch
1055 
1056  int natom; // = 5;
1057  cs_real_t *wmolat; // dim = natom
1058 
1059  int ngazem; // = 20;
1060  cs_real_t *wmole; // ngazem
1061  int *iym1;
1062 
1063  int ncharm; // cpincl
1064  cs_real_t *a1ch; // ncharm
1078 
1080 
1081 /*============================================================================
1082  * Global variables
1083  *============================================================================*/
1084 
1088 
1092 
1095 
1098 
1101 
1102 /* Lagrangian log output every frequency_n time steps */
1103 
1104 extern int cs_glob_lagr_log_frequency_n;
1105 
1106 /* Statisics on borders*/
1107 extern cs_real_t *bound_stat;
1108 
1109 extern int cs_glob_lagr_nzone_max;
1110 extern int cs_glob_lagr_nclass_max;
1111 
1122 
1125 
1129 
1130 /* Unit normals and offsets of boundary faces */
1132 
1133 /* Projection matrices for global to local coordinates on boundary faces */
1135 
1136 /*============================================================================
1137  * Public function prototypes
1138  *============================================================================*/
1139 
1140 /*----------------------------------------------------------------------------*/
1150 /* ----------------------------------------------------------------------------*/
1151 
1153 cs_lagr_get_zone_class_data(int iclass,
1154  int izone);
1155 
1156 /*----------------------------------------------------------------------------*/
1167 /*----------------------------------------------------------------------------*/
1168 
1169 void
1171  int izone,
1172  int number,
1173  int freq,
1174  int stat);
1175 
1176 /*----------------------------------------------------------------------------*/
1186 /*----------------------------------------------------------------------------*/
1187 
1188 void
1190  int izone,
1191  int profile,
1192  cs_real_t *temp,
1194 
1195 /*----------------------------------------------------------------------------*/
1203 /*----------------------------------------------------------------------------*/
1204 
1205 void
1206 cs_lagr_set_zone_class_cp(int iclass,
1207  int izone,
1208  cs_real_t cp);
1209 
1210 /*----------------------------------------------------------------------------*/
1226 /*----------------------------------------------------------------------------*/
1227 
1228 void
1229 cs_lagr_set_zone_class_coal(int iclass,
1230  int izone,
1231  int profile,
1232  int number,
1233  cs_real_t *temp,
1234  cs_real_t *coal_mf,
1235  cs_real_t *coke_mf,
1236  cs_real_t *coke_density,
1237  cs_real_t water_mf,
1238  cs_real_t shrink_diam,
1239  cs_real_t init_diam);
1240 
1241 /*----------------------------------------------------------------------------*/
1252 /*----------------------------------------------------------------------------*/
1253 
1254 void
1255 cs_lagr_set_zone_class_stat(int iclass,
1256  int izone,
1257  int profile,
1258  cs_real_t weight,
1259  cs_real_t flow);
1260 
1261 /*----------------------------------------------------------------------------*/
1271 /*----------------------------------------------------------------------------*/
1272 
1273 void
1274 cs_lagr_set_zone_class_diam(int iclass,
1275  int izone,
1276  int profile,
1277  cs_real_t diam,
1278  cs_real_t diam_dev);
1279 
1280 /*----------------------------------------------------------------------------*/
1288 /*----------------------------------------------------------------------------*/
1289 
1290 void
1292  int izone,
1293  cs_real_t density);
1294 
1295 /*----------------------------------------------------------------------------*/
1303 /*----------------------------------------------------------------------------*/
1304 
1305 void
1307  int izone,
1308  cs_real_t foul_index);
1309 
1310 /*----------------------------------------------------------------------------*/
1319 /*----------------------------------------------------------------------------*/
1320 
1321 void
1323  int izone,
1324  int profile,
1325  cs_real_t velocity[]);
1326 
1327 /*----------------------------------------------------------------------------
1328  * \brief Initialize Lagrangian module parameters for a given set of data
1329  *
1330  *
1331  *----------------------------------------------------------------------------*/
1332 
1334 cs_lagr_init_zone_class_new(int iclass,
1335  int izone);
1336 
1337 /*----------------------------------------------------------------------------*/
1344 /*----------------------------------------------------------------------------*/
1345 
1348 
1349 /*----------------------------------------------------------------------------*/
1358 /*----------------------------------------------------------------------------*/
1359 
1362 
1363 /*----------------------------------------------------------------------------
1364  * Provide access to cs_lagr_particle_counter_t
1365  *
1366  * needed to initialize structure with GUI
1367  *----------------------------------------------------------------------------*/
1368 
1371 
1372 /*----------------------------------------------------------------------------
1373  * Provide access to cs_lagr_reentrained_model_t
1374  *
1375  * needed to initialize structure with GUI
1376  *----------------------------------------------------------------------------*/
1377 
1380 
1381 /*----------------------------------------------------------------------------
1382  * Provide access to cs_lagr_precipitation_model_t
1383  *
1384  * needed to initialize structure with GUI
1385  *----------------------------------------------------------------------------*/
1386 
1389 
1390 /*----------------------------------------------------------------------------
1391  * Provide access to cs_lagr_clogging_model_t
1392  *
1393  * needed to initialize structure with GUI
1394  *----------------------------------------------------------------------------*/
1395 
1398 
1399 /*----------------------------------------------------------------------------
1400  * Provide access to cs_lagr_time_step_t
1401  *
1402  * needed to initialize structure with GUI
1403  *----------------------------------------------------------------------------*/
1404 
1406 cs_get_lagr_time_step(void);
1407 
1408 /*----------------------------------------------------------------------------
1409  * Provide access to cs_lagr_source_terms_t
1410  *
1411  * needed to initialize structure with GUI
1412  *----------------------------------------------------------------------------*/
1413 
1416 
1417 /*----------------------------------------------------------------------------
1418  * Provide access to cs_lagr_encrustation_t
1419  *
1420  * needed to initialize structure with GUI
1421  *----------------------------------------------------------------------------*/
1422 
1425 
1426 /*----------------------------------------------------------------------------
1427  * Provide access to cs_lagr_physico_chemical_t
1428  *
1429  * needed to initialize structure with GUI
1430  *----------------------------------------------------------------------------*/
1431 
1434 
1435 /*----------------------------------------------------------------------------
1436  * Provide access to cs_lagr_brownian_t
1437  *
1438  * needed to initialize structure with GUI
1439  *----------------------------------------------------------------------------*/
1440 
1442 cs_get_lagr_brownian(void);
1443 
1444 /*----------------------------------------------------------------------------*/
1451 /*----------------------------------------------------------------------------*/
1452 
1455 
1456 /*----------------------------------------------------------------------------*/
1463 /*----------------------------------------------------------------------------*/
1464 
1467 
1468 /*----------------------------------------------------------------------------
1469  * Destroy finalize the global cs_lagr_bdy_condition_t structure.
1470  *----------------------------------------------------------------------------*/
1471 
1472 void
1474 
1475 /*----------------------------------------------------------------------------
1476  * Destroy finalize the global cs_lagr_internal_condition_t structure.
1477  *----------------------------------------------------------------------------*/
1478 
1479 void
1481 
1482 /*----------------------------------------------------------------------------
1483  * Provide access to cs_lagr_boundary_interactions_t
1484  *
1485  * needed to initialize structure with GUI
1486  *----------------------------------------------------------------------------*/
1487 
1490 
1491 /*----------------------------------------------------------------------------
1492  * Provide access to cs_lagr_visualization_t
1493  *
1494  * needed to initialize structure with GUI
1495  *----------------------------------------------------------------------------*/
1496 
1499 
1500 /*----------------------------------------------------------------------------
1501  * Provide access to cs_lagr_extra_module_t
1502  *
1503  *----------------------------------------------------------------------------*/
1504 
1507 
1508 /*--------------------------------------------------------------------
1509  * Execute one time step of the Lagrangian model.
1510  *
1511  * This is the main function for that model.
1512  *
1513  * parameters:
1514  * itypfb <-- boundary face types
1515  * dt <-- time step (per cell)
1516  *-------------------------------------------------------------------- */
1517 
1518 void
1519 cs_lagr_solve_time_step(const int itypfb[],
1520  const cs_real_t *dt);
1521 
1522 /*----------------------------------------------------------------------------
1523  * Return pointers to lagrangian arrays
1524  *
1525  * This function is intended for use by Fortran wrappers.
1526  *
1527  * parameters:
1528  * dim_bound_stat --> dimensions for bound_stat pointer
1529  * p_bound_stat --> bound_stat pointer
1530  *----------------------------------------------------------------------------*/
1531 
1532 void
1533 cs_lagr_init_c_arrays(int dim_cs_glob_lagr_source_terms[2],
1534  cs_real_t **p_cs_glob_lagr_source_terms);
1535 
1536 
1537 /*----------------------------------------------------------------------------
1538  * Free lagrangian arrays
1539  *
1540  * This function is intended for use by Fortran wrappers.
1541  *----------------------------------------------------------------------------*/
1542 void
1543 cs_lagr_finalize(void);
1544 
1545 /*----------------------------------------------------------------------------*/
1546 
1548 
1549 #endif /* __CS_LAGR_H__ */
cs_real_t * bound_stat
const cs_lagr_particle_counter_t * cs_glob_lagr_particle_counter
cs_real_t cmu
Definition: cs_lagr.h:982
int itsli
Definition: cs_lagr.h:530
cs_gnum_t n_g_exit
Definition: cs_lagr.h:259
Definition: cs_lagr.h:77
cs_real_t csthpp
Definition: cs_lagr.h:380
cs_real_t emissivity
Definition: cs_lagr.h:474
unsigned long cs_gnum_t
global mesh entity number
Definition: cs_defs.h:280
Definition: cs_lagr.h:295
int ilapoi
Definition: cs_lagr.h:173
cs_real_t * visref
Definition: cs_lagr.h:684
int idiffl
Definition: cs_lagr.h:166
cs_real_t * enc2
Definition: cs_lagr.h:672
void cs_lagr_set_zone_class_coal(int iclass, int izone, int profile, int number, cs_real_t *temp, cs_real_t *coal_mf, cs_real_t *coke_mf, cs_real_t *coke_density, cs_real_t water_mf, cs_real_t shrink_diam, cs_real_t init_diam)
Set coal parameters for a given class and boundary zone.
Definition: cs_lagr.c:1147
int ivisv1
Definition: cs_lagr.h:915
cs_lagr_bdy_condition_t * cs_lagr_get_bdy_conditions(void)
Return pointer to the main boundary conditions structure.
Definition: cs_lagr.c:1675
cs_real_t * particle_flow_rate
Definition: cs_lagr.h:642
Definition: cs_lagr.h:57
int ivisch
Definition: cs_lagr.h:945
cs_real_t phi_p
Definition: cs_lagr.h:706
cs_lagr_boundary_interactions_t * cs_get_lagr_boundary_interactions(void)
Definition: cs_lagr.c:1739
cs_real_t diameter
Definition: cs_lagr.h:453
Definition: cs_lagr.h:58
Definition: cs_lagr.h:1046
int iclogh
Definition: cs_lagr.h:861
cs_field_t * cvar_r11
Definition: cs_lagr.h:1036
cs_lagr_zone_class_data_t * cs_lagr_init_zone_class_new(int iclass, int izone)
Initialize a new class/boundary zone parameters structure.
Definition: cs_lagr.c:1387
cs_real_t * enc1
Definition: cs_lagr.h:669
Definition: cs_lagr.h:386
int io2
Definition: cs_lagr.h:1049
int npstf
Definition: cs_lagr.h:743
Definition: cs_lagr.h:401
int cs_glob_lagr_nzone_max
int deposition
Definition: cs_lagr.h:217
Field descriptor.
Definition: cs_field.h:121
cs_lagr_physico_chemical_t * cs_get_lagr_physico_chemical(void)
Definition: cs_lagr.c:1619
cs_real_t * e2ch
Definition: cs_lagr.h:1068
cs_real_t rayasp
Definition: cs_lagr.h:349
int coal_profile
Definition: cs_lagr.h:435
cs_real_t w_exit
Definition: cs_lagr.h:280
integer, dimension(:), pointer itypfb
Definition: pointe.f90:122
cs_real_t * uetbor
Definition: cs_lagr.h:992
cs_real_t added_mass_const
Definition: cs_lagr.h:193
Definition: cs_lagr.h:60
int itsco
Definition: cs_lagr.h:587
int ires
Definition: cs_lagr.h:812
int ncharm2
Definition: cs_lagr.h:88
Definition: cs_lagr.h:68
int nstits
Definition: cs_lagr.h:604
int iencma
Definition: cs_lagr.h:821
int nvgaus
Definition: cs_lagr.h:85
cs_lagr_encrustation_t * cs_get_lagr_encrustation(void)
Definition: cs_lagr.c:1607
int clogging
Definition: cs_lagr.h:230
int temperature_profile
Definition: cs_lagr.h:425
Definition: cs_lagr.h:356
cs_real_t * thcdch
Definition: cs_lagr.h:1077
cs_real_t cppart
Definition: cs_lagr.h:329
cs_real_t velocity_magnitude
Definition: cs_lagr.h:444
int fouling
Definition: cs_lagr.h:233
int nvisbr
Definition: cs_lagr.h:99
int ngazem
Definition: cs_lagr.h:1059
void cs_lagr_set_zone_class_foul_index(int iclass, int izone, cs_real_t foul_index)
Set density for a given class and boundary zone.
Definition: cs_lagr.c:1339
int ivitbd
Definition: cs_lagr.h:779
int * b_face_zone_id
Definition: cs_lagr.h:638
int iclgst
Definition: cs_lagr.h:785
int ivisdk
Definition: cs_lagr.h:940
int n_b_zones
Definition: cs_lagr.h:631
cs_real_t denasp
Definition: cs_lagr.h:347
int ihdepm
Definition: cs_lagr.h:867
Definition: cs_lagr.h:71
int itsr13
Definition: cs_lagr.h:550
#define BEGIN_C_DECLS
Definition: cs_defs.h:448
Definition: cs_lagr.h:59
int iencdi
Definition: cs_lagr.h:824
Definition: cs_lagr.h:56
int itytur
Definition: cs_lagr.h:963
Definition: cs_lagr.h:500
int n_user_variables
Definition: cs_lagr.h:236
Definition: cs_lagr.h:49
int iencra
Definition: cs_lagr.h:661
int ireent
Definition: cs_lagr.h:339
int itsfp4
Definition: cs_lagr.h:590
cs_real_t * coke_density
Definition: cs_lagr.h:492
cs_real_t * rho0ch
Definition: cs_lagr.h:1074
int * iusb
Definition: cs_lagr.h:830
cs_real_t cp
Definition: cs_lagr.h:465
Definition: cs_lagr.h:67
cs_real_t * a2ch
Definition: cs_lagr.h:1067
int isuila
Definition: cs_lagr.h:138
int itsti
Definition: cs_lagr.h:573
void cs_lagr_set_zone_class_injection(int iclass, int izone, int number, int freq, int stat)
Set injection parameters for a given class and boundary zone.
Definition: cs_lagr.c:1258
int distribution_profile
Definition: cs_lagr.h:420
char ** nombrd
Definition: cs_lagr.h:894
cs_real_t foul_index
Definition: cs_lagr.h:462
Definition: cs_lagr.h:200
Definition: cs_lagr.h:54
int nusbrd
Definition: cs_lagr.h:79
int itsr23
Definition: cs_lagr.h:560
int coal_number
Definition: cs_lagr.h:438
cs_real_4_t * cs_glob_lagr_b_u_normal
Definition: cs_field_pointer.h:100
int isttio
Definition: cs_lagr.h:131
cs_lagr_particle_counter_t * cs_lagr_get_particle_counter(void)
Get read/write pointer to global particle counter.
Definition: cs_lagr.c:1467
cs_real_t * wmolat
Definition: cs_lagr.h:1057
cs_field_t * scal_t
Definition: cs_lagr.h:1001
Definition: cs_lagr.h:51
int idirla
Definition: cs_lagr.h:154
cs_field_t * cvar_k
Definition: cs_lagr.h:1027
int iflmbd
Definition: cs_lagr.h:767
void cs_lagr_set_zone_class_velocity(int iclass, int izone, int profile, cs_real_t velocity[])
Set velocity parameters for a given class and boundary zone.
Definition: cs_lagr.c:1284
cs_real_t * st_val
Definition: cs_lagr.h:620
int nlayer
Definition: cs_lagr.h:89
cs_real_t cs_real_4_t[4]
vector of 4 floating-point values
Definition: cs_defs.h:309
double cs_real_t
Floating-point value.
Definition: cs_defs.h:296
cs_real_t stat_weight
Definition: cs_lagr.h:468
Definition: cs_lagr.h:900
cs_real_t jamlim
Definition: cs_lagr.h:378
int iencnb
Definition: cs_lagr.h:818
Definition: cs_lagr.h:69
int dlvo
Definition: cs_lagr.h:218
cs_real_t lambda_vdw
Definition: cs_lagr.h:700
cs_real_t espasg
Definition: cs_lagr.h:346
int ivismp
Definition: cs_lagr.h:935
int iencck
Definition: cs_lagr.h:827
cs_real_t * xashch
Definition: cs_lagr.h:1076
Definition: cs_lagr.h:50
int idpvar
Definition: cs_lagr.h:307
int ihdiam
Definition: cs_lagr.h:873
int * imoybr
Definition: cs_lagr.h:849
void cs_lagr_set_zone_class_cp(int iclass, int izone, cs_real_t cp)
Set temperature parameters for a given class and boundary zone.
Definition: cs_lagr.c:1118
int iirayo
Definition: cs_lagr.h:973
void cs_lagr_set_zone_class_stat(int iclass, int izone, int profile, cs_real_t weight, cs_real_t flow)
Set coal parameters for a given class and boundary zone.
Definition: cs_lagr.c:1202
cs_lagr_coal_comb_t * cs_glob_lagr_coal_comb
int nozppm
Definition: cs_lagr.h:970
cs_real_t prefth
Definition: cs_lagr.h:1053
int t_order
Definition: cs_lagr.h:141
cs_real_t water_mass_fraction
Definition: cs_lagr.h:477
cs_field_t * x_eau
Definition: cs_lagr.h:1022
cs_real_t * tprenc
Definition: cs_lagr.h:678
int itsvz
Definition: cs_lagr.h:526
cs_real_t ttclag
Definition: cs_lagr.h:395
cs_real_t valen
Definition: cs_lagr.h:712
int ntersl
Definition: cs_lagr.h:98
cs_lagr_physico_chemical_t * cs_glob_lagr_physico_chemical
cs_field_t * cpro_cp
Definition: cs_lagr.h:1015
int iadded_mass
Definition: cs_lagr.h:190
Definition: cs_lagr.h:730
cs_real_t vmax
Definition: cs_lagr.h:614
int n_temperature_layers
Definition: cs_lagr.h:215
int ncharm
Definition: cs_lagr.h:1063
int inbr
Definition: cs_lagr.h:800
int iang
Definition: cs_lagr.h:806
int nor
Definition: cs_lagr.h:389
int inbrbd
Definition: cs_lagr.h:760
Definition: cs_lagr.h:48
cs_real_33_t * cs_glob_lagr_b_face_proj
cs_field_t * cvar_r33
Definition: cs_lagr.h:1040
int n_b_max_zones
Definition: cs_lagr.h:632
int ltsdyn
Definition: cs_lagr.h:505
int * itsmv2
Definition: cs_lagr.h:584
cs_real_t diftl0
Definition: cs_lagr.h:979
int npstft
Definition: cs_lagr.h:752
int iatc
Definition: cs_lagr.h:1052
cs_lagr_reentrained_model_t * cs_get_lagr_reentrained_model(void)
Provide access to cs_lagr_reentrained_model_t.
Definition: cs_lagr.c:1546
cs_field_t * x_m
Definition: cs_lagr.h:1023
int ih2o
Definition: cs_lagr.h:1048
int roughness
Definition: cs_lagr.h:222
int nusbor
Definition: cs_lagr.h:735
cs_field_t * x_oxyd
Definition: cs_lagr.h:1021
Definition: cs_lagr.h:648
int iscovc
Definition: cs_lagr.h:864
cs_lagr_zone_class_data_t * cs_lagr_get_zone_class_data(int iclass, int izone)
Provide access to class/boundary zone parameters structure.
Definition: cs_lagr.c:1366
cs_lagr_specific_physics_t * cs_get_lagr_specific_physics(void)
Provide access to cs_lagr_specific_physics_t.
Definition: cs_lagr.c:1533
cs_lagr_bc_type_t
Definition: cs_lagr.h:46
Definition: cs_field_pointer.h:174
Definition: cs_field_pointer.h:64
int idistu
Definition: cs_lagr.h:160
int itste
Definition: cs_lagr.h:569
int * itsmv1
Definition: cs_lagr.h:580
void cs_lagr_set_zone_class_density(int iclass, int izone, cs_real_t density)
Set density for a given class and boundary zone.
Definition: cs_lagr.c:1317
cs_real_t w_total
Definition: cs_lagr.h:274
int modcpl
Definition: cs_lagr.h:149
int itske
Definition: cs_lagr.h:535
cs_real_t diameter
Definition: cs_lagr.h:361
cs_real_t * a1ch
Definition: cs_lagr.h:1064
int ltsmas
Definition: cs_lagr.h:509
int itsr22
Definition: cs_lagr.h:555
cs_gnum_t n_g_cumulative_failed
Definition: cs_lagr.h:250
cs_real_t * ehetch
Definition: cs_lagr.h:1073
int iensi3
Definition: cs_lagr.h:910
cs_field_t * luminance
Definition: cs_lagr.h:1018
cs_real_t * solub
Definition: cs_lagr.h:367
void cs_lagr_set_zone_class_diam(int iclass, int izone, int profile, cs_real_t diam, cs_real_t diam_dev)
Set diameter parameters for a given class and boundary zone.
Definition: cs_lagr.c:1230
cs_field_t * cromf
Definition: cs_lagr.h:995
cs_gnum_t n_g_resuspended
Definition: cs_lagr.h:268
cs_real_t * e1ch
Definition: cs_lagr.h:1066
cs_real_t w_new
Definition: cs_lagr.h:277
cs_lagr_visualization_t * cs_get_lagr_visualization(void)
cs_real_t shrinking_diameter
Definition: cs_lagr.h:486
Definition: cs_lagr.h:55
void cs_lagr_set_zone_class_temperature(int iclass, int izone, int profile, cs_real_t *temp, cs_real_t emissivity)
Set temperature parameters for a given class and boundary zone.
Definition: cs_lagr.c:1089
int * iym1
Definition: cs_lagr.h:1061
int itsvy
Definition: cs_lagr.h:523
int iencckbd
Definition: cs_lagr.h:797
double precision, dimension(:,:,:), allocatable density
Definition: atimbr.f90:124
int ivistp
Definition: cs_lagr.h:923
bool steady_bndy_conditions
Definition: cs_lagr.h:640
int iclogt
Definition: cs_lagr.h:858
int npts
Definition: cs_lagr.h:607
cs_real_t trefth
Definition: cs_lagr.h:1054
cs_field_t * cvar_omg
Definition: cs_lagr.h:1033
int ntxerr
Definition: cs_lagr.h:611
int injection_frequency
Definition: cs_lagr.h:408
cs_lagr_deposition_state_t
Definition: cs_lagr.h:66
Definition: cs_lagr.h:959
Definition: cs_lagr.h:375
cs_gnum_t n_g_fouling
Definition: cs_lagr.h:265
int nflagm
Definition: cs_lagr.h:82
Definition: cs_lagr.h:629
int cs_glob_lagr_nclass_max
int ltsthe
Definition: cs_lagr.h:517
cs_real_t * xwatch
Definition: cs_lagr.h:1075
cs_real_t visls0
Definition: cs_lagr.h:985
int itsr12
Definition: cs_lagr.h:545
cs_real_t w_resuspended
Definition: cs_lagr.h:289
cs_field_t * temperature
Definition: cs_lagr.h:1002
cs_field_t * pressure
Definition: cs_lagr.h:998
int iencmabd
Definition: cs_lagr.h:791
cs_lagr_specific_physics_t * cs_glob_lagr_specific_physics
cs_lagr_time_scheme_t * cs_glob_lagr_time_scheme
int cluster
Definition: cs_lagr.h:441
cs_real_t cstham
Definition: cs_lagr.h:696
Definition: cs_lagr.h:721
int itsvx
Definition: cs_lagr.h:520
cs_lagr_internal_condition_t * cs_glob_lagr_internal_conditions
cs_real_t * cp2ch
Definition: cs_lagr.h:1071
int * nbprec
Definition: cs_lagr.h:365
int iangbd
Definition: cs_lagr.h:773
int * b_zone_natures
Definition: cs_lagr.h:636
cs_field_t * cvar_r22
Definition: cs_lagr.h:1038
cs_gnum_t n_g_failed
Definition: cs_lagr.h:271
cs_real_t * temperature
Definition: cs_lagr.h:450
cs_lagr_model_t * cs_glob_lagr_model
cs_real_t * y2ch
Definition: cs_lagr.h:1070
int ihdepv
Definition: cs_lagr.h:870
int ndlaim
Definition: cs_lagr.h:83
int * i_face_zone_id
Definition: cs_lagr.h:650
cs_field_t * cvar_ep
Definition: cs_lagr.h:1030
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:292
cs_lagr_clogging_model_t * cs_glob_lagr_clogging_model
void cs_lagr_solve_time_step(const int itypfb[], const cs_real_t *dt)
Definition: cs_lagr.c:1766
cs_real_t dtp
Definition: cs_lagr.h:392
int iturb
Definition: cs_lagr.h:962
int physical_model
Definition: cs_lagr.h:212
cs_real_t tstatp
Definition: cs_lagr.h:885
cs_real_t epseau
Definition: cs_lagr.h:703
int cs_glob_lagr_log_frequency_n
const cs_lagr_const_dim_t * cs_glob_lagr_const_dim
Definition: cs_lagr.h:242
int nbrclas
Definition: cs_lagr.h:359
void cs_lagr_finalize_bdy_cond(void)
Definition: cs_lagr.c:1701
cs_real_t * wmole
Definition: cs_lagr.h:1060
int iilagr
Definition: cs_lagr.h:120
int iflow
Definition: cs_lagr.h:343
cs_field_t * cpro_viscls
Definition: cs_lagr.h:1012
cs_lagr_clogging_model_t * cs_get_lagr_clogging_model(void)
Definition: cs_lagr.c:1571
int ivit
Definition: cs_lagr.h:809
int ncharm
Definition: cs_lagr.h:969
int precipitation
Definition: cs_lagr.h:232
cs_lagr_boundary_interactions_t * cs_glob_lagr_boundary_interactions
int itpvar
Definition: cs_lagr.h:302
#define END_C_DECLS
Definition: cs_defs.h:449
cs_lagr_source_terms_t * cs_get_lagr_source_terms(void)
Definition: cs_lagr.c:1595
int ico
Definition: cs_lagr.h:1050
cs_real_t * coke_mass_fraction
Definition: cs_lagr.h:483
int itsr11
Definition: cs_lagr.h:540
cs_real_t tpart
Definition: cs_lagr.h:320
int * b_zone_classes
Definition: cs_lagr.h:635
cs_gnum_t n_g_total
Definition: cs_lagr.h:253
cs_real_t mporos
Definition: cs_lagr.h:379
int ivisdm
Definition: cs_lagr.h:927
int iviswat
Definition: cs_lagr.h:953
int inclgt
Definition: cs_lagr.h:855
cs_real_t * ahetch
Definition: cs_lagr.h:1072
cs_real_t * h02ch
Definition: cs_lagr.h:1065
cs_lagr_particle_counter_t * cs_lagr_update_particle_counter(void)
Update global particle counter.
Definition: cs_lagr.c:1484
int nbrgau
Definition: cs_lagr.h:86
cs_real_t rho
Definition: cs_lagr.h:363
cs_real_t dnpenc
Definition: cs_lagr.h:687
cs_gnum_t n_g_cumulative_total
Definition: cs_lagr.h:246
cs_lagr_zone_class_data_t * lagr_zone_class_data
cs_real_t flow_rate
Definition: cs_lagr.h:471
cs_lagr_brownian_t * cs_glob_lagr_brownian
int iencdibd
Definition: cs_lagr.h:794
cs_real_t * y1ch
Definition: cs_lagr.h:1069
int resuspension
Definition: cs_lagr.h:226
cs_lagr_precipitation_model_t * cs_glob_lagr_precipitation_model
int itsr33
Definition: cs_lagr.h:565
int ihsum
Definition: cs_lagr.h:876
cs_lagr_dim_t * cs_glob_lagr_dim
cs_real_t cs_real_33_t[3][3]
3x3 matrix of floating-point values
Definition: cs_defs.h:314
int iflm
Definition: cs_lagr.h:803
cs_lagr_extra_module_t * cs_get_lagr_extra_module(void)
Definition: cs_lagr.c:1750
cs_lagr_brownian_t * cs_get_lagr_brownian(void)
Definition: cs_lagr.c:1631
cs_lnum_t * b_zone_id
Definition: cs_lagr.h:634
int ncharb
Definition: cs_lagr.h:966
int inclg
Definition: cs_lagr.h:852
int iflres
Definition: cs_lagr.h:815
cs_real_t fion
Definition: cs_lagr.h:715
int iviste
Definition: cs_lagr.h:931
cs_field_t * viscl
Definition: cs_lagr.h:1009
cs_lagr_time_step_t * cs_glob_lagr_time_step
cs_lagr_source_terms_t * cs_glob_lagr_source_terms
Definition: cs_lagr.h:656
cs_real_t w_deposited
Definition: cs_lagr.h:283
void cs_lagr_init_c_arrays(int dim_cs_glob_lagr_source_terms[2], cs_real_t **p_cs_glob_lagr_source_terms)
Definition: cs_lagr.c:992
cs_real_t w_fouling
Definition: cs_lagr.h:286
cs_field_t * t_gaz
Definition: cs_lagr.h:1003
cs_lagr_extra_module_t * cs_glob_lagr_extra_module
int n_stat_classes
Definition: cs_lagr.h:235
cs_real_t density
Definition: cs_lagr.h:459
Definition: cs_lagr.h:53
cs_real_t modyeq
Definition: cs_lagr.h:348
int velocity_profile
Definition: cs_lagr.h:415
int ivisck
Definition: cs_lagr.h:950
cs_lagr_encrustation_t * cs_glob_lagr_encrustation
cs_lagr_precipitation_model_t * cs_get_lagr_precipitation_model(void)
Provide access to cs_lagr_precipitation_model_t.
Definition: cs_lagr.c:1559
cs_lagr_time_step_t * cs_get_lagr_time_step(void)
Definition: cs_lagr.c:1583
cs_gnum_t n_g_deposited
Definition: cs_lagr.h:262
int impvar
Definition: cs_lagr.h:312
cs_real_t * coal_mass_fraction
Definition: cs_lagr.h:480
cs_field_t * vel
Definition: cs_lagr.h:1006
cs_real_t phi_s
Definition: cs_lagr.h:709
void cs_lagr_finalize(void)
Definition: cs_lagr.c:1019
Definition: cs_lagr.h:52
cs_real_t diameter_variance
Definition: cs_lagr.h:456
Definition: cs_lagr.h:693
int icp
Definition: cs_lagr.h:976
int iencnbbd
Definition: cs_lagr.h:788
int itsmas
Definition: cs_lagr.h:576
cs_lagr_bdy_condition_t * cs_glob_lagr_bdy_conditions
cs_lagr_internal_condition_t * cs_lagr_get_internal_conditions(void)
Return pointer to the main internal conditions structure.
Definition: cs_lagr.c:1646
int lamvbr
Definition: cs_lagr.h:724
cs_lagr_reentrained_model_t * cs_glob_lagr_reentrained_model
cs_real_t rayasg
Definition: cs_lagr.h:350
cs_gnum_t n_g_new
Definition: cs_lagr.h:256
int ivisv2
Definition: cs_lagr.h:919
Definition: cs_lagr.h:335
int npencr
Definition: cs_lagr.h:664
cs_real_t tmamax
Definition: cs_lagr.h:617
cs_real_t initial_diameter
Definition: cs_lagr.h:489
Definition: cs_lagr.h:70
cs_lnum_t nb_part
Definition: cs_lagr.h:404
void cs_lagr_finalize_internal_cond(void)
Definition: cs_lagr.c:1724
cs_real_t * mp_diss
Definition: cs_lagr.h:369
int diameter_profile
Definition: cs_lagr.h:430
Definition: cs_lagr.h:106
int natom
Definition: cs_lagr.h:1056
Definition: cs_lagr.h:96