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 POUR MATISSE
34
!-------------------------------------------------------------------------------
37
! --> PARAMETRES PHYSIQUES EN DUR
39
! TRFMAT : Temperature de reference pour l'air en degre
40
! RRFMAT : Masse volumique de reference de l'air a TRFMAT degres
41
! CRFMAT : CP de reference pour l'air
42
! XMUMAT : Viscosite moleculaire dynamique de l'air consideree
43
! RTURB0 : Intensite turbulente (selon k = (3/2)*(V_ref*RTURB0/100)**2)
45
double precision trfmat, rrfmat, crfmat, xmumat, rturb0
46
parameter(trfmat=20.d0)
47
parameter(rrfmat=1.177d0)
48
parameter(crfmat=1004.d0)
49
parameter(xmumat=1.85d-5)
50
parameter(rturb0=30.d0)
53
! --> INDICATEUR MATISSE
55
! IMATIS : Indicateur permettant de savoir si on utilise Matisse
57
common /imtmat/ imatis
60
! --> VARIABLES NUMERIQUES
62
! IMPMAT : Unite logique du fichier de resultats
64
common /imtfic/ impmat
66
! IICONR : "pointeur" sur ICONRA (connectivite rayonnement et panaches)
68
common /imtmem/ iiconr
70
! ICNROK : indique si la connectivite pour le rayonnement ICONRA
71
! a ete calculee ( = 1) ou non ( = 0)
73
common /imtray/ icnrok
76
! --> DONNEES GEOMETRIQUES
80
! NPTRAN : Nombre de pas d espace transversal
81
! NPLGRS : Nombre de pas d espace longitudinal
82
! NELGRS : Nombre d'elements par pas longitudinal
83
! NCHEST : Nombre de couche d element dans la zone stockage
84
! ITYPEN : Type d'entreposage
90
! EPREGI : Epaisseur des registres/cloisons amont et aval (en y)
91
! EPCHEM : Epaisseur des cheminees (en y)
92
! HCONVE : Hauteur du convergent eventuel
93
! RCONVE : Rapport du convergent eventuel sur le maillage (>=1)
94
! HCHALI : Hauteur de la cheminee d alimentation
95
! HCHEVA : Hauteur de la cheminee d evacuation
96
! HFTTOI : Hauteur du faite du toit
97
! PTRRES : Pas transversal du reseau de conteneur
98
! FRDTRA : Facteur de reduction transversal du modele/reel
99
! PLGRES : Pas longitudinal du reseau de conteneur
100
! EPCHEL : Epaisseur d une couche d element (zone stockage)
101
! DMCONT : Diametre des conteneurs
102
! HRESO : Hauteur du reseau de colis
103
! HPLEN : Hauteur du plenum inferieur (cas alveole uniquement)
105
integer nptran, nplgrs, nelgrs, nchest, netran
107
common /imtgeo/ nptran, nplgrs, nelgrs, nchest, netran, &
110
double precision epregi, epchem, hconve, rconve, hchali, hcheva
111
double precision hfttoi, ptrres, frdtra, plgres, epchel, dmcont
112
double precision hreso , hplen
113
common /rmtgeo/ epregi, epchem, hconve, rconve, hchali, hcheva, &
114
hfttoi, ptrres, frdtra, plgres, epchel, dmcont, &
118
! --> DONNEES PHYSIQUES
122
! IMDCNT (0 ou 1) : Modelisation des panaches de convection naturelle
123
! ICOFOR (0 ou 1) : Regime hydraulique de circulation forcee
124
! ICONLG (0 ou 1) : Reseau de conteneur en ligne (pas triangulaire sinon)
125
! IALVEO (0 ou 1) : Entreposage en alveole
129
! DTDTMX : Delta temperature max / pas de temps
130
! PUICON : Puissance d'un conteneur
131
! TINIT : Temperature d'air en entree en degres C
132
! TCRIT : Temperature d'air de sortie critique en degres C
133
! EMICON : Emissivite des conteneurs
134
! EMIMUR : Emissivite des murs
135
! HEPCNT : Hauteur d erosion des panaches de convection naturelle
136
! DHPCNT : Debit enthalpique des panaches de convection naturelle
137
! DEBMAS : Debit de circulation forcee
138
! PDCCHA : Perte de charge du diffuseur de cheminee d'ALIMENTATION
139
! PDCFCH : Perte de charge du filtre de cheminee d'ALIMENTATION
140
! DHCHEA : Diametre hydraulique de cheminee d'ALIMENTATION
141
! SDCHEA : Surface debitante de cheminee d'ALIMENTATION
142
! PDCCHE : Perte de charge du diffuseur de cheminee d'EVACUATION
143
! PDCCCH : Perte de charge du clapet de cheminee d'EVACUATION
144
! DHCHES : Diametre hydraulique de cheminee d'EVACUATION
145
! SDCHES : Surface debitante de cheminee d'EVACUATION
146
! PDCALG : Perte de charge porte d'entree AMONT longitudinal
147
! PDCATV : Perte de charge porte d'entree AMONT transversale (sur z)
148
! ARGAMT : Angle d inclinaison du registre AMONT (degre)
149
! PDCSLG : Perte de charge porte de sortie AVAL longitudinale
150
! PDCSTV : Perte de charge porte de sortie AVAL transversale (sur z)
151
! ARGAVL : Angle d inclinaison du registre AVAL (degre)
152
! AMPPDC : Amplification des pertes de charge de reseau
153
! DHALVE : Diametre hydraulique de l'alveole
154
! VITREF : Vitesse de reference pour calculer les pertes de charge
155
! PUITOT : Puissance totale de l'installation
156
! DPVENT : Differentiel de pression athmospherique entree/sortie
158
integer imdcnt, icofor, iconlg, ialveo
159
common /imtphy/ imdcnt, icofor, iconlg, ialveo
161
double precision dtdtmx, puicon, tinit , tcrit , emicon, emimur
162
double precision hepcnt, dhpcnt, debmas, pdccha, pdcfch, dhchea
163
double precision sdchea, pdcche, pdccch, dhches, sdches, pdcalg
164
double precision pdcatv, argamt, pdcslg, pdcstv, argavl, amppdc
165
double precision dhalve, vitref, puitot, debcon, cfecca, cfecma
166
double precision dpvent
167
common /rmtphy/ dtdtmx, puicon, tinit , tcrit , emicon, emimur, &
168
hepcnt, dhpcnt, debmas, pdccha, pdcfch, dhchea, &
169
sdchea, pdcche, pdccch, dhches, sdches, pdcalg, &
170
pdcatv, argamt, pdcslg, pdcstv, argavl, amppdc, &
171
dhalve, vitref, puitot, debcon, cfecca, cfecma, &
175
! --> CARTES 2D ET 3D
179
! NZONMX : nombre de zones maximum pour la definition des
180
! - Cartes 2D des pertes de charges de porte d'entree
181
! - Cartes 2D des pertes de charges de porte de sortie
182
! - Cartes 3D des pertes de charges reseau
183
! - Cartes 3D des puissances
184
! NCARTE : nombre de cartes 2D et 3D a definir
185
! NMTDIR : nombre de directions d'espace
187
! - Indicateurs de numero de carte (pour NZOCAR et VIZCAR)
190
! ICPDCE : indicateur carte 2D pdc de porte d'entree
191
! ICPDCS : indicateur carte 2D pdc de porte de sortie
192
! ICPDCR : indicateur carte 2D pdc de reseau
193
! ICPUIS : indicateur carte 3D puissance
195
! - Indicateurs de direction pour les cartes (NZOCAR, VIZCAR, VZOCAR)
198
! ILIGNE : ligne ( x variable, une ligne = (y;z) constant )
199
! IRANGE : rangee ( y variable, une rangee = (x;z) constant )
200
! IALTIT : altitude ( z variable, une altitude = (x;y) constant )
202
! - Tableaux de donnees
204
! NZOCAR(NMTDIR, NCARTE)
205
! : nombre de zones pour chaque direction de chaque carte
206
! VIZCAR(2, NZONMX, NMTDIR, NCARTE)
207
! : definition du debut et de la fin des zones pour
208
! chaque zone dans chaque direction de chaque carte
209
! c'est un reel puisque l'on peut definir une demi ligne
211
! VZOCAR(NZONMX, NMTDIR)
212
! : valeurs associees a la carte 3D de puissance pour
213
! chaque zone dans chaque direction
214
! en x et y, valeurs comprises entre 0 et 1, indiquant
215
! la fraction de colis representee sur une maille
216
! en z, les valeurs sont renormalisees a l'unite
218
! il n'est plus necessaire de stocker les valeurs pour les
219
! cartes de pertes de charges car il n'y a que 2 valeurs
220
! possibles : 1. perte de charge (colis present)
221
! 0. pas de perte de charge
222
! par defaut, on considere qu'il n'y a pas de perte de
223
! charge ; on repere avec VIZCAR le debut et la fin
224
! des zones dans lesquelles les pertes de charges doivent
227
! - Parametres des zones
230
parameter(nzonmx=100)
232
! - Parametres des types de carte
237
integer icpdce, icpdcs, icpdcr, icpuis
243
! - Parametres des directions
248
integer iligne, irange, ialtit
253
! - Numero des cartes
255
integer nzocar(nmtdir, ncarte)
256
common /imtcar/ nzocar
258
! - Valeur des cartes
260
double precision vizcar(2, nzonmx, nmtdir, ncarte)
261
double precision vcarth(nzonmx, nmtdir)
262
common /rmtcar/ vizcar, vcarth
265
! --> NUMERO DES SCALAIRES
267
! ITAAMT : temperature air ambiant
268
! ITPCMT : temperature de peau des colis
269
! ITPPMT : temperature de peau des parois
270
! (murs et alveoles eventuelles)
272
! ISCA(ITAAMT), ISCA(ITPCMT), ISCA(ITPPMT)
274
integer itaamt, itpcmt, itppmt
280
! --> NUMERO DES COULEURS
282
! - Couleurs d'elements
284
! ICMTDF : couleur par defaut ("df" pour "defaut")
285
! ICMTST : zone de stockage ("st" pour "stockage")
286
! ICMTCI : cheminee d'alimentation ("ci" pour "cheminee inlet")
287
! ICMTCO : cheminee d'evacuation ("co" pour "cheminee outlet")
288
! ICMTRI : registre amont ("ri" pour "registre inlet")
289
! ICMTRO : registre aval ("ro" pour "registre outlet")
290
! ICMTJI : jeu entre colis et registre amont
291
! ("ji" pour "jeu inlet")
292
! ICMTJO : jeu entre colis et registre aval
293
! ("jo" pour "jeu outlet")
295
! ICMTCI et ICMTCO designent la partie des cheminees situees
296
! au dessus des convergents enventuels
297
! ICMTJI et ICMTJO ne sont utilises que pour les configurations
298
! avec alveoles (et il n'y a pas de jeu aval en cathedrale)
300
integer icmtdf, icmtst
301
integer icmtci, icmtco, icmtri, icmtro, icmtji, icmtjo
302
parameter(icmtdf = 0)
303
parameter(icmtst = 8)
304
parameter(icmtci = 3)
305
parameter(icmtco = 6)
306
parameter(icmtri = 2)
307
parameter(icmtro = 4)
308
parameter(icmtji = 7)
309
parameter(icmtjo = 9)
311
! - Couleurs de faces (conditions aux limites)
313
! ICMTFI : entree ("i" pour "in")
314
! ICMTFO : sortie ("o" pour "out")
315
! ICMTFG : sol ("g" pour "ground")
316
! ICMTFC : plafond ("c" pour "ceiling")
317
! ICMTFS : symetrie ("s" pour "symmetry")
318
! ICMTFW : paroi ("w" pour "wall")
320
integer icmtfi, icmtfo, icmtfg, icmtfc, icmtfs, icmtfw
321
parameter(icmtfi = 5)
322
parameter(icmtfo = 1)
323
parameter(icmtfg = 10)
324
parameter(icmtfc = 11)
325
parameter(icmtfs = 12)
326
parameter(icmtfw = 13)
329
! --> VARIABLES DEDUITES
331
! HERCNT : Hauteur d'erosion HEPCNT reduite a un nombre entier de
332
! mailles en altitude (ici pour eviter de faire le calcul
333
! plusieurs fois, et donc de risquer des erreurs)
335
double precision hercnt
336
common /rmtbid/ hercnt