Afficher cette notice en
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