Télécharger retsat.procedur

Retour à la liste

Numérotation des lignes :

  1. * RETSAT PROCEDUR MAGN 18/06/08 21:15:13 9840
  2. 'DEBPROC' retsat ;
  3. 'ARGU' TBT*'TABLE' tic*'TABLE' $mod*'MMODEL' Pto*'FLOTTANT'
  4. ndl*'ENTIER' Rhomn*'FLOTTANT' ITER*'ENTIER' ;
  5. *-----------------------------------------------------------------------
  6. * LOIETAT objet de type 'MOT' indique le type de la loi d'état pour
  7. * la vapeur
  8. * LOIETAT='GP' gaz parfait -> ZVAP=1.
  9. * LOIETAT='GR' gaz reel -> ZVAP<1.
  10. * tic. objet de type TABLE
  11. *
  12. *Paramètres
  13. * TBT.'YBRO' : fraction massique de brouillard autorisée (Déf 1.E-4)
  14. * TBT.'V1alfa' : Coefficient de relaxation sur la pression (Déf 0.8)(wz)
  15. * TBT.'TAUXcm' : part de la vapeur sursaturée que l'on condense
  16. * Tols : Tolérance à la saturation en Pa = 0. Pascal
  17. *
  18. *Objets de type CHPOINT SCAL SOMMET
  19. * Ptc = kcht $mod scal sommet Pt
  20. * PSTF Pression de saturation PSTF = Psat(Tf °K)
  21. * RSTF Densité de vapeur à saturation RSTF = Psat/(Rvap Tf(°C)Zv)
  22. * tic.'TF' température mélange pas de temps n
  23. * tic.'RHO' densité du mélange pas de temps n
  24. * tic.'RVP' densité de vapeur + brouillard n
  25. * tic.'RAIR' densité air pas de temps n
  26. * . . pas de temps n
  27. * RVP = tic.'RVP' vapeur pure + RBRO pas de temps n
  28. * ROVP = tic.'RVP' - tic.'RBRO' vapeur pure pas de temps n
  29. * Rbro = tic.'RBRO' densité brouillard au pas de temps n-1
  30. * Rhot = densité du mélange gazeux -> Pt
  31. *
  32. *Objets de type LISTREEL
  33. * Ptn1 = extr tic.'PT' ndl (ndl correspond à n-1)
  34. * Rhomnm
  35. * tic.'Rhomn' densité du mélange gazeux (vapeur+incondensables) -> Pt
  36. * tic.'Rhomv'
  37. * tic.'Rhomvg'
  38. *
  39. *-----------------------------------------------------------------------
  40.  
  41.  
  42. Tols = 0. ;
  43.  
  44. mess ;
  45. mess '...............................................................';
  46.  
  47. * On initialise à gaz parfait
  48. LOIETAT = TBT.'LOI-ETAT' ;
  49. LDT = 'Gaz Parfait' ;
  50. Zv = 1. ;
  51. Si('EGA' LOIETAT 'GR') ;
  52. LDT = 'Gaz Reel' ;
  53. Finsi ;
  54.  
  55. Mess ' Procédure RETSAT' ;
  56. Mess ' Loi d état pour la vapeur: ' LDT ;
  57. *-----------------------------------------------------------------------
  58. *
  59. Diag = doma $mod 'XXDIAGSI' ;
  60. Vtotal= somt Diag ;
  61.  
  62. Tf = tic.'TF' ;
  63. mess ' Controle des temperatures: Mini/Maxi Tf' (mini Tf) (Maxi Tf);
  64.  
  65.  
  66. Si ((Maxi Tf) > 300.);
  67. mess ' Attention La température dépasse 300°C : On sort des tables';
  68. * La température est limitée à 300°C pour rester dans les tables !!
  69. a = 300.; al=0.99; ala=al*a ;
  70. * b=ala*al*((1.-al)**(-1.));
  71. * b2=ala*(1.+ (ala*al) - ala)*((1.-al)**(-1.));
  72. b3 = ala*(1.+ (ala*ala*al) - (ala*ala))*((1.-al)**(-1.));
  73. ik=masq Tf 'INFERIEUR' ala ;
  74. *Tf=(ik*Tf)+((1.-ik)*a*(Tf+b)*(inve(a+Tf+b)));
  75. *Tf=(ik*Tf)+((1.-ik)*a*((Tf*Tf+b2)*(inve(a+(Tf*Tf+b2)))));
  76. Tf=(ik*Tf)+((1.-ik)*a*((Tf*Tf*Tf+b3)*(inve(a+(Tf*Tf*Tf+b3)))));
  77. mess '%$$ Mini/Maxi Tf' (mini Tf) (Maxi Tf) ;
  78. Finsi;
  79.  
  80. Tfm = ('SOMT' (Diag * tic.'TF'))/Vtotal ;
  81. Tkfm = Tfm + 273.15 ;
  82. *......................................
  83.  
  84.  
  85. *-----------------------------------------------------------------------
  86. *-----------------------------------------------------------------------
  87. mess 'ZZZ On fait un bilan de la distribution des masses en entrée'
  88. ' RETSAT: Densités au pas de temps n' ;
  89.  
  90. * Constante des gaz parfaits au pas de temps n-1
  91. Rgp = tic.'Rgp' ;
  92.  
  93.  
  94. *-----------------------------------------------------------------------
  95. *-----------------------------------------------------------------------
  96. * On calcule les invariants de la part d'incondensable : Pinc et Roinc
  97.  
  98. * CALCP Calcul des Cp des incondensables
  99. Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair = CALCP Tfm ;
  100. Cpvap = TBT.'Cpvap' ;
  101.  
  102. Mair = somt (Diag * tic.'RAIR') ;
  103. Pairo = Mair * TBT.'Rgair' / vtotal ;
  104. Pinco = Pairo ;
  105. Roinc = tic.'RAIR' ;
  106. Rocpinc= tic.'RAIR'*Cpair ;
  107.  
  108. Mo2 = 0. ;
  109. Si TBT.'TO2' ;
  110. Mo2 = 'SOMT' (Diag*tic.'RO2') ;
  111. Po2o = Mo2 * TBT.'Rgo2' / vtotal ;
  112. Pinco = Pinco + Po2o ;
  113. Roinc = Roinc + tic.'RO2' ;
  114. Rocpinc= Rocpinc + (tic.'RO2'*Cpo2) ;
  115. Finsi ;
  116.  
  117. Mn2 = 0. ;
  118. Si TBT.'TN2' ;
  119. Mn2 = 'SOMT' (Diag*tic.'RN2') ;
  120. Pn2o = Mn2 * TBT.'Rgn2' / vtotal ;
  121. Pinco = Pinco + Pn2o ;
  122. Roinc = Roinc + tic.'RN2' ;
  123. Rocpinc= Rocpinc + (tic.'RN2'*Cpn2) ;
  124. Finsi ;
  125.  
  126. Mh2 = 0. ;
  127. Si TBT.'TH2' ;
  128. Mh2 = 'SOMT' (Diag*tic.'RH2') ;
  129. Ph2o = Mh2 * TBT.'Rgh2' / vtotal ;
  130. Pinco = Pinco + Ph2o ;
  131. Roinc = Roinc + tic.'RH2' ;
  132. Rocpinc= Rocpinc + (tic.'RH2'*Cph2) ;
  133. Finsi ;
  134.  
  135. Mhe = 0. ;
  136. Si TBT.'THE' ;
  137. Mhe = 'SOMT' (Diag*tic.'RHE') ;
  138. Pheo = Mhe * TBT.'Rghe' / vtotal ;
  139. Pinco = Pinco + Pheo ;
  140. Roinc = Roinc + tic.'RHE' ;
  141. Rocpinc= Rocpinc + (tic.'RHE'*Cphe) ;
  142. Finsi ;
  143.  
  144. Mco = 0. ;
  145. Si TBT.'TCO' ;
  146. Mco = 'SOMT' (Diag*tic.'RCO') ;
  147. Pcoo = Mco * TBT.'Rgco' / vtotal ;
  148. Pinco = Pinco + Pcoo ;
  149. Roinc = Roinc + tic.'RCO' ;
  150. Rocpinc= Rocpinc + (tic.'RCO'*Cpco) ;
  151. Finsi ;
  152.  
  153. Mco2 = 0. ;
  154. Si TBT.'TCO2' ;
  155. Mco2 = 'SOMT' (Diag*tic.'RCO2') ;
  156. Pco2o = Mco2 * TBT.'Rgco2' / vtotal ;
  157. Pinco = Pinco + Pco2o ;
  158. Roinc = Roinc + tic.'RCO2' ;
  159. Rocpinc= Rocpinc + (tic.'RCO2'*Cpco2) ;
  160. Finsi ;
  161.  
  162.  
  163. Pinc = Pinco * Tkfm ;
  164.  
  165. *--------------------------------------------------------------
  166. *A ce stade on a:
  167. * Mvap=Mrvp=Mvapb
  168. * Mtotg(n) = Mtotp(n) = Mtot(n+1)
  169. * Mtotp est déduite de Rhot déduit de la pression
  170. * Mtotg = Mvap + Minc
  171. * Mtot = est l'intégrale de Rho (dans laquelle on a remis le brouillard)
  172. * Pt=Pto # Ptn1 !!!
  173. *--------------------------------------------------------------
  174.  
  175. *--------------------------------------------------------------
  176. * On remet le brouillard en vapeur pour éventuellement
  177. * l'évaporer si les conditions sont réunies
  178. Rbro = tic.'RBRO' ;
  179. Rho = tic.'RHO' + Rbro ;
  180. Rvp = tic.'RVP' ;
  181. Rovp = Rvp ;
  182. **?? Mgini = somt (Diag*Rho) ;
  183. Rbro = 0. ;
  184. Rbrom = 0. ;
  185. Mliq = 0. ;
  186.  
  187. * Masses globales au pas de temps n
  188. Minc = Mair + Mo2 + Mn2 + Mh2 + Mhe + Mco + Mco2 ;
  189. Mtot = somt (Diag * tic.'RHO') ;
  190. Mrvp = somt (Diag * tic.'RVP') ;
  191. Mvap = Mrvp ;
  192. Mbro = 0. ;
  193. Mtotg= Mvap + Minc ;
  194.  
  195. mess 'ZZZ0 Masse(Rho) Mtot=' mtot ' Mvap+Minc Mtotg=' Mtotg
  196. ' Minc=' Minc ;
  197. mess 'ZZZ0 Mbro=' Mbro ' Mvap=' Mvap ' =? Mrvp=' Mrvp ;
  198.  
  199. dmt=(somt (Diag * Rho)) - mtot ;
  200. mtot = somt (Diag * Rho);
  201. mess '%$$ On remet le brouillard en suspension Mrhot='mtot
  202. ' dmt='dmt;
  203.  
  204. *-----------------------------------------------------------------------
  205. *-----------------------------------------------------------------------
  206. * Calcul de la pression 0D avec le champ Tf actualise
  207.  
  208. irtf = 'INVE' (rgp '*' (tic.'TF' '+' 273.15)) ;
  209. Pt = ((Rhomn-Rbrom) '*' Vtotal) '/' ('SOMT' (Diag '*' irtf));
  210. Ptc = kcht $mod scal sommet Pt ;
  211. rhot = Pt '*' ('INVE' (Rgp '*' (Tf '+' 273.15))) ;
  212. Pvap = Pt - Pinc ;
  213. Ptn1 = extr tic.'PT' ndl ;
  214. Mtotp= somt (Diag * rhot) ;
  215. mess 'ZZZ0 Mtotp est déduite de Rhot déduit de la pression' ;
  216. mess 'ZZZ0 Mtotg = Mvap + Minc ' ;
  217. mess 'ZZZ0 Mtot = est l intégrale de Rho (+ brouillard?)' ;
  218. mess 'ZZZ0 Mtotp=' Mtotp ' Mtotg=' Mtotg ' Mtot=' Mtot ;
  219. mess 'ZZZ0 Mbro=' Mbro ' Mvap=' Mvap 'Mrvp=' Mrvp ' Minc=' Minc;
  220.  
  221.  
  222. BILR=tic.'RHO' - tic.'RAIR' - tic.'RVP' ;
  223. BILT=rhot - tic.'RAIR' - tic.'RVP' ;
  224. mess 'ZZZ1 BILR='(mini BILR)(maxi BILR)' BILT='(mini BILT)(maxi BILT) ;
  225. mess 'ZZZ1 Rhomn= ' Rhomn ' Pto= ' Pto ' Ptn1=' Ptn1 ' Pt=' Pt;
  226. *-----------------------------------------------------------------------
  227. *-----------------------------------------------------------------------
  228.  
  229. *--------------------------------------------------------------
  230. * On Teste la saturation
  231. * Si non on sort
  232. PSTF = VARI 'PSATT' (Tf + 273.15) ;
  233. indsat = maxi (Pvap - PSTF - Tols) ;
  234. mess '1er RESP Pt=' Pt ' Rhomn=' Rhomn ' Mliq=' Mliq Mlq ;
  235.  
  236. Si(indsat < 0.) ;
  237. Mess 'ZZZ0 ******* On n est pas à saturation ** Pas de TPS ' ITER '**';
  238. Qcm = 0. ;
  239. Mliq = 0. ;
  240. Rebro= 0. ;
  241. Mlq = 0. ;
  242. yvap = Rovp * (inve Rhot) ;
  243. tic.'YVAP' = yvap ;
  244. Cpvap = TBT.'Cpvap' ;
  245. RESPROC Pt Rhot Qcm Rebro Mliq Mlq ;
  246. Mess '************ On quitte RETSAT tout de suite *******************';
  247. QUITTER RETSAT ;
  248. Finsi ;
  249.  
  250. Mess 'ZZZ2 ************* On est Sur saturé ** Pas de TPS ' ITER '****';
  251. *--------------------------------------------------------------
  252.  
  253. *-----------------------------------------------------------------------
  254. *-----------------------------------------------------------------------
  255. *2/ RETOUR A LA SATURATION
  256.  
  257. Tg = kcht $mod scal sommet Tfm ;
  258. Tg = Tf ;
  259. Tgm = ('SOMT' (Diag * tic.'TF'))/Vtotal ;
  260. Tkgm = Tgm + 273.15 ;
  261.  
  262. nbzv=10O ;
  263. wz = TBT.'V1alfa' ;
  264. wz=1.;
  265. wb=TBT.'TAUXcm' ;
  266. *wb=0.01 ;
  267. *wb=1. ;
  268.  
  269. REPETER BLOC nbzv ;
  270.  
  271. Mtotg=Mvap + Mair + Mo2 + Mn2 + Mh2 + Mhe + Mco + Mco2 ;
  272. Pinc = Pinco * Tkgm ;
  273. Pvap = Pt - Pinc ;
  274.  
  275. *--------------------------------------------------------------
  276. * CALCUL DE Rhosat On borne Rho par la saturation
  277. * On calcule d'abord Rhosat que l'on déduit de Psat
  278. * Si Gaz parfait Rhosat = Psat/(Rgvap (Tg + 273.15))
  279. * Sinon Rhosat = Psat/(Zv* Rgvap*(Tg + 273.15))
  280. * où Zv = f (Rhosat , Tgk) Il faut donc itérer
  281.  
  282. Si('EGA' LOIETAT 'GR') ;
  283. Zv = (VARI 'ZVAP' rstf (Tg + 273.15)) 'ABS' ;
  284. Sinon ;
  285. Zv = 1. ;
  286. Finsi ;
  287. * RSTF : Ro vapeur à saturation
  288. PSTF = VARI 'PSATT' (Tg + 273.15) ;
  289. RSTF = PSTF * (INVE (Zv * TBT.'Rgvap' * (Tg + 273.15))) ;
  290. Rstfm= ('SOMT' (Diag '*' rstf))/Vtotal ;
  291.  
  292. * ROVP : Ro vapeur borné par la saturation
  293. ROVP = 'KOPS' ROVP '>|' RSTF ;
  294. ROVP = 'KOPS' ROVP '|<' 1.e-10 ;
  295. *--------------------------------------------------------------
  296. *--------------------------------------------------------------
  297.  
  298.  
  299. *--------------------------------------------------------------
  300. *--------------------------------------------------------------
  301. * CALCUL de RBRO
  302. * On sépare vapeur + brouillard (tic.'RVP') et liquide
  303. * (tic.'RLIQ')
  304.  
  305. RDIF = (tic.'RVP' - ROVP) + 1.e-10 ;
  306. RDIF = wb*RDIF ;
  307. ROVP = (tic.'RVP' - RDIF) + 1.e-10 ;
  308. RVP = ROVP + RDIF ;
  309. YBRO = RDIF * (INVE rhot) ;
  310. YBRO = 'KOPS' YBRO '>|' TBT.'YBRO' ;
  311. RBRO = YBRO '*' Rhot ;
  312. RLIQ = (RDIF - (YBRO*rhot)) ;
  313. mess ' RLIQ=' (maxi RLIQ) (mini RLIQ) ;
  314.  
  315. * Masses globales au pas de temps n iteration &bloc
  316.  
  317. Minc = Mair + Mo2 + Mn2 + Mh2 + Mhe + Mco + Mco2 ;
  318. Mrvp = somt (Diag * tic.'RVP') ;
  319. Mbro = 0. ;
  320. Mliq = 'SOMT' (Diag '*' RLIQ) ;
  321. Mdif = 'SOMT' (Diag '*' RDIF) ;
  322. Mtot = somt (Diag * tic.'RHO') ;
  323. Mtotg= Mvap + Minc ;
  324. Mtotp= somt (Diag * rhot) ;
  325. Mtotl= Mvap + Minc + Mliq + Mbro ;
  326.  
  327. * On remet la chaleur latente de condensation dans le mélange
  328. * gazeux. On considère que le brouillard ou l'eau liquide se
  329. * mettent rapidement en équilibre thermique avec le gaz.
  330.  
  331. Cpmnm = 'EXTR' tic.'Cpm' ndl ;
  332. Cvmnm = 'EXTR' tic.'Cvm' ndl ;
  333.  
  334. Hls ='VARI' 'HLS' PTC (Tg + 273.15) ;
  335. Cpliq='VARI' 'CPLIQ' PTC Hls ;
  336.  
  337. dtg = (RLIQ*Cpliq)*(Tg+273.15)/((ROVP*Cpvap) + Rocpinc) ;
  338. **?? Tg = Tg - Dtg ;
  339.  
  340. mess ' DTG=' (mini DTG) (maxi dtg);
  341.  
  342.  
  343. Si ((Maxi Tg) > 300.);
  344. mess ' Attention La température dépasse 300°C : On sort des tables';
  345. * La température est limitée à 300°C pour rester dans les tables !!
  346. a = 300.; al=0.99; ala=al*a ;
  347. * b=ala*al*((1.-al)**(-1.));
  348. * b2=ala*(1.+ (ala*al) - ala)*((1.-al)**(-1.));
  349. b3 = ala*(1.+ (ala*ala*al) - (ala*ala))*((1.-al)**(-1.));
  350. ik=masq Tg 'INFERIEUR' ala ;
  351. Tg=(ik*Tg)+((1.-ik)*a*((Tg*Tg*Tg+b3)*(inve(a+(Tg*Tg*Tg+b3)))));
  352. mess '%$$ Mini/Maxi Tg' (mini Tg) (Maxi Tg) ;
  353. Finsi;
  354.  
  355. Tgm = ('SOMT' (Diag * tg))/Vtotal ;
  356. Tkgm = Tgm + 273.15 ;
  357.  
  358. * On recalcule la pression
  359. Si('EGA' LOIETAT 'GR') ;
  360. Zv = (VARI 'ZVAP' ROVP (Tg + 273.15)) 'ABS' ;
  361. Sinon ;
  362. Zv = 1. ;
  363. Finsi ;
  364.  
  365. sr = (Zv*ROVP*TBT.'Rgvap'*(Tg+273.15)) + Pinc ;
  366. sr = 'SOMT' (sr * Diag) ;
  367. a = sr *(1./Vtotal) ;
  368. da = (Pt - a)/a 'ABS' ;
  369. Pt = (wz*a) + ((1.-wz)*Pt) ;
  370. Ptc= kcht $mod scal sommet Pt ;
  371. rhot = Pt '*' ('INVE' (Rgp '*' (Tg '+' 273.15))) ;
  372.  
  373. mess 'Itér ' &bloc ' da=' da ' Pt=' Pt (mini Tg) (maxi Tg) ;
  374. Si(da < 1.e-3);Mess 'On a convergé sur la pression en' &bloc
  375. 'itérations ZZZ ' ITER ; quitter BLOC; finsi ;
  376.  
  377. FIN BLOC ;
  378.  
  379. Si(da > 1.e-3) ;
  380. Mess 'On n a pas convergé sur la pression en' &bloc
  381. 'itérations, da=' da ' ZZZ ' ITER ;
  382. finsi ;
  383.  
  384. mess 'ZZZ2 Iter=' &bloc ;
  385. mess 'ZZZ2 Mtotp est déduite de Rhot déduit de la pression' ;
  386. mess 'ZZZ2 Mtotg = Mvap + Minc ' ;
  387. mess 'ZZZ2 Mtot = est l intégrale de Rho (+ brouillard?)' ;
  388. mess 'ZZZ2 Mtotl = est la masse liquide condensée en masse' ;
  389. mess 'ZZZ2' ;
  390. mess 'ZZZ2 Mtotp=' mtotp ' Mtotg=Mvap+Minc=' Mtotg
  391. ' Mtot=' Mtot ' Mtotl=' Mtotl (Mtotg+Mliq+Mbro) ;
  392. mess 'ZZZ2 Mbro=' Mbro ' Mvap=' Mvap 'Mrvp=' Mrvp ' Minc=' Minc;
  393. mess 'ZZZ2 Mliq=' Mliq ' Mdif=' Mdif ' Minc=' Minc ;
  394.  
  395. tic.'TF' = Tg ;
  396. tic.'RVP' = Rvp ;
  397. tic.'RBRO' = Rbro ;
  398. tic.'RLIQ' = Rliq ;
  399. *--------------------------------------------------------------
  400. *--------------------------------------------------------------
  401.  
  402. * On s'occupe des pressions partielles
  403. Pvap = (Zv * ROVP * TBT.'Rgvap' * Diag * (Tf + 273.15)) ;
  404. Pvap = ('SOMT' Pvap) * (1./Vtotal) ;
  405. Pinc = Pinco * Tkfm ;
  406.  
  407. Mliq = 'SOMT' (Diag*Rliq) ;
  408. Mess 'Mgini=' Mgini ' rotfinal=' (somt(Diag*Rhot)) ' Mliq2=' mliq
  409. 'DT=' tic.'DT' ;
  410. Rliqm = Mliq / Vtotal ;
  411.  
  412. * On est obligé de reprendre la densité moyenne Rhom, la densité moyenne
  413. * de vapeur Rhomvn et la densité moyenne vapeur + gouttes et lui enlever
  414. * le drainage du liquide sans oublier bien sur l'énergie
  415.  
  416. **?? Rhomn = ( somt (Diag * Rhot)) / Vtotal ;
  417. Rhomn = ('EXTR' (tic.'Rhomn') (dime tic.'Rhomn'))
  418. -(( somt (Diag * Rdif))/Vtotal) ;
  419. tic.'Rhomn' = ('ENLEV' (tic.'Rhomn') (dime tic.'Rhomn'))
  420. 'ET' ('PROG' Rhomn) ;
  421.  
  422. **?? Rhomvn = ( somt (Diag * Rovp)) / Vtotal ;
  423. Rhomvn = ('EXTR' (tic.'Rhomv') (dime tic.'Rhomv'))
  424. -(( somt (Diag * Rdif))/Vtotal) ;
  425. tic.'Rhomv' = ('ENLEV' (tic.'Rhomv') (dime tic.'Rhomv'))
  426. 'ET' ('PROG' Rhomvn) ;
  427.  
  428. **?? Rhomvg = ( somt (Diag * (Rovp+Rbro))) / Vtotal ;
  429. Rhomvg = ('EXTR' (tic.'Rhomvg') (dime tic.'Rhomvg'))
  430. -(( somt (Diag * (Rdif+Rbro)))/Vtotal) ;
  431. tic.'Rhomvg' = ('ENLEV' (tic.'Rhomvg') (dime tic.'Rhomvg'))
  432. 'ET' ('PROG' Rhomvg) ;
  433.  
  434. *------------------------
  435. **??reliq = Cpmnm '*' ( 'SOMT' (Diag '*' RLIQ '*' (Tf '+' 273.15))) ;
  436. **??reliq = reliq + ('SOMT' (Lvv*Diag*tic.'RLIQ')) ;
  437. **?? reliq = somt (Diag*reliq) ;
  438. rebro = 0. ;
  439. *------------------------
  440.  
  441. **?? Remncm = 'SOMT' (Cpliq '*' Diag '*' Rliq '*' (Tf '+' 273.15)) ;
  442. Remncm = ('EXTR' tic.'Remn' (dime tic.'Remn'))
  443. -('SOMT' (Cpmnm '*' Diag '*' Rliq '*' (Tf '+' 273.15))) ;
  444.  
  445. Remn = ('EXTR' tic.'Remn' (dime tic.'Remn')) ;
  446. mess 'REmn =' Remn ' Remncm=' Remncm ;
  447. **?? tic.'Remn' = ('ENLEV' (tic.'Remn') (dime tic.'Remn'))
  448. **?? 'ET' ('PROG' Remncm) ;
  449.  
  450. Mbro = somt (Diag * RBRO) ;
  451. **?? Mvp = somt (Diag * ROVP) ;
  452. Rbrom = Mbro / Vtotal ;
  453.  
  454.  
  455. Cpvap = TBT.'Cpvap' ;
  456. Cvvap = Cpvap '-' TBT.'Rgvap' ;
  457. Qcm = Mliq / (tic.'DT') ;
  458. tic.'Qcm'= tic.'Qcm' 'ET' ('PROG' Qcm) ;
  459. *Econdm = Cvvap * (SOMT (Diag*RLIQ*(tic.'TF' + 273.15)))/tic.'DT' ;
  460. *Hcondm = Cpvap * (SOMT (Diag*RLIQ*(tic.'TF' + 273.15)))/tic.'DT' ;
  461.  
  462. * Détermination de la masse volumique moyenne de liquide en suspension :
  463. tic.'Rliqm' = tic.'Rliqm' et ('PROG' Rliqm) ;
  464. tic.'Rbrom' = tic.'Rbrom' et ('PROG' Rbrom) ;
  465. tic.'Mcondm'= tic.'Mcondm' 'ET' ('PROG' Mliq) ;
  466.  
  467.  
  468. mtot = somt (Diag * Rhot);
  469. **?? mair = somt (Diag * Rair);
  470. mrvp = somt (Diag * tic.'RVP');
  471. mvap = somt (Diag * (tic.'RVP' - tic.'RBRO')) ;
  472. mbro = somt (Diag * tic.'RBRO') ;
  473. mdif = somt (Diag * RDIF) ;
  474.  
  475. mess '%$$ RETSAT Masse totale (Rhot) APRES ' mtot
  476. 'Mair=' mair 'Mvap =' mvap ;
  477. mess '%$$ RETSAT APRES Mbrot=' Mbro ' Mrvp=' Mrvp ' Mdif='mdif;
  478.  
  479. mess '%$$ RETSAT Pression ' Pt ' Reliq=' reliq
  480. ' Flux condensés en Masse' Qcm 'Kg/s' ;
  481.  
  482. *-----------------------------------------------------------------------
  483. tic.'RHO' = Rhot ;
  484. tic.'RBRO' = Rbro ;
  485. tic.'RVP' = Rovp + Rbro ;
  486. tic.'ROVP' = ELNO (NOEL Rvp $mod) $mod ;
  487. yvap = Rovp * (inve Rhot) ;
  488. tic.'YVAP' = yvap ;
  489. tic.'YAIR' = tic.'RAIR' * (inve Rhot) ;
  490.  
  491. Si TBT.'TO2' ;
  492. tic.'YO2' = tic.'RO2' * (inve Rhot) ;
  493. Finsi ;
  494. Si TBT.'TN2' ;
  495. tic.'YN2' = tic.'RN2' * (inve Rhot) ;
  496. Finsi ;
  497. Si TBT.'TH2' ;
  498. tic.'YH2' = tic.'RH2' * (inve Rhot) ;
  499. Finsi ;
  500. Si TBT.'THE' ;
  501. tic.'YHE' = tic.'RHE' * (inve Rhot) ;
  502. Finsi ;
  503. Si TBT.'TCO' ;
  504. tic.'YCO' = tic.'RCO' * (inve Rhot) ;
  505. Finsi ;
  506. Si TBT.'TCO2' ;
  507. tic.'YCO2' = tic.'RCO2' * (inve Rhot) ;
  508. Finsi ;
  509.  
  510. * mess '%$$ Somme des Y ' (mini sy)(maxi sy) ;
  511. mess '...............................................................';
  512. mess ;
  513. irtf = 'INVE' (rgp '*' (tic.'TF' '+' 273.15)) ;
  514. Pt = ((Rhomn-Rbrom) '*' Vtotal) '/' ('SOMT' (Diag '*' irtf)) ;
  515. Ptc= kcht $mod scal sommet Pt ;
  516. rhot = Pt '*' ('INVE' (Rgp '*' (Tf '+' 273.15))) ;
  517. mess '2eme RESP Pt=' Pt ' Rhomn=' Rhomn ' Mliq=' Mliq Mlq ;
  518.  
  519. RESPROC Pt Rhot Qcm Rebro Mliq Mlq ;
  520.  
  521. *-----------------------------------------------------------------------
  522. * Fin du traitement condensation en masse
  523. *-----------------------------------------------------------------------
  524.  
  525. FINPROC;
  526.  
  527.  
  528.  

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