1 : * BROCHE PROCEDUR MAGN 12/03/20 21:15:00 7322 2 : 'DEBPROC' broche KAS*'MOT' nab*'ENTIER' rxt*'TABLE'; 3 : *mess '=================== DEBUT BROCHE ==================='; 4 : 5 : GEO = rxt.'GEO' ; 6 : nub = nab ; 7 : ibr=index (rxt.'Sorties') ; 8 : Briches =rxt.'Sorties'.(ibr.nub) ; 9 : rxtsce =Briches.'scenario' ; 10 : GEO =Briches ; 11 : CTinj =chai 'Tinj' ; 12 : CKsi =chai 'Ksi' ; 13 : CQlj =chai 'Qlj' ; 14 : CHlj =chai 'Hlj' ; 15 : 16 : TTOUT=FAUX ; 17 : TTINJ=FAUX ; 18 : THINJ=FAUX ; 19 : TTINJ=exist (rxtsce ) 'tinj' ; 20 : TTOUT=exist (rxtsce ) 'Pext' ; 21 : Briches.'TTOUT'=TTOUT ; 22 : 23 : 24 : Si ((NON TTOUT) et (NON TTINJ) ) ; 25 : Mess '*********************************************************' ; 26 : Mess ' ERREUR ERREUR ERREUR ERREUR ERREUR ERREUR ' ; 27 : Mess ' ' ; 28 : Mess 'il faut imposer au moins tinj ou pression ! ' ; 29 : Mess '*********************************************************' ; 30 : ERREUR 21 ; 31 : Finsi ; 32 : 33 : * Modifications ANTS du 06/01/2009 34 : * Evite d'avoir un debit d'injection nul 35 : 36 : 'SI' (('NON' ('EXIS' rxtsce Qmodif)) et (TTINJ ou THINJ)); 37 : 'MESS' '=============================================================' ; 38 : 'MESS' '= WARNING !! =' ; 39 : 'MESS' '= LES DEBITS D INJECTION NULS ONT ETE IMPOSES A 1E-10 =' ; 40 : 'MESS' '=============================================================' ; 41 : lstqxxx = 'MOTS' 'qeau' 'qliq' 'qair' 'qhe' 42 : 'qh2' 'qo2' 'qn2' 'qco' 'qco2' ; 43 : lmots1 = 'MOTS' ; 44 : 'REPE' bqxxx ('DIME' lstqxxx) ; 45 : lmots = 'EXTR' lstqxxx &bqxxx ; 46 : 'SI' ('EXIS' rxtsce lmots) ; 47 : lmots1 = lmots1 'ET' ('MOTS' lmots) ; 48 : 'FINSI' ; 49 : 'FIN' bqxxx ; 50 : mots1 = 'EXTR' lmots1 1 ; 51 : 'REPE' bvxxx ('DIME' rxtsce.mots1) ; 52 : indic = FAUX ; 53 : 'REPETER' bq1xxx ('DIME' lmots1) ; 54 : motsxx = 'EXTR' lmots1 &bq1xxx ; 55 : valxxx = 'EXTR' rxtsce.motsxx &bvxxx ; 56 : 'SI' ('>' valxxx 1.E-10) ; 57 : indic = VRAI ; 58 : 'QUIT' bq1xxx ; 59 : 'FINSI' ; 60 : 'FIN' bq1xxx ; 61 : 'SI' ('NON' indic) ; 62 : 'REMP' rxtsce.mots1 &bvxxx 1E-10 ; 63 : 'FINSI' ; 64 : 'FIN' bvxxx ; 65 : rxtsce.Qmodif = VRAI ; 66 : 'FINSI' ; 67 : * Fin modifications ANTS du 06/01/2009 68 : 69 : ************************************************************************ 70 : ********* DEBUT KAS1 *************************************************** 71 : ************************************************************************ 72 : 73 : 'SI KAS1' (EGA KAS 'KAS1'); 74 : 75 : * Initialisation des debits 76 : *? Mliqbr=0. ; 77 : aex =1. ; 78 : Ksi =1. ; 79 : qliq =0. ; hliq =0. ; Mliqj = 0.0 ; 80 : qeau =0. ; hvap =0. ; Yvapj = 0.0 ; 81 : qair =0. ; hair =0. ; Yairj = 0.0 ; 82 : qhe =0. ; hhe =0. ; Yhej = 0.0 ; 83 : qh2 =0. ; hh2 =0. ; Yh2j = 0.0 ; 84 : qo2 =0. ; ho2 =0. ; Yo2j = 0.0 ; 85 : qn2 =0. ; hn2 =0. ; Yn2j = 0.0 ; 86 : qco =0. ; hco =0. ; Ycoj = 0.0 ; 87 : qco2 =0. ; hco2 =0. ; Yco2j = 0.0 ; 88 : 89 : * Les donnees physiques 90 : Cpvap = TBT.'Cpvap' ; 91 : Lv = TBT.'Lv' ; 92 : Rabs = TBT.'Rabs' ; 93 : Rgh2 = TBT.'Rgh2' ; Mkh2 = TBT.'Mkh2' ; 94 : Rghe = TBT.'Rghe' ; Mkhe = TBT.'Mkhe' ; 95 : Rgo2 = TBT.'Rgo2' ; Mko2 = TBT.'Mko2' ; 96 : Rgn2 = TBT.'Rgn2' ; Mkn2 = TBT.'Mkn2' ; 97 : Rgco2 = TBT.'Rgco2' ; Mkco2 = TBT.'Mkco2' ; 98 : Rgco = TBT.'Rgco' ; Mkco = TBT.'Mkco' ; 99 : Rgvap = TBT.'Rgvap' ; Mkvap = TBT.'Mkvap' ; 100 : Rgair = TBT.'Rgair' ; Mkair = TBT.'Mkair' ; 101 : 102 : 103 : t = rxtsce .'t' ; 104 : 105 : *-- Formulation en Extraction ------------------------------------------ 106 : *------------------ TTOUT ---------------------------------------------- 107 : 'SI TTOUT KAS1' TTOUT ; 108 : 109 : $breche = Briches.'$sreche'; 110 : brechei = Briches.'srechei'; 111 : brechec = Briches.'srechec'; 112 : breche = Briches.'sreche' ; 113 : 114 : Rout=redu Breche tic.'RHO'; 115 : diru = Briches.'diru' ; 116 : 117 : Si(EGA (VALE DIME) 2); 118 : gjc = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0 119 : 'NATURE' 'DISCRET' ; 120 : Sinon ; 121 : gjc = 'MANU' 'CHPO' brechec 3 'UX' 0.0 'UY' 0.0 'UZ' 0. 122 : 'NATURE' 'DISCRET' ; 123 : Finsi; 124 : 125 : *============================================================== 126 : * C'est la qu'on definit l'extraction 127 : * 128 : Pt =extr tic.'PT' ndl ; 129 : Tfm =extr tic.'Tfm' ndl ; 130 : Rhom=extr tic.'Rhom' ndl ; 131 : Sout = rxtsce.'Section' ; 132 : Pext = 'IPOL' tps t (rxtsce .'Pext') ; 133 : Text = 'IPOL' tps t (rxtsce .'Text') ; 134 : Kext = 'IPOL' tps t (rxtsce .'Kext') ; 135 : Uout = ((Pt - Pext)*Kext)**(1./1.75) ; 136 : uus = 'KCHT' $breche 'VECT' 'SOMMET' (uout*diru) gjc ; 137 : Briches.'gj' = uus ; 138 : gamma=1.29; 139 : Pcr = Pext*((2./(gamma + 1.))**(gamma/(gamma - 1.))) ; 140 : Si (Pcr > Pt); 141 : ge=(gamma+1.)/2./(gamma-1.); 142 : Qcr = Pext*((gamma/rhom/Text)**0.5)*((2/(gamma+1))**ge)*Sout ; 143 : Uout = Qcr/Rhom/Sout ; 144 : uus = 'KCHT' $breche 'VECT' 'SOMMET' (uout*diru) gjc ; 145 : Briches.'gj' = uus ; 146 : 147 : mess ' Il y a un col sonique : Débit critique = ' Qcr ' Uout=' Uout ; 148 : Sinon ; 149 : mess ' Il n y a pas de col sonique : Uout=' Uout ; 150 : Finsi; 151 : * 152 : *============================================================== 153 : Tout = redu tic.'TF' breche ; 154 : Qout = 'DBIT' $breche uus ; 155 : 156 : Si(Qout > 0.) ; 157 : Tout = ('DBIT' $breche (Tout*uus)) / Qout ; 158 : Rout = ('DBIT' $breche (Rout*uus)) / Qout ; 159 : Roj = Rout ; 160 : Tinj= Tout ; 161 : Qj = Qout ; 162 : Sinon ; 163 : Rout=1.; 164 : Tout=0.; 165 : Roj = 1. ; 166 : Tinj= 0. ; 167 : Qj = 0. ; 168 : Finsi ; 169 : 170 : $menvf = rxt.'GEO'.'$menvf' ; 171 : src = 'DBIT' $breche uus 'IMPR' ; 172 : 173 : Si (src < 0.); 174 : mess; 175 : mess 'L Extraction est dans le mauvais sens ! '; 176 : mess; 177 : erre 2; 178 : quitter broche; 179 : Finsi; 180 : 181 : * ung=vect uus 0.5 ux uy uz jaune; 182 : * trace ung (doma $menvf maillage); 183 : 184 : RHO =tic.'RHO' redu Breche ; 185 : YAIR=tic.'YAIR' redu Breche ; 186 : Y = YAIR ; 187 : qair = 'DBIT' $breche (uus*RHO*YAIR) ; 188 : 189 : 'SI' TBT.'THE' ; 190 : YHE =tic.'YHE' redu Breche ; 191 : qhe = 'DBIT' $breche (uus*YHE*RHO) ; Y=Y+YHE ; Finsi ; 192 : 'SI' TBT.'TH2' ; 193 : YH2 =tic.'YH2' redu Breche ; 194 : qh2 = 'DBIT' $breche (uus*YH2*RHO) ; Y=Y+YH2 ; Finsi ; 195 : 'SI' TBT.'TO2' ; 196 : YO2 =tic.'YO2' redu Breche ; 197 : qo2 = 'DBIT' $breche (uus*YO2*RHO) ; Y=Y+YO2 ; Finsi ; 198 : 'SI' TBT.'TN2' ; 199 : YN2 =tic.'YN2' redu Breche ; 200 : qn2 = 'DBIT' $breche (uus*YN2*RHO) ; Y=Y+YN2 ; Finsi ; 201 : 'SI' TBT.'TCO'; 202 : YCO =tic.'YCO' redu Breche ; 203 : qco = 'DBIT' $breche (uus*YCO*RHO) ; Y=Y+YCO ; Finsi ; 204 : 'SI' TBT.'TCO2' ; 205 : YCO2 =tic.'YCO2' redu Breche ; 206 : qco2 = 'DBIT' $breche (uus*YCO2*RHO) ; Y=Y+YCO2 ; Finsi ; 207 : 208 : Qj= qair + qhe + qh2 + qo2 + qn2 + qco + qco2 ; 209 : 210 : 211 : 'SI' TBT.'VAPEUR' ; 212 : RVP =tic.'RVP' redu Breche ; 213 : qeau = 'DBIT' $breche (uus*RVP) ; 214 : YVP= RVP * (inve RHO) ; 215 : Y = Y + YVP ; 216 : Qj = Qj + qeau ; 217 : Pt =extr tic.'PT' ndl ; 218 : Tfm=extr tic.'Tfm' ndl ; 219 : 'FINSI' ; 220 : 221 : Qj = Qout ; 222 : 'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ; 223 : Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej = qhe '*' iQj ; 224 : Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj ; 225 : Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ; 226 : 227 : * On ne change pas le signe de Qj mais celui de qeau et autres 228 : aex=-1.; 229 : qeau = aex * qeau ; 230 : qliq = aex * qliq ; 231 : qair = aex * qair ; 232 : qhe = aex * qhe ; 233 : qh2 = aex * qh2 ; 234 : qo2 = aex * qo2 ; 235 : qn2 = aex * qn2 ; 236 : qco = aex * qco ; 237 : qco2 = aex * qco2 ; 238 : * mess ' MINI MAXI YSortie=' (mini Y) (maxi Y); 239 : 240 : TTINJ=FAUX; 241 : Tinj=Tout; 242 : 'FINSI TTOUT KAS1' ; 243 : *-------------- FIN TTOUT ---------------------------------------------- 244 : 245 : *-- Formulation en Extraction ------------------------------------------ 246 : *------------------ TTINJ ---------------------------------------------- 247 : 'SI TTINJ KAS1' TTINJ ; 248 : 249 : $breche = Briches.'$sreche'; 250 : brechei = Briches.'srechei'; 251 : brechec = Briches.'srechec'; 252 : breche = Briches.'sreche' ; 253 : 254 : Roj=redu Breche tic.'RHO'; 255 : diru =Briches.'diru' ; 256 : 257 : Si(EGA (VALE DIME) 2); 258 : gjc = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0 259 : 'NATURE' 'DISCRET' ; 260 : Sinon ; 261 : gjc = 'MANU' 'CHPO' brechec 3 'UX' 0.0 'UY' 0.0 'UZ' 0. 262 : 'NATURE' 'DISCRET' ; 263 : Finsi; 264 : 265 : $menvf = rxt.'GEO'.'$menvf' ; 266 : 267 : tinj = 'IPOL' tps t (rxtsce .'tinj') ; 268 : qair = 'IPOL' tps t (rxtsce .'qair') ; 269 : 'SI' TBT.'THE' ; 270 : qhe= 'IPOL' tps t (rxtsce .'qhe') ; Finsi ; 271 : 'SI' TBT.'TH2' ; 272 : qh2= 'IPOL' tps t (rxtsce .'qh2') ; Finsi ; 273 : 'SI' TBT.'TO2' ; 274 : qo2= 'IPOL' tps t (rxtsce .'qo2') ; Finsi ; 275 : 'SI' TBT.'TN2' ; 276 : qn2= 'IPOL' tps t (rxtsce .'qn2') ; Finsi ; 277 : 'SI' TBT.'TCO'; 278 : qco= 'IPOL' tps t (rxtsce .'qco') ; Finsi ; 279 : 'SI' TBT.'TCO2' ; 280 : qco2='IPOL' tps t (rxtsce .'qco2') ; Finsi ; 281 : 282 : Qj= qair + qhe + qh2 + qo2 + qn2 + qco + qco2 ; 283 : 284 : 'SI' TBT.'VAPEUR' ; 285 : qeau = 'IPOL' tps t (rxtsce .'qeau') ; 286 : Qj = Qj + qeau ; 287 : Pt =extr tic.'PT' ndl ; 288 : Tfm=extr tic.'Tfm' ndl ; 289 : 290 : 'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ; 291 : Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej = qhe '*' iQj ; 292 : Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj ; 293 : Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ; 294 : 295 : 296 : Mj=((Yvapj/Mkvap)+(Yairj/Mkair)+(Yhej/Mkhe)+(Yh2j/Mkh2)+(Yo2j/Mko2)+ 297 : (Yn2j*Mkn2)+(Yco2j/Mkco2)+(Ycoj/Mkco))**(-1.) ; 298 : Xvj=Yvapj*Mj/Mkvap ; 299 : Pvj = Pt * Xvj ; 300 : 301 : 'SI' (TBT.'CONDMAS') ; 302 : Pt =extr tic.'PT' ndl ; 303 : Tfm=extr tic.'Tfm' ndl ; 304 : Qj=qeau + qair + qhe + qh2 + qo2 + qn2 + qco + qco2; 305 : *a peu près Xvj = Yvj on fera mieux après 306 : Xvj = qeau/Qt ; 307 : Pvj = Pt * Xvj ; 308 : Hinj = VARI HVS Pvj (Tinj + 273.15) ; 309 : Tsat = VARI TSATP Pt ; 310 : Hliq = VARI HLS Pvj Tsat ; 311 : Hvap = VARI HVS Pvj Tsat ; 312 : mess ' Hinj=' Hinj ' Hliq=' Hliq ' Hvap=' Hvap ; 313 : mess 'Pvj=' pvj ' Tsat = ' Tsat ' Tinj=' (Tinj + 273.15); 314 : Ksi=(Hinj - Hliq)/(Hvap - Hliq) ; 315 : Ksi=MINI (PROG 1. Ksi) ; 316 : mess ' Ksi ' Ksi ' qeau= ' qeau ' qvap=' (qeau*Ksi); 317 : **? Mliqbr= Mliqbr + (qeau*(1. - Ksi)) ; 318 : qeau = qeau * Ksi ; 319 : 'FINSI' ; 320 : 321 : 'SINON' ; 322 : 323 : 'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ; 324 : Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej = qhe '*' iQj ; 325 : Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj ; 326 : Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ; 327 : 328 : 329 : 'FINSI' ; 330 : 331 : Si(EGA Qj 0.) ; 332 : Roj = 0. ; 333 : Sinon ; 334 : Pt = 'EXTR' (tic.'PT') ndl ; 335 : Rgj = (TBT.'Rgvap' * Yvapj) + (TBT.'Rgair' * Yairj) + 336 : (TBT.'Rghe' * Yhej) + (TBT.'Rgh2' * Yh2j) + 337 : (TBT.'Rgo2' * Yo2j) + (TBT.'Rgn2' * Yn2j) + 338 : (TBT.'Rgco' * Ycoj) + (TBT.'Rgco2' * Yco2j) ; 339 : 340 : Roj = Pt '/' (Tinj '+' 273.15) '/' Rgj ; 341 : Finsi ; 342 : 343 : 'FINSI TTINJ KAS1' ; 344 : *-------------- FIN TTINJ ---------------------------------------------- 345 : 346 : rxtsce .'Yvapj'= Yvapj ; 347 : rxtsce .'Yairj'= Yairj ; 348 : rxtsce .'Yhej' = Yhej ; 349 : rxtsce .'Yh2j' = Yh2j ; 350 : rxtsce .'Yo2j' = Yo2j ; 351 : rxtsce .'Yn2j' = Yn2j ; 352 : rxtsce .'Yco2j'= Yco2j ; 353 : rxtsce .'Ycoj' = Ycoj ; 354 : 355 : rxtsce .'aex' = aex ; 356 : rxtsce .'Qj' = Qj ; 357 : rxtsce .'Qlj' = Qlj ; 358 : rxtsce .'Tj' = tinj ; 359 : rxtsce .'Roj' = Roj ; 360 : 361 : Briches.CTinj= Briches.CTinj 'ET' ('PROG' Tinj) ; 362 : Briches.CKsi = Briches.CKsi 'ET' ('PROG' Ksi) ; 363 : Briches.CQlj = Briches.CQlj 'ET' ('PROG' qliq) ; 364 : Briches.CHlj = Briches.CHlj 'ET' ('PROG' hliq) ; 365 : *------------------------------------------------------------- 366 : *Caractéristiques de l'injection au pas de temps N 367 : * c'est a dire Dtn = tn - t(n-1) a l'instant n-1 368 : 369 : 'SI' ( TBT.'IMPR' >EG 2) ; 370 : 'MESS' '==> Détail des débits extraits (kg/s) Sortie ' nub ; 371 : 'MESS' 'Qeau=' qeau 'Qliq=' qliq 'Qair=' qair 'Qhe =' qhe 'Qh2 =' qh2 ; 372 : 'MESS' 'Qo2 =' qo2 'Qn2 =' qn2 'Qco =' qco 'Qco2=' qco2 ; 373 : 'FINSI' ; 374 : 375 : *mess '=================== FIN KAS1 ======================='; 376 : resproc qeau qliq qair qhe qh2 qo2 qn2 qco qco2 ; 377 : 378 : QUITTER BROCHE ; 379 : 'FINSI KAS1' ; 380 : ************************************************************************ 381 : ********* FIN KAS1 ***************************************************** 382 : ************************************************************************ 383 : 384 : ************************************************************************ 385 : ********* DEBUT KAS2 *************************************************** 386 : ************************************************************************ 387 : 388 : 'SI KAS2' (EGA KAS 'KAS2') ; 389 : 390 : CSbreche =chai 'Ssreche' ; 391 : C$breche =chai '$sreche' ; 392 : C$brechc =chai '$srechec' ; 393 : Cfacgeo =chai 'facgeo' ; 394 : Cbrechei =chai 'srechei' ; 395 : Cbrechec =chai 'srechec' ; 396 : Cbreche =chai 'sreche' ; 397 : 398 : 399 : diru =Briches.'diru' ; 400 : 401 : diru = diru/(norm diru) ; 402 : Sbreche = GEO.CSbreche ; 403 : $breche = GEO.C$breche ; 404 : facgeo = GEO.Cfacgeo ; 405 : brechei = GEO.Cbrechei ; 406 : brechec = GEO.Cbrechec ; 407 : breche = GEO.Cbreche ; 408 : 409 : $menvf = rxt.'GEO'.'$menvf' ; 410 : $vtf = rxt.'GEO'.'$vtf' ; 411 : 412 : Roj = rxtsce .'Roj' ; 413 : Qj = rxtsce .'Qj' ; 414 : Qlj = rxtsce .'Qlj' ; 415 : tinj= rxtsce .'Tj' ; 416 : 417 : ************************************************************** 418 : **** Evolution de la source en fct du temps ****************** 419 : 420 : t = rxtsce .'t' ; 421 : 422 : * Récupération des fractions massiques 423 : Yvapj= rxtsce .'Yvapj' ; 424 : Yairj= rxtsce .'Yairj' ; 425 : Yhej = rxtsce .'Yhej' ; 426 : Yh2j = rxtsce .'Yh2j' ; 427 : Yo2j = rxtsce .'Yo2j' ; 428 : Yn2j = rxtsce .'Yn2j' ; 429 : Yco2j= rxtsce .'Yco2j' ; 430 : Ycoj = rxtsce .'Ycoj' ; 431 : 432 : Qjf = Qj '*' facgeo ; 433 : 434 : 'SI' ('EGA' Qj 0.0) ; 435 : Roj = 0.0 ; 436 : d guj = 0.0 ; 437 : 'SINON' ; 438 : Roj = rxtsce .'Roj' ; 439 : guj = (Qjf '/' Sbreche '/' Roj ) ; 440 : 'FINSI' ; 441 : 442 : Rvapj = Roj '*' Yvapj ; 443 : Rairj = Roj '*' Yairj ; 444 : Rhej = Roj '*' Yhej ; 445 : Rh2j = Roj '*' Yh2j ; 446 : Ro2j = Roj '*' Yo2j ; 447 : Rn2j = Roj '*' Yn2j ; 448 : Rco2j = Roj '*' Yco2j ; 449 : Rcoj = Roj '*' Ycoj ; 450 : 451 : * Les donnees physiques 452 : Cpvap = TBT.'Cpvap' ; 453 : Lv = TBT.'Lv' ; 454 : Rabs = TBT.'Rabs' ; 455 : Rgh2 = TBT.'Rgh2' ; Mkh2 = TBT.'Mkh2' ; 456 : Rghe = TBT.'Rghe' ; Mkhe = TBT.'Mkhe' ; 457 : Rgo2 = TBT.'Rgo2' ; Mko2 = TBT.'Mko2' ; 458 : Rgn2 = TBT.'Rgn2' ; Mkn2 = TBT.'Mkn2' ; 459 : Rgco2 = TBT.'Rgco2' ; Mkco2 = TBT.'Mkco2' ; 460 : Rgco = TBT.'Rgco' ; Mkco = TBT.'Mkco' ; 461 : Rgvap = TBT.'Rgvap' ; Mkvap = TBT.'Mkvap' ; 462 : Rgair = TBT.'Rgair' ; Mkair = TBT.'Mkair' ; 463 : 464 : * Calcul Cp (Tj) CALCP 465 : Tj = tinj ; 466 : 467 : Cph2j Cphej Cpo2j Cpn2j Cpco2j Cpcoj Cpairj = CALCP Tj ; 468 : 469 : * Calcul des Cvs pour l'injection 470 : Cpvapj = TBT.'Cpvap' ; 471 : Cvvapj = Cpvapj '-' Rgvap ; 472 : Cvhej = Cphej '-' Rghe ; 473 : Cvh2j = Cph2j '-' Rgh2 ; 474 : Cvo2j = Cpo2j '-' Rgo2 ; 475 : Cvn2j = Cpn2j '-' Rgn2 ; 476 : Cvco2j = Cpco2j '-' Rgco2 ; 477 : Cvcoj = Cpcoj '-' Rgco ; 478 : Cvairj = Cpairj '-' Rgair ; 479 : 480 : Cvj = (Cvvapj * Yvapj)+(Cvairj * Yairj)+(Cvhej * Yhej)+ 481 : (Cvh2j * Yh2j)+(Cvo2j * Yo2j)+(Cvn2j * Yn2j)+ 482 : (Cvcoj * Ycoj)+(Cvco2j * Yco2j) ; 483 : Cpj = (Cpvapj * Yvapj)+(Cpairj * Yairj)+(Cphej * Yhej)+ 484 : (Cph2j * Yh2j)+(Cpo2j * Yo2j)+(Cpn2j * Yn2j)+ 485 : (Cpcoj * Ycoj)+(Cpco2j * Yco2j) ; 486 : Rgj = (Rgvap * Yvapj)+(Rgair * Yairj)+(Rghe * Yhej)+ 487 : (Rgh2 * Yh2j)+(Rgo2 * Yo2j)+(Rgn2 * Yn2j)+ 488 : (Rgco * Ycoj)+(Rgco2 * Yco2j) ; 489 : 490 : Hj = Cpj '*' (Tj '+' 273.15) ; 491 : Ej = Cvj '*' (Tj '+' 273.15) ; 492 : 493 : Lbreche=Sbreche**(1./IDIM) ; 494 : Knj=0.05 * guj ; 495 : Enj=(Knj**1.5)/Lbreche ; 496 : 497 : 'SI' ( TBT.'IMPR' >EG 2) ; 498 : 'MESS' '==> Caracteristiques de l extraction ' nub ; 499 : 'MESS' ' Débit total extrait Qj=' Qj 'kg/s' ; 500 : 'MESS' ' Vitesse extraction guj=' guj 'm/s' ; 501 : 'MESS' ' Densité extraction Roj=' Roj 'kg/m3' ; 502 : 'MESS' ' R/M extraction Rgj=' Rgj 'J/kg/K' ; 503 : 'FINSI' ; 504 : 505 : src=0.; 506 : 507 : 'SI' TBT.'FPAROI' ; 508 : brchp=breche ; 509 : 'SINON' ; 510 : brchp=brechei ; 511 : 'FINSI' ; 512 : 513 : 'SI' DIM3D ; 514 : $brechec= GEO.C$brechc ; 515 : gj = 'KCHT' $breche 'VECT' 'SOMMET' 'COMP' 516 : ('MOTS' 1UN 2UN 3UN) (diru * guj) ; 517 : gj = 'REDU' gj brchp ; 518 : gjc = 'KCHT' $brechec 'VECT' 'SOMMET' (0.0 0.0 0.0) ; 519 : gjx = 'KCHT' $breche 'VECT' 'SOMMET' (diru * guj) gjc ; 520 : 'SINON' ; 521 : gj = 'KCHT' $breche 'VECT' 'SOMMET' 'COMP' 522 : ('MOTS' 1UN 2UN) (diru * guj) ; 523 : gj = 'REDU' gj brchp ; 524 : gjc = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0 525 : 'NATURE' 'DISCRET' ; 526 : gjx = 'KCHT' $breche 'VECT' 'SOMMET' (diru * guj) gjc ; 527 : 'FINSI' ; 528 : 529 : unvf0 = 'KCHT' $menvf 'VECT' 'SOMMET' gjx ; 530 : src = 'DBIT' unvf0 $menvf ; 531 : 532 : *============================================================= 533 : 534 : Briches.'Rvapj'=Rvapj ; 535 : Briches.'Rairj'=Rairj ; 536 : Briches.'Rhej' =Rhej ; 537 : Briches.'Rh2j' =Rh2j ; 538 : Briches.'Ro2j' =Ro2j ; 539 : Briches.'Rn2j' =Rn2j ; 540 : Briches.'Rco2j'=Rco2j ; 541 : Briches.'Rcoj' =Rcoj ; 542 : Briches.'Rgj' =Rgj ; 543 : * Briches.'Hj' =Hj ; 544 : * Briches.'Ej' =Ej ; 545 : Briches.'Tj' =Tj ; 546 : * Briches.'Qj' =Qj ; 547 : Briches.'Roj' =Roj ; 548 : Briches.'gj' =gj ; 549 : * Briches.'guj' =guj ; 550 : Briches.'unvf0'=unvf0 ; 551 : Briches.'src' =src ; 552 : aex = rxtsce .'aex' ; 553 : *mess '=================== FIN KAS2 ======================='; 554 : resproc 555 : Rvapj Rairj Rhej Rh2j Ro2j Rn2j Rco2j Rcoj Rgj Hj Ej 556 : Tj aex Qj Roj gj guj unvf0 src ; 557 : 558 : QUITTER BROCHE ; 559 : 560 : 'FINSI KAS2' ; 561 : ************************************************************************ 562 : ********* FIN KAS2 ***************************************************** 563 : ************************************************************************ 564 : 565 : 566 : FINPROC ; 567 : 568 :
© Cast3M 2003 - All rights reserved.
Disclaimer