Afficher cette notice en
1 : $$$$ KONV NOTICE CHAT 11/09/12 21:16:46 7124 2 : DATE 11/09/12 3 : 4 : Operateur KONV Voir aussi : NAVI 5 : -------------- 6 : 7 : 8 : OBJET : 9 : 10 : I Formulation Volume Finis (OPTI VF) : 11 : ______________________________________ 12 : 13 : Discretise l'operateur de convection par des schemas volumes finis 14 : Il calcule un champ point FACE qui represente 15 : le flux de chaque arete du maillage selon 3 schemas. 16 : On ecrit le flux dans l'increment KIZG. 17 : 18 : SYNTAXE 1 : KONV VNF VNC option; 19 : VNF = U calcule aux sommets (par ksof) 20 : VNC = U calcule aux centres (par knol) 21 : Option = 'upwind' | 'quick' | 'muscl' 22 : 23 : SYNTAXE 2 : en utilisant EQEX 24 : RV = EQEX TABDOM ALFA .. ITMA .. 25 : 'ZONE' MOD1 'OPER' 'KONV' VNF VNC 'INCO' 'TN' 26 : CLIM 'TN' TIMP ENTREE ... ; 27 : 28 : on met l'option dans klop : klop ind rv 'option' ; 29 : 30 : Erreur possible dans le calcul de DT si VN est petit changer EPSILON 31 : 32 : 33 : COMMENTAIRES : 34 : 35 : 36 : Le flux est evalue sur chaque maille par : FI = VNF TK LGR 37 : ou TK est la temperature sur la face K calculee de facons differentes selon 38 : les schemas a partir des temperatures sur les mailles connexes: TIM ,TI ET TIP. 39 : (Voir le rapport pour le choix de ces trois temperatures) 40 : 41 : 42 : 43 : -Schema UPWIND : 44 : 45 : TK = TI 46 : 47 : 48 : -Schema QUICK : 49 : 50 : TK = (TI+TIP)/2 + (TIM+TIP-2*TI)/8 51 : 52 : 53 : -Schema MUSCL : 54 : 55 : 56 : TK = TI + DELTA_I/2 57 : 58 : 59 : ou DELTA_I est la pente entre deux mailles : 60 : 61 : GAMMA_I = SIGNE (TIP-TI) 62 : DELTA_I = GAMMA_I * MAX (0,MIN(GAMMA_I*(TIP-TI),GAMMA_I*(TI-TIM))) 63 : 64 : 65 : 66 : 67 : 68 : Le bilan des flux precedents est fait dans l'operateur AVCT 69 : 70 : 71 : II Formulation Element Finis (OPTI EF ou EFM1) : 72 : ________________________________________________ 73 : 74 : Discretise l'operateur de convection en Element Finis. 75 : Suivant l'option l'operateur est traites sous forme 76 : conservative ou non conservative. 77 : 78 : SYNTAXE - EQEX Cf operateur EQEX 79 : ______________ 80 : 81 : 'OPER' 'KONV' ROC UN LAM 'INCO' 'TN' : 82 : 83 : 84 : 1/ Formulation non conservative 85 : 86 : roc ( u Grad T ) 87 : 88 : 2/ Formulation conservative 89 : 90 : Div ( roc u T ) 91 : 92 : 93 : Commentaires : 94 : ______________ 95 : 96 : roc capacite calorique (J/M**3/°C) 97 : FLOTTANT ou CHPOINT SCAL CENTRE ou CHPOINT SCAL SOMMET ou MOT 98 : lam conductivite thermique (W/M/°C) 99 : cette donnee est necessaire pour l'evaluation 100 : du Peclet de maille. En general c'est le coefficent 101 : de l'operateur LAPN. Si c'est operateur est absent 102 : mettre lam=0. 103 : FLOTTANT ou CHPOINT SCAL CENTRE ou CHPOINT SCAL SOMMET ou MOT 104 : un Champ de vitesse transportant 105 : CHPOINT VECT SOMMET ou MOT 106 : tn Champ de temperature 107 : CHPOINT SCAL SOMMET ou MOT 108 : 109 : Un coefficient de type MOT indique que l'operateur va chercher le 110 : coefficient dans la table INCO a l'indice MOT. 111 : 112 : 113 : Options : (EQEX) 114 : _________ 115 : 116 : La discretisation des termes de convection peut etre : 117 : 118 : centree OPTION CENTREE 119 : decentree OPTION SUPG 120 : decentree avec capture de choc OPTION SUPGCC Option par defaut 121 : Crank Nicholson generalise OPTION CNG 122 : (ordre 4 en temps) 123 : 124 : 125 : Formulation non conservative OPTION NOCONS Option par defaut 126 : Formulation conservative OPTION CONS 127 : 128 : Formulation EF OPTION EF 129 : 130 : 131 : 132 : III Discretisation des Equations d'Euler 133 : _________________________________________ 134 : 135 : IIIa : gaz parfait mono-constituent polytropique 136 : ________________________________________________ 137 : 138 : 139 : Discretisation en VF "cell-centered" des equations d'Euler pour un gaz 140 : parfait mono-constituent polytropique 141 : 142 : Inconnues: 143 : 144 : densite, quantite de mouvement (qdm), energie totale par unite de volume 145 : (variables conservatives) 146 : 147 : ou 148 : 149 : densite, vitesse, pression (variables primitives) 150 : 151 : On peut calculer: 152 : 153 : IIIa.1. Le flux numerique 154 : IIIa.2 Le residu 155 : IIIa.3 La matrice jacobienne du residu par rapport aux variables 156 : conservatives 157 : IIIa.4 La matrice jacobienne du residu par rapport aux variables 158 : primitives 159 : IIIa.5 La matrice de preconditionnent bas Mach par rapport aux 160 : variables conservatives 161 : IIIa.6 La matrice de preconditionnent bas Mach par rapport aux 162 : variables primitives 163 : IIIa.7 La contribution de quelque condition limite au residu et 164 : a la matrice jacobienne 165 : 166 : 167 : IIIa.1 et IIIa.2 Le flux numerique et le residu 168 : ________________ 169 : 170 : 171 : RCHPO1 RFLOT1 = 'KONV' 'VF' 'PERFMONO' MOT1 MOT2 MOD1 LMOT1 MCHAM1 172 : MCHAM2 MCHAM3 MCHAM4 (CHPO5 CHPO6) (MAIL1) ; 173 : 174 : ENTRÉES 175 : 176 : MOT1 : objet de type MOT 177 : Il vaut 'RESI' si on veut calculer le residu 178 : Il vaut 'FLUX' si on veut calculer le flux 179 : 180 : MOT2 : objet de type MOT 181 : Il indique la methode: 182 : 'GODUNOV' = solveur exacte 183 : 'VANLEER' = solveur de van Leer 184 : 'VLH' = solveur de van Leer Hanel 185 : 'HUSVL' = HUS (van Leer + Osher) 186 : 'HUSVLH' = HUS (van Leer Hanel + Osher) 187 : 'AUSMPLUS' = AUSM+ 188 : 'ROE' = solveur de Roe 189 : 'SS' = solveur choc-choc 190 : 'AUSMPLM' = AUSM+ low Mach 191 : 'RUSANOV' = solveur de Rusanov 192 : 'RUSANOLM' = solveur de Rusanov pour le bas-Mach 193 : 'CENTERED' = solveur centre 194 : 'ROELM' = solveur de Roe-Turkel pour le bas-Mach 195 : 'HLLC' = solveur HLLC 196 : 'HLLCLM' = solveur HLLC-Turkel pour le bas-Mach 197 : 'AUSMPUP' = solveur AUSM+up low Mach 198 : 199 : 200 : LMOT1 : objet de type LISTMOTS 201 : Noms de composantes du resultat (RCHPO1) 202 : Il contient dans l'ordre suivant: le noms de la densite, 203 : de la qdm, de l'energie totale par unite de volume 204 : 205 : MOD1 : objet MODELE. 206 : 207 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 208 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 209 : (une composante, 'SCAL') 210 : (voir operateur PRET) 211 : 212 : MCHAM2 : MCHAML contenant la vitesse et les cosinus 213 : directeurs du repere locale (n,t) dans le repere 214 : global (x,y) (dans le cas 2D 6 composantes: 215 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 216 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 217 : * 'NX' = n.x (SPG = 'FACE') 218 : * 'NY' = n.y (SPG = 'FACE') 219 : * 'TX' = t.x (SPG = 'FACE') 220 : * 'TY' = t.y (SPG = 'FACE')). 221 : (voir operateur PRET) 222 : 223 : MCHAM3 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression du 224 : gaz (une seule composante, 'SCAL'). 225 : (voir operateur PRET) 226 : 227 : MCHAM4 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant le "gamma" du 228 : gaz (une seule composante, 'SCAL'). 229 : (voir operateur PRET) 230 : 231 : CHPO5 : CHPOINT contenant la vitesse de cut-off 232 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 233 : 'SCAL'). 234 : A donner dans le cas bas Mach 235 : 236 : CHPO6 : CHPOINT contenant la deuxieme vitesse de cut-off 237 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 238 : 'SCAL'). 239 : A donner dans le cas bas Mach 240 : 241 : (MAIL1) : MAILLAGE de POI1 ou, si appartenant a ('DOMA' MOD1 'FACE'), 242 : on ne calcule pas la contribution au flux ou au residu. 243 : 244 : SORTIES 245 : 246 : RCHPO1 : objet de type CHPOINT (composantes = LMOT1) 247 : Residu si MOT2 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE') 248 : Flux si MOT2 = 'FLUX' (SPG = 'DOMA' MOD1 'FACE') 249 : 250 : RFLOT1 : objet de type FLOTTANT 251 : Il est le temps caracteristique associe a l'onde la plus 252 : rapide (meme dans le cas bas Mach, on considere le systeme 253 : non preconditionne) 254 : 255 : Remarque 256 : -------- 257 : 258 : RCHPO1 est egal a: 259 : * la derive temporelle des inconnues si l'option 'RESI' est utilisee 260 : * la projection du flux sur ('DOMA' MOD1 'XXNORMAF') si l'option 261 : 'FLUX' est utilisee 262 : 263 : 264 : IIIa.3 La matrice jacobienne du residu par rapport aux variables 265 : _____ conservatives 266 : 267 : 268 : RMAT1 = 'KONV' 'VF' 'PERFMONO' 'JACOCONS' MOD1 LMOT1 (MAIL1) MOT1 269 : CHPO1 CHPO2 CHPO3 CHPO4 (CHPO5 CHPO6) ; 270 : 271 : ENTRÉES 272 : 273 : MOT1 : objet de type MOT 274 : 'VLH' : jacobien du residu pour la methode VLH 275 : 'AUSMPLUS' : jacobien du residu pour la methode AUSM+ 276 : 'AUSMPLM' : jacobien du residu pour la methode AUSM+ low Mach 277 : 'RUSANOLM' : jacobien du residu pour la methode de Rusanov low Mach 278 : 'CENTERED' : jacobien du residu pour la methode centree 279 : 280 : LMOT1 : objet de type LISTMOTS 281 : Noms des variables conservatives 282 : Il contient dans l'ordre suivant: le noms de la densite, 283 : de la qdm, de l'energie totale par unite de volume 284 : 285 : MOD1 : objet MODELE. 286 : 287 : CHPO1 : CHPOINT contenant la masse volumique 288 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 289 : 'SCAL'). 290 : 291 : CHPO2 : CHPOINT contenant la vitesse 292 : (SPG = 'DOMA' DOM1 'CENTRE', deux/trois composantes 293 : 'UX', 'UY', 'UZ') 294 : 295 : CHPO3 : CHPOINT contenant la pression du gaz 296 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 297 : 'SCAL'). 298 : 299 : CHPO4 : CHPOINT contenant le "gamma" du gaz 300 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 301 : 'SCAL'). 302 : 303 : CHPO5 : CHPOINT contenant la vitesse de cut-off 304 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 305 : 'SCAL'). 306 : A donner dans le cas bas Mach 307 : 308 : CHPO6 : CHPOINT contenant la deuxieme vitesse de cut-off 309 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 310 : 'SCAL'). 311 : A donner dans le cas bas Mach 312 : 313 : (MAIL1) : MAILLAGE de POI1 ou, si appartenant a ('DOMA' DOM1 'FACE'), 314 : on ne calcule pas la contribution au jacobien. 315 : 316 : SORTIES 317 : 318 : RMAT1 : objet de type MATRIK 319 : (SPG = 'DOMA' DOM1 'CENTRE') 320 : (inconnues primales = inconnues duales = LMOT1) 321 : Il contient le jacobien du residu par rapport aux variables 322 : conservatives. 323 : 324 : 325 : IIIa.4 La matrice jacobienne du residu par rapport aux variables 326 : ______ primitives 327 : 328 : RMAT1 = 'KONV' 'VF' 'PERFMONO' 'JACOPRIM' MOD1 LMOT1 LMOT2 (MAIL1) MOT1 329 : CHPO1 CHPO2 CHPO3 CHPO4 (CHPO5 CHPO6) ; 330 : 331 : ENTREES 332 : 333 : MOT1 : objet de type MOT 334 : 'VLH' : jacobien du residu pour la methode VLH 335 : 'AUSMPLUS' : jacobien du residu pour la methode AUSM+ 336 : 'AUSMPLM' : jacobien du residu pour la methode AUSM+ low Mach 337 : 338 : LMOT1 : objet de type LISTMOTS 339 : Noms des variables conservatives 340 : Il contient dans l'ordre suivant: le noms de la densite, 341 : de la qdm, de l'energie totale par unite de volume 342 : 343 : LMOT2 : objet de type LISTMOTS 344 : Noms des variables primitives 345 : Il contient dans l'ordre suivant: le noms de la densite, 346 : de la vitesse, de la pression 347 : 348 : MOD1 : objet MODELE. 349 : 350 : CHPO1 : CHPOINT contenant la masse volumique 351 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 352 : 'SCAL'). 353 : 354 : CHPO2 : CHPOINT contenant la vitesse 355 : (SPG = 'DOMA' MOD1 'CENTRE', deux/trois composantes 356 : 'UX', 'UY', 'UZ') 357 : 358 : CHPO3 : CHPOINT contenant la pression du gaz 359 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 360 : 'SCAL'). 361 : 362 : CHPO4 : CHPOINT contenant le "gamma" du gaz 363 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 364 : 'SCAL'). 365 : 366 : CHPO5 : CHPOINT contenant la vitesse de cut-off 367 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 368 : 'SCAL'). 369 : A donner dans le cas bas Mach 370 : 371 : CHPO6 : CHPOINT contenant la deuxieme vitesse de cut-off 372 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 373 : 'SCAL'). 374 : A donner dans le cas bas Mach 375 : 376 : (MAIL1) : MAILLAGE de POI1 ou, si appartenant a ('DOMA' MOD1 'FACE'), 377 : on ne calcule pas la contribution au jacobien. 378 : 379 : SORTIES 380 : 381 : RMAT1 : objet de type MATRIK 382 : (SPG = 'DOMA' MOD1 'CENTRE') 383 : (inconnues primales = variables primitives = LMOT2) 384 : (inconnues duales = variables conservatives = LMOT1) 385 : Il contient le jacobien du residu par rapport aux variables 386 : primitives. 387 : 388 : 389 : IIIa.5 La matrice de preconditionnement bas Mach par rapport aux 390 : ______ variables conservatives (divise par le pas de temps locale; 391 : le pas de temps est calcule pour le systeme preconditionnne 392 : (intervalle de temps pendant lequel le signale le plus rapide 393 : traverse une cellule reguliere)) 394 : 395 : RMAT1 = 'KONV' 'VF' 'PERFMONO' 'GAMMCONS' MAIL1 LMOT1 396 : CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 ; 397 : 398 : LMOT1 : objet de type LISTMOTS 399 : Noms des variables conservatives 400 : Il contient dans l'ordre suivant: le noms de la densite, 401 : de la qdm, de l'energie totale par unite de volume 402 : 403 : MAIL1 : SPG des CHPOINTs 404 : 405 : CHPO0 : CHPOINT contenant le diametre des elts 406 : (une seule composante, 'SCAL'). 407 : 408 : CHPO1 : CHPOINT contenant la masse volumique 409 : (une seule composante, 'SCAL'). 410 : 411 : CHPO2 : CHPOINT contenant la vitesse 412 : (deux/trois composantes 413 : 'UX', 'UY', 'UZ') 414 : 415 : CHPO3 : CHPOINT contenant la pression du gaz 416 : (une seule composante, 417 : 'SCAL'). 418 : 419 : CHPO4 : CHPOINT contenant le "gamma" du gaz 420 : (une seule composante, 421 : 'SCAL'). 422 : 423 : CHPO5 : CHPOINT contenant la vitesse de cut-off 424 : (SPG = 'CENTRE', une seule composante, 425 : 'SCAL'). 426 : A donner dans le cas bas Mach 427 : 428 : CHPO6 : CHPOINT contenant la deuxieme vitesse de cut-off 429 : (SPG = 'CENTRE', une seule composante, 430 : 'SCAL'). 431 : A donner dans le cas bas Mach 432 : 433 : RMAT1 : objet de type MATRIK 434 : (SPG = MAIL1) 435 : (inconnues primales = variables duales = LMOT1) 436 : 437 : IIIa.6 La matrice de preconditionnement bas Mach par rapport aux 438 : ______ variables primitives (divise par le pas de temps locale; 439 : le pas de temps est calcule pour le systeme preconditionnne 440 : (intervalle de temps pendant lequel le signale le plus rapide 441 : traverse une cellule reguliere)) 442 : 443 : RMAT1 = 'KONV' 'VF' 'PERFMONO' 'GAMMPRIM' MAIL1 LMOT1 LMOT2 444 : CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 ; 445 : 446 : LMOT1 : objet de type LISTMOTS 447 : Noms des variables conservatives 448 : Il contient dans l'ordre suivant: le noms de la densite, 449 : de la qdm, de l'energie totale par unite de volume 450 : 451 : LMOT2 : objet de type LISTMOTS 452 : Noms des variables primitives 453 : Il contient dans l'ordre suivant: le noms de la densite, 454 : de la vitesse, de la pression 455 : 456 : MAIL1 : SPG des CHPOINTs 457 : 458 : CHPO0 : CHPOINT contenant le diametre des elts 459 : (une seule composante, 'SCAL'). 460 : 461 : CHPO1 : CHPOINT contenant la masse volumique 462 : (une seule composante, 'SCAL'). 463 : 464 : CHPO2 : CHPOINT contenant la vitesse 465 : (deux/trois composantes 466 : 'UX', 'UY', 'UZ') 467 : 468 : CHPO3 : CHPOINT contenant la pression du gaz 469 : (une seule composante, 470 : 'SCAL'). 471 : 472 : CHPO4 : CHPOINT contenant le "gamma" du gaz 473 : (une seule composante, 474 : 'SCAL'). 475 : 476 : CHPO5 : CHPOINT contenant la vitesse de cut-off 477 : (SPG = 'CENTRE', une seule composante, 478 : 'SCAL'). 479 : A donner dans le cas bas Mach 480 : 481 : CHPO6 : CHPOINT contenant la deuxieme vitesse de cut-off 482 : (SPG = 'CENTRE', une seule composante, 483 : 'SCAL'). 484 : A donner dans le cas bas Mach 485 : 486 : RMAT1 : objet de type MATRIK 487 : (SPG = 'CENTRE') 488 : (inconnues primales = variables primitives = LMOT2) 489 : (inconnues duales = variables conservatives = LMOT1) 490 : 491 : 492 : IIIa.7 La contribution de quelque condition limite au residu et 493 : ------ a la matrice jacobienne 494 : 495 : RCHPLI RCHPRE = 'KONV' 'VF' 'PERFMONO' 'CLIM' 496 : 'RESI' $MOD1 $MOD2 LMOTC LMOTP CHPRN CHPVN CHPPN 497 : CHPGN MOT1 CHPLI ; 498 : 499 : ou 500 : 501 : RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' 502 : $MOD1 $MOD2 LMOTC LMOTP CHPRN CHPVN CHPPN 503 : CHPGN MOT1 CHPLI ; 504 : 505 : RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' 506 : $MOD1 $MOD2 LMOTC LMOTP CHPRN CHPVN CHPPN 507 : CHPGN MOT1 CHPLI ; 508 : 509 : $MOD1 : l'objet modele du domaine total 510 : 511 : $MOD2 : l'objet modele du domaine du bord 512 : 513 : LMOTC : LISTMOTS, noms des variables conservatives 514 : 515 : LMOTP : LISTMOTS, noms des variables primitives 516 : 517 : CHPRN : densite (SPG = 'DOMA' $MOD1 'CENTRE', une seule 518 : composante, 'SCAL') 519 : 520 : CHPVN : vitesse (SPG = 'DOMA' $MOD1 'CENTRE', 521 : composantes: 'UX', 'UY', ('UZ')) 522 : 523 : CHPPN : pression (SPG = 'DOMA' $MOD1 'CENTRE', une seule 524 : composante, 'SCAL') 525 : 526 : CHPGN : gamma (SPG = 'DOMA' $MOD1 'CENTRE', une seule 527 : composante, 'SCAL') 528 : 529 : MOT1 : MOT, type de condition limite: 530 : 'INRI' : entree subsonique. On utilise les invariants de 531 : Riemann. On donne densite 'RN', vitesse 'UX' 'UY' 532 : ('UZ'), pression 'PN'. 533 : 'INSU' : entree subsonique. On donne l'enthalpie totale 'HT' 534 : par unite de masse, l'entropie 'S' (pression divisee 535 : par la densite a la puissance gamma). On impose que 536 : la vitesse tangentielle est nulle. On recupere la 537 : pression a l'interieur. La contribution au residu et 538 : a la matrice jacobienne sont calculee via la methode 539 : 'AUSMPLUS'. 540 : 'OUTRI': sortie subsonique. On utilise les invariants de 541 : Riemann. On donne densite 'RN', vitesse 'UX' 'UY' 542 : ('UZ'), pression 'PN'. 543 : 'OUTP' : sortie subsonique. On donne la pression 'PN', 544 : on recupere la vitesse et la densite a l'interieur. 545 : La contribution au residu et a la matrice jacobienne 546 : sont calculee via la methode 'AUSMPLUS'. 547 : 'INSS' : entree supersonique. On donne densite 'RN', vitesse 548 : 'UX' 'UY' ('UX'), pression 'PN'. 549 : 'OUTSS': sortie supersonique. On donne un CHPOINT vide. 550 : 'INJE' : condition limite d'injection compressible. 551 : On donne le flux de masse 'MOME' et 'RT' (temperature fois 552 : la constante de gaz); on impose que la vitesse 553 : tangentielle est nulle, on recupere la pression. 554 : 'INJELM': condition limite d'injection faiblement compressible. 555 : On donne le flux de masse et RT (temperature fois 556 : la constante de gaz); on impose que la vitesse 557 : tangentielle est nulle, on recupere la pression. 558 : 559 : CHPLI : condition limite impose (SPG ='DOMA' $MOD2 'CENTRE') 560 : Les composantes dependent de MOT1 561 : 562 : Resultats 563 : 564 : RCHPLI : CHPOINT qui contient les valeur de la densite, de la vitesse et 565 : de pression associe a la condition limite. 566 : (SPG ='DOMA' $MOD2 'CENTRE', composantes = LMOTP) 567 : 568 : RCHPRE : CHPOINT qui contient la contribution de la condition limite au 569 : residu (SPG en 'DOMA' $MOD1 'CENTRE', composantes = LMOTC) 570 : 571 : RJACO : MATRIK qui contient la contribution de la condition limite a la 572 : matrice jacobienne du residu. 573 : 574 : 575 : IIIb :gaz thermiquement parfait 576 : ____ 577 : 578 : Discretisation en VF "cell-centered" des equations d'Euler pour 579 : un melange de gaz parfaits. 580 : Inconnues: densite (totale) du melange, quantite de mouvement, 581 : energie totale par unite de volume, densites des constituants du gaz, 582 : scalaires passifs (multiplies par la densite du melange) 583 : 584 : RCHPO1 RFLOT1 = 'KONV' 'VF' 'PERFTEMP' MOT1 MOT2 585 : MOD1 TAB2 LMOT1 MCHAM1 MCHAM2 MCHAM3 586 : (MCHAM4) (MCHAM5) ; 587 : 588 : ENTREES 589 : 590 : 591 : MOT1 : objet de type MOT 592 : Il vaut 'RESI' si on veut calculer le residu 593 : Il vaut 'FLUX' si on veut calculer le flux 594 : 595 : MOT2 : objet de type MOT 596 : Il indique la methode de decentrement: 597 : 'VLH' = solveur de van Leer Hanel 598 : 'SS' = solveur choc-choc 599 : 600 : MOD1 : objet MODELE 601 : 602 : TAB2 : table qui contient les proprietes du gaz. Plus precisement: 603 : * le nom de l'espece qui n'est pas dans les Equations d'Euler 604 : en TAB2 . 'ESPNEULE' (MOT); 605 : * les noms des especes qui apparaissent explicitement dans les 606 : equations d'Euler en TAB2 . 'ESPEULE' (LISTMOTS); dans le 607 : cas mono-espece cet indice n'existe pas. 608 : * le degre des polynomes cv_i=cv_i(T), en TAB2 . 'NORD' 609 : (ENTIER >= 0) 610 : * les proprietes de chaque gaz 'ESPI', dans la table 611 : TAB2 . 'ESPI': 612 : - TAB2 . 'ESPI' . 'A' (LISTREEL) 613 : qui contient les (TAB2.'NORD')+1 coefficients des polynomes 614 : cv(T), (A0,A1,...); si le cv(T) sont supposes etre en 615 : J/kg/K, les Ai doivent etre en unites coherentes; 616 : - TAB2 . 'ESPI' . 'R' (J/kg/K dans le SI, FLOTTANT) qui 617 : contient la constante du gaz parfait 618 : - TAB2 . 'ESPI' . 'H0K' (J/kg, FLOTTANT) qui contient 619 : l'enthalpie de formation du gaz a 0K 620 : * TAB2 . 'SCALPASS' = si existante, noms des scalaires passifs 621 : a transporter (LISTMOTS) 622 : 623 : LMOT1 : objet de type LISTMOTS 624 : Noms de composantes du resultat (RCHPO1) 625 : Il contient dans l'ordre suivant: le nom de la densite 626 : totale, de la quantite de mouvement, de l'energie totale 627 : par unite de volume, des densites des constituants du gaz 628 : (en TAB2 . 'ESPEULE'), des scalaires passifs 629 : (en TAB2 . 'SCALPASS') fois la densite 630 : 631 : MCHAM1: MCHAML contenant la masse volumique, qui a comme SPG (support 632 : geometrique) 'DOMA' MOD1 'FACEL' (une composante, 'SCAL') 633 : (voir operateur 'PRET') 634 : 635 : MCHAM2: MCHAML contenant la vitesse et les cosinus 636 : directeurs du repere locale (n,t) dans le repere 637 : global (x,y) (dans le cas 2D 6 composantes: 638 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 639 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 640 : * 'NX' = n.x (SPG = 'FACE') 641 : * 'NY' = n.y (SPG = 'FACE') 642 : * 'TX' = t.x (SPG = 'FACE') 643 : * 'TY' = t.y (SPG = 'FACE')). 644 : (voir operateur 'PRET') 645 : 646 : MCHAM3: MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression du 647 : gaz (une seule composante, 'SCAL'). 648 : (voir operateur 'PRET') 649 : 650 : (MCHAM4):MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant les fractions 651 : massiques des especes qui apparaissent explicitement dans 652 : les equations d'Euler (voir TAB2 . 'ESPEULE') ; 653 : Dans le cas mono-espece, MCHAM4 n'est pas a donner 654 : (voir operateur 'PRET') 655 : 656 : (MCHAM5):MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant les scalaires 657 : transportes (noms des composantes en TAB2 . 'SCALPASS'') ; 658 : (voir operateur 'PRET') 659 : 660 : SORTIES 661 : 662 : 663 : RCHPO1: objet de type CHPOINT (composantes = LMOT1) 664 : Residu si MOT1 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE') 665 : Flux si MOT1 = 'FLUX' (SPG = 'DOMA' MOD1 'FACE') 666 : 667 : RFLOT1: objet de type FLOTTANT 668 : Il est le temps caracteristique associe a l'onde la plus 669 : rapide. 670 : 671 : Remarque 672 : -------- 673 : 674 : RCHPO1 est egal a: 675 : * la derive temporelle des inconnues si l'option 'RESI' est utilisee 676 : * la composante normale a la face du flux convectif (les normales 677 : etant donnees par 'DOMA' MOD1 'XXNORMAF') si l'option 'FLUX' 678 : est utilisee 679 : 680 : 681 : 682 : IIIc : gaz parfait, "Free matrix method" 683 : ____ 684 : 685 : Discretisation en VF "cell-centered" des equations d'Euler pour un gaz 686 : parfait. 687 : Inconnues: U (densite, quantite de mouvement, energie totale par unite 688 : de volume). 689 : Implicitation via une methode "sans matrice" (a stoker), ou la 690 : convergence a l'etat stationnaire est effectue en utilisant la methode 691 : de Rusanov. Dans la cellule i-eme on doit calculer U_i^{n+1} t.q. 692 : 693 : (U_i^{n+1} - U_i^{n}) * AN_i(U^{n}) = 694 : RES_i(U^{n}) + (BN_i(U^{n}) - BN_i(U_i^{n+1})) 695 : 696 : ou 697 : AN_i contient les contributions liees au reciproque du pas de temps 698 : (local) et aux diffusivites numeriques (Rusanov) aux interfaces; 699 : RES_i contient le residu calcule avec n'import quelle methode 700 : numerique (voir operateur KONV, cas IIIa). 701 : BN_i contient les contributions aux interfaces liees au flux 702 : numerique centres et la diffusivite numerique (Rusanov) multiplie par 703 : l'etat voisin. 704 : On remarque que a l'etat stationnaire on a 705 : 706 : RES_i(U)=0 707 : 708 : Pour calculer le CHPOINT RES on utilise l'operateur KONV, cas IIIa. 709 : Pour calculer les CHPOINTs AN et BN: 710 : 711 : AN = 'KONV' 'VF' 'PMONOFMM' 'AN' LMOT1 MOD1 712 : CHPO1 CHPO2 CHPO3 CHPO4 FLOT1 ('CLIM' LMOT2 CHPO5) ; 713 : 714 : BN = 'KONV' 'VF' 'PMONOFMM' 'BN' LMOT1 MOD1 715 : CHPO1 CHPO2 CHPO3 CHPO4 ('CLIM' LMOT2 CHPO5) ; 716 : 717 : ENTREES 718 : 719 : MOD1 : objet MODELE 'EULER' 720 : 721 : LMOT1 : objet de type LISTMOTS 722 : Noms de composantes du resultat (RCHPO1) 723 : Il contient dans l'ordre suivant: le nom de la densite 724 : totale, de la quantite de mouvement, de l'energie totale 725 : par unite de volume. 726 : 727 : CHPO1 : CHPOINT contenant la masse volumique (en kg/m3; une 728 : composante, 'SCAL'). 729 : 730 : CHPO2 : CHPOINT contenant les debits (en kg/s/m2; deux 731 : composantes en 2D, 'UX ','UY ', trois composantes 732 : en 3D, 'UX ','UY ', 'UZ '). 733 : 734 : CHPO3 : CHPOINT contenant l'energie totale par unite de volume 735 : (en J/m3; une composante, 'SCAL'). 736 : 737 : CHPO4 : CHPOINT contenant le "gamma" du gaz (une composante, 738 : 'SCAL'). 739 : 740 : FLOT1 : flottant, le double de la CFL 741 : 742 : LMOT2 : objet de type LISTMOTS 743 : Noms de composantes de conditions aux bords (CHPO5) 744 : Il contient dans l'ordre suivant: le nom de la densite 745 : totale, de la vitesse, de la pression. 746 : 747 : CHPO5 : CHPOINT contenant les conditions aux bords (densite, 748 : vitesse et pression sur le bord). 749 : 750 : SORTIES 751 : 752 : AN : objet de type CHPOINT (SPG = 'DOMA' MOD1 'CENTRE', 753 : une composante, 'SCAL'). 754 : 755 : BN : objet de type CHPOINT (SPG = 'DOMA' MOD1 'CENTRE', 756 : composantes = LMOT1) 757 : 758 : IIId : gaz parfait, "Free matrix method", Euler/NS bas Mach 759 : ____ 760 : 761 : Discretisation en VF "cell-centered" des equations d'Euler/NS pour un gaz 762 : parfait. 763 : Inconnues: U (densite, quantite de mouvement, energie totale par unite 764 : de volume). 765 : Implicitation via une methode "sans matrice" (a stoker), ou la 766 : convergence a l'etat stationnaire est effectue en utilisant la methode 767 : de Rusanov (voir IIIc). Dans la cellule i-eme on doit calculer 768 : 769 : 770 : (U_i^{n+1} - U_i^{n}) = DUN 771 : 772 : Syntaxe: 773 : 774 : DUN IPRO = 'KONV' 'VF' 'PMON1FMM' MOT1 LMOT1 MOD1 775 : CHPORE CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 776 : FLOT1 FLOT2 777 : NJAC 'CLIM' LMOT2 CHPO6 CHPO7 ; 778 : 779 : ENTREES 780 : 781 : MOT1 : objet de type MOT, methode d'inversion ('PJACO', point Jacobi, 782 : 'LJACOF', 'LJACOB', 'LJACOFB' Gauss-Seidel Forward, Backward et 783 : Symmetric) 784 : 785 : LMOT1 : objet de type LISTMOTS 786 : Noms de composantes du resultat (RCHPO1) 787 : Il contient dans l'ordre suivant: le nom de la densite 788 : totale, de la quantite de mouvement, de l'energie totale 789 : par unite de volume. 790 : 791 : MOD1 : objet MODELE 'EULER' 792 : 793 : CHPORE: CHPOINT contenant le residu explicite (Euler/NS) 794 : 795 : CHPO1 : CHPOINT contenant la masse volumique (en kg/m3; une 796 : composante, 'SCAL'). 797 : 798 : CHPO2 : CHPOINT contenant les debits (en kg/s/m2; deux 799 : composantes en 2D, 'UX ','UY ', trois composantes 800 : en 3D, 'UX ','UY ', 'UZ '). 801 : 802 : CHPO3 : CHPOINT contenant l'energie totale par unite de volume 803 : (en J/m3 en SI, une composante, 'SCAL'). 804 : 805 : CHPO4 : CHPOINT contenant le "gamma" du gaz (une composante, 806 : 'SCAL'). 807 : 808 : CHPO5 : CHPOINT contenat la vitesse de cut-off (une composante, 809 : 'SCAL'). 810 : 811 : FLOT1 : flottant, le pas de temps physique 812 : 813 : FLOT2 : flottant, le double de la CFL "duale" 814 : 815 : LMOT2 : objet de type LISTMOTS 816 : Noms de composantes de conditions aux bords (CHPO6) 817 : Il contient dans l'ordre suivant: le nom de la densite 818 : totale, de la vitesse, de la pression. 819 : 820 : CHPO6 : CHPOINT contenant les conditions aux bords (densite, 821 : vitesse et pression sur le bord). 822 : 823 : CHPO7 : CHPOINT contenant le rayon spectral visquex (une composante, 824 : 'SCAL'). 825 : 826 : SORTIES 827 : 828 : DUN : objet de type CHPOINT (SPG = 'DOMA' MOD1 'CENTRE', 829 : composantes = LMOT1). 830 : 831 : IPRO : entier, 0 si tout s'est bien passe 832 : 833 : 834 : IIIg. Gas parfait multi-constituant 835 : ___________________________________ 836 : 837 : IIIg.1 Le residu 838 : __________________ 839 : 840 : 841 : RESIDU DELTAT = 'KONV' 'VF' 'PERFMULT' 'RESI' MOT1 LMOT1 MOD1 842 : MCHAM1 MCHAM2 MCHAM3 MCHAM4 MCHAM5 843 : TABGAS (CHPO6) (CHPO7) (MAIL1) ; 844 : 845 : ENTRÉES 846 : 847 : 848 : MOT1 : objet de type MOT 849 : Il indique la methode: 850 : 'GODUNOV' = solveur exacte 851 : 'VANLEER' = solveur de van Leer 852 : 'VLH' = solveur de van Leer Hanel 853 : 'HUSVL' = HUS (van Leer + Osher) 854 : 'HUSVLH' = HUS (van Leer Hanel + Osher) 855 : 'AUSMPLUS' = AUSM+ 856 : 'ROE' = solveur de Roe 857 : 'SS' = solveur choc-choc 858 : 'AUSMPLM' = AUSM+ low Mach 859 : 'RUSANOV' = solveur de Rusanov 860 : 'RUSANOLM' = solveur de Rusanov pour le bas-Mach 861 : 'CENTERED' = schema centre 862 : 'ROELM' = solveur de Roe-Turkel pour le bas-Mach 863 : 'HLLC' = solveur HLLC 864 : 'HLLCLM' = solveur HLLC-Turkel pour le bas-Mach 865 : 'AUSMPUP' = solveur AUSM+up low Mach 866 : 867 : 868 : LMOT1 : objet de type LISTMOTS 869 : Noms de composantes du resultat (RCHPO1) 870 : Il contient dans l'ordre suivant: le noms de la densite, 871 : de la qdm, de l'energie totale par unite de 872 : volume, de l'especes qui sont dans TABGAS.'ESPEULE'. 873 : 874 : MOD1 : objet MODELE. 875 : 876 : MCHAM1 : MCHAML contenant la masse volumique, qui a comme 877 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 878 : (une composante, 'SCAL') 879 : (voir operateur PRET) 880 : 881 : MCHAM2 : MCHAML contenant la vitesse et les cosinus 882 : directeurs du repere locale (n,t) dans le repere 883 : global (x,y) (dans le cas 2D 6 composantes: 884 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 885 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 886 : * 'NX' = n.x (SPG = 'FACE') 887 : * 'NY' = n.y (SPG = 'FACE') 888 : * 'TX' = t.x (SPG = 'FACE') 889 : * 'TY' = t.y (SPG = 'FACE')). 890 : (voir operateur PRET) 891 : 892 : MCHAM3 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression du 893 : gaz (une seule composante, 'SCAL'). 894 : (voir operateur PRET) 895 : 896 : MCHAM4 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant le "gamma" du 897 : gaz (une seule composante, 'SCAL'). 898 : (voir operateur PRET) 899 : 900 : MCHAM5 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant les 901 : fractions massiques (nombre des composants egal 902 : celui dans TABGAS.'ESPEULE') 903 : 904 : TABGAS : la table contenant les properties de gas 905 : (voir PRET ou PRIM) 906 : 907 : CHPO6 : CHPOINT contenant la vitesse de cut-off 908 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 909 : 'SCAL'). 910 : A donner dans le cas bas Mach 911 : 912 : CHPO7 : CHPOINT contenant la deuxieme vitesse de cut-off 913 : (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 914 : 'SCAL'). 915 : A donner dans le cas bas Mach 916 : 917 : (MAIL1) : MAILLAGE de POI1 ou, si appartenant a ('DOMA' MOD1 'FACE'), 918 : on ne calcule pas la contribution au flux ou au residu. 919 : 920 : SORTIES 921 : 922 : RESIDU : objet de type CHPOINT (composantes = LMOT1) 923 : Residu si MOT2 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE') 924 : 925 : DELTAT : objet de type FLOTTANT 926 : Il est le temps caracteristique associe a l'onde la plus 927 : rapide (meme dans le cas bas Mach, on considere le systeme 928 : non preconditionne) 929 : 930 : Remarque 931 : -------- 932 : 933 : RESIDU est egal a: 934 : * la derive temporelle des inconnues si l'option 'RESI' est utilisee 935 : 936 : 937 : 938 : IIIg.2 La matrice jacobienne du residu par rapport aux variables 939 : _____ conservatives 940 : 941 : 942 : RMAT1 = 'KONV' 'VF' 'PERFMULT' 'JACOCONS' MOT1 MOD1 943 : TABGAS LMOT1 (MAIL1) 944 : CHPO1 CHPO2 CHPO3 CHPO4 (CHPO5 CHPO6) ; 945 : 946 : 947 : ENTRÉES 948 : 949 : MOT1 : objet de type MOT 950 : 'VLH' : jacobien du residu pour la methode VLH 951 : 'AUSMPLUS' : jacobien du residu pour la methode AUSM+ 952 : 'AUSMPLM' : jacobien du residu pour la methode AUSM+ low Mach 953 : 954 : TABGAS : la table contenant les properties de gas 955 : (voir PRET ou PRIM) 956 : 957 : LMOT1 : objet de type LISTMOTS 958 : Noms des variables conservatives 959 : Il contient dans l'ordre suivant: le noms de la densite, 960 : de la qdm, de l'energie totale par unite de volume, 961 : de l'especes qui sont dans TABGAS.'ESPEULE'. 962 : 963 : MOD1 : objet MODELE. 964 : 965 : CHPO1 : CHPOINT contenant la masse volumique 966 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 967 : 'SCAL'). 968 : 969 : CHPO2 : CHPOINT contenant la vitesse 970 : (SPG = 'DOMA' DOM1 'CENTRE', deux/trois composantes 971 : 'UX', 'UY', 'UZ') 972 : 973 : CHPO3 : CHPOINT contenant la pression du gaz 974 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 975 : 'SCAL'). 976 : 977 : CHPO4 : CHPOINT contenant les especes 978 : (SPG = 'DOMA' DOM1 'CENTRE', nombre des composantes 979 : egal a celui dans TABGAS.'ESPEULE'). 980 : 981 : CHPO5 : CHPOINT contenant la vitesse de cut-off 982 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 983 : 'SCAL'). 984 : A donner dans le cas bas Mach 985 : 986 : CHPO6 : CHPOINT contenant la deuxieme vitesse de cut-off 987 : (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 988 : 'SCAL'). 989 : A donner dans le cas bas Mach 990 : 991 : (MAIL1) : MAILLAGE de POI1 ou, si appartenant a ('DOMA' DOM1 'FACE'), 992 : on ne calcule pas la contribution au jacobien. 993 : 994 : SORTIES 995 : 996 : RMAT1 : objet de type MATRIK 997 : (SPG = 'DOMA' DOM1 'CENTRE') 998 : (inconnues primales = inconnues duales = LMOT1) 999 : Il contient le jacobien du residu par rapport aux variables 1000 : conservatives. 1001 : 1002 : IIIg.3 La contribution de quelque condition limite au residu et 1003 : ------ a la matrice jacobienne 1004 : 1005 : 1006 : RCHPLI RCHPRE = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' 1007 : $MOD1 $MOD2 TABGAS LMOTC LMOTP 1008 : CHPRN CHPVN CHPPN CHPYN 1009 : CHPLI MOT1 ; 1010 : 1011 : ou 1012 : 1013 : RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' 1014 : $MOD1 $MOD2 TABGAS LMOTC LMOTP 1015 : CHPRN CHPVN CHPPN CHPYN 1016 : CHPLI MOT1 ; 1017 : 1018 : 1019 : $MOD1 : l'objet modele du domaine total 1020 : 1021 : $MOD2 : l'objet modele du domaine du bord 1022 : 1023 : TABGAS : la table contenant les properties de gas 1024 : (voir PRET ou PRIM) 1025 : 1026 : LMOTC : LISTMOTS, noms des variables conservatives 1027 : 1028 : LMOTP : LISTMOTS, noms des variables primitives 1029 : 1030 : CHPRN : densite (SPG = 'DOMA' $MOD1 'CENTRE', une seule 1031 : composante, 'SCAL') 1032 : 1033 : CHPVN : vitesse (SPG = 'DOMA' $MOD1 'CENTRE', 1034 : composantes: 'UX', 'UY', ('UZ')) 1035 : 1036 : CHPPN : pression (SPG = 'DOMA' $MOD1 'CENTRE', une seule 1037 : composante, 'SCAL') 1038 : 1039 : CHPYN : les especes (SPG = 'DOMA' $MOD1 'CENTRE', le nombre 1040 : des composantes egal a celui dans TABGAS.'ESPEULE') 1041 : 1042 : MOT1 : MOT, type de condition limite: 1043 : 'INRI' : entree subsonique. On utilise les invariants de 1044 : Riemann. On donne densite 'RN', vitesse 'UX' 'UY' 1045 : ('UZ'), pression 'PN', les fractions massiques 1046 : 'INSU' : entree subsonique. On donne l'enthalpie totale 'HT' 1047 : par unite de masse, l'entropie 'S' (pression divisee 1048 : par la densite a la puissance gamma),les fractions massiques 1049 : On impose que 1050 : la vitesse tangentielle est nulle. On recupere la 1051 : pression a l'interieur. La contribution au residu et 1052 : a la matrice jacobienne sont calculee via la methode 1053 : 'AUSMPLUS'. 1054 : 'OUTP' : sortie subsonique. On donne la pression 'PN', 1055 : on recupere la vitesse, la densite 1056 : et les fractions massiques a l'interieur. 1057 : La contribution au residu et a la matrice jacobienne 1058 : sont calculee via la methode 'AUSMPLUS'. 1059 : 'INJE' : condition limite d'injection compressible. 1060 : On donne le flux de masse 'MOME', 'RT' (temperature fois 1061 : la constante de gaz) et les fractions massiques; 1062 : on impose que la vitesse 1063 : tangentielle est nulle, on recupere la 1064 : pression. 1065 : 'INSS': entree supersonique. On impose les 1066 : conditions limites dans le flux. 1067 : On donne la densite, la vitesse, 1068 : la pression et les fractions massiques 1069 : 'OUTSS': sortie supersonique. On donn le 1070 : Champ par points vide. 1071 : Toute l'information est recoupere de 1072 : l'interieur. 1073 : 'RESE' : condition limite de type reservoir. 1074 : On donne la pression 'PN' la densite 'RN' et les 1075 : fractions massiques ; on impose que la vitesse 1076 : tangentielle est nulle, on recupere la 1077 : pression interieure et, selon sa valeur, 1078 : on impose les conditions de col sonique ou 1079 : subsonique. 1080 : 1081 : CHPLI : condition limite impose (SPG ='DOMA' $MOD2 'CENTRE') 1082 : Les composantes dependent de MOT1 1083 : 1084 : Resultats 1085 : 1086 : RCHPLI : CHPOINT qui contient les valeur de la densite, de la vitesse, 1087 : de pression et des fractions massiques associe a la condition limite. 1088 : (SPG ='DOMA' $MOD2 'CENTRE', composantes = LMOTP) 1089 : 1090 : RCHPRE : CHPOINT qui contient la contribution de la condition limite au 1091 : residu (SPG en 'DOMA' $MOD1 'CENTRE', composantes = LMOTC) 1092 : 1093 : RJACO : MATRIK qui contient la contribution de la condition limite a la 1094 : matrice jacobienne du residu. 1095 : 1096 : 1097 : IV Transport des scalaires 1098 : __________________________ 1099 : 1100 : 1101 : Discretisation en VF "cell-centered" du transport des scalaires 1102 : --> --> 1103 : d/dt S + div . (u S) = 0 1104 : 1105 : RCHPO1 RFLOT1 = 'KONV' 'VF' 'CLAUDEIS' 'FACE' MOT1 MOT2 1106 : MOD1 CHPO1 MCHAM1 ; 1107 : 1108 : or 1109 : 1110 : RMAT1 = 'KONV' 'VF' 'CLAUDEIS' 'FACE' 'JACO' MOT2 1111 : MOD1 CHPO1 MCHAM1 ; 1112 : 1113 : ENTRÉES 1114 : 1115 : MOT1 : objet de type MOT 1116 : Il vaut 'RESI' si on veut calculer le residu, 1117 : (i.e. l' increment des scalaires passives par unite 1118 : de temps) 1119 : Il vaut 'FLUX' si on veut calculer le flux aux 1120 : interfaces 1121 : 1122 : MOT2 : objet de type MOT 1123 : Il indique la methode de calcul du flux 1124 : Pour l'instant 1125 : 'UPWIND' 1126 : 'CENTERED' 1127 : 1128 : MOD1 : objet MODELE. 1129 : 1130 : CHPO1 : CHPOINT contenant la vitesse aux FACEs 1131 : (SPG = 'DOMA' MOD1 'FACE') 1132 : 2/3 composantes, 'UX', 'UY', 'UZ' 1133 : 1134 : MCHAM1 : MCHAML contenant les scalaires a transporter 1135 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 1136 : 1137 : SORTIES 1138 : 1139 : RCHPO1 : objet de type CHPOINT (memes composantes que MCHAM1) 1140 : Residu si MOT2 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE') 1141 : Flux si MOT2 = 'FLUX' (SPG = 'DOMA' MOD1 'FACE') 1142 : 1143 : RFLOT1 : objet de type FLOTTANT 1144 : Il est le temps caracteristique associe a l'onde la plus 1145 : rapide. 1146 : 1147 : RMAT1 : objet de type MATRIK 1148 : (SPG = 'DOMA' MOD1 'CENTRE') 1149 : (inconnues primales = inconnues duales; memes composantes 1150 : que MCHAM1) 1151 : 1152 : 1153 : V Flux bifluide 1154 : ___________________________________________________ 1155 : 1156 : 1157 : Discretisation en VF "cell-centered" des equations bifluides pour un 1158 : melange d'eau et air. 1159 : 1160 : Inconnues: 1161 : 1162 : mass volumique, quantite de mouvement (qdm) et energie totale par 1163 : unite de volume (variables conservatives) de chaque fluide 1164 : 1165 : ou 1166 : 1167 : taux de vide, vitesses du gaz et du liquide, pression, 1168 : temperature du gaz et du liquide (variables primitives) 1169 : 1170 : On peut calculer: 1171 : 1172 : RCHPO2 RCHPO1 RFLOT1 = 'KONV' 'VF' 'TWOFLUID' MOT1 MOT2 MOD1 LMOT1 1173 : MCHAM1 MCHAM2 MCHAM3 MCHAM4 MCHAM5 MCHAM6 MCHAM7 MCHAM8; 1174 : 1175 : ENTRÉES 1176 : 1177 : MOT1 : objet de type MOT 1178 : Il vaut 'RESI' si on veut calculer le residu 1179 : Il vaut 'FLUX' si on veut calculer le flux 1180 : 1181 : MOT2 : object de type MOT 1182 : Il indique la methode: 1183 : 'AUSMP1' = AUSM+ 1184 : 'AUSMP2' = AUSM+ preconditione 1185 : 'AUSMDV1' = AUSMDV 1186 : 'AUSMDV2' = AUSMDV preconditione 1187 : 1188 : LMOT1 : object de type LISTMOTS 1189 : Noms de composantes du resultat (RCHPO1) 1190 : Il contient dans l'ordre suivant: le noms de la densite, 1191 : de la qdm, 1192 : de la masse volumique du gaz, 1193 : de la masse volumique du liquide, 1194 : de la quantite de mouvement du gaz, 1195 : de la quantite de mouvement du liquide, 1196 : de l'energie totale du gaz par unite de volume 1197 : de l'energie totale du liquide par unite de volume 1198 : 1199 : MOD1 : object MODELE 1200 : 1201 : MCHAM1 : MCHAML contenant la taux de vide, qui a comme 1202 : SPG (support geometrique) 'DOMA' MOD1 'FACEL' 1203 : (une composante, 'SCAL') 1204 : (voir operateur PRET) 1205 : 1206 : MCHAM2 : MCHAML contenant la vitesse du gaz et les cosinus 1207 : directeurs du repere locale (n,t) dans le repere 1208 : global (x,y) (dans le cas 2D 6 composantes: 1209 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 1210 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 1211 : * 'NX' = n.x (SPG = 'FACE') 1212 : * 'NY' = n.y (SPG = 'FACE') 1213 : * 'TX' = t.x (SPG = 'FACE') 1214 : * 'TY' = t.y (SPG = 'FACE')). 1215 : (voir operateur PRET) 1216 : 1217 : MCHAM3 : MCHAML contenant la vitesse du liquide et les cosinus 1218 : directeurs du repere locale (n,t) dans le repere 1219 : global (x,y) (dans le cas 2D 6 composantes: 1220 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 1221 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 1222 : * 'NX' = n.x (SPG = 'FACE') 1223 : * 'NY' = n.y (SPG = 'FACE') 1224 : * 'TX' = t.x (SPG = 'FACE') 1225 : * 'TY' = t.y (SPG = 'FACE')). 1226 : (voir operateur PRET) 1227 : 1228 : MCHAM4 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression 1229 : (une seule composante, 'SCAL'). 1230 : (voir operateur PRET) 1231 : 1232 : MCHAM5 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la temperature du 1233 : gaz (une seule composante, 'SCAL'). 1234 : (voir operateur PRET) 1235 : 1236 : MCHAM6 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la temperature du 1237 : liquide (une seule composante, 'SCAL'). 1238 : (voir operateur PRET) 1239 : 1240 : MCHAM7 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du 1241 : gaz (une seule composante, 'SCAL'). 1242 : (voir operateur PRET) 1243 : 1244 : MCHAM8 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du 1245 : liquide (une seule composante, 'SCAL'). 1246 : (voir operateur PRET) 1247 : 1248 : 1249 : SORTIES 1250 : ------- 1251 : 1252 : RCHPO1 : objet de type CHPOINT (composantes = LMOT1) 1253 : Residu si MOT2 = 'RESI' ou flux si MOT2 = 'FLUX' 1254 : (SPG = 'DOMA' MOD1 'CENTRE') 1255 : 1256 : RFLOT1 : objet de type FLOTTANT 1257 : Il est le temps caracteristique associe a l'onde la plus 1258 : rapide 1259 : 1260 : 1261 : Vl Ecoulement Reactif avec resolution "Discrete Equation Method". 1262 : ___________________________________________________________________ 1263 : 1264 : 1265 : 1266 : Discretisation en VF "cell-centered" des equations d'Euler pour 1267 : un gaz parfait multi-constituent. Methode "DEM". 1268 : 1269 : 1270 : RCHPO1 RFLOT1 = 'KONV' 'VF' 'DEM' MOT1 MOT2 MOT3 MOD1 TABG LMOT1 1271 : CHPA1 CHPA2 1272 : MCHAA1 MCHAA2 1273 : MCHAR1 MCHAR2 1274 : MCHAV1 MCHAV2 1275 : MCHAP1 MCHAP2 1276 : GRALP1 K0 EPS MAILLIM 1277 : (CHPV1 CHPV2) ; 1278 : 1279 : ARGUMENTS : 1280 : ---------- 1281 : 1282 : MOT1 : objet de type MOT 1283 : Il vaut 'RESI' si on veut calculer le 1284 : residu 1285 : 1286 : MOT2 : objet de type MOT 1287 : Il indique la methode de decentrement (partie non reactive) 1288 : 'SS' = solveur choc-choc 1289 : 'VLH' = solver de Van Leer Hanel 1290 : 'AUSMPUP'= solver AUSM+up (bas Mach) 1291 : 1292 : MOT3 : objet de type MOT 1293 : 'CONS' = la vitesse fondamentale est constante (objet REEL) 1294 : 'VARI' = la vitesse fondamentale est variable (objet CHAMPOINT) 1295 : 1296 : MOD1 : objet modele, type EULER 1297 : 1298 : TABG : objet de type TABLE; 1299 : la table contenant les properties de gas 1300 : (voir PRET ou PRIM) 1301 : 1302 : LMOT1 : objet de type LISTMOTS 1303 : contenant les noms des composants du vecteur RCHPO1. 1304 : Ils sont donnes dans l'ordre suivant : noms de alpha, 1305 : de la densite, de la vitesse, de 1306 : l'energie specifique totale. 1307 : 1308 : CHPA1 : CHPOINT contenant la fraction volumique \alpha de 1309 : la phase 1, qui a comme SPG (support geometric) 1310 : l'index 'CENTRE' de la table MOD1 (une 1311 : composante, 'SCAL') 1312 : 1313 : CHPA2 : CHPOINT contenant la fraction volumique \alpha de 1314 : la phase 2, qui a comme SPG (support geometric) 1315 : l'index 'CENTRE' de la table MOD1 (une 1316 : composante, 'SCAL') 1317 : 1318 : MCHAA1 : MCHAML contenant la fraction volumique \alpha de 1319 : la phase 1, qui a comme SPG (support geometric) 1320 : l'index 'FACEL' de la table MOD1 (une 1321 : composante, 'SCAL') 1322 : 1323 : MCHAA2 : MCHAML contenant la fraction volumique \alpha de 1324 : la phase 2, le meme SPG que pour MCHAA1, (une 1325 : composante, 'SCAL') 1326 : 1327 : MCHAR1 : MCHAML contenant la densite de la 1328 : phase 1, le meme SPG que pour MCHAA1, (une 1329 : composante, 'SCAL') 1330 : 1331 : MCHAR2 : MCHAML contenant la densite de la 1332 : phase 2, le meme SPG que pour MCHAA1, (une 1333 : composante, 'SCAL') 1334 : 1335 : 1336 : MCHAV1 : MCHAML contenant la vitesse du liquide de la 1337 : phase 1 et les cosinus 1338 : directeurs du repere locale (n,t) dans le repere 1339 : global (x,y) (dans le cas 2D 6 composantes: 1340 : * 'UN' = vitesse normale (SPG = 'DOMA' MOD1 'FACEL') 1341 : * 'UT' = vitesse tangentielle (SPG = 'DOMA' MOD1 FACEL') 1342 : * 'NX' = n.x (SPG = 'FACE') 1343 : * 'NY' = n.y (SPG = 'FACE') 1344 : * 'TX' = t.x (SPG = 'FACE') 1345 : * 'TY' = t.y (SPG = 'FACE')). 1346 : (voir operateur PRET) 1347 : 1348 : 1349 : MCHAV2 : MCHAML contenant la vitesse du liquide de la 1350 : phase 2 et les cosinus 1351 : directeurs du repere locale (n,t) dans le repere 1352 : global (x,y) (la meme structure que pour MCHAV1) 1353 : 1354 : 1355 : 1356 : MCHAP1 : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant la 1357 : la pression de la phase 1 1358 : (une composante, 'SCAL'). 1359 : 1360 : MCHAP2 : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant la 1361 : la pression de la phase 2 1362 : (une composante, 'SCAL'). 1363 : 1364 : K0 : FLOTTANT/CHAMPOINT (voir MOT3), vitesse fondamentale de la 1365 : flamme 1366 : 1367 : GRALP1 : CHPOINT, grad(alp1)/|grad(alp1)| 1368 : 1369 : EPSILON : FLOTTANT t.q. a < EPSILON => a = 0 1370 : 1371 : MAILLIM : MAILLAGE -- la partie du maillage ou le flux 1372 : n'est pas determine ; il est 1373 : calcule via les procedures pour 1374 : les Conditions Limites. 1375 : 1376 : CHPV1 : CHPOINT, cut-off speed dans la methode AUSM+up, 1377 : phase 1 1378 : 1379 : CHPV2 : CHPOINT, cut-off speed dans la methode AUSM+up, 1380 : phase 2 1381 : 1382 : RESULTS : 1383 : ----------- 1384 : 1385 : RCHPO1 : objet de type CHPOINT (le noms des composantes 1386 : sont dans LMOT1) 1387 : Residu MOT2 = 'RESI' (SPG =('DOMA' MOD1 'CENTRE')) 1388 : 1389 : RFLOT1 : objet de type FLOTTANT 1390 : Il est le temps caracteristique associe a l'onde la plus 1391 : rapide 1392 : 1393 : VII Ghost fluid method for the poor. 1394 : ___________________________________________________ 1395 : 1396 : 1397 : Discretisation en VF "cell-centered" des equations bifluides pour le 1398 : transport d'une interface. 1399 : 1400 : Inconnues: phi, densite, quantite de mouvement, energie totale par 1401 : unite de volume, densites et fractions volumiques des constituants. 1402 : 1403 : 1404 : RCHPO1 RFLOT1 = 'KONV' 'VF' 'GFMP' MOT1 MOT2 MOD1 1405 : TABG 1406 : LMOT1 MCHAPH MCHAR MCHAV MCHAP (MCHAY 1407 : MCHAA) MCHPPH LOG1 MAILLIM ; 1408 : ENTREES 1409 : 1410 : 1411 : MOT1 : objet de type MOT 1412 : Il vaut 'RESI' si on veut calculer le residu 1413 : 1414 : MOT2 : objet de type MOT 1415 : Il indique la methode de decentrement: 1416 : 'GODUNOV' 1417 : 1418 : MOD1 : objet modele de type Euler 1419 : 1420 : TABG : objet de type TABLE 1421 : il contient les proprietes du gaz 1422 : (voir operateur PRIM) 1423 : 1424 : LMOT1 : objet de type LISTMOTS 1425 : Noms de composantes du resultat (RCHPO1) 1426 : Il contient dans l'ordre suivant les noms de 1427 : phi, de la densites, de la vitesse, de l'energie totale 1428 : par unite de volume, des densites des especes, des 1429 : fractions volumiques. 1430 : 1431 : MCHAPH : MCHAML contenant la fonction phi, 1432 : qui a comme SPG (support geometrique) l'indice 'FACEL' de la 1433 : table associee a MOD1 (une composante, 'SCAL') 1434 : 1435 : MCHAR : MCHAML contenant la masse volumique, 1436 : meme SPG que MCHAA1, une composante, 'SCAL' 1437 : 1438 : MCHAV : MCHAML contenant la vitesse et les cosinus 1439 : directeurs du repere locale (n,t) dans le repere 1440 : global (x,y), meme SPG que MCHAA1, 1441 : (dans le cas 2D 6 composantes: 1442 : * 'UN' = vitesse normale (SPG =('DOMA' MOD1 'FACEL')) 1443 : * 'UT' = vitesse tangentielle (SPG =('DOMA' MOD1 'FACEL')) 1444 : * 'NX' = n.x (SPG = 'FACE') 1445 : * 'NY' = n.y (SPG = 'FACE') 1446 : * 'TX' = t.x (SPG = 'FACE') 1447 : * 'TY' = t.y (SPG = 'FACE')). 1448 : 1449 : MCHAP : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant la pression 1450 : (une seule composante, 'SCAL'). 1451 : 1452 : MCHAY : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant les fractions 1453 : massiques (composantes dans TABG . 'ESPEULE'). 1454 : 1455 : MCHAA : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant les fractions 1456 : volumiques (composantes dans TABG . 'ESPEULE'). 1457 : 1458 : MCHPPH : CHPO contenant phi, 1459 : qui a comme SPG (support geometrique) l'indice 'CENTRE' de la 1460 : table associee a MOD1 (une composante, 'SCAL') 1461 : 1462 : LOG1 : LOGIQUE, si VRAI phi et les fractions volumiques sont 1463 : traites avec une approche conservative 1464 : 1465 : MAILLIM : MAILLAGE -- points faces ou le flux n'est pas calcule. 1466 : 1467 : SORTIES 1468 : 1469 : RCHPO1 : objet de type CHPOINT (composantes = LMOT1) 1470 : Residu si MOT2 = 'RESI' (SPG =('DOMA' MOD1 'CENTRE')) 1471 : 1472 : RFLOT1 : objet de type FLOTTANT 1473 : Il est le temps caracteristique associe a l'onde la plus 1474 : rapide. 1475 : 1476 : 1477 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales