Télécharger fla010.eso

Retour à la liste

Numérotation des lignes :

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

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