Afficher cette notice en
1 : $$$$ PRET NOTICE CHAT 11/09/12 21:17:43 7124 2 : DATE 11/09/12 3 : 4 : Operateur PRET Voir aussi : PRIM 5 : -------------- KONV 6 : 7 : 8 : I) ECOULEMENT COMPRESSIBLE 9 : _______________________ 10 : 11 : Dans le cadre de la modelisation d'un ecoulement compressible 12 : en discretisation volumes finis des Équations d'Euler, cet 13 : operateur permet de calculer les variables primitives aux 14 : interfaces (i.e. etats 'gauche' et 'droite' de chaque face), a 15 : partir des variables primitives aux centres. 16 : 17 : 18 : 19 : 20 : ----------------- 21 : | 1ere modele | 22 : ----------------- 23 : 24 : Gaz parfait (mono-espece); les chaleurs specifiques cp et cv 25 : sont independantes de la temperature. 26 : 27 : 28 : 29 : Cas a: premier ordre en espace, premier ordre en temps 30 : ______ 31 : 32 : MCHAM1 MCHAM2 MCHAM3 MCHAM4 = 'PRET' MCLE1 ENTI1 ENTI2 33 : MOD1 CHPO1 CHPO2 CHPO3 CHPO4 ; 34 : 35 : Commentaire : 36 : _____________ 37 : 38 : MCLE1 : MOT; 'PERFMONO'. 39 : 40 : ENTI1 : ENTIER; ordre en espace (=1). 41 : 42 : ENTI2 : ENTIER; ordre en temps (=1). 43 : 44 : MOD1 : Objet MODELE. 45 : 46 : CHPO1 : CHPOINT "CENTRE" contenant la masse volumique (en 47 : kg/m^3; une composante, 'SCAL'). 48 : 49 : CHPO2 : CHPOINT "CENTRE" contenant la vitesse (en m/s; 50 : deux composantes en 2D, 'UX ','UY '). 51 : 52 : CHPO3 : CHPOINT "CENTRE" contenant la pression du gaz 53 : (en Pa; une composante, 'SCAL'). 54 : 55 : CHPO4 : CHPOINT "CENTRE" contenant le "gamma" du gaz (une 56 : composante, 'SCAL'). 57 : 58 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 59 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 60 : (une composante, 'SCAL'; en kg/m^3) 61 : 62 : MCHAM2 : MCHAML contenant la vitesse (m/s) et les cosinus 63 : directeurs du repere locale (n,t) dans le repere 64 : global (x,y) (dans le cas 2D 6 composantes: 65 : * 'UN' = vitesse normale (SPG = 'FACEL') 66 : * 'UT' = vitesse tangentielle (SPG = 'FACEL') 67 : * 'NX' = n.x (SPG = 'FACE') 68 : * 'NY' = n.y (SPG = 'FACE') 69 : * 'TX' = t.x (SPG = 'FACE') 70 : * 'TY' = t.y (SPG = 'FACE')). 71 : 72 : MCHAM3 : MCHAML (SPG = "FACEL") contenant la pression du 73 : gaz (en Pa, une seule composante, 'SCAL'). 74 : 75 : MCHAM4 : MCHAML (SPG = "FACEL") contenant le "gamma" du 76 : gaz (une seule composante, 'SCAL'). 77 : 78 : 79 : 80 : Cas b: deuxieme ordre en espace, premier ou deuxieme ordre en temps 81 : ______ 82 : 83 : MCHAM1 MCHAM2 MCHAM3 MCHAM4 = 'PRET' MCLE1 ENTI1 ENTI2 84 : MOD1 CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 85 : CHPO7 CHPO8 CHPO9 CHPO10 (FLOT1) ; 86 : 87 : Commentaire : 88 : _____________ 89 : 90 : MCLE1 : MOT; 'PERFMONO' 91 : 92 : ENTI1 : ENTIER; ordre en espace (=2) 93 : 94 : ENTI2 : ENTIER; ordre en temps (=1 ou 2) 95 : 96 : MOD1 : Objet MODELE. 97 : 98 : CHPO1 : CHPOINT "CENTRE" contenant la masse volumique (en 99 : kg/m^3; une composante, 'SCAL'). 100 : 101 : CHPO2 : CHPOINT "CENTRE" contenant le gradient de la 102 : masse volumique (en kg/m^4; 2 composantes en 2D, 103 : 'P1DX', 'P1DY'). 104 : 105 : CHPO3 : CHPOINT "CENTRE" contenant le limiteur du 106 : gradient de la masse volumique (une seule 107 : composante 'P1 ') 108 : 109 : CHPO4 : CHPOINT "CENTRE" contenant la vitesse (en m/s; 110 : deux composantes en 2D, 'UX ','UY '). 111 : 112 : CHPO5 : CHPOINT "CENTRE" contenant le gradient de la 113 : vitesse (en s^-1; 4 composantes en 2D, 'P1DX', 114 : 'P1DY', 'P2DX', 'P2DY'). 115 : 116 : CHPO6 : CHPOINT "CENTRE" contenant le limiteur du 117 : gradient de la vitesse (2 composantes en 2D, 118 : 'P1', 'P2'). 119 : 120 : CHPO7 : CHPOINT "CENTRE" contenant la pression du gaz 121 : (en Pa; une composante, 'SCAL'). 122 : 123 : CHPO8 : CHPOINT "CENTRE" contenant le gradient de la 124 : pression (en Pa/m; 2 composantes en 2D, 'P1DX', 125 : 'P1DY'). 126 : 127 : CHPO9 : CHPOINT "CENTRE" contenant le limiteur du 128 : gradient de la pression (une composante, 'P1' ). 129 : 130 : CHPO10 : CHPOINT "CENTRE" contenant le "gamma" du gaz 131 : (une composante, 'SCAL'). 132 : 133 : FLOT1 : FLOTTANT (a specifier dans le cas ENTI2 = 2) qui 134 : contient l'increment en temps (en s) pour l'etape 135 : de prediction (valeur conseillee = Dt / 2); 136 : 137 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 138 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 139 : (une composante, 'SCAL'; en kg/m^3) 140 : 141 : MCHAM2 : MCHAML contenant la vitesse (m/s) et les cosinus 142 : directeurs du repere locale (n,t) dans le repere 143 : global (x,y) (dans le cas 2D 6 composantes: 144 : * 'UN' = vitesse normale (SPG = 'FACEL') 145 : * 'UT' = vitesse tangentielle (SPG = 'FACEL') 146 : * 'NX' = n.x (SPG = 'FACE') 147 : * 'NY' = n.y (SPG = 'FACE') 148 : * 'TX' = t.x (SPG = 'FACE') 149 : * 'TY' = t.y (SPG = 'FACE')). 150 : 151 : MCHAM3 : MCHAML (SPG = "FACEL") contenant la pression du 152 : gaz (en Pa, une seule composante, 'SCAL'). 153 : 154 : MCHAM4 : MCHAML (SPG = "FACEL") contenant le "gamma" du 155 : gaz (une seule composante, 'SCAL'). 156 : 157 : 158 : 159 : 160 : ----------------- 161 : | 2eme modele | 162 : ----------------- 163 : 164 : Melange de gaz parfaits (cp et cv independants de la temperature) 165 : 166 : 167 : 168 : Cas a: premier ordre en espace, premier ordre en temps 169 : ______ 170 : 171 : MCHAM1 MCHAM2 MCHAM3 MCHAM4 MCHAM5 = 'PRET' MCLE1 ENTI1 172 : ENTI2 MOD1 CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 ; 173 : 174 : Commentaire : 175 : _____________ 176 : 177 : MCLE1 : MOT; 'PERFMULT' 178 : 179 : ENTI1 : ENTIER; ordre en espace (=1) 180 : 181 : ENTI2 : ENTIER; ordre en temps (=1) 182 : 183 : MOD1 : Objet MODELE. 184 : 185 : CHPO1 : CHPOINT "CENTRE" contenant la masse volumique (en 186 : kg/m^3; une composante, 'SCAL'). 187 : 188 : CHPO2 : CHPOINT "CENTRE" contenant la vitesse (en m/s; 189 : deux composantes en 2D, 'UX ','UY '). 190 : 191 : CHPO3 : CHPOINT "CENTRE" contenant la pression du gaz 192 : (en Pa; une composante, 'SCAL'). 193 : 194 : CHPO4 : CHPOINT "CENTRE" contenant les fractions 195 : massiques du gaz (une ou plusieurs composantes) 196 : 197 : CHPO5 : CHPOINT "CENTRE" contenant le "gamma" du gaz 198 : (une composante, 'SCAL'). 199 : 200 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 201 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 202 : (une composante, 'SCAL'; en kg/m^3) 203 : 204 : MCHAM2 : MCHAML contenant la vitesse (m/s) et les cosinus 205 : directeurs du repere locale (n,t) dans le repere 206 : global (x,y) (dans le cas 2D 6 composantes: 207 : * 'UN' = vitesse normale (SPG = 'FACEL') 208 : * 'UT' = vitesse tangentielle (SPG = 'FACEL') 209 : * 'NX' = n.x (SPG = 'FACE') 210 : * 'NY' = n.y (SPG = 'FACE') 211 : * 'TX' = t.x (SPG = 'FACE') 212 : * 'TY' = t.y (SPG = 'FACE')). 213 : 214 : MCHAM3 : MCHAML (SPG = "FACEL") contenant la pression du 215 : gaz (en Pa, une seule composante, 'SCAL'). 216 : 217 : MCHAM4 : MCHAML (SPG = "FACEL") contenant les fractions 218 : massiques du gaz (meme nombre de composantes et 219 : memes nom que CHPO4). 220 : 221 : MCHAM5 : MCHAML (SPG = "FACEL") contenant le "gamma" du 222 : gaz (une seule composante, 'SCAL'). 223 : 224 : 225 : 226 : Cas b: deuxieme ordre en espace, premier ou deuxieme ordre en temps 227 : ______ 228 : 229 : MCHAM1 MCHAM2 MCHAM3 MCHAM4 MCHAM5 = 'PRET' MCLE1 ENTI1 ENTI2 230 : MOD1 TAB2 CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 231 : CHPO7 CHPO8 CHPO9 CHPO10 CHPO11 CHPO12 CHPO13 (FLOT1) ; 232 : 233 : Commentaire : 234 : _____________ 235 : 236 : MCLE1 : MOT; 'PERFMULT' 237 : 238 : ENTI1 : ENTIER; ordre en espace (=2) 239 : 240 : ENTI2 : ENTIER; ordre en temps (=1 ou 2) 241 : 242 : MOD1 : Objet MODELE. 243 : 244 : TAB2 : TABLE, qui contient les propriete du gaz 245 : (voir operateur PRIM) 246 : 247 : CHPO1 : CHPOINT "CENTRE" contenant la masse volumique (en 248 : kg/m^3; une composante, 'SCAL'). 249 : 250 : CHPO2 : CHPOINT "CENTRE" contenant le gradient de la 251 : masse volumique (en kg/m^4; 2 composantes en 2D, 252 : 'P1DX', 'P1DY'). 253 : 254 : CHPO3 : CHPOINT "CENTRE" contenant le limiteur du 255 : gradient de la masse volumique (une seule 256 : composante 'P1 ') 257 : 258 : CHPO4 : CHPOINT "CENTRE" contenant la vitesse (en m/s; 259 : deux composantes en 2D, 'UX ','UY '). 260 : 261 : CHPO5 : CHPOINT "CENTRE" contenant le gradient de la 262 : vitesse (en s^-1; 4 composantes en 2D, 'P1DX', 263 : 'P1DY', 'P2DX','P2DY'). 264 : 265 : CHPO6 : CHPOINT "CENTRE" contenant le limiteur du 266 : gradient de la vitesse (2 composantes en 2D, 267 : 'P1', 'P2'). 268 : 269 : CHPO7 : CHPOINT "CENTRE" contenant la pression du gaz 270 : (en Pa; une composante, 'SCAL'). 271 : 272 : CHPO8 : CHPOINT "CENTRE" contenant le gradient de la 273 : pression (en Pa/m; 2 composantes en 2D, 'P1DX', 274 : 'P1DY'). 275 : 276 : CHPO9 : CHPOINT "CENTRE" contenant le limiteur du 277 : gradient de la pression (une composante, 'P1' ). 278 : 279 : CHPO10 : CHPOINT "CENTRE" contenant les fractions 280 : massiques (n composantes (n >= 1), dont le nom 281 : est dans l'objet LISTMOTS en TAB2 . 'ESPEULE') 282 : 283 : CHPO11 : CHPOINT "CENTRE" contenant le gradient des 284 : fractions massiques (en m^-1, 2n composantes en 285 : 2D) 286 : 287 : CHPO12 : CHPOINT "CENTRE" contenant le limiteur du 288 : gradient des fractions massiques (n 289 : composantes). 290 : 291 : CHPO13 : CHPOINT "CENTRE" contenant le "gamma" du gaz 292 : (une composante, 'SCAL'). 293 : 294 : FLOT1 : FLOTTANT (a specifier dans le cas ENTI2 = 2) qui 295 : contient l'increment en temps (en s) pour l'etape 296 : de prediction (valeur conseillee = Dt / 2); 297 : 298 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 299 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 300 : (une composante, 'SCAL'; en kg/m^3) 301 : 302 : MCHAM2 : MCHAML contenant la vitesse (m/s) et les cosinus 303 : directeurs du repere locale (n,t) dans le repere 304 : global (x,y) (dans le cas 2D 6 composantes: 305 : * 'UN' = vitesse normale (SPG = 'FACEL') 306 : * 'UT' = vitesse tangentielle (SPG = 'FACEL') 307 : * 'NX' = n.x (SPG = 'FACE') 308 : * 'NY' = n.y (SPG = 'FACE') 309 : * 'TX' = t.x (SPG = 'FACE') 310 : * 'TY' = t.y (SPG = 'FACE')). 311 : 312 : MCHAM3 : MCHAML (SPG = "FACEL") contenant la pression du 313 : gaz (en Pa, une seule composante, 'SCAL'). 314 : 315 : MCHAM4 : MCHAML (SPG = "FACEL") contenant les fractions 316 : massiques du gaz (meme nombre de composantes et 317 : meme noms que CHPO10). 318 : 319 : MCHAM5 : MCHAML (SPG = "FACEL") contenant le "gamma" du 320 : gaz (une seule composante, 'SCAL'). 321 : 322 : 323 : 324 : ----------------- 325 : | 3eme modele | 326 : ----------------- 327 : 328 : Gaz parfait / melange de gaz parfaits (cv dependants de la temperature 329 : selon une fonction polynomiale) 330 : 331 : 332 : Cas a: premier ordre en espace, premier ordre en temps 333 : ______ 334 : 335 : MCHAM1 MCHAM2 MCHAM3 (MCHAM4) (MCHAM5) = 'PRET' MCLE1 ENTI1 336 : ENTI2 MOD1 TAB2 CHPO1 CHPO2 CHPO3 (CHPO4) (CHPO5) ; 337 : 338 : Commentaire : 339 : _____________ 340 : 341 : MCLE1 : MOT; 'PERFTEMP' 342 : 343 : ENTI1 : ENTIER; ordre en espace (=1) 344 : 345 : ENTI2 : ENTIER; ordre en temps (=1) 346 : 347 : MOD1 : Objet MODELE. 348 : 349 : TAB2 : TABLE, qui contient les propriete du gaz 350 : (voir operateur PRIM) 351 : 352 : CHPO1 : CHPOINT "CENTRE" contenant la masse volumique (en 353 : kg/m^3; une composante, 'SCAL'). 354 : 355 : CHPO2 : CHPOINT "CENTRE" contenant la vitesse (en m/s; 356 : deux composantes en 2D, 'UX ','UY '). 357 : 358 : CHPO3 : CHPOINT "CENTRE" contenant la pression du gaz 359 : (en Pa; une composante, 'SCAL'). 360 : 361 : (CHPO4) : CHPOINT "CENTRE" contenant les fractions 362 : massiques du gaz (uniquement dans le cas 363 : d'un melange; ses composantes sont aussi en 364 : TAB2 . 'ESPEULE') 365 : 366 : (CHPO5) : CHPOINT "CENTRE" contenant les scalaires 367 : passifs du gaz (uniquement dans le cas de 368 : transport de scalaires passifs; ses composantes 369 : sont en TAB2 . 'SCALPASS') 370 : 371 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 372 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 373 : (une composante, 'SCAL'; en kg/m^3) 374 : 375 : MCHAM2 : MCHAML contenant la vitesse (m/s) et les cosinus 376 : directeurs du repere locale (n,t) dans le repere 377 : global (x,y) (dans le cas 2D 6 composantes: 378 : * 'UN' = vitesse normale (SPG = 'FACEL') 379 : * 'UT' = vitesse tangentielle (SPG = 'FACEL') 380 : * 'NX' = n.x (SPG = 'FACE') 381 : * 'NY' = n.y (SPG = 'FACE') 382 : * 'TX' = t.x (SPG = 'FACE') 383 : * 'TY' = t.y (SPG = 'FACE')). 384 : 385 : MCHAM3 : MCHAML (SPG = "FACEL") contenant la pression du 386 : gaz (en Pa, une seule composante, 'SCAL'). 387 : 388 : (MCHAM4): MCHAML (SPG = "FACEL") contenant les fractions 389 : massiques du gaz (uniquement dans le cas d'un 390 : melange; meme composantes que CHPO4). 391 : 392 : (MCHAM5): MCHAML (SPG = "FACEL") contenant les scalaires 393 : passifs (uniquement dans le cas de transport de 394 : scalaires passifs). 395 : 396 : 397 : 398 : Cas b: deuxieme ordre en espace, premier ou deuxieme ordre en temps 399 : ______ 400 : 401 : Premier ordre en temps 402 : 403 : MCHAM1 MCHAM2 MCHAM3 (MCHAM4)(MCHAM5)= 'PRET' MCLE1 ENTI1 ENTI2 404 : MOD1 TAB2 405 : CHPO1 CHPO2 CHPO3 406 : CHPO4 CHPO5 CHPO6 407 : CHPO7 CHPO8 CHPO9 408 : (CHPO10 CHPO11 CHPO12) 409 : (CHPO13 CHPO14 CHPO15) ; 410 : 411 : Deuxieme ordre en temps 412 : 413 : MCHAM1 MCHAM2 MCHAM3 (MCHAM4)(MCHAM5)= 'PRET' MCLE1 ENTI1 ENTI2 414 : MOD1 TAB2 415 : CHPO1 CHPO2 CHPO3 416 : CHPO4 CHPO5 CHPO6 417 : CHPO7 CHPO8 CHPO9 418 : (CHPO10 CHPO11 CHPO12) 419 : (CHPO13 CHPO14 CHPO15) 420 : CHPO16 FLOT1 ; 421 : 422 : Commentaire : 423 : _____________ 424 : 425 : MCLE1 : MOT; 'PERFTEMP' 426 : 427 : ENTI1 : ENTIER; ordre en espace (=2) 428 : 429 : ENTI2 : ENTIER; ordre en temps (=1 ou 2) 430 : 431 : MOD1 : Objet MODELE. 432 : 433 : TAB2 : TABLE, qui contient les propriete du gaz 434 : (voir operateur PRIM) 435 : 436 : CHPO1 : CHPOINT "CENTRE" contenant la masse volumique (en 437 : kg/m^3; une composante, 'SCAL'). 438 : 439 : CHPO2 : CHPOINT "CENTRE" contenant le gradient de la 440 : masse volumique (en kg/m^4; 2 composantes en 2D, 441 : 'P1DX', 'P1DY'). 442 : 443 : CHPO3 : CHPOINT "CENTRE" contenant le limiteur du 444 : gradient de la masse volumique (une seule 445 : composante 'P1 ') 446 : 447 : CHPO4 : CHPOINT "CENTRE" contenant la vitesse (en m/s; 448 : deux composantes en 2D, 'UX ','UY '). 449 : 450 : CHPO5 : CHPOINT "CENTRE" contenant le gradient de la 451 : vitesse (en s^-1; 4 composantes en 2D, 'P1DX', 452 : 'P1DY', 'P2DX','P2DY'). 453 : 454 : CHPO6 : CHPOINT "CENTRE" contenant le limiteur du 455 : gradient de la vitesse (2 composantes en 2D, 456 : 'P1', 'P2')." 457 : 458 : CHPO7 : CHPOINT "CENTRE" contenant la pression du gaz 459 : (en Pa; une composante, 'SCAL'). 460 : 461 : CHPO8 : CHPOINT "CENTRE" contenant le gradient de la 462 : pression (en Pa/m; 2 composantes en 2D, 'P1DX', 463 : 'P1DY'). 464 : 465 : CHPO9 : CHPOINT "CENTRE" contenant le limiteur du 466 : gradient de la pression (une composante, 'P1' ). 467 : 468 : (CHPO10) : CHPOINT "CENTRE" contenant les fractions 469 : massiques (uniquement dans le cas d'un melange; 470 : n composantes, qui sont aussi en 471 : TAB2 . 'ESPEULE') 472 : 473 : (CHPO11) : CHPOINT "CENTRE" contenant le gradient des 474 : fractions massiques (en m^-1, uniquement dans le 475 : cas d'un melange; 2n composantes en 2D) 476 : 477 : (CHPO12) : CHPOINT "CENTRE" contenant le limiteur du 478 : gradient des fractions massiques (uniquement dans 479 : le cas d'un melange; n composantes). 480 : 481 : (CHPO13) : CHPOINT "CENTRE" contenant les scalaires 482 : passifs (uniquement dans le cas de transport de 483 : scalaires passifs; k composantes, qui sont aussi 484 : en TAB2 . 'SCALPASS') 485 : 486 : (CHPO14) : CHPOINT "CENTRE" contenant le gradient des 487 : scalaires passifs (uniquement dans le cas de 488 : transport de scalaires passifs; 2k composantes 489 : en 2D) 490 : 491 : (CHPO15) : CHPOINT "CENTRE" contenant le limiteur du gradient 492 : des scalaires passifs (uniquement dans le cas de 493 : transport de scalaires passifs; k composantes) 494 : 495 : CHPO16 : CHPOINT "CENTRE" contenant le "gamma" du gaz 496 : (une composante, 'SCAL'; uniquement dans le cas 497 : ENTI2 = 2). 498 : 499 : FLOT1 : FLOTTANT qui contient l'increment en temps (en s) 500 : pour l'etape de prediction (uniquement dans le cas 501 : ENTI2 = 2; valeur conseillee = Dt / 2); 502 : 503 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 504 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 505 : (une composante, 'SCAL'; en kg/m^3) 506 : 507 : MCHAM2 : MCHAML contenant la vitesse (m/s) et les cosinus 508 : directeurs du repere locale (n,t) dans le repere 509 : global (x,y) (dans le cas 2D 6 composantes: 510 : * 'UN' = vitesse normale (SPG = 'FACEL') 511 : * 'UT' = vitesse tangentielle (SPG = 'FACEL') 512 : * 'NX' = n.x (SPG = 'FACE') 513 : * 'NY' = n.y (SPG = 'FACE') 514 : * 'TX' = t.x (SPG = 'FACE') 515 : * 'TY' = t.y (SPG = 'FACE')). 516 : 517 : MCHAM3 : MCHAML (SPG = "FACEL") contenant la pression du 518 : gaz (en Pa, une seule composante, 'SCAL'). 519 : 520 : (MCHAM4): MCHAML (SPG = "FACEL") contenant les fractions 521 : massiques du gaz (uniquement dans le cas d'un 522 : melange; meme composantes que CHPO10). 523 : 524 : (MCHAM5): MCHAML (SPG = "FACEL") contenant les scalaires 525 : passifs (uniquement dans le cas de transport de 526 : scalaires passifs). 527 : 528 : Remarque: 529 : _________ 530 : 531 : Le SPG des MCHAMLs " 'DOMA' MOD1 'FACEL' " est un MAILLAGE de SEG3. 532 : Si la face n'est pas a la frontiere, le premier et le troisieme 533 : POINT sont les centres des ELTs s'appuyant sur la face 534 : (pour convention le premier point est le centre a gauche de la 535 : face et le troisieme est le centre a droite de la face); le 536 : deuxieme point est le centre de la face. Les MCHAMLs 'FACEL' 537 : sont definis sur le premier et le troisieme POINT et 538 : representent la reconstruction des variables a gauche et a 539 : droite de la face. 540 : Si la face est sur la frontiere, le premier POINT contient 541 : l'etat reconstruit et le troisieme contient son etat 'miroir' 542 : (i.e. on impose des conditions limites de type 'mur'). 543 : 544 : 545 : 546 : II) TRANSPORT DE SCALAIRES 547 : _______________________ 548 : 549 : Dans le cadre de la resolution en volumes finis de 550 : 551 : --> ---> 552 : d/dt S + u . grad (S) = 0 553 : --> 554 : avec div( u ) = 0 555 : 556 : cet operateur permet de calculer les variables S aux 557 : interfaces (i.e. etats 'gauche' et 'droite' de chaque face), a 558 : partir des valeurs aux centres. 559 : 560 : 561 : a) 1er ordre en espace 562 : 563 : CHAM1S = 'PRET' 'CLAUDEIS' 'FACE' ENTI1 MOD1 CHPO1 ; 564 : 565 : ENTI1 : ordre en espace (=1) 566 : 567 : MOD1 : Objet MODELE. 568 : 569 : CHPO1 : CHPOINT contenant les scalaires a transporter 570 : (SPG = 'DOMA' MOD1 'CENTRE') 571 : 572 : CHAM1S : MCHAML contenant les S scalaires a transporter 573 : (a gauche et a droite de chaque face, SPG = 'DOMA' MOD1 574 : 'CENTRE', autant de composantes que CHPO1) 575 : 576 : Remarque: 577 : _________ 578 : 579 : Le SPG des MCHAMLs " 'DOMA' MO1 'FACEL' " est un MAILLAGE de SEG3. 580 : Si la face n'est pas a la frontiere, le premier et le troisieme 581 : POINT sont les centres des ELTs s'appuyant sur la face 582 : (pour convention le premier point est le centre a gauche de la 583 : face et le troisieme est le centre a droite de la face); le 584 : deuxieme point est le centre de la face. Les MCHAMLs 'FACEL' 585 : sont definis sur le premier et le troisieme POINT et 586 : representent la reconstruction des variables a gauche et a 587 : droite de la face. 588 : Si la face est sur la frontiere, le premier POINT contient 589 : l'etat reconstruit et le troisieme contient son etat 'miroir'. 590 : 591 : 592 : 593 : III) ECOULEMENT COMPRESSIBLE, EQUATIONS BIFLUIDES 594 : ____________________________________________ 595 : 596 : Dans le cadre de la modelisation d'un ecoulement compressible 597 : en discretisation volumes finis des Équations du Flux Bifluide a 598 : une pression, cet operateur permet de calculer les variables 599 : primitives aux interfaces (i.e. etats 'gauche' et 'droite' de 600 : chaque face), a partir des variables primitives aux centres. 601 : 602 : Cas a: Premier ordre en espace et en temps 603 : ______ 604 : 605 : MCHAM8 MCHAM7 MCHAM6 MCHAM5 MCHAM4 MCHAM3 MCHAM2 MCHAM1 = 606 : 'PRET' MCLE1 ENTI1 ENTI2 MOD1 607 : CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 CHPO7 CHPO8 ; 608 : 609 : Commentaire: 610 : ____________ 611 : 612 : MCLE1 : MOT; 'TWOFLUID'. 613 : 614 : ENTI1 : ENTIER; ordre en espace (=1). 615 : 616 : ENTI2 : ENTIER; ordre en temps (=1). 617 : 618 : MOD1 : Objet MODELE. 619 : 620 : CHPO1 : CHPOINT contenant la taux de vide 621 : (une composante, 'SCAL'). 622 : 623 : CHPO2 : CHPOINT contenant la vitesse du gaz 624 : (m/s; deux composantes en 2D, 'UVX ','UVY ', trois 625 : composantes en 3D, 'UVX ','UVY ','UVZ '). 626 : 627 : CHPO3 : CHPOINT contenant la vitesse du liquide 628 : (m/s; deux composantes en 2D, 'ULX ','ULY ', trois 629 : composantes en 3D, 'ULX ','ULY ','ULZ '). 630 : 631 : CHPO4 : CHPOINT contenant la pression du gaz (Pa; une 632 : composante, 'SCAL'). 633 : 634 : CHPO5 : CHPOINT contenant la temperature du gaz 635 : (K; one component, 'SCAL'). 636 : 637 : CHPO6 : CHPOINT contenant la temperature du liquide 638 : (K; une composante, 'SCAL'). 639 : 640 : CHPO7 : CHPOINT contenant la densite du gaz 641 : (kg/m^3; une componsamte, 'SCAL'). 642 : 643 : CHPO8 : CHPOINT contenant la densite du liquide 644 : (kg/m^3; une composante, 'SCAL'). 645 : 646 : MCHAM1 : MCHAML contenant la taux de vide, qui a comme 647 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 648 : (une composante, 'SCAL') 649 : (voir operateur PRET) 650 : 651 : MCHAM2 : MCHAML contenant la vitesse du gaz et les cosinus 652 : directeurs du repere locale (n,t) dans le repere 653 : global (x,y) (dans le cas 2D 6 composantes: 654 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 655 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 656 : * 'NX' = n.x (SPG = 'FACE') 657 : * 'NY' = n.y (SPG = 'FACE') 658 : * 'TX' = t.x (SPG = 'FACE') 659 : * 'TY' = t.y (SPG = 'FACE')). 660 : (voir operateur PRET) 661 : 662 : MCHAM3 : MCHAML contenant la vitesse du liquide et les cosinus 663 : directeurs du repere locale (n,t) dans le repere 664 : global (x,y) (dans le cas 2D 6 composantes: 665 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 666 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 667 : * 'NX' = n.x (SPG = 'FACE') 668 : * 'NY' = n.y (SPG = 'FACE') 669 : * 'TX' = t.x (SPG = 'FACE') 670 : * 'TY' = t.y (SPG = 'FACE')). 671 : (voir operateur PRET) 672 : 673 : MCHAM4 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression 674 : (une seule composante, 'SCAL'). 675 : (voir operateur PRET) 676 : 677 : MCHAM5 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la 678 : temperature du gaz (une seule composante, 'SCAL'). 679 : (voir operateur PRET) 680 : 681 : MCHAM6 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la 682 : temperature du liquide (une seule composante, 'SCAL'). 683 : (voir operateur PRET) 684 : 685 : MCHAM7 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du 686 : gaz (une seule composante, 'SCAL'). 687 : (voir operateur PRET) 688 : 689 : MCHAM8 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du 690 : liquide (une seule composante, 'SCAL'). 691 : (voir operateur PRET) 692 : 693 : Cas b: deuxieme ordre en espace et en temps 694 : ______ 695 : 696 : MCHAM8 MCHAM7 MCHAM6 MCHAM5 MCHAM4 MCHAM3 MCHAM2 MCHAM1 = 697 : 'PRET' MCLE1 ENTI1 ENTI2 MOD1 698 : CHPO1 CHPO2 CHPO3 699 : CHPO4 CHPO5 CHPO6 700 : CHPO7 CHPO8 CHPO9 701 : CHPO10 CHPO11 CHPO12 702 : CHPO13 CHPO14 CHPO15 703 : CHPO16 CHPO17 CHPO18 704 : CHPO19 CHPO20 FLOT1 ; 705 : 'PRET' MCLE1 ENTI1 ENTI2 MOD1 706 : CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 CHPO7 CHPO8 ; 707 : 708 : Contents : 709 : _________ 710 : 711 : MCLE1 : MOT; 'TWOFLUID'. 712 : 713 : ENTI1 : ENTIER; ordre en espace (=1). 714 : 715 : ENTI2 : ENTIER; ordre en temps (=1). 716 : 717 : MOD1 : Objet MODELE. 718 : 719 : CHPO1 : CHPOINT contenant la taux de vide 720 : (une composante, 'SCAL'). 721 : 722 : CHPO2 : CHPOINT "CENTRE" contenant le gradient de la 723 : taux de vide (en 1/m; 2 composantes en 2D, 'P1DX', 724 : 'P1DY'). 725 : 726 : CHPO3 : CHPOINT "CENTRE" contenant le limiteur du 727 : gradient de la taux de vide (one component, 'P1 ') 728 : 729 : CHPO4 : CHPOINT "CENTRE" contenant la vitesse du gaz 730 : (m/s; deux composantes en 2D, 'UVX ','UVY ', trois 731 : composantes en 3D, 'UVX ','UVY ','UVZ '). 732 : 733 : CHPO5 : CHPOINT "CENTRE" contenant le gradient de la 734 : vitesse du gaz (en s^-1; 4 composantes en 2D, 'P1DX', 735 : 'P1DY', 'P2DX','P2DY'). 736 : 737 : CHPO6 : CHPOINT "CENTRE" contenant le limiteur du 738 : gradient de la vitesse du gaz (2 composantes en 2D, 739 : 'P1', 'P2'). 740 : 741 : CHPO7 : CHPOINT "CENTRE" contenant la vitesse du liquide 742 : (m/s; deux composantes en 2D, 'ULX ','ULY ', trois 743 : composantes en 3D, 'ULX ','ULY ','ULZ '). 744 : 745 : CHPO8 : CHPOINT "CENTRE" contenant le gradient de la 746 : vitesse du liquide (en s^-1; 4 composantes en 2D, 747 : 'P1DX', 'P1DY', 'P2DX','P2DY'). 748 : 749 : CHPO9 : CHPOINT "CENTRE" contenant le limiteur du 750 : gradient de la vitesse du liquid (2 composantes en 2D, 751 : 'P1', 'P2'). 752 : 753 : CHPO10 : CHPOINT "CENTRE" contenant la pression (Pa; une 754 : composante, 'SCAL'). 755 : 756 : CHPO11 : CHPOINT "CENTRE" contenant le gradient de la 757 : pression (en Pa/m; 2 composantes en 2D, 'P1DX', 758 : 'P1DY'). 759 : 760 : CHPO12 : CHPOINT "CENTRE" contenant le limiteur du 761 : gradient de la pression (une composante, 'P1' ). 762 : 763 : CHPO13 : CHPOINT contenant la temperature du gaz 764 : (K; une composante component, 'SCAL'). 765 : 766 : CHPO14 : CHPOINT "CENTRE" contenant le gradient de la 767 : temperature du gaz (en K/m; 2 composantes 768 : en 2D, 'P1DX', 'P1DY'). 769 : 770 : CHPO15 : CHPOINT "CENTRE" contenant le limiteur du 771 : gradient de la temperature du gaz 772 : (une composante, 'P1' ). 773 : 774 : CHPO16 : CHPOINT "CENTRE" contenant la temperature du liquide 775 : (K; une composante, 'SCAL'). 776 : 777 : CHPO17 : CHPOINT "CENTRE" contenant le gradient de la 778 : temperature du liquide (en K/m; 2 composantes 779 : en 2D, 'P1DX', 'P1DY'). 780 : 781 : CHPO18 : CHPOINT "CENTRE" contenant le limiteur du 782 : gradient de la temperature du liquide 783 : (une composante, 'P1' ). 784 : 785 : CHPO19 : CHPOINT "CENTRE" contenant la densite du gaz 786 : (kg/m^3; une componsamte, 'SCAL'). 787 : 788 : CHPO20 : CHPOINT "CENTRE" contenant la densite du liquide 789 : (kg/m^3; une composante, 'SCAL'). 790 : 791 : FLOT1 : FLOTTANT qui contient l'increment en temps (en s) 792 : pour l'etape de prediction (uniquement dans le cas 793 : ENTI2 = 2; valeur conseillee = Dt); 794 : 795 : MCHAM1 : MCHAML contenant la taux de vide, qui a comme 796 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 797 : (une composante, 'SCAL') 798 : (voir operateur PRET) 799 : 800 : MCHAM2 : MCHAML contenant la vitesse du gaz et les cosinus 801 : directeurs du repere locale (n,t) dans le repere 802 : global (x,y) (dans le cas 2D 6 composantes: 803 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 804 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 805 : * 'NX' = n.x (SPG = 'FACE') 806 : * 'NY' = n.y (SPG = 'FACE') 807 : * 'TX' = t.x (SPG = 'FACE') 808 : * 'TY' = t.y (SPG = 'FACE')). 809 : (voir operateur PRET) 810 : 811 : MCHAM3 : MCHAML contenant la vitesse du liquide et les cosinus 812 : directeurs du repere locale (n,t) dans le repere 813 : global (x,y) (dans le cas 2D 6 composantes: 814 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 815 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 816 : * 'NX' = n.x (SPG = 'FACE') 817 : * 'NY' = n.y (SPG = 'FACE') 818 : * 'TX' = t.x (SPG = 'FACE') 819 : * 'TY' = t.y (SPG = 'FACE')). 820 : (voir operateur PRET) 821 : 822 : MCHAM4 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression 823 : (une seule composante, 'SCAL'). 824 : (voir operateur PRET) 825 : 826 : MCHAM5 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la 827 : temperature du gaz (une seule composante, 'SCAL'). 828 : (voir operateur PRET) 829 : 830 : MCHAM6 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la 831 : temperature du liquide (une seule composante, 'SCAL'). 832 : (voir operateur PRET) 833 : 834 : MCHAM7 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du 835 : gaz (une seule composante, 'SCAL'). 836 : (voir operateur PRET) 837 : 838 : MCHAM8 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du 839 : liquide (une seule composante, 'SCAL'). 840 : (voir operateur PRET) 841 : 842 : 843 : IV) ECOULEMENT COMPRESSIBLE, 'REACTIVE DISCRETE EQUATION METHOD' 844 : ______________________________________ 845 : 846 : Dans le cadre de la modelisation d'un ecoulement compressible 847 : en discretisation volumes finis des Équations 848 : d'Euler avec REACTIVE DISCRETE EQUATION METHOD, 849 : cet operateur permet de calculer les variables 850 : primitives aux interfaces (i.e. etats 'gauche' et 'droite' de 851 : chaque face), a partir des variables primitives aux 852 : centres. Pour une description des variables veillez regarder 853 : la notice d'operateur 'PRIM'. 854 : 855 : 856 : 857 : a) deuxieme ordre en espace (pour appliquer premier 858 : ordre en espace veillez multiplier CHPGRA* et CHPLIA* 859 : par zero) 860 : 861 : CHAMAL1 CHAMAL2 CHAMRN1 CHAMRN2 CHAMV1 CHAMV2 CHAMP1 CHAMP2 = 862 : 'PRET' 'DEM' MOD1 TABPGAS 863 : CHPAL1 CHPGRA1 CHPLIA1 864 : CHPAL2 CHPGRA1 CHPLIA1 865 : CHPRN1 CHPGRR1 CHPLIR1 866 : CHPRN2 CHPGRR2 CHPLIR2 867 : CHPVN1 CHPGRV1 CHPLIV1 868 : CHPVN2 CHPGRV2 CHPLIV2 869 : CHPP1 CHPGRP1 CHPLIP1 870 : CHPP2 CHPGRP2 CHPLIP2 871 : EPS ; 872 : 873 : ARGUMENTS : 874 : 875 : 'DEM' : le mot cle pour "Reactive Discrete 876 : Equation Method" 877 : 878 : MOD1 : OBJET MODELE de TYPE 'EULER' 879 : 880 : TABPGAS : TABLE des proprietes des gaz 881 : (voir la notice de l'operateur 'PRIM') 882 : 883 : CHPAL1 CHPGRA1 CHPLIA1 : CHPOINTs "CENTRE": 884 : alpha1, son gradient et limiteur 885 : 886 : CHPAL2 CHPGRA2 CHPLIA2 : CHPOINTs "CENTRE": 887 : alpha2, son gradient et limiteur 888 : 889 : CHPRN1 CHPGRR1 CHPLIR1 : CHPOINTs "CENTRE": 890 : rho1 - la densite, son gradient et limiteur 891 : 892 : CHPRN2 CHPGRR2 CHPLIR2 : CHPOINTs "CENTRE": 893 : rho2 - la densite, son gradient et limiteur 894 : 895 : CHPVN1 CHPGRV1 CHPLIV1 : CHPOINTs "CENTRE": 896 : v1 - la vitesse, son gradient et limiteur 897 : 898 : CHPVN2 CHPGRV2 CHPLIV2 : CHPOINTs "CENTRE": 899 : v2 - la vitesse, son gradient et limiteur 900 : 901 : CHPVP1 CHPGRP1 CHPLIP1 : CHPOINTs "CENTRE": 902 : p1 - la pression, son gradient et limiteur 903 : 904 : CHPVP2 CHPGRP2 CHPLIP2 : CHPOINTs "CENTRE": 905 : p2 - la pression, son gradient et limiteur 906 : 907 : EPS : FLOTTANT; alpha_i < eps => no phase_i 908 : 909 : RESULTS : 910 : 911 : CHAMAL1 CHAMAL2 CHAMRN1 CHAMRN2 CHAMV1 CHAMV2 CHAMP1 912 : CHAMP2 : MCHAML "FACEL" qui contiennent alpha 913 : (1 et 2), densite (1 et 2), vitesse (1 et 2) et 914 : pression (1 et 2) a gauche et a droite de 915 : chaque interface. 916 : 917 : 918 : 919 : V) ECOULEMENT COMPRESSIBLE, GHOST FLUID METHOD FOR THE POOR (GFMP) 920 : _______________________________________________________________ 921 : 922 : Dans le cadre de la modelisation d'un ecoulement compressible 923 : en discretisation volumes finis des Équations d'Euler, methode 924 : GFMP, cet operateur permet de calculer les variables 925 : primitives aux interfaces (i.e. etats 'gauche' et 'droite' de 926 : chaque face), a partir des variables primitives aux centres. 927 : 928 : V.a) 2eme ordre en espace, une espece 929 : 930 : CHAMPH CHAMRN CHAMV CHAMP = 931 : 'PRET' 'GFMP' LMOT1 MOD1 932 : CHPPH CHPGRPH CHPLIP 933 : CHPRN CHPGRR CHPLIR 934 : CHPVN CHPGRV CHPLIV 935 : CHPP CHPGRP CHPLIP ; 936 : ou 937 : 938 : V.b) 2eme ordre en espace, multi-espece 939 : 940 : CHAMPH CHAMRN CHAMV CHAMP CHAMY CHAMAL = 941 : 'PRET' 'GFMP' LMOT1 MOD1 942 : CHPPH CHPGRPH CHPLIP 943 : CHPRN CHPGRR CHPLIR 944 : CHPVN CHPGRV CHPLIV 945 : CHPP CHPGRP CHPLIP 946 : CHPYN CHPGRY CHPLIY 947 : CHPAN CHPGRA CHPLIA ; 948 : 949 : 950 : INPUT 951 : 952 : 'GFMP' : OBJET MOT 953 : 954 : LMOT1 : OBJET LISTMOTS, noms des especes 955 : dans les equations d'Euler 956 : 957 : MOD1 : OBJET MODELE de TYPE EULER 958 : 959 : CHPPH CHPGRP CHPLIP : CHPOINTs "CENTRE": 960 : phi, gradient, limiteur 961 : 962 : CHPRN CHPGRR CHPLIR : CHPOINTs "CENTRE": 963 : rho, gradient, limiteur 964 : 965 : CHPVN CHPGRV CHPLIV : CHPOINTs "CENTRE": 966 : v, gradient, limiteur 967 : 968 : CHPVP CHPGRP CHPLIP : CHPOINTs "CENTRE": 969 : p gradient, limiteur 970 : 971 : CHPYN CHPGRY CHPLIY : CHPOINTs "CENTRE": 972 : y, gradient, limiteur 973 : 974 : CHPAN CHPGRA CHPLIA : CHPOINTs "CENTRE": 975 : alpha, gradient, limiteur 976 : 977 : OUTPUT 978 : 979 : CHAMPH CHAMRN CHAMV CHAMP CHAMY CHAMAL: 980 : MCHAML "FACEL" 981 : contenant phi, la densite, 982 : la velocite, la pression, les fractions 983 : massiques et volumiques 984 : a gauche et a droite des faces 985 : 986 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales