Télécharger fla019.eso

Retour à la liste

Numérotation des lignes :

  1. C FLA019 SOURCE TTMF3 12/07/05 21:15:42 7425
  2. SUBROUTINE FLA019(T,TC,MP,TS,XS,IDEB,TRAV)
  3. C---------------------------------------------------------------------
  4. C Calcul du débit et des conditions en sortie de recombineur
  5. C (concentrations et température)
  6. C---------------------------------------------------------------------
  7. C
  8. C---------------------------
  9. C Parametres Entree/Sortie :
  10. C---------------------------
  11. C
  12. C E/ T : flottant : Temps courant (s)
  13. C E/ TC : flottant : Température des plaques (K)
  14. C /S MP : flottant : Debit total au travers du PAR (kg/s)
  15. C /S TS : flottant : Température du gaz à la sortie du PAR (K)
  16. C /S XS : flottant : Fractions molaires à la sortie du PAR
  17. C /S IDEB : entier : Flag identifiant la situation
  18. C (négatif si le débit est négatif !)
  19. C
  20. C------------------------------
  21. C Variables de TRAV utilisées :
  22. C------------------------------
  23. C
  24. C E/ PRESSION : Pression à l'entrée du PAR (Pa)
  25. C E/ TEMPENT : Température à l'entrée du PAR (K)
  26. C E/ XiMOY : Fraction molaire à l'entrée du PAR
  27. C E/ M(7) : Masse molaire des constituants du mélange (kg/mol)
  28. C E/ EPS_MH2 : Seuil en deça duquel on suppose que la consommation
  29. C d'H2 max est mise à 0. (kg/s)
  30. C E/ EPS_DT : Ecart entre la température des plaques et celle du
  31. C gaz (Tc-T) en deça duquel on suppose qu'il n'y a pas
  32. C d'échange convectif (K)
  33. C E/ U : Vitesse minimale dans le PAR (= 0.01 m/s)
  34. C E/ S : Surface des plaques (m2)
  35. C E/ SP : Section de passage entrée/sortie fluide du PAR (m2)
  36. C E/ MC : Masse des plaques (kg)
  37. C E/ CPC : Chaleur spécifique des plaques (J/kg/K)
  38. C
  39. C-------------------
  40. C Autres variables :
  41. C-------------------
  42. C
  43. C MH2MAX : Consommation maximale par les plaques catalytique (kg/s)
  44. C MH2 : Consommation d'H2 par les plaques pour le débit MP (kg/s)
  45. C MH2EMAX : Débit d'hydrogène disponible à l'entrée du PAR (kg/s)
  46. C MO2EMAX : Débit d'hydrogène disponible à l'entrée du PAR limité par
  47. C le débit d'O2 insuffisant (kg/s)
  48. C
  49. C---------------------------------------------------------------------
  50. C
  51. C Langage : ESOPE + FORTRAN 77
  52. C
  53. C Mise en oeuvre : H. Paillère (1997, TTMF)
  54. C
  55. C---------------------------------------------------------------------
  56. IMPLICIT INTEGER(I-N)
  57. IMPLICIT REAL*8 (A-H,O-Z)
  58. REAL*8 XE(7),OME(7),XS(7),OMS(7)
  59. REAL*8 MH2MAX,MH2,MP,MH2EMAX,MO2EMAX
  60. segment trav
  61. integer iKALP,iMODEL
  62. real*8 e,L,Lch,Dh,S,sp,Ck
  63. real*8 mc,Cpc
  64. real*8 g,R,deltah
  65. real*8 M(nbesp),cstmod(ncst)
  66. real*8 Cpi(nbesp),al
  67. real*8 eps_mh2,eps_dt,eps_con,u
  68. real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  69. real*8 XHEMOY,XCO2MOY,XCOMOY
  70. endsegment
  71. C
  72. C OME(7) : Fraction massique
  73. PE = PRESSION
  74. TE = TEMPENT
  75. XE(1) = XN2MOY
  76. XE(2) = XO2MOY
  77. XE(3) = XH2MOY
  78. XE(4) = XH2OMOY
  79. XE(5) = XHEMOY
  80. XE(6) = XCO2MOY
  81. XE(7) = XCOMOY
  82. CALL FLA012(XE,OME,TRAV)
  83. ROE = FLA001(PE,TE,XE,TRAV)
  84. CPE = FLA003(OME,TRAV)
  85. C
  86. C Calcul de dmh2/dt sans limitation par le débit d'H2
  87. C et/ou d'O2 à l'entrée du PAR (kg/s)
  88. IF (IMODEL .EQ. 1) THEN
  89. MH2MAX = FLA008(PE,XE,CSTMOD(1),CSTMOD(2),CSTMOD(4),CSTMOD(5))
  90. ELSEIF (IMODEL .EQ. 2) THEN
  91. MH2MAX = FLA008(PE,XE,CSTMOD(1),CSTMOD(2),CSTMOD(4),CSTMOD(5))
  92. ENDIF
  93. C
  94. C 1. Cas sans comsommation d'hydrogène
  95. IF (MH2MAX.LT.EPS_MH2) THEN
  96. IF ((TC-TE) .LT. EPS_DT) THEN
  97. C 1.1 et sans echange convectif
  98. MH2 = 0.D0
  99. MP = U*ROE*SP
  100. TC = TE
  101. IDEB = 11
  102. ELSE
  103. C 1.2 et avec echange convectif (bilan d'énergie des plaques)
  104. MH2 = 0.D0
  105. MP = FLA015(T,TC,MH2,TRAV)
  106. IDEB = 12
  107. ENDIF
  108. ELSE
  109. C 2. Cas avec comsommation d'hydrogène
  110. IF ((TC-TE) .LT. EPS_DT) THEN
  111. C 2.1 et sans echange convectif mais prise en compte d'un déficit
  112. C éventuel en H2/O2 et débit minimum correspondant à la vitesse U
  113. BOBO = MH2MAX
  114. MH2EMAX = U * ROE * SP * OME(3)
  115. MO2EMAX = U * ROE * SP * 2.D0*M(3)*OME(2)/M(2)
  116. IF (MH2EMAX .LT. BOBO) THEN
  117. BOBO = MH2EMAX
  118. ENDIF
  119. IF (MO2EMAX .LT. BOBO) THEN
  120. BOBO = MO2EMAX
  121. ENDIF
  122. MH2 = BOBO
  123. MP = U*ROE*SP
  124. IDEB = 21
  125. ELSE
  126. C 2.2 et avec echange convectif (bilan d'énergie des plaques)
  127. C et prise en compte d'un déficit éventuel en H2/O2 par rapport
  128. C au débit mp traversant le PAR et correspondant à MH2MAX
  129. MP = FLA015(T,TC,MH2MAX,TRAV)
  130. C 2.2.1 Déficit en O2 (indic=2) ou déficit en H2 (indic=3)
  131. BNERF = OME(3)
  132. IF (BNERF .GT. 2.D0*M(3)*OME(2)/M(2)) THEN
  133. INDIC = 2
  134. BNERF = 2.D0*M(3)*OME(2)/M(2)
  135. ENDIF
  136. IF (MP*BNERF .LT. MH2MAX) THEN
  137. MP = FLA014(T,TC,INDIC,TRAV)
  138. MH2 = MP*BNERF
  139. IDEB = 221
  140. ELSE
  141. C 2.2.2 Pas de déficit en O2/H2
  142. MH2 = MH2MAX
  143. IDEB = 222
  144. ENDIF
  145. ENDIF
  146. ENDIF
  147. C
  148. C Température en sortie du PAR (bilan d'énergie du gaz)
  149. IF ((TC-TE).LT.EPS_DT) THEN
  150. TS = TE
  151. ELSE
  152. TS = TE + S*FLA009(PE,TE,TC,TE,XE,TRAV)*(TC-TE)
  153. & / (MP*CPE - AL*MH2)
  154. ENDIF
  155. C
  156. C Composition du mélange en sortie
  157. OMS(1) = OME(1)
  158. OMS(2) = OME(2) - M(2)*MH2/(2.D0*M(3)*MP)
  159. OMS(3) = OME(3) - MH2/MP
  160. OMS(4) = OME(4) + M(4)*MH2/(M(3)*MP)
  161. OMS(5) = OME(5)
  162. OMS(6) = OME(6)
  163. OMS(7) = OME(7)
  164. CALL FLA013(XS,OMS,TRAV)
  165. C
  166. C Conditions de réalisabilité
  167. XS(1) = ABS(XS(1))
  168. XS(2) = ABS(XS(2))
  169. XS(3) = ABS(XS(3))
  170. XS(5) = ABS(XS(5))
  171. XS(6) = ABS(XS(6))
  172. XS(7) = ABS(XS(7))
  173. XS(4) = 1.D0 - XS(1) - XS(2) - XS(3) - XS(5) - XS(6) - XS(7)
  174. C
  175. RETURN
  176. END
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  

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