Télécharger execrxt.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ EXECRXT  NOTICE  MAGN      17/02/24    21:15:09     9323           
   2 :                                              DATE     17/02/24
   3 :                                                                                 
   4 :    Procedure EXECRXT                        Voir aussi :
   5 :     ----------------  
   6 :       EXECRXT NDT RXT ;
   7 : 
   8 : 
   9 :                                                                                 
  10 :     Commentaires
  11 :     ____________
  12 : 
  13 :     NDT   ENTIER : Nombre de pas de temps
  14 :     RXT   TABLE  : contenant les informations decrites ci-apres
  15 : 
  16 : 
  17 :     OBJET :
  18 :     -----
  19 : 
  20 :   La procedure EXECRXT calcule, a partir d'un etat initial, l'evolution
  21 :  au cours du temps d'un melange gazeux dans une enceinte fermee.
  22 : 
  23 :   L'etat initial peut etre soit un etat uniforme dans l'espace donne par
  24 :  l'utilisateur soit le resultat d'un calcul precedent.
  25 : 
  26 :   Le melange gazeux a la base est compose d'air mais il peut contenir
  27 :  un ou plusieurs des constituants suivants: vapeur,H2,N2,He,O2,CO,CO2.
  28 : 
  29 : Deux modèles sont prévu pour la vapeur d'eau: un modèle 'gaz parfait'
  30 : comme pour l'ensemble des incondensables, version V0 (modèle par défaut)
  31 : et un modèle gaz réel, spécifique à la vapeur d'eau, version V1.
  32 : Dans le cas du modèle gaz réel pour la vapeur d'eau (version V1), les
  33 : conditions limites à l'injection changent. Il faut donner les
  34 : enthalpies d'injection et non la température d'injection, comme pour
  35 : le gaz parfait.
  36 : 
  37 :   En presence de vapeur la condensation en paroi peut apparaitre si
  38 :  les conditions locales sont reunies: Pvap > Psat.
  39 :   Le modele de condensation est de type Chilton-Colburn associe a une
  40 :  correlation d'echange de type convection naturelle le long d'une plaque
  41 :  verticale.
  42 : 
  43 :   On distingue 3 types de conditions limites: les zones d'injection
  44 :  (breches), les parois thermiques et les parois inertes.
  45 :   - Les zones d'injection (possibilite de trois breches) sont des
  46 :  conditions aux limites de type valeur imposee. Elles concernent la
  47 :  vitesse,la temperature du melange, l'air et les differents constituants
  48 :  selon leur presence: vapeur,H2,N2,He,O2,CO,CO2.
  49 :  Il faut alors preciser le debit massique du melange (kg/s),
  50 :  la temperature d'entree (°C) le debit massique d'air et des autres
  51 :  constituants en fonction du temps  (Voir entree 'scenario').
  52 :   - Les parois thermiques : La vitesse est nulle (ou possibilite de
  53 :  fonctions de paroi (entree FPAROI)). La température est soit imposée au
  54 :  cours du temps (entree TIMP1) soit imposée a une temperature par défaut
  55 :  susceptible d'être modifiée par une procedure perso (entree ECHANP)
  56 :  soit le resultat d'un calcul thermique de paroi (entree THERMP). Il est 
  57 :  possible de coupler les equations de th des parois et du fluide de façon 
  58 :  implicite (entree THERCO).
  59 :  En presence de vapeur ces parois sont susceptibles
  60 :  de condenser. Elles sont impermeables pour tous les incondensables.
  61 :   - Les parois inertes. La vitesse est nulle, elles sont impermeables
  62 :  pour toutes les autres inconnues (temperature vapeur incondensables).
  63 :  Pour l'instant elles font office de plan de symetrie. Elles
  64 :  correspondent au maillage obtenu par difference :
  65 :  Enveloppe - paroi - breche(s).
  66 : 
  67 :   La turbulence des mouvements de gaz est modelisee soit par une
  68 :  viscosite tourbillonnaire constante, soit par un modele de longueur de
  69 :  melange soit par un modele K-Epsilon (entree MODTURB).
  70 :   L'absence de l'entree MODTURB stipule que l'ecoulement est laminaire.
  71 : 
  72 :   Un modele d'aspersion est disponible.
  73 : 
  74 :   Un modele de condensation en masse est en test.
  75 : 
  76 :   La constitution d'un jeu de donnees revient a remplir la table RXT et
  77 :  ensuite a appeler la procedure EXECRXT avec cette table.
  78 : 
  79 :   Les indices et leur significations sont les suivants :
  80 : 
  81 :  rxt          = TABLE ; Definition de la table de travail.
  82 : 
  83 :             Cette table sera completee au moment de l'execution par
  84 :                         trois tables: GEO TBT TIC.
  85 :       la sous table GEO contient tous les objets geometriques crees a
  86 :                         partir des donnees de base.
  87 :       la sous table TBT est la table de travail proprement dite et
  88 :                         contient les autres objets crees necessaires au
  89 :                         calcul hormis les inconnues.
  90 :       la sous table TIC contient les inconnues calculees et sert a
  91 :                         l'initialisation du calcul.
  92 : 
  93 :  !!!  Les entrees de la table RXT fournis par l'utilisateur ne sont donc
  94 :       pas modifiees.
  95 : 
  96 :  Définition du modèle de gaz pour la vapeur:
  97 :   modèle de gaz parfait (modèle par défaut):
  98 :   rxt.'VERSION' = 'V0' ;
  99 :   modèle de gaz réel :
 100 :   rxt.'VERSION' = 'V1' ;
 101 : 
 102 : I/ Definition des objets geometriques
 103 : -------------------------------------
 104 : 
 105 :  rxt.'vtf'    = MAILLAGE ;   Definition du maillage fluide (OBLIGATOIRE)
 106 : 
 107 : 
 108 :  rxt.'breche' = MAILLAGE ;   Definition de la localisation de la breche
 109 :                              (OBLIGATOIRE pour le moment).
 110 :  rxt.'diru1'  = POINT    ;   Direction du champ de vitesse a l'injection
 111 :  rxt.'pi'     =  POIN    ;   Definition d'un point interieur au domaine
 112 :                              fluide ce qui sera utilise pour imposer la
 113 :                              pression en ce point (OBLIGATOIRE).
 114 :  rxt.'axe'    = MAILLAGE ;   axe si 2D AXI
 115 : 
 116 : < rxt.'breche2' = MAILLAGE ; > Seconde breche
 117 : < rxt.'diru2'   = POINT    ; > Direction du champ de vitesse a l'injection
 118 : 
 119 : < rxt.'breche3' = MAILLAGE ; > Troisieme breche.
 120 : < rxt.'diru3'   = POINT    ; >
 121 : 
 122 : II/ Modelisation Physique.
 123 : --------------------------
 124 : 
 125 :  A/ Thermique paroi
 126 :  ------------------
 127 : 
 128 :   Definition de la thermique paroi
 129 : 
 130 :   A-1 Ancienne methode
 131 :   --------------------
 132 :  rxt.'THERMP' = 'LOGIQUE'  ;   VRAI si thermique paroi (mur).
 133 :  rxt.'THERCO' = 'LOGIQUE'  ;   VRAI si couplage implicite th. paroi/fluide 
 134 : 
 135 :  rxt.'vtp'    = 'MAILLAGE' ;  maillage de la paroi.
 136 :  rxt.'LAMBDA' = 'FLOTTANT' ;  conductivite thermique de la paroi
 137 :                                                            (W/m/°C).
 138 :  rxt.'ROCP'   = 'FLOTTANT' ;  Rho*Cp de la paroi (J/m3/°C).
 139 :  rxt.'Tp0'    = 'FLOTTANT' ;  Temperature initiale de la paroi (°C).
 140 :  rxt.'ECHAN'  = 'FLOTTANT' ;  Coefficient d'echange fluide / paroi
 141 :                                                           (W/m**2/°C).
 142 :  Remarque : Il est possible de faire varier le coefficient d'échange
 143 :  fluide / paroi au cours du temps en réactualisant dans une procédure
 144 :  PERSO la valeur de rxt.'TIC'.'KHW' .
 145 : 
 146 :   A-2 Nouvelle methode
 147 :   --------------------
 148 :   Cette nouvelle methode permet de definir les differentes parties
 149 :  constituant la paroi et de leur affecter des proprietes differentes.
 150 :  On procede de la maniere suivante:
 151 :  rxt.'PAROIS'=table                                         ;
 152 :  rxt.'PAROIS'.'Peinture'=table                              ;
 153 :  rxt.'PAROIS'.'Peinture'.'vtp'    =  mpeinture              ;
 154 :  rxt.'PAROIS'.'Peinture'.'LAMBDA' = 0.2                     ;
 155 :  rxt.'PAROIS'.'Peinture'.'ROCP'   = 1.E6                    ;
 156 :  rxt.'PAROIS'.'Peinture'.'Tp0'    = 60.                     ;
 157 :  rxt.'PAROIS'.'Acier'=table                                 ;
 158 :  rxt.'PAROIS'.'Acier'.'vtp'    =   macier                   ;
 159 :  rxt.'PAROIS'.'Acier'.'LAMBDA' = 15.                        ;
 160 :  rxt.'PAROIS'.'Acier'.'ROCP'   = 3.9E6                      ;
 161 :  rxt.'PAROIS'.'Acier'.'Tp0'    = 60.                        ;
 162 :  Il reste necessaire de definir le coefficient d'echange de la
 163 :  surface en contact avec le fluide comme precedemment.
 164 :  rxt.'ECHAN'  = 'FLOTTANT' ;  Coefficient d'echange fluide / paroi
 165 :                                                           (W/m**2/°C).
 166 : 
 167 :  La suite est identique pour les deux methodes
 168 : 
 169 : <rxt.'ECHEXT' = 'LOGIQUE'  ;>  VRAI si echange externe paroi/exterieure.
 170 :  rxt.'parext' = 'MAILLAGE' ;  Maillage de la surface exterieure de la
 171 :                               paroi.
 172 :  rxt.'HEXT'   = 'FLOTTANT' ;  Coefficient d'echange avec l'exterieur
 173 :                                                          (W/m**2/°C).
 174 :  rxt.'TPEXT'  = 'FLOTTANT' ;  Temperature exterieure (°C).
 175 : 
 176 :   Il est possible de definir directement les CHPO des proprietes des murs pour,
 177 :   par exemple modeliser differentes couches de materiaux. 
 178 :    1) On initialise tous les champs (avec 1 valeur reelle) en appelant 
 179 :    EXECRXT sans iteration (il faut penser desactiver le trace des graphs)   
 180 :      rxt.'GRAPH' = FAUX;
 181 :      EXECRXT 0 RXT ;
 182 :    2) On surcharge les champs de RXT en question rxt.'TIC'.'LAMBDA' et 
 183 :       rxt.'TIC'.'ROCP'
 184 :      * creation des CHPO sur les supports concernes (ici Mg et Md)
 185 :      * Ces maillages doivent etre issus du spg des murs
 186 :      cvtp = 'DOMA' RXT.'GEO'.'$vtp' 'CENTRE' ;
 187 :      cWg = 'CONT' Wg;
 188 :      cWd = 'CONT' Wd;
 189 :      * Dans cet exemple, Wg et Wd etaient inclus dans RXT.'vtp' lors de la
 190 :      creation du maillage. 
 191 :      Mg = incl cvtp cWg 'BARY' ;
 192 :      Md = incl cvtp cWd 'BARY' ;
 193 :      l_g = 'MANU' 'CHPO' Mg 1 'SCAL' 15.0 'NATURE' 'DISCRET' ;
 194 :      l_d = 'MANU' 'CHPO' Md 1 'SCAL' 30.0 'NATURE' 'DISCRET' ;
 195 :      * surchargement du CHPO rxt.'TIC'.'LAMBDA'
 196 :      rxt.'TIC'.'LAMBDA' = 'KCHT' rxt.'GEO'.'$vtp' 'SCAL' 'CENTRE'
 197 :           15.0 l_d l_g ;
 198 :      * idem pour rxt.'TIC'.'ROCP'
 199 :      r_g = 'MANU' 'CHPO' Mg 1 'SCAL' (7800.0 '*' 500.0) 'NATURE' 
 200 :             'DISCRET' ;
 201 :      r_d = 'MANU' 'CHPO' Md 1 'SCAL' (780.0 '*' 500.0) 'NATURE' 
 202 :             'DISCRET' ;
 203 :      rxt.'TIC'.'ROCP' = 'KCHT' rxt.'GEO'.'$vtp' 'SCAL' 'CENTRE'
 204 :          (7800.0 '*' 500.0) r_d r_g ;
 205 :     3) On reprend le calcul avec les nouveaux CHPO
 206 :       EXECRXT NBIT RXT
 207 : 
 208 :  B/ Temperatures de paroi imposees.
 209 :  ----------------------------------
 210 : 
 211 :   On peut imposer une temperature de paroi via un coefficient d'echange.
 212 :   Cela permet de modeliser des condenseurs. Trois zones sont prevues
 213 : 
 214 : 
 215 :  rxt.'TIMP1'             = table                             ;
 216 :  rxt.'TIMP1'.'MAILLAGE'  = bidon                             ;
 217 :  rxt.'TIMP1'.'t'         = prog   0.0    19620.              ;
 218 :  rxt.'TIMP1'.'TIMP'      = prog   110.0  110.0               ;
 219 :  rxt.'TIMP1'.'ECHAN'     = 1.e1                              ;
 220 : 
 221 : <rxt.'TIMP2'             = table                             ;>
 222 :            .
 223 :            .
 224 : 
 225 : <rxt.'TIMP3'             = table                             ;>
 226 :            .
 227 :            .
 228 :  La valeur du coefficient d'echange donne est la valeur par defaut s'il n'y
 229 :  a pas de condensation ou est la valeur minoree en cas de condensation.
 230 : 
 231 :    Une nouvelle possibilite est offerte pour definir un nombre
 232 :   quelconque de surface a temperature imposee. On procede de la
 233 :   maniere suivante: (TIMP1 est dans l'exemple, le nom donne a la
 234 :   premiere surface).
 235 :  rxt.'TIMP'              = table                              ;
 236 :  rxt.'TIMP'.'TIMP1'             = table                       ;
 237 :  rxt.'TIMP'.'TIMP1'.'MAILLAGE'  = mailsurf                    ;
 238 :  rxt.'TIMP'.'TIMP1'.'t'         = prog   0.0 900000.0         ;
 239 :  rxt.'TIMP'.'TIMP1'.'TIMP'      = prog  21.     21.           ;
 240 :  rxt.'TIMP'.'TIMP1'.'ECHAN'     = 1.e2                        ;
 241 :  rxt.'TIMP'.'TIMP2'             = table                       ;
 242 :  rxt.'TIMP'.'TIMP2'.'MAILLAGE'  = (bas et paroid et plaf)     ;
 243 :  rxt.'TIMP'.'TIMP2'.'t'         = prog   0.0 900000.0         ;
 244 :            .
 245 :            .
 246 :  La valeur du coefficient d'echange donne est la valeur par defaut s'il n'y
 247 :  a pas de condensation ou est la valeur minoree en cas de condensation.
 248 :  
 249 :  On peut egalement imposer une temperature de paroi constante au cours du temps
 250 :  via un coefficient d'echange.
 251 :  
 252 :  rxt.'ECHANP'            = table                             ;
 253 :  rxt.'ECHANP'.'MAILLAGE' = bidon                             ;
 254 :  rxt.'ECHANP'.'TMUR'     = 110.                              ;
 255 :  rxt.'ECHANP'.'ECHAN'    = 1.e1                              ;
 256 :  
 257 :  L'avantage de ceci est de pouvoir ensuite surcharger les champs crees dans la
 258 :  table TIC afin d'imposer un echange thermique issu de procedures personnelles.
 259 : 
 260 :  C/ Definition des constituants
 261 :  ------------------------------
 262 : 
 263 :  L'air est traite par defaut
 264 : 
 265 :  rxt.'VAPEUR' = LOGIQUE ; Presence de vapeur
 266 :  rxt.'H2'     = LOGIQUE ; Presence d'hydrogene
 267 :  rxt.'HE'     = LOGIQUE ; Presence d'helium
 268 :  rxt.'N2'     = LOGIQUE ; Presence d'azote
 269 :  rxt.'O2'     = LOGIQUE ; Presence d'oxygene
 270 :  rxt.'CO'     = LOGIQUE ; Presence d'oxyde de carbone
 271 :  rxt.'CO2'    = LOGIQUE ; Presence de gaz carbonique
 272 : 
 273 : - Conditions initiales dans le volume fluide
 274 : 
 275 :  rxt.'TF0'    = FLOTTANT ; Temperature initiale du melange en °C.
 276 :  rxt.'PT0'    = FLOTTANT ; Pression totale initiale en Pascal.
 277 : 
 278 :  rxt.'Yvap0'  = FLOTTANT ; Fraction massique initiale de vapeur d'eau dans
 279 :                            le volume fluide (Si VAPEUR).
 280 :  rxt.MOT      = FLOTTANT ; Fraction massique des incondensables initialement
 281 :                             presents dans le volume
 282 :  rxt. MOT peut etre egal a 'Yhe0', 'Yh20', 'Yo20', 'Yn20', 'Yco0', 'Yco20'.
 283 : 
 284 :  REMARQUE :
 285 :  Il est a remarquer que l'air est considere dans ce modele comme un composant
 286 :  a part entiere. Il faudra faire attention lors de la poursuite par un
 287 :  eventuel calcul de combustion.
 288 : 
 289 : * Proprietes physiques du melange gazeux prises par defaut
 290 :  Constantes des gaz
 291 : Rgh2  = 4156.5
 292 : Rghe  = 2078.25
 293 : Rgo2  = 259.8
 294 : Rgn2  = 296.9
 295 : Rgco2 = 188.9
 296 : Rgco  = 296.9
 297 : Rgvap = 461.513
 298 : Rgair = 287.1
 299 : 
 300 : muair =  1.800e-5 :  viscosite dynamique (air)  (kg/m/s)
 301 : alf   =  1.800e-5 :  diffusivite thermique      (m**2/s)
 302 : lamair=  2.580e-2 :  conductivite thermique     (W/m/°C)
 303 : db    =  1.000e-5 :  diffusion Browniene
 304 : Cpvap =  1700.0   :  Chaleur specifique a pression constante
 305 :                      pour la vapeur d'eau       (J/°C/kg)
 306 : Lv    =     2.3E6 :  Chaleur latente            (J/kg)
 307 : 
 308 : Les chaleurs specifiques des autres gaz sont donnees par la
 309 : procedure CALCP
 310 : 
 311 : 
 312 :  D/ Modeles de turbulence.
 313 :  -------------------------
 314 : 
 315 :  L'absence de cette entree signifie que l'on est en laminaire.
 316 :  rxt.'MODTURB'    =  MOT ; Type du modele de turbulence
 317 : 
 318 :   Les possibilites sont :
 319 : 
 320 :  a/  rxt.'MODTURB'    =  'NUTURB'    ; Nu turbulent constant
 321 :      rxt.'NUT'        =  'FLOTTANT'  ; Valeur du Nu turbulent
 322 : 
 323 :  b/  rxt.'MODTURB'    =  'LMEL'      ; Modele de longueur de melange
 324 :      rxt.'LMEL'       =  'FLOTTANT'  ; Valeur de la longueur de melange
 325 : 
 326 :  c/  rxt.'MODTURB'    = 'KEPSILON';
 327 : 
 328 : < rxt.'FPAROI'  = 'LOGIQUE'   ; >
 329 : < rxt.'YP'      = 'FLOTTANT'  ; >
 330 : 
 331 :   Le modele KEPSILON n'est disponible qu'avec l'algorithme
 332 :   semi-explicite (ALGO EFM1). Voir entree ALGO.
 333 : 
 334 :  E/ Aspersion.
 335 :  -------------
 336 : 
 337 :  rxt.'ASPER'  = LOGIQUE ; Variable logique definissant la
 338 :  presence d'aspersion dans le cas test.
 339 : 
 340 :  En cas d'aspersion le logique rxt.'VAPEUR' doit etre = VRAI
 341 : 
 342 :  rxt.'aspinj' = MAILLAGE ; Maillage de la surface d'injection
 343 :                            de la phase dispersee
 344 :  rxt.'toitf'  = MAILLAGE ; Maillage de la surface superieure (par
 345 :                            rapport a la verticale) du volume
 346 :                            fluide
 347 :  rxt.'rod'    = MAILLAGE ; Densite de la phase dispersee
 348 :  rxt.'Cpd'    = MAILLAGE ; Chaleur specifique de la phase dispersee
 349 :  rxt.'scenario'.'vzinj'  = LISTREEL ; Liste des vitesses (verticales)
 350 :                                       d'injection de la phase disp.
 351 :  rxt.'scenario'.'xdinj'  = LISTREEL ; Liste de la fraction volumique
 352 :                                       de la phase disp. a l'inject.
 353 :  rxt.'scenario'.'tdinj'  = LISTREEL ; Liste de la temperature de la
 354 :                                       phase disp. a l'injection.
 355 :  rxt.'scenario'.'ddinj'  = LISTREEL ; Liste du diametre de la phase
 356 :                                       disp. a l'injection.
 357 : 
 358 :  F/ Condensation en masse (en test)
 359 :  ----------------------------------
 360 : 
 361 :  rxt.'CONDMAS' =  MOT ; Type du modele de condensation en masse.
 362 : 
 363 :   La seule possibilite est : 'CMAS0'
 364 : 
 365 :  <rxt.'HMAS' = 'FLOTTANT'>  ; Valeur du coefficient d'echange
 366 :  pour l'elimination du brouillard. Par defaut HMAS = 0.
 367 : 
 368 :  G/ Definition du scenario de calcul
 369 :  -----------------------------------
 370 : 
 371 : - Conditions aux limites pour l'injection et les parois
 372 :   pilotees thermiquement
 373 : 
 374 :   G-1 Ancienne methode
 375 :   --------------------
 376 :   Les tables qui contiennent les donnees d'injection pour chacune 
 377 :   des breches (3 breches possibles) sont rxt.'scenario', rxt.'scenario2'
 378 :   et rxt.'scenario3'.
 379 : 
 380 :  rxt.'scenario'          = TABLE ; Table definissant le scenario.
 381 :  rxt.'scenario'.'t'      = LISTREEL ; Liste des temps en secondes.
 382 :  rxt.'scenario'.'qeau'   = LISTREEL ; Liste des debits de
 383 :                                       vapeur injectes en kg/s.
 384 : 
 385 :  rxt.'scenario'.MOT   = LISTREEL ; Liste des debits d'incondensables
 386 :                injectes en kg/s. MOT peut etre egal a 'qair', 'qhe',
 387 :                'qh2', 'qo2', 'qn2', 'qco', 'qco2'.
 388 :  qair est obligatoire, presence d'air dans tous les cas.
 389 :  Les autres debits sont necessaires si presence de l'incondensable.
 390 :  Si l'utilisateur souhaite stopper une injection (tous les debits 
 391 :  par espece nuls), il convient de redefinir les maillages des breches.
 392 : 
 393 :  rxt.'scenario'.'tinj'   = LISTREEL ; Liste des temperatures
 394 :                           d'injection en fonction du temps.
 395 : 
 396 :  Si Vapeur gaz réel (VERSION V1)
 397 : 
 398 :  rxt.'scenario'.'model'  = mot 'A'        ;
 399 :  rxt.'scenario'.'qeau'   = lqm            ; (en kg/s)
 400 :  La liste des enthalpies est en Joule par kg
 401 :  rxt.'scenario'.'hvap'   = Lhj                                 ;
 402 :  rxt.'scenario'.'qliq'   = prog     nbL*0.                     ;
 403 :  rxt.'scenario'.'hliq'   = prog     nbL*0.                     ;
 404 :  rxt.'scenario'.'qair'   = prog     nbL*0.                     ;
 405 :  rxt.'scenario'.'hair'   = prog     nbL*0.                     ;
 406 :  rxt.'Mliq0' = 0.;
 407 : 
 408 :  Sinon  Vapeur gaz parfait on aura (VERSION V0 par défaut)
 409 : 
 410 :  rxt.'scenario'.'qeau'   = lqm*1.         ;
 411 :  rxt.'scenario'.'qair'   = prog nbL*0.    ;
 412 :  rxt.'scenario'.'tinj'   = prog nbL*180.  ;
 413 : 
 414 : Commentaires :
 415 :  Les modeles d'injection sont:
 416 :  modele A : Temperature d'ebulition
 417 :             on effectue la separation des phases à la temperature de
 418 :             saturation de la vapeur, a la pression totale de l'enceinte.
 419 : 
 420 :  modele B : Temperature moyenne de l'enceinte
 421 : 
 422 :  modele C : Temperature de saturation (pression de vapeur à l'injection)
 423 : 
 424 :  La partie liquide va directement dans le puisard dont la masse est
 425 :  initialisée avec rxt.'Mliq0'
 426 : 
 427 : 
 428 :   G-2 Nouvelle methode
 429 :   --------------------
 430 :    Une autre possibilité a été implementée pour definir un nombre
 431 :   quelconque de breches. On procede de la maniere suivante: (A est le
 432 :   le nom donne a la premiere breche).
 433 :  rxt.'Breches'                =   table                      ;
 434 : 
 435 :  rxt.'Breches'.'A'            =   table                      ;
 436 :  rxt.'Breches'.'A'.'scenario' =   table                      ;
 437 :  rxt.'Breches'.'A'.'Maillage' =   ma                         ;
 438 :  rxt.'Breches'.'A'.'diru'     =  (0. 0. 1.)                  ;
 439 :  rxt.'Breches'.'A'.'scenario'.'t'      = prog   0.0 1000.0   ;
 440 :  rxt.'Breches'.'A'.'scenario'.'qeau'   = prog   0.050  0.050 ;
 441 :  rxt.'Breches'.'A'.'scenario'.'qair'   = prog   0.000  0.000 ;
 442 :  rxt.'Breches'.'A'.'scenario'.'tinj'   = prog   150.0  150.0 ;
 443 :  puis pour la breche B .....
 444 :  rxt.'Breches'.'B'            =   table                      ;
 445 :  rxt.'Breches'.'B'.'scenario' =   table                      ;
 446 :  ....
 447 :  à repeter autant que necessaire
 448 : 
 449 :    De la même manière on peut définir un nombre quelconque de
 450 :  sorties. Pour les sorties deux modèles sont disponibles.
 451 :  a/ On impose un débit de sortie: 'Qext'.
 452 :  b/ On modélise un clapet de décharge dont le débit est 
 453 :     déterminé par la différence de pression enceinte/extérieur: 'Pext'.
 454 :     Il faut alors renseigner la pression extérieure et le coefficient
 455 :     de perte de charge (Kext U ** Bext). Un débit en est déduit et
 456 :     est imposé. Si les conditions sont réunies on peut faire apparaître
 457 :     un col sonique qui limitera de débit.
 458 : 
 459 : REMARQUE: Si la condition de sortie n'est pas activée pendant tout le
 460 :  scénario, il est néanmoins conseillé de la définir durant tous le 
 461 :  scénario, de la maintenir non active en donnant une section nulle
 462 :  (<1.e-10), puis de l'activer en lui donnant une section non nulle.
 463 :  Ceci évite de réinitialiser la matrice de pression.
 464 : 
 465 :   On procède ensuite de la manière suivante: (S est le
 466 :   le nom donné à la premiere sortie).
 467 :  rxt.'Sorties'                =   table                      ;
 468 : 
 469 :  rxt.'Sorties'.'S'            =   table                      ;
 470 :  rxt.'Sorties'.'S'.'scenario' =   table                      ;
 471 :  rxt.'Sorties'.'S'.'Maillage' =   ms                         ;
 472 :  rxt.'Sorties'.'S'.'diru'     =  (0. 0. 1.)                  ;
 473 :  rxt.'Sorties'.'S'.'scenario'.'t'      = prog   0.0 1000.0   ;
 474 :  rxt.'Sorties'.'S'.'scenario'.'Section'= 4.68e-4             ;
 475 : 
 476 : Cas a : 'Qext'
 477 :  rxt.'Sorties'.'S'.'scenario'.'Qext'   = prog   1.     1.    ;
 478 : 
 479 : Cas b : 'Pext'
 480 :  rxt.'Sorties'.'S'.'scenario'.'Pext'   = prog   1.e5   1.e5  ;
 481 :  rxt.'Sorties'.'S'.'scenario'.'Kext'   = prog   1.000  1.000 ;
 482 :  rxt.'Sorties'.'S'.'scenario'.'Bext'   = prog   0.75   0.75  ;
 483 : 
 484 :  puis pour la sortie T .....
 485 :  rxt.'Sorties'.'T'            =   table                      ;
 486 :  rxt.'Sorties'.'T'.'scenario' =   table                      ;
 487 :  ....
 488 :  a repeter autant que necessaire
 489 : 
 490 :   G-3 Flux de masse imposes
 491 :   -------------------------
 492 : 
 493 :   On peut imposer les debits de masse sous forme de flux diffusif.
 494 : Cette possibilite concerne les faibles flux (diffusif) de masse.
 495 : 
 496 : rxt.'XFIMP'=table;
 497 : rxt.'XFIMP'.'SURF1'      = table                             ;
 498 : rxt.'XFIMP'.'SURF1'.'MAILLAGE'=permali                       ;
 499 : rxt.'XFIMP'.'SURF1'.'t'         = listreel    ;
 500 : rxt.'XFIMP'.'SURF1'.'qh2'       = listreel    ;
 501 : 
 502 :  Le mots cles caracterisant les incondensables et la vapeur d'eau sont
 503 : les memes que ceux pour decrire les debits breches.
 504 : 'qeau' , 'qair', 'qhe', 'qh2', 'qo2', 'qn2', 'qco', 'qco2'.
 505 : 
 506 :  H/ Possibilite de definir des recombineurs
 507 :  ------------------------------------------
 508 :  Les recombineurs ne peuvent etre defini uniquement si 
 509 :  les constituants suivants sont presents :
 510 :        H2, N2, O2, VAPEUR
 511 :  *-- Recombineur
 512 :  RXT.'RECOMB' = 'TABLE' ;
 513 :  RXT.'RECOMB' . 1 = 'TABLE' ;
 514 :  RXT.'RECOMB' . 1 . 'PAREXT' = MAILLAGE ;
 515 :  RXT.'RECOMB' . 1 . 'ENTREE' = MAILLAGE ;
 516 :  RXT.'RECOMB' . 1 . 'SORTIE' = MAILLAGE ;
 517 :  RXT.'RECOMB' . 1 . 'direntr' = POINT ;
 518 :  RXT.'RECOMB' . 1 . 'dirsort' = POINT ;
 519 :  RXT.'RECOMB' . 1 . 'NSECTION' = ENTIER ;
 520 :  RXT.'RECOMB' . 1 . 'A'        = FLOTTANT ;
 521 :  RXT.'RECOMB' . 1 . 'B'        = FLOTTANT ;
 522 :  RXT.'RECOMB' . 1 . 'XH2ON'    = FLOTTANT ;
 523 :  RXT.'RECOMB' . 1 . 'XH2OFF'   = FLOTTANT ;
 524 :  RXT.'RECOMB' . 1 . 'BETA'     = FLOTTANT ;
 525 : 
 526 :  On peut definir N recombineurs, il suffit de les rajouter 
 527 :  a la table RXT.'RECOMB' . Pour modeliser un recombineur,
 528 :  plusieurs objets doivent etre definis.
 529 :   - PAREXT : Objet de type maillage permettant de decrire
 530 :     la paroi exterieur du recombineur (sans l'entree ni la
 531 :     sortie).
 532 :   - ENTREE : Maillage representant l'entree du recombineur
 533 :   - SORTIE : Maillage representant la sortie du recombineur
 534 :   - direntr : Point definissant la direction du fluide au
 535 :     niveau de l'entree du recombineur
 536 :   - dirsort : Point definissant la direction du fluide au
 537 :     niveau de la sortie du recombineur
 538 :   - NSECTION : Nombre d'unites inserees dans le bati
 539 :     (1 par defaut)
 540 :  Le taux de consommation d'H2 est calcule par la loi SIEMENS
 541 :  (constructeur des recombineurs modelises, par defaut FR90/1-150)
 542 :   - A et B   : Coefficient de la loi lineaire en pression
 543 :     (0.48D-8 kg/s/Pa et 0.58D-3 kg/s par defaut)
 544 :   - XH2ON    : Fraction molaire d'H2 au dela de laquelle le
 545 :     recombineur entre en fonctionnement (0.005 par defaut)
 546 :   - XH2OFF   : Fraction molaire d'H2 en deça de laquelle le 
 547 :     recombineur passe en arret (0.005 par defaut)
 548 :   - BETA     : Zone de transition affectee par la correction
 549 :     pour les faibles concentrations en H2 (0.2 par defaut)  
 550 :  
 551 :  I/ Choix du modele de flux de condensation
 552 :  ------------------------------------------
 553 :  rxt.'MODCOND' = MOT: choix de l'expression du flux de condensation
 554 : 
 555 :   Les possibilites sont :
 556 :  
 557 :  a/ rxt.'MODCOND' = 'CHIL0' : Jv = Jstand = kc rho (Yv - Yvsat)
 558 :  valeur par defaut.
 559 : 
 560 :  * Modele de type Chilton-Colburn valable uniquement dans le cas de tres
 561 :  faible fraction massique de vapeur.
 562 :  La correlation utilisee est une correlation de convection naturelle.
 563 :  
 564 :  b/ rxt.'MODCOND' = 'CHIL1' : Jv = Jstand / (1-Yvsat)
 565 :   (tire de loi de Fick en masse volumique avec diffusivite 
 566 :   de vapeur dans melange de RALOC)
 567 :  
 568 : III/ Parametres de fonctionnement.
 569 : ----------------------------------
 570 : 
 571 :  rxt.'DT0'    =  FLOTTANT ; valeur du pas de temps.
 572 : <rxt.'IMPR'   =  ENTIER   ;>  0 ou 1 impression supplementaires ou pas
 573 : 
 574 :  rxt.'epsi'   =  FLOTTANT ; precision des maillages lors de
 575 :  l'utilisation de l'operateur ELIM.
 576 : 
 577 :  rxt.'GRAPH' = LOGIQUE ; Variable logique permettant
 578 :  l'obtention de graphes de controle sur les variables moyennes.
 579 : 
 580 :   Il y a la possibilite d'executer une procedure personnelle
 581 :  en tete de la boucle en temps de l'algorithme.
 582 : 
 583 : <rxt.'PERSO'   = 'LOGIQUE' ;>
 584 :  rxt.'PRCPERSO'= 'MOT'     ;  nom de la procedure a executer
 585 :  rxt.'TABPERSO'= 'TABLE'   ;  table associee
 586 : 
 587 : <rxt.'FRPREC'= ENTIER  ;> frequence de recalcul du
 588 :  preconditionnement imposee par l'utilisateur exprimee en
 589 :  terme de nombre de pas de temps. 
 590 :  rxt.'FCPRECT'= ENTIER ;  frequence de recalcul du preconditionneur 
 591 :  en fonction de l'indice de boucle sur les pas de temps (utilise par
 592 :  les procedures PREPAENC et EXEC).
 593 :  Par defaut : 1 
 594 :  rxt.'FCPRECI'= ENTIER ;  frequence de recalcul du preconditionneur
 595 :  en fonction de l'indice de boucle sur la boucle d'iterations utilisee
 596 :  pour resoudre les non-linearites (utilise par les procedures PREPAENC
 597 :  et EXEC).
 598 :  Par defaut : 1 
 599 :  Attention ! : si FCPRECT = 1 ou FCPRECI = 1 (valeurs par defaut), le
 600 :  preconditionnement est effectue a chaque iteration, quelle que soit la
 601 :  valeur de FRPREC. Dans le cas ou FRPREC > 1, il faut donc egalement 
 602 :  renseigner FCPRECT et FCPRECI avec FRPREC.
 603 :  
 604 : <rxt.'STAB' = FLOTTANT ;> Ce reel correspond au facteur de
 605 : stabilisation des elements MACRO d'ancienne generation.
 606 : 
 607 :  rxt.'DISCR' = MOT      ; type de la discretisation (vitesse
 608 :                      et grandeurs scalaires hormis la pression)
 609 :                         (voir NAVI ou MODEL).
 610 : 
 611 :  rxt.'KPRE'  = MOT      ; type de la discretisation pour la
 612 :                           pression (voir NAVI ou MODEL).
 613 : 
 614 :  rxt.'MAXELPRE'= ENTIER ; Nombre maximum d'elements pour lequel on 
 615 :                  utilise une methode directe (Cholesky) pour inverser la
 616 :                  matrice de pression. Par defaut 20000.
 617 : 
 618 :  rxt.'TYPINV'= ENTIER ; Type d'inversion (1) utilise une methode directe
 619 :   (Cholesky) pour inverser les matrices autres que celle de la pression
 620 :   par défaut 3 (Bi CG). Pour les autres possibilités voir notice de KRES
 621 : 
 622 :  rxt.'RESIDU'= FLOTTANT ; résidu toléré pour les méthodes itératives y
 623 :  compris la pression. Par défaut 1.E-20.
 624 : 
 625 :  rxt.'NTIT'= ENTIER ; Nombre d'itérations maximum pour les méthodes
 626 :  itératives y compris la pression. Par défaut 800 .
 627 : 
 628 :  rxt.'ALGO'  = MOT      ; Type d'algorithme
 629 :                           'IMPL' Algorithme semi implicite
 630 :                           'EFM1' Algorithme semi explicite
 631 : 
 632 : <rxt.'DETMAT' = LOGIQUE ;> indique que l'on detruit ou non
 633 :  les matrices a la fin du calcul.
 634 : 
 635 : <rxt.'REINIT' = LOGIQUE ;> indique que l'on reinitialise le
 636 :  calcul pour tenir compte d'eventuels changements dans le modele
 637 :  (ajout d'un composant). Si une modification dans les entrees
 638 :  de la table RXT (ajout ou retrait) est operee et que cette
 639 :  variable n'est pas positionnee a VRAI, le calcul s'arrete.
 640 : 
 641 : <rxt.'CORTEMP' = LOGIQUE ;> indique que l'on utilise ou non
 642 :  le controle du bilan d'energie. Le defaut est VRAI.
 643 : 
 644 : 
 645 : 
 646 :  LIMITATIONS DU MODELE ACTUEL
 647 : 
 648 : * modele gaz parfait
 649 : * modele designe pour les grands pas de temps.
 650 : * Pas de possibilite de plan de symetrie
 651 : 
 652 : Quelques recettes :
 653 : 
 654 : I/ Changement de nature des conditions aux limites
 655 :    Si on repart d'un calcul precedent
 656 :   1/ On recupere la table rxt.'TIC'
 657 :      opti rest 'MONFIC.sauv';
 658 :      rest ;
 659 :      tic = rxt.'TIC' ;
 660 :   2/ On recupere eventuellemnt le maillage
 661 :      opti rest 'MONMAIL.sauv';
 662 :      rest ;
 663 :   3/ On decrit entierement le nouveau probleme
 664 :      rxt=table ;
 665 :      rxt.'vtf'= mon maillage ;
 666 :      rxt.  .... etc
 667 :   4/ On initialise avec la table tic precedente
 668 :      rxt.'TIC'=tic ;
 669 : 
 670 : II/ Reduction de la taille du fichier de sauvetage
 671 :      voir entree DETMAT pour detruire les objets MATRIK
 672 : 
 673 : 
 674 : IV/ Resultats du calcul.
 675 : -----------------------
 676 : 
 677 :  La TABLE  'TIC'  Contient les inconnues calculees 0D et multi-D
 678 :  -----------------------------------------------------------------
 679 : 
 680 :          Indice                          Objet
 681 :       Type    Valeur                  Type   Valeur
 682 :  MOT       SOUSTYPE              MOT       INCO                
 683 : 
 684 : * Cvm Cpm Gamm : capacite calorifique du melange (J/kg/K) et rapport
 685 : * Roj : densite a l'injection (kg/m3)
 686 : 
 687 : * A/ Suivi temporel 0D
 688 : *---------------------
 689 :  MOT    Tps      FLOTTANT : Temps physique
 690 :  MOT    NUPADT   ENTIER   : Numero du pas de temps.
 691 :  MOT    LTPS     LISTREEL : liste des instants calcules.
 692 :  MOT    DT       FLOTTANT : pas de temps
 693 : 
 694 :  MOT    LMAXU    LISTREEL : Liste Temporelle (LT) du module max de la
 695 :                                vitesse (m/s)
 696 :  MOT    Rhom     LISTREEL : LT densite du melange (kg/m3)
 697 :  MOT    Rhomv    LISTREEL : LT densite de vapeur dans le melange (kg/m3)
 698 :                                     (Si VAPEUR VRAI)
 699 :  MOT    Rhomhe   LISTREEL : LT densite de helium dans le melange (kg/m3)
 700 :                                     (Si THE    VRAI)
 701 :  MOT    Rhomh2   LISTREEL : LT densite de h2     dans le melange (kg/m3)
 702 :                                     (Si TH2    VRAI)
 703 :  MOT    Rhomo2   LISTREEL : LT densite de o2     dans le melange (kg/m3)
 704 :                                     (Si TO2    VRAI)
 705 :  MOT    Rhomn2   LISTREEL : LT densite de n2     dans le melange (kg/m3)
 706 :                                     (Si Tn2    VRAI)
 707 :  MOT    Rhomco   LISTREEL : LT densite de CO     dans le melange (kg/m3)
 708 :                                     (Si Tco    VRAI)
 709 :  MOT    Rhomco2  LISTREEL : LT densite de CO2    dans le melange (kg/m3)
 710 :                                     (Si Tco2   VRAI)
 711 :  MOT    Tfm      LISTREEL : LT Temperature moyenne du melange (°C)
 712 :  MOT    Tpm      LISTREEL : LT Temperature moyenne de la paroi (°C)
 713 :                                     (Si THERMP VRAI)
 714 :  MOT    Ltbp1    LISTREEL : LT temperature imposee (TTIMP1)
 715 :  MOT    Ltbp2    LISTREEL : LT temperature imposee (TTIMP2)
 716 :  MOT    Ltbp3    LISTREEL : LT temperature imposee (TTIMP3)
 717 : 
 718 :  MOT    Qc       LISTREEL : LT Debit total de condensation (kg/s)
 719 :  MOT    Qcw      LISTREEL : LT Debit de condensation (THERMP) (kg/s)
 720 :  MOT    Qc1      LISTREEL : LT Debit de condensation (TTIMP1) (kg/s)
 721 :  MOT    Qc2      LISTREEL : LT Debit de condensation (TTIMP2) (kg/s)
 722 :  MOT    Qc3      LISTREEL : LT Debit de condensation (TTIMP3) (kg/s)
 723 :  MOT    Qc0      LISTREEL : LT Debit de condensation (ECHANP) (kg/s)
 724 :  MOT    Econd    LISTREEL : LT energie extraite par condensation (J/m3)
 725 :  MOT    Hcond    LISTREEL : LT enthalpie extraite par condensation(J/m3)
 726 :  MOT    Econv    LISTREEL : LT energie de convection (J/m3)
 727 : 
 728 :  MOT    Easpe    LISTREEL : LT energie extraite par l'aspersion (J/m3) 
 729 :                                     (Si ASPER  VRAI)
 730 :  MOT    Haspe    LISTREEL : LT enthalpie extraite par l'aspersion (J/m3)
 731 :  MOT    Qaspe    LISTREEL : LT debit aspersion  (kg/s)
 732 : 
 733 :  MOT    Remn     LISTREEL : LT energie interne specifique melange (J/m3)
 734 : 
 735 :  MOT    Rgpm     LISTREEL : LT constante gaz parfaits melange (J/kg/K)
 736 :  MOT    Cvm      LISTREEL : LT }capacite calorifique du melange (J/kg/K)
 737 :  MOT    Cpm      LISTREEL : LT }
 738 :  MOT    Gamm     LISTREEL : LT }et rapport
 739 : 
 740 :  MOT    PT       LISTREEL : LT pression thermodynamique
 741 :  MOT    Minj     LISTREEL : LT masses injectees (kg)
 742 :  MOT    Mcond    LISTREEL : LT masses condensees (kg)
 743 :  MOT    Mrest    LISTREEL : LT masses restantes  (kg)
 744 : 
 745 :  MOT    guj      LISTREEL : LT vitesse a l'injection (breche 1) (m/s)
 746 :  MOT    Qj       LISTREEL : LT debit d'injection (breche 1) (kg/s)
 747 :  MOT    Hj       LISTREEL : LT enthalpie a l'injection (breche 1) (J/kg)
 748 :  MOT    Ej       LISTREEL : LT energie a l'injection (breche 1) (J/kg)
 749 : 
 750 :  MOT    guj2     LISTREEL : LT 
 751 :  MOT    Qj2      LISTREEL : LT 
 752 :  MOT    Hj2      LISTREEL : LT    breche 2
 753 :  MOT    Ej2      LISTREEL : LT 
 754 : 
 755 :  MOT    TBP1     FLOTTANT : Temperature de paroi imposee (TTIMP1)
 756 :  MOT    TBP2     FLOTTANT : Temperature de paroi imposee (TTIMP2)
 757 :  MOT    TBP3     FLOTTANT : Temperature de paroi imposee (TTIMP3)
 758 :  MOT    TBP0     CHPOINT  : Temperature de paroi constante (ECHANP)
 759 :  MOT    KHEW     FLOTTANT : coefficient d'echange (convectif)
 760 :                             paroi/fluide  THERMP
 761 :  MOT    KHE1     FLOTTANT : coefficient d'echange (convectif) TIMP1
 762 :  MOT    KHE2     FLOTTANT : coefficient d'echange (convectif) TIMP2
 763 :  MOT    KHE3     FLOTTANT : coefficient d'echange (convectif) TIMP3
 764 :  MOT    KHE0     CHPOINT  : coefficient d'echange (convectif) ECHANP
 765 : 
 766 :  *--- RECOMBINEUR ---*
 767 :  
 768 :  MOT    QIN_H2   LISTREEL : LT debit d'hydrogene total a l'entree 
 769 :                             des recombineurs
 770 :  MOT    QIN_H2O  LISTREEL : LT debit de vapeur total a l'entree 
 771 :                             des recombineurs
 772 :  MOT    QIN_O2   LISTREEL : LT debit d'oxygene total a l'entree 
 773 :                             des recombineurs
 774 :  MOT    QIN_N2   LISTREEL : LT debit d'azote total a l'entree 
 775 :                             des recombineurs
 776 :  MOT    QOUT_H2  LISTREEL : LT debit d'hydrogene total a la sortie
 777 :                             des recombineurs
 778 :  MOT    QOUT_H2O LISTREEL : LT debit de vapeur total a la sortie
 779 :                             des recombineurs
 780 :  MOT    QOUT_O2  LISTREEL : LT debit d'oxygene total a la sortie
 781 :                             des recombineurs
 782 :  MOT    QOUT_N2  LISTREEL : LT debit d'azote total a la sortie
 783 :                             des recombineurs
 784 :  MOT    RECOMB   TABLE    : Table contenant differentes informations
 785 :                             a propos des recombineurs
 786 :  RXT.TIC.RECOMB. 1        : Table du premier recombineur
 787 :   (...)
 788 :  RXT.TIC.RECOMB. 1 . 'DEB' : LT contenant le debit du recombineur 1
 789 :  RXT.TIC.RECOMB. 1 . 'Uin' : LT contenant la vitesse a l'entree
 790 :                              du recombineur 1
 791 :  RXT.TIC.RECOMB. 1 . 'Uout': LT contenant la vitesse a la sortie
 792 :                              du recombineur 1
 793 :  RXT.TIC.RECOMB. 1 . 'Tin' : LT contenant la temperature a l'entree
 794 :                              du recombineur 1
 795 :  RXT.TIC.RECOMB. 1 . 'Tout': LT contenant la temperature a la sortie
 796 :                              du recombineur 1
 797 :  RXT.TIC.RECOMB. 1 . 'Tpla': LT contenant la temperature des plaques
 798 :                              du recombineur 1
 799 :  RXT.TIC.RECOMB. 1 . 'EFF' : LT contenant l'efficacite du recombineur 1
 800 :  RXT.TIC.RECOMB. 1 . 'Hin' : LT contenant l'enthalpie a l'entree
 801 :                              du recombineur 1
 802 :  RXT.TIC.RECOMB. 1 . 'Hout': LT contenant l'enthalpie a la sortie
 803 :                              du recombineur 1
 804 :  RXT.TIC.RECOMB. 1 . 'R?E' : LT contenant la densite de l'espece ? a 
 805 :                              l'entree du recombineur 1
 806 :  RXT.TIC.RECOMB. 1 . 'R?S' : LT contenant la densite de l'espece ? a 
 807 :                              la sortie du recombineur 1
 808 :  RXT.TIC.RECOMB. 1 . 'Y?E' : LT contenant la fraction massique de 
 809 :                              l'espece ? a l'entree du recombineur 1
 810 :  RXT.TIC.RECOMB. 1 . 'Y?S' : LT contenant la fraction massique de 
 811 :                              l'espece ? a la sortie du recombineur 1
 812 : 
 813 : * B/ Grandeurs multi-D
 814 : *---------------------
 815 :  MOT    UN       CHPOINT  : Champ de vitesse au pas de temps courant
 816 :  MOT    UNM      CHPOINT  : Champ de vitesse au pas de temps precedent
 817 :  MOT    ROG      CHPOINT  :
 818 :  MOT    PRES     CHPOINT  : Champ de pression
 819 :  MOT    TF       CHPOINT  : Champ de temperature au pas de temps courant
 820 :  MOT    TFNM     CHPOINT  : temperature au pas de temps precedent
 821 :  MOT    RHO      CHPOINT  : Champ de densite au pas de temps courant
 822 :  MOT    RHONM    CHPOINT  : Champ de densite au pas de temps precedent
 823 :  MOT    Mu       CHPOINT  : viscosite dynamique du melange
 824 :  MOT    NU       CHPOINT  : viscosite cinematique du melange
 825 :  MOT    NUEFF    CHPOINT  : viscosite cinematique effective
 826 : 
 827 : Aspersion
 828 :  MOT    VN       CHPOINT  : Champ de vitesse de la phase dispersee au
 829 :                              pas de temps courant
 830 :  MOT    TD       CHPOINT  : Champ de temperature de la phase dispersee
 831 :                              au pas de temps courant
 832 :  MOT    XD       CHPOINT  : Fraction volumique de la phase dispersee
 833 :                              au pas de temps courant
 834 :  MOT    DD       CHPOINT  : Diametre de la phase dispersee
 835 :                              au pas de temps courant
 836 : 
 837 : 
 838 :  MOT    RAIR     CHPOINT  : densite de l'air (kg)
 839 :  MOT    RVP      CHPOINT  : densite de vapeur (kg)
 840 :  MOT    RHE      CHPOINT  : densite de He     (kg)
 841 :  MOT    RH2      CHPOINT  : densite de H2     (kg)
 842 :  MOT    RO2      CHPOINT  : densite de O2     (kg)
 843 :  MOT    RN2      CHPOINT  : densite de N2     (kg)
 844 :  MOT    RCO2     CHPOINT  : densite de CO2    (kg)
 845 :  MOT    RCO      CHPOINT  : densite de CO     (kg)
 846 : 
 847 :  MOT    YVAP     FLOTTANT / CHPOINT : Fraction (F.M.) massique vapeur
 848 :  MOT    YHE      FLOTTANT / CHPOINT :
 849 :  MOT    YH2      FLOTTANT / CHPOINT : FLOTTANT si l'espece est absente
 850 :  MOT    YO2      FLOTTANT / CHPOINT :
 851 :  MOT    YN2      FLOTTANT / CHPOINT :
 852 :  MOT    YCO2     FLOTTANT / CHPOINT :
 853 :  MOT    YCO      FLOTTANT / CHPOINT :
 854 :  MOT    YAIR     FLOTTANT / CHPOINT :
 855 : 
 856 : 
 857 :  La TABLE GEO contient les informations geometriques en particulier
 858 :  les objets MMODEL.
 859 : 
 860 :      Indice               Objet
 861 :  Type      Valeur   Type        Commentaires
 862 : 
 863 :  MOT    epsi     FLOTTANT : tolerance pour confondre deux points dont
 864 :                                distance est inferieure a epsi
 865 :  MOT    $vtf     MMODEL   :
 866 :  MOT    $menvf   MMODEL   :
 867 :  MOT    $axe     MMODEL   : axe ou plan de symetrie
 868 :  MOT    $vtp     MMODEL   : paroi thermique
 869 :  MOT    $mtp1    MMODEL   : paroi temperature imposee (TIMP1)
 870 :  MOT    $mtp2    MMODEL   : paroi temperature imposee (TIMP2)
 871 :  MOT    $mtp3    MMODEL   : paroi temperature imposee (TIMP3)
 872 :  MOT    $mtp0    MMODEL   : paroi temperature constante (ECHANP)
 873 :  MOT    Pimp     MAILLAGE : maillage type POI1 contenant le point ou
 874 :                                la pression est imposee.
 875 : 
 876 : 
 877 : 

© Cast3M 2003 - Tous droits réservés.
Mentions légales