1 : * ENCEINTE PROCEDUR CHAT 12/08/07 21:15:02 7481 2 : 'DEBPROC' ENCEINTE nbit*'ENTIER ' rxt*'TABLE ' KENC/'MOT'; 3 : *---------------------------------------------------------- 4 : *---- Suivi des modifications 5 : *--- 20/09/99 : ES prise en compte des incondensables 6 : * dans la procédure CONDENS 7 : *--- 18/05/00 : FP prises en comptes de la condensation en masse par 8 : * l'ajout d'une nouvelle équation sur la masse volumique 9 : * de liquide en suspension. Les équations sur rhovap et 10 : * sur T sont également complétées par l'introduction 11 : * d'une nouveau terme source. La détermination des 12 : * termes sources est réalisée dans une procedure appelée 13 : * CONDENSM. 14 : *--- 04/11/05 : Ajout d'un choix de formulation pour le flux condensé 15 : * (Chilton-Colburn) dans la procédure CONDENS 16 : *---------------------------------------------------------- 17 : 18 : *============================================================== 19 : Si(Exist KENC) ; 20 : MKENC=KENC ; 21 : Sinon ; 22 : MKENC='????' ; 23 : Finsi ; 24 : *----------------------------------------------------------------------- 25 : * Liste des entrées nécessaires ; 26 : *----------------------------------------------------------------------- 27 : * Géométrie 28 : 29 : listand = MOTS ; 30 : listchp = MOTS ; 31 : listentr= table ; 32 : listentr.'VERSION' ='' ; 33 : listentr.'vtf' ='' ; 34 : listentr.'TF0' ='' ; 35 : listentr.'PT0' ='' ; 36 : 37 : listentr.'THERMP' ='' ; 38 : listentr.'THERCO' ='' ; 39 : listentr.'vtp' ='' ; 40 : listentr.'ROCP' ='' ; 41 : listentr.'LAMBDA' ='' ; 42 : listentr.'Tp0' ='' ; 43 : listentr.'ECHAN' ='' ; 44 : 45 : listentr.'ECHEXT' ='' ; 46 : listentr.'parext' ='' ; 47 : listentr.'HEXT' ='' ; 48 : listentr.'TPEXT' ='' ; 49 : 50 : listentr.'TIMP1' ='' ; 51 : listentr.'TIMP2' ='' ; 52 : listentr.'TIMP3' ='' ; 53 : 54 : listentr.'ECHANP' ='' ; 55 : 56 : listentr.'breche' ='' ; 57 : listentr.'breche2' ='' ; 58 : listentr.'breche3' ='' ; 59 : listentr.'sortie' ='' ; 60 : listentr.'axe' ='' ; 61 : listentr.'epsi' ='' ; 62 : listentr.'pi' ='' ; 63 : 64 : * Table utilisée pour la définition de recombineurs 65 : listentr.'RECOMB' ='' ; 66 : 67 : * Table de travail rajoutées par la procédure 68 : listentr.'REINIT' ='' ; 69 : listentr.'TIC' ='' ; 70 : listentr.'GEO' ='' ; 71 : listentr.'TBT' ='' ; 72 : 73 : * paramètres numériques / algorithmes 74 : listentr.'DISCR' ='' ; 75 : listentr.'KPRE' ='' ; 76 : listentr.'STAB' ='' ; 77 : listentr.'ALGO' ='' ; 78 : 79 : * paramètres de fonctionnement 80 : listentr.'DT0' ='' ; 81 : listentr.'PERSO' ='' ; 82 : listentr.'GRAPH' ='' ; 83 : listentr.'RENU' ='' ; 84 : listentr.'TYPINV' ='' ; 85 : listentr.'IMPR' ='' ; 86 : listentr.'FRPREC' ='' ; 87 : listentr.'PPDT' ='' ; 88 : 89 : 90 : ** Modèles physiques 91 : * Constituant (air par défaut ) vapeur/incondensables 92 : listentr.'VAPEUR' ='' ; 93 : listentr.'HE' ='' ; 94 : listentr.'H2' ='' ; 95 : listentr.'O2' ='' ; 96 : listentr.'N2' ='' ; 97 : listentr.'CO' ='' ; 98 : listentr.'CO2' ='' ; 99 : listentr.'AIR' ='' ; 100 : 101 : * Modèles de turbulence 102 : listentr.'MODTURB' ='' ; 103 : listentr.'FPAROI' ='' ; 104 : 105 : * Formulation pour le modèle de condensation 106 : listentr.'MODCOND' ='' ; 107 : 108 : idim = vale 'DIME' ; 109 : DIM3D=FAUX ; 110 : Si(EGA idim 3) ; 111 : DIM3D=VRAI ; 112 : Finsi ; 113 : 114 : *----------------------------------------------------------------------- 115 : 116 : REINIT= (NON ('EXIST' rxt 'TBT')) ; 117 : 'SI' ('EXIST' rxt 'REINIT') ; 118 : vertytab rxt 'REINIT' 'LOGIQUE' ; 119 : REINIT=rxt.'REINIT' ; 120 : 'FINSI' ; 121 : 122 : 'SI'(('EXIST' rxt 'TBT') et (NON REINIT)) ; 123 : vertytab rxt 'TBT' 'TABLE ' ; 124 : vertytab rxt 'GEO' 'TABLE ' ; 125 : TPREPA=FAUX ; 126 : TBT=rxt.'TBT' ; 127 : GEO=rxt.'GEO' ; 128 : 'SINON' ; 129 : TPREPA=VRAI ; 130 : TBT='TABLE' ; 131 : GEO='TABLE' ; 132 : rxt.'GEO'=GEO ; 133 : 'FINSI' ; 134 : 135 : 'SI' ('EXIST' rxt 'TIC') ; 136 : vertytab rxt 'TIC' 'TABLE ' ; 137 : vertytab (rxt.'TIC') 'LTPS' 'LISTREEL' ; 138 : tic=rxt.'TIC' ; 139 : 'SINON' ; 140 : tic= 'TABLE' 'INCO' ; 141 : rxt.'TIC'=tic ; 142 : tic.'LTPS' = 'PROG' 0.0 ; 143 : 'FINSI' ; 144 : 145 : ndl = 'DIME' (tic.'LTPS') ; 146 : Tps='EXTR' (tic.'LTPS') ndl ; 147 : MESS ' Le calcul demarre a partir du temps ' Tps ; 148 : 149 : * Initialisation 150 : * Liste des Gaz qui peuvent etre calcule 151 : LINCOND=MOTS 'N2' 'H2' 'HE' 'O2' 'CO' 'CO2' 'VAPE' ; 152 : 153 : Si (non(exist tic 'Lg0D')) ; 154 : tic.'Lg0D'= Mots ; 155 : tic.'Li0D'= Prog ; 156 : Finsi ; 157 : 158 : *------------------------------------------------------------- 159 : * On vérifie que l'utilisateur a bien crée un volume FLUIDE 160 : * Donnée obligatoire 161 : *------------------------------------------------------------- 162 : 163 : vertytab rxt 'vtf' 'MAILLAGE' ; 164 : vertytab rxt 'TF0' 'FLOTTANT' ; 165 : vertytab rxt 'PT0' 'FLOTTANT' ; 166 : 167 : *------------------------------------------------------------- 168 : * DT0 Le pas de temps : Donnée obligatoire 169 : *------------------------------------------------------------- 170 : 171 : vertytab rxt 'DT0' 'FLOTTANT' ; 172 : tic.'DT'=rxt.'DT0' ; 173 : 174 : *---------------------------------------------------------- 175 : * VERSION= 'MOT' indiquant le niveau de la version 176 : * Niveau disponible V0 (standard) et V1 177 : *---------------------------------------------------------- 178 : 179 : TBT.'VERSION' = 0 ; 180 : 'SI' (exist rxt 'VERSION') ; 181 : lxt = mots 'V0' 'V1' ; 182 : tbvs= table; tbvs.'V0'=0; tbvs.'V1'=1; ; 183 : vertytab lxt (rxt.'VERSION') 'MOT' ; 184 : TBT.'VERSION'= tbvs.(rxt.'VERSION') ; 185 : 'SI' (TBT.'VERSION' > 0) ; 186 : TBT.'CONDMAS' = VRAI ; 187 : 'FINSI' ; 188 : 'FINSI' ; 189 : 190 : *---------------------------------------------------------- 191 : * GRAPH = Logique permettant les dessins 192 : *---------------------------------------------------------- 193 : 194 : 195 : TBT.'GRAPH' = FAUX ; 196 : 'SI' (exist rxt 'GRAPH') ; 197 : vertytab rxt 'GRAPH' 'LOGIQUE ' ; 198 : TBT.'GRAPH'= rxt.'GRAPH' ; 199 : 'FINSI' ; 200 : 201 : *--------------------------------------------------------- 202 : * STAB = coefficient de Stabilisation pour les MACRO 203 : * pour l'algorithme de Projection 204 : *--------------------------------------------------------- 205 : 206 : TBT.'STAB' = 1. ; 207 : 'SI' ('EXIST' rxt 'STAB') ; 208 : vertytab rxt 'STAB' 'FLOTTANT' ; 209 : TBT.'STAB' = rxt.'STAB' ; 210 : listand=listand et (MOTS 'STAB') ; 211 : 'FINSI' ; 212 : 213 : *--------------------------------------------------------- 214 : * DETMAT Logique indiquant qu'on détruit toutes les 215 : * matrices à la fin du calcul. 216 : *--------------------------------------------------------- 217 : 218 : TBT.'DETMAT' = VRAI ; 219 : 'SI' ('EXIST' rxt 'DETMAT') ; 220 : vertytab rxt 'DETMAT' 'LOGIQUE' ; 221 : listentr.'DETMAT' ='' ; 222 : TBT.'DETMAT' = rxt.'DETMAT' ; 223 : 'FINSI' ; 224 : 225 : *--------------------------------------------------------- 226 : * CORTEMP Logique indiquant qu'on supprime ou non le controle 227 : * du bilan d'énergie. 228 : *--------------------------------------------------------- 229 : 230 : TBT.'CORTEMP' = VRAI ; 231 : 'SI' ('EXIST' rxt 'CORTEMP') ; 232 : vertytab rxt 'CORTEMP' 'LOGIQUE' ; 233 : listentr.'CORTEMP' ='' ; 234 : TBT.'CORTEMP' = rxt.'CORTEMP' ; 235 : listand=listand et (MOTS 'CORTEMP') ; 236 : 'FINSI' ; 237 : 238 : *--------------------------------------------------------- 239 : * RENU = Type de renumérotation pour les méthodes d'inversion 240 : * RIEN SLOA 241 : *--------------------------------------------------------- 242 : 243 : TBT.'RENU' = 'SLOA' ; 244 : 'SI' ('EXIST' rxt 'RENU') ; 245 : vertytab rxt 'RENU' 'MOT' ; 246 : lxt=mots 'RIEN' 'SLOA' 'GIBR' 'GIBA' ; 247 : vertytab lxt (rxt.'RENU') 'MOT' ; 248 : TBT.'RENU' = rxt.'RENU' ; 249 : listand=listand et (MOTS 'RENU') ; 250 : 'FINSI' ; 251 : 252 : *--------------------------------------------------------- 253 : * TYPINV Type d'inversion 254 : * RIEN BICG sinon voir notice KRES 255 : *--------------------------------------------------------- 256 : 257 : TBT.'TYPINV'= 3 ; 258 : 'SI' ('EXIST' rxt 'TYPINV') ; 259 : vertytab rxt 'TYPINV' 'ENTIER' ; 260 : TBT.'TYPINV' = rxt.'TYPINV' ; 261 : listand=listand et (MOTS 'TYPINV') ; 262 : 'FINSI' ; 263 : 264 : *--------------------------------------------------------- 265 : * PPDT formulation designé petits/grand pas de temps 266 : * conservative ou non conservative 267 : *--------------------------------------------------------- 268 : 269 : TBT.'PPDT'= FAUX ; 270 : 'SI' ('EXIST' rxt 'PPDT') ; 271 : vertytab rxt 'PPDT' 'LOGIQUE' ; 272 : TBT.'PPDT' = rxt.'PPDT' ; 273 : listand=listand et (MOTS 'PPDT') ; 274 : 'FINSI' ; 275 : 276 : Si TBT.'PPDT' ; 277 : TBT.'FCONS'= 'CONS2' ; 278 : Sinon ; 279 : TBT.'FCONS'= 'NOCONS' ; 280 : Finsi ; 281 : *-------------------------------------------------------- 282 : * ALGO = choix de l'algorithme implicite / explicite 283 : * Les choix possibles sont : 284 : * IMPL = semi implicite 285 : * EFM1 = semi explicite 286 : *-------------------------------------------------------- 287 : 288 : TBT.'ALGO' ='IMPL' ; 289 : TBT.'SPGNU'='CHAI' sommet ; 290 : TBT.'FEF' ='EF' ; 291 : TBT.'FIMPL'='IMPL' ; 292 : 'SI' ('EXIST' rxt 'ALGO') ; 293 : vertytab rxt 'ALGO' 'MOT' ; 294 : lxt=mots 'IMPL' 'EFM1' ; 295 : vertytab lxt (rxt.'ALGO') 'MOT' ; 296 : TBT.'ALGO' = rxt.'ALGO' ; 297 : Si (EGA TBT.'ALGO' 'EFM1') ; 298 : TBT.'SPGNU'='CHAI' centre ; 299 : TBT.'FEF' ='EFM1' ; 300 : TBT.'FIMPL'='EXPL' ; 301 : Finsi ; 302 : listand=listand et (MOTS 'ALGO') ; 303 : 'FINSI' ; 304 : 'MESS' '==> Algorithme choisi ' TBT.'ALGO' ; 305 : 306 : *-------------------------------------------------------- 307 : * DISCR = choix de la discrétisation 308 : * Les dicrétisations possibles sont : 309 : * LINE = éléments linéaires 310 : * MACRO = éléments MACRO (par DEFAUT) 311 : * QUAF = éléments quadratiques fluide 312 : *-------------------------------------------------------- 313 : 314 : TBT.'discr'='MACRO' ; 315 : 'SI' ('EXIST' rxt 'DISCR') ; 316 : vertytab rxt 'DISCR' 'MOT' ; 317 : lxt=mots 'MACRO' 'LINE' 'QUAF' ; 318 : vertytab lxt (rxt.'DISCR') 'MOT' ; 319 : TBT.'discr' = rxt.'DISCR' ; 320 : listand=listand et (MOTS 'DISCR') ; 321 : 'FINSI' ; 322 : 'MESS' '==> Discrétisation choisie pour la vitesse ' TBT.'discr' ; 323 : TBT.'discrH'=TBT.'discr' ; 324 : Si ((Ega TBT.'discr' 'QUAF') ET (EGA (VALE DIME) 3)) ; 325 : TBT.'discrH'='MACRO' ; 326 : Finsi ; 327 : 328 : *--------------------------------------------------------- 329 : * KPRE = position de l'inconnue de pression 330 : * possibilités : CENTRE avec algorithme de projection 331 : * CENTREP1 avec algorithme IMPLICIT 332 : *--------------------------------------------------------- 333 : 334 : TBT.'KPRE'='CENTRE' ; 335 : TBT.'MAXELPRE'=20000 ; 336 : 'SI' ('EXIST' rxt 'KPRE') ; 337 : vertytab rxt 'KPRE' 'MOT' ; 338 : lxt=table; lxt.'CENTRE'=' '; lxt.'CENTREP1'=' ' ; 339 : lxt.'MSOMMET'=' ' ; 340 : vertytab lxt (rxt.'KPRE') 'MOT' ; 341 : TBT.'KPRE' = rxt.'KPRE' ; 342 : listand=listand et (MOTS 'KPRE') ; 343 : 'FINSI' ; 344 : 'MESS' '==> Discrétisation choisie pour la pression ' TBT.'KPRE'; 345 : 346 : *--------------------------------------------------------- 347 : * Definition des grandeurs Physiques : 348 : * Rgi = constante des gaz pour chaque especes J/Kg/K 349 : * 350 : * Modif FP 26/06/00 351 : Rabs = 6.022*1.38066 ; 352 : TBT.'Rabs' = Rabs ; 353 : TBT.'Rgh2' = Rabs/2.E-3 ; TBT.'Mkh2' = 2.016E-3 ; 354 : TBT.'Rghe' = Rabs/4.E-3 ; TBT.'Mkhe' = 4.003E-3 ; 355 : TBT.'Rgo2' = Rabs/3.2E-2 ; TBT.'Mko2' = 32.00E-3 ; 356 : TBT.'Rgn2' = Rabs/2.8E-2 ; TBT.'Mkn2' = 28.02E-3 ; 357 : TBT.'Rgco2' = Rabs/4.4E-2 ; TBT.'Mkco2' = 44.01E-3 ; 358 : TBT.'Rgco' = Rabs/2.8E-2 ; TBT.'Mkco' = 28.01E-3 ; 359 : TBT.'Rgvap' = Rabs/1.8E-2 ; TBT.'Mkvap' = 18.02E-3 ; 360 : TBT.'Rgair' = 287.1 ; TBT.'Mkair' = 28.97E-3 ; 361 : TBT.'Rgar' = Rabs/4.E-2 ; TBT.'Mkar' = 39.94E-3 ; 362 : TBT.'Rgch4' = Rabs/16.E-3 ; TBT.'Mkch4' = 16.04E-3 ; 363 : TBT.'lamair'= 2.580e-2 ; 364 : 365 : TBT.'B0' = -4512.70756 ; TBT.'B1' = 4245585.9 ; 366 : TBT.'B2' = -1584624342.; 367 : 368 : * Ancienne mouture 369 : **TBT.'Rgh2' = 4156.5 ; 370 : **TBT.'Rghe' = 2078.25 ; 371 : **TBT.'Rgo2' = 259.8 ; 372 : **TBT.'Rgn2' = 296.9 ; 373 : **TBT.'Rgco2' = 188.9 ; 374 : **TBT.'Rgco' = 296.9 ; 375 : **TBT.'Rgvap' = 461.513 ; 376 : **TBT.'Rgair' = (0.8 '*' TBT.'Rgn2') '+' (0.2 '*' TBT.'Rgo2'); 377 : *-------------------------------------------------------- 378 : * Cpvap = Chaleur specifique a pression constante(J/Kg/K) 379 : * Hypothese de gaz parfait pour la vapeur d'eau 380 : * 381 : 'SI' ('EXIST' rxt 'Cpvap') ; 382 : vertytab rxt 'Cpvap' 'FLOTTANT' ; 383 : Cpvap = rxt.'Cpvap' ; 384 : 'SINON' ; 385 : Cpvap = 1700. ; 386 : 'MESS' '==> Cpvap mis par defaut a 1700 J/kg/K ' ; 387 : 'FINSI' ; 388 : TBT.'Cpvap'=Cpvap ; 389 : TBT.'Cvvap'=TBT.'Cpvap' - TBT.'Rgvap' ; 390 : *-------------------------------------------------------- 391 : * Lv = Chaleur Latente supposée constante (J/kg) 392 : * 393 : 'SI' ('EXIST' rxt 'Lv') ; 394 : vertytab rxt 'Lv' 'FLOTTANT' ; 395 : Lv = rxt.'Lv' ; 396 : 'SINON' ; 397 : Lv = 2.3E6 ; 398 : 'MESS' '==> Chaleur latente Lv mise par defaut a 2.3E6 J/kg' ; 399 : 'FINSI' ; 400 : TBT.'Lv'=Lv ; 401 : *-------------------------------------------------------- 402 : * DB = coefficient de diffusion Brownienne 403 : * 404 : 'SI' ('EXIST' rxt 'db') ; 405 : vertytab rxt 'db' 'FLOTTANT' ; 406 : db = rxt.'db' ; 407 : 'SINON' ; 408 : db = 1.000e-5 ; 409 : 'FINSI' ; 410 : TBT.'db'=db ; 411 : *---------------------------------------------------------- 412 : 413 : TBT.'Taxe' = FAUX ; 414 : Si (exist rxt 'axe') ; 415 : vertytab rxt 'axe' 'MAILLAGE' ; 416 : TBT.'Taxe' = VRAI ; 417 : Finsi ; 418 : 419 : *---------------------------------------------------------- 420 : * SCENARIO = on verifie existence de la TABLE necessaire 421 : * pour faire un calcul avec ENCEINTE 422 : * si il y a au moins une brèche 423 : 424 : TBT.'Tbreche' = FAUX ; 425 : Si (exist rxt 'breche') ; 426 : 427 : vertytab rxt 'breche' 'MAILLAGE' ; 428 : TBT.'Tbreche' = VRAI ; 429 : 430 : * listentr= listentr et (mots 'diru1' 'scenario') ; 431 : listentr.'diru1' ='' ; 432 : listentr.'scenario' ='' ; 433 : vertytab rxt 'diru1' 'POINT ' ; 434 : TBT.'diru1'=rxt.'diru1' ; 435 : 436 : vertytab rxt 'scenario' 'TABLE ' ; 437 : scn = rxt.'scenario' ; 438 : TBT.'Modinj'=mot 'A' ; 439 : Si(EXIST rxt.'scenario' 'Modinj') ; 440 : TBT.'Modinj'=rxt.'scenario'.'Modinj' ; 441 : Finsi ; 442 : Finsi ; 443 : 444 : TBT.'Tbreche2' = FAUX ; 445 : Si (exist rxt 'breche2') ; 446 : 447 : * pas de breche2 sans breche sinon pb ds les Qo 448 : vertytab rxt 'breche' 'MAILLAGE' ; 449 : vertytab rxt 'breche2' 'MAILLAGE' ; 450 : TBT.'Tbreche2' = VRAI ; 451 : 452 : * listentr= listentr et (mots 'diru2''scenario2') ; 453 : listentr.'diru2' ='' ; 454 : listentr.'scenario2' ='' ; 455 : vertytab rxt 'diru2' 'POINT ' ; 456 : TBT.'diru2'=rxt.'diru2' ; 457 : 458 : vertytab rxt 'scenario2' 'TABLE ' ; 459 : scn2 = rxt.'scenario2' ; 460 : TBT.'Modinj2'=mot 'A' ; 461 : Si(EXIST rxt.'scenario2' 'Modinj') ; 462 : TBT.'Modinj2'=rxt.'scenario2'.'Modinj' ; 463 : Finsi ; 464 : Finsi ; 465 : 466 : TBT.'Tbreche3' = FAUX ; 467 : Si (exist rxt 'breche3') ; 468 : 469 : * pas de breche3 sans breche sinon pb ds les Qo 470 : vertytab rxt 'breche' 'MAILLAGE' ; 471 : vertytab rxt 'breche2' 'MAILLAGE' ; 472 : vertytab rxt 'breche3' 'MAILLAGE' ; 473 : TBT.'Tbreche3' = VRAI ; 474 : 475 : * listentr= listentr et (mots 'diru3''scenario3') ; 476 : listentr.'diru3' ='' ; 477 : listentr.'scenario3' ='' ; 478 : vertytab rxt 'diru3' 'POINT ' ; 479 : TBT.'diru3'=rxt.'diru3' ; 480 : 481 : vertytab rxt 'scenario3' 'TABLE ' ; 482 : scn3 = rxt.'scenario3' ; 483 : TBT.'Modinj3'=mot 'A' ; 484 : Si(EXIST rxt.'scenario3' 'Modinj') ; 485 : TBT.'Modinj3'=rxt.'scenario3'.'Modinj' ; 486 : Finsi ; 487 : Finsi ; 488 : 489 : TBT.'Tsortie' = FAUX ; 490 : 491 : Si (exist rxt 'sortie') ; 492 : vertytab rxt 'sortie' 'MAILLAGE' ; 493 : TBT.'Tsortie' = VRAI ; 494 : Si (NON (exist rxt 'CORTEMP')) ; 495 : TBT.'CORTEMP' = FAUX ; 496 : Finsi ; 497 : Finsi ; 498 : 499 : *---------------------------------------------------------- 500 : * vtp : thermique paroi 501 : * 502 : 503 : TBT.'THERMP' = FAUX ; 504 : TBT.'THERCO' = FAUX ; 505 : TBT.'ECHEXT' = FAUX ; 506 : 507 : Si (Exist rxt 'THERMP') ; 508 : vertytab rxt 'THERMP' 'LOGIQUE' ; 509 : TBT.'THERMP' = rxt.'THERMP' ; 510 : Finsi ; 511 : 512 : Si (Exist rxt 'THERCO') ; 513 : vertytab rxt 'THERCO' 'LOGIQUE' ; 514 : TBT.'THERCO' = rxt.'THERCO' ; 515 : Si (TBT.'THERCO') ; 516 : TBT.'THERMP' = VRAI ; 517 : listand=listand et (MOTS 'THCO') ; 518 : Finsi ; 519 : Finsi ; 520 : 521 : Si TBT.'THERMP' ; 522 : * On vérifie que le modèle est bien renseigné 523 : vertytab rxt 'vtp' 'MAILLAGE' ; 524 : vertytab rxt 'ROCP' 'FLOTTANT' 'CHPOINT' ; 525 : vertytab rxt 'LAMBDA' 'FLOTTANT' 'CHPOINT' ; 526 : vertytab rxt 'Tp0' 'FLOTTANT' 'CHPOINT' ; 527 : vertytab rxt 'ECHAN' 'FLOTTANT' 'CHPOINT' ; 528 : 529 : listchp=listchp et (MOTS 'ROCP') ; 530 : listchp=listchp et (MOTS 'LAMBDA') ; 531 : listchp=listchp et (MOTS 'Tp0' ) ; 532 : listchp=listchp et (MOTS 'ECHAN' ) ; 533 : 534 : 'SI' ('NON' ('EXIST' tic 'ROCP')) ; 535 : tic.'ROCP'=rxt.'ROCP' ; 536 : 'FINSI' ; 537 : 'SI' ('NON' ('EXIST' tic 'LAMBDA')) ; 538 : tic.'LAMBDA'=rxt.'LAMBDA' ; 539 : 'FINSI' ; 540 : 'SI' ('NON' ('EXIST' tic 'Tp0')) ; 541 : tic.'Tp0' =rxt.'Tp0' ; 542 : 'FINSI' ; 543 : *? 'SI' ('NON' ('EXIST' tic 'KHW')) ; 544 : tic.'KHW' =rxt.'ECHAN' ; 545 : *? 'FINSI' ; 546 : 547 : Si (Exist rxt 'ECHEXT') ; 548 : vertytab rxt 'ECHEXT' 'LOGIQUE' ; 549 : TBT.'ECHEXT' = rxt.'ECHEXT' ; 550 : Finsi ; 551 : 552 : Si (TBT.'ECHEXT') ; 553 : vertytab rxt 'parext' 'MAILLAGE' ; 554 : vertytab rxt 'HEXT' 'FLOTTANT' 'CHPOINT' ; 555 : vertytab rxt 'TPEXT' 'FLOTTANT' 'CHPOINT' ; 556 : 557 : listchp=listchp et (MOTS 'HEXT' ) ; 558 : listchp=listchp et (MOTS 'TPEXT' ) ; 559 : 560 : 'SI' ('NON' ('EXIST' tic 'HEXT')) ; 561 : tic.'HEXT'=rxt.'HEXT' ; 562 : 'FINSI' ; 563 : 'SI' ('NON' ('EXIST' tic 'TPEXT')) ; 564 : tic.'TPEXT'=rxt.'TPEXT' ; 565 : 'FINSI' ; 566 : 567 : Finsi ; 568 : 569 : 'MESS' ' : ==> Scénario avec thermique paroi ' ; 570 : Sinon ; 571 : 'MESS' ' : ==> Scénario sans thermique paroi ' ; 572 : Finsi ; 573 : *---------------------------------------------------------- 574 : * TIMPi : Températures imposées 575 : * 576 : TBT.'TTIMP1' = FAUX ; 577 : TBT.'TTIMP2' = FAUX ; 578 : TBT.'TTIMP3' = FAUX ; 579 : Si (Exist rxt 'TIMP1') ; 580 : TBT.'TTIMP1' = VRAI ; 581 : vertytab rxt 'TIMP1' 'TABLE' ; 582 : vertytab (rxt.'TIMP1') 'MAILLAGE' 'MAILLAGE' ; 583 : vertytab (rxt.'TIMP1') 'ECHAN' 'FLOTTANT' ; 584 : vertytab (rxt.'TIMP1') 't' 'LISTREEL' ; 585 : vertytab (rxt.'TIMP1') 'TIMP' 'LISTREEL' ; 586 : Finsi ; 587 : 588 : Si (Exist rxt 'TIMP2') ; 589 : TBT.'TTIMP2' = VRAI ; 590 : vertytab rxt 'TIMP2' 'TABLE' ; 591 : vertytab (rxt.'TIMP2') 'MAILLAGE' 'MAILLAGE' ; 592 : vertytab (rxt.'TIMP2') 'ECHAN' 'FLOTTANT' ; 593 : vertytab (rxt.'TIMP2') 't' 'LISTREEL' ; 594 : vertytab (rxt.'TIMP2') 'TIMP' 'LISTREEL' ; 595 : Finsi ; 596 : 597 : Si (Exist rxt 'TIMP3') ; 598 : TBT.'TTIMP3' = VRAI ; 599 : vertytab rxt 'TIMP3' 'TABLE' ; 600 : vertytab (rxt.'TIMP3') 'MAILLAGE' 'MAILLAGE' ; 601 : vertytab (rxt.'TIMP3') 'ECHAN' 'FLOTTANT' 'CHPOIN' ; 602 : vertytab (rxt.'TIMP3') 't' 'LISTREEL' ; 603 : vertytab (rxt.'TIMP3') 'TIMP' 'LISTREEL' ; 604 : Finsi ; 605 : 606 : *---------------------------------------------------------- 607 : * ECHANP : Température du mur constante ou calculée via une procédure 608 : * personnelle 609 : TBT.'TECHANP' = FAUX ; 610 : Si (Exist rxt 'ECHANP') ; 611 : TBT.'TECHANP' = VRAI ; 612 : vertytab rxt 'ECHANP' 'TABLE' ; 613 : vertytab (rxt.'ECHANP') 'MAILLAGE' 'MAILLAGE' ; 614 : vertytab (rxt.'ECHANP') 'ECHAN' 'FLOTTANT' ; 615 : vertytab (rxt.'ECHANP') 'TMUR' 'FLOTTANT' ; 616 : Finsi ; 617 : 618 : * On verifie les données de la table RECOMB 619 : * pour définir des recombineurs 620 : TBT.'TRECOMB' = FAUX ; 621 : 'SI' ('EXIST' rxt 'RECOMB') ; 622 : vertytab rxt 'RECOMB' 'TABLE ' ; 623 : NBREC = 'DIME' rxt.'RECOMB' ; 624 : 625 : 'REPETER' BCL1 NBREC ; 626 : vertytab (rxt.'RECOMB' . &BCL1) 'PAREXT' 'MAILLAGE' ; 627 : vertytab (rxt.'RECOMB' . &BCL1) 'ENTREE' 'MAILLAGE' ; 628 : vertytab (rxt.'RECOMB' . &BCL1) 'SORTIE' 'MAILLAGE' ; 629 : vertytab (rxt.'RECOMB' . &BCL1) 'direntr' 'POINT ' ; 630 : vertytab (rxt.'RECOMB' . &BCL1) 'dirsort' 'POINT ' ; 631 : 'FIN' BCL1 ; 632 : TBT.'TRECOMB' = VRAI ; 633 : TBT.'RECOMB' = rxt.'RECOMB' ; 634 : 'FINSI' ; 635 : *---------------------------------------------------------- 636 : * MODTURB = Mot identifiant le modèle de turbulence 637 : * NUTURB ou LMEL ou KEPSILON ou NUTCONS 638 : 639 : Si (Exist rxt 'MODTURB') ; 640 : 641 : vertytab rxt 'MODTURB' 'MOT' ; 642 : TBT.'MODTURB'=rxt.'MODTURB' ; 643 : 644 : * On vérifie que le modèle est bien identifié 645 : lxt=table ; 646 : lxt.'NUTURB' = ' ' ; 647 : lxt.'LMEL' = ' ' ; 648 : lxt.'KEPSILON' = ' ' ; 649 : vertytab lxt rxt.'MODTURB' 'MOT' ; 650 : 651 : * Nu turbulent constant 652 : Si(Ega rxt.'MODTURB' 'NUTURB') ; 653 : * listentr=listentr et (mots 'NUT') ; 654 : listentr.'NUT' ='' ; 655 : vertytab rxt 'NUT' 'FLOTTANT' ; 656 : TBT.'NUT' = rxt.'NUT' ; 657 : Finsi ; 658 : 659 : * Longueur de mélange 660 : Si(Ega rxt.'MODTURB' 'LMEL' ) ; 661 : * listentr=listentr et (mots 'LMEL') ; 662 : listentr.'LMEL' ='' ; 663 : vertytab rxt 'LMEL' 'FLOTTANT' ; 664 : TBT.'LMEL' = rxt.'LMEL' ; 665 : Finsi ; 666 : 667 : 'MESS' ' : ==> Scenario avec modèle de turbulence ' 668 : rxt.'MODTURB' ; 669 : Sinon ; 670 : 'MESS' ' : ==> Scenario sans modèle de turbulence '; 671 : Finsi ; 672 : *---------------------------------------------------------- 673 : * FPAROI = logique identifiant l'utilisation de fonctions 674 : * de paroi 675 : 676 : TBT.'FPAROI' = FAUX ; 677 : 'SI' ('EXIST' rxt 'FPAROI') ; 678 : vertytab rxt 'FPAROI' 'LOGIQUE' ; 679 : TBT.'FPAROI' = rxt.'FPAROI' ; 680 : 'FINSI' ; 681 : 682 : 'SI' TBT.'FPAROI' ; 683 : 'MESS' ' : ==> Scenario avec fonction de paroi '; 684 : * listentr=listentr et (mots 'YP') ; 685 : listentr.'YP' ='' ; 686 : vertytab rxt 'YP' 'FLOTTANT' ; 687 : TBT.'YP'=rxt.'YP' ; 688 : 'FINSI' ; 689 : 690 : *---------------------------------------------------------- 691 : * MODCOND = mot identifiant la formulation du flux de 692 : * condensation 693 : TBT.'MODCOND' = 'CHIL0' ; 694 : 'SI' ('EXIST' rxt 'MODCOND') ; 695 : vertytab rxt 'MODCOND' 'MOT' ; 696 : TBT.'MODCOND' = rxt.'MODCOND' ; 697 : 698 : * On vérifie que le modèle est bien identifié 699 : lxt=table ; 700 : lxt.'CHIL0' = ' ' ; 701 : lxt.'CHIL1' = ' ' ; 702 : vertytab lxt rxt.'MODCOND' 'MOT' ; 703 : 704 : 'FINSI' ; 705 : 'MESS' ' : ==> Scénario avec flux de condensation ' 706 : TBT.'MODCOND' ; 707 : *---------------------------------------------------------- 708 : * PERSO = logique identifiant une procédure PERSO 709 : * dans le cas de calcul 710 : 711 : TBT.'PERSO' = FAUX ; 712 : 'SI' ('EXIST' rxt 'PERSO') ; 713 : vertytab rxt 'PERSO' 'LOGIQUE' ; 714 : TBT.'PERSO' = rxt.'PERSO' ; 715 : 'FINSI' ; 716 : 717 : 'SI' TBT.'PERSO' ; 718 : 'MESS' ' : ==> Scénario avec exécution d une ' ; 719 : 'MESS' ' procédure personnelle ' ; 720 : 721 : * listentr = listentr et (mots 'PRCPERSO' 'TABPERSO') ; 722 : listentr.'PRCPERSO' ='' ; 723 : listentr.'TABPERSO' ='' ; 724 : 725 : vertytab rxt 'PRCPERSO' 'MOT' ; 726 : TBT.'PRCPERSO' = rxt.'PRCPERSO' ; 727 : vertytab rxt 'TABPERSO' 'TABLE' ; 728 : TBT.'TABPERSO' = rxt.'TABPERSO' ; 729 : 'FINSI' ; 730 : 731 : *---------------------------------------------------------- 732 : * VAPEUR = logique identifiant l'utilisation de vapeur 733 : * dans le cas de calcul 734 : 735 : TBT.'VAPEUR' = FAUX ; 736 : 'SI' ('EXIST' rxt 'VAPEUR') ; 737 : vertytab rxt 'VAPEUR' 'LOGIQUE' ; 738 : TBT.'VAPEUR' = rxt.'VAPEUR' ; 739 : 'FINSI' ; 740 : 741 : 'SI' TBT.'VAPEUR' ; 742 : vertytab rxt 'Yvap0' 'FLOTTANT' ; 743 : listentr.'Yvap0' ='' ; 744 : listentr.'Mliq0' ='' ; 745 : listentr.'Lv' ='' ; 746 : listentr.'Cpvap' ='' ; 747 : 'FINSI' ; 748 : 749 : *----------------------------------------------------------- 750 : * Modif. FP 18/05/00 751 : * Introduction d'une nouvelle variable logique qui est vrai lorsque la 752 : * condensation en masse est prise en compte. 753 : 754 : * CONDMAS = Mot identifiant le modèle de condensation en masse 755 : * CMAS0 756 : 757 : TBT.'CONDMAS' = FAUX ; 758 : TBT.'MODCMAS' ='????' ; 759 : listentr.'CONDMAS' ='' ; 760 : listentr.'HMAS' ='' ; 761 : listentr.'WCM' ='' ; 762 : *'SI' ('EXIST' rxt 'CONDMAS') ; 763 : * vertytab rxt 'CONDMAS' 'MOT' ; 764 : * TBT.'CONDMAS' = VRAI ; 765 : ** On vérifie que le modèle est bien identifié 766 : * lxt=table ; 767 : * lxt.'CMAS0' = ' ' ; 768 : * vertytab lxt rxt.'CONDMAS' 'MOT' ; 769 : * TBT.'MODCMAS' = rxt.'CONDMAS' ; 770 : * Si ('EGA' rxt.'CONDMAS' 'CMAS0') ; 771 : * Si ('EXIST' rxt 'HMAS') ; 772 : * vertytab rxt 'HMAS' 'FLOTTANT' ; 773 : * TBT.'HMAS'=rxt.'HMAS' ; 774 : * Sinon ; 775 : * TBT.'HMAS'=0. ; 776 : * Finsi ; 777 : * Si ('EXIST' rxt 'WCM') ; 778 : * vertytab rxt 'WCM' 'FLOTTANT' ; 779 : * TBT.'WCM'=rxt.'WCM' ; 780 : * Sinon ; 781 : * TBT.'WCM'=0. ; 782 : * Finsi ; 783 : * Finsi ; 784 : *'FINSI' ; 785 : *----------------------------------------------------------- 786 : * ASPER logique indiquant qu'on traite l'aspersion 787 : 788 : 'SI' TBT.'VAPEUR' ; 789 : listentr.'ASPER' ='' ; 790 : 'FINSI' ; 791 : 792 : TBT.'ASPER' = FAUX ; 793 : 'SI' ('EXIST' rxt 'ASPER') ; 794 : vertytab rxt 'ASPER' 'LOGIQUE' ; 795 : TBT.'ASPER' = rxt.'ASPER' ; 796 : 'FINSI' ; 797 : 798 : 'SI' (TBT.'ASPER') ; 799 : listentr.'scenasp' ='' ; 800 : vertytab rxt 'scenasp' 'TABLE ' ; 801 : scn = rxt.'scenasp' ; 802 : 803 : vertytab rxt 'aspinj' 'MAILLAGE' ; 804 : listentr.'aspinj' ='' ; 805 : vertytab rxt 'toitf' 'MAILLAGE' ; 806 : listentr.'toitf' ='' ; 807 : vertytab rxt 'rod' 'FLOTTANT' ; 808 : listentr.'rod' ='' ; 809 : vertytab rxt 'Cpd' 'FLOTTANT' ; 810 : listentr.'Cpd' ='' ; 811 : 'FINSI' ; 812 : 813 : *---------------------------------------------------------- 814 : * Vérification des incondensables utilisés 815 : * Les incondensables possibles sont H2, HE, O2, N2, CO, CO2 816 : * 817 : *--------------------------------------------------------- 818 : LINCOND = 'MOTS' ; 819 : TBT.'TH2' = FAUX ; 820 : TBT.'THE' = FAUX ; 821 : TBT.'TO2' = FAUX ; 822 : TBT.'TN2' = FAUX ; 823 : TBT.'TCO2'= FAUX ; 824 : *? TCO2AIR 825 : TBT.'TAIR'= FAUX ; 826 : TBT.'TCO' = FAUX ; 827 : 828 : 'SI' ('EXIST' rxt 'H2' ) ; 829 : vertytab rxt 'H2' 'LOGIQUE' ; 830 : TBT.'TH2' = rxt.'H2' ; 831 : 'SI' TBT.'TH2' ; 832 : LINCOND=LINCOND 'ET' ('MOTS' 'H2') ; 833 : vertytab rxt 'Yh20' 'FLOTTANT' ; 834 : listentr.'Yh20' ='' ; 835 : 'SI' ('EXIST' rxt 'scenario') ; 836 : vertytab (rxt.'scenario') 'qh2' 'LISTREEL' ; 837 : 'FINSI' ; 838 : 'FINSI' ; 839 : 'FINSI' ; 840 : 841 : 'SI' ('EXIST' rxt 'HE' ) ; 842 : vertytab rxt 'HE' 'LOGIQUE' ; 843 : TBT.'THE' = rxt.'HE' ; 844 : 'SI' TBT.'THE' ; 845 : LINCOND=LINCOND 'ET' ('MOTS' 'HE') ; 846 : vertytab rxt 'Yhe0' 'FLOTTANT' ; 847 : listentr.'Yhe0' ='' ; 848 : 'SI' ('EXIST' rxt 'scenario') ; 849 : vertytab (rxt.'scenario') 'qhe' 'LISTREEL' ; 850 : 'FINSI' ; 851 : 'FINSI' ; 852 : 'FINSI' ; 853 : 854 : 'SI' ('EXIST' rxt 'O2' ) ; 855 : vertytab rxt 'O2' 'LOGIQUE' ; 856 : TBT.'TO2' = rxt.'O2' ; 857 : 'SI' TBT.'TO2' ; 858 : LINCOND=LINCOND 'ET' ('MOTS' 'O2') ; 859 : vertytab rxt 'Yo20' 'FLOTTANT' ; 860 : listentr.'Yo20' ='' ; 861 : 'SI' ('EXIST' rxt 'scenario') ; 862 : vertytab (rxt.'scenario') 'qo2' 'LISTREEL' ; 863 : 'FINSI' ; 864 : 'FINSI' ; 865 : 'FINSI' ; 866 : 867 : 'SI' ('EXIST' rxt 'N2' ) ; 868 : vertytab rxt 'N2' 'LOGIQUE' ; 869 : TBT.'TN2' = rxt.'N2' ; 870 : 'SI' TBT.'TN2' ; 871 : LINCOND=LINCOND 'ET' ('MOTS' 'N2') ; 872 : vertytab rxt 'Yn20' 'FLOTTANT' ; 873 : listentr.'Yn20' ='' ; 874 : 'SI' ('EXIST' rxt 'scenario') ; 875 : vertytab (rxt.'scenario') 'qn2' 'LISTREEL' ; 876 : 'FINSI' ; 877 : 'FINSI' ; 878 : 'FINSI' ; 879 : 880 : 'SI' ('EXIST' rxt 'CO2' ) ; 881 : vertytab rxt 'CO2' 'LOGIQUE' ; 882 : TBT.'TCO2' = rxt.'CO2' ; 883 : 'SI' TBT.'TCO2' ; 884 : LINCOND=LINCOND 'ET' ('MOTS' 'CO2') ; 885 : vertytab rxt 'Yco20' 'FLOTTANT' ; 886 : listentr.'Yco20' ='' ; 887 : 'SI' ('EXIST' rxt 'scenario') ; 888 : vertytab (rxt.'scenario') 'qco2' 'LISTREEL' ; 889 : 'FINSI' ; 890 : 'FINSI' ; 891 : 'FINSI' ; 892 : 893 : *? TCO2AIR 894 : 'SI' ('EXIST' rxt 'AIR' ) ; 895 : vertytab rxt 'AIR' 'LOGIQUE' ; 896 : TBT.'TAIR' = rxt.'AIR' ; 897 : 'SI' TBT.'TAIR' ; 898 : LINCOND=LINCOND 'ET' ('MOTS' 'AIR') ; 899 : *? non vertytab rxt 'Yair0' 'FLOTTANT' ; 900 : *? non listentr.'Yair0' ='' ; 901 : 'SI' ('EXIST' rxt 'scenario') ; 902 : vertytab (rxt.'scenario') 'qair' 'LISTREEL' ; 903 : 'FINSI' ; 904 : 'FINSI' ; 905 : 'FINSI' ; 906 : 907 : 'SI' ('EXIST' rxt 'CO' ) ; 908 : vertytab rxt 'CO' 'LOGIQUE' ; 909 : TBT.'TCO' = rxt.'CO' ; 910 : 'SI' TBT.'TCO' ; 911 : LINCOND=LINCOND 'ET' ('MOTS' 'CO') ; 912 : vertytab rxt 'Yco0' 'FLOTTANT' ; 913 : listentr.'Yco0' ='' ; 914 : 'SI' ('EXIST' rxt 'scenario') ; 915 : vertytab (rxt.'scenario') 'qco' 'LISTREEL' ; 916 : 'FINSI' ; 917 : 'FINSI' ; 918 : 'FINSI' ; 919 : 920 : *--------------------------------------------------------- 921 : * Affichage de la liste des incondensables 922 : *--------------------------------------------------------- 923 : 924 : TINCOND= 'NON' ( 'EGA' ('DIME' LINCOND ) 0 ) ; 925 : 'SI' TINCOND ; 926 : 'MESS' ' Liste des incondensables : ' ; 927 : 'LIST' LINCOND ; 928 : 'SINON' ; 929 : 'MESS' ' Il n y a pas d incondensable hormis AIR ' ; 930 : 'FINSI' ; 931 : 932 : 'SI' TBT.'TRECOMB' ; 933 : LOG0 = TBT.'TH2' 'ET' TBT.'VAPEUR' 'ET' TBT.'TN2' 'ET' TBT.'TO2' ; 934 : 'SI' ('NON' LOG0) ; 935 : 'MESS' 'RECOMBINEUR ==> H2, H20, O2 et N2 par defaut !' ; 936 : 'ERRE' 5 ; 937 : 'FINSI' ; 938 : 'FINSI' ; 939 : 940 : *--------------------------------------------------------- 941 : * Initialisation des paramètres du modèle 942 : * nimpr = Niveau impression lors de la résolution 943 : * 0 = Pas d'impression 944 : * epsi = tolérance pour les maillages (OPERATEUR ELIM) 945 : * FRPREC : ENTIER fréquence de recalcul du préconditionnement 946 : *--------------------------------------------------------- 947 : 948 : TBT.'nimpr'=0 ; 949 : TBT.'IMPR' =2 ; 950 : 'SI' ('EXIST' rxt 'IMPR' ) ; 951 : vertytab rxt 'IMPR' 'ENTIER' ; 952 : TBT.'IMPR' =rxt.'IMPR' ; 953 : 'FINSI' ; 954 : 955 : GEO.'epsi'= 1.e-4 ; 956 : 'SI' ('EXIST' rxt 'epsi') ; 957 : vertytab rxt 'epsi' 'FLOTTANT' ; 958 : GEO.'epsi'= rxt.'epsi' ; 959 : 'FINSI' ; 960 : 961 : 'SI' ( 'EXIST' rxt 'FRPREC' ) ; 962 : vertytab rxt 'FRPREC' 'ENTIER' ; 963 : TBT.'FRPREC' = rxt.'FRPREC' ; 964 : 'FINSI' ; 965 : *----------------------------------------------------------------------- 966 : * Controle des entrées de la table RXT et de la compatibilité 967 : * entre les entrées. 968 : *----------------------------------------------------------------------- 969 : * On vérifie que les entrées de RXT sont autorisées 970 : 971 : 972 : lindex=index rxt ; 973 : n=dime lindex; 974 : 975 : repeter Bver n ; 976 : *mess lindex.&bver ; 977 : tae=type lindex.&bver ; 978 : Si (non (ega tae 'MOT')); 979 : Mess ' Il y a un indice qui n est pas de type mot dans la table '; 980 : erreur 21 ; 981 : Finsi ; 982 : Si (non (exist listentr lindex.&bver)) ; 983 : Mess ; 984 : Mess '**************************************************************'; 985 : Mess ; 986 : Mess ' Il y a un indice non reconnu : ' lindex.&bver 987 : ' dans la table '; 988 : Mess ; 989 : Mess '**************************************************************'; 990 : 991 : erreur 21 ; 992 : Finsi ; 993 : Fin Bver ; 994 : 995 : n=dime listand ; 996 : Si (non (ega n 0)); 997 : Mess ; 998 : Mess '**************************************************************'; 999 : Mess ; 1000 : Mess ' Il y a des options surchargées dont voici la liste : ' ; 1001 : list listand ; 1002 : Si (exist listand 'THCO') ; 1003 : Mess ; 1004 : Mess 'Commentaire:' ; 1005 : mess 'THCO: le couplage thermique paroi/thermohydraulique' 1006 : ' est implicite' ; 1007 : Finsi ; 1008 : Mess ; 1009 : Mess '**************************************************************'; 1010 : Finsi ; 1011 : 1012 : * Vérification de Compatibilités 1013 : 1014 : 'SI'(('EGA' TBT.'ALGO' 'EFM1') et ('EGA' TBT.'discr' 'QUAF')) ; 1015 : vertytab 'ERREUR' 'ERREUR' 1016 : (CHAI 'Options QUAF et EFM1 incompatibles'); 1017 : 'FINSI' ; 1018 : 1019 : 1020 : 'SI'(Exist TBT 'MODTURB') ; 1021 : 'SI'( (Ega TBT.'MODTURB' 'KEPSILON') et (non('EGA' TBT.'ALGO' 'EFM1'))); 1022 : vertytab 'ERREUR' 'ERREUR' 1023 : (CHAI 'L option EFM1 est nécessaire pour le modèle K-epsilon') ; 1024 : 'FINSI' ; 1025 : 'FINSI' ; 1026 : 1027 : 'SI'(TBT.'Tsortie' et TBT.'CORTEMP') ; 1028 : vertytab 'ERREUR' 'ERREUR' 1029 : (CHAI 'Ecoulement ouvert (Sortie) incompatible avec CORTEMP VRAI') ; 1030 : 'FINSI' ; 1031 : 1032 : 1033 : 1034 : 1035 : *----------------------------------------------------------------------- 1036 : *----------------------------------------------------------------------- 1037 : *----------------------------------------------------------------------- 1038 : * ------------- FIN DES INITIALISATIONS ET VERIFICATIONS --------------- 1039 : *----------------------------------------------------------------------- 1040 : *----------------------------------------------------------------------- 1041 : *----------------------------------------------------------------------- 1042 : 1043 : 'SI' (TPREPA); 1044 : PREPAENC rxt TBT GEO TIC ; 1045 : 'FINSI' ; 1046 : 1047 : rxt.'TBT'=TBT; 1048 : 1049 : *============================================================== 1050 : GEO=rxt.'GEO' ; 1051 : 1052 : * On depile un certain nombre de renseignements sur le 1053 : * cas test a traiter 1054 : 1055 : idim = vale 'DIME'; 1056 : DIM3D=FAUX ; 1057 : Si(EGA idim 3); 1058 : DIM3D=VRAI ; 1059 : Finsi ; 1060 : 1061 : * B/ Les maillages ou modeles 1062 : 1063 : vtf = rxt.'vtf' ; 1064 : 'SI' (TBT.'Tbreche'); 1065 : breche = rxt.'breche' ; 1066 : brechec=GEO.'brechec' ; 1067 : brechei=GEO.'brechei' ; 1068 : $breche=GEO.'$breche' ; 1069 : Sbreche = GEO.'Sbreche' ; 1070 : 'FINSI'; 1071 : 'SI' (TBT.'Tbreche2'); 1072 : breche2=rxt.'breche2' ; 1073 : brech2c=GEO.'brech2c' ; 1074 : brech2i=GEO.'brech2i' ; 1075 : $breche2=GEO.'$breche2' ; 1076 : Sbreche2= GEO.'Sbreche2' ; 1077 : 'FINSI'; 1078 : 'SI' (TBT.'Tbreche3'); 1079 : breche3=rxt.'breche3' ; 1080 : brech3c=GEO.'brech3c' ; 1081 : brech3i=GEO.'brech3i' ; 1082 : $breche3=GEO.'$breche3' ; 1083 : Sbreche3= GEO.'Sbreche3' ; 1084 : 'FINSI'; 1085 : 1086 : $vtf = GEO.'$vtf' ; 1087 : 'SI' TBT.'THERMP' ; 1088 : $vtp = GEO.'$vtp' ; 1089 : $paroif = GEO.'$paroif' ; 1090 : paroif = GEO.'paroif' ; 1091 : 'FINSI' ; 1092 : 1093 : Diag = 'DOMA' $vtf 'XXDIAGSI' ; 1094 : 1095 : Volvtf = 'DOMA' $vtf 'VOLUME' ; 1096 : VTotal = GEO.'VTotal' ; 1097 : 1098 : 1099 : 'MESS' ' ************* GEOMETRIE *********************' ; 1100 : 'MESS' ' **********************************************' ; 1101 : 'DOMA' $vtf 'IMPR' ; 1102 : 'SI' TBT.'THERMP' ; 1103 : 'DOMA' $vtp 'IMPR' ; 1104 : 'FINSI' ; 1105 : 1106 : 1107 : * C/ Les chaines de caracteres 1108 : 1109 : *? DISCR = TBT.'discr' ; 1110 : tic = rxt.'TIC' ; 1111 : 1112 : * D/ Les données physiques 1113 : 1114 : Cpvap = TBT.'Cpvap' ; 1115 : Lv = TBT.'Lv' ; 1116 : Rgh2 = TBT.'Rgh2' ; 1117 : Rghe = TBT.'Rghe' ; 1118 : Rgo2 = TBT.'Rgo2' ; 1119 : Rgn2 = TBT.'Rgn2' ; 1120 : Rgco2 = TBT.'Rgco2' ; 1121 : Rgco = TBT.'Rgco' ; 1122 : Rgvap = TBT.'Rgvap' ; 1123 : Rgair = TBT.'Rgair' ; 1124 : 1125 : 'MESS' ' ************* DISCRETISATION *****************' ; 1126 : 'MESS' ' **********************************************' ; 1127 : 'MESS' ' Discr = ' tbt.'discr' ' Pression = ' tbt.'KPRE' ; 1128 : 'MESS' ' Algo = ' algo ' Niveau impression = ' nimpr ; 1129 : 'MESS' ' Pas de temps DT = ' rxt.'DT0' ; 1130 : 1131 : * Calcul Cp (Tfm) 1132 : ndl = 'DIME' (tic.'Tfm') ; 1133 : Tfm = 'EXTR' (tic.'Tfm') ndl ; 1134 : * CALCP initialisation ? 1135 : Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair = CALCP Tfm ; 1136 : Cpvap = TBT.'Cpvap' ; 1137 : 1138 : * Calcul Mu (Tkm) 1139 : Tkm=Tfm '+' 273.15 ; 1140 : Muh2 Muhe Muo2 Mun2 Muco2 Muco Muvap Muair = CALMU Tkm ; 1141 : 1142 : 'MESS' ' ******* PROPRIETES PHYSIQUES *****************' ; 1143 : 'MESS' ' **********************************************' ; 1144 : 1145 : 'MESS' ' ***** Air ***** ' ; 1146 : Cvair = Cpair '-' Rgair ; 1147 : 'MESS' '==> Cpair = ' Cpair 'J/kg/K' ; 1148 : 'MESS' '==> Cvair = ' Cvair 'J/kg/K' ; 1149 : 'MESS' '==> Rgair = ' Rgair 'J/kg/K' ; 1150 : 1151 : 'SI' TBT.'VAPEUR' ; 1152 : 'MESS' ' ***** Vapeur H2o ***** ' ; 1153 : Cvvap = Cpvap '-' Rgvap ; 1154 : 'MESS' '==> Cpvap = ' Cpvap 'J/kg/K' ; 1155 : 'MESS' '==> Cvvap = ' Cvvap 'J/kg/K' ; 1156 : 'MESS' '==> Rgvap = ' Rgvap 'J/kg/K' ; 1157 : 'MESS' '==> Lv = ' Lv 'J/kg/K' ; 1158 : 'FINSI' ; 1159 : 1160 : 'SI' TBT.'THE' ; 1161 : 'MESS' ' ***** He ***** ' ; 1162 : Cvhe = Cphe '-' Rghe ; 1163 : 'MESS' '==> Cphe = ' Cphe 'J/kg/K' ; 1164 : 'MESS' '==> Cvhe = ' Cvhe 'J/kg/K' ; 1165 : 'MESS' '==> Rghe = ' Rghe 'J/kg/K' ; 1166 : 'FINSI' ; 1167 : 1168 : 'SI' TBT.'TH2' ; 1169 : 'MESS' ' ***** H2 ***** ' ; 1170 : Cvh2 = Cph2 '-' Rgh2 ; 1171 : 'MESS' '==> Cph2 = ' Cph2 'J/kg/K' ; 1172 : 'MESS' '==> Cvh2 = ' Cvh2 'J/kg/K' ; 1173 : 'MESS' '==> Rgh2 = ' Rgh2 'J/kg/K' ; 1174 : 'FINSI' ; 1175 : 1176 : 'SI' TBT.'TO2' ; 1177 : 'MESS' ' ***** O2 ***** ' ; 1178 : Cvo2 = Cpo2 '-' Rgo2 ; 1179 : 'MESS' '==> Cpo2 = ' Cpo2 'J/kg/K' ; 1180 : 'MESS' '==> Cvo2 = ' Cvo2 'J/kg/K' ; 1181 : 'MESS' '==> Rgo2 = ' Rgo2 'J/kg/K' ; 1182 : 'FINSI' ; 1183 : 1184 : 'SI' TBT.'TN2' ; 1185 : 'MESS' ' ***** N2 ***** ' ; 1186 : Cvn2 = Cpn2 '-' Rgn2 ; 1187 : 'MESS' '==> Cpn2 = ' Cpn2 'J/kg/K' ; 1188 : 'MESS' '==> Cvn2 = ' Cvn2 'J/kg/K' ; 1189 : 'MESS' '==> Rgn2 = ' Rgn2 'J/kg/K' ; 1190 : 'FINSI' ; 1191 : 1192 : 'SI' TBT.'TCO2' ; 1193 : 'MESS' ' ***** CO2 ***** ' ; 1194 : Cvco2 = Cpco2 '-' Rgco2 ; 1195 : 'MESS' '==> Cpco2 = ' Cpco2 'J/kg/K' ; 1196 : 'MESS' '==> Cvco2 = ' Cvco2 'J/kg/K' ; 1197 : 'MESS' '==> Rgco2 = ' Rgco2 'J/kg/K' ; 1198 : 'FINSI' ; 1199 : 1200 : *? TCO2AIR 1201 : 'SI' TBT.'TAIR' ; 1202 : 'MESS' ' ***** AIR ***** ' ; 1203 : Cvair = Cpair '-' Rgair ; 1204 : 'MESS' '==> Cpair = ' Cpair 'J/kg/K' ; 1205 : 'MESS' '==> Cvair = ' Cvair 'J/kg/K' ; 1206 : 'MESS' '==> Rgair = ' Rgair 'J/kg/K' ; 1207 : 'FINSI' ; 1208 : 1209 : 'SI' TBT.'TCO' ; 1210 : 'MESS' ' ***** CO ***** ' ; 1211 : Cvco = Cpco '-' Rgco ; 1212 : 'MESS' '==> Cpco = ' Cpco 'J/kg/K' ; 1213 : 'MESS' '==> Cvco = ' Cvco 'J/kg/K' ; 1214 : 'MESS' '==> Rgco = ' Rgco 'J/kg/K' ; 1215 : 'FINSI' ; 1216 : 1217 : 'MESS' ' *************** ' ; 1218 : 'MESS' ' Diffusion browniene db = ' db ; 1219 : 'SI' TBT.'THERMP' ; 1220 : 'SI' ('EGA' ('TYPE' tic.'KHW') 'CHPOINT') ; 1221 : 'MESS' ' Coefficient d echange ECHAN (Min/Max)= ' 1222 : ('MINI' tic.'KHW') ' ' ('MAXI' tic.'KHW') 'W/m2/K' ; 1223 : 'SINON' ; 1224 : 'MESS' ' Coefficient d echange ECHAN= ' tic.'KHW' 'W/m2/K' ; 1225 : 'FINSI' ; 1226 : 'FINSI' ; 1227 : 1228 : 1229 : 'MESS' '*** Bilans initiaux 0D-MultiD ***' ; 1230 : 'MESS' '_____________________________________________________________'; 1231 : 'MESS' '---------Bilans 0D ____________________Bilans MultiD ________' 1232 : '_____________ Erreur Relative (%)'; 1233 : *---- Bilan masse 1234 : Rhom = 'EXTR' (tic.'Rhom') ndl ; 1235 : M0D = Rhom '*' Vtotal ; 1236 : MMD = 'SOMT' (Diag '*' rxt.'TIC'.'RHO') ; 1237 : 1238 : 'SI' TBT.'VAPEUR' ; 1239 : Rhomv = 'EXTR' (tic.'Rhomv') ndl ; 1240 : MV0D = Rhomv '*' Vtotal ; 1241 : MVMD = 'SOMT' (Diag '*' rxt.'TIC'.'ROVP') ; 1242 : Rhomvg = 'EXTR' (tic.'Rhomvg') ndl ; 1243 : *? Rliqm = 'EXTR' (tic.'Rliqm') ndl ; 1244 : *? Rbrom = 'EXTR' (tic.'Rbrom') ndl ; 1245 : MVG0D = Rhomvg '*' Vtotal ; 1246 : MVGMD = 'SOMT' (Diag '*' rxt.'TIC'.'RVP') ; 1247 : 'FINSI' ; 1248 : 'SI' TBT.'THE' ; 1249 : Rhomhe = 'EXTR' (tic.'Rhomhe') ndl ; 1250 : MHE0D = Rhomhe '*' Vtotal ; 1251 : MHEMD = 'SOMT' (Diag '*' rxt.'TIC'.'RHE') ; 1252 : 'FINSI' ; 1253 : 'SI' TBT.'TH2' ; 1254 : Rhomh2 = 'EXTR' (tic.'Rhomh2') ndl ; 1255 : MH20D = Rhomh2 '*' Vtotal ; 1256 : MH2MD = 'SOMT' (Diag '*' rxt.'TIC'.'RH2') ; 1257 : 'FINSI' ; 1258 : 'SI' TBT.'TO2' ; 1259 : Rhomo2 = 'EXTR' (tic.'Rhomo2') ndl ; 1260 : MO20D = Rhomo2 '*' Vtotal ; 1261 : MO2MD = 'SOMT' (Diag '*' rxt.'TIC'.'RO2') ; 1262 : 'FINSI' ; 1263 : 'SI' TBT.'TN2' ; 1264 : Rhomn2 = 'EXTR' (tic.'Rhomn2') ndl ; 1265 : MN20D = Rhomn2 '*' Vtotal ; 1266 : MN2MD = 'SOMT' (Diag '*' rxt.'TIC'.'RN2') ; 1267 : 'FINSI' ; 1268 : 'SI' TBT.'TCO' ; 1269 : Rhomco = 'EXTR' (tic.'Rhomco') ndl ; 1270 : MCO0D = Rhomco '*' Vtotal ; 1271 : MCOMD = 'SOMT' (Diag '*' rxt.'TIC'.'RCO') ; 1272 : 'FINSI' ; 1273 : 'SI' TBT.'TCO2' ; 1274 : Rhomco2 = 'EXTR' (tic.'Rhomco2') ndl ; 1275 : MCO20D = Rhomco2 '*' Vtotal ; 1276 : MCO2MD = 'SOMT' (Diag '*' rxt.'TIC'.'RCO2') ; 1277 : 'FINSI' ; 1278 : *? TCO2AIR 1279 : 'SI' TBT.'TAIR' ; 1280 : Rhomair = 'EXTR' (tic.'Rhomair') ndl ; 1281 : MAIR0D = Rhomair '*' Vtotal ; 1282 : MAIRMD = 'SOMT' (Diag '*' rxt.'TIC'.'RAIR') ; 1283 : 'FINSI' ; 1284 : 1285 : *---- Bilan energie 1286 : 1287 : Remn = 'EXTR' (tic.'Remn') ndl ; 1288 : E0D = Remn ; 1289 : 'SI' ( 'EGA' ndl 1 ) ; 1290 : Cvm = 'EXTR' (tic.'Cvm') ndl ; 1291 : 'SINON' ; 1292 : Cvm = 'EXTR' (tic.'Cvm') (ndl-1) ; 1293 : 'FINSI' ; 1294 : EMD = ( 'SOMT' (Diag '*' Cvm '*' rxt.'TIC'.'RHO' '*' 1295 : (rxt.'TIC'.'TF' '+' 273.15)) ) '/' Vtotal ; 1296 : 'MESS' ' Masse (Kg) ' M0D '_______' MMD '________' 1297 : (100. '*'(M0D '-' MMD) '/'M0D) ; 1298 : 'MESS' ' Energie (J/m3)' E0D '_______' EMD '________' 1299 : (100. '*' (E0D '-' EMD) '/' E0D) ; 1300 : 'SI' (TBT.'VAPEUR' 'ET' ('NON' ('EGA' MV0D 0.0))) ; 1301 : 'MESS' ' M Vapeur(Kg) ' MV0D '_______' MVMD '_______' 1302 : (100. '*' (MV0D '-' MVMD) '/' MV0D) ; 1303 : 'FINSI' ; 1304 : 'SI' (TBT.'THE' 'ET' ('NON' ('EGA' MHE0D 0.0))) ; 1305 : 'MESS' ' M Helium (Kg) ' MHE0D '_______' MHEMD '_______' 1306 : (100. '*' (MHE0D '-' MHEMD) '/' MHE0D) ; 1307 : 'FINSI' ; 1308 : 'SI' (TBT.'TH2' 'ET' ('NON' ('EGA' MH20D 0.0))) ; 1309 : 'MESS' ' M Hydrogene (Kg) ' MH20D '_______' MH2MD '_______' 1310 : (100. '*' (MH20D '-' MH2MD) '/' MH20D) ; 1311 : 'FINSI' ; 1312 : 'SI' (TBT.'TO2' 'ET' ('NON' ('EGA' MO20D 0.0))) ; 1313 : 'MESS' ' M Oxygene (Kg) ' MO20D '_______' MO2MD '_______' 1314 : (100. '*' (MO20D '-' MO2MD) '/' MO20D) ; 1315 : 'FINSI' ; 1316 : 'SI' (TBT.'TN2' 'ET' ('NON' ('EGA' MN20D 0.0))) ; 1317 : 'MESS' ' M Azote (Kg) ' MN20D '_______' MN2MD '_______' 1318 : (100. '*' (MN20D '-' MN2MD) '/' MN20D) ; 1319 : 'FINSI' ; 1320 : 'SI' (TBT.'TCO' 'ET' ('NON' ('EGA' MCO0D 0.0))) ; 1321 : 'MESS' ' M CO (Kg) ' MCO0D '_______' MCOMD '_______' 1322 : (100. '*' (MCO0D '-' MCOMD) '/' MCO0D) ; 1323 : 'FINSI' ; 1324 : 'SI' (TBT.'TCO2' 'ET' ('NON' ('EGA' MCO20D 0.0))) ; 1325 : 'MESS' ' M CO2 (Kg) ' MCO20D '_______' MCO2MD '_______' 1326 : (100. '*' (MCO20D '-' MCO2MD) '/' MCO20D) ; 1327 : 'FINSI' ; 1328 : *? TCO2AIR 1329 : 'SI' (TBT.'TAIR' 'ET' ('NON' ('EGA' MAIR0D 0.0))) ; 1330 : 'MESS' ' M AIR (Kg) ' MAIR0D '_______' MAIRMD '_______' 1331 : (100. '*' (MAIR0D '-' MAIRMD) '/' MAIR0D) ; 1332 : 'FINSI' ; 1333 : 'MESS' '_____________________________________________' 1334 : '________________' ; 1335 : 'SI' TBT.'THERMP' ; 1336 : Emurn = 'EXTR' (tic.'Emur') ndl ; 1337 : 'MESS' ' Energie Stockée Mur (J) ' Emurn ; 1338 : 'MESS' '_____________________________________________' 1339 : '________________' ; 1340 : 'FINSI' ; 1341 : 1342 : 'MESS' ' *** Divers **** ' ; 1343 : 1344 : 'MESS' ' Dimension espace: ' idim ; 1345 : 'MESS' ' Nombre d elements/noeuds ' ('NBEL' vtf) 1346 : ('NBNO' vtf) ; 1347 : 1348 : 'MESS' '***********************************************' ; 1349 : 1350 : 1351 : Econv = 0.0 ; 1352 : Qc = 0.0 ; 1353 : Econd = 0.0 ; 1354 : Hcond = 0.0 ; 1355 : Rgp = tic.'Rgp' ; 1356 : 1357 : dRhov = 0. ; 1358 : dRhohe = 0. ; 1359 : dRhoh2 = 0. ; 1360 : dRhoo2 = 0. ; 1361 : dRhon2 = 0. ; 1362 : dRhoco = 0. ; 1363 : dRhoco2= 0. ; 1364 : 1365 : *---------------------------------------------------------- 1366 : *----- BOUCLE PRINCIPALE DE RESOLUTION SUR LE NB ---------- 1367 : *----- TOTAL DE PAS DE TEMPS A REALISER ------------------- 1368 : *---------------------------------------------------------- 1369 : 1370 : Si (NON (EXIST tic 'DT')) ; 1371 : *?DT = rxt.'DT0' ; 1372 : tic.'DT' = rxt.'DT0' ; 1373 : Finsi ; 1374 : 1375 : *>1 'REPETER' BCLTPS nbit 1376 : 'REPETER' BCLTPS nbit ; 1377 : 1378 : * Gestion de la frequence de preconditionnement 1379 : 'SI' ('OU' ('EGA' &BCLTPS 1) 1380 : ('MULT' &BCLTPS RXT.'FRPREC')) ; 1381 : calprec = VRAI ; 1382 : 'SINON' ; 1383 : calprec = FAUX ; 1384 : 'FINSI' ; 1385 : 1386 : 'SI' ('EGA' nbit 0) ; 1387 : 'QUITTER' BCLTPS ; 1388 : 'FINSI' ; 1389 : 1390 : nupadt = tic.'NUPADT' ; 1391 : nupadt = nupadt '+' 1 ; 1392 : tic.'NUPADT' = nupadt ; 1393 : *>1 ndl = 'DIME' (tic.'LTPS') 1394 : ndl = 'DIME' (tic.'LTPS') ; 1395 : 1396 : Tps='EXTR' (tic.'LTPS') ndl ; 1397 : 1398 : Tps = tic.'Tps' '+' tic.'DT' ; 1399 : tic.'Tps' = Tps ; 1400 : tic.'LTPS' = tic.'LTPS' 'ET' ('PROG' Tps ) ; 1401 : 1402 : 'SI' ( TBT.'IMPR' >EG 1) ; 1403 : 'MESS' '*************************************************************' 1404 : '***********************'; 1405 : 'MESS' '*************************************************************' 1406 : '***********************'; 1407 : 'MESS' '==> Pas de temps ' nupadt 1408 : ' Temps = ' tps 'Secondes Pas de temps : DT= ' tic.'DT' ; 1409 : 'FINSI' ; 1410 : 1411 : 1412 : **** Initialisations en début de bas de temps **************** 1413 : 1414 : **** Evolution de la source en fct du temps ****************** 1415 : ******** Scenario transitoire ******************************** 1416 : 1417 : Ksi=1. ; Ksi2=1. ; Ksi3=1. ; 1418 : Qlj=0. ; Qlj2=0. ; Qlj3=0. ; 1419 : 1420 : Si (TBT.'Tbreche') ; 1421 : qeau qlj qair qhe qh2 qo2 qn2 qco qco2 = bruche 'KAS1' 1 rxt; 1422 : Finsi; 1423 : 1424 : Si (TBT.'Tbreche2') ; 1425 : qeau2 qlj2 qair2 qhe2 qh22 qo22 qn22 q2co q2co2 = bruche 'KAS1' 2 rxt; 1426 : Finsi; 1427 : 1428 : Si (TBT.'Tbreche3') ; 1429 : qeau3 qlj3 qair3 qhe3 qh23 qo23 qn23 q3co q3co2 = bruche 'KAS1' 3 rxt; 1430 : Finsi; 1431 : 1432 : 1433 : *-------------------------------------------------------------- 1434 : * Interpolation des conditions aux limites en température 1435 : * imposée sur une ou plusieurs parois et chargement coef d'échange 1436 : * TTIMP1 -> tic.'TBP1' etc 1437 : 1438 : 'SI' TBT.'TTIMP1' ; 1439 : tic.'TBP1' = 'IPOL' tps (rxt.'TIMP1'.'t') (rxt.'TIMP1'.'TIMP'); 1440 : Si (NON ((TBT.'VAPEUR') et (exist tic 'KH1'))) ; 1441 : tic.'KH1' = (rxt.'TIMP1'.'ECHAN') ; 1442 : Finsi ; 1443 : tic.'Ltbp1'= tic.'Ltbp1' et (prog (tic.'TBP1')) ; 1444 : 'FINSI' ; 1445 : 1446 : 'SI' TBT.'TTIMP2' ; 1447 : tic.'TBP2' = 'IPOL' tps (rxt.'TIMP2'.'t') (rxt.'TIMP2'.'TIMP'); 1448 : Si (NON ((TBT.'VAPEUR') et (exist tic 'KH2'))) ; 1449 : tic.'KH2' = (rxt.'TIMP2'.'ECHAN') ; 1450 : Finsi ; 1451 : tic.'Ltbp2'= tic.'Ltbp2' et (prog (tic.'TBP2')) ; 1452 : 'FINSI' ; 1453 : 1454 : 'SI' TBT.'TTIMP3' ; 1455 : tic.'TBP3' = 'IPOL' tps (rxt.'TIMP3'.'t') (rxt.'TIMP3'.'TIMP'); 1456 : Si (NON ((TBT.'VAPEUR') et (exist tic 'KH3'))) ; 1457 : tic.'KH3' = (rxt.'TIMP3'.'ECHAN') ; 1458 : Finsi ; 1459 : tic.'Ltbp3'= tic.'Ltbp3' et (prog (tic.'TBP3')) ; 1460 : 'FINSI' ; 1461 : 1462 : 'SI' TBT.'TECHANP' ; 1463 : 'SI' ('NON' ('EXIST' TIC 'TBP0')) ; 1464 : tic.'TBP0' = 'KCHT' RXT.'GEO'.'$mtp0' 'SCAL' 'SOMMET' 1465 : rxt.'ECHANP'.'TMUR' ; 1466 : 'FINSI' ; 1467 : Si (NON ((TBT.'VAPEUR') et (exist tic 'KH0'))) ; 1468 : tic.'KH0' = 'KCHT' RXT.'GEO'.'$mtp0' 'SCAL' 'CENTRE' 1469 : rxt.'ECHANP'.'ECHAN' ; 1470 : Finsi ; 1471 : 'FINSI' ; 1472 : 1473 : *-------------------------------------------------------------- 1474 : * Traitement de l'aspersion 1475 : 1476 : 'SI' TBT.'ASPER ' ; 1477 : t = rxt.'scenasp'.'t' ; 1478 : vzinj = 'IPOL' tps t (rxt.'scenasp'.'vzinj') ; 1479 : xdinj = 'IPOL' tps t (rxt.'scenasp'.'xdinj') ; 1480 : tdinj = 'IPOL' tps t (rxt.'scenasp'.'tdinj') ; 1481 : ddinj = 'IPOL' tps t (rxt.'scenasp'.'ddinj') ; 1482 : $aspinj = GEO.'$aspinj' ; 1483 : $toitf = GEO.'$toitf' ; 1484 : * vnj = 'KCHT' $aspinj 'VECT' 'SOMMET' 'COMP' 1485 : * ('MOTS' 1VN 2VN 3VN) (0.0 0.0 vzinj) ; 1486 : 1487 : 'SI' (DIM3D) ; 1488 : vnj = 'KCHT' $aspinj 'VECT' 'SOMMET' 'COMP' 1489 : ('MOTS' 1VN 2VN 3VN) (0.0 0.0 vzinj) ; 1490 : 'SINON' ; 1491 : *-------------------------------------------------------------- 1492 : * modif Tojo 12/07/2007 * 1493 : * 1494 : * mise en place de vecteurs inclinés conformément 1495 : * à l'angle d'ouverture du jet 1496 : * ici l'inclinaison est fixée à 30° 1497 : * une prochaine évolution pourrait consister à mettre 1498 : * l'angle d'ouverture comme donnée d'entrée 1499 : * dans la table scenasp 1500 : * 1501 : *-------------------------------------------------------------- 1502 : * 1503 : cc1 = coor 1 rxt.'aspinj' ; 1504 : z0 = (sin(300. '*' cc1)) '*' ((cos(300. '*' cc1)) '**' (-1.0)); 1505 : z0 = z0 '*' (vzinj '*' -1.) ; 1506 : z0 = chan attribut z0 nature diffus ; 1507 : z1 = z0 '*' 0. ; 1508 : ux0 = nomc '1VN' z0 ; 1509 : z1 = z1 '+' vzinj ; 1510 : z1 = chan attribut z1 nature diffus ; 1511 : uy0 = nomc '2VN' z1 ; 1512 : vnj1 = ux0 et uy0 ; 1513 : vnj = 'KCHT' $aspinj 'VECT' 'SOMMET' 'COMP' 1514 : ('MOTS' '1VN' '2VN') vnj1 ; 1515 : * 1516 : *-------------------------------------------------------------- 1517 : 'FINSI' ; 1518 : *--------------------fin modif Tojo Juillet 2007--------------- 1519 : 1520 : xdj = 'KCHT' $aspinj 'SCAL' 'SOMMET' 'COMP' 1521 : ('MOT' 'XD') xdinj ; 1522 : tdj = 'KCHT' $aspinj 'SCAL' 'SOMMET' 'COMP' 1523 : ('MOT' 'TD') tdinj ; 1524 : ddj = 'KCHT' $toitf 'SCAL' 'SOMMET' 'COMP' 1525 : ('MOT' 'DD') ddinj ; 1526 : 'FINSI' ; 1527 : 1528 : *-------------------------------------------------------------- 1529 : *>1********************************************************** 1530 : *>1* Grandeurs moyennes au debut du pas de temps N ********** 1531 : *>1* i.e. au temps N-1 ************************************** 1532 : 1533 : Ptnm = 'EXTR' (tic.'PT') ndl ; 1534 : dPdtnm = 'EXTR' (tic.'dPdt') ndl ; 1535 : Rgpmnm = 'EXTR' (tic.'Rgpm') ndl ; 1536 : Cvmnm = 'EXTR' (tic.'Cvm') ndl ; 1537 : Gamnm = 'EXTR' (tic.'Gamm') ndl ; 1538 : Cpmnm = 'EXTR' (tic.'Cpm') ndl ; 1539 : Rhomnm = 'EXTR' (tic.'Rhom') ndl ; 1540 : 1541 : Rhomvn = 0. ; 1542 : Rhovgn = 0. ; 1543 : Rbrom = 0. ; 1544 : Rmhen = 0. ; 1545 : Rmh2n = 0. ; 1546 : Rmn2n = 0. ; 1547 : Rmo2n = 0. ; 1548 : Rmcon = 0. ; 1549 : Rmco2n = 0. ; 1550 : Rmairn = 0. ; 1551 : 1552 : 'SI' TBT.'VAPEUR' ; 1553 : *>1 Rhomvnm,Rhovgnm,Rliqm = EXTR tic.'Rhomv','Rhomvg','Rliqm' ndl 1554 : Rhomvnm = 'EXTR' (tic.'Rhomv') ndl ; 1555 : Rhovgnm = 'EXTR' (tic.'Rhomvg') ndl ; 1556 : 'SI' TBT.'CONDMAS' ; 1557 : Rliqm = 'EXTR' (tic.'Rliqm') ndl ; 1558 : Rbrom = 'EXTR' (tic.'Rbrom') ndl ; 1559 : 'FINSI' ; 1560 : 'FINSI' ; 1561 : 'SI' TBT.'THE' ; 1562 : Rmhenm = 'EXTR' (tic.'Rhomhe') ndl ; 1563 : 'FINSI' ; 1564 : 'SI' TBT.'TH2' ; 1565 : Rmh2nm = 'EXTR' (tic.'Rhomh2') ndl ; 1566 : 'FINSI' ; 1567 : 'SI' TBT.'TO2' ; 1568 : Rmo2nm = 'EXTR' (tic.'Rhomo2') ndl ; 1569 : 'FINSI' ; 1570 : 'SI' TBT.'TN2' ; 1571 : Rmn2nm = 'EXTR' (tic.'Rhomn2') ndl ; 1572 : 'FINSI' ; 1573 : 'SI' TBT.'TCO' ; 1574 : Rmconm = 'EXTR' (tic.'Rhomco') ndl ; 1575 : 'FINSI' ; 1576 : 'SI' TBT.'TCO2' ; 1577 : Rmco2nm = 'EXTR' (tic.'Rhomco2') ndl ; 1578 : 'FINSI' ; 1579 : 1580 : Rmairnm = 'EXTR' (tic.'Rhomair') ndl ; 1581 : Remnm = 'EXTR' (tic.'Remn') ndl ; 1582 : 'SI' TBT.'THERMP' ; 1583 : Emurnm = 'EXTR' (tic.'Emur') ndl ; 1584 : 'FINSI' ; 1585 : 'SI' TBT.'ASPER ' ; 1586 : Qasp = extr (tic.'Qaspe') ndl ; 1587 : Easp = extr (tic.'Easpe') ndl ; 1588 : Hasp = extr (tic.'Haspe') ndl ; 1589 : 'SINON' ; 1590 : Qasp = 0.0 ; 1591 : Easp = 0.0 ; 1592 : Hasp = 0.0 ; 1593 : 'FINSI' ; 1594 : 1595 : Qcnm = 'EXTR' (tic.'Qc') ndl ; 1596 : Qcm = 0. ; 1597 : * Qcm = 'EXTR' (tic.'Qcm') ndl ; 1598 : * Pour faire marcher la version V0 1599 : Econdnm = 'EXTR' (tic.'Econd') ndl ; 1600 : Econvnm = 'EXTR' (tic.'Econv') ndl ; 1601 : Hcondnm = 'EXTR' (tic.'Hcond') ndl ; 1602 : MinjT = 'EXTR' (tic.'Minj') ndl ; 1603 : McondT = 'EXTR' (tic.'Mcond') ndl ; 1604 : 1605 : 'SI' (TBT.'VAPEUR') ; 1606 : *? 'SI' (TBT.'CONDMAS') ; 1607 : 'SI' (TBT.'VERSION' > 0) ; 1608 : 'SI' (EXIST TIC 'Mcondm') ; 1609 : Mcondm = 'EXTR' (tic.'Mcondm') ndl ; 1610 : 'SINON' ; 1611 : tic.'Mcondm' = 'PROG' 0.0 ; 1612 : Mcondm = 0.0 ; 1613 : 'FINSI' ; 1614 : 'SI' (NON (EXIST TIC 'Qcm')) ; 1615 : tic.'Qcm' = 'PROG' 0.0 ; 1616 : 'FINSI' ; 1617 : 'FINSI' ; 1618 : 1619 : 'SI' (TBT.'TPAROIF') ; 1620 : 'SI' (EXIST TIC 'Mcondw') ; 1621 : McondTw = 'EXTR' (tic.'Mcondw') ndl ; 1622 : 'SINON' ; 1623 : tic.'Mcondw' = 'PROG' 0.0 ; 1624 : McondTw = 0.0 ; 1625 : 'FINSI' ; 1626 : 'SI' (NON (EXIST TIC 'Qcw')) ; 1627 : tic.'Qcw' = 'PROG' 0.0 ; 1628 : 'FINSI' ; 1629 : 'FINSI' ; 1630 : 'SI' (TBT.'TTIMP1') ; 1631 : 'SI' (EXIST TIC 'Mcond1') ; 1632 : McondT1 = 'EXTR' (tic.'Mcond1') ndl ; 1633 : 'SINON' ; 1634 : tic.'Mcond1' = 'PROG' 0.0 ; 1635 : McondT1 = 0.0 ; 1636 : 'FINSI' ; 1637 : 'SI' (NON (EXIST TIC 'Qc1')) ; 1638 : tic.'Qc1' = 'PROG' 0.0 ; 1639 : 'FINSI' ; 1640 : 'FINSI' ; 1641 : 'SI' (TBT.'TTIMP2') ; 1642 : 'SI' (EXIST TIC 'Mcond2') ; 1643 : McondT2 = 'EXTR' (tic.'Mcond2') ndl ; 1644 : 'SINON' ; 1645 : tic.'Mcond2' = 'PROG' 0.0 ; 1646 : McondT2 = 0.0 ; 1647 : 'FINSI' ; 1648 : 'SI' (NON (EXIST TIC 'Qc2')) ; 1649 : tic.'Qc2' = 'PROG' 0.0 ; 1650 : 'FINSI' ; 1651 : 'FINSI' ; 1652 : 'SI' (TBT.'TTIMP3') ; 1653 : 'SI' (EXIST TIC 'Mcond3') ; 1654 : McondT3 = 'EXTR' (tic.'Mcond3') ndl ; 1655 : 'SINON' ; 1656 : tic.'Mcond3' = 'PROG' 0.0 ; 1657 : McondT3 = 0.0 ; 1658 : 'FINSI' ; 1659 : 'SI' (NON (EXIST TIC 'Qc3')) ; 1660 : tic.'Qc3' = 'PROG' 0.0 ; 1661 : 'FINSI' ; 1662 : 'FINSI' ; 1663 : 'SI' (TBT.'TECHANP') ; 1664 : 'SI' (EXIST TIC 'Mcond0') ; 1665 : McondT0 = 'EXTR' (tic.'Mcond0') ndl ; 1666 : 'SINON' ; 1667 : tic.'Mcond0' = 'PROG' 0.0 ; 1668 : McondT0 = 0.0 ; 1669 : 'FINSI' ; 1670 : 'SI' (NON (EXIST TIC 'Qc0')) ; 1671 : tic.'Qc0' = 'PROG' 0.0 ; 1672 : 'FINSI' ; 1673 : 'FINSI' ; 1674 : 'FINSI' ; 1675 : MrestT = 'EXTR' (tic.'Mrest') ndl ; 1676 : 1677 : *------------------------------------------------------------- 1678 : * Calcul KH/Ro Cp pour l'équation en température TF 1679 : 1680 : 'SI' TBT.'THERMP' ; 1681 : tic.'KHEW'=tic.'KHW' / (Rhomnm '*' Cpmnm) ; 1682 : tic.'LHEW' = -1.0 * tic.'KHEW' ; 1683 : 'FINSI' ; 1684 : 1685 : 'SI' TBT.'TTIMP1' ; 1686 : tic.'KHE1' = (tic.'KH1') '/' (Rhomnm '*' Cpmnm) ; 1687 : 'FINSI' ; 1688 : 1689 : 'SI' TBT.'TTIMP2' ; 1690 : tic.'KHE2' = (tic.'KH2') '/' (Rhomnm '*' Cpmnm) ; 1691 : 'FINSI' ; 1692 : 1693 : 'SI' TBT.'TTIMP3' ; 1694 : tic.'KHE3' = (tic.'KH3') '/' (Rhomnm '*' Cpmnm) ; 1695 : 'FINSI' ; 1696 : 1697 : 'SI' TBT.'TECHANP' ; 1698 : tic.'KHE0' = (tic.'KH0') '/' (Rhomnm '*' Cpmnm) ; 1699 : 'FINSI' ; 1700 : 1701 : *============================================================= 1702 : *BRECHE 1703 : *============================================================= 1704 : Qj = 0. ; 1705 : Roj = 1.e-10 ; 1706 : src = 0. ; 1707 : guj = 0. ; 1708 : Ej=0. ; 1709 : Hj=0. ; 1710 : Tj=0. ; 1711 : 1712 : Si (TBT.'Tbreche') ; 1713 : Rvapj Rairj Rhej Rh2j Ro2j Rn2j Rco2j Rcoj Rgj Hj Ej 1714 : Tj Qj Roj gj guj unvf0 src = bruche 'KAS2' 1 rxt ; 1715 : asrc = Roj '*' ('ABS' src) ; 1716 : Finsi; 1717 : 1718 : *---------------------------------------------------------- 1719 : * STOCKAGE des valeurs de l'injection pour le temps courant 1720 : * 1ere Breche 1721 : * La brèche peut avoir eu une existance évanescente (Tbreche=FAUX) 1722 : * Si tic.'Qj' existe on complete avec des valeurs nulles 1723 : * tic.'Qj',tic.'guj',tic.'Hj',tic.'Ej' 1724 : * 1725 : Si(Exist tic 'Qj') ; 1726 : tic.'guj'= tic.'guj' 'ET' ('PROG' guj) ; 1727 : asrc = Roj '*' ('ABS' src) ; 1728 : tic.'Qj' = tic.'Qj' 'ET' ('PROG' asrc) ; 1729 : 1730 : tic.'Hj'= tic.'Hj' 'ET' ('PROG' Hj) ; 1731 : tic.'Ej'= tic.'Ej' 'ET' ('PROG' Ej) ; 1732 : Finsi ; 1733 : 1734 : *============================================================= 1735 : *BRECHE 2 1736 : *============================================================= 1737 : Qj2 = 0. ; 1738 : Roj2 = 1.e-10 ; 1739 : src2 = 0. ; 1740 : guj2 = 0. ; 1741 : Ej2=0. ; 1742 : Hj2=0. ; 1743 : Tj2=0. ; 1744 : 1745 : Si (TBT.'Tbreche2') ; 1746 : Rvapj2 Rairj2 Rhej2 Rh2j2 Ro2j2 Rn2j2 Rco2j2 Rcoj2 Rgj2 Hj2 Ej2 1747 : Tj2 Qj2 Roj2 gj2 guj2 unvf0 src2 = bruche 'KAS2' 2 rxt ; 1748 : 1749 : asrc2 = Roj2 '*' ('ABS' src2) ; 1750 : Finsi; 1751 : 1752 : *---------------------------------------------------------- 1753 : * STOCKAGE des valeurs de l'injection pour le temps courant 1754 : * 2eme Breche 1755 : * La brèche peut avoir eu une existance évanescente (Tbreche2=FAUX) 1756 : * Si tic.'Qj2' existe on complete avec des valeurs nulles 1757 : * tic.'Qj2',tic.'guj2',tic.'Hj2',tic.'Ej2' 1758 : * 1759 : Si(Exist tic 'Qj2') ; 1760 : tic.'guj2'= tic.'guj2' 'ET' ('PROG' guj2) ; 1761 : asrc2 = Roj2 '*' ('ABS' src2) ; 1762 : tic.'Qj2' = tic.'Qj2' 'ET' ('PROG' asrc2) ; 1763 : tic.'Hj2'= tic.'Hj2' 'ET' ('PROG' Hj2) ; 1764 : tic.'Ej2'= tic.'Ej2' 'ET' ('PROG' Ej2) ; 1765 : Finsi ; 1766 : *============================================================= 1767 : *BRECHE 3 1768 : *============================================================= 1769 : Qj3 = 0. ; 1770 : Roj3 = 1.e-10 ; 1771 : src3 = 0.; 1772 : guj3 = 0. ; 1773 : Ej3=0. ; 1774 : Hj3=0. ; 1775 : Tj3=0. ; 1776 : 1777 : Si (TBT.'Tbreche3') ; 1778 : Rvapj3 Rairj3 Rhej3 Rh2j3 Ro2j3 Rn2j3 Rco2j3 Rcoj3 Rgj3 Hj3 Ej3 1779 : Tj3 Qj3 Roj3 gj3 guj3 unvf0 src3 = bruche 'KAS2' 3 rxt ; 1780 : 1781 : asrc3 = Roj3 '*' ('ABS' src3) ; 1782 : Finsi; 1783 : 1784 : *---------------------------------------------------------- 1785 : * STOCKAGE des valeurs de l'injection pour le temps courant 1786 : * 3ème Breche 1787 : * La brèche peut avoir eu une existance évanescente (Tbreche3=FAUX) 1788 : * Si tic.'Qj3' existe on complète avec des valeurs nulles 1789 : * tic.'Qj3',tic.'guj3',tic.'Hj3',tic.'Ej3' 1790 : * 1791 : Si(Exist tic 'Qj3') ; 1792 : tic.'guj3'= tic.'guj3' 'ET' ('PROG' guj3) ; 1793 : asrc3 = Roj3 '*' ('ABS' src3) ; 1794 : tic.'Qj3' = tic.'Qj3' 'ET' ('PROG' asrc3) ; 1795 : tic.'Hj3'= tic.'Hj3' 'ET' ('PROG' Hj3) ; 1796 : tic.'Ej3'= tic.'Ej3' 'ET' ('PROG' Ej3) ; 1797 : Finsi ; 1798 : *============================================================= 1799 : 1800 : 1801 : 'SI' TBT.'TRECOMB' ; 1802 : 1803 : MDRECOMB RXT ; 1804 : 1805 : 'FINSI' ; 1806 : 1807 : *------------------------------------------------------------ 1808 : **** Avancement des grandeurs moyennes 0D ******************* 1809 : **** compte tenu des injections ******************* 1810 : * Pour le calcul de RHO, P 1811 : * 1812 : 1813 : 1814 : * Calcul de l energie moyenne au pas de temps N 1815 : Htot = Hcondnm '-' Hasp ; 1816 : Etot = Econvnm '-' Easp ; 1817 : dr = 0. '-' Etot '-' Htot ; 1818 : 1819 : *mess ' Etot ' Etot ' Htot ' Htot 'dr=' dr; 1820 : *mess 'Econv pour calcul de dr ' Econv; 1821 : 1822 : * Calcul de la densité moyenne au pas de temps N 1823 : Qtot = Qasp '-' Qcnm ; 1824 : * Source de vapeur liée à la première brèche : 1825 : Si (TBT.'Tbreche') ; 1826 : Qtot = Qtot '+' asrc ; 1827 : dr = dr + (Qj '*' Hj) ; 1828 : Finsi ; 1829 : * Source de vapeur liée à la deuxième brèche : 1830 : Si (TBT.'Tbreche2') ; 1831 : Qtot = Qtot '+' asrc2 ; 1832 : dr = dr + (Qj2 '*' Hj2) ; 1833 : Finsi ; 1834 : * Source de vapeur liée à la troisième brèche : 1835 : Si (TBT.'Tbreche3') ; 1836 : Qtot = Qtot '+' asrc3 ; 1837 : dr = dr + (Qj3 '*' Hj3) ; 1838 : Finsi ; 1839 : 1840 : * Prise en compte des données du ou des recombineurs eventuels 1841 : 'SI' TBT.'TRECOMB' ; 1842 : 1843 : NBREC = 'DIME' rxt.'RECOMB' ; 1844 : 'REPETER' BCLREC NBREC ; 1845 : 1846 : DIRNAME = 'CHAIN' 'DREC' 'ent' &BCLREC ; 1847 : direntr = RXT.'GEO'. DIRNAME ; 1848 : DIRNAME = 'CHAIN' 'DREC' 'sor' &BCLREC ; 1849 : dirsort = RXT.'GEO'. DIRNAME ; 1850 : 1851 : $NAME = 'CHAIN' '$REC' 'ent' &BCLREC ; 1852 : $ENTREE = RXT.'GEO'. $NAME ; 1853 : $NAME = 'CHAIN' '$REC' 'sor' &BCLREC ; 1854 : $SORTIE = RXT.'GEO'. $NAME ; 1855 : 1856 : 'SI' TBT.'FPAROI' ; 1857 : NAME = 'CHAIN' 'REC' 'ent' &BCLREC ; 1858 : ENTREE = RXT.'GEO'.NAME ; 1859 : NAME = 'CHAIN' 'REC' 'sor' &BCLREC ; 1860 : SORTIE = RXT.'GEO'.NAME ; 1861 : 'SINON' ; 1862 : NAME = 'CHAIN' 'REC' 'entI' &BCLREC ; 1863 : ENTREE = RXT.'GEO'.NAME ; 1864 : NAME = 'CHAIN' 'REC' 'sorI' &BCLREC ; 1865 : SORTIE = RXT.'GEO'.NAME ; 1866 : 'FINSI' ; 1867 : 1868 : LMOT1 = 'MOTS' 1UN 2UN ; 1869 : LMOT2 = 'MOTS' UX UY ; 1870 : 1871 : 'SI' DIM3D; 1872 : LMOT1 = LMOT1 'ET' ('MOTS' 3UN ) ; 1873 : LMOT2 = LMOT2 'ET' ('MOTS' UZ ) ; 1874 : 'FINSI' ; 1875 : NRECT = 'DIME' RXT.'TIC'.'RECOMB'.&BCLREC.'Uin' ; 1876 : Uin1 = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'Uin' NRECT ; 1877 : Uout1 = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'Uout' NRECT ; 1878 : DEB1 = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'DEB' NRECT ; 1879 : Hin1 = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'Hin' NRECT ; 1880 : Hout1 = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'Hout' NRECT ; 1881 : toto1 = 'KCHT' $ENTREE 'VECT' 'SOMMET' 'COMP' 1882 : LMOT1 (direntr * Uin1); 1883 : toto1 = 'REDU' toto1 ENTREE ; 1884 : toto1 = 'KCHT' $ENTREE 'VECT' 'SOMMET' 'COMP' 1885 : LMOT1 toto1; 1886 : toto2 = 'NOMC' LMOT1 LMOT2 toto1 ; 1887 : toto3 = 'KCHT' RXT.'GEO'.'$menvf' 'VECT' 'SOMMET' toto2 ; 1888 : toto3 = 'DBIT' toto3 RXT.'GEO'.'$menvf' ; 1889 : 1890 : tutu1 = 'KCHT' $SORTIE 'VECT' 'SOMMET' 'COMP' 1891 : LMOT1 (dirsort * Uout1); 1892 : tutu1 = 'REDU' tutu1 SORTIE ; 1893 : tutu1 = 'KCHT' $SORTIE 'VECT' 'SOMMET' 'COMP' 1894 : LMOT1 tutu1; 1895 : tutu2 = 'NOMC' LMOT1 LMOT2 tutu1 ; 1896 : tutu3 = 'KCHT' RXT.'GEO'.'$menvf' 'VECT' 'SOMMET' tutu2 ; 1897 : tutu3 = 'DBIT' tutu3 RXT.'GEO'.'$menvf' ; 1898 : 1899 : 'SI' ('EGA' &BCLREC 1) ; 1900 : Urec0 = toto1 'ET' tutu1 ; 1901 : * Srho0 = tutu3 '-' toto3 ; 1902 : 'SINON' ; 1903 : Urec0 = Urec0 'ET' toto1 'ET' tutu1 ; 1904 : * Srho0 = Srho0 '+' tutu3 '-' toto3 ; 1905 : 'FINSI' ; 1906 : Srho0 = 0.0 ; 1907 : 1908 : ** Dans le recombineur, on suppose que le debit massique entrant 1909 : ** est egale au debit massique sortant 1910 : ** ==> auncune correction n'est a effectué sur la densité totale ! 1911 : * Qtot = Qtot '+' 0.0 ; 1912 : 1913 : dr = dr '+' (DEB1 '*' (Hout1 '-' Hin1)) ; 1914 : 1915 : 'FIN' BCLREC ; 1916 : 1917 : 'FINSI' ; 1918 : * Modif. FP 19/05/00 1919 : * Retranchement de la masse et de l'énergie perdues par condensation 1920 : * en masse : 1921 : 1922 : *µµµµµµµµµµµµµ calcul Rhomn et Remn µµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµ 1923 : Si TBT.'Tsortie'; 1924 : 1925 : Tf=tic.'TF' ; 1926 : Rho=tic.'RHO' ; 1927 : Cvm = 'EXTR' (tic.'Cvm') ndl ; 1928 : Rhomn = Rhom ; 1929 : Remn = (Somt(Diag * Rho * cvm * (Tf+273.15))) '/' VTotal ; 1930 : tic.'Rhom' = tic.'Rhom' 'ET' ('PROG' Rhomn) ; 1931 : tic.'Remn' = tic.'Remn' 'ET' ('PROG' Remn) ; 1932 : 1933 : Sinon; 1934 : 1935 : * Détermination de la masse volumique et de l'énergie interne moyenne 1936 : 1937 : Rhomn = MAXI (PROG (Rhomnm '+' (tic.'DT' '*' Qtot '/' VTotal )) 0.0); 1938 : Remn = MAXI (PROG (Remnm '+' (tic.'DT' '*' dr '/' VTotal ) ) 0.0); 1939 : 1940 : tic.'Rhom' = tic.'Rhom' 'ET' ('PROG' Rhomn) ; 1941 : tic.'Remn' = tic.'Remn' 'ET' ('PROG' Remn) ; 1942 : 1943 : Finsi ; 1944 : *µµµµµµµµµµµµµ calcul Rhomn et Remn µµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµ 1945 : 1946 : * Calcul de l'energie dans le mur au pas de temps N 1947 : * Il manque la partie convective !!!! 1948 : 'SI' TBT.'THERMP' ; 1949 : Hmur = Qcnm '*' TBT.'Lv' ; 1950 : Emurn = Emurnm '+' ( tic.'DT' '*' Hmur) ; 1951 : tic.'Emur' = tic.'Emur' 'ET' ('PROG' Emurn) ; 1952 : 'FINSI' ; 1953 : 1954 : * Calcul de la densité moyenne de vapeur au pas de temps N 1955 : 'SI' TBT.'VAPEUR' ; 1956 : Qo=0. ; Qlo=0.; 1957 : Si (TBT.'Tbreche') ; Qo=Qeau ; Qlo=Qlj ; Finsi ; 1958 : Si (TBT.'Tbreche2') ; Qo=Qo + Qeau2; Qlo=Qlo + Qlj2 ; Finsi ; 1959 : Si (TBT.'Tbreche3') ; Qo=Qo + Qeau3; Qlo=Qlo + Qlj3 ; Finsi ; 1960 : 1961 : 'SI' TBT.'TRECOMB' ; 1962 : QIN_H2O = 'EXTR' RXT.'TIC'.'QIN_H2O' ('DIME' RXT.'TIC'.'QIN_H2O') ; 1963 : QOUT_H2O = 'EXTR' RXT.'TIC'.'QOUT_H2O' ('DIME' RXT.'TIC'.'QOUT_H2O') ; 1964 : Qo=Qo + QOUT_H2O - QIN_H2O ; 1965 : 'FINSI' ; 1966 : *>2 Calcul Rhovgn Rhomvn -> tic.'Rhomvg' et tic.'Rhomv' 1967 : *>2 Rhovgn= MAXI (PROG (Rhovgnm '+' 1968 : *>2 (tic.'DT' '*' (Qo '-' Qcnm '-' Qcm '+' Qasp)'/' VTotal)) 0.); 1969 : *>2 Rhomvn= Rhovgn - Rbrom ; 1970 : 1971 : Rhovgn= MAXI (PROG (Rhovgnm '+' 1972 : (tic.'DT' '*' (Qo '-' Qcnm '+' Qasp - Qcm) '/' VTotal )) 0.0) ; 1973 : Rhomvn= Rhovgn - Rbrom ; 1974 : 1975 : tic.'Rhomv' = tic.'Rhomv' 'ET' ('PROG' Rhomvn) ; 1976 : tic.'Rhomvg'= tic.'Rhomvg' 'ET' ('PROG' Rhovgn) ; 1977 : 1978 : Mlq = (tic.'DT' '*' (Qcnm + Qlo)) ; 1979 : 1980 : 'FINSI' ; 1981 : 1982 : 'SI' TBT.'THE' ; 1983 : Qo=0. ; 1984 : Si (TBT.'Tbreche') ; Qo=Qhe ; Finsi ; 1985 : Si (TBT.'Tbreche2') ; Qo=Qo + Qhe2; Finsi ; 1986 : Si (TBT.'Tbreche3') ; Qo=Qo + Qhe3; Finsi ; 1987 : Rmhen= MAXI (PROG (Rmhenm '+' (tic.'DT' '*' Qo '/' VTotal )) 0.0) ; 1988 : tic.'Rhomhe' = tic.'Rhomhe' 'ET' ('PROG' Rmhen) ; 1989 : 'FINSI' ; 1990 : 1991 : 'SI' TBT.'TH2' ; 1992 : Qo=0. ; 1993 : Si (TBT.'Tbreche') ; Qo=Qh2 ; Finsi ; 1994 : Si (TBT.'Tbreche2') ; Qo=Qo + Qh22; Finsi ; 1995 : Si (TBT.'Tbreche3') ; Qo=Qo + Qh23; Finsi ; 1996 : 1997 : 'SI' TBT.'TRECOMB' ; 1998 : QIN_H2 = 'EXTR' RXT.'TIC'.'QIN_H2' ('DIME' RXT.'TIC'.'QIN_H2' ) ; 1999 : QOUT_H2 = 'EXTR' RXT.'TIC'.'QOUT_H2' ('DIME' RXT.'TIC'.'QOUT_H2') ; 2000 : Qo=Qo + QOUT_H2 - QIN_H2 ; 2001 : 'FINSI' ; 2002 : 2003 : Rmh2n= MAXI (PROG (Rmh2nm '+' (tic.'DT' '*' Qo '/' VTotal )) 0.0) ; 2004 : tic.'Rhomh2' = tic.'Rhomh2' 'ET' ('PROG' Rmh2n) ; 2005 : 'FINSI' ; 2006 : 2007 : 'SI' TBT.'TN2' ; 2008 : Qno= 0. ; 2009 : Si (TBT.'Tbreche') ; Qno=Qn2; Finsi ; 2010 : Si (TBT.'Tbreche2') ; Qno=Qno + Qn22; Finsi ; 2011 : Si (TBT.'Tbreche3') ; Qno=Qno + Qn23; Finsi ; 2012 : 2013 : 'SI' TBT.'TRECOMB' ; 2014 : QIN_N2 = 'EXTR' RXT.'TIC'.'QIN_N2' ('DIME' RXT.'TIC'.'QIN_N2' ) ; 2015 : QOUT_N2 = 'EXTR' RXT.'TIC'.'QOUT_N2' ('DIME' RXT.'TIC'.'QOUT_N2') ; 2016 : Qo=Qo + QOUT_N2 - QIN_N2 ; 2017 : 'FINSI' ; 2018 : 2019 : Rmn2n= MAXI (PROG (Rmn2nm '+' (tic.'DT' '*' Qno '/' VTotal )) 0.0); 2020 : tic.'Rhomn2' = tic.'Rhomn2' 'ET' ('PROG' Rmn2n) ; 2021 : 'FINSI' ; 2022 : 2023 : 'SI' TBT.'TO2' ; 2024 : Qo=0. ; 2025 : Si (TBT.'Tbreche') ; Qo=Qo2 ; Finsi ; 2026 : Si (TBT.'Tbreche2') ; Qo=Qo + Qo22; Finsi ; 2027 : Si (TBT.'Tbreche3') ; Qo=Qo + Qo23; Finsi ; 2028 : 2029 : 'SI' TBT.'TRECOMB' ; 2030 : QIN_O2 = 'EXTR' RXT.'TIC'.'QIN_O2' ('DIME' RXT.'TIC'.'QIN_O2' ) ; 2031 : QOUT_O2 = 'EXTR' RXT.'TIC'.'QOUT_O2' ('DIME' RXT.'TIC'.'QOUT_O2') ; 2032 : Qo=Qo + QOUT_O2 - QIN_O2 ; 2033 : 'FINSI' ; 2034 : 2035 : Rmo2n= MAXI (PROG (Rmo2nm '+' (tic.'DT' '*' Qo '/' VTotal )) 0.0) ; 2036 : tic.'Rhomo2' = tic.'Rhomo2' 'ET' ('PROG' Rmo2n) ; 2037 : 'FINSI' ; 2038 : 2039 : 'SI' TBT.'TCO' ; 2040 : Qo=0. ; 2041 : Si (TBT.'Tbreche') ; Qo=Qco ; Finsi ; 2042 : Si (TBT.'Tbreche2') ; Qo=Qo + Q2co; Finsi ; 2043 : Si (TBT.'Tbreche3') ; Qo=Qo + Q3co; Finsi ; 2044 : Rmcon= MAXI (PROG (Rmconm '+' (tic.'DT' '*' Qo '/' VTotal )) 0.0) ; 2045 : tic.'Rhomco' = tic.'Rhomco' 'ET' ('PROG' Rmcon) ; 2046 : 'FINSI' ; 2047 : 2048 : 'SI' TBT.'TCO2' ; 2049 : Qo=0. ; 2050 : Si (TBT.'Tbreche') ; Qo=Qco2 ; Finsi ; 2051 : Si (TBT.'Tbreche2') ; Qo=Qo + Q2co2; Finsi ; 2052 : Si (TBT.'Tbreche3') ; Qo=Qo + Q3co2; Finsi ; 2053 : Rmco2n= MAXI (PROG (Rmco2nm '+' (tic.'DT' '*' Qo '/' VTotal )) 0.0); 2054 : tic.'Rhomco2' = tic.'Rhomco2' 'ET' ('PROG' Rmco2n) ; 2055 : 'FINSI' ; 2056 : 2057 : *D'où la masse d'air totale 2058 : *Rmairn = Rhomn - Rhovgn - Rmhen - Rmh2n - Rmn2n - 2059 : Rmairn = Rhomn - Rhomvn - Rmhen - Rmh2n - Rmn2n - 2060 : Rmo2n - Rmcon - Rmco2n ; 2061 : 'SI' (Rmairn '<' 0.0) ; 2062 : 'MESS' 'Pb avec la masse d air 0D' ; 2063 : 'FINSI'; 2064 : tic.'Rhomair' = tic.'Rhomair' 'ET' ('PROG' Rmairn) ; 2065 : * 2066 : *? d_ME = tic.'DT' * (Easp + Hasp) ; 2067 : *? d_RhoE = tic.'DT' * (Easp + Hasp) / Vtotal ; 2068 : * 2069 : *µµµµµµµµ calcul de la Pression 0D au pas de temps N µµµµµµµµµµµµµµµµµµµ 2070 : * Si Cavité ouverte Pression nulle ! 2071 : Si TBT.'Tsortie' ; 2072 : Pt = 0. ; 2073 : DPDT = 0. ; 2074 : Sinon ; 2075 : * Calcul de la Pression 0D au pas de temps N 2076 : * Calcul du terme DPDT 0D au pas de temps N 2077 : * Estimation de la pression Pt au temps N 2078 : 2079 : Pt = (Gamnm '-' 1.) '*' Remn ; 2080 : * mess 'Estimation de Pt =' Pt; 2081 : DPDT = (Pt '-' Ptnm) '/' (tic.'DT') ; 2082 : Finsi ; 2083 : tic.'dPdt' = tic.'dPdt' 'ET' ('PROG' DPDT) ; 2084 : 2085 : *µµµµµµµµ calcul de la Pression 0D au pas de temps N µµµµµµµµµµµµµµµµµµµ 2086 : *----------------------------------------------------------------------- 2087 : * Execution d'une eventuelle procedure PERSO 2088 : Si TBT.'PERSO' ; 2089 : 'SI' ( TBT.'IMPR' >EG 2) ; 2090 : mess ' ***********************************'; 2091 : mess ' * Execution d une procédure PERSO *'; 2092 : mess ' * '(rxt.'PRCPERSO') 2093 : ' *'; 2094 : mess ' ***********************************'; 2095 : 'FINSI' ; 2096 : (TEXT (rxt.'PRCPERSO')) Tps rxt ; 2097 : 'SI' ( TBT.'IMPR' >EG 2) ; 2098 : mess ' ***********************************'; 2099 : mess ' * Fin execution de la procédure *'; 2100 : mess ' * '(rxt.'PRCPERSO') 2101 : ' *'; 2102 : mess ' ***********************************'; 2103 : 'FINSI' ; 2104 : Finsi ; 2105 : 2106 : *----------------------------------------------------------------------- 2107 : * RESOLUTION DE LA QDM 2108 : * Calcul du terme source de l'equation de QDM 2109 : * Ajout de la CLIM injection 2110 : 2111 : tic.'DSRC' = ((-1.) '*' (src + src2 + src3) '/' VTotal) ; 2112 : 2113 : ** contribution source densité des recombineurs --- 2114 : 'SI' (TBT.'TRECOMB') ; 2115 : tic.'DSRC' = tic.'DSRC' '+' Srho0 ; 2116 : 'FINSI' ; 2117 : ** ------------------------------------------------ 2118 : 2119 : rv = TBT.'RV' ; 2120 : 2121 : gjo matxx= 'KOPS' 'MATRIK' ; 2122 : Si (TBT.'Tbreche') ; 2123 : gjo= gjo et gj ; 2124 : Finsi ; 2125 : Si (TBT.'Tbreche2') ; 2126 : gjo= gjo et gj2 ; 2127 : Finsi ; 2128 : Si (TBT.'Tbreche3') ; 2129 : gjo= gjo et gj3 ; 2130 : Finsi ; 2131 : 2132 : 'SI'((EGA TBT.'MODTURB' 'KEPSILON') et (EGA TBT.'FEF' 'EFM1')); 2133 : *---- Modèle K-Epsilon EFM1 -------------------------------- 2134 : *mess 'Modèle K-Epsilon EFM1' ; 2135 : Si (TBT.'Tbreche') ; 2136 : lcj=extr gj 'COMP' ; 2137 : mdj=gj 'PSCA' gj lcj lcj ; 2138 : gk=((0.05**2.)*mdj) nomc 'KN' ; 2139 : gk=CHAN 'ATTRIBUT' gk 'NATURE' 'DISCRET' ; 2140 : ge=((gk**1.5)/(GEO.'Dbreche')) nomc 'EN' ; 2141 : gjo= gjo et gk ; 2142 : gjo= gjo et ge ; 2143 : Finsi ; 2144 : Si (TBT.'Tbreche2') ; 2145 : lcj=extr gj2 'COMP' ; 2146 : mdj=gj2 'PSCA' gj2 lcj lcj ; 2147 : gk=((0.05**2.)*mdj) nomc 'KN' ; 2148 : gk=CHAN 'ATTRIBUT' gk 'NATURE' 'DISCRET' ; 2149 : ge=((gk**1.5)/(GEO.'Dbreche2')) nomc 'EN' ; 2150 : gjo= gjo et gk ; 2151 : gjo= gjo et ge ; 2152 : Finsi ; 2153 : Si (TBT.'Tbreche3') ; 2154 : lcj=extr gj3 'COMP' ; 2155 : mdj=gj3 'PSCA' gj3 lcj lcj ; 2156 : gk=((0.05**2.)*mdj) nomc 'KN' ; 2157 : gk=CHAN 'ATTRIBUT' gk 'NATURE' 'DISCRET' ; 2158 : ge=((gk**1.5)/(GEO.'Dbreche3')) nomc 'EN' ; 2159 : gjo= gjo et gk ; 2160 : gjo= gjo et ge ; 2161 : Finsi ; 2162 : 2163 : *---- Modèle K-Epsilon EFM1 -- FIN ------------------------- 2164 : 'FINSI' ; 2165 : 2166 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')) ; 2167 : rv.'CLIM' = (rv.'CLIM' '-' rv.'CLIM') '+' gjo ; 2168 : Finsi ; 2169 : 2170 : ** contribution CL vitesse des recombineurs --- 2171 : 'SI' (TBT.'TRECOMB') ; 2172 : 'SI' ((TBT.'Tbreche') 'OU' (TBT.'Tbreche2') 'OU' (TBT.'Tbreche3')) ; 2173 : rv.'CLIM' = rv.'CLIM' '+' Urec0 ; 2174 : 'SINON' ; 2175 : rv.'CLIM' = (rv.'CLIM' '-' rv.'CLIM') '+' Urec0 ; 2176 : 'FINSI' ; 2177 : 'FINSI' ; 2178 : ** -------------------------------------------- 2179 : 2180 : 2181 : Rgp = tic.'Rgp' ; 2182 : 2183 : 'SI' ( TBT.'IMPR' >EG 2) ; 2184 : 'MESS' '==> RESOLUTION QDM' ; 2185 : 'FINSI' ; 2186 : 2187 : *--- Actualisation de NU 2188 : tic.'NU' = (tic.'Mu' '/' Rhomnm) ; 2189 : tic.'NUm' = somt (GEO.'Diag' * tic.'NU') / GEO.'VTotal' ; 2190 : 2191 : *--- On actualise YP (si Fonction de paroi) 2192 : Si (TBT.'FPAROI'); 2193 : tic.'YP'=rxt.'YP'; 2194 : Finsi ; 2195 : 2196 : *--- Actualisation de NUEFF 2197 : Si (Exist TBT 'MODTURB') ; 2198 : 2199 : *---- Modèle viscosité tourbillonnaire constante ----------- 2200 : Si ( Ega TBT.'MODTURB' 'NUTURB') ; 2201 : tic.'NUEFF'=tic.'NU' + rxt.'NUT' ; 2202 : Finsi ; 2203 : *---- Modèle viscosité tourbillonnaire constante -- FIN ---- 2204 : 2205 : *---- Modèle de longueur de mélange ------------------------ 2206 : Si ( Ega TBT.'MODTURB' 'LMEL') ; 2207 : un = tic.'UN' ; 2208 : P = PRODT UN $vtf ; 2209 : lm=TBT.'LMEL' ; 2210 : NUT= (lm*lm*(P**0.5)) ; 2211 : *Réalisabilité sur Nut 2212 : Nu=tic.'NU' ; 2213 : a = 5000.*Nu ; 2214 : al=0.8 ; ala=al*a ; b=ala*((2.*(1.-al))**(-1.)) ; 2215 : ik=masq Nut 'INFERIEUR' ala ; 2216 : Nut=(ik*Nut)+( (1.-ik)*a*(Nut + b)*(inve (a + Nut + b))); 2217 : Nutsnu=Nut*(inve Nu) ; 2218 : 'SI' ( TBT.'IMPR' >EG 2) ; 2219 : MESS ' Mini Maxi Nut/Nu ' (mini Nutsnu) (maxi Nutsnu) ; 2220 : 'FINSI' ; 2221 : tic.'NUEFF'=NUT + tic.'NU' ; 2222 : 2223 : Finsi ; 2224 : *---- Modèle de longueur de mélange -- FIN ----------------- 2225 : 2226 : *---- Modèle K-Epsilon Implicite --------------------------- 2227 : Si (( Ega TBT.'MODTURB' 'KEPSILON') et (non(EGA TBT.'FEF' 'EFM1'))); 2228 : * Calcul KEPSILON 2229 : ************** Constantes du K epsilon ****** 2230 : cnu=0.09; 2231 : c1=1.44; 2232 : c2=1.92; 2233 : sgk=1. ; 2234 : sge=1.3; 2235 : 2236 : Dg = doma $vtf 'XXDIAGSI' ; 2237 : Lma=(somt Dg)**(1./(vale DIME)); 2238 : Lmi=(mini Dg)**((1.)/(vale DIME)); 2239 : un = tic.'UN' ; 2240 : P = PRODT UN $vtf ; 2241 : * P= PRODT UN $vtf m4 m5 ; 2242 : 2243 : Rrk= TBT.'Rrk' ; 2244 : Rre= TBT.'Rre' ; 2245 : 2246 : En=tic.'EN' ; 2247 : Kn=tic.'KN' ; 2248 : *Nut=Cnu * Kn * Kn * (inve En) ; 2249 : Nut = (Lmi*Lma*(P**0.5)) ; 2250 : * Réalisabilité sur Nut 2251 : * Nut > 1.e-15 2252 : ik=masq Nut 'SUPERIEUR' 1.e-15 ; 2253 : Nut=(ik*Nut)+((1.-ik)*1.e-15); 2254 : * Nut P < En 2255 : *a= En * (inve P) * 0.133; 2256 : a= En * (inve P) * 0.133; 2257 : al=0.8; ala=al*a; b=ala*((2.*(1.-al))**(-1.)); 2258 : ik=masq Nut 'INFERIEUR' ala ; 2259 : *Nut=(ik*Nut)+( (1.-ik)*a*(Nut + b)*(inve (a + Nut + b))); 2260 : *Nut = kcht $vtf scal sommet Nut ; 2261 : mess 'mini maxi Nut ' (mini Nut ) (maxi Nut ) ; 2262 : Sk=Nut * P - ((KN**1.5)*(1./Lma)) ; 2263 : Kn = 'KOPS' Kn '|<' 1.e-5 ; 2264 : EsK= En * (inve Kn) ; 2265 : Se=EsK*(Nut * P - (c2*En)); 2266 : tic.'Se'=Se; 2267 : tic.'Sk'=Sk; 2268 : 2269 : tic.'NUEFF'=NUT + tic.'NU' ; 2270 : *-- On traite l'injection comme une C.L. en température 2271 : cpv matxx= 'KOPS' 'MATRIK' ; 2272 : Si (TBT.'Tbreche') ; 2273 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'KN' 2274 : Knj 'NATURE' 'DISCRET') ; 2275 : Finsi ; 2276 : Si (TBT.'Tbreche2') ; 2277 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'KN' 2278 : Knj2 'NATURE' 'DISCRET') ; 2279 : Finsi ; 2280 : Si (TBT.'Tbreche3') ; 2281 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'KN' 2282 : Knj3 'NATURE' 'DISCRET') ; 2283 : Finsi ; 2284 : 2285 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2286 : Rrk.'CLIM' = cpv; 2287 : Finsi ; 2288 : 2289 : 2290 : Rrk.'METHINV'.'CALPREC' = calprec ; 2291 : EXEC Rrk ; 2292 : * Filtre 0 '<' Kn 2293 : Kn = tic.'KN' ; 2294 : Kn = 'KOPS' Kn '|<' 1.e-5 ; 2295 : tic.'KN' = Kn ; 2296 : 2297 : 2298 : 2299 : *-- On traite l'injection comme une C.L. en température 2300 : cpv matxx= 'KOPS' 'MATRIK' ; 2301 : Si (TBT.'Tbreche') ; 2302 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'EN' 2303 : Enj 'NATURE' 'DISCRET') ; 2304 : Finsi ; 2305 : Si (TBT.'Tbreche2') ; 2306 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'EN' 2307 : Enj2 'NATURE' 'DISCRET') ; 2308 : Finsi ; 2309 : Si (TBT.'Tbreche3') ; 2310 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'EN' 2311 : Enj3 'NATURE' 'DISCRET') ; 2312 : Finsi ; 2313 : 2314 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2315 : Rre.'CLIM' = cpv; 2316 : Finsi ; 2317 : 2318 : 2319 : Rre.'METHINV'.'CALPREC' = calprec ; 2320 : EXEC Rre ; 2321 : * Filtre 0 '<' En 2322 : En = tic.'EN' ; 2323 : En = 'KOPS' En '|<' 1.e-15 ; 2324 : tic.'EN' = En ; 2325 : 2326 : * Réalisabilité sur Nut 2327 : Nut = cnu * Kn * Kn * (inve En) ; ; 2328 : Nu= tic.'NU' ; 2329 : mess 'Réalisabilité sur Nut'; 2330 : a = 1500.*Nu ; al=0.8 ; ala=al*a ; b=ala*((2.*(1.-al))**(-1.)); 2331 : tic.'ala'=ala ; 2332 : tic.'nut'=nut ; 2333 : ik=masq Nut 'INFERIEUR' ala ; 2334 : Nut=(ik*Nut)+( (1.-ik)*a*(Nut + b)*(inve (a + Nut + b))); 2335 : Nut = kcht $vtf scal sommet Nut ; 2336 : Nutsnu=(1./Nu)*Nut; 2337 : 'SI' ( TBT.'IMPR' >EG 2) ; 2338 : MESS ' Mini Maxi Nut/Nu ' (mini Nutsnu) (maxi Nutsnu) ; 2339 : 'FINSI' ; 2340 : 2341 : tic.'NUEFF'=NUT + tic.'NU' ; 2342 : 2343 : 2344 : Finsi ; 2345 : *---- Modèle K-Epsilon Implicite -- FIN -------------------- 2346 : 2347 : *---- Modèle K-Epsilon EFM1 -------------------------------- 2348 : 'SI'(( Ega TBT.'MODTURB' 'KEPSILON') et (EGA TBT.'FEF' 'EFM1')); 2349 : tic.'NUEFF'=tic.'NUT' + tic.'NU' ; 2350 : *---- Modèle K-Epsilon EFM1 -- FIN ------------------------- 2351 : Finsi ; 2352 : 2353 : Sinon ; 2354 : *---- Il n'y a pas de modèle de turbulence ----------------- 2355 : tic.'NUEFF'=tic.'NU' ; 2356 : Finsi ; 2357 : 2358 : *--- Actualisation de MUT 2359 : tic.'MUT' = (tic.'NUEFF'* Rhomnm) ; 2360 : *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2361 : 2362 : * Mess ' Avt EXEC RV tic.DT=' tic.'DT' ; 2363 : 2364 : rv.'METHINV'.'CALPREC' = calprec ; 2365 : 'SI' ('EXIS' rv 'PROJ') ; 2366 : rv.'PROJ'.'METHINV'.'CALPREC' = calprec ; 2367 : 'FINSI' ; 2368 : EXEC rv ; 2369 : 2370 : * Modif. FP 18/05/00 2371 : * Détermination des termes sources liés à la condensation en masse : 2372 : * ------------------------------------------------------------------ 2373 : 2374 : *--------------------------------------------------- 2375 : * RESOLUTION DE LA DENSITE DE VAPEUR 2376 : * Affichage des masses totales vapeur AVANT et APRES 2377 : * Creation de la CLIM pour injection 2378 : 2379 : 'SI' TBT.'VAPEUR' ; 2380 : Rrvap = TBT.'Rrvap' ; 2381 : 2382 : *-- On traite l'injection comme une C.L. en densite 2383 : 2384 : cpv matxx= 'KOPS' 'MATRIK' ; 2385 : Si (TBT.'Tbreche') ; 2386 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RVP' 2387 : Rvapj 'NATURE' 'DISCRET') ; 2388 : Finsi ; 2389 : Si (TBT.'Tbreche2') ; 2390 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RVP' 2391 : Rvapj2 'NATURE' 'DISCRET') ; 2392 : Finsi ; 2393 : Si (TBT.'Tbreche3') ; 2394 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RVP' 2395 : Rvapj3 'NATURE' 'DISCRET') ; 2396 : Finsi ; 2397 : 2398 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2399 : Rrvap.'CLIM' = cpv; 2400 : Finsi ; 2401 : 2402 : 'SI' TBT.'TRECOMB' ; 2403 : cpvrec matxx= 'KOPS' 'MATRIK' ; 2404 : 2405 : NBREC = 'DIME' rxt.'RECOMB' ; 2406 : 'REPETER' BCLREC NBREC ; 2407 : NAME = 'CHAIN' 'REC' 'sorI' &BCLREC ; 2408 : SORTIEi = GEO.NAME ; 2409 : NRECT = 'DIME' RXT.'TIC'.'RECOMB'.&BCLREC.'RH2OS' ; 2410 : RH2Os = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'RH2OS' NRECT ; 2411 : 2412 : cpvrec = cpvrec 'ET' ('MANU' 'CHPO' SORTIEi 1 'RVP' 2413 : RH2Os 'NATURE' 'DISCRET') ; 2414 : 'FIN' BCLREC ; 2415 : 2416 : 'SI' ((TBT.'Tbreche') 'OU' (TBT.'Tbreche2') 'OU' (TBT.'Tbreche3')) ; 2417 : Rrvap.'CLIM' = Rrvap.'CLIM' '+' cpvrec ; 2418 : 'SINON' ; 2419 : Rrvap.'CLIM' = cpvrec ; 2420 : 'FINSI' ; 2421 : 'FINSI' ; 2422 : 2423 : 2424 : 'SI' ( TBT.'IMPR' >EG 2) ; 2425 : 'MESS' '==> RESOLUTION Rrvap' ; 2426 : 'FINSI' ; 2427 : 2428 : Rrvap.'METHINV'.'CALPREC' = calprec ; 2429 : EXEC Rrvap ; 2430 : 2431 : 'SINON' ; 2432 : Rvp = 0. ; 2433 : tic.'RVP' = Rvp ; 2434 : 'FINSI' ; 2435 : tic.'ROVP' = tic.'RVP' ; 2436 : 2437 : *------------------------------------------------ 2438 : * RESOLUTION DE LA DENSITE HELIUM 2439 : * 2440 : 'SI' TBT.'THE' ; 2441 : Rrhe = TBT.'Rrhe' ; 2442 : 2443 : *-- On traite l'injection comme une C.L. en densite 2444 : cpv matxx= 'KOPS' 'MATRIK' ; 2445 : Si (TBT.'Tbreche') ; 2446 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RHE' 2447 : Rhej 'NATURE' 'DISCRET') ; 2448 : Finsi ; 2449 : Si (TBT.'Tbreche2') ; 2450 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RHE' 2451 : Rhej2 'NATURE' 'DISCRET') ; 2452 : Finsi ; 2453 : Si (TBT.'Tbreche3') ; 2454 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RHE' 2455 : Rhej3 'NATURE' 'DISCRET') ; 2456 : Finsi ; 2457 : 2458 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2459 : Rrhe.'CLIM' = cpv; 2460 : Finsi ; 2461 : 2462 : 'SI' ( TBT.'IMPR' >EG 2) ; 2463 : 'MESS' '==> RESOLUTION Rrhe' ; 2464 : 'FINSI' ; 2465 : 2466 : Rrhe.'resmn' = tic.'RHE' NOMC 'RHE' ; 2467 : 2468 : Rrhe.'METHINV'.'CALPREC' = calprec ; 2469 : EXEC Rrhe ; 2470 : 2471 : 2472 : * Filtre 0 '<' Rhe 2473 : rhe = tic.'RHE' ; 2474 : rhe = 'KOPS' rhe '|<' 1.e-10 ; 2475 : rhe = 'KOPS' rhe '>|' (maxi tic.'RHO') ; 2476 : 'SINON' ; 2477 : rhe = 0. ; 2478 : 'FINSI' ; 2479 : 2480 : *------------------------------------------------ 2481 : * RESOLUTION DE LA DENSITE HYDROGENE 2482 : * 2483 : 'SI' TBT.'TH2' ; 2484 : Rrh2 = TBT.'Rrh2' ; 2485 : 2486 : *-- On traite l'injection comme une C.L. en densite 2487 : cpv matxx= 'KOPS' 'MATRIK' ; 2488 : Si (TBT.'Tbreche') ; 2489 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RH2' 2490 : Rh2j 'NATURE' 'DISCRET') ; 2491 : Finsi ; 2492 : Si (TBT.'Tbreche2') ; 2493 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RH2' 2494 : Rh2j2 'NATURE' 'DISCRET') ; 2495 : Finsi ; 2496 : Si (TBT.'Tbreche3') ; 2497 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RH2' 2498 : Rh2j3 'NATURE' 'DISCRET') ; 2499 : Finsi ; 2500 : 2501 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2502 : Rrh2.'CLIM' = cpv ; 2503 : Finsi ; 2504 : 2505 : 'SI' TBT.'TRECOMB' ; 2506 : cpvrec matxx= 'KOPS' 'MATRIK' ; 2507 : 2508 : NBREC = 'DIME' rxt.'RECOMB' ; 2509 : 'REPETER' BCLREC NBREC ; 2510 : NAME = 'CHAIN' 'REC' 'sorI' &BCLREC ; 2511 : SORTIEi = GEO.NAME ; 2512 : NRECT = 'DIME' RXT.'TIC'.'RECOMB'.&BCLREC.'RH2S' ; 2513 : RH2s = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'RH2S' NRECT ; 2514 : 2515 : cpvrec = cpvrec 'ET' ('MANU' 'CHPO' SORTIEi 1 'RH2' 2516 : RH2s 'NATURE' 'DISCRET') ; 2517 : 'FIN' BCLREC ; 2518 : 2519 : 'SI' ((TBT.'Tbreche') 'OU' (TBT.'Tbreche2') 'OU' (TBT.'Tbreche3')) ; 2520 : Rrh2.'CLIM' = Rrh2.'CLIM' '+' cpvrec ; 2521 : 'SINON' ; 2522 : Rrh2.'CLIM' = cpvrec ; 2523 : 'FINSI' ; 2524 : 'FINSI' ; 2525 : 2526 : 'SI' ( TBT.'IMPR' >EG 2) ; 2527 : 'MESS' '==> RESOLUTION Rrh2' ; 2528 : 'FINSI' ; 2529 : 2530 : Rrh2.'resmn' = tic.'RH2' NOMC 'RH2' ; 2531 : 2532 : Rrh2.'METHINV'.'CALPREC' = calprec ; 2533 : EXEC Rrh2 ; 2534 : 2535 : * Filtre 0 '<' Rh2 2536 : rh2 = tic.'RH2' ; 2537 : rh2 = 'KOPS' rh2 '|<' 1.e-10 ; 2538 : rh2 = 'KOPS' rh2 '>|' (maxi tic.'RHO') ; 2539 : 'SINON' ; 2540 : rh2 = 0. ; 2541 : 'FINSI' ; 2542 : 2543 : *------------------------------------------------ 2544 : * RESOLUTION DE LA DENSITE OXYGENE 2545 : * 2546 : 'SI' TBT.'TO2' ; 2547 : Rro2 = TBT.'Rro2' ; 2548 : 2549 : *-- On traite l'injection comme une C.L. en densite 2550 : 2551 : cpv matxx= 'KOPS' 'MATRIK' ; 2552 : Si (TBT.'Tbreche') ; 2553 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RO2' 2554 : Ro2j 'NATURE' 'DISCRET') ; 2555 : Finsi ; 2556 : Si (TBT.'Tbreche2') ; 2557 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RO2' 2558 : Ro2j2 'NATURE' 'DISCRET') ; 2559 : Finsi ; 2560 : Si (TBT.'Tbreche3') ; 2561 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RO2' 2562 : Ro2j3 'NATURE' 'DISCRET') ; 2563 : Finsi ; 2564 : 2565 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2566 : Rro2.'CLIM' = cpv ; 2567 : Finsi ; 2568 : 2569 : 'SI' TBT.'TRECOMB' ; 2570 : cpvrec matxx= 'KOPS' 'MATRIK' ; 2571 : 2572 : NBREC = 'DIME' rxt.'RECOMB' ; 2573 : 'REPETER' BCLREC NBREC ; 2574 : NAME = 'CHAIN' 'REC' 'sorI' &BCLREC ; 2575 : SORTIEi = GEO.NAME ; 2576 : NRECT = 'DIME' RXT.'TIC'.'RECOMB'.&BCLREC.'RO2S' ; 2577 : RO2s = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'RO2S' NRECT ; 2578 : 2579 : cpvrec = cpvrec 'ET' ('MANU' 'CHPO' SORTIEi 1 'RO2' 2580 : RO2s 'NATURE' 'DISCRET') ; 2581 : 'FIN' BCLREC ; 2582 : 2583 : 'SI' ((TBT.'Tbreche') 'OU' (TBT.'Tbreche2') 'OU' (TBT.'Tbreche3')) ; 2584 : Rro2.'CLIM' = Rro2.'CLIM' '+' cpvrec ; 2585 : 'SINON' ; 2586 : Rro2.'CLIM' = cpvrec ; 2587 : 'FINSI' ; 2588 : 'FINSI' ; 2589 : 2590 : 'SI' ( TBT.'IMPR' >EG 2) ; 2591 : 'MESS' '==> RESOLUTION Rro2' ; 2592 : 'FINSI' ; 2593 : 2594 : Rro2.'resmn' = tic.'RO2' NOMC 'RO2' ; 2595 : 2596 : Rro2.'METHINV'.'CALPREC' = calprec ; 2597 : EXEC Rro2 ; 2598 : 2599 : * Filtre 0 '<' Ro2 2600 : ro2 = tic.'RO2' ; 2601 : ro2 = 'KOPS' ro2 '|<' 1.e-10 ; 2602 : ro2 = 'KOPS' ro2 '>|' (maxi tic.'RHO') ; 2603 : 'SINON' ; 2604 : ro2 = 0. ; 2605 : 'FINSI' ; 2606 : 2607 : *------------------------------------------------ 2608 : * RESOLUTION DE LA DENSITE AZOTE 2609 : * 2610 : 'SI' TBT.'TN2' ; 2611 : Rrn2 = TBT.'Rrn2' ; 2612 : 2613 : *-- On traite l'injection comme une C.L. en densite 2614 : 2615 : cpv matxx= 'KOPS' 'MATRIK' ; 2616 : Si (TBT.'Tbreche') ; 2617 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RN2' 2618 : Rn2j 'NATURE' 'DISCRET') ; 2619 : Finsi ; 2620 : Si (TBT.'Tbreche2') ; 2621 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RN2' 2622 : Rn2j2 'NATURE' 'DISCRET') ; 2623 : Finsi ; 2624 : Si (TBT.'Tbreche3') ; 2625 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RN2' 2626 : Rn2j3 'NATURE' 'DISCRET') ; 2627 : Finsi ; 2628 : 2629 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2630 : Rrn2.'CLIM' = cpv ; 2631 : Finsi ; 2632 : 2633 : 'SI' TBT.'TRECOMB' ; 2634 : cpvrec matxx= 'KOPS' 'MATRIK' ; 2635 : 2636 : NBREC = 'DIME' rxt.'RECOMB' ; 2637 : 'REPETER' BCLREC NBREC ; 2638 : NAME = 'CHAIN' 'REC' 'sorI' &BCLREC ; 2639 : SORTIEi = GEO.NAME ; 2640 : NRECT = 'DIME' RXT.'TIC'.'RECOMB'.&BCLREC.'RN2S' ; 2641 : RN2s = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'RN2S' NRECT ; 2642 : 2643 : cpvrec = cpvrec 'ET' ('MANU' 'CHPO' SORTIEi 1 'RN2' 2644 : RN2s 'NATURE' 'DISCRET') ; 2645 : 'FIN' BCLREC ; 2646 : 2647 : 'SI' ((TBT.'Tbreche') 'OU' (TBT.'Tbreche2') 'OU' (TBT.'Tbreche3')) ; 2648 : Rrn2.'CLIM' = Rrn2.'CLIM' '+' cpvrec ; 2649 : 'SINON' ; 2650 : Rrn2.'CLIM' = cpvrec ; 2651 : 'FINSI' ; 2652 : 'FINSI' ; 2653 : 2654 : 'SI' ( TBT.'IMPR' >EG 2) ; 2655 : 'MESS' '==> RESOLUTION Rrn2' ; 2656 : 'FINSI' ; 2657 : 2658 : Rrn2.'resmn' = tic.'RN2' NOMC 'RN2' ; 2659 : 2660 : Rrn2.'METHINV'.'CALPREC' = calprec ; 2661 : EXEC Rrn2 ; 2662 : 2663 : * Filtre 0 '<' Rn2 2664 : rn2 = tic.'RN2' ; 2665 : rn2 = 'KOPS' rn2 '|<' 1.e-10 ; 2666 : rn2 = 'KOPS' rn2 '>|' (maxi tic.'RHO') ; 2667 : 'SINON' ; 2668 : rn2 = 0. ; 2669 : 'FINSI' ; 2670 : 2671 : *------------------------------------------------ 2672 : * RESOLUTION DE LA DENSITE DIOXYDE DE CARBONE 2673 : * 2674 : 'SI' TBT.'TCO2' ; 2675 : Rrco2 = TBT.'Rrco2' ; 2676 : 2677 : *-- On traite l'injection comme une C.L. en densite 2678 : 2679 : cpv matxx= 'KOPS' 'MATRIK' ; 2680 : Si (TBT.'Tbreche') ; 2681 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RCO2' 2682 : Rco2j 'NATURE' 'DISCRET') ; 2683 : Finsi ; 2684 : Si (TBT.'Tbreche2') ; 2685 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RCO2' 2686 : Rco2j2 'NATURE' 'DISCRET') ; 2687 : Finsi ; 2688 : Si (TBT.'Tbreche3') ; 2689 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RCO2' 2690 : Rco2j3 'NATURE' 'DISCRET') ; 2691 : Finsi ; 2692 : 2693 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2694 : Rrco2.'CLIM' = cpv ; 2695 : Finsi ; 2696 : 2697 : 'SI' ( TBT.'IMPR' >EG 2) ; 2698 : 'MESS' '==> RESOLUTION Rrco2' ; 2699 : 'FINSI' ; 2700 : 2701 : Rrco2.'resmn' = tic.'RCO2' NOMC 'RCO2' ; 2702 : 2703 : Rrco2.'METHINV'.'CALPREC' = calprec ; 2704 : EXEC Rrco2 ; 2705 : 2706 : * Filtre 0 '<' Rco2 2707 : rco2 = tic.'RCO2' ; 2708 : rco2 = 'KOPS' rco2 '|<' 1.e-10 ; 2709 : rco2 = 'KOPS' rco2 '>|' (maxi tic.'RHO') ; 2710 : 'SINON' ; 2711 : rco2 = 0. ; 2712 : 'FINSI' ; 2713 : 2714 : *? TCO2AIR 2715 : *------------------------------------------------ 2716 : * RESOLUTION DE LA DENSITE D'AIR 2717 : * 2718 : 'SI' TBT.'TAIR' ; 2719 : Rrair = TBT.'Rrair' ; 2720 : 2721 : *-- On traite l'injection comme une C.L. en densite 2722 : 2723 : cpv matxx= 'KOPS' 'MATRIK' ; 2724 : Si (TBT.'Tbreche') ; 2725 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RAIR' 2726 : Rairj 'NATURE' 'DISCRET') ; 2727 : Finsi ; 2728 : Si (TBT.'Tbreche2') ; 2729 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RAIR' 2730 : Rairj2 'NATURE' 'DISCRET') ; 2731 : Finsi ; 2732 : Si (TBT.'Tbreche3') ; 2733 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RAIR' 2734 : Rairj3 'NATURE' 'DISCRET') ; 2735 : Finsi ; 2736 : 2737 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2738 : Rrair.'CLIM' = cpv ; 2739 : Finsi ; 2740 : 2741 : 'SI' ( TBT.'IMPR' >EG 2) ; 2742 : 'MESS' '==> RESOLUTION Rrair' ; 2743 : 'FINSI' ; 2744 : 2745 : Rrair.'resmn' = tic.'RAIR' NOMC 'RAIR' ; 2746 : 2747 : Rrair.'METHINV'.'CALPREC' = calprec ; 2748 : EXEC Rrair ; 2749 : 2750 : * Filtre 0 '<' Rair 2751 : rair = tic.'RAIR' ; 2752 : rair = 'KOPS' rair '|<' 1.e-10 ; 2753 : rair = 'KOPS' rair '>|' (maxi tic.'RHO') ; 2754 : 'FINSI' ; 2755 : 2756 : *------------------------------------------------ 2757 : * RESOLUTION DE LA DENSITE MONOXYDE DE CARBONE 2758 : * 2759 : 'SI' TBT.'TCO' ; 2760 : Rrco = TBT.'Rrco' ; 2761 : 2762 : *-- On traite l'injection comme une C.L. en densite 2763 : 2764 : cpv matxx= 'KOPS' 'MATRIK' ; 2765 : Si (TBT.'Tbreche') ; 2766 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'RCO' 2767 : Rcoj 'NATURE' 'DISCRET') ; 2768 : Finsi ; 2769 : Si (TBT.'Tbreche2') ; 2770 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'RCO' 2771 : Rcoj2 'NATURE' 'DISCRET') ; 2772 : Finsi ; 2773 : Si (TBT.'Tbreche3') ; 2774 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'RCO' 2775 : Rcoj3 'NATURE' 'DISCRET') ; 2776 : Finsi ; 2777 : 2778 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2779 : Rrco.'CLIM' = cpv ; 2780 : Finsi ; 2781 : 2782 : 'SI' ( TBT.'IMPR' >EG 2) ; 2783 : 'MESS' '==> RESOLUTION Rrco' ; 2784 : 'FINSI' ; 2785 : 2786 : Rrco.'resmn' = tic.'RCO' NOMC 'RCO' ; 2787 : 2788 : Rrco.'METHINV'.'CALPREC' = calprec ; 2789 : EXEC Rrco ; 2790 : 2791 : * Filtre 0 '<' Rco 2792 : rco = tic.'RCO' ; 2793 : rco = 'KOPS' rco '|<' 1.e-10 ; 2794 : rco = 'KOPS' rco '>|' (maxi tic.'RHO') ; 2795 : 'SINON' ; 2796 : rco = 0. ; 2797 : 'FINSI' ; 2798 : 2799 : *------------------------------------------------ 2800 : * RESOLUTION DE LA TEMPERATURE DU FLUIDE TF 2801 : * Definition des températures TBPW et TBP pour ECHI 2802 : * TBPW température paroi (THERMP) 2803 : * Definition des CLIM pour l'injection et ajout 2804 : * des CLIM pour température paroi imposée 2805 : * Le terme source STF est calcule avant (fin pas precedent) 2806 : 2807 : rtf = TBT.'RTF' ; 2808 : 2809 : 'SI' TBT.'TPAROIF' ; 2810 : TBPW = 'REDU' (tic.'TP') GEO.'paroif' ; 2811 : tic.'TBPW' = 'KCHT' GEO.'$paroif' 'SCAL' 'SOMMET' TBPW ; 2812 : TBFW = 'REDU' (tic.'TF') GEO.'paroif' ; 2813 : tic.'TBFW' = 'KCHT' GEO.'$paroif' 'SCAL' 'SOMMET' TBFW ; 2814 : 'FINSI' ; 2815 : 2816 : *-- On traite l'injection comme une C.L. en température 2817 : cpv matxx= 'KOPS' 'MATRIK' ; 2818 : Si (TBT.'Tbreche') ; 2819 : cpv = cpv et ('MANU' 'CHPO' brechei 1 'TF' 2820 : Tj 'NATURE' 'DISCRET') ; 2821 : Finsi ; 2822 : Si (TBT.'Tbreche2') ; 2823 : cpv = cpv et ('MANU' 'CHPO' brech2i 1 'TF' 2824 : Tj2 'NATURE' 'DISCRET') ; 2825 : Finsi ; 2826 : Si (TBT.'Tbreche3') ; 2827 : cpv = cpv et ('MANU' 'CHPO' brech3i 1 'TF' 2828 : Tj3 'NATURE' 'DISCRET') ; 2829 : Finsi ; 2830 : 2831 : Si ((TBT.'Tbreche') ou (TBT.'Tbreche2') ou (TBT.'Tbreche3')); 2832 : rtf.'CLIM' = cpv; 2833 : Finsi ; 2834 : 2835 : 'SI' TBT.'TRECOMB' ; 2836 : cpvrec matxx= 'KOPS' 'MATRIK' ; 2837 : 2838 : NBREC = 'DIME' rxt.'RECOMB' ; 2839 : 'REPETER' BCLREC NBREC ; 2840 : NAME = 'CHAIN' 'REC' 'sorI' &BCLREC ; 2841 : SORTIEi = GEO.NAME ; 2842 : NRECT = 'DIME' RXT.'TIC'.'RECOMB'.&BCLREC.'Tout' ; 2843 : RTS0 = 'EXTR' RXT.'TIC'.'RECOMB'.&BCLREC.'Tout' NRECT ; 2844 : 2845 : cpvrec = cpvrec 'ET' ('MANU' 'CHPO' SORTIEi 1 'TF' 2846 : RTS0 'NATURE' 'DISCRET') ; 2847 : 'FIN' BCLREC ; 2848 : 2849 : 'SI' ((TBT.'Tbreche') 'OU' (TBT.'Tbreche2') 'OU' (TBT.'Tbreche3')) ; 2850 : rtf.'CLIM' = rtf.'CLIM' '+' cpvrec ; 2851 : 'SINON' ; 2852 : rtf.'CLIM' = cpvrec ; 2853 : 'FINSI' ; 2854 : 'FINSI' ; 2855 : 2856 : 'SI' ( TBT.'IMPR' >EG 2) ; 2857 : 'MESS' '==> RESOLUTION rtf ' ; 2858 : 'FINSI' ; 2859 : 2860 : *? TF=tic.'TF' ; 2861 : *? TFmoy=(somt(Diag*TF))/Vtotal; 2862 : *? mess '%$$ TFmoy avant exec rtf ' TFmoy; 2863 : 2864 : rtf.'METHINV'.'CALPREC' = calprec ; 2865 : EXEC rtf ; 2866 : 2867 : *? TF=tic.'TF' ; 2868 : *? TFmoy=(somt(Diag*TF))/Vtotal; 2869 : *? mess '%$$ TFmoy apres exec rtf ' TFmoy; 2870 : 2871 : *-------------------------------------------------------- 2872 : * RESOLUTION DE LA TEMPERATURE PAROI TP 2873 : * Le terme SOURCE FHP1 est calulé dans la procédure CONDENS 2874 : 2875 : 'SI' TBT.'THERMP' ; 2876 : 'SI' ('NON' TBT.'THERCO') ; 2877 : rtp = TBT.'RTP' ; 2878 : 'SI' ( TBT.'IMPR' >EG 2) ; 2879 : 'MESS' '==> RESOLUTION rtp ' ; 2880 : 'FINSI' ; 2881 : rtp.'METHINV'.'CALPREC' = calprec ; 2882 : EXEC rtp ; 2883 : 'FINSI' ; 2884 : Tp = tic.'TP' ; 2885 : 'FINSI' ; 2886 : 2887 : Tfgd = 1.E10 ; 2888 : * Filtre Tfmini '<' Tf 2889 : CHrtf1 = rtf.'CLIM' ; 2890 : MArtf1 = EXTR CHrtf1 'MAILLAGE' ; 2891 : 'SI' ( MArtf1 'EGA' 0 ) ; 2892 : tfmi0 = Tfgd ; 2893 : 'SINON' ; 2894 : tfmi0 = mini rtf.'CLIM' ; 2895 : 'FINSI' ; 2896 : 'SI' TBT.'ASPER ' ; 2897 : tfmiij= tdinj ; 2898 : 'SINON' ; 2899 : tfmiij= Tfgd ; 2900 : 'FINSI' ; 2901 : 'SI' TBT.'THERMP' ; 2902 : tfmiw = mini tic.'TP' ; 2903 : 'SINON' ; 2904 : tfmiw = Tfgd ; 2905 : 'FINSI' ; 2906 : 'SI' TBT.'TTIMP1' ; 2907 : tfmi1 = tic.'TBP1' ; 2908 : 'SINON' ; 2909 : tfmi1 = Tfgd ; 2910 : 'FINSI' ; 2911 : 'SI' TBT.'TTIMP2' ; 2912 : tfmi2 = tic.'TBP2' ; 2913 : 'SINON' ; 2914 : tfmi2 = Tfgd ; 2915 : 'FINSI' ; 2916 : 'SI' TBT.'TTIMP3' ; 2917 : tfmi3 = tic.'TBP3' ; 2918 : 'SINON' ; 2919 : tfmi3 = Tfgd ; 2920 : 'FINSI' ; 2921 : 'SI' TBT.'TECHANP' ; 2922 : tfmi4 = 'MINI' tic.'TBP0' ; 2923 : 'SINON' ; 2924 : tfmi4 = Tfgd ; 2925 : 'FINSI' ; 2926 : Tfmi = MINI (PROG tfmi0 tfmiw tfmi1 tfmi2 tfmi3 tfmi4 tfmiij 2927 : (rxt.'TF0')); 2928 : Tf = tic.'TF' ; 2929 : Tf = 'KOPS' Tf '|<' tfmi ; 2930 : tic.'TF' = Tf ; 2931 : * Tf est bornée inférieurement par la température minimale 2932 : * physique 2933 : 2934 : 'SI' TBT.'ASPER ' ; 2935 : *** RVN *** Velocita gocce 2936 : rvn = TBT.'RVN' ; 2937 : rvn.'CLIM' = (rvn.'CLIM' '-' rvn.'CLIM') '+' vnj ; 2938 : rvn.'METHINV'.'CALPREC' = calprec ; 2939 : EXEC rvn ; 2940 : 2941 : *** RXD *** Frazione volume fase dispersa 2942 : rxd = TBT.'RXD' ; 2943 : rxd.'CLIM' = (rxd.'CLIM' '-' rxd.'CLIM') '+' xdj ; 2944 : rxd.'METHINV'.'CALPREC' = calprec ; 2945 : EXEC rxd ; 2946 : tic.'XD' = 'KOPS' (tic.'XD') '|<' 1.e-10 ; 2947 : 2948 : *** RTD *** Temperature fase dispersa 2949 : rtd = TBT.'RTD' ; 2950 : rtd.'CLIM' = (rtd.'CLIM' '-' rtd.'CLIM') '+' tdj ; 2951 : rtd.'METHINV'.'CALPREC' = calprec ; 2952 : EXEC rtd ; 2953 : 2954 : *** RDD *** Diametro fase dispersa 2955 : rdd = TBT.'RDD' ; 2956 : rdd.'CLIM' = (rdd.'CLIM' '-' rdd.'CLIM') '+' ddj ; 2957 : rdd.'METHINV'.'CALPREC' = calprec ; 2958 : EXEC rdd ; 2959 : tic.'DD' = 'KOPS' (tic.'DD') '|<' 1.e-10 ; 2960 : 2961 : 'FINSI' ; 2962 : 2963 : **** Fin EDP ********************** 2964 : 2965 : *tc mise en commentaire de ce finsi 2966 : *'FINSI' ; 2967 : 2968 : *------------------------------------------------------------ 2969 : *>7 MISE EN PLACE DES CONTRAINTES POUR LES BILANS 0D 2970 : *>7 Contraintes sur : 2971 : *>7 - la masse totale 2972 : *>7 - l'énergie et donc la température 2973 : *>7 - les masses des differents constituants incondensables 2974 : *>7 - la masse de vapeur d'eau 2975 : * 2976 : *------------------------------------------------------------ 2977 : 'SI' TBT.'VAPEUR' ; 2978 : Si (exist tic 'RAIR'); 2979 : Mvp=somt (diag*tic.'RVP'); 2980 : Mair=somt (diag*tic.'RAIR'); 2981 : mess '%$$ RETSAT avt CORMASSE Mvp=' mvp ' Mair='Mair; 2982 : finsi; 2983 : 2984 : dRhov rvp Rhomvnt = CORMASSE 'RVP' Rhovgn 'Mdrvap' ; 2985 : Mvp=somt (diag*tic.'RVP'); 2986 : Mdrov= drhov*Vtotal; 2987 : mess '%$$ RETSAT apr CORMASSE Mvp=' mvp ' Mdrov=' Mdrov; 2988 : 'FINSI' ; 2989 : *>7 D'où RVP 2990 : 2991 : *------------------------------------------------------------ 2992 : 'SI' TBT.'THE' ; 2993 : dRhohe rhe Rhomhet = CORMASSE 'RHE' Rmhen 'Mdrhe' ; 2994 : 'FINSI' ; 2995 : 2996 : *-------------------------------------------------------- 2997 : 'SI' TBT.'TH2' ; 2998 : dRhoh2 rh2 Rhomh2t = CORMASSE 'RH2' Rmh2n 'Mdrh2' ; 2999 : 'FINSI' ; 3000 : 3001 : *-------------------------------------------------------- 3002 : 'SI' TBT.'TO2' ; 3003 : dRhoo2 ro2 Rhomo2t = CORMASSE 'RO2' Rmo2n 'Mdro2' ; 3004 : 'FINSI' ; 3005 : 3006 : *--------------------------------------------------------- 3007 : 'SI' TBT.'TN2' ; 3008 : dRhon2 rn2 Rhomn2t = CORMASSE 'RN2' Rmn2n 'Mdrn2' ; 3009 : 'FINSI' ; 3010 : 3011 : *--------------------------------------------------------- 3012 : 'SI' TBT.'TCO' ; 3013 : dRhoco rco Rhomcot = CORMASSE 'RCO' Rmcon 'Mdrco' ; 3014 : 'FINSI' ; 3015 : 3016 : *--------------------------------------------------------- 3017 : 'SI' TBT.'TCO2' ; 3018 : dRhoco2 rco2 Rhomco2t = CORMASSE 'RCO2' Rmco2n 'Mdrco2'; 3019 : 'FINSI' ; 3020 : 3021 : *--------------------------------------------------------- 3022 : *µµµµµµµµ Calcul de la Densité au pas de temps N µµµµµµµµµµµµµµµµµµµµµµµ 3023 : * Si Cavité ouverte Pression nulle ! 3024 : Si TBT.'Tsortie' ; 3025 : Pt = 0. ; 3026 : tic.'PT' = tic.'PT' 'ET' ('PROG' Pt) ; 3027 : rho = Rhom ; 3028 : Sinon ; 3029 : * Calcul de la pression 0D avec le champ Tf actualise 3030 : 3031 : irtf = 'INVE' (rgp '*' (tic.'TF' '+' 273.15)) ; 3032 : Pt = ((Rhomn-Rbrom) '*' Vtotal) '/' ('SOMT' (Diag '*' irtf)); 3033 : *mess ' Rhomn= ' Rhomn ' Pression : Pt = ' Pt ; 3034 : 3035 : * Calcul de la Densité au pas de temps N 3036 : 3037 : Si((EGA TBT.'VERSION' 0) ou (NON (TBT.'VAPEUR'))) ; 3038 : 3039 : Mess ' Version V0 ==> Loi d état gaz parfait pour la vapeur '; 3040 : rho = Pt '*' ('INVE' (Rgp '*' (Tf '+' 273.15))) ; 3041 : rhot = Pt '*' ('INVE' (Rgp '*' (Tf '+' 273.15))) ; 3042 : tic.'PT' = tic.'PT' 'ET' ('PROG' Pt) ; 3043 : Rebro = 0. ; 3044 : Mliq = 0. ; 3045 : Mlq = 0. ; 3046 : Sinon ; 3047 : *----------------------------------------------------------------------- 3048 : * Mess ' Version V1 ==> Loi d état gaz réel pour la vapeur '; 3049 : 3050 : 3051 : Si (NON (EXIST tic 'RAIR')) ; 3052 : Rho = tic.'RHO' ; 3053 : Rair = tic.'YAIR' * tic.'RHO' ; 3054 : tic.'RAIR' = Rair ; 3055 : Finsi ; 3056 : Rbro = tic.'RBRO' ; 3057 : Rair = tic.'RAIR' ; 3058 : Rhot = Rair + Rvp + Rhe + Rh2 + Ro2 + Rn2 + Rco2 + Rco - Rbro; 3059 : 3060 : Qcm Pt yvap cpvap Rhot Mliq Rebro 3061 : = RETSAT tic ndl $vtf TBT 0.5 1.e-4 Rhot Tfmi Cvmnm ; 3062 : 3063 : *?? au cas où on veut tracer Econd = Econd + Econdm ; 3064 : *?? au cas où on veut tracer Hcond = Hcond + Hcondm ; 3065 : Mair = somt (Diag*tic.'RAIR') ; 3066 : Rmairn = Mair/Vtotal ; 3067 : 3068 : tic.'PT' = tic.'PT' 'ET' ('PROG' Pt) ; 3069 : 3070 : Mlq = Mlq + Mliq ; 3071 : 3072 : 3073 : *....................................................................... 3074 : Rbro = tic.'RBRO' ; 3075 : *Rovp = tic.'RVP' - tic.'RBRO' ; 3076 : Tf = tic.'TF' ; 3077 : 3078 : * Fin traitement Gaz réel version V1 3079 : Finsi ; 3080 : 3081 : Mlq = Mlq + (extr (dime tic.'Mliqpuis') tic.'Mliqpuis') ; 3082 : tic.'Mliqpuis'= tic.'Mliqpuis' 'ET' ('PROG' Mlq) ; 3083 : 3084 : * Fin traitement enceinte fermée 3085 : Finsi ; 3086 : *µµµµµµµµ Calcul de la Densité au pas de temps N µµµµµµµµµµµµµµµµµµµµµµµ 3087 : rho = 'KCHT' $vtf 'SCAL' 'SOMMET' rhot ; 3088 : Rhomnt = ('SOMT' (Diag '*' rho)) + (Rbrom*Vtotal) ; 3089 : *mess ' Rhomn=' Rhomn ' Rhomnt=' Rhomn; 3090 : *----------------------------- 3091 : dRho = ((Rhomn '*' VTotal) '-' Rhomnt ) '/' Vtotal ; 3092 : tic.'drho' = tic.'drho' 'ET' ( 'PROG' dRho ) ; 3093 : 3094 : *---------- Recalcul pour bilan nul 3095 : *?rho1 = rho ; 3096 : drot=rho - rhot; 3097 : mess '%$$ Mini / Maxi drot=' (mini drot) (maxi drot); 3098 : rho = 'KCHT' $vtf 'SCAL' 'SOMMET' (rho '+' dRho) ; 3099 : rhot = 'KCHT' $vtf 'SCAL' 'SOMMET' rho ; 3100 : *?mess ' Recalcul pour Bilan Nul' (somt (Diag'*' rho1)) 3101 : *?(somt (Diag'*' rho)) (somt (Diag'*' tic.'RHO')) ; 3102 : 3103 : tic.'RHO' = rho ; 3104 : Rhomnt = 'SOMT' (Diag '*' rho) ; 3105 : dRho = dRho '*' Vtotal ; 3106 : 3107 : *--------------------------------------------------------- 3108 : * On corrige l'air apres le calcul de la masse total car 3109 : * dans le cas où on ne transporte pas l'air on en a besoin !! 3110 : 'SI'((NON TBT.'TAIR') et (TBT.'VERSION' < 1)) ; 3111 : * On calcule la masse d'air 3112 : *mess 'Rho ' (type Rho) 'Rvp ' (type Rvp) 'rhe ' (type rhe); 3113 : *mess 'rh2 ' (type rh2) 'ro2 ' (type ro2) 'rco2 ' (type rco2); 3114 : *mess 'rco ' (type rco) ; 3115 : Rair = Rho - Rvp - Rhe - Rh2 - Ro2 - Rn2 - Rco2 - Rco ; 3116 : tic.'RAIR' = Rair ; 3117 : 'SINON' ; 3118 : Rair = tic.'RAIR' ; 3119 : 'FINSI' ; 3120 : 3121 : *? Mair = somt (Diag*Rair) ; 3122 : 3123 : dRhoair rair Rhomairt = CORMASSE 'RAIR' Rmairn 'Mdrair'; 3124 : 3125 : *? mess '%$$ masse d air attendue ' (Rmairn '*' VTotal) 3126 : *? 'masse d air présente avant correction ' Mair 3127 : *? 'masse d air présente apres correction ' Rhomairt ; 3128 : 3129 : *--------------------------------------------------------- 3130 : * Correction de la température pour vérifier l'énergie 3131 : * ATTENTION : Cv n'est pas reactualise pour l'instant 3132 : * N 3133 : *mess 'Cvmnm=' cvmnm 'mini maxi Tf ' (mini tf) (maxi tf) ; 3134 : *mess 'mini maxi rho ' (mini rho) (maxi rho) ; 3135 : 3136 : Remnt = Cvmnm '*' ( 'SOMT' (Diag '*' rho '*' 3137 : (Tf '+' 273.15)) ) ; 3138 : 3139 : * On rajoute à l'énergie du mélange l'énergie du Brouillard 3140 : * Si il y en a 3141 : Remnt = Remnt - rebro ; 3142 : 3143 : dE = (Remn '*' Vtotal) '-' Remnt ; 3144 : 3145 : *MESS '==> Erreur Energie = ' dE ' Joules' ; 3146 : 3147 : Cvmn=Cvmnm ; 3148 : IndE = rho '*' Cvmn '*' (Tf '+' 273.15) '*' Diag 3149 : '/' Remnt ; 3150 : dE = IndE '*' dE ; 3151 : 3152 : dTf= ('INVE' (Diag '*' rho)) '*' (dE '/' Cvmn) ; 3153 : dTf = 'KCHT' $vtf 'SCAL' 'SOMMET' dTf ; 3154 : 3155 : Si (NON (TBT.'CORTEMP')) ; 3156 : *mess ' On supprime la correction sur Tf' ; 3157 : dTf = 'KCHT' $vtf 'SCAL' 'SOMMET' 0. ; 3158 : FINsi ; 3159 : 3160 : mess '%$$ dtf ' (mini dtf) (maxi dtf) ; 3161 : Tf = 'KCHT' $vtf 'SCAL' 'SOMMET' (Tf '+' dTf) ; 3162 : Tf = KOPS Tf '|<' Tfmi ; 3163 : tic.'TF' = Tf ; 3164 : tic.'MdTf' = tic.'MdTf' 'ET' ( 'PROG' ('MAXI' dTf) ) ; 3165 : tic.'mdTf' = tic.'mdTf' 'ET' ( 'PROG' ('MINI' dTf) ) ; 3166 : 3167 : *'MESS' '==>>> Mini dTf = ' ('MINI' dTf) ' Maxi dTf ' 3168 : * ( 'MAXI' dTf) ; 3169 : 3170 : *--- Recalcul pour bilan nul 3171 : Remnt = Cvmn '*' ( 'SOMT' (Diag '*' rho '*' 3172 : (Tf '+' 273.15)) ) '/' Vtotal ; 3173 : dE = 'SOMT' ( Diag '*' dE ) '/' VTotal ; 3174 : dTf= dE '/' Cvmn '/' ('SOMT'(Diag '*' rho) ) ; 3175 : 3176 : *--------------------------------------------------------- 3177 : **** Test sur la conservation energie dans le mur 3178 : 3179 : 'SI' TBT.'THERMP' ; 3180 : Diagp1 = 'DOMA' $vtp 'VOLUME' ; 3181 : Vtotalp1 = 'SOMT' Diagp1 ; 3182 : tp01 = 'NOEL' $vtp Tp ; 3183 : Emurm = ( 'SOMT' (tic.'ROCP' '*' Diagp1 '*' 3184 : (tp01 '+' 273.15 ) ) ) 3185 : '/' Vtotalp1 ; 3186 : tic.'TP' = Tp ; 3187 : 'FINSI' ; 3188 : 3189 : *--------------------------------------------------------- 3190 : * Calcul de la Température moyenne fluide 3191 : 3192 : Tfm= (Remn '*' Vtotal) 3193 : / ('SOMT' (Diag '*' (tic.'RHO') * Cvmn)) ; 3194 : Tfm=Tfm - 273.15 ; 3195 : Tfm = maxi (prog Tfm Tfmi) ; 3196 : 3197 : tic.'Tfm'= tic.'Tfm' 'ET' ('PROG' Tfm) ; 3198 : 3199 : *--------------------------------------------------------- 3200 : * Calcul de la Température moyenne paroi 3201 : 3202 : 'SI' TBT.'THERMP' ; 3203 : Tpm= 'SOMT' (GEO.'Diagp' '*' Tp) '/' GEO.'VTotalp' ; 3204 : tic.'Tpm'= tic.'Tpm' 'ET' ('PROG' Tpm) ; 3205 : 'FINSI' ; 3206 : 3207 : *--------------------------------------------------------- 3208 : **** Actualisation des proprietes physiques ************** 3209 : **** Calcul de la pression ******************************* 3210 : **** et avancement en temps ****************************** 3211 : 3212 : * Calcul Cp (Tfm) 3213 : Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair = CALCP Tfm ; 3214 : * CALCP on est en fin de boucle 3215 : 3216 : * Calcul Mu (Tkm) 3217 : Tkm=Tfm '+' 273.15 ; 3218 : Muh2 Muhe Muo2 Mun2 Muco2 Muco Muvap Muair = CALMU Tkm ; 3219 : 3220 : 3221 : *Dans ces conditions la masse totale est 3222 : * Ceci en vue de calculer des fractions massiques 3223 : * dont la somme vaut 1 qqst x 3224 : * 0<= Yi <= 1 parceque les rair Rvp etc sont >0 3225 : * Rhot diffère de Rho dans sa distribution mais 3226 : * l'intégrale doit être la même 3227 : 3228 : Si (TBT.'VERSION' '<' 1) ; 3229 : Cpvap = TBT.'Cpvap' ; 3230 : Rhot=rair + Rvp + rhe + rh2 + ro2 + rn2 + rco + rco2 ; 3231 : Finsi ; 3232 : * Sinon Rhot est donné par RETSAT (un peu plus haut) 3233 : 3234 : ermas=somt ((Rho - Rhot)*Diag) ; 3235 : 3236 : 'SI' ( TBT.'IMPR' >EG 2) ; 3237 : MESS 'ERMAS=' ermas ; 3238 : 'FINSI' ; 3239 : 3240 : * Calcul fraction massique eau 3241 : 'SI' TBT.'VAPEUR' ; 3242 : Yvap = 'KOPS' Rvp '/' Rhot ; 3243 : 'SINON' ; 3244 : Yvap = 0. ; 3245 : 'FINSI' ; 3246 : tic.'YVAP'= Yvap ; 3247 : 3248 : * Calcul fraction massique de He 3249 : 'SI' TBT.'THE' ; 3250 : Yhe = 'KOPS' rhe '/' Rhot ; 3251 : 'SINON' ; 3252 : Yhe = 0. ; 3253 : 'FINSI' ; 3254 : tic.'YHE'= Yhe ; 3255 : 3256 : * Calcul fraction massique de H2 3257 : 'SI' TBT.'TH2' ; 3258 : Yh2 = 'KOPS' rh2 '/' Rhot ; 3259 : 'SINON' ; 3260 : Yh2 = 0. ; 3261 : 'FINSI' ; 3262 : tic.'YH2'= Yh2 ; 3263 : 3264 : * Calcul fraction massique de O2 3265 : 'SI' TBT.'TO2' ; 3266 : YO2 = 'KOPS' ro2 '/' Rhot ; 3267 : 'SINON' ; 3268 : Yo2 = 0. ; 3269 : 'FINSI' ; 3270 : tic.'YO2'= Yo2 ; 3271 : 3272 : * Calcul fraction massique de N2 3273 : 'SI' TBT.'TN2' ; 3274 : YN2 = 'KOPS' rn2 '/' Rhot ; 3275 : 'SINON' ; 3276 : Yn2 = 0. ; 3277 : 'FINSI' ; 3278 : tic.'YN2'= Yn2 ; 3279 : 3280 : * Calcul fraction massique de CO2 3281 : 'SI' TBT.'TCO2' ; 3282 : YCO2 = 'KOPS' rco2 '/' Rhot ; 3283 : 'SINON' ; 3284 : Yco2 = 0. ; 3285 : 'FINSI' ; 3286 : tic.'YCO2'= Yco2 ; 3287 : 3288 : * Calcul fraction massique de CO 3289 : 'SI' TBT.'TCO' ; 3290 : YCO = 'KOPS' rco '/' Rhot ; 3291 : 'SINON' ; 3292 : Yco = 0. ; 3293 : 'FINSI' ; 3294 : tic.'YCO'= Yco ; 3295 : 3296 : * Calcul de la fraction massique d'air 3297 : Yair = 'KOPS' rair '/' Rhot ; 3298 : Yair = KCHT $vtf 'SCAL' 'SOMMET' Yair ; 3299 : tic.'YAIR'= Yair ; 3300 : 3301 : * a= yair + yvap + yh2 + yhe + yn2 + yo2 + yco2 + yco ; 3302 : * mess ' Controle de la somme des fractions massiques ' 3303 : *(mini a) (maxi a) ; 3304 : 3305 : * Calcul de la constante des gaz du mélange et du gamma 3306 : * du mélange 3307 : * Calcul du Cp du mélange 3308 : Rgp = (Rgvap '*' Yvap) '+' (Rgair '*' Yair) 3309 : '+' (Rghe '*' Yhe ) '+' (Rgh2 '*' Yh2 ) 3310 : '+' (Rgo2 '*' Yo2 ) '+' (Rgn2 '*' Yn2 ) 3311 : '+' (Rgco '*' Yco ) '+' (Rgco2 '*' Yco2 ) ; 3312 : 3313 : Rgp = 'KCHT' $vtf 'SCAL' 'SOMMET' Rgp ; 3314 : Rgpm = 'SOMT' (Diag '*' Rgp) '/' VTotal ; 3315 : * Mess ' Rgpm = ' Rgpm ; 3316 : 3317 : Cp = (Cpvap '*' Yvap) '+' (Cpair '*' Yair) 3318 : '+' (Cphe '*' Yhe ) '+' (Cph2 '*' Yh2 ) 3319 : '+' (Cpo2 '*' Yo2 ) '+' (Cpn2 '*' Yn2 ) 3320 : '+' (Cpco '*' Yco ) '+' (Cpco2 '*' Yco2 ) ; 3321 : 3322 : Cp = 'KCHT' $vtf 'SCAL' 'SOMMET' Cp ; 3323 : Cpm = 'SOMT' (Diag '*' Cp) '/' VTotal ; 3324 : Cvm = Cpm '-' Rgpm ; 3325 : Gamm = Cpm '/' Cvm ; 3326 : 3327 : Mu = (Muvap '*' Yvap) '+' (Muair '*' Yair) 3328 : '+' (Muhe '*' Yhe ) '+' (Muh2 '*' Yh2 ) 3329 : '+' (Muo2 '*' Yo2 ) '+' (Mun2 '*' Yn2 ) 3330 : '+' (Muco '*' Yco ) '+' (Muco2 '*' Yco2 ) ; 3331 : 3332 : tic.'Mu' = Mu ; 3333 : 3334 : tic.'Rgp' = Rgp ; 3335 : tic.'Rgpm' = tic.'Rgpm' 'ET' ('PROG' Rgpm) ; 3336 : tic.'Cvm' = tic.'Cvm' 'ET' ('PROG' Cvm) ; 3337 : tic.'Gamm' = tic.'Gamm' 'ET' ('PROG' Gamm) ; 3338 : tic.'Cpm' = tic.'Cpm' 'ET' ('PROG' Cpm) ; 3339 : *--------------------------------------------------------- 3340 : * Calcul du terme source de l'equation energie 3341 : 3342 : stf = DPDT '/' (Cpm '*' Rhomn) ; 3343 : tic.'STF' = stf ; 3344 : 3345 : *--------------------------------------------------------- 3346 : * Calcul de l'energie convective 3347 : 3348 : Tf=tic.'TF' ; 3349 : Econvw = 0.0 ; 3350 : Econv1 = 0.0 ; 3351 : Econv2 = 0.0 ; 3352 : Econv3 = 0.0 ; 3353 : Econv0 = 0.0 ; 3354 : 3355 : 3356 : 'SI'((TBT.'THERMP') et (TBT.'TPAROIF')) ; 3357 : Tp=tic.'TP' ; 3358 : * Energie extraite du gaz par convection 3359 : Tfpw = 'REDU' tic.'TF' GEO.'paroif' ; 3360 : Tbpw = 'REDU' tic.'TP' GEO.'paroif' ; 3361 : Hconv = KCHT GEO.'$paroif' SCAL CENTRE tic.'KHW' ; 3362 : Hconv = elno GEO.'$paroif' Hconv ; 3363 : Econvw = 'SOMT' ( GEO.'Diagpf' '*' 3364 : (Hconv '*' (Tfpw '-' Tbpw))) ; 3365 : 'FINSI' ; 3366 : 3367 : 'SI' TBT.'TTIMP1' ; 3368 : * Energie extraite du gaz par convection 3369 : Tfp1 = 'REDU' Tf (GEO.'mtp1') ; 3370 : Tbp1 = tic.'TBP1' ; 3371 : Hconv = KCHT GEO.'$mtp1' SCAL CENTRE tic.'KH1' ; 3372 : Hconv = elno GEO.'$mtp1' Hconv ; 3373 : * Diagmp1 = 'DOMA' (GEO.'$mtp1') 'XXDIAGSI' ; 3374 : Econv1 = 'SOMT' ( GEO.'Diag1f' '*' (Hconv '*' 3375 : (Tfp1 '-' Tbp1))) ; 3376 : 'FINSI' ; 3377 : 3378 : 'SI' TBT.'TTIMP2' ; 3379 : * Energie extraite du gaz par convection 3380 : Tfp2 = 'REDU' Tf (GEO.'mtp2') ; 3381 : Tbp2 = tic.'TBP2' ; 3382 : Hconv = KCHT GEO.'$mtp2' SCAL CENTRE tic.'KH2' ; 3383 : Hconv = elno GEO.'$mtp2' Hconv ; 3384 : * Diagmp2 = 'DOMA' (GEO.'$mtp2') 'XXDIAGSI' ; 3385 : Econv2 = 'SOMT' ( GEO.'Diag2f' '*' (Hconv '*' 3386 : (Tfp2 '-' Tbp2))) ; 3387 : 'FINSI' ; 3388 : 3389 : 'SI' TBT.'TTIMP3' ; 3390 : * Energie extraite du gaz par convection 3391 : Tfp3 = 'REDU' Tf (GEO.'mtp3') ; 3392 : Tbp3 = tic.'TBP3' ; 3393 : Hconv = KCHT GEO.'$mtp3' SCAL CENTRE tic.'KH3' ; 3394 : Hconv = elno GEO.'$mtp3' Hconv ; 3395 : * Diagmp3 = 'DOMA' (GEO.'$mtp3') 'XXDIAGSI' ; 3396 : Econv3 = 'SOMT' ( GEO.'Diag3f' '*' (Hconv '*' 3397 : (Tfp3 '-' Tbp3))) ; 3398 : 'FINSI' ; 3399 : 3400 : 'SI' TBT.'TECHANP' ; 3401 : * Energie extraite du gaz par convection 3402 : Tfp0 = 'REDU' Tf (GEO.'mtp0') ; 3403 : Tbp0 = tic.'TBP0' ; 3404 : Hconv = KCHT GEO.'$mtp0' SCAL CENTRE tic.'KH0' ; 3405 : Hconv = elno GEO.'$mtp0' Hconv ; 3406 : * Diagmp3 = 'DOMA' (GEO.'$mtp3') 'XXDIAGSI' ; 3407 : Econv0 = 'SOMT' ( GEO.'Diag0f' '*' (Hconv '*' 3408 : (Tfp0 '-' Tbp0))) ; 3409 : 'FINSI' ; 3410 : 3411 : Econv = Econvw '+' Econv1 + Econv2 '+' Econv3 '+' Econv0 ; 3412 : tic.'Econv' = tic.'Econv' 'ET' ('PROG' Econv) ; 3413 : 3414 : Easp = 0.0 ; 3415 : Hasp = 0.0 ; 3416 : Qasp = 0.0 ; 3417 : 3418 : 'SI' TBT.'ASPER ' ; 3419 : Tf = tic.'TF' ; 3420 : Td = tic.'TD' ; 3421 : Hgas = tic.'HGAS' ; 3422 : Diagvtf = doma $vtf 'XXDIAGSI' ; 3423 : Easp = (somt ( Diagvtf * (Hgas * (Td - Tf)))) *Rhomnm*Cpmnm; 3424 : srvp = tic.'SRVP' ; 3425 : Qasp = somt (srvp * Volvtf) ; 3426 : Tfe = noel $vtf tic.'TFNM' ; 3427 : Hasp = somt ((srvp * (tfe + 273.15) * Volvtf) * Cpvap) ; 3428 : 3429 : *****modif 2005 3430 : * température goutte max et min 3431 : Td = tic.'TD' ; 3432 : maTd = 'MAXI' Td ; 3433 : tic.'Tdma' = tic.'Tdma' et (prog maTd) ; 3434 : miTd = 'MINI' Td ; 3435 : tic.'Tdmi' = tic.'Tdmi' et (prog miTd) ; 3436 : ***************** 3437 : 3438 : 3439 : 'FINSI' ; 3440 : 3441 : tic.'Easpe' = tic.'Easpe' et (prog Easp ) ; 3442 : tic.'Haspe' = tic.'Haspe' et (prog Hasp ) ; 3443 : tic.'Qaspe' = tic.'Qaspe' et (prog Qasp ) ; 3444 : 3445 : *--------------------------------------------------------- 3446 : *>9* Prediction de la condensation pour le pas suivant *** 3447 : 3448 : Qc=0. ; Econd=0. ; Hcond=0. ; 3449 : 3450 : 3451 : 3452 : *tc mise en commentaires du sinon et du finsi la ligne au milieu 3453 : * l'etant deja 3454 : *'SINON' ; 3455 : *???? Qcm=0. ; Econdm=0. ; Hcondm=0. ; 3456 : *'FINSI' ; 3457 : 3458 : 'SI' (TBT.'VAPEUR') ; 3459 : tic.'Qcm'= tic.'Qcm' 'ET' ('PROG' Qcm) ; 3460 : 'FINSI' ; 3461 : 3462 : 3463 : *>9-2 Calcul des flux condenses thermique paroi 3464 : 'SI' (TBT.'VAPEUR' 'ET' (TBT.'TPAROIF')) ; 3465 : KHcu = rxt.'ECHAN' ; 3466 : Qcw Fcondw Econdw Hcondw KKC ROVI FHP HT 3467 : = CONDENS rxt (GEO.'$paroif') 'TP' 'TF' 'RVP' KHcu ; 3468 : tic.'KKCW' = KKC ; 3469 : tic.'ROVIW'= ROVI ; 3470 : tic.'FHPW' = FHP ; 3471 : tic.'KHW' = HT ; 3472 : tic.'Fcondw'=Fcondw ; 3473 : *'MESS' 'Calcul des flux condenses ' Qcw 'Kg/s' ; 3474 : Qc=Qc+Qcw ; 3475 : Econd= Econd + Econdw ; 3476 : Hcond= Hcond + Hcondw ; 3477 : 'SINON' ; 3478 : Qcw=0. ; Econdw=0. ; Hcondw=0. ; 3479 : 'FINSI' ; 3480 : 3481 : * Calcul des flux condenses Température imposée 1 3482 : 'SI' (TBT.'VAPEUR' 'ET' TBT.'TTIMP1') ; 3483 : tic.'TPI1'=kcht (GEO.'$mtp1') scal sommet (tic.'TBP1') ; 3484 : KHcu=(rxt.'TIMP1'.'ECHAN') ; 3485 : Qc1 Fcond1 Econd1 Hcond1 KKC ROVI FHP HT 3486 : = CONDENS rxt (GEO.'$mtp1') 'TPI1' 'TF' 'RVP' KHcu ; 3487 : tic.'KKC1' = KKC ; 3488 : tic.'ROVI1'= ROVI ; 3489 : tic.'FHP1' = FHP ; 3490 : tic.'KH1' = HT ; 3491 : tic.'Fcond1'=Fcond1 ; 3492 : 'SI' ( TBT.'IMPR' >EG 2) ; 3493 : 'FINSI' ; 3494 : 'SI' ( TBT.'IMPR' >EG 2) ; 3495 : *'MESS' 'Calcul des flux condenses ' Qc1 'Kg/s' ; 3496 : 'FINSI' ; 3497 : Qc=Qc+Qc1 ; 3498 : Econd= Econd + Econd1 ; 3499 : Hcond= Hcond + Hcond1 ; 3500 : 'SINON' ; 3501 : Qc1=0. ; Econd1=0. ; Hcond1=0. ; 3502 : 'FINSI' ; 3503 : 3504 : * Calcul des flux condenses Température imposée 2 3505 : 'SI' (TBT.'VAPEUR' 'ET' TBT.'TTIMP2') ; 3506 : tic.'TPI2'=kcht (GEO.'$mtp2') scal sommet (tic.'TBP2') ; 3507 : KHcu=(rxt.'TIMP2'.'ECHAN') ; 3508 : Qc2 Fcond2 Econd2 Hcond2 KKC ROVI FHP HT 3509 : = CONDENS rxt (GEO.'$mtp2') 'TPI2' 'TF' 'RVP' KHcu ; 3510 : tic.'KKC2' = KKC ; 3511 : tic.'ROVI2'= ROVI ; 3512 : tic.'FHP2' = FHP ; 3513 : tic.'KH2' = HT ; 3514 : tic.'Fcond2'=Fcond2 ; 3515 : 'SI' ( TBT.'IMPR' >EG 2) ; 3516 : *'MESS' 'Calcul des flux condenses ' Qc2 'Kg/s' ; 3517 : 'FINSI' ; 3518 : Qc=Qc+Qc2 ; 3519 : Econd= Econd + Econd2 ; 3520 : Hcond= Hcond + Hcond2 ; 3521 : 'SINON' ; 3522 : Qc2=0. ; Econd2=0. ; Hcond2=0. ; 3523 : 'FINSI' ; 3524 : 3525 : * Calcul des flux condenses Température imposée 3 3526 : 'SI' (TBT.'VAPEUR' 'ET' TBT.'TTIMP3') ; 3527 : tic.'TPI3'=kcht (GEO.'$mtp3') scal sommet (tic.'TBP3') ; 3528 : KHcu=(rxt.'TIMP3'.'ECHAN') ; 3529 : Qc3 Fcond3 Econd3 Hcond3 KKC ROVI FHP HT 3530 : = CONDENS rxt (GEO.'$mtp3') 'TPI3' 'TF' 'RVP' KHcu ; 3531 : tic.'KKC3' = KKC ; 3532 : tic.'ROVI3'= ROVI ; 3533 : tic.'FHP3' = FHP ; 3534 : tic.'KH3' = HT ; 3535 : tic.'Fcond3'=Fcond3 ; 3536 : 'SI' ( TBT.'IMPR' >EG 2) ; 3537 : *'MESS' 'Calcul des flux condenses ' Qc3 'Kg/s' ; 3538 : 'FINSI' ; 3539 : Qc=Qc+Qc3 ; 3540 : Econd= Econd + Econd3 ; 3541 : Hcond= Hcond + Hcond3 ; 3542 : 'SINON' ; 3543 : Qc3=0. ; Econd3=0. ; Hcond3=0. ; 3544 : 'FINSI' ; 3545 : 3546 : * Calcul des flux condenses 3547 : 'SI' (TBT.'VAPEUR' 'ET' TBT.'TECHANP') ; 3548 : tic.'TPI0'=kcht (GEO.'$mtp0') scal sommet (tic.'TBP0') ; 3549 : KHcu=(rxt.'ECHANP'.'ECHAN') ; 3550 : Qc0 Fcond0 Econd0 Hcond0 KKC ROVI FHP HT 3551 : = CONDENS rxt (GEO.'$mtp0') 'TPI0' 'TF' 'RVP' KHcu ; 3552 : tic.'KKC0' = KKC ; 3553 : tic.'ROVI0'= ROVI ; 3554 : tic.'FHP0' = FHP ; 3555 : tic.'KH0' = HT ; 3556 : tic.'Fcond0'=Fcond0 ; 3557 : 'SI' ( TBT.'IMPR' >EG 2) ; 3558 : *'MESS' 'Calcul des flux condenses ' Qc0 'Kg/s' ; 3559 : 'FINSI' ; 3560 : Qc=Qc+Qc0 ; 3561 : Econd= Econd + Econd0 ; 3562 : Hcond= Hcond + Hcond0 ; 3563 : 'SINON' ; 3564 : Qc0=0. ; Econd0=0. ; Hcond0=0. ; 3565 : 'FINSI' ; 3566 : 3567 : ********************************************************** 3568 : 'SI' TBT.'ASPER ' ; 3569 : asparam rxt ; 3570 : 'FINSI' ; 3571 : 3572 : 3573 : ********************************************************** 3574 : **** Avancement en temps ************************************ 3575 : *mess 'Avancement en temps' ; 3576 : 3577 : * Calcul densité drho/rho pour le terme de flotabilite 3578 : 3579 : Mrho=somt (Diag*Rho); 3580 : Mrhot=somt (Diag*Rhot); 3581 : mess '%$$ Mrho=' Mrho ' Mrhot=' Mrhot ; 3582 : mess '%$$ Mini/Maxi rho=' (mini rho)(maxi rho); 3583 : mess '%$$ Mini/Maxi rhot=' (mini rhot)(maxi rhot); 3584 : 3585 : *Rhop = (-1.) '*' (Rhomn '-' Rho) '*' (Rho **(-1.0)) ; 3586 : Rhop = (-1.) '*' (Rhomn '-' Rhot) '*' (Rhot **(-1.0)) ; 3587 : Rhop = 'KCHT' $vtf 'SCAL' 'SOMMET' Rhop ; 3588 : rmc = 'NOEL' $vtf Rhop ; 3589 : tic.'Rhop' = Rhop ; 3590 : 'SI' DIM3D ; 3591 : rogx = 'KCHT' $vtf 'SCAL' 'CENTRE' 0.0 ; 3592 : rogy = 'KCHT' $vtf 'SCAL' 'CENTRE' 0.0 ; 3593 : rogz = 'KCHT' $vtf 'SCAL' 'CENTRE' (-9.81 '*' rmc); 3594 : rogx = 'NOMC' 'UX' rogx 'NATU' 'DISCRET' ; 3595 : rogy = 'NOMC' 'UY' rogy 'NATU' 'DISCRET' ; 3596 : rogz = 'NOMC' 'UZ' rogz 'NATU' 'DISCRET' ; 3597 : tic.'ROG'=kcht $vtf 'VECT' 'CENTRE' 3598 : (rogx 'ET' rogy 'ET' rogz) ; 3599 : 'SINON' ; 3600 : rogx = 'KCHT' $vtf 'SCAL' 'CENTRE' 0.0 ; 3601 : rogy = 'KCHT' $vtf 'SCAL' 'CENTRE' (-9.81 '*' rmc); 3602 : rogx = 'NOMC' 'UX' rogx 'NATU' 'DISCRET' ; 3603 : rogy = 'NOMC' 'UY' rogy 'NATU' 'DISCRET' ; 3604 : tic.'ROG'=kcht $vtf 'VECT' 'CENTRE' 3605 : (rogx 'ET' rogy) ; 3606 : 'FINSI' ; 3607 : 3608 : * Vitesse max 3609 : un = tic.'UN' ; 3610 : modu = ( 'KOPS' un 'PSCA' un ) '**' 0.5 ; 3611 : mmu = 'MAXI' modu ; 3612 : tic.'LMAXU' = tic.'LMAXU' 'ET' ('PROG' mmu) ; 3613 : 3614 : tic.'Qc' = tic.'Qc' 'ET' ('PROG' QC) ; 3615 : Minj = (Qj '+' Qj2 '+' Qj3 '+' Qlj '+' Qlj2 '+' Qlj3) 3616 : '*' tic.'DT' ; 3617 : MinjT = MinjT '+' Minj ; 3618 : Mcond = Qc '*' tic.'DT' ; 3619 : McondT = McondT '+' Mcond + Mliq ; 3620 : Mrest =(Qj '+' Qj2 '+' Qj3 '-' Qc '+' Qasp) 3621 : '*' tic.'DT' ; 3622 : MrestT = MrestT '+' Mrest ; 3623 : tic.'Minj' = tic.'Minj' 'ET' ('PROG' MinjT) ; 3624 : tic.'Mcond' = tic.'Mcond' 'ET' ('PROG' McondT) ; 3625 : 'SI' TBT.'VAPEUR' ; 3626 : 3627 : 'SI' TBT.'TPAROIF' ; 3628 : Mcondw = Qcw '*' tic.'DT' ; 3629 : McondTw = McondTw '+' Mcondw ; 3630 : tic.'Mcondw'= tic.'Mcondw' 'ET' ('PROG' McondTw) ; 3631 : tic.'Qcw'= tic.'Qcw' 'ET' ('PROG' Qcw) ; 3632 : 'FINSI' ; 3633 : 3634 : 'SI' TBT.'TTIMP1' ; 3635 : Mcond1 = Qc1 '*' tic.'DT' ; 3636 : McondT1 = McondT1 '+' Mcond1 ; 3637 : tic.'Mcond1'= tic.'Mcond1' 'ET' ('PROG' McondT1) ; 3638 : tic.'Qc1'= tic.'Qc1' 'ET' ('PROG' Qc1) ; 3639 : 'FINSI' ; 3640 : 3641 : 'SI' TBT.'TTIMP2' ; 3642 : Mcond2 = Qc2 '*' tic.'DT' ; 3643 : McondT2 = McondT2 '+' Mcond2 ; 3644 : tic.'Mcond2'= tic.'Mcond2' 'ET' ('PROG' McondT2) ; 3645 : tic.'Qc2'= tic.'Qc2' 'ET' ('PROG' Qc2) ; 3646 : 'FINSI' ; 3647 : 3648 : 'SI' TBT.'TTIMP3' ; 3649 : Mcond3 = Qc3 '*' tic.'DT' ; 3650 : McondT3 = McondT3 '+' Mcond3 ; 3651 : tic.'Mcond3'= tic.'Mcond3' 'ET' ('PROG' McondT3) ; 3652 : tic.'Qc3'= tic.'Qc3' 'ET' ('PROG' Qc3) ; 3653 : 'FINSI' ; 3654 : 3655 : 'SI' TBT.'TECHANP' ; 3656 : Mcond0 = Qc0 '*' tic.'DT' ; 3657 : McondT0 = McondT0 '+' Mcond0 ; 3658 : tic.'Mcond0'= tic.'Mcond0' 'ET' ('PROG' McondT0) ; 3659 : tic.'Qc0'= tic.'Qc0' 'ET' ('PROG' Qc0) ; 3660 : 'FINSI' ; 3661 : 3662 : 'FINSI' ; 3663 : 3664 : tic.'Mrest' = tic.'Mrest' 'ET' ('PROG' MrestT) ; 3665 : tic.'Econd' = tic.'Econd' 'ET' ('PROG' Econd) ; 3666 : * tic.'Econdw'= tic.'Econdw' 'ET' ('PROG' Econdw) ; 3667 : * tic.'Econd1'= tic.'Econd1' 'ET' ('PROG' Econd1) ; 3668 : * tic.'Econd2'= tic.'Econd2' 'ET' ('PROG' Econd2) ; 3669 : * tic.'Econd3'= tic.'Econd3' 'ET' ('PROG' Econd3) ; 3670 : tic.'Hcond' = tic.'Hcond' 'ET' ('PROG' Hcond) ; 3671 : * tic.'Hcondw'= tic.'Hcondw' 'ET' ('PROG' Hcondw) ; 3672 : * tic.'Hcond1'= tic.'Hcond1' 'ET' ('PROG' Hcond1) ; 3673 : * tic.'Hcond2'= tic.'Hcond2' 'ET' ('PROG' Hcond2) ; 3674 : * tic.'Hcond3'= tic.'Hcond3' 'ET' ('PROG' Hcond3) ; 3675 : 3676 : 3677 : Si( TBT.'IMPR' >EG 2) ; 3678 : 'MESS' ' ************* BILAN au temps ' Tps ' -> ' 3679 : &bcltps ' DT=' tic.'DT' 3680 : ' ****' ; 3681 : 'SI' TBT.'Tbreche' ; 3682 : 'MESS' '________________________________________________'; 3683 : 'MESS' '-------- Paramètres Injection 1 ---------------' ; 3684 : 'MESS' ' Débit Masse (Kg/s) Qj=' Qj ; 3685 : 'SI' TBT.'VAPEUR' ; 3686 : 'MESS' ' Ksi (taux de vapeur injectée) Ksi=' Ksi ; 3687 : 'FINSI' ; 3688 : 'MESS' ' Vitesse (m/s) guj=' guj ; 3689 : 'MESS' ' Densité à l injection (Kg/m**3) Roj=' Roj ; 3690 : 'MESS' ' R/M a l Injection (J/kg/K) Rgj=' Rgj ; 3691 : 'MESS' ' Température d injection (oC) Tj=' Tj ; 3692 : 'MESS' ' Enthalpie d injection (J/Kg) Hj=' Hj ; 3693 : 'MESS' ' Energie injection (J/Kg) Ej=' Ej ; 3694 : 'MESS' ' Puissance à l injection (Watt) Pj=' (Qj '*'Ej); 3695 : 'MESS' '________________________________________________'; 3696 : 'FINSI' ; 3697 : 'SI' TBT.'Tbreche2' ; 3698 : 'MESS' '-------- Paramètres Injection 2 ---------------' ; 3699 : 'MESS' ' Débit Masse (Kg/s) Qj=' Qj2 ; 3700 : 'SI' TBT.'VAPEUR' ; 3701 : 'MESS' ' Ksi2 (taux de vapeur injectée) Ksi2=' Ksi2 ; 3702 : 'FINSI' ; 3703 : 'MESS' ' Vitesse (m/s) guj=' guj2 ; 3704 : 'MESS' ' Densité à l injection (Kg/m**3) Roj=' Roj2 ; 3705 : 'MESS' ' R/M a l Injection (J/kg/K) Rgj=' Rgj2 ; 3706 : 'MESS' ' Température d injection (oC) Tj=' Tj2 ; 3707 : 'MESS' ' Enthalpie d injection (J/Kg) Hj=' Hj2 ; 3708 : 'MESS' ' Energie injection (J/Kg) Ej=' Ej2 ; 3709 : 'MESS' ' Puissance à l injection (Watt) Pj='(Qj2'*'Ej2); 3710 : 'MESS' '________________________________________________'; 3711 : 'FINSI' ; 3712 : 'SI' TBT.'Tbreche3' ; 3713 : 'MESS' '-------- Paramètres Injection 2 ---------------' ; 3714 : 'MESS' ' Débit Masse (Kg/s) Qj=' Qj3 ; 3715 : 'SI' TBT.'VAPEUR' ; 3716 : 'MESS' ' Ksi3 (taux de vapeur injectée) Ksi3=' Ksi3 ; 3717 : 'FINSI' ; 3718 : 'MESS' ' Vitesse (m/s) guj=' guj3 ; 3719 : 'MESS' ' Densité à l injection (Kg/m**3) Roj=' Roj3 ; 3720 : 'MESS' ' R/M a l Injection (J/kg/K) Rgj=' Rgj3 ; 3721 : 'MESS' ' Température d injection (oC) Tj=' Tj3 ; 3722 : 'MESS' ' Enthalpie d injection (J/Kg) Hj=' Hj3 ; 3723 : 'MESS' ' Energie injection (J/Kg) Ej=' Ej3 ; 3724 : 'MESS' ' Puissance à l injection (Watt) Pj='(Qj3'*'Ej3); 3725 : 'MESS' '________________________________________________'; 3726 : 'FINSI' ; 3727 : 3728 : 'MESS' ' Pression corrigée (Pa) Pt =' Ptnm ; 3729 : 'MESS' '________________________________________________'; 3730 : 'MESS' ' ' ; 3731 : 'SI' TBT.'VAPEUR' ; 3732 : 'MESS' '________________________________________________'; 3733 : 'MESS' '--------Bilan sur la vapeur--------------------' ; 3734 : 'MESS' ' Débit masse condensée (kg/s) QC=' Qc ; 3735 : 'MESS' ' Masse injectée (kg) Minj=' Minj ; 3736 : 'MESS' ' Masse injectée totale (kg) MinjT=' MinjT ; 3737 : 'MESS' ' Masse condensée (kg) Mcond=' Mcond ; 3738 : 'SI' TBT.'TPAROIF' ; 3739 : 'MESS' ' Masse condensée (mur vtp) (kg) Mcondw=' Mcondw ; 3740 : 'MESS' ' Masse condensée totale (mur vtp) (kg) McondTw=' McondTw; 3741 : 'FINSI'; 3742 : 'SI' TBT.'TTIMP1' ; 3743 : 'MESS' ' Masse condensée (mur Timp1) (kg) Mcond1=' Mcond1; 3744 : 'MESS' ' Masse condensée totale (mur Timp1) (kg) McondT1=' McondT1; 3745 : 'FINSI'; 3746 : 'SI' TBT.'TTIMP2' ; 3747 : 'MESS' ' Masse condensée (mur Timp2) (kg) Mcond2=' Mcond2; 3748 : 'MESS' ' Masse condensée totale (mur Timp2) (kg) McondT2=' McondT2; 3749 : 'FINSI'; 3750 : 'SI' TBT.'TTIMP3' ; 3751 : 'MESS' ' Masse condensée (mur Timp3) (kg) Mcond3=' Mcond3; 3752 : 'MESS' ' Masse condensée totale (mur Timp3) (kg) McondT3=' McondT3; 3753 : 'FINSI'; 3754 : 'MESS' ' Masse condensée totale (kg) McondT=' McondT; 3755 : 'MESS' ' Masse vapeur restante (kg) Mrest=' Mrest ; 3756 : 'MESS' ' Masse vapeur restante totale(kg) MrestT=' MrestT; 3757 : 'MESS' ' Puissance condensation (Watt) QC*Lv=' 3758 : (QC '*' Lv) ; 3759 : 'MESS' ' Energie extraite du gaz par condensation (J)' 3760 : 'Econd=' Econd ; 3761 : 'MESS' '________________________________________________'; 3762 : 'MESS' ' ' ; 3763 : 'FINSI' ; 3764 : 'MESS' '________________________________________________'; 3765 : 'MESS' '----- Grandeurs Moyennes ----------------------' ; 3766 : 'MESS' ' Energie extraite du gaz par convection (J) Econv= ' Econv ; 3767 : 'MESS' ' Densité moyenne (kg/m3) Rhomn= ' Rhomn ; 3768 : 'SI' TBT.'VAPEUR' ; 3769 : 'MESS' ' Densité moyenne de vapeur(kg/m3) Rhomvn= ' Rhomvn; 3770 : 'FINSI' ; 3771 : 'SI' TBT.'THE' ; 3772 : 'MESS' ' Densité moyenne de helium(kg/m3) Rmhen= ' Rmhen ; 3773 : 'FINSI' ; 3774 : 'SI' TBT.'TH2' ; 3775 : 'MESS' ' Densité moyenne de Hydrogene(kg/m3) Rmh2n= ' Rmh2n ; 3776 : 'FINSI' ; 3777 : 'SI' TBT.'TO2' ; 3778 : 'MESS' ' Densité moyenne de oxygene (kg/m3) Rmo2n= ' Rmo2n ; 3779 : 'FINSI' ; 3780 : 'SI' TBT.'TN2' ; 3781 : 'MESS' ' Densité moyenne de azote(kg/m3) Rmn2n= ' Rmn2n ; 3782 : 'FINSI' ; 3783 : 'SI' TBT.'TCO' ; 3784 : 'MESS' ' Densité moyenne de CO (kg/m3) Rmcon= ' Rmcon ; 3785 : 'FINSI' ; 3786 : 'SI' TBT.'TCO2' ; 3787 : 'MESS' ' Densité moyenne de CO2 (kg/m3) Rmco2n= ' Rmco2n; 3788 : 'FINSI' ; 3789 : *? TCO2AIR Faut il un test ??? 3790 : 'SI' TBT.'TAIR' ; 3791 : 'MESS' ' Densité moyenne de AIR (kg/m3) Rmairn= ' Rmairn; 3792 : 'FINSI' ; 3793 : 'SI' TBT.'TPAROIF' ; 3794 : 'MESS' ' Température moyenne des murs THERMP (oC) Tpm= ' Tpm ; 3795 : 'SI' ('EGA' ('TYPE' tic.'KHW') 'FLOTTANT') ; 3796 : KHWm=tic.'KHW' ; 3797 : 'SINON' ; 3798 : KHW1 = 'ELNO' GEO.'$paroif' (tic.'KHW') ; 3799 : KHWm = 'SOMT' (GEO.'Diagpf' '*' KHW1) '/' GEO.'Sparoif' ; 3800 : 'FINSI' ; 3801 : 'MESS' ' Coefficient echange convectif moyen (Watt/m**2/°C) KHWm=' 3802 : KHWm ; 3803 : 'FINSI' ; 3804 : 'SI' TBT.'TTIMP1' ; 3805 : 'MESS' ' Température paroi 1 (oC) (TIMP1) TBP1= ' tic.'TBP1'; 3806 : 'SI' TBT.'VAPEUR' ; 3807 : KH11 = 'ELNO' GEO.'$mtp1' (tic.'KH1') ; 3808 : KH1m = 'SOMT' (GEO.'Diag1f' '*' KH11) '/' GEO.'Smtp1' ; 3809 : 'SINON' ; 3810 : KH1m = tic.'KH1' ; 3811 : 'FINSI' ; 3812 : 'MESS' ' Coefficient echange convectif (Watt/m**2/°C) KH1m=' 3813 : KH1m ; 3814 : 'FINSI' ; 3815 : 'SI' TBT.'TTIMP2' ; 3816 : 'MESS' ' Température paroi 2 (oC) (TIMP2) TBP2= ' tic.'TBP2'; 3817 : 'SI' TBT.'VAPEUR' ; 3818 : KH21 = 'ELNO' GEO.'$mtp2' (tic.'KH2') ; 3819 : KH2m = 'SOMT' (GEO.'Diag2f' '*' KH21) '/' GEO.'Smtp2' ; 3820 : 'SINON' ; 3821 : KH2m = tic.'KH2' ; 3822 : 'FINSI' ; 3823 : 'MESS' ' Coefficient echange convectif (Watt/m**2/°C) KH2m=' 3824 : KH2m ; 3825 : 'FINSI' ; 3826 : 'SI' TBT.'TTIMP3' ; 3827 : 'MESS' ' Température paroi 3 (oC) (TIMP3) TBP3= ' tic.'TBP3'; 3828 : 'SI' TBT.'VAPEUR' ; 3829 : KH31 = 'ELNO' GEO.'$mtp3' (tic.'KH3') ; 3830 : KH3m = 'SOMT' (GEO.'Diag3f' '*' KH31) '/' GEO.'Smtp3' ; 3831 : 'SINON' ; 3832 : KH3m = tic.'KH3' ; 3833 : 'FINSI' ; 3834 : 'MESS' ' Coefficient echange convectif (Watt/m**2/°C) KH3m=' 3835 : KH3m ; 3836 : 'FINSI' ; 3837 : 'MESS' ' Température fluide moyenne (oC) Tfm= ' Tfm ; 3838 : 'MESS' ' Pression moyenne (Pa) Pt= ' Pt ; 3839 : 'MESS' ' Gamma moyen Gamm= ' Gamm ; 3840 : 'MESS' ' Rgaz moyen (J/kg/K) Rgpm= ' Rgpm ; 3841 : 'MESS' ' Cv moyen (J/kg/K) Cvm= ' Cvm ; 3842 : 'MESS' ' Cp moyen (J/kg/K) Cpm= ' Cpm ; 3843 : Si(EGA (TYPE Cpvap) 'FLOTTANT') ; 3844 : 'MESS' ' Cp vapeur moyen (J/kg/K) Cpvap= ' Cpvap ; 3845 : Sinon ; 3846 : 'MESS' ' Cp vapeur moyen (J/kg/K) Cpvap= ' 3847 : (mini Cpvap) (maxi Cpvap) ; 3848 : Finsi ; 3849 : 'MESS' '____________________________________________________________'; 3850 : 'MESS' ' ' ; 3851 : 'MESS' '____________________________________________________________'; 3852 : 'MESS' '-------------------- Corrections -------------------------- '; 3853 : 'MESS' ' dTf =' dTf ' K dRho=' dRho ' Kg dRhov = ' dRhov ' Kg' ; 3854 : 'MESS' ' dRhohe=' dRhohe ' Kg dRhoh2=' dRhoh2 ' Kg'; 3855 : 'SI' TBT.'TO2' ; 3856 : 'MESS' ' dRhoo2=' dRhoo2 ' Kg'; 'FINSI' ; 3857 : 'SI' TBT.'TN2' ; 3858 : 'MESS' ' dRhon2=' dRhon2 ' Kg'; 'FINSI' ; 3859 : 'SI' TBT.'TCO' ; 3860 : 'MESS' ' dRhoco=' dRhoco ' Kg'; 'FINSI' ; 3861 : 'SI' TBT.'TCO2' ; 3862 : 'MESS' ' dRhoco2=' dRhoco2 ' Kg'; 'FINSI' ; 3863 : *? TCO2AIR Faut t il un test ?? 3864 : 'SI' TBT.'TAIR' ; 3865 : 'MESS' ' dRhoair=' dRhoair ' Kg'; 'FINSI' ; 3866 : 'MESS' '____________________________________________________________'; 3867 : 'MESS' '___Bilans 0D_________Bilans MultiD _____ Erreur Relative (%)'; 3868 : 3869 : *---- Bilan masse 3870 : M0D = Rhomn '*' Vtotal ; 3871 : MMD = Rhomnt ; 3872 : 'SI' TBT.'VAPEUR' ; 3873 : MV0D = Rhomvn '*' Vtotal ; 3874 : MVMD = Rhomvnt ; 3875 : 'FINSI' ; 3876 : 'SI' TBT.'THE' ; 3877 : MHE0D = Rmhen '*' Vtotal ; 3878 : MHEMD = Rhomhet ; 3879 : 'FINSI' ; 3880 : 'SI' TBT.'TH2' ; 3881 : MH20D = Rmh2n '*' Vtotal ; 3882 : MH2MD = Rhomh2t ; 3883 : 'FINSI' ; 3884 : 'SI' TBT.'TO2' ; 3885 : MO20D = Rmo2n '*' Vtotal ; 3886 : MO2MD = Rhomo2t ; 3887 : 'FINSI' ; 3888 : 'SI' TBT.'TN2' ; 3889 : MN20D = Rmn2n '*' Vtotal ; 3890 : MN2MD = Rhomn2t ; 3891 : 'FINSI' ; 3892 : 'SI' TBT.'TCO' ; 3893 : MCO0D = Rmcon '*' Vtotal ; 3894 : MCOMD = Rhomcot ; 3895 : 'FINSI' ; 3896 : 'SI' TBT.'TCO2' ; 3897 : MCO20D = Rmco2n '*' Vtotal ; 3898 : MCO2MD = Rhomco2t ; 3899 : 'FINSI' ; 3900 : *? TCO2AIR Faut il un test ?? 3901 : 'SI' TBT.'TAIR' ; 3902 : MAIR0D = Rmairn '*' Vtotal ; 3903 : MAIRMD = Rhomairt ; 3904 : 'FINSI' ; 3905 : 3906 : *---- Bilan energie 3907 : E0D = Remn ; 3908 : EMD = Remnt ; 3909 : 'MESS' ' Masse (Kg) ' M0D '_______' MMD '________' 3910 : (100. '*' (M0D '-' MMD) '/' M0D) ; 3911 : 'MESS' ' Energie (J/m3)' E0D '_______' EMD '________' 3912 : (100. '*' (E0D '-' EMD) '/' E0D) ; 3913 : 3914 : 'SI' (TBT.'VAPEUR' 'ET' ('NON' ('EGA' MV0D 0.0))) ; 3915 : 'MESS' ' M Vapeur(Kg) ' MV0D '_______' MVMD '_______' 3916 : (100. '*' (MV0D '-' MVMD) '/' MV0D) ; 3917 : 'FINSI' ; 3918 : 'SI' (TBT.'THE' 'ET' ('NON' ('EGA' MHE0D 0.0))) ; 3919 : 'MESS' ' M Helium (Kg) ' MHE0D '_______' MHEMD '_______' 3920 : (100. '*' (MHE0D '-' MHEMD) '/' MHE0D) ; 3921 : 'FINSI' ; 3922 : 'SI' (TBT.'TH2' 'ET' ('NON' ('EGA' MH20D 0.0))) ; 3923 : 'MESS' ' M Hydrogene (Kg) ' MH20D '_______' MH2MD '_______' 3924 : (100. '*' (MH20D '-' MH2MD) '/' MH20D) ; 3925 : 'FINSI' ; 3926 : 'SI' (TBT.'TO2' 'ET' ('NON' ('EGA' MO20D 0.0))) ; 3927 : 'MESS' ' M Oxygene (Kg) ' MO20D '_______' MO2MD '_______' 3928 : (100. '*' (MO20D '-' MO2MD) '/' MO20D) ; 3929 : 'FINSI' ; 3930 : 'SI' (TBT.'TN2' 'ET' ('NON' ('EGA' MN20D 0.0))) ; 3931 : 'MESS' ' M Azote (Kg) ' MN20D '_______' MN2MD '_______' 3932 : (100. '*' (MN20D '-' MN2MD) '/' MN20D) ; 3933 : 'FINSI' ; 3934 : 'SI' (TBT.'TCO' 'ET' ('NON' ('EGA' MCO0D 0.0))) ; 3935 : 'MESS' ' M CO (Kg) ' MCO0D '_______' MCOMD '_______' 3936 : (100. '*' (MCO0D '-' MCOMD) '/' MCO0D) ; 3937 : 'FINSI' ; 3938 : 'SI' (TBT.'TCO2' 'ET' ('NON' ('EGA' MCO20D 0.0))) ; 3939 : 'MESS' ' M CO2 (Kg) ' MCO20D '_______' MCO2MD '_______' 3940 : (100. '*' (MCO20D '-' MCO2MD) '/' MCO20D) ; 3941 : 'FINSI' ; 3942 : *? TCO2AIR 3943 : 'SI' (TBT.'TAIR' 'ET' ('NON' ('EGA' MAIR0D 0.0))) ; 3944 : 'MESS' ' M AIR (Kg) ' MAIR0D '_______' MAIRMD '_______' 3945 : (100. '*' (MAIR0D '-' MAIRMD) '/' MAIR0D) ; 3946 : 'FINSI' ; 3947 : 'MESS' '___________________________________________________________'; 3948 : 'SI' TBT.'THERMP' ; 3949 : 'MESS' ' Energie Stockée Mur (J) ' Emurn Emurm 3950 : (100. '*' (Emurn '-' Emurm) '/' Emurn) ; 3951 : 'MESS' '___________________________________________________________'; 3952 : 'FINSI' ; 3953 : 'MESS' ' *** Volumes' ; 3954 : 'MESS' ' VTotal=' VTotal ' m3' ; 3955 : Si (TBT.'Tbreche'); 3956 : 'MESS' ' Facteur Geometrique BRECHE (>1) =' GEO.'facgeo' ; 3957 : Finsi ; 3958 : Si (TBT.'Tbreche2'); 3959 : 'MESS' ' Facteur Geometrique BRECHE2 (>1) =' GEO.'facgeo2' ; 3960 : Finsi ; 3961 : Si (TBT.'Tbreche3'); 3962 : 'MESS' ' Facteur Geometrique BRECHE3 (>1) =' GEO.'facgeo3' ; 3963 : Finsi ; 3964 : 'SI' TBT.'THERMP' ; 3965 : Sparoif = GEO.'Sparoif' ; 3966 : 'MESS' ' VTotal paroi =' GEO.'VTotalp' ' m3' ' Sparoif=' Sparoif ' m2'; 3967 : 'FINSI' ; 3968 : 'SI' TBT.'TTIMP1' ; 3969 : 'MESS' ' Surface T imposée (TIMP1) =' GEO.'Smtp1' ' m2' ; 3970 : 'FINSI' ; 3971 : 'SI' TBT.'TTIMP2' ; 3972 : 'MESS' ' Surface T imposée (TIMP2) =' GEO.'Smtp2' ' m2' ; 3973 : 'FINSI' ; 3974 : 'SI' TBT.'TTIMP3' ; 3975 : 'MESS' ' Surface T imposée (TIMP3) =' GEO.'Smtp3' ' m2' ; 3976 : 'FINSI' ; 3977 : 'SI' (TBT.'Tbreche') ; 3978 : 'MESS' ' Sbreche=' Sbreche ' m2' ; 3979 : 'FINSI' ; 3980 : 'MESS' ' **********************************************************'; 3981 : 3982 : 'FINSI' ; 3983 : 3984 : ************************************************************* 3985 : 3986 : 3987 : *----------------------------------------------------------------- 3988 : 3989 : tic.'UNM' = un ; 3990 : tic.'RHONM' = Rho ; 3991 : *---- Modèle K-Epsilon EFM1 -------------------------------- 3992 : 'SI'((EGA TBT.'MODTURB' 'KEPSILON') et (EGA TBT.'FEF' 'EFM1')); 3993 : tic.'KNM' = tic.'KN' ; 3994 : tic.'ENM' = tic.'EN' ; 3995 : 'FINSI' ; 3996 : *---- Modèle K-Epsilon EFM1 -- FIN ------------------------- 3997 : 'SI' TBT.'VAPEUR' ; 3998 : tic.'RVPNM' = Rvp ; 3999 : 'FINSI' ; 4000 : 'SI' TBT.'THE' ; 4001 : tic.'RHEM' = Rhe ; 4002 : 'FINSI' ; 4003 : 'SI' TBT.'TH2' ; 4004 : tic.'RH2M' = Rh2 ; 4005 : 'FINSI' ; 4006 : 'SI' TBT.'TO2' ; 4007 : tic.'RO2M' = Ro2 ; 4008 : 'FINSI' ; 4009 : 'SI' TBT.'TN2' ; 4010 : tic.'RN2M' = Rn2 ; 4011 : 'FINSI' ; 4012 : 'SI' TBT.'TCO2' ; 4013 : tic.'RCO2M' = Rco2 ; 4014 : 'FINSI' ; 4015 : 'SI' TBT.'TAIR' ; 4016 : tic.'RAIRM' = Rair ; 4017 : 'FINSI' ; 4018 : 'SI' TBT.'TCO' ; 4019 : tic.'RCOM' = Rco ; 4020 : 'FINSI' ; 4021 : tic.'TFNM' = Tf ; 4022 : 'SI' TBT.'THERMP' ; 4023 : tic.'TPNM' = Tp ; 4024 : 'FINSI' ; 4025 : 'SI' TBT.'ASPER ' ; 4026 : tic.'VNM' = tic.'VN' ; 4027 : tic.'XDM' = tic.'XD' ; 4028 : tic.'TDM' = tic.'TD' ; 4029 : tic.'DDM' = tic.'DD' ; 4030 : 'FINSI' ; 4031 : *----------------------------------------------------------- 4032 : 4033 : 4034 : 'SI' (EGA TBT.'FEF' 'EFM1') ; 4035 : DTDIFU=rv.'PASDETPS'.'DTDIFU' ; 4036 : DTCONV=rv.'PASDETPS'.'DTCONV' ; 4037 : DTEXPL= (MINI (PROG DTDIFU DTCONV)) * 0.5 ; 4038 : * mess 'DTEXPL=' DTEXPL ; 4039 : DTM1 = tic.'DT' ; 4040 : 4041 : ect=(DTEXPL - tic.'DT') / DTEXPL ; 4042 : Si (ect > 0.1 ) ; 4043 : tic.'DT' = (0.3 * DTEXPL) + (0.7 * tic.'DT') ; 4044 : Finsi ; 4045 : 4046 : Si (ect < 0.1 ) ; 4047 : tic.'DT' = DTEXPL ; 4048 : Finsi ; 4049 : 4050 : Si(rxt.'DT0' < DTEXPL) ; 4051 : tic.'DT' = rxt.'DT0' ; 4052 : Finsi ; 4053 : 4054 : rv.'PASDETPS'.'NUPASDT' = tic.'NUPADT' ; 4055 : rv.'PASDETPS'.'DELTAT-1'= DTM1 ; 4056 : rv.'PASDETPS'.'DELTAT' = 1.e30 ; 4057 : rv.'PASDETPS'.'TPS' = tic.'Tps' ; 4058 : 4059 : 'FINSI' ; 4060 : * Mess ' Fin BclTPS tic.DT=' tic.'DT' ; 4061 : 'FIN' BCLTPS ; 4062 : *>10'FIN' BCLTPS 4063 : 4064 : 'SI' ('ET' (TBT.'DETMAT') ('>' nbit 0)) ; 4065 : 4066 : matvide = 0 ; 4067 : 4068 : rvm = TBT.'RV'.'METHINV' ; 4069 : rvm.'TABRES' = matvide ; 4070 : rvm.'MATASS' = matvide ; 4071 : rvm.'MAPREC' = matvide ; 4072 : rvm = 'ENLE' rvm 'TABRES' ; 4073 : rvm = 'ENLE' rvm 'MATASS' ; 4074 : rvm = 'ENLE' rvm 'MAPREC' ; 4075 : TBT.'RV'.'METHINV' = rvm ; 4076 : 4077 : 'SI' ('EXIS' TBT.'RV' 'PROJ') ; 4078 : rvprm = TBT.'RV'.'PROJ'.'METHINV' ; 4079 : rvprm.'TABRES' = matvide ; 4080 : rvprm.'MATASS' = matvide ; 4081 : rvprm.'MAPREC' = matvide ; 4082 : rvprm = 'ENLE' rvprm 'TABRES' ; 4083 : rvprm = 'ENLE' rvprm 'MATASS' ; 4084 : rvprm = 'ENLE' rvprm 'MAPREC' ; 4085 : TBT.'RV'.'PROJ'.'METHINV' = rvprm ; 4086 : rvpr = TBT.'RV'.'PROJ' ; 4087 : rvpr.'MATP' = matvide ; 4088 : rvpr.'MATC' = matvide ; 4089 : rvpr = 'ENLE' rvpr 'MATP' ; 4090 : rvpr = 'ENLE' rvpr 'MATC' ; 4091 : TBT.'RV'.'PROJ' = rvpr ; 4092 : rv = TBT.'RV' ; 4093 : rv = 'ENLE' rv 'rvpd' ; 4094 : TBT.'RV' = rv ; 4095 : 'FINSI' ; 4096 : 4097 : rtfm = TBT.'RTF'.'METHINV' ; 4098 : rtfm.'TABRES' = matvide ; 4099 : rtfm.'MATASS' = matvide ; 4100 : rtfm.'MAPREC' = matvide ; 4101 : rtfm = 'ENLE' rtfm 'TABRES' ; 4102 : rtfm = 'ENLE' rtfm 'MATASS' ; 4103 : rtfm = 'ENLE' rtfm 'MAPREC' ; 4104 : TBT.'RTF'.'METHINV' = rtfm ; 4105 : 4106 : 'SI' TBT.'THERMP' ; 4107 : rtpm = TBT.'RTP'.'METHINV' ; 4108 : rtpm.'TABRES' = matvide ; 4109 : rtpm.'MATASS' = matvide ; 4110 : rtpm.'MAPREC' = matvide ; 4111 : rtpm = 'ENLE' rtpm 'TABRES' ; 4112 : rtpm = 'ENLE' rtpm 'MATASS' ; 4113 : rtpm = 'ENLE' rtpm 'MAPREC' ; 4114 : TBT.'RTP'.'METHINV' = rtpm ; 4115 : 'FINSI' ; 4116 : 4117 : 'SI' TBT.'VAPEUR' ; 4118 : Rrvapm = TBT.'Rrvap'.'METHINV' ; 4119 : Rrvapm.'TABRES' = matvide ; 4120 : Rrvapm.'MATASS' = matvide ; 4121 : Rrvapm.'MAPREC' = matvide ; 4122 : Rrvapm = 'ENLE' Rrvapm 'TABRES' ; 4123 : Rrvapm = 'ENLE' Rrvapm 'MATASS' ; 4124 : Rrvapm = 'ENLE' Rrvapm 'MAPREC' ; 4125 : TBT.'Rrvap'.'METHINV' = Rrvapm ; 4126 : 'FINSI' ; 4127 : 4128 : 'SI' TBT.'THE' ; 4129 : Rrhem = TBT.'Rrhe'.'METHINV' ; 4130 : Rrhem.'TABRES' = matvide ; 4131 : Rrhem.'MATASS' = matvide ; 4132 : Rrhem.'MAPREC' = matvide ; 4133 : Rrhem = 'ENLE' Rrhem 'TABRES' ; 4134 : Rrhem = 'ENLE' Rrhem 'MATASS' ; 4135 : Rrhem = 'ENLE' Rrhem 'MAPREC' ; 4136 : TBT.'Rrhe'.'METHINV' = Rrhem ; 4137 : 'FINSI' ; 4138 : 4139 : 'SI' TBT.'TH2' ; 4140 : Rrh2m = TBT.'Rrh2'.'METHINV' ; 4141 : Rrh2m.'TABRES' = matvide ; 4142 : Rrh2m.'MATASS' = matvide ; 4143 : Rrh2m.'MAPREC' = matvide ; 4144 : Rrh2m = 'ENLE' Rrh2m 'TABRES' ; 4145 : Rrh2m = 'ENLE' Rrh2m 'MATASS' ; 4146 : Rrh2m = 'ENLE' Rrh2m 'MAPREC' ; 4147 : TBT.'Rrh2'.'METHINV' = Rrh2m ; 4148 : 'FINSI' ; 4149 : 4150 : 'SI' TBT.'TO2' ; 4151 : Rro2m = TBT.'Rro2'.'METHINV' ; 4152 : Rro2m.'TABRES' = matvide ; 4153 : Rro2m.'MATASS' = matvide ; 4154 : Rro2m.'MAPREC' = matvide ; 4155 : Rro2m = 'ENLE' Rro2m 'TABRES' ; 4156 : Rro2m = 'ENLE' Rro2m 'MATASS' ; 4157 : Rro2m = 'ENLE' Rro2m 'MAPREC' ; 4158 : TBT.'Rro2'.'METHINV' = Rro2m ; 4159 : 'FINSI' ; 4160 : 4161 : 'SI' TBT.'TN2' ; 4162 : Rrn2m = TBT.'Rrn2'.'METHINV' ; 4163 : Rrn2m.'TABRES' = matvide ; 4164 : Rrn2m.'MATASS' = matvide ; 4165 : Rrn2m.'MAPREC' = matvide ; 4166 : Rrn2m = 'ENLE' Rrn2m 'TABRES' ; 4167 : Rrn2m = 'ENLE' Rrn2m 'MATASS' ; 4168 : Rrn2m = 'ENLE' Rrn2m 'MAPREC' ; 4169 : TBT.'Rrn2'.'METHINV' = Rrn2m ; 4170 : 'FINSI' ; 4171 : 4172 : 'SI' TBT.'TCO2' ; 4173 : Rrco2m = TBT.'Rrco2'.'METHINV' ; 4174 : Rrco2m.'TABRES' = matvide ; 4175 : Rrco2m.'MATASS' = matvide ; 4176 : Rrco2m.'MAPREC' = matvide ; 4177 : Rrco2m = 'ENLE' Rrco2m 'TABRES' ; 4178 : Rrco2m = 'ENLE' Rrco2m 'MATASS' ; 4179 : Rrco2m = 'ENLE' Rrco2m 'MAPREC' ; 4180 : TBT.'Rrco2'.'METHINV' = Rrco2m ; 4181 : 'FINSI' ; 4182 : 4183 : 'SI' TBT.'TAIR' ; 4184 : Rrairm = TBT.'Rrair'.'METHINV' ; 4185 : Rrairm.'TABRES' = matvide ; 4186 : Rrairm.'MATASS' = matvide ; 4187 : Rrairm.'MAPREC' = matvide ; 4188 : Rrairm = 'ENLE' Rrairm 'TABRES' ; 4189 : Rrairm = 'ENLE' Rrairm 'MATASS' ; 4190 : Rrairm = 'ENLE' Rrairm 'MAPREC' ; 4191 : TBT.'Rrair'.'METHINV' = Rrairm ; 4192 : 'FINSI' ; 4193 : 4194 : 'SI' TBT.'TCO' ; 4195 : Rrcom = TBT.'Rrco'.'METHINV' ; 4196 : Rrcom.'TABRES'= matvide ; 4197 : Rrcom.'MATASS'= matvide ; 4198 : Rrcom.'MAPREC'= matvide ; 4199 : Rrcom = 'ENLE' Rrcom 'TABRES' ; 4200 : Rrcom = 'ENLE' Rrcom 'MATASS' ; 4201 : Rrcom = 'ENLE' Rrcom 'MAPREC' ; 4202 : TBT.'Rrco'.'METHINV' = Rrcom ; 4203 : 'FINSI' ; 4204 : 4205 : 'SI' TBT.'ASPER' ; 4206 : rvnm = TBT.'RVN'.'METHINV' ; 4207 : rvnm.'TABRES'= matvide ; 4208 : rvnm.'MATASS'= matvide ; 4209 : rvnm.'MAPREC'= matvide ; 4210 : rvnm = 'ENLE' rvnm 'TABRES' ; 4211 : rvnm = 'ENLE' rvnm 'MATASS' ; 4212 : rvnm = 'ENLE' rvnm 'MAPREC' ; 4213 : TBT.'RVN'.'METHINV' = rvnm ; 4214 : rxdm = TBT.'RXD'.'METHINV' ; 4215 : rxdm.'TABRES'= matvide ; 4216 : rxdm.'MATASS'= matvide ; 4217 : rxdm.'MAPREC'= matvide ; 4218 : rxdm = 'ENLE' rxdm 'TABRES' ; 4219 : rxdm = 'ENLE' rxdm 'MATASS' ; 4220 : rxdm = 'ENLE' rxdm 'MAPREC' ; 4221 : TBT.'RXD'.'METHINV' = rxdm ; 4222 : rddm = TBT.'RDD'.'METHINV' ; 4223 : rddm.'TABRES'= matvide ; 4224 : rddm.'MATASS'= matvide ; 4225 : rddm.'MAPREC'= matvide ; 4226 : rddm = 'ENLE' rddm 'TABRES' ; 4227 : rddm = 'ENLE' rddm 'MATASS' ; 4228 : rddm = 'ENLE' rddm 'MAPREC' ; 4229 : TBT.'RDD'.'METHINV' = rddm ; 4230 : rtdm = TBT.'RTD'.'METHINV' ; 4231 : rtdm.'TABRES'= matvide ; 4232 : rtdm.'MATASS'= matvide ; 4233 : rtdm.'MAPREC'= matvide ; 4234 : rtdm = 'ENLE' rtdm 'TABRES' ; 4235 : rtdm = 'ENLE' rtdm 'MATASS' ; 4236 : rtdm = 'ENLE' rtdm 'MAPREC' ; 4237 : TBT.'RTD'.'METHINV' = rtdm ; 4238 : 'FINSI' ; 4239 : 4240 : 'FINSI' ; 4241 : 4242 : *----------------------------------------------------------- 4243 : *----- FIN de la boucle principale sur les pas de temps ---- 4244 : *----------------------------------------------------------- 4245 : 4246 : *============================================================== 4247 : Si TBT.'GRAPH' ; 4248 : ******************************************************************* 4249 : ******** Post Graphique de controle 4250 : ******************************************************************* 4251 : 4252 : tic = rxt.'TIC' ; 4253 : $vtf = GEO.'$vtf' ; 4254 : Si (Exist tic 'LTPSC') ; 4255 : ltpsc= tic.'LTPSC' ; 4256 : ltps = extr (tic.'LTPS') (tic.'LTPSC') ; 4257 : Sinon ; 4258 : ltps = tic.'LTPS' ; 4259 : ltpsc= lect 1 pas 1 (dime ltps) ; 4260 : Finsi ; 4261 : 4262 : tab1 = table ; 4263 : tab1.1 = 'MARQ CROI REGU' ; 4264 : tab1.2 = 'MARQ TRIA REGU' ; 4265 : tab1.3 = 'MARQ LOSA REGU' ; 4266 : tab1.'TITRE' = table ; 4267 : 4268 : mvtf= 'DOMA' $vtf maillage ; 4269 : 4270 : un = tic.'UN' ; 4271 : modu= ( 'KOPS' un 'PSCA' un ) ** 0.5 ; mmu= 'MAXI' modu ; 4272 : 'MESS' ' MAXI module de U ********* ' mmu ; 4273 : 4274 : ************************************************************** 4275 : * Débit injection / Débit condensation 4276 : 4277 : evsrcv=EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/s' 4278 : (extr tic.'Qj' ltpsc); 4279 : 4280 : evqc = EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/s' 4281 : (extr tic.'Qc' ltpsc); 4282 : 4283 : 4284 : tab1.'TITRE'. 1 = 'MOT' 'Q_cond Kg/s' ; 4285 : tab1.'TITRE'. 2 = 'MOT' 'Q_inj Kg/s' ; 4286 : 4287 : evsrcv=evsrcv et (EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/s' 4288 : (extr (tic.'Qj2') ltpsc)) ; 4289 : tab1.'TITRE'. 3 = 'MOT' 'Q_inj2 Kg/s' ; 4290 : 4291 : evsrcv=evsrcv et (EVOL TURQ 'MANU' 'Temps en sec' ltps 'Kg/s' 4292 : (extr (tic.'Qj3') ltpsc)) ; 4293 : tab1.'TITRE'. 4 = 'MOT' 'Q_inj3 Kg/s' ; 4294 : 4295 : evtq=evqc et evsrcv; 4296 : 4297 : Si TBT.'ASPER ' ; 4298 : * ASPER dessin 4299 : lqQas = tic.'Qaspe' ; 4300 : evQas = EVOL JAUNE 'MANU' 'Temps en sec' ltps 'Kg/s' 4301 : (extr lqQas ltpsc) ; 4302 : tab1.'TITRE'. 5 = 'MOT' 'Q_aspe Kg/s' ; 4303 : evtq=evtq et evqas ; 4304 : dess evtq 4305 : LEGE TAB1 DATE MIMA TITR ' Débit injection/condensé/aspersion' ; 4306 : Sinon ; 4307 : dess evtq 4308 : LEGE TAB1 DATE MIMA TITR ' Débit injection/condensé' ; 4309 : Finsi ; 4310 : 4311 : ************************************************************** 4312 : ************************************************************** 4313 : * Taux d'injection de la vapeur 4314 : 'SI' (TBT.'VAPEUR') ; 4315 : 4316 : 4317 : *............................................................. 4318 : Si ('EGA'(dime tic.'Tinj') (dime ltpsc)); 4319 : evtinj =EVOL VERT 'MANU' 'Temps en sec' ltps '%' 4320 : (extr tic.'Tinj' ltpsc); 4321 : 4322 : Si ('EGA'(dime tic.'Tinj2') (dime ltpsc)); 4323 : evtinj=evtinj et (EVOL BLEU 'MANU' 'Temps en sec' ltps '%' 4324 : (extr tic.'Tinj2' ltpsc)); 4325 : Finsi ; 4326 : 4327 : Si ('EGA'(dime tic.'Tinj3') (dime ltpsc)); 4328 : evtinj=evtinj et (EVOL BLEU 'MANU' 'Temps en sec' ltps '%' 4329 : (extr tic.'Tinj3' ltpsc)); 4330 : Finsi ; 4331 : 4332 : tab1.'TITRE'. 1 = 'MOT' 'Tinj ' ; 4333 : tab1.'TITRE'. 2 = 'MOT' 'Tinj2' ; 4334 : tab1.'TITRE'. 3 = 'MOT' 'Tinj3' ; 4335 : 4336 : dess evtinj 4337 : LEGE TAB1 DATE MIMA TITR ' Température d injection aux brèches' ; 4338 : Finsi ; 4339 : 4340 : *............................................................. 4341 : 4342 : *............................................................. 4343 : Si ('EGA'(dime tic.'Ksi') (dime ltpsc)); 4344 : evksi =EVOL VERT 'MANU' 'Temps en sec' ltps '%' 4345 : (extr tic.'Ksi' ltpsc); 4346 : 4347 : Si ('EGA'(dime tic.'Ksi2') (dime ltpsc)); 4348 : evksi=evksi et (EVOL BLEU 'MANU' 'Temps en sec' ltps '%' 4349 : (extr tic.'Ksi2' ltpsc)); 4350 : Finsi ; 4351 : 4352 : Si ('EGA'(dime tic.'Ksi3') (dime ltpsc)); 4353 : evksi=evksi et (EVOL BLEU 'MANU' 'Temps en sec' ltps '%' 4354 : (extr tic.'Ksi3' ltpsc)); 4355 : Finsi ; 4356 : 4357 : tab1.'TITRE'. 1 = 'MOT' 'Ksi ' ; 4358 : tab1.'TITRE'. 2 = 'MOT' 'Ksi2' ; 4359 : tab1.'TITRE'. 3 = 'MOT' 'Ksi3' ; 4360 : 4361 : dess evksi 4362 : LEGE TAB1 DATE MIMA TITR ' Taux d injection de la vapeur aux brèches' ; 4363 : Finsi ; 4364 : 4365 : *............................................................. 4366 : 4367 : * Trace des Debits liquide vapeur et total aux breches 4368 : 4369 : tab1.1 = 'MARQ CROI REGU' ; 4370 : tab1.2 = 'MARQ CROI REGU TIRM' ; 4371 : tab1.3 = 'MARQ CROI REGU TIRC' ; 4372 : tab1.4 = 'MARQ TRIA REGU' ; 4373 : tab1.5 = 'MARQ TRIA REGU TIRM' ; 4374 : tab1.6 = 'MARQ TRIA REGU TIRC' ; 4375 : tab1.7 = 'MARQ LOSA REGU' ; 4376 : tab1.8 = 'MARQ LOSA REGU TIRM' ; 4377 : tab1.9 = 'MARQ LOSA REGU TIRC' ; 4378 : 4379 : Si ('EGA'(dime tic.'Qlj') (dime ltpsc)); 4380 : evql =EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/s' 4381 : (extr tic.'Qlj' ltpsc); 4382 : Sinon ; 4383 : evql=(EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/s' 4384 : (prog (dime ltpsc)*0.)); 4385 : Finsi ; 4386 : qt= (extr evql 'ORDO'); 4387 : qt= qt + (extr tic.'Qj' ltpsc); 4388 : evqlj=(EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/s' qt); 4389 : evqlj=evqlj et evql; 4390 : evqlj=evqlj et (EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/s' 4391 : (extr tic.'Qj' ltpsc)); 4392 : 4393 : tab1.'TITRE'. 1 = 'MOT' 'Qjtot' ; 4394 : tab1.'TITRE'. 2 = 'MOT' 'Qlj ' ; 4395 : tab1.'TITRE'. 3 = 'MOT' 'Qj '; 4396 : 4397 : Si ('EGA'(dime tic.'Qlj2') (dime ltpsc)); 4398 : evql=(EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/s' 4399 : (extr tic.'Qlj2' ltpsc)); 4400 : Sinon ; 4401 : evql=(EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/s' 4402 : (prog (dime ltpsc)*0.)); 4403 : Finsi ; 4404 : qt= (extr evql 'ORDO'); 4405 : qt= qt + (extr tic.'Qj2' ltpsc); 4406 : evqlj=evqlj et (EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/s' qt); 4407 : evqlj=evqlj et evql ; 4408 : evqlj=evqlj et (EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/s' 4409 : (extr tic.'Qj2' ltpsc)); 4410 : tab1.'TITRE'. 4 = 'MOT' 'Qj2tot'; 4411 : tab1.'TITRE'. 5 = 'MOT' 'Qlj2'; 4412 : tab1.'TITRE'. 6 = 'MOT' 'Qj2'; 4413 : 4414 : Si ('EGA'(dime tic.'Qlj3') (dime ltpsc)); 4415 : evql=(EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/s' 4416 : (extr tic.'Qlj3' ltpsc)); 4417 : Sinon ; 4418 : evql=(EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/s' 4419 : (prog (dime ltpsc)*0.)); 4420 : Finsi ; 4421 : qt= (extr evql 'ORDO'); 4422 : qt= qt + (extr tic.'Qj3' ltpsc); 4423 : evqlj=evqlj et (EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/s' qt); 4424 : evqlj=evqlj et evql ; 4425 : evqlj=evqlj et (EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/s' 4426 : (extr tic.'Qj3' ltpsc)); 4427 : tab1.'TITRE'. 7 = 'MOT' 'Qj3tot' ; 4428 : tab1.'TITRE'. 8 = 'MOT' 'Qlj3' ; 4429 : tab1.'TITRE'. 9 = 'MOT' 'Qj3' ; 4430 : 4431 : dess evqlj 4432 : LEGE TAB1 DATE MIMA 4433 : TITR ' Debits injection liquide/vapeur aux brèches' ; 4434 : 4435 : tab1.1 = 'MARQ CROI REGU' ; 4436 : tab1.2 = 'MARQ TRIA REGU' ; 4437 : tab1.3 = 'MARQ LOSA REGU' ; 4438 : 4439 : *............................................................. 4440 : 4441 : 'FINSI' ; 4442 : ************************************************************** 4443 : ************************************************************** 4444 : * Débits condensations Total et détails 4445 : Si (TBT.'VAPEUR'); 4446 : 4447 : tab2 = table ; 4448 : tab2.1 = 'MARQ CROI REGU' ; 4449 : tab2.2 = 'MARQ TRIA REGU' ; 4450 : tab2.3 = 'MARQ LOSA REGU' ; 4451 : tab2.4 = 'MARQ CARR REGU' ; 4452 : tab2.5 = 'MARQ TRIB REGU' ; 4453 : tab2.'TITRE' = table ; 4454 : 4455 : evsrcv=EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/s' 4456 : (extr tic.'Qc' ltpsc); 4457 : evtq = evsrcv ; 4458 : tab2.'TITRE'. 1 = 'MOT' 'Q_cond Kg/s' ; 4459 : 4460 : Numero = 1 ; 4461 : Si('EXIST' tic 'Qc1') ; 4462 : Numero = Numero + 1 ; 4463 : evqc1 = EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/s' 4464 : (extr tic.'Qc1' ltpsc); 4465 : evtq=evtq et evqc1 ; 4466 : tab2.'TITRE'. Numero = 'MOT' 'Qc1 Kg/s' ; 4467 : Finsi ; 4468 : 4469 : Si('EXIST' tic 'Qc2') ; 4470 : Numero = Numero + 1 ; 4471 : evqc2 = EVOL JAUNE 'MANU' 'Temps en sec' ltps 'Kg/s' 4472 : (extr tic.'Qc2' ltpsc); 4473 : evtq=evtq et evqc2 ; 4474 : tab2.'TITRE'. Numero = 'MOT' 'Qc2 Kg/s' ; 4475 : Finsi ; 4476 : 4477 : Si('EXIST' tic 'Qc3') ; 4478 : Numero = Numero + 1 ; 4479 : evqc3 = EVOL ROSE 'MANU' 'Temps en sec' ltps 'Kg/s' 4480 : (extr tic.'Qc3' ltpsc); 4481 : evtq=evtq et evqc3 ; 4482 : tab2.'TITRE'. Numero = 'MOT' 'Qc3 Kg/s' ; 4483 : Finsi ; 4484 : 4485 : Si('EXIST' tic 'Qcw') ; 4486 : Numero = Numero + 1 ; 4487 : evqcw = EVOL TURQ 'MANU' 'Temps en sec' ltps 'Kg/s' 4488 : (extr tic.'Qcw' ltpsc); 4489 : evtq=evtq et evqcw ; 4490 : tab2.'TITRE'. Numero = 'MOT' 'Qcw Kg/s' ; 4491 : Finsi ; 4492 : 4493 : Si('EXIST' tic 'Qc0') ; 4494 : Numero = Numero + 1 ; 4495 : evqc0 = EVOL TURQ 'MANU' 'Temps en sec' ltps 'Kg/s' 4496 : (extr tic.'Qc0' ltpsc); 4497 : evtq=evtq et evqc0 ; 4498 : tab2.'TITRE'. Numero = 'MOT' 'Qc0 Kg/s' ; 4499 : Finsi ; 4500 : 4501 : Si('EXIST' tic 'Qcm') ; 4502 : Numero = Numero + 1 ; 4503 : evqcm = EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/s' 4504 : (extr tic.'Qcm' ltpsc); 4505 : evtq=evtq et evqcm ; 4506 : tab2.'TITRE'. Numero = 'MOT' 'Qcm Kg/s' ; 4507 : Finsi ; 4508 : 4509 : dess evtq 4510 : 4511 : LEGE TAB2 DATE MIMA TITR ' Débits condensés détails' ; 4512 : Finsi ; 4513 : ************************************************************** 4514 : 4515 : ******************************************************************* 4516 : 4517 : evsrhv=EVOL VERT 'MANU' 'Temps en sec' ltps 'Htot' 4518 : (extr tic.'Hj' ltpsc); 4519 : lqe = tic.'Remn' ; 4520 : evqe = EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Hj-Em ' 4521 : (extr lqe ltpsc) ; 4522 : tab1.'TITRE'. 2 = 'MOT' 'Hj J/Kg' ; 4523 : tab1.'TITRE'. 1 = 'MOT' 'Em J/kg' ; 4524 : Si('EXIST' tic 'Hj2') ; 4525 : evsrhv= evsrhv et (EVOL BLEU 'MANU' 'Temps en sec' ltps 'HVap' 4526 : (extr (tic.'Hj2') ltpsc) ) ; 4527 : tab1.'TITRE'. 3 = 'MOT' 'Hj2 J/kg' ; 4528 : Finsi ; 4529 : Si('EXIST' tic 'Hj3') ; 4530 : evsrhv= evsrhv et (EVOL TURQ 'MANU' 'Temps en sec' ltps 'HVap' 4531 : (extr (tic.'Hj3') ltpsc) ) ; 4532 : tab1.'TITRE'. 4 = 'MOT' 'Hj3 J/kg' ; 4533 : Finsi ; 4534 : 4535 : dess (evqe et evsrhv) 4536 : LEGE TAB1 DATE MIMA TITR 4537 : ' Enthalpie(s) à la (aux) brèche(s) / énergie moyenne' ; 4538 : ******************************************************************* 4539 : 4540 : 4541 : lqEcd = tic.'Econd' ; 4542 : evEcd = EVOL VERT 'MANU' 'Temps en sec' ltps 'Ecd' (extr lqEcd ltpsc); 4543 : 4544 : lqEcv = tic.'Econv' ; 4545 : evEcv = EVOL BLEU 'MANU' 'Temps en sec' ltps 'Ecv' (extr lqEcv ltpsc); 4546 : 4547 : 4548 : 'SI' TBT.'VAPEUR' ; 4549 : tab1.'TITRE'. 1 = 'MOT' 'Ecv Joule' ; 4550 : dess evEcd 4551 : LEGE TAB1 DATE MIMA TITR ' Energie de condensation ' ; 4552 : dess evEcv 4553 : LEGE TAB1 DATE MIMA TITR ' Energie echange convectif'; 4554 : 'FINSI' ; 4555 : 4556 : lqMinj = extr tic.'Minj' ltpsc ; 4557 : evMinjv = EVOL BLEU 'MANU' 'Temps en sec' ltps 'Minj-Mcond' lqMinj ; 4558 : lqMcond = extr tic.'Mcond' ltpsc ; 4559 : evMcond = EVOL VERT 'MANU' 'Temps en sec' ltps 'Minj-Mcond' lqMcond ; 4560 : lqMrest = extr tic.'Mrest' ltpsc ; 4561 : evMrest = EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Minj-Mcond' lqMrest ; 4562 : Si (EGA (dime tic.'Mliqpuis') (dime ltpsc)); 4563 : lqMpuis = extr tic.'Mliqpuis' ltpsc ; 4564 : Sinon; 4565 : lqMpuis=prog (dime ltpsc)*0.; 4566 : Finsi; 4567 : evMpuis = EVOL JAUNE 'MANU' 'Temps en sec' ltps 'Minj-Mcond' lqMpuis ; 4568 : 4569 : ******************************************************************** 4570 : tab1.'TITRE'. 1 = 'MOT' 'M_inj Kg' ; 4571 : tab1.'TITRE'. 2 = 'MOT' 'M_cond Kg' ; 4572 : tab1.'TITRE'. 3 = 'MOT' 'M_rest Kg' ; 4573 : tab1.'TITRE'. 4 = 'MOT' 'M_puis Kg' ; 4574 : dess (evMinjv et evMcond et evMrest et evMpuis) 4575 : LEGE TAB1 DATE MIMA 4576 : TITR ' Masses injectées/ condensées/ restantes/ Puisard' ; 4577 : ******************************************************************** 4578 : 4579 : Si ( 'EXIST' tic 'guj' ) ; 4580 : ******************************************************************* 4581 : * Vitesse max / Vitesse injection brèche 1 4582 : evmaxu= EVOL VERT 'MANU' 'Temps en sec' ltps 'M / s ' 4583 : (extr tic.'LMAXU' ltpsc); 4584 : evgj = EVOL ROSE 'MANU' 'Temps en sec' ltps 'M / s ' 4585 : (extr tic.'guj' ltpsc) ; 4586 : tab1.'TITRE'. 1 = 'MOT' 'Max u m/s' ; 4587 : tab1.'TITRE'. 2 = 'MOT' 'U inj m/s' ; 4588 : dess (evmaxu et evgj) 4589 : LEGE TAB1 DATE MIMA TITR ' Vitesse max / Vitesse inject brèche 1 ' ; 4590 : ******************************************************************* 4591 : Finsi ; 4592 : 4593 : Si ( 'EXIST' tic 'guj2' ) ; 4594 : ******************************************************************* 4595 : * Vitesse max / Vitesse injection brèche 2 4596 : evgj2 = EVOL ROSE 'MANU' 'Temps en sec' ltps 'M / s ' 4597 : (extr tic.'guj2' ltpsc); 4598 : tab1.'TITRE'. 1 = 'MOT' 'Max u m/s' ; 4599 : tab1.'TITRE'. 2 = 'MOT' 'U inj m/s' ; 4600 : dess (evmaxu et evgj2) 4601 : LEGE TAB1 DATE MIMA TITR ' Vitesse max / Vitesse inject brèche 2 ' ; 4602 : ******************************************************************* 4603 : Finsi ; 4604 : 4605 : Si ( 'EXIST' tic 'guj3' ) ; 4606 : ******************************************************************* 4607 : * Vitesse max / Vitesse injection brèche 3 4608 : evgj3 = EVOL ROSE 'MANU' 'Temps en sec' ltps 'M / s ' 4609 : (extr tic.'guj3' ltpsc); 4610 : tab1.'TITRE'. 1 = 'MOT' 'Max u m/s' ; 4611 : tab1.'TITRE'. 2 = 'MOT' 'U inj m/s' ; 4612 : dess (evmaxu et evgj3) 4613 : LEGE TAB1 DATE MIMA TITR ' Vitesse max / Vitesse inject brèche 3 ' ; 4614 : ******************************************************************* 4615 : Finsi ; 4616 : 4617 : ******************************************************************* 4618 : * Densités 4619 : evrhom= EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4620 : (extr tic.'Rhom' ltpsc) ; 4621 : tab1.'TITRE'. 1 = 'MOT' 'Rhom Kg/m3' ; 4622 : evdens = evrhom ; 4623 : 4624 : ID = 1 ; 4625 : 'SI' TBT.'VAPEUR' ; 4626 : ID = ID + 1 ; 4627 : evrhomv= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4628 : (extr tic.'Rhomv' ltpsc) ; 4629 : tab1.'TITRE'. ID = 'MOT' 'Rhom_vap Kg/m3' ; 4630 : evdens = evdens et evrhomv ; 4631 : evrhomvg= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4632 : (extr tic.'Rhomvg' ltpsc) ; 4633 : ID = ID + 1 ; 4634 : tab1.'TITRE'. ID = 'MOT' 'Rhom_vg Kg/m3' ; 4635 : evdens = evdens et evrhomvg ; 4636 : 'FINSI' ; 4637 : 4638 : 'SI' TBT.'THE'; 4639 : ID = ID + 1 ; 4640 : tab1.'TITRE'. ID = 'MOT' 'Rhom_He Kg/m3' ; 4641 : evrhomhe=EVOL JAUNE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4642 : (extr tic.'Rhomhe' ltpsc) ; 4643 : evdens = evdens et evrhomhe ; 4644 : 'FINSI' ; 4645 : 4646 : 'SI' TBT.'TH2' ; 4647 : ID = ID + 1 ; 4648 : tab1.'TITRE'. ID = 'MOT' 'Rhom_H2 Kg/m3' ; 4649 : evrhomh2=EVOL BLEU 'MANU' 'Temps en sec' ltps 'Kg/m3' 4650 : (extr tic.'Rhomh2' ltpsc) ; 4651 : evdens = evdens et evrhomh2 ; 4652 : 'FINSI' ; 4653 : 4654 : 'SI' TBT.'TO2'; 4655 : ID = ID + 1 ; 4656 : tab1.'TITRE'. ID = 'MOT' 'Rhom_O2 Kg/m3' ; 4657 : evrhomo2=EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4658 : (extr tic.'Rhomo2' ltpsc) ; 4659 : evdens = evdens et evrhomo2 ; 4660 : 'FINSI' ; 4661 : 4662 : 'SI' TBT.'TN2'; 4663 : ID = ID + 1 ; 4664 : tab1.'TITRE'. ID = 'MOT' 'Rhom_N2 Kg/m3' ; 4665 : evrhomn2=EVOL ROSE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4666 : (extr tic.'Rhomn2' ltpsc); 4667 : evdens = evdens et evrhomn2 ; 4668 : 'FINSI' ; 4669 : 4670 : 'SI' TBT.'TCO' ; 4671 : ID = ID + 1 ; 4672 : tab1.'TITRE'. ID = 'MOT' 'Rhom_CO Kg/m3' ; 4673 : evrhomco=EVOL JAUN 'MANU' 'Temps en sec' ltps 'Kg/m3' 4674 : (extr tic.'Rhomco' ltpsc); 4675 : evdens = evdens et evrhomco ; 4676 : 'FINSI' ; 4677 : 4678 : 'SI' TBT.'TCO2' ; 4679 : ID = ID + 1 ; 4680 : tab1.'TITRE'. ID = 'MOT' 'Rhom_CO2 Kg/m3' ; 4681 : evrco2 =EVOL ROSE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4682 : (extr tic.'Rhomco2' ltpsc); 4683 : evdens = evdens et evrco2 ; 4684 : 'FINSI' ; 4685 : 4686 : *? TCO2AIR 4687 : 'SI' TBT.'TAIR' ; 4688 : ID = ID + 1 ; 4689 : tab1.'TITRE'. ID = 'MOT' 'Rhom_AIR Kg/m3' ; 4690 : evrair =EVOL JAUNE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4691 : (extr tic.'Rhomair' ltpsc); 4692 : evdens = evdens et evrair ; 4693 : 'FINSI' ; 4694 : 4695 : dess evdens 4696 : LEGE TAB1 DATE MIMA TITR ' Densités mélange / vapeur / incondensables'; 4697 : ******************************************************************* 4698 : 4699 : ******************************************************************* 4700 : * Densités et correction 0D 4701 : evrhom= EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4702 : (extr tic.'Rhom' ltpsc); 4703 : tab1.'TITRE'. 1 = 'MOT' 'Rhom Kg/m3' ; 4704 : evrho = evrhom ; 4705 : 4706 : evdrho = EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4707 : (extr tic.'drho' ltpsc); 4708 : tab1.'TITRE'. 2 = 'MOT' 'dRho Kg/m3' ; 4709 : evrho = evrho et evdrho ; 4710 : 4711 : dess evrho 4712 : LEGE TAB1 DATE MIMA TITR ' Densité mélange et correction ' ; 4713 : 4714 : 'SI' TBT.'VAPEUR' ; 4715 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4716 : (extr tic.'Rhomv' ltpsc); 4717 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_VAP Kg/m3' ; 4718 : evrhov = evrhomv; 4719 : 4720 : evrhomvg = EVOL TURQ 'MANU' 'Temps en sec' ltps 'Kg/m3' 4721 : (extr tic.'Rhomvg' ltpsc); 4722 : tab1.'TITRE'. 2 = 'MOT' 'Rhom_H2O Kg/m3' ; 4723 : evrhov = evrhov et evrhomvg; 4724 : 4725 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4726 : (extr tic.'Mdrvap' ltpsc); 4727 : tab1.'TITRE'. 3 = 'MOT' 'Maxi dRho_H2O Kg/m3' ; 4728 : evrhov = evrhov et evdrhov1 ; 4729 : 4730 : dess evrhov 4731 : LEGE TAB1 DATE MIMA TITR ' Densité vapeur et correction ' ; 4732 : 'FINSI' ; 4733 : 4734 : 'SI' TBT.'THE'; 4735 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4736 : (extr tic.'Rhomhe' ltpsc); 4737 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_He Kg/m3' ; 4738 : evrhov = evrhomv; 4739 : 4740 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4741 : (extr tic.'Mdrhe' ltpsc); 4742 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dRho_He Kg/m3' ; 4743 : evrhov = evrhov et evdrhov1 ; 4744 : 4745 : dess evrhov 4746 : LEGE TAB1 DATE MIMA TITR ' Densité Helium et correction ' ; 4747 : 'FINSI' ; 4748 : 4749 : 'SI' TBT.'TH2' ; 4750 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4751 : (extr tic.'Rhomh2' ltpsc); 4752 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_H2 Kg/m3' ; 4753 : evrhov = evrhomv; 4754 : 4755 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4756 : (extr tic.'Mdrh2' ltpsc); 4757 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dRho_H2 Kg/m3' ; 4758 : evrhov = evrhov et evdrhov1 ; 4759 : 4760 : dess evrhov 4761 : LEGE TAB1 DATE MIMA TITR ' Densité Hydrogène et correction ' ; 4762 : 'FINSI' ; 4763 : 4764 : 'SI' TBT.'TO2'; 4765 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4766 : (extr tic.'Rhomo2' ltpsc); 4767 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_o2 Kg/m3' ; 4768 : evrhov = evrhomv; 4769 : 4770 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4771 : (extr tic.'Mdro2' ltpsc); 4772 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dRho_o2 Kg/m3' ; 4773 : evrhov = evrhov et evdrhov1 ; 4774 : 4775 : dess evrhov 4776 : LEGE TAB1 DATE MIMA TITR ' Densité Oxygène et correction ' ; 4777 : 'FINSI' ; 4778 : 4779 : 4780 : 'SI' TBT.'TN2'; 4781 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4782 : (extr tic.'Rhomn2' ltpsc); 4783 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_N2 Kg/m3' ; 4784 : evrhov = evrhomv; 4785 : 4786 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4787 : (extr tic.'Mdrn2' ltpsc); 4788 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dRho_N2 Kg/m3' ; 4789 : evrhov = evrhov et evdrhov1 ; 4790 : 4791 : dess evrhov 4792 : LEGE TAB1 DATE MIMA TITR ' Densité Azote et correction ' ; 4793 : 'FINSI' ; 4794 : 4795 : 'SI' TBT.'TCO'; 4796 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4797 : (extr tic.'Rhomco' ltpsc); 4798 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_Co Kg/m3' ; 4799 : evrhov = evrhomv; 4800 : 4801 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4802 : (extr tic.'Mdrco' ltpsc); 4803 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dRho_Co Kg/m3' ; 4804 : evrhov = evrhov et evdrhov1 ; 4805 : 4806 : dess evrhov 4807 : LEGE TAB1 DATE MIMA TITR ' Densité Co et correction ' ; 4808 : 'FINSI' ; 4809 : 4810 : *? TCO2AIR Pas fait !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 4811 : 'SI' TBT.'TCO2'; 4812 : evrhomv = EVOL VERT 'MANU' 'Temps en sec' ltps 'Kg/m3' 4813 : (extr tic.'Rhomco2' ltpsc); 4814 : tab1.'TITRE'. 1 = 'MOT' 'Rhom_Co2 Kg/m3' ; 4815 : evrhov = evrhomv; 4816 : 4817 : evdrhov1= EVOL ROUGE 'MANU' 'Temps en sec' ltps 'Kg/m3' 4818 : (extr tic.'Mdrco2' ltpsc); 4819 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dRho_Co2 Kg/m3' ; 4820 : evrhov = evrhov et evdrhov1 ; 4821 : 4822 : dess evrhov 4823 : LEGE TAB1 DATE MIMA TITR ' Densité Co2 et correction ' ; 4824 : 'FINSI' ; 4825 : 4826 : ******************************************************************* 4827 : 4828 : ******************************************************************* 4829 : * Pression moyenne 4830 : evPT=EVOL BLEU 'MANU' 'Temps en sec' ltps 'Pt' 4831 : (extr tic.'PT' ltpsc) ; 4832 : tab1.'TITRE'. 1 = 'MOT' 'Pmoy Pa' ; 4833 : dess evPT 4834 : LEGE TAB1 DATE MIMA TITR ' Pression moyenne '; 4835 : ******************************************************************* 4836 : 4837 : ******************************************************************* 4838 : * Températures Fluides paroi etc 4839 : evtpf = EVOL ROUGE 'MANU' 'Temps en sec' ltps 'oC ' 4840 : (extr tic.'Tfm' ltpsc); 4841 : tab1.'TITRE'. 1 = 'MOT' 'Tfm oC' ; 4842 : evtpf1= EVOL BLEU 'MANU' 'Temps en sec' ltps 'oC' 4843 : (extr tic.'MdTf' ltpsc); 4844 : tab1.'TITRE'. 2 = 'MOT' 'Maxi dTf oC' ; 4845 : evtpf = evtpf et evtpf1 ; 4846 : evtpf2= EVOL ROSE 'MANU' 'Temps en sec' ltps 'oC' 4847 : (extr tic.'mdTf' ltpsc); 4848 : tab1.'TITRE'. 3 = 'MOT' 'Mini dTf oC' ; 4849 : evtpf = evtpf et evtpf2 ; 4850 : 4851 : ID = 3 ; 4852 : 'SI' TBT.'THERMP' ; 4853 : ID = ID + 1 ; 4854 : tab1.'TITRE'. ID= 'MOT' 'Tfp oC' ; 4855 : evTfp = EVOL VERT 'MANU' 'Temps en sec' ltps 'oC ' 4856 : (extr tic.'Tpm' ltpsc); 4857 : evtpf = evtpf et evTfp ; 4858 : 'FINSI' ; 4859 : 'SI' TBT.'TTIMP1' ; 4860 : ID = ID + 1 ; 4861 : tab1.'TITRE'. ID= 'MOT' 'Tfpc1 oC' ; 4862 : evTfpc1= EVOL TURQ 'MANU' 'Temps en sec' ltps 'oC ' 4863 : (extr tic.'Ltbp1' ltpsc); 4864 : evtpf = evtpf et evTfpc1 ; 4865 : 'FINSI' ; 4866 : 'SI' TBT.'TTIMP2' ; 4867 : ID = ID + 1 ; 4868 : tab1.'TITRE'. ID= 'MOT' 'Tfpc2 oC' ; 4869 : evTfpc2= EVOL TURQ 'MANU' 'Temps en sec' ltps 'oC ' 4870 : (extr tic.'Ltbp2' ltpsc); 4871 : evtpf = evtpf et evTfpc2 ; 4872 : 'FINSI' ; 4873 : 'SI' TBT.'TTIMP3' ; 4874 : ID = ID + 1 ; 4875 : tab1.'TITRE'. ID= 'MOT' 'Tfpc3 oC' ; 4876 : evTfpc3= EVOL TURQ 'MANU' 'Temps en sec' ltps 'oC ' 4877 : (extr tic.'Ltbp3' ltpsc); 4878 : evtpf = evtpf et evTfpc3 ; 4879 : 'FINSI' ; 4880 : dess evtpf 4881 : LEGE TAB1 DATE MIMA TITR ' Températures moyennes Fluide / Parois'; 4882 : 4883 : 4884 : Finsi ; 4885 : *============================================================== 4886 : 'FINPROC' ; 4887 : 4888 : 4889 : 4890 : 4891 : 4892 : 4893 : 4894 :
© Cast3M 2003 - All rights reserved.
Disclaimer