Télécharger enceinte.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ ENCEINTE NOTICE  CHAT      11/09/12    21:16:01     7124           
   2 :                                              DATE     11/09/12
   3 :                                                                                 
   4 :   Procedure ENCEINTE                       Voir aussi :
   5 :     ------------------  
   6 :         ENCEINTE NDT RXT ;
   7 :        
   8 : 
   9 : 
  10 :                                                                                 
  11 :     Commentaires
  12 :     ____________
  13 : 
  14 :     NDT   ENTIER : Nombre de pas de temps
  15 :     RXT   TABLE  : contenant les informations decrites ci-apres
  16 : 
  17 :                                                                                 
  18 :     OBJET :
  19 :     -----
  20 : 
  21 :   La procedure ENCEINTE calcule, a partir d'un etat initial, l'evolution
  22 :  au cours du temps d'un melange gazeux dans une enceinte fermee.
  23 : 
  24 :   L'etat initial peut etre soit un etat uniforme dans l'espace donne par
  25 :  l'utilisateur soit le resultat d'un calcul precedent.
  26 : 
  27 :   Le melange gazeux a la base est compose d'air mais il peut contenir
  28 :  un ou plusieurs des constituants suivants: vapeur,H2,N2,He,O2,CO,CO2.
  29 : 
  30 :   En presence de vapeur la condensation en paroi peut apparaitre si
  31 :  les conditions locales sont reunies: Pvap > Psat.
  32 :   Le modele de condensation est de type Chilton-Colburn associe a une
  33 :  correlation d'echange de type convection naturelle le long d'une plaque
  34 :  verticale.
  35 : 
  36 :   On distingue 3 types de conditions limites: les zones d'injection
  37 :  (breches), les parois thermiques et les parois inertes.
  38 :   - Les zones d'injection (possibilite de trois breches) sont des
  39 :  conditions aux limites de type valeur imposee. Elles concernent la
  40 :  vitesse,la temperature du melange, l'air et les differents constituants
  41 :  selon leur presence: vapeur,H2,N2,He,O2,CO,CO2.
  42 :  Il faut alors preciser le debit massique du melange (Kg/s),
  43 :  la temperature d'entree (°C) le debit massique d'air et des autres
  44 :  constituants en fonction du temps  (Voir entree 'scenario').
  45 :   - Les parois thermiques : La vitesse est nulle (ou possibilite de
  46 :  fonctions de paroi (entree FPAROI)). La temperature est soit imposee au
  47 :  cours du temps (entree TIMP1) soit imposee a une temperture par defaut 
  48 :  susceptible d'etre modiofiee par une procedure perso (entree ECHANP)
  49 :  soit le resultat d'un calcul thermique de paroi (entree THERMP). Il est 
  50 :  possible de coupler les equations de th des parois et du fluide de façon 
  51 :  implicite (entree THERCO).
  52 :  En presence de vapeur ces parois sont susceptibles
  53 :  de condenser. Elles sont impermeables pour tous les incondensables.
  54 :   - Les parois inertes. La vitesse est nulle, elles sont impermeables
  55 :  pour toutes les autres inconnues (temperature vapeur incondensables).
  56 :  Pour l'instant elles font office de plan de symetrie. Elles
  57 :  correspondent au maillage obtenu par difference :
  58 :  Enveloppe - paroi - breche(s).
  59 : 
  60 :   La turbulence des mouvements de gaz est modelisee soit par une
  61 :  viscosite tourbillonnaire constante, soit par un modele de longueur de
  62 :  melange soit par un modele K-Epsilon (entree MODTURB).
  63 :   L'absence de l'entree MODTURB stipule que l'ecoulement est laminaire.
  64 : 
  65 :   Un modele d'aspersion est disponible.
  66 : 
  67 :   La constitution d'un jeu de donnees revient a remplir la table RXT et
  68 :  ensuite a appeler la procedure ENCEINTE avec cette table.
  69 : 
  70 :   Les indices et leur significations sont les suivants :
  71 : 
  72 :  rxt          = TABLE ; Definition de la table de travail.
  73 : 
  74 :             Cette table sera completee au moment de l'execution par
  75 :                         trois tables: GEO TBT TIC.
  76 :       la sous table GEO contient tous les objets geometriques crees a
  77 :                         partir des donnees de base.
  78 :       la sous table TBT est la table de travail proprement dite et
  79 :                         contient les autres objets crees necessaires au
  80 :                         calcul hormis les inconnues.
  81 :       la sous table TIC contient les inconnues calculees et sert a
  82 :                         l'initialisation du calcul.
  83 : 
  84 :  !!!  Les entrees de la table RXT fournis par l'utilisateur ne sont donc
  85 :       pas modifiees.
  86 : 
  87 : I/ Definition des objets geometriques
  88 : -------------------------------------
  89 : 
  90 :  rxt.'vtf'    = MAILLAGE ;   Definition du maillage fluide (OBLIGATOIRE)
  91 : 
  92 : 
  93 :  rxt.'breche' = MAILLAGE ;   Definition de la localisation de la breche
  94 :                              (OBLIGATOIRE pour le moment).
  95 :  rxt.'diru1'  = POINT    ;   Direction du champ de vitesse a l'injection
  96 :  rxt.'pi'     =  POIN    ;   Definition d'un point interieur au domaine
  97 :                              fluide ce qui sera utilise pour imposer la
  98 :                              pression en ce point (OBLIGATOIRE).
  99 :  rxt.'axe'    = MAILLAGE ;   axe si 2D AXI
 100 : 
 101 : < rxt.'breche2' = MAILLAGE ; > Seconde breche
 102 : < rxt.'diru2'   = POINT    ; > Direction du champ de vitesse a l'injection
 103 : 
 104 : < rxt.'breche3' = MAILLAGE ; > Troisieme breche.
 105 : < rxt.'diru3'   = POINT    ; >
 106 : 
 107 : II/ Modelisation Physique.
 108 : --------------------------
 109 : 
 110 :  A/ Thermique paroi
 111 :  ------------------
 112 : 
 113 :   Definition de la thermique paroi
 114 : 
 115 :  rxt.'THERMP' = 'LOGIQUE'  ;   VRAI si thermique paroi (mur).
 116 :  rxt.'THERCO' = 'LOGIQUE'  ;   VRAI si couplage implicite th. paroi/fluide 
 117 : 
 118 :  rxt.'vtp'    = 'MAILLAGE' ;  maillage de la paroi.
 119 :  rxt.'LAMBDA' = 'FLOTTANT' ;  conductivite thermique de la paroi
 120 :                                                            (W/m/°C).
 121 :  rxt.'ROCP'   = 'FLOTTANT' ;  Rho*Cp de la paroi (J/m3/°C).
 122 :  rxt.'Tp0'    = 'FLOTTANT' ;  Temperature initiale de la paroi (°C).
 123 :  rxt.'ECHAN'  = 'FLOTTANT' ;  Coefficient d'echange fluide / paroi
 124 :                                                           (W/m**2/°C).
 125 : <rxt.'ECHEXT' = 'LOGIQUE'  ;>  VRAI si echange externe paroi/exterieure.
 126 :  rxt.'parext' = 'MAILLAGE' ;  Maillage de la surface exterieure de la
 127 :                               paroi.
 128 :  rxt.'HEXT'   = 'FLOTTANT' ;  Coefficient d'echange avec l'exterieur
 129 :                                                          (W/m**2/°C).
 130 :  rxt.'TPEXT'  = 'FLOTTANT' ;  Temperature exterieure (°C).
 131 : 
 132 :   Il est possible de definir directement les CHPO des proprietes des murs pour,
 133 :   par exemple modeliser differentes couches de materiaux. 
 134 :    1) On initialise tous les champs (avec 1 valeur reelle) en appelant 
 135 :    ENCEINTE sans iteration (il faut penser desactiver le trace des graphs)   
 136 :      rxt.'GRAPH' = FAUX;
 137 :      ENCEINTE 0 RXT ;
 138 :    2) On surcharge les champs de RXT en question rxt.'TIC'.'LAMBDA' et 
 139 :       rxt.'TIC'.'ROCP'
 140 :      * creation des CHPO sur les supports concernes (ici Mg et Md)
 141 :      * Ces maillages doivent etre issus du spg des murs
 142 :      cvtp = 'DOMA' RXT.'GEO'.'$vtp' 'CENTRE' ;
 143 :      cWg = 'CONT' Wg;
 144 :      cWd = 'CONT' Wd;
 145 :      * Dans cet exemple, Wg et Wd etaient inclus dans RXT.'vtp' lors de la
 146 :      creation du maillage. 
 147 :      Mg = incl cvtp cWg 'BARY' ;
 148 :      Md = incl cvtp cWd 'BARY' ;
 149 :      l_g = 'MANU' 'CHPO' Mg 1 'SCAL' 15.0 'NATURE' 'DISCRET' ;
 150 :      l_d = 'MANU' 'CHPO' Md 1 'SCAL' 30.0 'NATURE' 'DISCRET' ;
 151 :      * surchargement du CHPO rxt.'TIC'.'LAMBDA'
 152 :      rxt.'TIC'.'LAMBDA' = 'KCHT' rxt.'GEO'.'$vtp' 'SCAL' 'CENTRE'
 153 :           15.0 l_d l_g ;
 154 :      * idem pour rxt.'TIC'.'ROCP'
 155 :      r_g = 'MANU' 'CHPO' Mg 1 'SCAL' (7800.0 '*' 500.0) 'NATURE' 
 156 :             'DISCRET' ;
 157 :      r_d = 'MANU' 'CHPO' Md 1 'SCAL' (780.0 '*' 500.0) 'NATURE' 
 158 :             'DISCRET' ;
 159 :      rxt.'TIC'.'ROCP' = 'KCHT' rxt.'GEO'.'$vtp' 'SCAL' 'CENTRE'
 160 :          (7800.0 '*' 500.0) r_d r_g ;
 161 :     3) On reprend le calcul avec les nouveaux CHPO
 162 :       ENCEINTE NBIT RXT
 163 : 
 164 :  B/ Temperatures de paroi imposees.
 165 :  ----------------------------------
 166 : 
 167 :   On peut imposer une temperature de paroi via un coefficient d'echange.
 168 :   Cela permet de modeliser des condenseurs. Trois zones sont prevues
 169 : 
 170 : 
 171 :  rxt.'TIMP1'             = table                             ;
 172 :  rxt.'TIMP1'.'MAILLAGE'  = bidon                             ;
 173 :  rxt.'TIMP1'.'t'         = prog   0.0    19620.              ;
 174 :  rxt.'TIMP1'.'TIMP'      = prog   110.0  110.0               ;
 175 :  rxt.'TIMP1'.'ECHAN'     = 1.e1                              ;
 176 : 
 177 : <rxt.'TIMP2'             = table                             ;>
 178 :            .
 179 :            .
 180 : 
 181 : <rxt.'TIMP3'             = table                             ;>
 182 :            .
 183 :            .
 184 : 
 185 :  La valeur du coefficient d'echange donnee est la valeur par defaut s'il n'y
 186 :  a pas de condensation ou est la valeur minoree en cas de condensation.
 187 :  
 188 :  On peut egalement imposer une temperature de paroi constante au cours du temps
 189 :  via un coefficient d'echange.
 190 :  
 191 :  rxt.'ECHANP'            = table                             ;
 192 :  rxt.'ECHANP'.'MAILLAGE' = bidon                             ;
 193 :  rxt.'ECHANP'.'TMUR'     = 110.                              ;
 194 :  rxt.'ECHANP'.'ECHAN'    = 1.e1                              ;
 195 :  
 196 :  L'avantage de ceci est de pouvoir ensuite surcharger les champs crees dans la
 197 :  table TIC afin d'imposer un echange thermique issu de procedures personnelles.
 198 : 
 199 :  C/ Definition des constituants
 200 :  ------------------------------
 201 : 
 202 :  L'air est traite par defaut
 203 : 
 204 :  rxt.'VAPEUR' = LOGIQUE ; Presence de vapeur
 205 :  rxt.'H2'     = LOGIQUE ; Presence d'hydrogene
 206 :  rxt.'HE'     = LOGIQUE ; Presence d'helium
 207 :  rxt.'N2'     = LOGIQUE ; Presence d'azote
 208 :  rxt.'O2'     = LOGIQUE ; Presence d'oxygene
 209 :  rxt.'CO'     = LOGIQUE ; Presence d'oxyde de carbone
 210 :  rxt.'CO2'    = LOGIQUE ; Presence de gaz carbonique
 211 : 
 212 : - Conditions initiales dans le volume fluide
 213 : 
 214 :  rxt.'TF0'    = FLOTTANT ; Temperature initiale du melange en °C.
 215 :  rxt.'PT0'    = FLOTTANT ; Pression totale initiale en Pascal.
 216 : 
 217 :  rxt.'Yvap0'  = FLOTTANT ; Fraction massique initiale de vapeur d'eau dans
 218 :                            le volume fluide (Si VAPEUR).
 219 :  rxt.MOT      = FLOTTANT ; Fraction massique des incondensables initialement
 220 :                             presents dans le volume
 221 :  rxt. MOT peut etre egal a 'Yhe0', 'Yh20', 'Yo20', 'Yn20', 'Yco0', 'Yco20'.
 222 : 
 223 :  REMARQUE :
 224 :  Il est a remarquer que l'air est considere dans ce modele comme un composant
 225 :  a part entiere. Il faudra faire attention lors de la poursuite par un
 226 :  eventuel calcul de combustion.
 227 : 
 228 : * Proprietes physiques du melange gazeux prises par defaut
 229 :  Constantes des gaz
 230 : Rgh2  = 4156.5
 231 : Rghe  = 2078.25
 232 : Rgo2  = 259.8
 233 : Rgn2  = 296.9
 234 : Rgco2 = 188.9
 235 : Rgco  = 296.9
 236 : Rgvap = 461.513
 237 : Rgair = 287.1
 238 : 
 239 : muair =  1.800e-5 :  viscosite dynamique (air)  (Kg/m/s)
 240 : alf   =  1.800e-5 :  diffusivite thermique      (m**2/s)
 241 : lamair=  2.580e-2 :  conductivite thermique     (W/m/°C)
 242 : db    =  1.000e-5 :  diffusion Browniene
 243 : Cpvap =  1700.0   :  Chaleur specifique a pression constante
 244 :                      pour la vapeur d'eau       (J/°C/Kg)
 245 : Lv    =     2.3E6 :  Chaleur latente            (J/Kg)
 246 : 
 247 : Les chaleurs specifiques des autres gaz sont donnees par la
 248 : procedure CALCP
 249 : 
 250 : 
 251 :  D/ Modeles de turbulence.
 252 :  -------------------------
 253 : 
 254 :  L'absence de cette entree signifie que l'on est en laminaire.
 255 :  rxt.'MODTURB'    =  MOT ; Type du modele de turbulence
 256 : 
 257 :   Les possibilites sont :
 258 : 
 259 :  a/  rxt.'MODTURB'    =  'NUTURB'    ; Nu turbulent constant
 260 :      rxt.'NUT'        =  'FLOTTANT'  ; Valeur du Nu turbulent
 261 : 
 262 :  b/  rxt.'MODTURB'    =  'LMEL'      ; Modele de longueur de melange
 263 :      rxt.'LMEL'       =  'FLOTTANT'  ; Valeur de la longueur de melange
 264 : 
 265 :  c/  rxt.'MODTURB'    = 'KEPSILON'; (Non disponible pour l'instant)
 266 : 
 267 : < rxt.'FPAROI'  = 'LOGIQUE'   ; >
 268 : < rxt.'YP'      = 'FLOTTANT'  ; >
 269 : 
 270 :  E/ Aspersion.
 271 :  -------------
 272 : 
 273 :  rxt.'ASPER'  = LOGIQUE ; Variable logique definissant la
 274 :  presence d'aspersion dans le cas test.
 275 : 
 276 :  En cas d'aspersion le logique rxt.'VAPEUR' doit etre = VRAI
 277 : 
 278 :  rxt.'aspinj' = MAILLAGE ; Maillage de la surface d'injection
 279 :                            de la phase dispersee
 280 :  rxt.'toitf'  = MAILLAGE ; Maillage de la surface superieure (par
 281 :                            rapport a la verticale) du volume
 282 :                            fluide
 283 :  rxt.'rod'    = MAILLAGE ; Densite de la phase dispersee
 284 :  rxt.'Cpd'    = MAILLAGE ; Chaleur specifique de la phase dispersee
 285 :  rxt.'scenario'.'vzinj'  = LISTREEL ; Liste des vitesses (verticales)
 286 :                                       d'injection de la phase disp.
 287 :  rxt.'scenario'.'xdinj'  = LISTREEL ; Liste de la fraction volumique
 288 :                                       de la phase disp. a l'inject.
 289 :  rxt.'scenario'.'tdinj'  = LISTREEL ; Liste de la temperature de la
 290 :                                       phase disp. a l'injection.
 291 :  rxt.'scenario'.'ddinj'  = LISTREEL ; Liste du diametre de la phase
 292 :                                       disp. a l'injection.
 293 : 
 294 :  F/ Condensation en masse (Non disponible pour l'instant)
 295 :  --------------------------------------------------------
 296 : 
 297 : 
 298 :  G/ Definition du scenario de calcul
 299 :  -----------------------------------
 300 : 
 301 : 
 302 : - Conditions aux limites pour l'injection et les parois
 303 :   pilotees thermiquement
 304 :   Les tables qui contiennent les donnees d'injection pour chacune 
 305 :   des breches (3 breches possibles) sont rxt.'scenario', rxt.'scenario2'
 306 :   et rxt.'scenario3'.
 307 : 
 308 :  rxt.'scenario'          = TABLE ; Table definissant le scenario.
 309 :  rxt.'scenario'.'t'      = LISTREEL ; Liste des temps en secondes.
 310 :  rxt.'scenario'.'qeau'   = LISTREEL ; Liste des debits de
 311 :                                       vapeur injectes en kg/s.
 312 : 
 313 :  rxt.'scenario'.MOT   = LISTREEL ; Liste des debits d'incondensables
 314 :                injectes en kg/s. MOT peut etre egal a 'qair', 'qhe',
 315 :                'qh2', 'qo2', 'qn2', 'qco', 'qco2'.
 316 :  qair est obligatoire, presence d'air dans tous les cas.
 317 :  Les autres debits sont necessaires si presence de l'incondensable.
 318 :  Si l'utilisateur souhaite stopper une injection (tous les debits 
 319 :  par espece nuls), il convient de redefinir les maillages des breches.
 320 : 
 321 :  rxt.'scenario'.'tinj'   = LISTREEL ; Liste des temperatures
 322 :                           d'injection en fonction du temps.
 323 : 
 324 :  H/ Possibilite de definir des recombineurs
 325 :  ------------------------------------------
 326 :  Les recombineurs ne peuvent etre defini uniquement si 
 327 :  les constituants suivants sont presents :
 328 :        H2, N2, O2, VAPEUR
 329 :  *-- Recombineur
 330 :  RXT.'RECOMB' = 'TABLE' ;
 331 :  RXT.'RECOMB' . 1 = 'TABLE' ;
 332 :  RXT.'RECOMB' . 1 . 'PAREXT' = MAILLAGE ;
 333 :  RXT.'RECOMB' . 1 . 'ENTREE' = MAILLAGE ;
 334 :  RXT.'RECOMB' . 1 . 'SORTIE' = MAILLAGE ;
 335 :  RXT.'RECOMB' . 1 . 'direntr' = POINT ;
 336 :  RXT.'RECOMB' . 1 . 'dirsort' = POINT ;
 337 : 
 338 :  On peut definir N recombineurs, il suffit de les rajouter 
 339 :  a la table RXT.'RECOMB' . Pour modeliser un recombineur,
 340 :  plusieurs objets doivent etre definis.
 341 :   - PAREXT : Objet de type maillage permettant de decrire
 342 :     la paroi exterieur du recombineur (sans l'entree ni la
 343 :     sortie).
 344 :   - ENTREE : Maillage representant l'entree du recombineur
 345 :   - SORTIE : Maillage representant la sortie du recombineur
 346 :   - direntr : Point definissant la direction du fluide au
 347 :     niveau de l'entree du recombineur
 348 :   - dirsort : Point definissant la direction du fluide au
 349 :     niveau de la sortie du recombineur
 350 :  
 351 :  I/ Choix du modele de flux de condensation
 352 :  ------------------------------------------
 353 :  rxt.'MODCOND' = MOT: choix de l'expression du flux de condensation
 354 : 
 355 :   Les possibilites sont :
 356 :  
 357 :  a/ rxt.'MODCOND' = 'CHIL0' : Jv = Jstand = kc rho (Yv - Yvsat) 
 358 :  valeur par defaut.
 359 : 
 360 :  * Modele de type Chilton-Colburn valable uniquement dans le cas de tres
 361 :  faible fraction massique de vapeur.
 362 :  La correlation utilisee est une correlation de convection naturelle.
 363 :  
 364 :  b/ rxt.'MODCOND' = 'CHIL1' : Jv = Jstand / (1-Yvsat) 
 365 :   (tire de loi de Fick en masse volumique avec diffusivite 
 366 :   de vapeur dans melange de RALOC)
 367 : 
 368 : III/ Parametres de fonctionnement.
 369 : ----------------------------------
 370 : 
 371 :  rxt.'DT0'    =  FLOTTANT ; valeur du pas de temps.
 372 : <rxt.'IMPR'   =  ENTIER   ;>  0 ou 1 impression supplementaires ou pas
 373 : 
 374 :  rxt.'epsi'   =  FLOTTANT ; precision des maillages lors de
 375 :  l'utilisation de l'operateur ELIM.
 376 : 
 377 :  rxt.'GRAPH' = LOGIQUE ; Variable logique permettant
 378 :  l'obtention de graphes de controle sur les variables moyennes.
 379 : 
 380 :   Il y a la possibilite d'executer une procedure personnelle
 381 :  en tete de la boucle en temps de l'algorithme.
 382 : 
 383 : <rxt.'PERSO'   = 'LOGIQUE' ;>
 384 :  rxt.'PRCPERSO'= 'MOT'     ;  nom de la procedure a executer
 385 :  rxt.'TABPERSO'= 'TABLE'   ;  table associee
 386 : 
 387 :  rxt.'FRPREC'= ENTIER ; frequence de recalcul du
 388 :  preconditionnement imposee par l'utilisateur exprimee en
 389 :  terme de nombre de pas de temps. Si on veut que la
 390 :  procedure calcule elle-meme la frequence de
 391 :  repreconditionnement, il ne faut pas specifier cette valeur
 392 :  et definir l'indice suivant.
 393 : 
 394 : <rxt.'STAB' = FLOTTANT ;> Ce reel correspond au facteur de
 395 : stabilisation des elements MACRO d'ancienne generation.
 396 : 
 397 :  rxt.'DISCR' = MOT      ; type de la discretisation (vitesse
 398 :                      et grandeurs scalaires hormis la pression)
 399 :                         (voir NAVI ou MODEL).
 400 :  rxt.'KPRE'  = MOT      ; type de la discretisation pour la
 401 :                           pression (voir NAVI ou MODEL).
 402 : 
 403 :  rxt.'ALGO'  = MOT      ; Type d'algorithme
 404 :                           'IMPL' Algorithme semi implicite
 405 :                           'EFM1' Algorithme semi explicite
 406 : 
 407 : <rxt.'DETMAT' = LOGIQUE ;> indique que l'on detruit ou non
 408 :  les matrices a la fin du calcul.
 409 : 
 410 : <rxt.'REINIT' = LOGIQUE ;> indique que l'on reinitialise le
 411 :  calcul pour tenir compte d'eventuels changements dans le modele
 412 :  (ajout d'un composant). Si une modification dans les entrees
 413 :  de la table RXT (ajout ou retrait) est operee et que cette
 414 :  variable n'est pas positionnee a VRAI, le calcul s'arrete.
 415 : 
 416 : <rxt.'CORTEMP' = LOGIQUE ;> indique que l'on utilise ou non
 417 :  le controle du bilan d'energie. Le defaut est VRAI.
 418 : 
 419 : 
 420 : 
 421 :  LIMITATIONS DU MODELE ACTUEL
 422 : 
 423 : * modele gaz parfait
 424 : * modele designe pour les grands pas de temps.
 425 : * Pas de possibilite de plan de symetrie
 426 : 
 427 : Quelques recettes :
 428 : 
 429 : I/ Changement de nature des conditions aux limites
 430 :    Si on repart d'un calcul precedent
 431 :   1/ On recupere la table rxt.'TIC'
 432 :      opti rest 'MONFIC.sauv';
 433 :      rest ;
 434 :      tic = rxt.'TIC' ;
 435 :   2/ On recupere eventuellemnt le maillage
 436 :      opti rest 'MONMAIL.sauv';
 437 :      rest ;
 438 :   3/ On decrit entierement le nouveau probleme
 439 :      rxt=table ;
 440 :      rxt.'vtf'= mon maillage ;
 441 :      rxt.  .... etc
 442 :   4/ On initialise avec la table tic precedente
 443 :      rxt.'TIC'=tic ;
 444 : 
 445 : II/ Reduction de la taille du fichier de sauvetage
 446 :      voir entree DETMAT pour detruire les objets MATRIK
 447 : 
 448 : 
 449 : IV/ Resultats du calcul.
 450 : -----------------------
 451 : 
 452 :  La TABLE  'TIC'  Contient les inconnues calculees 0D et multi-D
 453 :  -----------------------------------------------------------------
 454 : 
 455 :          Indice                          Objet
 456 :       Type    Valeur                  Type   Valeur
 457 :  MOT       SOUSTYPE              MOT       INCO                
 458 : 
 459 : * Cvm Cpm Gamm : capacite calorifique du melange (J/kg/K) et rapport
 460 : * Roj : densite a l'injection (kg/m3)
 461 : 
 462 : * A/ Suivi temporel 0D
 463 : *---------------------
 464 :  MOT    Tps      FLOTTANT : Temps physique
 465 :  MOT    NUPADT   ENTIER   : Numero du pas de temps.
 466 :  MOT    LTPS     LISTREEL : liste des instants calcules.
 467 :  MOT    DT       FLOTTANT : pas de temps
 468 : 
 469 :  MOT    LMAXU    LISTREEL : Liste Temporelle (LT) du module max de la
 470 :                                vitesse (m/s)
 471 :  MOT    Rhom     LISTREEL : LT densite du melange (kg/m3)
 472 :  MOT    Rhomv    LISTREEL : LT densite de vapeur dans le melange (kg/m3)
 473 :                                     (Si VAPEUR VRAI)
 474 :  MOT    Rhomhe   LISTREEL : LT densite de helium dans le melange (kg/m3)
 475 :                                     (Si THE    VRAI)
 476 :  MOT    Rhomh2   LISTREEL : LT densite de h2     dans le melange (kg/m3)
 477 :                                     (Si TH2    VRAI)
 478 :  MOT    Rhomo2   LISTREEL : LT densite de o2     dans le melange (kg/m3)
 479 :                                     (Si TO2    VRAI)
 480 :  MOT    Rhomn2   LISTREEL : LT densite de n2     dans le melange (kg/m3)
 481 :                                     (Si Tn2    VRAI)
 482 :  MOT    Rhomco   LISTREEL : LT densite de CO     dans le melange (kg/m3)
 483 :                                     (Si Tco    VRAI)
 484 :  MOT    Rhomco2  LISTREEL : LT densite de CO2    dans le melange (kg/m3)
 485 :                                     (Si Tco2   VRAI)
 486 :  MOT    Tfm      LISTREEL : LT Temperature moyenne du melange (°C)
 487 :  MOT    Tpm      LISTREEL : LT Temperature moyenne de la paroi (°C)
 488 :                                     (Si THERMP VRAI)
 489 :  MOT    Ltbp1    LISTREEL : LT temperature imposee (TTIMP1)
 490 :  MOT    Ltbp2    LISTREEL : LT temperature imposee (TTIMP2)
 491 :  MOT    Ltbp3    LISTREEL : LT temperature imposee (TTIMP3)
 492 : 
 493 :  MOT    Qc       LISTREEL : LT Debit total de condensation (Kg/s)
 494 :  MOT    Qcw      LISTREEL : LT Debit de condensation (THERMP) (Kg/s)
 495 :  MOT    Qc1      LISTREEL : LT Debit de condensation (TTIMP1) (Kg/s)
 496 :  MOT    Qc2      LISTREEL : LT Debit de condensation (TTIMP2) (Kg/s)
 497 :  MOT    Qc3      LISTREEL : LT Debit de condensation (TTIMP3) (Kg/s)
 498 :  MOT    Qc0      LISTREEL : LT Debit de condensation (ECHANP) (Kg/s)
 499 :  MOT    Econd    LISTREEL : LT energie extraite par condensation (J/s)
 500 :  MOT    Hcond    LISTREEL : LT enthalpie extraite par condensation(J/s)
 501 :  MOT    Econv    LISTREEL : LT energie de convection (J/s)
 502 : 
 503 :  MOT    Easpe    LISTREEL : LT energie extraite par l'aspersion (W=J/s) 
 504 :                                     (Si ASPER  VRAI)
 505 :  MOT    Haspe    LISTREEL : LT enthalpie extraite par l'aspersion(W=J/s)
 506 :  MOT    Qaspe    LISTREEL : LT debit aspersion  (kg/s)
 507 : 
 508 :  MOT    Remn     LISTREEL : LT energie interne specifique melange (J/m3)
 509 : 
 510 :  MOT    Rgpm     LISTREEL : LT constante gaz parfaits melange (J/kg/K)
 511 :  MOT    Cvm      LISTREEL : LT }capacite calorifique du melange (J/kg/K)
 512 :  MOT    Cpm      LISTREEL : LT }
 513 :  MOT    Gamm     LISTREEL : LT }et rapport
 514 : 
 515 :  MOT    PT       LISTREEL : LT pression thermodynamique
 516 :  MOT    Minj     LISTREEL : LT masses injectees (Kg)
 517 :  MOT    Mcond    LISTREEL : LT masses condensees (Kg)
 518 :  MOT    Mrest    LISTREEL : LT masses restantes  (Kg)
 519 : 
 520 :  MOT    guj      LISTREEL : LT vitesse a l'injection (breche 1) (m/s)
 521 :  MOT    Qj       LISTREEL : LT debit d'injection (breche 1) (Kg/s)
 522 :  MOT    Hj       LISTREEL : LT enthalpie a l'injection (breche 1) (J/kg)
 523 :  MOT    Ej       LISTREEL : LT energie a l'injection (breche 1) (J/kg)
 524 : 
 525 :  MOT    guj2     LISTREEL : LT 
 526 :  MOT    Qj2      LISTREEL : LT 
 527 :  MOT    Hj2      LISTREEL : LT    breche 2
 528 :  MOT    Ej2      LISTREEL : LT 
 529 : 
 530 :  MOT    TBP1     FLOTTANT : Temperature de paroi imposee (TTIMP1)
 531 :  MOT    TBP2     FLOTTANT : Temperature de paroi imposee (TTIMP2)
 532 :  MOT    TBP3     FLOTTANT : Temperature de paroi imposee (TTIMP3)
 533 :  MOT    TBP0     CHPOINT  : Temperature de paroi constante (ECHANP)
 534 :  MOT    KHEW     FLOTTANT : coefficient d'echange (convectif)
 535 :                             paroi/fluide  THERMP
 536 :  MOT    KHE1     FLOTTANT : coefficient d'echange (convectif) TIMP1
 537 :  MOT    KHE2     FLOTTANT : coefficient d'echange (convectif) TIMP2
 538 :  MOT    KHE3     FLOTTANT : coefficient d'echange (convectif) TIMP3
 539 :  MOT    KHE0     CHPOINT  : coefficient d'echange (convectif) ECHANP
 540 : 
 541 :  *--- RECOMBINEUR ---*
 542 :  
 543 :  MOT    QIN_H2   LISTREEL : LT debit d'hydrogene total a l'entree 
 544 :                             des recombineurs
 545 :  MOT    QIN_H2O  LISTREEL : LT debit de vapeur total a l'entree 
 546 :                             des recombineurs
 547 :  MOT    QIN_O2   LISTREEL : LT debit d'oxygene total a l'entree 
 548 :                             des recombineurs
 549 :  MOT    QIN_N2   LISTREEL : LT debit d'azote total a l'entree 
 550 :                             des recombineurs
 551 :  MOT    QOUT_H2  LISTREEL : LT debit d'hydrogene total a la sortie
 552 :                             des recombineurs
 553 :  MOT    QOUT_H2O LISTREEL : LT debit de vapeur total a la sortie
 554 :                             des recombineurs
 555 :  MOT    QOUT_O2  LISTREEL : LT debit d'oxygene total a la sortie
 556 :                             des recombineurs
 557 :  MOT    QOUT_N2  LISTREEL : LT debit d'azote total a la sortie
 558 :                             des recombineurs
 559 :  MOT    RECOMB   TABLE    : Table contenant differentes informations
 560 :                             a propos des recombineurs
 561 :  RXT.TIC.RECOMB. 1        : Table du premier recombineur
 562 :   (...)
 563 :  RXT.TIC.RECOMB. 1 . 'DEB' : LT contenant le debit du recombineur 1
 564 :  RXT.TIC.RECOMB. 1 . 'Uin' : LT contenant la vitesse a l'entree
 565 :                              du recombineur 1
 566 :  RXT.TIC.RECOMB. 1 . 'Uout': LT contenant la vitesse a la sortie
 567 :                              du recombineur 1
 568 :  RXT.TIC.RECOMB. 1 . 'Tin' : LT contenant la temperature a l'entree
 569 :                              du recombineur 1
 570 :  RXT.TIC.RECOMB. 1 . 'Tout': LT contenant la temperature a la sortie
 571 :                              du recombineur 1
 572 :  RXT.TIC.RECOMB. 1 . 'Tpla': LT contenant la temperature des plaques
 573 :                              du recombineur 1
 574 :  RXT.TIC.RECOMB. 1 . 'EFF' : LT contenant l'efficacite du recombineur 1
 575 :  RXT.TIC.RECOMB. 1 . 'Hin' : LT contenant l'enthalpie a l'entree
 576 :                              du recombineur 1
 577 :  RXT.TIC.RECOMB. 1 . 'Hout': LT contenant l'enthalpie a la sortie
 578 :                              du recombineur 1
 579 :  RXT.TIC.RECOMB. 1 . 'R?E' : LT contenant la densite de l'espece ? a 
 580 :                              l'entree du recombineur 1
 581 :  RXT.TIC.RECOMB. 1 . 'R?S' : LT contenant la densite de l'espece ? a 
 582 :                              la sortie du recombineur 1
 583 :  RXT.TIC.RECOMB. 1 . 'Y?E' : LT contenant la fraction massique de 
 584 :                              l'espece ? a l'entree du recombineur 1
 585 :  RXT.TIC.RECOMB. 1 . 'Y?S' : LT contenant la fraction massique de 
 586 :                              l'espece ? a la sortie du recombineur 1
 587 : 
 588 : * B/ Grandeurs multi-D
 589 : *---------------------
 590 :  MOT    UN       CHPOINT  : Champ de vitesse au pas de temps courant
 591 :  MOT    UNM      CHPOINT  : Champ de vitesse au pas de temps precedent
 592 :  MOT    ROG      CHPOINT  :
 593 :  MOT    PRES     CHPOINT  : Champ de pression
 594 :  MOT    TF       CHPOINT  : Champ de temperature au pas de temps courant
 595 :  MOT    TFNM     CHPOINT  : temperature au pas de temps precedent
 596 :  MOT    RHO      CHPOINT  : Champ de densite au pas de temps courant
 597 :  MOT    RHONM    CHPOINT  : Champ de densite au pas de temps precedent
 598 :  MOT    Mu       CHPOINT  : viscosite dynamique du melange
 599 :  MOT    NU       CHPOINT  : viscosite cinematique du melange
 600 :  MOT    NUEFF    CHPOINT  : viscosite cinematique effective
 601 : 
 602 : Aspersion
 603 :  MOT    VN       CHPOINT  : Champ de vitesse de la phase dispersee au
 604 :                              pas de temps courant
 605 :  MOT    TD       CHPOINT  : Champ de temperature de la phase dispersee
 606 :                              au pas de temps courant
 607 :  MOT    XD       CHPOINT  : Fraction volumique de la phase dispersee
 608 :                              au pas de temps courant
 609 :  MOT    DD       CHPOINT  : Diametre de la phase dispersee
 610 :                              au pas de temps courant
 611 : 
 612 : 
 613 :  MOT    RAIR     CHPOINT  : densite de l'air (Kg)
 614 :  MOT    RVP      CHPOINT  : densite de vapeur (Kg)
 615 :  MOT    RHE      CHPOINT  : densite de He     (Kg)
 616 :  MOT    RH2      CHPOINT  : densite de H2     (Kg)
 617 :  MOT    RO2      CHPOINT  : densite de O2     (Kg)
 618 :  MOT    RN2      CHPOINT  : densite de N2     (Kg)
 619 :  MOT    RCO2     CHPOINT  : densite de CO2    (Kg)
 620 :  MOT    RCO      CHPOINT  : densite de CO     (Kg)
 621 : 
 622 :  MOT    YVAP     FLOTTANT / CHPOINT : Fraction (F.M.) massique vapeur
 623 :  MOT    YHE      FLOTTANT / CHPOINT :
 624 :  MOT    YH2      FLOTTANT / CHPOINT : FLOTTANT si l'espece est absente
 625 :  MOT    YO2      FLOTTANT / CHPOINT :
 626 :  MOT    YN2      FLOTTANT / CHPOINT :
 627 :  MOT    YCO2     FLOTTANT / CHPOINT :
 628 :  MOT    YCO      FLOTTANT / CHPOINT :
 629 :  MOT    YAIR     FLOTTANT / CHPOINT :
 630 : 
 631 : 
 632 :  La TABLE GEO contient les informations geometriques en particulier
 633 :  les objets MMODEL.
 634 : 
 635 :      Indice               Objet
 636 :  Type      Valeur   Type        Commentaires
 637 : 
 638 :  MOT    epsi     FLOTTANT : tolerance pour confondre deux points dont
 639 :                                distance est inferieure a epsi
 640 :  MOT    $vtf     MMODEL   :
 641 :  MOT    $menvf   MMODEL   :
 642 :  MOT    $axe     MMODEL   : axe ou plan de symetrie
 643 :  MOT    $vtp     MMODEL   : paroi thermique
 644 :  MOT    $mtp1    MMODEL   : paroi temperature imposee (TIMP1)
 645 :  MOT    $mtp2    MMODEL   : paroi temperature imposee (TIMP2)
 646 :  MOT    $mtp3    MMODEL   : paroi temperature imposee (TIMP3)
 647 :  MOT    $mtp0    MMODEL   : paroi temperature constante (ECHANP)
 648 :  MOT    Pimp     MAILLAGE : maillage type POI1 contenant le point ou
 649 :                                la pression est imposee.
 650 : 
 651 : 
 652 : 

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