1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
|
/*============================================================================
* Code_Saturne documentation page
*============================================================================*/
/*
This file is part of Code_Saturne, a general-purpose CFD tool.
Copyright (C) 1998-2021 EDF S.A.
This program is free software; you can redistribute it and/or modify it under
the terms of the GNU General Public License as published by the Free Software
Foundation; either version 2 of the License, or (at your option) any later
version.
This program is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
details.
You should have received a copy of the GNU General Public License along with
this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
/*----------------------------------------------------------------------------*/
/*!
\page cs_var_dico Variables and structures reference (C and Fortran)
\section cs_var_dico_intro Introduction
This page is meant to help users find their way through when implementing
user C functions or Fortran subroutines or even developing inside the code
kernel. It provides cross-reference tables containing the names of Fortran
variables and their C counterparts as well as guidelines about
how to manage the mesh entities and the fields (variables, properties, ...).
In addition, some naming conventions are described.
Note: variables with the same name in Fortran and in C are most of the time not listed.
\section cs_var_dico_namingcontent Contents
Cross-reference tables and guidelines are organized as follows:
- \subpage mesh
- \subpage field
- \subpage local
- \subpage field_kw
*/
// _____________________________________________________________________________
/*!
\page mesh How to access and manage the main mesh entities and mesh quantities ?
\section cs_var_dico_mesh_vars Mesh variables
These variables are defined in the files \ref mesh.f90 and \ref cs_mesh.h.
- The Fortran variables are global in the code, accessible directly by their
name.
- Members of the global C structure \c cs_glob_mesh are accessed as: \n
<tt>cs_glob_mesh->name</tt>, \n
\e i.e. adding <tt> cs_glob_mesh-> </tt> in front of the name of the
variable.
- For a more concise syntax, defining a local \c m or \c mesh variable is recommended:
\code{.c}
const cs_mesh_t *m = cs_glob_mesh;
\endcode
so that the shorter name of the local variable may be used.
- In functions which have a <tt> domain </tt> argument,
using <tt> domain->mesh </tt> is recommended, and may be combined
with the use of a local argument (as above) for conciseness.
- Note that the number of elements given here are *local* to a given
rank when running in parallel. Global sizes should only be used in specific
instances, though they are available in \ref cs_mesh_t structures.
Fortran code | C code | Description
------------- | ------------------------------------- | ------------
<tt> ndim | cs_glob_mesh->dim </tt> | Space dimension (always 3)
<tt> ncelet | cs_glob_mesh->n_cells_with_ghosts </tt> | Total number of cells on the local rank \n (n_cells + n_ghost_cells)
<tt> ncel | cs_glob_mesh->n_cells </tt> | Number of cells
<tt> nfac | cs_glob_mesh->n_i_faces </tt> | Number of interior faces
<tt> nfabor | cs_glob_mesh->n_b_faces </tt> | Number of boundary faces
<tt> nnod | cs_glob_mesh->n_vertices </tt> | Number of vertices
- | <tt> cs_glob_mesh->n_b_cells </tt> | Number of boundary cells
<tt> lndfac | cs_glob_mesh->i_face_vtx_connect_size </tt> | Size of the connectivity \n interior faces -> vertices
<tt> lndfbr | cs_glob_mesh->b_face_vtx_connect_size </tt>| Size of the connectivity \n boundary faces -> vertices
<tt> ifacel | cs_glob_mesh->i_face_cells </tt> | Interior faces -> cells connectivity
<tt> ifabor | cs_glob_mesh->b_face_cells </tt> | Boundary faces -> cells connectivity
<tt> ipnfac | cs_glob_mesh->i_face_vtx_idx </tt> | Interior faces -> vertices index
<tt> nodfac | cs_glob_mesh->i_face_vtx_lst </tt> | Interior faces -> vertices connectivity
<tt> ipnfbr | cs_glob_mesh->b_face_vtx_idx </tt> | Boundary faces -> vertices index
<tt> nodfbr | cs_glob_mesh->b_face_vtx_lst </tt> | Boundary faces -> vertices connectivity
- | <tt> cs_glob_mesh->b_cells </tt> | Boundary cell list
<tt> xyznod | cs_glob_mesh->vtx_coord </tt> | Vertex coordinates
\section cs_var_dico_mesh_q_vars Mesh quantity variables
These variables are defined in the files \ref mesh.f90 and
\ref cs_mesh_quantities.h.
- The Fortran variables are global in the code, accessible directly by their
name.
- Members of the global C structure \c cs_glob_mesh_quantities are accessed
as: \n
<tt>cs_glob_mesh_quantities->name</tt>, \n
\e i.e. adding <tt> cs_glob_mesh_quantities-> </tt> in front of the name of
the variable.
- For a more concise syntax, defining a local \c mq variable is recommended:
\code{.c}
const cs_mesh_quantities_t *mq = cs_glob_mesh_quantities;
\endcode
so that the shorter name of the local variable may be used.
- In functions which have a <tt> domain </tt> argument,
using <tt> domain->mesh_quantities </tt> is recommended, and may be combined
with the use of a local argument (as above) for conciseness.
Fortran code | C code | Description
------------- | --------------------------------------------- |-------------
<tt> isympa | cs_glob_mesh_quantities->b_sym_flag </tt> | Symmetry flag for boundary faces
<tt> xyzcen | cs_glob_mesh_quantities->cell_cen </tt> | Cell center coordinates
<tt> surfac | cs_glob_mesh_quantities->i_face_normal </tt> | Surface normal of interior faces
<tt> surfbo | cs_glob_mesh_quantities->b_face_normal </tt> | Surface normal of border faces
<tt> cdgfac | cs_glob_mesh_quantities->i_face_cog </tt> | Center of gravity of interior faces
<tt> cdgfbo | cs_glob_mesh_quantities->b_face_cog </tt> | Center of gravity of border faces
<tt> volume | cs_glob_mesh_quantities->cell_vol </tt> | Cell volume
<tt> surfan | cs_glob_mesh_quantities->i_face_surf </tt> | Surface of interior faces
<tt> surfbn | cs_glob_mesh_quantities->b_face_surf </tt> | Surface of boundary faces
<tt> dist | cs_glob_mesh_quantities->i_dist </tt> | Distance between the cell center and \n the center of gravity of interior faces
<tt> distb | cs_glob_mesh_quantities->b_dist </tt> | Distance between the cell center and \n the center of gravity of border faces
<tt> pond | cs_glob_mesh_quantities->weight </tt> | Interior faces weighting factor
*/
// _____________________________________________________________________________
/*!
\page field How to access and manage variables and properties using the cs_field API?
\ref cs_var_dico_vars "Variables" and \ref cs_var_dico_props "properties" can be accessed both in Fortran and in C using the \ref cs_field.h "cs_field" API.
First, let us specify a few conventions used in code_saturne:
- The *dynamic variables* designation includes velocity, pressure, and
variables related to the turbulence model (<em>k</em>, <em>ε</em>,
<em>R<sub>ij</sub></em>, <em>ω</em>, <em>ϕ</em>, \f$ \overline{f} \f$,
<em>α</em>, <em>ν<sub>t</sub></em>). \n
- The designation “scalar” refers to (usually scalar) variables which are
solution of an advection equation, apart from the *dynamic* variables listed
above; for instance the temperature, scalars which may be passive or not,
user-defined or model-defined.
The mean value of the square of the fluctuations of a “scalar” is a
“scalar”, too.
Scalars may be divided into two groups: **user-defined** scalars
and **model-defined** (sometimes referred to as “specific physics”) scalars.
In Fortran, there are \c nscaus user-defined scalars and \c nscapp
“specific physics” (i.e. model-defined) scalars, with
<tt>nscal = nscaus + nscapp</tt>.
- Depending on the thermal model used; the solved thermal scalar variable may be
the temperature, the enthalpy, or in the case of the compressible module,
the total energy. When the temperature is not the solved thermal variable,
it is usually available as a property.
- When a user scalar represents the mean of the
square of the fluctuations of another scalar (<em>i.e.</em> the average
\f$ \overline{\varphi^\prime\varphi^\prime} \f$ for a fluctuating scalar
<em>ϕ</em>). This can be made either {\em via} the
interface or by declaring that scalar using
\texttt{cs\_parameters\_add\_variable\_variance} in
\texttt{cs\_user\_parameters.f90} (if the scalar in question is not a ``user''
scalar, the selection is made automatically). For instance, if \texttt{j}
and \texttt{k} are ``user'' scalars, the variable $\varphi$ corresponding
to \texttt{k} is the variable number
\texttt{isca(k)=isca(iscavr(j))}.
\end{list}
\par Accessing variables and properties in Fortran:
- The most general way of accessing a field is to use its name: \n\n
<tt>call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("pressure", cvar_pr)
\n pres = cvar_pr(iel)</tt> \n\n
- Both \ref cs_var_dico_vars "variables" and \ref cs_var_dico_props "properties" can be accessed via the
\ref field.f90 "cs_field" API, using the global field indices and indirections arrays, as in the following examples:
- For one-dimensional arrays:\n\n
<tt>call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref numvar::ipr "ipr"), cvar_pr)
\n pres = cvar_pr(iel)</tt>, \n\n
<tt>call \ref field::field_get_val_s "field_get_val_s"(\ref cstphy::icp "icp", cpro_cp) \n
cp = cpro_cp(iel)</tt> \n\n
The values of scalar variable can be accessed as follows:\n\n
<tt>call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref isca "isca"(iscalt)), cvar_scalt) \n
temp = cvar_scalt(iel)</tt> \n\n
- For interleaved multidimensional arrays:\n\n
<tt>call \ref field::field_get_val_v "field_get_val_v"(ivarfl(\ref numvar::iu "iu"), cvar_vel)
\n ux = cvar_vel(1,iel)</tt> \n\n
\ref numvar::ipr "ipr", \ref numvar::iu "iu" are here variable indices and the array ivarfl allows to get the corresponding field indices.\n
\ref isca "isca"(iscalt) is also a variable index (\ref optcal::iscalt "iscalt" is the scalar index of the thermal scalar).\n
\ref cstphy::icp "icp" is directly a field index (there is no equivalent to the array ivarfl for field of type properties).\n\n
- the solved variable number for scalar <tt>j</tt>, where
<tt>1 <= j <= nscal</tt>, is given by the <tt>iscal</tt> array.
So to access the associated field, use:\n
<tt>ivarfl(iscal(j))</tt>). \n\n
- For model-defined scalars, the total scalar number for model scalar <tt>j</tt>, where
<tt>1 <= j <= nscapp</tt>, is given by the <tt>iscapp</tt> array.
So to access the associated field, use:\n
<tt>ivarfl(iscal(iscapp(j)))</tt>). \n\n
- The thermal scalar can be accessed using the <tt>iscalt</tt> index in the list of the
scalars. It's variable number is thus <tt> isca(iscalt) </tt>.\n
if there is no thermal scalar, <tt>iscalt</tt> is equal to -1. \n\n
\par Accessing variables and properties in C:
- The most general way of accessing a field is to use its name: \n\n
<tt>cs_real_t *cvar_pr = cs_field_by_name "field_get_val_s_by_name"("pressure")->val; </tt>\n
then: \n
<tt> cvar_pr(cell_id) </tt> \n\n
- The main \ref cs_var_dico_vars "variables" and \ref cs_var_dico_props "properties" can be accessed using the \ref CS_F_ macro:\n\n
- For one-dimensional arrays:\n
<tt>press = CS_F_(p)->val[cell_id]</tt>, \n
<tt>cp = CS_F_(cp)->val[cell_id]</tt> \n
<tt>temp = CS_F_(t)->val[cell_id]</tt> \n\n
- For multidimensional arrays:\n
<tt>uz = CS_F_(vel)->val[3*cell_id + 2]</tt>\n\n
These arrays can also be cast as follows (for a 3-D array):\n
<tt>\ref cs_real_3_t *cvar_vel = (\ref cs_real_3_t *)CS_F_(vel)->val</tt> \n\n
The cell value can then be accessed as \n
<tt>ux = cvar_vel[cell_id][0]</tt>\n\n
<b>\c vel, \c p, or \c cp </b> are defined in \ref cs_field_pointer.h. \n\n
- Indexed variables (such as user scalars) and indexed properties
are accessed as: \n
<tt>\ref CS_FI_(name,ii-1)->val[cell_id]</tt>. \n\n
- The thermal scalar can be accessed using the \ref cs_thermal_model_field
function: \n
<tt>field_t *tf = cs_thermal_model_field();</tt>\n
then: \n
<tt>tf->val[cell_id]</tt>. \n
if there is no thermal scalar, cs_thermal_model_field returns NULL. \n\n
- In any case, all variables can be accessed using the function \ref cs_field_by_name:\n
<tt>\ref cs_real_t *cvar_pr = \ref cs_field_by_name("pressure")->val</tt> \n\n
- A field's \ref scalar_id key-word should be ≥ 0 if it is a scalar. \n\n
- When a field is user-defined (rather than model-defined),
its \ref cs_field_t::type "type" flag
should also match the \ref CS_FIELD_USER mask: \n
<tt> if (f->type & CS_FIELD_USER) ... </tt>
\remark Note that indexes in C begin at 0, while indexes in Fortran begin at 1.
Thus, Fortran and C loop counters are related in the following as:\n
<tt>cell_id = iel-1</tt>.
Cross-reference tables are available for the variables and properties of the
standard solver and the specific physics features:
- \ref cs_var_dico_vars
- \ref cs_var_dico_props
- \ref cs_var_dico_part
- \ref cs_var_dico_atmo
- \ref cs_var_dico_comb
- \ref cs_var_dico_cfbl
- \ref cs_var_dico_elec
- \ref cs_var_dico_cogz
- \ref cs_var_dico_rayt
\section cs_var_dico_vars Variables
The Fortran variables indexes are defined in the files \ref numvar.f90 (with
the exception of \c ihm and \c iscal, which are respectively defined in \ref
ppincl.f90 and \ref optcal.f90) and the C variables names are defined in
\ref cs_field_pointer.h. \n Note that \c dt is accessible by its name
(using the \ref cs_field_by_name family of functions), but not through
a permanent index..
Fortran code | C code | Description
------------------------------------------------ | ---------------------------- | ------------
<tt> call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ipr), cvar_pr) | CS_F_(p)->val | Pressure
call \ref field::field_get_val_v "field_get_val_v"(ivarfl(\ref iu), cvar_vel) | CS_F_(vel)->val | Velocity
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ivolf2), cvar_voidf) | CS_F_(void_f)->val | Void fraction for Volume of Fluid model
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ik ), cvar_k ) | CS_F_(k)->val | Turbulent kinetic energy \f$ k \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref iep ), cvar_eps) | CS_F_(eps)->val | Turbulent dissipation \f$ \varepsilon \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir11), cvar_r11) | CS_F_(r11)->val | Reynolds stress component \f$ R_{xx} \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir22), cvar_r22) | CS_F_(r22)->val | Reynolds stress component \f$ R_{yy} \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir33), cvar_r33) | CS_F_(r33)->val | Reynolds stress component \f$ R_{zz} \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir12), cvar_r12) | CS_F_(r12)->val | Reynolds stress component \f$ R_{xy} \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir23), cvar_r23) | CS_F_(r23)->val | Reynolds stress component \f$ R_{yz} \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir13), cvar_r13) | CS_F_(r13)->val | Reynolds stress component \f$ R_{xz} \f$
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref iphi), cvar_phi) | CS_F_(phi)->val | \f$ \phi \f$ for \f$ \phi-f_b \f$ model
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ifb ), cvar_fb ) | CS_F_(f_bar)->val | \f$ f_b \f$ for \f$ \phi-f_b \f$ model
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ial ), cvar_al ) | CS_F_(alp_bl)->val | \f$ \alpha \f$ for \f$ Bl-v^2-k \f$ \n or EBRSM model
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref iomg), cvar_omg) | CS_F_(omg)->val | \f$ \omega \f$ for \f$ k-\omega \f$ SST model
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref inusa), cvar_nusa) | CS_F_(nusa)->val | \f$ \widetilde{\nu}_T \f$ for Spalart-Allmaras
call \ref field::field_get_val_v "field_get_val_v"(ivarfl(\ref iuma), cvar_mesh_v) | CS_F_(mesh_u)->val | Mesh velocity
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(isca(\ref ppincl::ihm "ihm")), cvar_hm) | CS_F_(h)->val | Enthalpy
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(isca(\ref optcal::iscalt "iscalt")), cvar_scalt) | CS_F_(t)->val | Temperature </tt>
\section cs_var_dico_props Properties
These properties are defined in the files \ref numvar.f90 and
\ref cs_field_pointer.h.
Fortran code |C code | Description
------------------------------------------------------------------------------------------------|---------------------------------- | ------------
<tt> dt |CS_F_(dt)->val | Local time step
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::iviscl "iviscl", cpro_viscl) |CS_F_(mu)->val | Molecular viscosity
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ivisct "ivisct", cpro_visct) |CS_F_(mu_t)->val | Turbulent dynamic viscosity
call \ref field::field_get_val_s "field_get_val_s"(\ref cstphy::icp "icp", cpro_cp) |CS_F_(cp)->val | Specific heat
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::icrom "icrom", cpro_crom) |CS_F_(rho)->val | Density (at cells)
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ibrom "ibrom", bpro_rho) |CS_F_(rho_b)->val[face_id] | Density (at boundary faces)
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ismago "ismago", cpro_smago) |\ref cs_real_t *cpro_smago = \ref cs_field_by_name("smagorinsky_constant^2")->val| Field id of the anisotropic turbulent viscosity
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::icour "icour", cpro_cour) |\ref cs_real_t *cpro_cour = \ref cs_field_by_name("courant_number")->val | Courant number
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ifour "ifour", cpro_four) |\ref cs_real_t *cpro_four = \ref cs_field_by_name("fourier_number")->val | Fourier number
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::iprtot "iprtot", cpro_prtot) |\ref cs_real_t *cpro_prtot = \ref cs_field_by_name("total_pressure")->val | Total pressure at cell centers
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ivisma "ivisma", cpro_visma_s) |\ref cs_real_t *cpro_visma_s = \ref cs_field_by_name("mesh_viscosity")->val | Mesh velocity viscosity (scalar) for the ALE module
call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::ivisma "ivisma", cpro_visma_s) |\ref cs_real_t *cpro_visma_v = \ref cs_field_by_name("mesh_viscosity")->val | Mesh velocity viscosity (vector) for the ALE module
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::itsrho "itsrho"), cpro_tsrho ) |\ref cs_real_t *cpro_tsrho = \ref cs_field_by_name("dila_st")->val | Global dilatation source terms
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ibeta "ibeta"), cpro_beta ) |\ref cs_real_t *cpro_beta = \ref cs_field_by_name("thermal_expansion")->val | Thermal expansion coefficient
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ipori "ipori", cpro_ipori) |CS_F_(poro)->val | Porosity
call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::iporf "iporf", cpro_iporf) |CS_F_(t_poro)->val | Tensorial porosity
call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::iforbr "iforbr", bpro_forbr) |\ref cs_real_t *bpro_forbr = \ref cs_field_by_name("boundary_forces")->val| Field id of the stresses at boundary
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::iyplbr "iyplbr", bpro_yplus) |\ref cs_real_t *bpro_yplus = \ref cs_field_by_name("yplus")->val | Field id of \f$y^+\f$ at boundary
call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::idtten "idtten", dttens) |\ref cs_real_t *dttens = \ref cs_field_by_name("dttens")->val | Field id for the dttens tensor
call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::itempb "itempb", t_b) |CS_F_(t_b)->val | Boundary temperature </tt>
Some physical properties such as specific heat or diffusivity are often
constant (depending on the model or user parameters).
In that case, these properties are stored as a simple real numbers
rather than in a field over all mesh cells.
- This is the case for the specific heat <em>C<sub>p</sub></em>.
- If <em>C<sub>p</sub></em> is constant, it is based on the reference
value \ref cs_glob_fluid_properties->cp0 (<tt>cp0</tt> in Fortran).\n
When this is the case, \ref cs_glob_fluid_properties->icp (mapped as <tt>icp</tt> in
Fortran) should remain set to 0.\n
When <tt>icp</tt> is initialized to 1 (by the GUI, or in
\ref cs_user_parameters.c, it is automatically reset to the id of
the cell-based property field referenced in the above table.
- This is the same for the diffusivity <em>K</em> of each scalar.
- When <em>K</em> is constant, its value is based on the field's reference
diffusivity, accessible through the scalar field's \ref diffusivity_ref
keyword.
When it is variable, the matching field can be specified and accessed using
the base scalar field's \ref diffusivity_id key (accessible using
\ref cs_field_key_id("diffusivity_id") in C, or <tt>kivisl</tt> in Fortran).\n
For example, for a scalar field <tt>f</tt>:\n
<tt> int k_f_id = \ref cs_field_get_key_int(f, \ref cs_field_key_id("diffusivity_id")); \n
cs_field_t *kf = \ref cs_field_by_id(k_f_id); </tt>
\section cs_var_dico_part Specific physical models
\subsection cs_var_dico_atmo Atmospheric
Defined in \ref optcal.f90, \ref atincl.f90, \ref atvarp.f90 and
\ref cs_field_pointer.h.
Fortran code | C code | Description
------------ | ------ | ------------
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(isca(iscalt)), cvar_scalt) | <tt> CS_F_(pot_t)->val </tt> | Potential temperature
- | <tt> CS_F_(totwt)->val </tt> | Total water content
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref isca(\ref atincl::intdrp "intdrp")), cvar_intdrp) | <tt> CS_F_(ntdrp)->val </tt> | Total number of droplets
call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref isca(\ref atchem::isca_chem "isca_chem"(iesp))), cvar_sc) | <tt> CS_FI_(chemistry,iesp-1)->val </tt>| Chemistry species (indexed)
\subsection cs_var_dico_comb Coal combustion
Defined in \ref ppincl.f90, \ref ppcpfu.f90 and \ref cs_field_pointer.h.
Fortran code | C code | Description
------------------------------------------------------------------------------- | -------------------------------- | ------------
<tt> call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::inp "inp"(iesp)), cvar_inpcl) | \ref CS_FI_(np,iesp-1)->val | Particles per kg for coal class
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ixch "ixch"(iesp)), cvar_xchcl) | \ref CS_FI_(xch,iesp-1)->val | Reactive coal mass fraction for coal class
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ixck "ixck"(iesp)), cvar_xckcl) | \ref CS_FI_(xck,iesp-1)->val | Coke mass fraction for coal class
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ixwt "ixwt"(iesp)), cvar_xwtcl) | \ref CS_FI_(xwt,iesp-1)->val | Water mass fraction for coal class
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ih2 "ih2"(iesp)), cvar_h2cl) | \ref CS_FI_(h2,iesp-1)->val | Mass enthalpy for coal class (permeatic case)
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if1m "if1m"(iesp)), cvar_f1mcl) | \ref CS_FI_(f1m,iesp-1)->val | Mean value light volatiles for coal class
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if2m "if2m"(iesp)), cvar_f2mcl) | \ref CS_FI_(f2m,iesp-1)->val | Mean value heavy volatiles for coal class
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if4m "if4m"), cvar_f4m) | CS_F_(f4m)->val | Oxydant 2 mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if5m "if5m"), cvar_f5m)) | CS_F_(f5m)->val | Oxydant 3 mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if6m "if6m"), cvar_f6m)) | CS_F_(f6m)->val | Water from coal drying mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if7m "if7m"), cvar_f7m)) | CS_F_(f7m)->val | Carbon from coal oxidyzed by O2 mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if8m "if8m"), cvar_f8m)) | CS_F_(f8m)->val | Carbon from coal gasified by CO2 mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if9m "if9m"), cvar_f9m)) | CS_F_(f9m)->val | Carbon from coal gasified by H2O mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifvp2m "ifvp2m"), cvar_fvp2m) | CS_F_(fvp2m)->val | f1f2 variance
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iyco2 "iyco2"), cvar_yco2) | CS_F_(yco2)->val | CO2 fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iyhcn "iyhcn"), cvar_yhnc) | CS_F_(yhcn)->val | HCN fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iyno "iyno"), cvar, yno) | CS_F_(yno)->val | NO fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iynh3 "iynh3"), cvar_ynh3) | CS_F_(ynh3)->val | NH3 enthalpy
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::ihox "ihox"), cvar_hox) | CS_F_(hox)->val | Ox enthalpy </tt>
\subsection cs_var_dico_cfbl Compressible
Defined in \ref ppincl.f90 and \ref cs_field_pointer.h.
Fortran code | C code | Description
------------------------------------------------------------------------------ | ----------------------------- | ------------
<tt> call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ienerg "ienerg"), cvar_energ) | CS_F_(e_tot)->val | Total energy
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::itempk "itempk"), cvar_tempk) | CS_F_(t_kelvin)->val | Temperature, in Kelvin </tt>
\subsection cs_var_dico_elec Electric arcs
Defined in \ref ppincl.f90 and \ref cs_field_pointer.h.
Fortran code | C code | Description
----------------------------------------------------------------------------- | ---------------------------------- | ------------
<tt> call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("elec_pot_r", cvar_potr) | CS_F_(potr)->val | Electric potential, real part
call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("elec_pot_i", cvar_poti) | CS_F_(poti)->val | Electric potential, imaginary part
call \ref field::field_get_val_v_by_name "field_get_val_v_by_name"("vec_potential", cvar_potva) | CS_F_(potva)->val | Vector potential
call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("esl_fraction_01", cvar_ycoel_01) | \ref CS_FI_(ycoel,iesp-1)->val | Constituent mass fraction </tt>
\subsection cs_var_dico_cogz Gas combustion
Defined in \ref ppincl.f90 and \ref cs_field_pointer.h.
Fortran code | C code | Description
------------------------------------------------------------------------- | -------------------------- | ------------
<tt> call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifm "ifm"), cvar_fm) | CS_F_(fm)->val | Mixture fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifp2m "ifp2m"), cvar_fp2m) | CS_F_(fp2m)->val | Mixture fraction variance
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifsm "ifsm"), cvar_fsm) | CS_F_(fsm)->val | Soot mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::inpm "inpm"), cvar_npm) | CS_F_(npm)->val | Soot precursor number
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::iygfm "iygfm"), cvar_ygfm) | CS_F_(ygfm)->val | Fresh gas fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::iyfm "iyfm"), cvar_yfm) | CS_F_(yfm)->val | Mass fraction
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::iyfp2m "iyfp2m"), cvar_yfp2m) | CS_F_(yfp2m)->val | Mass fraction variance
call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::icoyfp "icoyfp"), cvar_coyfp) | CS_F_(coyfp)->val | Mass fraction covariance </tt>
\subsection cs_var_dico_rayt Radiative transfer
Defined in \ref cs_field_pointer.h.
C code | Description
--------------------------------------- | ------------
<tt> CS_F_(rad_lumin)->val </tt> | Radiative luminance
<tt> CS_F_(rad_q)->val </tt> | Radiative flux
<tt> CS_FI_(rad_ets,iesp-1)->val </tt> | Radiative flux explicit source term
<tt> CS_FI_(rad_its,iesp-1)->val </tt> | Radiative flux implicit source term
<tt> CS_FI_(rad_abs,iesp-1)->val </tt> | Radiative absorption
<tt> CS_FI_(rad_emi,iesp-1)->val </tt> | Radiative emission
<tt> CS_FI_(rad_cak,iesp-1)->val </tt> | Radiative absorption coefficient
<tt> CS_F_(qinci)->val </tt> | Radiative incident radiative flux density
<tt> CS_F_(xlam)->val </tt> | Wall thermal conductivity
<tt> CS_F_(epa)->val </tt> | Wall thickness
<tt> CS_F_(emissivity)->val </tt> | Wall emissivity
<tt> CS_F_(fnet)->val </tt> | Boundary radiative flux
<tt> CS_F_(fconv)->val </tt> | Boundary radiative convective flux
<tt> CS_F_(hconv)->val </tt> | Radiative exchange coefficient
\subsection cs_var_dico_multiphase Eulerian-Eulerian multiphase flows
Defined in \ref cs_field_pointer.h.
C code | Description
--------------------------------------- | ------------
<tt> CS_FI_(yf_ncond,inc)->val </tt> | Non-condensable gas mass fractions
<tt> CS_FI_(qp,ip)->val </tt> | Particles turbulent kinetic energy Q2
<tt> CS_FI_(qfp,ip)->val </tt> | Covariance of the turbulent Q12
<tt> CS_FI_(qfpxx,ip)->val </tt> | XX component of qfp
<tt> CS_FI_(qfpxy,ip)->val </tt> | XY component of qfp
<tt> CS_FI_(qfpxz,ip)->val </tt> | XZ component of qfp
<tt> CS_FI_(qfpyx,ip)->val </tt> | YX component of qfp
<tt> CS_FI_(qfpyy,ip)->val </tt> | YY component of qfp
<tt> CS_FI_(qfpyz,ip)->val </tt> | YZ component of qfp
<tt> CS_FI_(qfpzx,ip)->val </tt> | ZX component of qfp
<tt> CS_FI_(qfpzy,ip)->val </tt> | ZY component of qfp
<tt> CS_FI_(qfpzz,ip)->val </tt> | ZZ component of qfp
<tt> CS_FI_(gamma,ip)->val </tt> | Interfacial mass transfer
<tt> CS_FI_(ia,ip)->val </tt> | Interfacial area
<tt> CS_FI_(x2,ip)->val </tt> | Droplets x2
<tt> CS_FI_(d32,ip)->val </tt> | Droplets Sauter mean diameter
<tt> CS_FI_(drag,ipcpl)->val </tt> | Drag between phases
<tt> CS_FI_(ad_mass,ip)->val </tt> | Added mass
<tt> CS_FI_(th_diff,ip)->val </tt> | Thermal diffusivity
<tt> CS_FI_(th_diff_t,ip)->val </tt> | Turbulent thermal diffusivity
<tt> CS_FI_(drho_dp,ip)->val </tt> | dRho over dP
<tt> CS_FI_(drho_dh,ip)->val </tt> | dRho over dH
<tt> CS_FI_(tau12_t,ip)->val </tt> | Turbulent tau12 for particles
<tt> CS_FI_(lift,ip)->val </tt> | Particles lift
<tt> CS_FI_(disp_t,ip)->val </tt> | Particles turbulent dispersion
<tt> CS_FI_(drift_vel,ip)->val </tt> | Particles drift velocity
*/
// _____________________________________________________________________________
/*!
\page local How to name common local variables ?
The following table provides a non-exhaustive list of local variables which
are used in the code in a recurring manner.
Fortran code | C code | Description
------------------ | ---------- | ------------
<tt> iel | cell_id </tt> | Cell index
<tt> ifac | face_id </tt> | Face index
<tt> ig | g_id </tt> | Interior face number of associated groups (OpenMP)
<tt> it | t_id </tt> | Interior face number of threads (OpenMP)
<tt> idimtr | tr_dim </tt> | Indicator for tensor perodicity of rotation
<tt> flumas | i_massflux </tt> | Mass flux at interior faces
<tt> flumab | b_massflux </tt> | Mass flux at boundary faces
<tt> viscf | i_visc </tt> | \f$ \mu_\fij \dfrac{S_\fij}{\ipf \jpf} \f$ \n at interior faces for the r.h.s.
<tt> viscb | b_visc </tt> | \f$ \mu_\fib \dfrac{S_\fib}{\ipf \centf} \f$ \n at border faces for the r.h.s.
<tt> smbrp | rhs </tt> | Right hand side \f$ \vect{Rhs} \f$
\section cs_var_dico_conv Local naming convention for fields (Fortran and C)
Rules have been stablished for local names denoting fields, depending on their nature. The convention, applying both in Fortran and in C, is as follows:
- The first letter of the name indicates the location at which the field values are defined:
- \b c for values at the cell centers.
- \b i for values at the interior faces.
- \b b for values at the boundary faces.
- The next three letters indicate if the field is a variable (at the current time step or the previous time step) or a property:
- \b var for variables at the current time step.
- \b vara for variables at the previous time step.
- \b pro for properties.
- An underscore \b _ follows.
- Finally, the <b> short name </b> of the variable/property is specified. This short name is built from the variable/property Fortran index, removing the \c i at the beginning of the word.
The following examples ilustrate this convention:
\c cvar_pr: Values of the variable pressure field defined at the cell centers, at the current time step. \n
\c cvara_pr: Values of the variable pressure field defined at the cell centers, at the previous time step. \n
\c cpro_cp: Values of the property specific heat defined field at the cell centers. \n
*/
// _____________________________________________________________________________
/*!
\page field_kw List of main field keys
A non-exhaustive \ref field_keywords "list of field keys" which are used in the code is provided.
*/
|
