Télécharger bruche.procedur

Retour à la liste

Numérotation des lignes :

  1. * BRUCHE PROCEDUR MAGN 17/12/12 21:15:04 9663
  2. 'DEBPROC' bruche KAS*'MOT' nab*'ENTIER' rxt*'TABLE';
  3.  
  4. GEO = rxt.'GEO' ;
  5. Si(nab <EG 3) ;
  6. nub = nab ;
  7. rxtsce =rxt.'scenario' ;
  8. Briches =rxt.'TIC' ;
  9. CTinj =chai 'Tinj' ;
  10. CKsi =chai 'Ksi' ;
  11. CQlj =chai 'Qlj' ;
  12. CHlj =chai 'Hlj' ;
  13. ;
  14. Si (nub > 1) ;
  15. rxtsce=rxt.(chai 'scenario' nub) ;
  16. CTinj = chai CTinj nub ;
  17. CKsi = chai CKsi nub ;
  18. CQlj = chai CQlj nub ;
  19. CHlj = chai CHlj nub ;
  20. Finsi ;
  21. Sinon ;
  22. nub = nab - 3 ;
  23. ibr=index (rxt.'Breches') ;
  24. Briches =rxt.'Breches'.(ibr.nub) ;
  25. rxtsce =Briches.'scenario' ;
  26. GEO =Briches ;
  27. CTinj =chai 'Tinj' ;
  28. CKsi =chai 'Ksi' ;
  29. CQlj =chai 'Qlj' ;
  30. CHlj =chai 'Hlj' ;
  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)) ;
  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. ************************************************************************
  71. ********* DEBUT KAS1 ***************************************************
  72. ************************************************************************
  73.  
  74. Si (EGA KAS 'KAS1');
  75.  
  76. * Initialisation des debits
  77. *? Mliqbr=0. ;
  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. TTINJ=FAUX ;
  106. THINJ=FAUX ;
  107. TTINJ=exist (rxtsce ) 'tinj' ;
  108. THINJ=exist (rxtsce ) 'hair' ;
  109.  
  110. Si ((TTINJ et THINJ) ou ((NON TTINJ) et (NON THINJ))) ;
  111. Mess '*********************************************************' ;
  112. Mess ' ERREUR ERREUR ERREUR ERREUR ERREUR ERREUR ' ;
  113. Mess ' ' ;
  114. Mess 'On ne peut imposer enthalpie et température ' ;
  115. Mess 'et il faut imposer au moins l une des deux !! ' ;
  116. Mess '*********************************************************' ;
  117. ERREUR 21 ;
  118. Finsi ;
  119.  
  120. *-- Formulation en Température imposée ---------------------------------
  121. *------------------ TTINJ ----------------------------------------------
  122. 'SI' TTINJ ;
  123. mess ' Formulation en Température imposée -';
  124. 'SI' (TBT.'CATHARE2') ;
  125. mess ' CATHAR2 ';
  126. mess ' HELLO ';
  127. titi='EXTE' 'cat < CASCAR' ;
  128. tinj = titi.1 ; qair = titi.2 ; qeau = titi.3 ;
  129. mess 'tinj qair qeau = ' tinj qair qeau ;
  130. * Finsi ;
  131. 'SINON' ;
  132.  
  133. tinj = 'IPOL' tps t (rxtsce .'tinj') ;
  134. qair = 'IPOL' tps t (rxtsce .'qair') ;
  135.  
  136. 'SI' TBT.'THE' ;
  137. qhe= 'IPOL' tps t (rxtsce .'qhe') ; Finsi ;
  138. 'SI' TBT.'TH2' ;
  139. qh2= 'IPOL' tps t (rxtsce .'qh2') ; Finsi ;
  140. 'SI' TBT.'TO2' ;
  141. qo2= 'IPOL' tps t (rxtsce .'qo2') ; Finsi ;
  142. 'SI' TBT.'TN2' ;
  143. qn2= 'IPOL' tps t (rxtsce .'qn2') ; Finsi ;
  144. 'SI' TBT.'TCO' ;
  145. qco= 'IPOL' tps t (rxtsce .'qco') ; Finsi ;
  146. 'SI' TBT.'TCO2' ;
  147. qco2='IPOL' tps t (rxtsce .'qco2') ; Finsi ;
  148. 'SI' TBT.'VAPEUR' ;
  149. qeau = 'IPOL' tps t (rxtsce .'qeau'); Finsi ;
  150. 'FINSI' ;
  151.  
  152. Qj= qair + qhe + qh2 + qo2 + qn2 + qco + qco2 + qeau ;
  153.  
  154.  
  155. 'SI' TBT.'VAPEUR' ;
  156. Pt =extr tic.'PT' ndl ;
  157. Tfm=extr tic.'Tfm' ndl ;
  158.  
  159. 'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ;
  160. Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej = qhe '*' iQj ;
  161. Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj ;
  162. Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ;
  163.  
  164.  
  165. Mj=((Yvapj/Mkvap)+(Yairj/Mkair)+(Yhej/Mkhe)+(Yh2j/Mkh2)+(Yo2j/Mko2)+
  166. (Yn2j*Mkn2)+(Yco2j/Mkco2)+(Ycoj/Mkco))**(-1.) ;
  167. Xvj=Yvapj*Mj/Mkvap ;
  168. Pvj = Pt * Xvj ;
  169.  
  170. 'SI' (TBT.'CONDMAS') ;
  171. Pt =extr tic.'PT' ndl ;
  172. Tfm=extr tic.'Tfm' ndl ;
  173. Qj=qeau + qair + qhe + qh2 + qo2 + qn2 + qco + qco2;
  174. *a peu près Xvj = Yvj on fera mieux après
  175. Xvj = qeau/Qt ;
  176. Pvj = Pt * Xvj ;
  177. Hinj = VARI HVS Pvj (Tinj + 273.15) ;
  178. Tsat = VARI TSATP Pt ;
  179. Hliq = VARI HLS Pvj Tsat ;
  180. Hvap = VARI HVS Pvj Tsat ;
  181. mess ' Hinj=' Hinj ' Hliq=' Hliq ' Hvap=' Hvap ;
  182. mess 'Pvj=' pvj ' Tsat = ' Tsat ' Tinj=' (Tinj + 273.15);
  183. Ksi=(Hinj - Hliq)/(Hvap - Hliq) ;
  184. Ksi=MINI (PROG 1. Ksi) ;
  185. mess ' Ksi ' Ksi ' qeau= ' qeau ' qvap=' (qeau*Ksi);
  186. **? Mliqbr= Mliqbr + (qeau*(1. - Ksi)) ;
  187. qeau = qeau * Ksi ;
  188. 'FINSI' ;
  189.  
  190. 'SINON' ;
  191.  
  192. 'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ;
  193. Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej = qhe '*' iQj ;
  194. Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj ;
  195. Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ;
  196.  
  197.  
  198. 'FINSI' ;
  199.  
  200. Si(EGA Qj 0.) ;
  201. Roj = 0. ;
  202. Sinon ;
  203. Pt = 'EXTR' (tic.'PT') ndl ;
  204. Rgj = (TBT.'Rgvap' * Yvapj) + (TBT.'Rgair' * Yairj) +
  205. (TBT.'Rghe' * Yhej) + (TBT.'Rgh2' * Yh2j) +
  206. (TBT.'Rgo2' * Yo2j) + (TBT.'Rgn2' * Yn2j) +
  207. (TBT.'Rgco' * Ycoj) + (TBT.'Rgco2' * Yco2j) ;
  208.  
  209. Roj = Pt '/' (Tinj '+' 273.15) '/' Rgj ;
  210. Finsi ;
  211.  
  212. Finsi ;
  213. *-------------- FIN TTINJ ----------------------------------------------
  214.  
  215. *-- Formulation en Enthalpie imposée -----------------------------------
  216. *------------------ THINJ ----------------------------------------------
  217. Si THINJ ;
  218. * mess ' Enthalpie Imposee';
  219. qair = 'IPOL' tps t (rxtsce .'qair') ;
  220. hair = 'IPOL' tps t (rxtsce .'hair') ;
  221. 'SI' TBT.'THE' ;
  222. qhe= 'IPOL' tps t (rxtsce .'qhe') ;
  223. hhe= 'IPOL' tps t (rxtsce .'hhe') ; Finsi ;
  224. 'SI' TBT.'TH2' ;
  225. qh2= 'IPOL' tps t (rxtsce .'qh2') ;
  226. hh2= 'IPOL' tps t (rxtsce .'hh2') ; Finsi ;
  227. 'SI' TBT.'TO2' ;
  228. qo2= 'IPOL' tps t (rxtsce .'qo2') ;
  229. ho2= 'IPOL' tps t (rxtsce .'ho2') ; Finsi ;
  230. 'SI' TBT.'TN2' ;
  231. qn2= 'IPOL' tps t (rxtsce .'qn2') ;
  232. hn2= 'IPOL' tps t (rxtsce .'hn2') ; Finsi ;
  233. 'SI' TBT.'TCO';
  234. qco= 'IPOL' tps t (rxtsce .'qco') ;
  235. hco= 'IPOL' tps t (rxtsce .'hco') ; Finsi ;
  236. 'SI' TBT.'TCO2' ;
  237. qco2='IPOL' tps t (rxtsce .'qco2') ;
  238. hco2='IPOL' tps t (rxtsce .'hco2') ; Finsi ;
  239.  
  240. Qji = qair + qhe + qh2 + qo2 + qn2 + qco + qco2 ;
  241. Qj = Qji ;
  242.  
  243. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  244. * Calcul Cpj, Cp du melange : CALCP -- CAS SANS VAPEUR
  245. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  246. 'SI' (NON (TBT.'VAPEUR')) ;
  247.  
  248. * Calcul des fractions massiques des incondensables a l'injection
  249.  
  250. 'SI'('NON'('EGA' Qji 0.)) ; iQj=1./Qji ; 'SINON' ; iQj=0. ; 'FINSI' ;
  251. Yairj = qair '*' iQj ; Yhej = qhe '*' iQj ;
  252. Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj ;
  253. Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ;
  254.  
  255. * On traite ici les incondensables. On se donne l'enthalpie.
  256. * Pour remonter a la temperature on considere le Cp constant
  257. * pour chaque incondensable et on le choisit a 50degreC
  258. * Cette hypothese est justifiee dans l'intervalle de temperature
  259. * 0. a 250 degresC
  260.  
  261. Cph2j Cphej Cpo2j Cpn2j Cpco2j Cpcoj Cpairj = CALCP tinj;
  262.  
  263. * On calcule le Cp du melange a l'injection
  264. Cpj = (Cpairj * Yairj)+(Cphej * Yhej)+
  265. (Cph2j * Yh2j )+(Cpo2j * Yo2j)+
  266. (Cpn2j * Yn2j )+(Cpcoj * Ycoj)+
  267. (Cpco2j * Yco2j) ;
  268.  
  269. * On etablit l'enthalpie a l'injection pour les incondensables
  270. Hj = hair + hhe + hh2 + ho2 + hn2 + hco + hco2 ;
  271.  
  272. * On en deduit la temperature a l'injection dans le cas ou il n'y
  273. * a que des incondensables
  274. tinj=(Hj/Cpj) - 273.15 ;
  275. 'FINSI' ;
  276. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  277. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  278.  
  279. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  280. * Calcul Cpj, Cp du melange : CALCP -- CAS AVEC VAPEUR
  281. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  282. 'SI' TBT.'VAPEUR' ;
  283.  
  284. qeau = 'IPOL' tps t (rxtsce .'qeau') ;
  285. hvap = 'IPOL' tps t (rxtsce .'hvap') ;
  286.  
  287. * On etablit l'enthalpie a l'injection
  288. Hj = hair + hhe + hh2 + ho2 + hn2 + hco + hco2 + hvap ;
  289.  
  290. * Initialisations
  291. * On prend comme temperature d'injection la temperature au
  292. * pas de temps precedent.
  293.  
  294. Pt = extr tic.'PT' ndl ;
  295. Si(EXIST rxt.'scenario' 'Tj') ;
  296. tinj = rxt.'scenario'.'Tj' ;
  297. Sinon ;
  298. tinj = extr tic.'Tfm' ndl ;
  299. Finsi ;
  300. tinja= tinj ;
  301. qvapj = qeau ;
  302.  
  303. * Repeter CCP 10 ;
  304. Repeter CCP 1 ;
  305. * Calcul des fractions massiques a l'injection
  306. Qj = qvapj + Qji ;
  307. 'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ;
  308. Yvapj = qvapj '*' iQj ; Yairj = qair '*' iQj ; Yhej = qhe '*' iQj;
  309. Yh2j = qh2 '*' iQj ; Yo2j = qo2 '*' iQj ; Yn2j = qn2 '*' iQj;
  310. Yco2j = qco2 '*' iQj ; Ycoj = qco '*' iQj ;
  311.  
  312. * On calcule d'abord les Cp des incondensables a la temperature d injection
  313. * on considere le Cp constant pour chaque incondensable et on le choisit a 50deg
  314. * Cette hypothese est justifiee dans l'intervalle de temperature
  315. * 0. a 250 degresC
  316.  
  317. Cph2j Cphej Cpo2j Cpn2j Cpco2j Cpcoj Cpairj = CALCP tinj ;
  318.  
  319. * On va chercher le Cp de la vapeur iterativement
  320. Cpvapj='VARI' 'DHVDT' Pt (tinj + 273.15) ;
  321. Cpvakj='VARI' 'DHVDT' Pt (tinj + 273.15) ;
  322. * Cpvapj=2200.;
  323.  
  324. * On calcule le Cp du melange a l'injection
  325. Cpj = (Cpvapj * Yvapj)+(Cpairj * Yairj)+
  326. (Cphej * Yhej )+(Cph2j * Yh2j )+
  327. (Cpo2j * Yo2j )+(Cpn2j * Yn2j )+
  328. (Cpcoj * Ycoj )+(Cpco2j * Yco2j) ;
  329. mess ' Iteration ' &ccp ' Cpj=' Cpj 'Tinj=' Tinj ' Cpvakj=' Cpvakj;
  330.  
  331. * Mais Toute la vapeur ne peut pas rentrer sous forme de vapeur
  332.  
  333. Mj=((Yvapj/Mkvap)+(Yairj/Mkair)+(Yhej/Mkhe)+(Yh2j/Mkh2)+(Yo2j/Mko2)+
  334. (Yn2j*Mkn2)+(Yco2j/Mkco2)+(Ycoj/Mkco))**(-1.) ;
  335. Xvj=Yvapj*Mj/Mkvap ;
  336. Pvj = Pt * Xvj ;
  337. * La il s'agit de la pression partielle de vapeur a l'injection
  338. * Si on injeccte de la vapeur pure c'est la pression de l'enceinte Pt
  339. Tsat = 'VARI' 'TSATP' Pvj ;
  340. Hliqs= 'VARI' 'HLS' Pvj Tsat ;
  341. Hvaps= 'VARI' 'HVS' Pvj Tsat ;
  342.  
  343. * On en deduit le taux de vide a l'injection
  344. Si(Hj > Hvaps);
  345. Ksi = 1. ;
  346. Sinon ;
  347. Si(Hj < Hliqs);
  348. Ksi = 0. ;
  349. Sinon;
  350. Ksi=(Hj - Hliqs)/(Hvaps - Hliqs) ;
  351. Ksi=MINI (PROG 1. Ksi) ;
  352. Finsi;
  353. Finsi ;
  354.  
  355.  
  356. Si ('EXIST' TBT 'Modinj') ; Modinj = TBT.'Modinj' ;Finsi;
  357. Si ('EXIST' TBT 'Modinj2'); Modinj = TBT.'Modinj2';Finsi;
  358. Si ('EXIST' TBT 'Modinj3'); Modinj = TBT.'Modinj3';Finsi;
  359. Si ('EXIST' Briches 'Modinj'); Modinj = Briches.'Modinj';Finsi;
  360.  
  361.  
  362.  
  363. * On en deduit la temperature a l'injection
  364. Si (EGA TBT.'Modinj' 'A') ;
  365. * Temperature d'ebullition
  366. mess ' Injection - modèle A : Température d ébullition' ;
  367. tinj=(VARI 'TSATP' Pt) - 273.15 ;
  368. Finsi ;
  369. Si (EGA TBT.'Modinj' 'B') ;
  370. * Temperature moyenne enceinte
  371. mess ' Injection - modèle B : Température moyenne' ;
  372. tinj=extr tic.'Tfm' ndl ;
  373. Finsi ;
  374. Si (EGA TBT.'Modinj' 'C') ;
  375. * Temperature Saturation a Pvj
  376. mess ' Injection - modèle C : Température de saturation';
  377. tinj =('VARI' 'TSATP' Pvj) - 273.15 ;
  378. Finsi ;
  379.  
  380. Si(Hj > Hvaps);
  381. * Temperature H/Cp
  382. mess ' Injection surchauffée : Température inj. =' tinj ;
  383. tinj =(Hj - Hvaps)/(Cpj + 1.e-5) + Tsat - 273.15 ;
  384. Finsi ;
  385.  
  386. dtinjr=(tinj - tinja)/(tinj + tinja)*2. ;
  387. dtinjr=abs dtinjr ;
  388. tinja=tinj ;
  389.  
  390. mess ' Xvj=' Xvj ' Yvapj=' Yvapj ' Tsat=' (Tsat - 273.15) ;
  391. mess ' Hj =' Hj ' hvap=' hvap ' Hliq=' Hliq ' Hvaps=' Hvaps ;
  392. mess ' Pvj=' pvj ' Pt =' Pt ;
  393. mess ' Ksi ' Ksi ' qeau= ' qeau ' qvapj=' (qeau*Ksi) ;
  394. * Mliqbr= Mliqbr + (qeau*(1. - Ksi)) ;
  395.  
  396. qvapj = qeau * Ksi ;
  397. Si (dtinjr < 1.e-3); quitter CCP ; Finsi ;
  398. Fin CCP;
  399.  
  400.  
  401. Si(EGA Xvj 0.) ;
  402. Rovap = 0. ;
  403. Sinon ;
  404. tknj=tinj + 273.15 ;
  405. tk2j=tknj * tknj ;
  406. a=(TBT.'B0'*tk2j)+(TBT.'B1'*tknj)+TBT.'B2' ;
  407. a=(abs a)/TBT.'Rgvap'/tk2j/tknj ;
  408. dlt=((4.*a*Pvj/TBT.'Rgvap'/tknj) + 1.)**0.5 ;
  409. Rovap=(dlt - 1.)/2./a. ;
  410. Finsi ;
  411.  
  412. Rginc = (TBT.'Rgair' * Yairj)+(TBT.'Rghe' * Yhej)+
  413. (TBT.'Rgh2' * Yh2j) +(TBT.'Rgo2' * Yo2j)+
  414. (TBT.'Rgn2' * Yn2j) +(TBT.'Rgco' * Ycoj)+
  415. (TBT.'Rgco2' * Yco2j) ;
  416. Si(EGA Rginc 0.) ;
  417. Roinc = 0. ;
  418. Sinon ;
  419. Roinc=(1. - Xvj)*Pt/Rginc/tknj ;
  420. Finsi ;
  421. Roj = Rovap + Roinc ;
  422.  
  423. qliq = 'IPOL' tps t (rxtsce .'qliq') ;
  424. qliq = qliq + (qeau*(1. - Ksi)) ;
  425. hliq = 'IPOL' tps t (rxtsce .'hliq') ;
  426.  
  427. qeau = qvapj ;
  428. Qlj = qliq ;
  429.  
  430. * Mliqbr= Mliqbr + (qeau*(1. - Ksi)) ;
  431. * Mliq = Mliq + (qeau*(1. - Ksi)) + qliq ;
  432.  
  433. 'FINSI' ;
  434. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  435. *%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  436.  
  437. Finsi ;
  438. *-------------- FIN THINJ --------------------
  439.  
  440. rxtsce .'Yvapj'= Yvapj ;
  441. rxtsce .'Yairj'= Yairj ;
  442. rxtsce .'Yhej' = Yhej ;
  443. rxtsce .'Yh2j' = Yh2j ;
  444. rxtsce .'Yo2j' = Yo2j ;
  445. rxtsce .'Yn2j' = Yn2j ;
  446. rxtsce .'Yco2j'= Yco2j ;
  447. rxtsce .'Ycoj' = Ycoj ;
  448.  
  449. Qj=qeau + qair + qhe + qh2 + qo2 + qn2 + qco + qco2;
  450. rxtsce .'Qj' = Qj ;
  451. rxtsce .'Qlj' = Qlj ;
  452. rxtsce .'Tj' = tinj ;
  453. rxtsce .'Roj' = Roj ;
  454.  
  455. Briches.CTinj= Briches.CTinj 'ET' ('PROG' Tinj) ;
  456. Briches.CKsi = Briches.CKsi 'ET' ('PROG' Ksi) ;
  457. Briches.CQlj = Briches.CQlj 'ET' ('PROG' qliq) ;
  458. Briches.CHlj = Briches.CHlj 'ET' ('PROG' hliq) ;
  459. *-------------------------------------------------------------
  460. *Caractéristiques de l'injection au pas de temps N
  461. * c'est a dire Dtn = tn - t(n-1) a l'instant n-1
  462.  
  463. 'SI' ( TBT.'IMPR' >EG 2) ;
  464. 'MESS' '==>tail des débits injectés (kg/s) Brèche ' nub ;
  465. 'MESS' 'Qeau=' qeau 'Qliq=' qliq 'Qair=' qair 'Qhe =' qhe 'Qh2 =' qh2 ;
  466. 'MESS' 'Qo2 =' qo2 'Qn2 =' qn2 'Qco =' qco 'Qco2=' qco2 ;
  467. 'FINSI' ;
  468.  
  469. resproc qeau qliq qair qhe qh2 qo2 qn2 qco qco2 ;
  470.  
  471. QUITTER BRUCHE ;
  472. Finsi ;
  473. ************************************************************************
  474. ********* FIN KAS1 *****************************************************
  475. ************************************************************************
  476.  
  477. ************************************************************************
  478. ********* DEBUT KAS2 ***************************************************
  479. ************************************************************************
  480.  
  481. Si (EGA KAS 'KAS2');
  482.  
  483. CSbreche =chai 'Sbreche' ;
  484. C$breche =chai '$breche' ;
  485. C$brechc =chai '$brechec' ;
  486. Cfacgeo =chai 'facgeo' ;
  487. Cbrechei =chai 'brechei' ;
  488. Cbrechec =chai 'brechec' ;
  489. Cbreche =chai 'breche' ;
  490.  
  491. Si ((nub > 1) et (nab < 4)) ;
  492. CSbreche = chai 'Sbreche' nub ;
  493. C$breche = chai '$breche' nub ;
  494. C$brechc = chai '$brech' nub 'c' ;
  495. Cfacgeo = chai 'facgeo' nub ;
  496. Cbrechei = chai 'brech' nub 'i' ;
  497. Cbrechec = chai 'brech' nub 'c' ;
  498. Cbreche = chai 'breche' nub ;
  499.  
  500. Finsi ;
  501.  
  502. Si ((nub > 0) et (nab < 4)) ;
  503. diru =TBT.(chai 'diru' nub) ;
  504. Sinon ;
  505. diru =Briches.'diru' ;
  506. Finsi ;
  507.  
  508. diru = diru/(norm diru) ;
  509. Sbreche = GEO.CSbreche ;
  510. $breche = GEO.C$breche ;
  511. facgeo = GEO.Cfacgeo ;
  512. brechei = GEO.Cbrechei ;
  513. brechec = GEO.Cbrechec ;
  514. breche = GEO.Cbreche ;
  515.  
  516. $menvf = rxt.'GEO'.'$menvf' ;
  517. $vtf = rxt.'GEO'.'$vtf' ;
  518.  
  519. Qj = rxtsce .'Qj' ;
  520. Qlj = rxtsce .'Qlj' ;
  521. tinj= rxtsce .'Tj' ;
  522. **************************************************************
  523. **** Evolution de la source en fct du temps ******************
  524.  
  525. t = rxtsce .'t' ;
  526.  
  527. * Initialisation des fractions massiques
  528. Yvapj= rxtsce .'Yvapj' ;
  529. Yairj= rxtsce .'Yairj' ;
  530. Yhej = rxtsce .'Yhej' ;
  531. Yh2j = rxtsce .'Yh2j' ;
  532. Yo2j = rxtsce .'Yo2j' ;
  533. Yn2j = rxtsce .'Yn2j' ;
  534. Yco2j= rxtsce .'Yco2j' ;
  535. Ycoj = rxtsce .'Ycoj' ;
  536.  
  537.  
  538. * Les donnees physiques
  539. Cpvap = TBT.'Cpvap' ;
  540. Lv = TBT.'Lv' ;
  541. Rabs = TBT.'Rabs' ;
  542. Rgh2 = TBT.'Rgh2' ; Mkh2 = TBT.'Mkh2' ;
  543. Rghe = TBT.'Rghe' ; Mkhe = TBT.'Mkhe' ;
  544. Rgo2 = TBT.'Rgo2' ; Mko2 = TBT.'Mko2' ;
  545. Rgn2 = TBT.'Rgn2' ; Mkn2 = TBT.'Mkn2' ;
  546. Rgco2 = TBT.'Rgco2' ; Mkco2 = TBT.'Mkco2' ;
  547. Rgco = TBT.'Rgco' ; Mkco = TBT.'Mkco' ;
  548. Rgvap = TBT.'Rgvap' ; Mkvap = TBT.'Mkvap' ;
  549. Rgair = TBT.'Rgair' ; Mkair = TBT.'Mkair' ;
  550.  
  551. * Calcul Cp (Tj) CALCP
  552. Tj = tinj ;
  553.  
  554. Cph2j Cphej Cpo2j Cpn2j Cpco2j Cpcoj Cpairj = CALCP Tj ;
  555.  
  556. * Calcul des Cvs pour l'injection
  557. Cpvapj = TBT.'Cpvap' ;
  558. Cvvapj = Cpvapj '-' Rgvap ;
  559. Cvhej = Cphej '-' Rghe ;
  560. Cvh2j = Cph2j '-' Rgh2 ;
  561. Cvo2j = Cpo2j '-' Rgo2 ;
  562. Cvn2j = Cpn2j '-' Rgn2 ;
  563. Cvco2j = Cpco2j '-' Rgco2 ;
  564. Cvcoj = Cpcoj '-' Rgco ;
  565. Cvairj = Cpairj '-' Rgair ;
  566.  
  567. Cvj = (Cvvapj * Yvapj)+(Cvairj * Yairj)+(Cvhej * Yhej)+
  568. (Cvh2j * Yh2j)+(Cvo2j * Yo2j)+(Cvn2j * Yn2j)+
  569. (Cvcoj * Ycoj)+(Cvco2j * Yco2j) ;
  570. Cpj = (Cpvapj * Yvapj)+(Cpairj * Yairj)+(Cphej * Yhej)+
  571. (Cph2j * Yh2j)+(Cpo2j * Yo2j)+(Cpn2j * Yn2j)+
  572. (Cpcoj * Ycoj)+(Cpco2j * Yco2j) ;
  573. Rgj = (Rgvap * Yvapj)+(Rgair * Yairj)+(Rghe * Yhej)+
  574. (Rgh2 * Yh2j)+(Rgo2 * Yo2j)+(Rgn2 * Yn2j)+
  575. (Rgco * Ycoj)+(Rgco2 * Yco2j) ;
  576.  
  577. Qjf = Qj '*' facgeo ;
  578.  
  579. 'SI' ('EGA' Qj 0.0) ;
  580. Roj = 0.0 ;
  581. guj = 0.0 ;
  582. 'SINON' ;
  583. *µµµµµ Roj µµµµµµµµµµµµµµµµµµµµ
  584. Si TBT.'Tsortie' ;
  585. Roj = Rhom ;
  586. guj = (Qjf '/' Sbreche '/' Roj ) ;
  587. Sinon ;
  588. Roj = rxtsce .'Roj' ;
  589. guj = (Qjf '/' Sbreche '/' Roj ) ;
  590. Finsi ;
  591. * mess 'Roj=' Roj ' Qjf=' Qjf 'Sbreche=' Sbreche 'Tps=' Tps ;
  592. *µµµµµ Roj µµµµµµµµµµµµµµµµµµµµ
  593. 'FINSI' ;
  594.  
  595.  
  596. Qjf = Qj '*' facgeo ;
  597. Hj = Cpj '*' (Tj '+' 273.15) ;
  598. Ej = Cvj '*' (Tj '+' 273.15) ;
  599.  
  600. Lbreche=Sbreche**(1./IDIM) ;
  601. Knj=0.05 * guj ;
  602. Enj=(Knj**1.5)/Lbreche ;
  603.  
  604. Rvapj = Roj '*' Yvapj ;
  605. Rairj = Roj '*' Yairj ;
  606. Rhej = Roj '*' Yhej ;
  607. Rh2j = Roj '*' Yh2j ;
  608. Ro2j = Roj '*' Yo2j ;
  609. Rn2j = Roj '*' Yn2j ;
  610. Rco2j = Roj '*' Yco2j ;
  611. Rcoj = Roj '*' Ycoj ;
  612.  
  613. 'SI' ( TBT.'IMPR' >EG 2) ;
  614. 'MESS' '==> Caracteristiques de l injection ' nub ;
  615. 'MESS' ' Débit total Injecté Qj=' Qj 'kg/s' ;
  616. 'MESS' ' Vitesse Injection guj=' guj 'm/s' ;
  617. 'MESS' ' Densité Injection Roj=' Roj 'kg/m3' ;
  618. 'MESS' ' R/M Injection Rgj=' Rgj 'J/kg/K' ;
  619. 'FINSI' ;
  620.  
  621. *=============================================================
  622. src=0.;
  623.  
  624. 'SI' TBT.'FPAROI' ;
  625. brchp=breche ;
  626. 'SINON' ;
  627. brchp=brechei ;
  628. 'FINSI' ;
  629.  
  630. 'SI' DIM3D ;
  631. $brechec= GEO.C$brechc ;
  632. gj = 'KCHT' $breche 'VECT' 'SOMMET' 'COMP'
  633. ('MOTS' 1UN 2UN 3UN) (diru * guj) ;
  634. gj = 'REDU' gj brchp ;
  635. gjc = 'KCHT' $brechec 'VECT' 'SOMMET' (0.0 0.0 0.0) ;
  636. gjx = 'KCHT' $breche 'VECT' 'SOMMET' (diru * guj) gjc ;
  637. 'SINON' ;
  638. gj = 'KCHT' $breche 'VECT' 'SOMMET' 'COMP'
  639. ('MOTS' 1UN 2UN) (diru * guj) ;
  640. gj = 'REDU' gj brchp ;
  641. gjc = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0
  642. 'NATURE' 'DISCRET' ;
  643. gjx = 'KCHT' $breche 'VECT' 'SOMMET' (diru * guj) gjc ;
  644. 'FINSI' ;
  645.  
  646. unvf0 = 'KCHT' $menvf 'VECT' 'SOMMET' gjx ;
  647. src = 'DBIT' unvf0 $menvf ;
  648.  
  649. *=============================================================
  650.  
  651. Si(nab > 3) ;
  652. Briches.'Rvapj'=Rvapj ;
  653. Briches.'Rairj'=Rairj ;
  654. Briches.'Rhej' =Rhej ;
  655. Briches.'Rh2j' =Rh2j ;
  656. Briches.'Ro2j' =Ro2j ;
  657. Briches.'Rn2j' =Rn2j ;
  658. Briches.'Rco2j'=Rco2j ;
  659. Briches.'Rcoj' =Rcoj ;
  660. Briches.'Rgj' =Rgj ;
  661. * Briches.'Hj' =Hj ;
  662. * Briches.'Ej' =Ej ;
  663. Briches.'Tj' =Tj ;
  664. * Briches.'Qj' =Qj ;
  665. Briches.'Roj' =Roj ;
  666. Briches.'gj' =gj ;
  667. * Briches.'guj' =guj ;
  668. Briches.'unvf0'=unvf0 ;
  669. Briches.'src' =src ;
  670. Finsi ;
  671. resproc
  672. Rvapj Rairj Rhej Rh2j Ro2j Rn2j Rco2j Rcoj Rgj Hj Ej
  673. Tj Qj Roj gj guj unvf0 src ;
  674.  
  675. QUITTER BRUCHE ;
  676.  
  677. ************************************************************************
  678. ********* FIN KAS2 *****************************************************
  679. ************************************************************************
  680. 'FINSI';
  681.  
  682. FINPROC ;
  683.  
  684.  
  685.  
  686.  
  687.  

İ Cast3M 2003 - Tous droits réservés.
Mentions légales