1
!-------------------------------------------------------------------------------
3
! This file is part of the Code_Saturne Kernel, element of the
4
! Code_Saturne CFD tool.
6
! Copyright (C) 1998-2009 EDF S.A., France
8
! contact: saturne-support@edf.fr
10
! The Code_Saturne Kernel is free software; you can redistribute it
11
! and/or modify it under the terms of the GNU General Public License
12
! as published by the Free Software Foundation; either version 2 of
13
! the License, or (at your option) any later version.
15
! The Code_Saturne Kernel is distributed in the hope that it will be
16
! useful, but WITHOUT ANY WARRANTY; without even the implied warranty
17
! of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18
! GNU General Public License for more details.
20
! You should have received a copy of the GNU General Public License
21
! along with the Code_Saturne Kernel; if not, write to the
22
! Free Software Foundation, Inc.,
23
! 51 Franklin St, Fifth Floor,
24
! Boston, MA 02110-1301 USA
26
!-------------------------------------------------------------------------------
30
!===============================================================================
32
! INCLUDE THERMOCHIMIE POUR LA PHYSIQUE PARTICULIERE
34
!-------------------------------------------------------------------------------
38
!--> CONSTANTES THERMOCHIMIE
40
! RR --> Constante des gaz parfaits en J/mol/K
41
! TREFTH --> Temperature de reference (K)
42
! VOLMOL --> Volume molaire dans les conditions NTP
43
! T = 0 C et P = 1 atm
46
double precision trefth, prefth, volmol
47
parameter ( rr = 8.31434d0 , &
48
trefth = 25.d0 + tkelvi , &
49
prefth = 1.01325d5 , &
54
! NRGAZ --> Nb de reactions globales en phase gaz
55
! NRGAZM --> Nb maximal de reactions globales en phase gaz
56
! NATO --> Nb d especes atomiques (C,H,..)
57
! NATOM --> Nb maximal d especes atomiques (C,H,..)
58
! NGAZE --> Nb de constituants gazeux elementaires
59
! NGAZEM --> Nb maximal de constituants gazeux elementaires
60
! NGAZG --> Nb d especes globales (ex:Fuel,Oxyd,Prod1,Prod2)
61
! NGAZGM --> Nb maximal d especes globales
62
! NPO --> Nb de points de tabulation
63
! NPOT --> Nb maximal de points de tabulation
64
! TH --> Temperature en Kelvin
65
! EHGAZG(G,IT) --> Enthalpie massique (J/kg) de l espece globale
66
! no G a la temperature T(IT)
67
! WMOLG(G) --> Masse molaire de l espece globale
68
! EHGAZE(G) --> Enthalpie massique (J/kg) constituant gazeux
69
! elementaire no E a la temperature T(IT)
70
! WMOLE(G) --> Masse molaire du constituant gazeux elementaire
71
! WMOLAT(E) --> Masse molaire des atomes (C,H,..)
72
! IATC, IATH --> Pointeur dans WMOLEL pour les ecpeces
73
! IATO, IATN, IATS elementaires (C,H,..)
74
! FS(R) --> Taux de melange pour la reaction gloable R
75
! STOEG(G,R) --> Stoechio en especes globales des reactions
76
! pour l espece no G et pour la reaction no R
77
! CKABSG(G) --> Coefficient d'absorption des especes globales
78
! CKABS1 --> Coefficient d'absorption du melange gazeux
80
! DIFTL0 --> Diffusivite dynamique en kg/(m s)
82
integer ngazgm, ngazem, npot, natom, nrgazm
83
parameter( ngazgm = 25 , ngazem = 20 , &
84
npot = 500 , natom = 5 , nrgazm = 1 )
85
integer iatc, iath, iato, iatn , iats
86
parameter( iatc = 1, iath = 2, iato = 3, iatn = 4 , iats = 5 )
88
integer npo, ngaze, ngazg, nato, nrgaz
89
common / tchppi / npo, ngaze, ngazg, nato, nrgaz
91
double precision th(npot), &
92
ehgaze(ngazem,npot), ehgazg(ngazgm,npot), &
93
wmole(ngazem), wmolg(ngazgm), wmolat(natom), &
94
stoeg(ngazgm,nrgazm), fs(nrgazm), &
95
ckabsg(ngazgm), ckabs1, &
97
! ..v.7..1....v ....2....v....3....v....4....v....5....v....6....v....7.I
98
common / tchppr / th, ehgaze, ehgazg, wmole, wmolg, wmolat, &
99
stoeg, fs, ckabsg, ckabs1, diftl0, xco2, xh2o