Télécharger la plaquette Cast3M


fr en

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

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