Télécharger prim.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ PRIM     NOTICE  CHAT      11/09/12    21:17:44     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Operateur PRIM                           Voir aussi :
   5 : 
   6 :        
   7 :     a)   EVOL2  = PRIM    EVOL1 ; 
   8 : 
   9 :     b1)  RCHPO1 RCHPO2 =  'PRIM' 'PERFMONO' CHPO1 CHPO2 CHPO3 CHPO4 ;
  10 : 
  11 :     b2)  RCHPO1 RCHPO2 RCHPO3 RCHPO4 RCHPO5 = 'PRIM' 'PERFMULT' TAB1   
  12 :           CHPO1 CHPO2  CHPO3  CHPO4  ;
  13 : 
  14 :     b3)  RCHPO1 RCHPO2 RCHPO3 RCHPO4 (RCHPO5) RCHPO6 = 
  15 :           'PRIM' 'PERFTEMP'  TAB1  CHPO1  CHPO2  CHPO3  CHPO4  
  16 :           (CHPO5) (CHPO6) ;
  17 :        ou
  18 : 
  19 :         RCHPO1 RCHPO2 RCHPO3 (RCHPO5) RCHPO6 = 'PRIM' 'PERFTEMP'
  20 :           TAB1 CHPO1 CHPO2  CHPO3 (CHPO5) (CHPO6) ;
  21 : 
  22 :     c)  RMAT1 =  'PRIM' 'CONSPRIM' MAIL1 LMOT1 LMOT2
  23 :            CHPO1 CHPO2 CHPO3 CHPO4 ;
  24 : 
  25 :     d)  RCHPO8 RCHPO7 RCHPO6 RCHPO5 RCHPO4 RCHPO3 RCHPO2 RCHPO1 =  
  26 :         'PRIM' 'TWOFLUID'
  27 :         CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 
  28 :         CHPO6 CHPO7 CHPO8 CHPO9 CHPO10 CHPO11;
  29 : 
  30 :     e)  RCHD1 RCHD2 RCHV1 RCHV2 RCHP1 RCHP2 RCHT1 RCHT2 = 
  31 :                                'PRIM' 'DEM' TABPGAS
  32 :         CHPAL1 CHPAL2 CHPARN1 CHPARN2 CHPAGN1 CHPAGN2 
  33 :         CHPARET1 CHPARET2 CHPTGUE1 CHPTGUE2 EPS ;
  34 : 
  35 :     f)  RCHPO0 RCHPO1 (RCHPO2)  = 'PRIM' 'GFMP' TAB1  
  36 :           CHPO0 CHPO1 CHPO2 CHPO3 (CHPO4 CHPO5) ;
  37 : 
  38 : 
  39 :  a) L'operateur PRIMITIVE calcule la primitive d'un objet EVOLUTION 
  40 :     pouvant representer des fonctions dont on connait la valeur pour 
  41 :     des abscisses croissantes. 
  42 :     La valeur de la primitive pour la premiere abscisse de chaque
  43 :     fonction est 0.
  44 :     
  45 :      ATTENTION : les valeurs sont rendues aux abscisses  de l'evolution
  46 :     initiale ( l'evolution  resultat ne represente pas la primitive 
  47 :     en tous points de l'intervale) 
  48 : 
  49 :  b) Dans le cadre de la modelisation d'un ecoulement compressible en 
  50 :     discretisation volumes finis (equations d'Euler ou Navier-Stokes), 
  51 :     cet operateur permet de calculer les variables primitives (i.e. 
  52 :     pression, vitesse, temperature, ...) a partir des variables 
  53 :     conservatives (i.e. masse volumique, quantite de mouvement, energie 
  54 :     volumique totale, masse volumique de chaque espece). 
  55 : 
  56 : 
  57 :  b1 -----------------
  58 :     |  1ere modele  |
  59 :     -----------------
  60 : 
  61 :     Gaz parfait (mono-espece); les chaleurs specifiques cp et cv 
  62 :     sont independantes de la temperature.
  63 : 
  64 :     RCHPO1 RCHPO2 =  'PRIM' MCLE1  CHPO1 CHPO2 CHPO3 CHPO4 ;
  65 : 
  66 :     Commentaire :
  67 :     _____________
  68 : 
  69 : 
  70 :     MCLE1   : MOT,  'PERFMONO'.
  71 :                     
  72 :     CHPO1   : CHPOINT  contenant la masse volumique (en kg/m^3; une 
  73 :               composante, 'SCAL').
  74 : 
  75 :     CHPO2   : CHPOINT  contenant les debits (en  kg/s/m^2; deux 
  76 :               composantes en 2D, 'UX  ','UY  ', trois composantes
  77 :               en 3D, 'UX  ','UY  ', 'UZ  ').
  78 :      
  79 :     CHPO3   : CHPOINT  contenant l'energie totale par unite de volume 
  80 :               (en J/m^3; une composante, 'SCAL').
  81 : 
  82 :     CHPO4   : CHPOINT  contenant le "gamma" du gaz (une composante, 
  83 :               'SCAL').
  84 :      
  85 :     RCHPO1  : CHPOINT  contenant la vitesse (m/s; deux  composantes 
  86 :               en 2D, 'UX  ','UY  ', trois  composantes
  87 :               en 3D, 'UX  ','UY  ', 'UZ  ').
  88 : 
  89 :     RCHPO2  : CHPOINT contenant la pression du gaz (Pa; une composante,
  90 :               'SCAL').
  91 :       
  92 : 
  93 :     Remarques :
  94 :     ___________
  95 : 
  96 :     1) On controle que
  97 :        * la pression est positive, 
  98 :        * 1 < gamma < 3
  99 :        * les CHPOINTs sont definis sur le meme support geometrique
 100 : 
 101 :     2) CHPO1, CHPO2, CHPO3  sont les variables conservatives des 
 102 :        Équations d'Euler.      
 103 : 
 104 : 
 105 :  b2 -----------------
 106 :     |  2eme modele  |
 107 :     -----------------
 108 : 
 109 :     Melange de gaz parfaits (cp et cv independants de la temperature)
 110 : 
 111 :     RCHPO1 RCHPO2 RCHPO3 RCHPO4 RCHPO5 = 'PRIM' MCLE1  TAB1 CHPO1  
 112 :            CHPO2  CHPO3  CHPO4  ;
 113 : 
 114 :     Commentaire :
 115 :     _____________
 116 : 
 117 : 
 118 :     MCLE1   : MOT, 'PERFMULT'.
 119 :                   
 120 :     TAB1    : TABLE qui contient :
 121 :               * les noms des especes qui apparaissent explicitement
 122 :                 dans les equations d'Euler en TAB1 . 'ESPEULE' 
 123 :                 (LISTMOTS);
 124 :               * le nom de l'espece qui n'y est pas dans  
 125 :                 TAB1 . 'ESPNEULE' (MOT);
 126 :               * les CP et les CV des gaz qui apparaissent en
 127 :                 TAB1 . 'ESPEULE' et en TAB1 . 'ESPNEULE' 
 128 :                 TAB1 . 'CP' (TABLE)
 129 :                 TAB1 . 'CV' (TABLE).
 130 : 
 131 :     CHPO1   : CHPOINT  contenant la masse volumique (en kg/m^3; une
 132 :               composante, 'SCAL').
 133 : 
 134 :     CHPO2   : CHPOINT  contenant les debits (en kg/s/m^2; deux 
 135 :               composantes en 2D, 'UX  ','UY  ', trois composantes
 136 :               en 3D, 'UX  ','UY  ', 'UZ  ').
 137 :      
 138 :     CHPO3   : CHPOINT  contenant l'energie totale par unite de volume 
 139 :               (en J/m^3; une composante, 'SCAL').
 140 : 
 141 :     CHPO4   : CHPOINT  contenant la masse volumique des especes qui sont
 142 :               explicitement "splitted" dans les equations  d'Euler 
 143 :               (en kg/m^3; leurs noms sont dans  TAB1 . 'ESPEULE').
 144 : 
 145 :     RCHPO1  : CHPOINT  contenant la vitesse (en m/s; deux composantes 
 146 :               en 2D, 'UX  ','UY  ', trois composantes
 147 :               en 3D, 'UX  ','UY  ', 'UZ  ').
 148 : 
 149 :     RCHPO2  : CHPOINT  contenant la pression du gaz (en Pa; une
 150 :               composante, 'SCAL').
 151 :       
 152 :     RCHPO3  : CHPOINT  contenant la temperature du gaz (en K; une
 153 :               composante, 'SCAL').
 154 :       
 155 :     RCHPO4  : CHPOINT  contenant les fractions massiques des 
 156 :               differentes especes (autant composantes que CHPO4;
 157 :               meme nom de composantes).
 158 : 
 159 :     RCHPO5  : CHPOINT  contenant les "gamma" du gaz (une composante,
 160 :               'SCAL').
 161 :       
 162 : 
 163 :     Remarques :
 164 :     ___________
 165 : 
 166 :     1) On controle que
 167 :        * la pression soit positive
 168 :        * 1 < gamma < 3
 169 :        * 0 < Y_i < 1
 170 :        * sum Y_i < 1
 171 :        * les CHPOINTs sont definis sur le meme support geometrique.
 172 : 
 173 :     2) CHPO1, CHPO2, CHPO3, CHPO4 sont les variables conservatives des
 174 :        Équations d'Euler.
 175 : 
 176 : 
 177 :     Exemple : 
 178 :     _________
 179 : 
 180 :     Pour montrer la structure de TAB1 on donne un exemple
 181 : 
 182 : 
 183 :     PGAZ = 'TABLE' ;
 184 : 
 185 :     *
 186 :     *** GAZ: H_2, O_2, H_2O, N_2
 187 :     *
 188 :     *   CP, CV en J/Kg/K @ T = 3000
 189 :     *
 190 : 
 191 :     *
 192 :     **** Especes qui sont dans les equations d'Euler
 193 :     *
 194 : 
 195 :     PGAZ . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
 196 : 
 197 :     *
 198 :     * Attention: 4 LETTRES    |____| |____| |____|
 199 :     *
 200 : 
 201 :     *
 202 :     **** Espece qui n'y est pas
 203 :     *
 204 :     
 205 :     
 206 :     PGAZ . 'ESPNEULE' = 'N2  ';
 207 :     
 208 :     * 4 LETTRES         |____| 
 209 : 
 210 :     
 211 :     PGAZ . 'CP' = 'TABLE' ;
 212 :     PGAZ . 'CP' . 'H2  '  = .18729066D+05 ;
 213 :     PGAZ . 'CP' . 'O2  '  = .11886820D+04 ;
 214 :     PGAZ . 'CP' . 'H2O '  = .31209047D+04 ;
 215 :     PGAZ . 'CP' . 'N2  '  = .12993995D+04 ;
 216 :     
 217 :     PGAZ . 'CV' = 'TABLE' ;
 218 :     PGAZ . 'CV' . 'H2  '  = .14571861D+05 ;
 219 :     PGAZ . 'CV' . 'O2  '  = .92885670D+03 ;
 220 :     PGAZ . 'CV' . 'H2O '  = .26589930D+04 ;
 221 :     PGAZ . 'CV' . 'N2  '  = .10024563D+04 ;
 222 :     
 223 : 
 224 :  b3  -----------------
 225 :     |  3eme modele  |
 226 :     -----------------
 227 : 
 228 :     Gaz parfait / melange de gaz parfaits avec cv dependants de la 
 229 :     temperature selon une loi polynomial.
 230 :     
 231 :     Cas multi-especes:
 232 :     RCHPO1 RCHPO2 RCHPO3 RCHPO4 (RCHPO5) RCHPO6 = 'PRIM' MCLE1  
 233 :           TAB1  CHPO1  CHPO2  CHPO3  CHPO4  (CHPO5) (CHPO6) ;
 234 : 
 235 :     Cas mono-espece:
 236 :     RCHPO1 RCHPO2 RCHPO3 (RCHPO5) RCHPO6 = 'PRIM' MCLE1  TAB1 
 237 :           CHPO1 CHPO2  CHPO3 (CHPO5) (CHPO6) ;
 238 : 
 239 :     Commentaire :
 240 :     _____________
 241 : 
 242 : 
 243 :     MCLE1   : MOT, 'PERFTEMP'.
 244 :                   
 245 :     TAB1    : TABLE qui contient :
 246 :               * le nom de l'espece qui n'est pas dans les Equations
 247 :                 d'Euler en
 248 :                 TAB1 . 'ESPNEULE' (MOT);
 249 :               * les noms des especes qui apparaissent explicitement
 250 :                 dans les equations d'Euler en 
 251 :                 TAB1 . 'ESPEULE' (LISTMOTS); 
 252 :                 dans le cas mono-espece  cet indice n'existe pas.
 253 :               * le degre de polynoms cv_i=cv_i(T), en
 254 :                 TAB1 . 'NORD' (ENTIER >= 0)
 255 :               * les proprietes de chaque gaz 'ESPI', dans la table
 256 :                 TAB1 . 'ESPI':
 257 :                 - TAB1 . 'ESPI' . 'A' (LISTREEL)
 258 :                   qui contient les (TAB1.'NORD')+1 coefficients des 
 259 :                   polynoms cv(T), (A0,A1,...); le cv(T) sont supposes 
 260 :                   etre en J/kg/K, donc les Ai doivent etre en unites
 261 :                   coherentes; 
 262 :                 - TAB1 . 'ESPI' . 'R' (J/kg/K, FLOTTANT) 
 263 :                   qui contient la constante du gaz parfait 
 264 :                 - TAB1 . 'ESPI' . 'H0K' (J/kg, FLOTTANT)
 265 :                   qui contient l'enthalpie de formation du gaz a 0K
 266 :                   (quantite numerique et pas physique)
 267 :               * TAB1 . 'SCALPASS' = si existante, noms des scalaires
 268 :                 passifs a transporter (LISTMOTS)
 269 : 
 270 :     CHPO1   : CHPOINT  contenant la masse volumique (en kg/m^3; une
 271 :               composante, 'SCAL').
 272 : 
 273 :     CHPO2   : CHPOINT  contenant les debits (en kg/s/m^2; deux 
 274 :               composantes en 2D, 'UX  ','UY  ', trois composantes
 275 :               en 3D, 'UX  ','UY  ', 'UZ  ').
 276 :      
 277 :     CHPO3   : CHPOINT  contenant l'energie totale par unite de volume 
 278 :               (en J/m^3; une composante, 'SCAL').
 279 : 
 280 :     CHPO4   : CHPOINT  contenant la masse volumique des especes qui 
 281 :               sont explicitement "splitted" dans les equations  
 282 :               d'Euler 
 283 :               (en kg/m^3; leurs noms sont dans  TAB1 . 'ESPEULE').
 284 : 
 285 :    (CHPO5)  : CHPOINT contenant les produits entre rho et les scalaires 
 286 :               passifs a transporter (leur noms sont en TAB1 . 'SCALPASS').
 287 : 
 288 :    (CHPO6)  : CHPOINT contenant une temperature de reference 
 289 :               (initialisation de la methode de Newton permettent de 
 290 :               calculer T a partir de l'energie)
 291 :               (en K, une composante, 'SCAL')
 292 : 
 293 :     RCHPO1  : CHPOINT  contenant la vitesse (en m/s; deux composantes 
 294 :               en 2D, 'UX  ','UY  ', trois composantes
 295 :               en 3D, 'UX  ','UY  ', 'UZ  ').
 296 : 
 297 :     RCHPO2  : CHPOINT  contenant la pression du gaz (en Pa; une
 298 :               composante, 'SCAL').
 299 :       
 300 :     RCHPO3  : CHPOINT contenant la temperature du gaz (en K; une
 301 :               composante, 'SCAL').
 302 :       
 303 :     RCHPO4  : CHPOINT  contenant les fractions massiques des 
 304 :               differentes especes (autant de composantes que CHPO4;
 305 :               meme nom de composantes).
 306 : 
 307 :    (RCHPO5) : CHPOINT contenant les scalaires passifs a transporter 
 308 :               (leur noms sont en TAB1 . 'SCALPASS').
 309 :       
 310 :     RCHPO6  : CHPOINT  contenant les "gamma" du gaz (une composante,
 311 :               'SCAL').
 312 : 
 313 : 
 314 :     Remarques :
 315 :     ___________
 316 : 
 317 :     1) On controle que
 318 :        * la pression et la temperature soient positives
 319 :        * 1 < gamma < 3
 320 :        * 0 < Y_i < 1
 321 :        * sum Y_i < 1
 322 :        * les CHPOINTs sont definis sur le meme support geometrique.
 323 : 
 324 :     2) CHPO1, CHPO2, CHPO3, CHPO4 CHPO5 sont les variables conservatives des
 325 :        Équations d'Euler.
 326 : 
 327 : 
 328 :     Exemple : 
 329 :     _________
 330 : 
 331 :     Pour montrer la structure de TAB1 on donne un exemple dans le cas 
 332 :     multi-especes:
 333 : 
 334 : 
 335 :     PGAZ = 'TABLE' ;
 336 : 
 337 :     *
 338 :     *** GAZ: H_2, O_2, H_2O, N_2
 339 :     *
 340 :     * 
 341 : 
 342 :     *
 343 :     **** Especes qui sont dans les equations d'Euler
 344 :     *
 345 : 
 346 :     PGAZ . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
 347 : 
 348 :     *
 349 :     * Attention: 4 LETTRES    |____| |____| |____|
 350 :     *
 351 : 
 352 :     *
 353 :     **** Espece qui n'y est pas
 354 :     *
 355 :     
 356 :     
 357 :     PGAZ . 'ESPNEULE' = 'N2  ';
 358 :     
 359 :     * 4 LETTRES         |____| 
 360 :     
 361 :     *
 362 :     **** Degre des polynoms
 363 :     *
 364 : 
 365 :     PGAZ . 'NORD' = 4 ;
 366 : 
 367 :     *
 368 :     **** Les tables qui contiennent les proprietes de chaque gaz
 369 :     *
 370 : 
 371 :     PGAZ .  'H2  ' = 'TABLE'  ;
 372 :     PGAZ .  'H2O ' = 'TABLE'  ;
 373 :     PGAZ .  'N2  ' = 'TABLE'  ;
 374 :     PGAZ .  'O2  ' = 'TABLE'  ;
 375 : 
 376 :     *
 377 :     **** R (J/Kg/K)
 378 :     *
 379 : 
 380 :     PGAZ .  'H2  ' . 'R' = 4130.0 ;
 381 :     PGAZ .  'H2O ' . 'R' = 461.4 ;
 382 :     PGAZ .  'N2  ' . 'R' = 296.8 ;
 383 :     PGAZ .  'O2  ' . 'R' = 259.8 ;
 384 : 
 385 :     *
 386 :     **** Regressions polynomials
 387 :     *
 388 : 
 389 :     PGAZ .  'H2  ' . 'A' = 'PROG'  9834.91866 0.54273926 0.000862203836
 390 :                                -2.37281455E-07 1.84701105E-11 ;
 391 :     PGAZ .  'H2O ' . 'A' = 'PROG' 1155.95625 0.768331151 -5.73129958E-05
 392 :                               -1.82753232E-08 2.44485692E-12 ;
 393 :     PGAZ .  'N2  ' . 'A' = 'PROG' 652.940766 0.288239099 -7.80442298E-05
 394 :                              8.78233606E-09 -3.05514485E-13 ;
 395 :     PGAZ .  'O2  ' . 'A' = 'PROG' 575.012333  0.350522002 
 396 :                       -0.000128294865    2.33636971E-08 -1.53304905E-12;
 397 :                               
 398 :     *
 399 :     **** "Enthalpies" (ou energies) de formations a 0K (J/Kg)
 400 :     *
 401 : 
 402 :     PGAZ .  'H2  ' . 'H0K' = -4.195D6 ;
 403 :     PGAZ .  'H2O ' . 'H0K' = -1.395D7 ;
 404 :     PGAZ .  'N2  ' . 'H0K' = -2.953D5 ;
 405 :     PGAZ .  'O2  ' . 'H0K' = -2.634D5 ;
 406 : 
 407 : 
 408 :  c) Dans le cadre de la modelisation d'un ecoulement compressible en 
 409 :     discretisation volumes finis (equations d'Euler ou Navier-Stokes), 
 410 :     cet operateur calcule la matrice jacobienne des variables 
 411 :     conservatives par rapport aux variables primitives (gaz parfait 
 412 :     mono-espece, les chaleurs specifiques cp et cv sont independantes 
 413 :     de la temperature).
 414 : 
 415 : 
 416 :  RMAT1 =  'KONV' 'VF' 'PERFMONO' 'CONSPRIM' MAIL1 LMOT1 LMOT2
 417 :            CHPO1 CHPO2 CHPO3 CHPO4 ;
 418 :    
 419 :  LMOT1   : objet de type LISTMOTS
 420 :            Noms des variables conservatives
 421 :            Il contient dans l'ordre suivant: le noms de la densite,
 422 :            de la qdm, de l'energie totale par unite de volume
 423 : 
 424 :  LMOT2   : objet de type LISTMOTS
 425 :            Noms des variables primitives
 426 :            Il contient dans l'ordre suivant: le noms de la densite,
 427 :            de la vitesse, de la pression
 428 : 
 429 :  MAIL1   : SPG des CHPOINTs
 430 : 
 431 :  CHPO1   : CHPOINT contenant la masse volumique
 432 :            (une seule composante, 'SCAL').
 433 : 
 434 :  CHPO2   : CHPOINT contenant la vitesse
 435 :            (deux/trois composantes
 436 :             'UX', 'UY', 'UZ')
 437 : 
 438 :  CHPO3   : CHPOINT contenant la pression du gaz
 439 :            (une seule composante, 
 440 :             'SCAL').
 441 : 
 442 :  CHPO4   : CHPOINT contenant le "gamma" du gaz
 443 :            (une seule composante, 
 444 :             'SCAL').
 445 : 
 446 :  RMAT1   : objet de type MATRIK
 447 :            (SPG =  MAIL1)
 448 :            (inconnues primales = variables primitives  = LMOT2)
 449 :            (inconnues duales = variables conservatives = LMOT1)
 450 :            Il contient le jacobien des variables conservatives 
 451 :            par rapport aux variables primitives.
 452 : 
 453 : 
 454 : d) Calcul des variables primitives (taux du vide, vitesse du gaz, 
 455 :    vitesse du liquide, pression du gaz, temperature du gaz, 
 456 :    temperature du liquide) si on connait les valeurs des 
 457 :    variables conservatives (masse volumique, quantite de mouvement, 
 458 :    energie volumique totale, masse volumique de chaque compossant)
 459 :    sur la modele bifluide a un pression et six equations,
 460 :    seulement valide por melanges d'eau et air
 461 : 
 462 :     -----------------
 463 : 
 464 :     L'air est considere un gaz parfait et l'eau est traitee
 465 :     avec l'equation de etat "Stiffened gas", les chaleurs specifiques 
 466 :     cp et cv sont independantes de la temperature.
 467 : 
 468 :     RCHPO8 RCHPO7 RCHPO6 RCHPO5 RCHPO4 RCHPO3 RCHPO2 RCHPO1 =  
 469 :     'PRIM' MCLE1  
 470 :     CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 CHPO7 CHPO8 CHPO9 CHPO10 CHPO11;
 471 : 
 472 :     Commentaire :
 473 :     _____________
 474 : 
 475 : 
 476 :     MCLE1   : MOT object,  'TWOFLUID'.
 477 :                     
 478 :     CHPO1   : CHPOINT contenant la masse volumique du gaz
 479 :               (kg/m^3; un composante 'SCAL').
 480 : 
 481 :     CHPO2   : CHPOINT contenant la masse volumique du liquide 
 482 :               (kg/m^3; un composante 'SCAL').
 483 : 
 484 :     CHPO3   : CHPOINT contenant la quantite de mouvement du gaz (kg/s/m^2; 
 485 :               deux composantes en 2D, 'UVX  ','UVY  ', trois composantes
 486 :               en 3D, 'UVX  ','UVY  ','UVZ  ').
 487 : 
 488 :     CHPO4   : CHPOINT contenant la quantite de mouvement du liquide 
 489 :               (kg/s/m^2; deux composantes en 2D, 'ULX  ','ULY  ', 
 490 :               trois composantes en 3D, 'ULX  ','ULY  ','ULZ  ').
 491 :      
 492 :     CHPO5   : CHPOINT contenant l'energie totale du gaz par unite de volume
 493 :               (J/m^3; une composante, 'SCAL').
 494 : 
 495 :     CHPO6   : CHPOINT contenant l'energie totale du gaz par unite de volume
 496 :               (J/m^3; une composante, 'SCAL').
 497 : 
 498 :     CHPO7   : CHPOINT contenant la taux de vide evalue le dernier pas du temp
 499 :               (une composante, 'SCAL').
 500 : 
 501 :     CHPO8   : CHPOINT contenant la temperature du gaz le dernier pas du temp
 502 :               (K; one component, 'SCAL').
 503 : 
 504 :     CHPO9   : CHPOINT contenant la temperature du liquide le dernier pas 
 505 :               du temp (K; une composante, 'SCAL').
 506 :     
 507 :     CHPO10  : CHPOINT contenant le parametre de correccion de pression
 508 :               selon le model du CATHARE (une composante, 'SCAL').
 509 : 
 510 :     CHPO11  : CHPOINT contenant le parametre de masse virtual correction
 511 :               selon le model du CATHARE (une composante, 'SCAL').
 512 : 
 513 :     RCHPO1  : CHPOINT contenant la taux de vide  
 514 :               (une composante, 'SCAL').
 515 :  
 516 :     RCHPO2  : CHPOINT contenant la vitesse du gaz 
 517 :               (m/s; deux composantes en 2D, 'UVX  ','UVY  ', trois 
 518 :               composantes en 3D, 'UVX  ','UVY  ','UVZ  ').
 519 : 
 520 :     RCHPO3  : CHPOINT  contenant la vitesse du liquide
 521 :               (m/s; deux composantes en 2D, 'ULX  ','ULY  ', trois 
 522 :               composantes en 3D, 'ULX  ','ULY  ','ULZ  ').
 523 : 
 524 :     RCHPO4  : CHPOINT contenant la pression du gaz (Pa; une
 525 :               composante, 'SCAL').
 526 : 
 527 :     RCHPO5  : CHPOINT contenant la temperature du gaz
 528 :               (K; one component, 'SCAL').
 529 : 
 530 :     RCHPO6  : CHPOINT contenant la temperature du liquide 
 531 :               (K; une composante, 'SCAL').
 532 : 
 533 :     RCHPO7  : CHPOINT contenant la densite du gaz
 534 :               (kg/m^3; une componsamte, 'SCAL').
 535 : 
 536 :     RCHPO8  : CHPOINT contenant la densite du liquide
 537 :               (kg/m^3; une composante, 'SCAL').
 538 : 
 539 :     RCHPO9  : CHPOINT contenant la correccion de pression selon la 
 540 :               modele du CATHARE (Pa; une composant, 'SCAL').
 541 :       
 542 : 
 543 :     Remarks :
 544 :     ___________
 545 : 
 546 :     1) On controle que
 547 :        * la pression soit positive
 548 :        * les CHPOINTs sont definis sur le meme support geometrique.
 549 : 
 550 :     2) CHPO1, CHPO2, CHPO3, CHPO4, CHPO5, CHPO6  sont les variables 
 551 :        conservatives de la modele bifluide a une pression et six 
 552 :        equationes
 553 : 
 554 : 
 555 : e) Evaluation des variables primitives (i.e. les masses volumiques
 556 :    des melanges brule et non-brules,  les vitesses
 557 :    des melanges  brule et non-brules, etc.) en
 558 :    utilisant les variables conservatives
 559 :    dans le cadre de la modelisation de la combustion en
 560 :    discretisation volumes finis + Reactive Discrete Equation Method.
 561 : 
 562 :    Ici le melange non-brule (en amont de la flamme) et
 563 :    le melange brule (en aval de la flamme) sont
 564 :    notes respectivement 1 et 2, et les variables
 565 :    primitives et conservatives ont attribuees a chaque melange.
 566 :  
 567 :    la fraction volumique \alpha_1 est egale a 1 dans
 568 :    le melange non-brule et 0 - dans le melange brule;
 569 :    \alpha_2 est egale a 0 dans le melange
 570 :    non-brule, et 1 - dans le melange brule.
 571 : 
 572 : 
 573 :         RCHD1 RCHD2 RCHV1 RCHV2 RCHP1 RCHP2 RCHT1 RCHT2 = 
 574 :                                'PRIM' 'DEM' TABPGAS
 575 :         CHPAL1 CHPAL2 CHPARN1 CHPARN2 CHPAGN1 CHPAGN2 
 576 :         CHPARET1 CHPARET2 CHPTGUE1 CHPTGUE2 EPS ;
 577 :  
 578 :   
 579 :   Arguments:
 580 :  ----------- 
 581 : 
 582 :       TABPGAS   : TABLE qui contient :
 583 :                       * 'SPECIES'  - les noms des especes 
 584 :                       * 'CHEM_COEF' - les coefficients dans
 585 :                            la reaction chimique consideree
 586 :                       * 'MASSFRA' - les fraction massiques 
 587 :                              initiale et finale de la
 588 :                              premiere espece dans 'SPECIES',
 589 :                              les fractions massiques finales
 590 :                              des autres especes ayants
 591 :                              de coefficients positifs dans 'CHEM_COEF',
 592 :                              les fractions massiques
 593 :                              initiales des especes ayants de coefficients
 594 :                              negatifs dans 'CHEM_COEF'                          
 595 :                       * 'RUNIV' = la constante universelle des gaz parfaits , 
 596 :                       * ESPi = la table qui contient les proprietes des 
 597 :                                des especes EPSi
 598 :                       * 'TMAX' la temperature maximale
 599 :                                pour expansion de cv, i.e. 
 600 :                         pour T>'TMAX', cv(T)=cv('TMAX') 
 601 :                       * ESPI . 'A' 
 602 :                         CV_i = \sum_{j=0,k} A_{i,j} T^j 
 603 :                       * ESPI . 'W' (Kg/mole)
 604 :                       * ESPI . 'H0K' 
 605 :                         e_{0,i} = h_{0,i} = h_{T_0,i} - {R_i * T_0 +
 606 :                                   {\sum_{j=0,k} A_{i,j} / (j+1) T_0^(j+1)}}; 
 607 : 
 608 :             CHPAL1  : CHPOINT qui contient la fraction volumique alpha_1 
 609 :                       de 1 (une composante, 'SCAL').
 610 : 
 611 :             CHPAL2  : CHPOINT qui contient la fraction volumique alpha_2
 612 :                       de 2 (une composante, 'SCAL').
 613 : 
 614 :             CHPARN1 : CHPOINT qui contient la alpha_1 * densite  de 1
 615 :                       (une composante, 'SCAL').
 616 : 
 617 :             CHPARN2 : CHPOINT qui contient la alpha_2 * density  de 2
 618 :                       (une composante, 'SCAL').
 619 : 
 620 :             CHPAGN1 : CHPOINT qui contient la alpha_1 * momentum de 1
 621 :                    (deux composantes, 'UX', 'UY').
 622 : 
 623 :             CHPAGN2 : CHPOINT qui contient la alpha_2 * momentum de 2
 624 :                    (deux composantes, 'UX', 'UY').
 625 : 
 626 :             CHPARET1: CHPOINT qui contient la alpha_1 * energie_totale 
 627 :                       de 1 (une composante, 'SCAL') .
 628 : 
 629 :             CHPARET2: CHPOINT qui contient la alpha_2 * energie_totale 
 630 :                       de 2 (une composante, 'SCAL').
 631 : 
 632 :             CHPTGUE1: CHPOINT qui contient la valeur 'guess' 
 633 :                       pour la temperature de 1 (une composante, 'SCAL').
 634 :              
 635 :             CHPTGUE2: CHPOINT qui contient la valeur 'guess'
 636 :                       pour la temperature de 2 (une composante, 'SCAL').
 637 : 
 638 : 
 639 :             EPS     : FLOTTANT tel que si ALPHA_i < EPS, on assume 
 640 :                       que lespece i n'existe plus.
 641 : 
 642 :     Results:
 643 :    ---------
 644 :             RCHD1   : CHPOINT qui contient la densite de 1
 645 :    
 646 :             RCHD2   : CHPOINT qui contient la densite de 2      
 647 : 
 648 :             RCHV1   : CHPOINT qui contient la vitesse de 1
 649 : 
 650 :             RCHV2   : CHPOINT qui contient la vitesse de 2
 651 :       
 652 :             RCHP1   : CHPOINT qui contient la pression de 1
 653 : 
 654 :             RCHP2   : CHPOINT qui contient la pression de 2
 655 :       
 656 :             RCHT1   : CHPOINT qui contient la temperature
 657 :                       de 1
 658 : 
 659 :             RCHT2   : CHPOINT qui contient la temperature
 660 :                       de 2
 661 : 
 662 : 
 663 :  f) Évaluation des variables primitives (i.e. masse volumique,
 664 :     vitesse, fractions massiques) en utilisant les variables 
 665 :     conservatives dans le cadre de la modelisation de la 
 666 :     propagation d'interfaces dans un milieu diphasique,
 667 :     via la methode GFMP ("the Ghost Fluid method for the poor", 
 668 :     equation d'etat "Stiffened gas").
 669 : 
 670 :     RCHPO0 RCHPO1 (RCHPO2)  = 'PRIM' 'GFMP' TAB1  CHPO0 CHPO1 CHPO2 
 671 :        CHPO3 (CHPO4 CHPO5) ;
 672 : 
 673 :     TAB1 : TABLE qui contient :
 674 :          * les noms des especes qui apparaissent explicitement
 675 :            dans les equations d'Euler en TAB1 . 'ESPEULE' 
 676 :            (LISTMOTS) ;
 677 :          * le nom de l'espece qui n'y est pas 
 678 :            (TAB1 . 'ESPNEULE' (MOT)) ;
 679 :          * les gamma et les pinf dans la zone phi < 0
 680 :            TAB1 . 'MGAM' (LISTREEL) ;
 681 :            TAB1 . 'MPIN' (LISTREEL) ;
 682 :          * les gamma et les pinf dans la zone phi > 0
 683 :            TAB1 . 'PGAM' (LISTREEL) ;
 684 :            TAB1 . 'PPIN' (LISTREEL) ;
 685 :            NB 
 686 :            La premiere valeur dans les objets LISTEEL
 687 :            TAB1 . 'MGAM', ... est celle de l'espece
 688 :            TAB1 . 'ESPNEULE'; les autres sont celles des
 689 :            especes TAB1 . 'ESPEULE'.
 690 :                    
 691 :            CHPO0   : CHPOINT  contenant la fonction phi
 692 :                      (une composante, 'SCAL').
 693 : 
 694 :            CHPO1   : CHPOINT  contenant la masse volumique
 695 :                      (une composante, 'SCAL').
 696 : 
 697 :            CHPO2   : CHPOINT  contenant les debits
 698 :                     (2 composantes en 2D, 'UX  ','UY  ');
 699 :     
 700 :            CHPO3   : CHPOINT  contenant l'energie totale per
 701 :                      unite de volume (RHO Et),
 702 :                      (une composante, 'SCAL').
 703 : 
 704 :            CHPO4   : CHPOINT  contenant les masses volumiques
 705 :                      des especes en (TAB1. 'ESPEULE')
 706 :                      (composante dans TAB1. 'ESPEULE').
 707 : 
 708 :            CHPO5   : CHPOINT  contenant les fractions volumiques
 709 :                      des especes en (TAB1. 'ESPEULE')
 710 :                      (composante dans TAB1. 'ESPEULE').
 711 : 
 712 : 
 713 :            RCHPO0  : CHPOINT  contenant la vitesse ;
 714 : 
 715 :            RCHPO1  : CHPOINT  contenant la pression du gaz ;
 716 :      
 717 :            RCHPO2  : CHPOINT  contenant les fractions massiques
 718 :                      des especes en (TAB1. 'ESPEULE')
 719 :                      (composante dans TAB1. 'ESPEULE').
 720 : 
 721 : 

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