~ubuntu-branches/ubuntu/precise/code-saturne/precise

« back to all changes in this revision

Viewing changes to include/mati/matiss.h

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-11-01 17:43:32 UTC
  • mto: (6.1.7 sid)
  • mto: This revision was merged to the branch mainline in revision 11.
  • Revision ID: package-import@ubuntu.com-20111101174332-tl4vk45no0x3emc3
Tags: upstream-2.1.0
ImportĀ upstreamĀ versionĀ 2.1.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
!-------------------------------------------------------------------------------
2
 
 
3
 
!     This file is part of the Code_Saturne Kernel, element of the
4
 
!     Code_Saturne CFD tool.
5
 
 
6
 
!     Copyright (C) 1998-2009 EDF S.A., France
7
 
 
8
 
!     contact: saturne-support@edf.fr
9
 
 
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.
14
 
 
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.
19
 
 
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
25
 
 
26
 
!-------------------------------------------------------------------------------
27
 
 
28
 
!                             matiss.h
29
 
 
30
 
!===============================================================================
31
 
 
32
 
!            INCLUDE POUR MATISSE
33
 
 
34
 
!-------------------------------------------------------------------------------
35
 
 
36
 
 
37
 
! --> PARAMETRES PHYSIQUES EN DUR
38
 
 
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)
44
 
 
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)
51
 
 
52
 
 
53
 
! --> INDICATEUR MATISSE
54
 
 
55
 
!     IMATIS : Indicateur permettant de savoir si on utilise Matisse
56
 
integer          imatis
57
 
common /imtmat/  imatis
58
 
 
59
 
 
60
 
! --> VARIABLES NUMERIQUES
61
 
 
62
 
!     IMPMAT : Unite logique du fichier de resultats
63
 
integer          impmat
64
 
common /imtfic/  impmat
65
 
 
66
 
!     IICONR : "pointeur" sur ICONRA (connectivite rayonnement et panaches)
67
 
integer          iiconr
68
 
common /imtmem/  iiconr
69
 
 
70
 
!     ICNROK : indique si la connectivite pour le rayonnement ICONRA
71
 
!              a ete calculee ( = 1) ou non ( = 0)
72
 
integer          icnrok
73
 
common /imtray/  icnrok
74
 
 
75
 
 
76
 
! --> DONNEES GEOMETRIQUES
77
 
 
78
 
!   - Entiers
79
 
 
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
85
 
!                * 1 : Emm
86
 
!                * 0 : Vault
87
 
 
88
 
!   - Reels
89
 
 
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)
104
 
 
105
 
integer          nptran, nplgrs, nelgrs, nchest, netran
106
 
integer          itypen
107
 
common /imtgeo/  nptran, nplgrs, nelgrs, nchest, netran,          &
108
 
                 itypen
109
 
 
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,  &
115
 
                 hreso , hplen
116
 
 
117
 
 
118
 
! --> DONNEES PHYSIQUES
119
 
 
120
 
!   - Entiers
121
 
 
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
126
 
 
127
 
!   - Reels
128
 
 
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
157
 
 
158
 
integer          imdcnt, icofor, iconlg, ialveo
159
 
common /imtphy/  imdcnt, icofor, iconlg, ialveo
160
 
 
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,  &
172
 
                 dpvent
173
 
 
174
 
 
175
 
! --> CARTES 2D ET 3D
176
 
 
177
 
!   - Dimensions
178
 
 
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
186
 
 
187
 
!   - Indicateurs de numero de carte (pour NZOCAR et VIZCAR)
188
 
!       de 1 a NCARTE
189
 
 
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
194
 
 
195
 
!   - Indicateurs de direction pour les cartes (NZOCAR, VIZCAR, VZOCAR)
196
 
!       de 1 a NMTDIR
197
 
 
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 )
201
 
 
202
 
!   - Tableaux de donnees
203
 
 
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
210
 
!                par exemple
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
217
 
 
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
225
 
!                etre activees
226
 
 
227
 
!   - Parametres des zones
228
 
 
229
 
integer   nzonmx
230
 
parameter(nzonmx=100)
231
 
 
232
 
!   - Parametres des types de carte
233
 
 
234
 
integer   ncarte
235
 
parameter(ncarte=4)
236
 
 
237
 
integer   icpdce, icpdcs, icpdcr, icpuis
238
 
parameter(icpdce=1)
239
 
parameter(icpdcs=2)
240
 
parameter(icpdcr=3)
241
 
parameter(icpuis=4)
242
 
 
243
 
!   - Parametres des directions
244
 
 
245
 
integer   nmtdir
246
 
parameter(nmtdir=3)
247
 
 
248
 
integer   iligne, irange, ialtit
249
 
parameter(iligne=1)
250
 
parameter(irange=2)
251
 
parameter(ialtit=3)
252
 
 
253
 
!   - Numero des cartes
254
 
 
255
 
integer          nzocar(nmtdir, ncarte)
256
 
common /imtcar/  nzocar
257
 
 
258
 
!   - Valeur des cartes
259
 
 
260
 
double precision vizcar(2, nzonmx, nmtdir, ncarte)
261
 
double precision vcarth(nzonmx, nmtdir)
262
 
common /rmtcar/  vizcar, vcarth
263
 
 
264
 
 
265
 
! --> NUMERO DES SCALAIRES
266
 
 
267
 
!     ITAAMT : temperature air ambiant
268
 
!     ITPCMT : temperature de peau des colis
269
 
!     ITPPMT : temperature de peau des parois
270
 
!                             (murs et alveoles eventuelles)
271
 
!     A utiliser comme
272
 
!       ISCA(ITAAMT), ISCA(ITPCMT), ISCA(ITPPMT)
273
 
 
274
 
integer   itaamt, itpcmt, itppmt
275
 
parameter(itaamt=1)
276
 
parameter(itpcmt=2)
277
 
parameter(itppmt=3)
278
 
 
279
 
 
280
 
! --> NUMERO DES COULEURS
281
 
 
282
 
!   - Couleurs d'elements
283
 
 
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")
294
 
 
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)
299
 
 
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)
310
 
 
311
 
!   - Couleurs de faces (conditions aux limites)
312
 
 
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")
319
 
 
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)
327
 
 
328
 
 
329
 
! --> VARIABLES DEDUITES
330
 
 
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)
334
 
 
335
 
double precision hercnt
336
 
common /rmtbid/  hercnt
337
 
 
338
 
! FIN