Télécharger fla015.eso

Retour à la liste

Numérotation des lignes :

  1. C FLA015 SOURCE TTMF3 12/07/05 21:15:37 7425
  2. FUNCTION FLA015(T,TC,MH2,TRAV)
  3. C---------------------------------------------------------------------
  4. C Calcul du débit au travers du PAR pour une température de plaque de
  5. C Tc et une consommation d'hydrogène mh2=mh2max
  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 E/ MH2 : flottant : Hydrogène consommé (kg/s)
  15. C /S FLA015 : flottant : Débit total traversant le PAR (kg/s)
  16. C
  17. C------------------------------
  18. C Variables de TRAV utilisées :
  19. C------------------------------
  20. C
  21. C E/ PRESSION : Pression à l'entrée du PAR (Pa)
  22. C E/ TEMPENT : Température à l'entrée du PAR (K)
  23. C E/ XiMOY : Fraction molaire à l'entrée du PAR
  24. C E/ M(7) : Masse molaire des constituants du mélange (kg/mol)
  25. C E/ U : Vitesse minimale dans le PAR (= 0.01 m/s)
  26. C E/ S : Surface des plaques (m2)
  27. C E/ SP : Section de passage entrée/sortie fluide du PAR (m2)
  28. C E/ L : Hauteur des plaques (m)
  29. C E/ LCH : Hauteur de la cheminée (m)
  30. C E/ DH : Diamètre hydraulique entre les plaques (m)
  31. C E/ CK : Demi somme des coefficients de perte de charge à
  32. C l'entrée et à la sortie du PAR
  33. C E/ G : Gravité (=9.81 m/s2)
  34. C E/ AL : Cste fonction de CP(i)
  35. C E/ EPS_CON : Seuil de convergence pour Newton (|f(x)|<EPS_CON)
  36. C
  37. C--------------------
  38. C Algorithme utilisé
  39. C--------------------
  40. C
  41. C L'équation de Bernouilli (qdm dans le PAR) se met sous la forme
  42. C F(mp,ros) = a mp^2 + b mp + c (ro^2 - ros^2) = 0
  43. C ros(mp) = ro / (1 + D/(mp-E))/ (1 - F/mp)
  44. C On recherche la plus grande racine de ce polynome par une méthode
  45. C de dichotomie (c'est pourquoi on privilégie le coté Xmax).
  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. *
  55. IMPLICIT INTEGER(I-N)
  56. IMPLICIT REAL*8 (A-H,O-Z)
  57. REAL*8 FLA015
  58. REAL*8 XE(7),OME(7)
  59. REAL*8 MUE,ME,MH2,MP
  60. -INC CCOPTIO
  61. segment trav
  62. integer iKALP,iMODEL
  63. real*8 e,L,Lch,Dh,S,sp,Ck
  64. real*8 mc,Cpc
  65. real*8 g,R,deltah
  66. real*8 M(nbesp),cstmod(ncst)
  67. real*8 Cpi(nbesp),al
  68. real*8 eps_mh2,eps_dt,eps_con,u
  69. real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  70. real*8 XHEMOY,XCO2MOY,XCOMOY
  71. endsegment
  72. C
  73. MAXNEWT = 10000
  74. COFNEWT = 1.D0
  75. C
  76. PE = PRESSION
  77. TE = TEMPENT
  78. XE(1) = XN2MOY
  79. XE(2) = XO2MOY
  80. XE(3) = XH2MOY
  81. XE(4) = XH2OMOY
  82. XE(5) = XHEMOY
  83. XE(6) = XCO2MOY
  84. XE(7) = XCOMOY
  85. CALL FLA012(XE,OME,TRAV)
  86. ROE = FLA001(PE,TE,XE,TRAV)
  87. MUE = FLA004(TE,XE,TRAV)
  88. CPE = FLA003(OME,TRAV)
  89. ibou = cpi(/1)
  90. sum = 0.D0
  91. do 5 i=1,ibou
  92. sum = sum + xe(i)*m(i)
  93. 5 continue
  94. ME = sum
  95. C
  96. ROEXT = ROE
  97. C
  98. C Coefficient de l'équation de Bernouilli et de la densité en sortie
  99. A = Ck*COFNEWT*COFNEWT
  100. B = COFNEWT*48.D0*mue*sp*L/(Dh*Dh)
  101. C = -0.5D0*sp*sp*g*(L/2.D0+Lch)
  102. DD = S*fla009(pe,Te,Tc,Te,xe,trav)*(Tc-Te)/(Te*Cpe)*COFNEWT
  103. EE = al*mh2/Cpe*COFNEWT
  104. FF = mh2*Me/(2.D0*M(3))*COFNEWT
  105. C
  106. C Méthode de Dichotomie : la plus grande des racines est au-delà
  107. C des poles du denominateur (EE-DD) et FF
  108. I0 = 0
  109. XMIN = MAX(EE-DD,FF)
  110. IF (XMIN .LE. 0.D0) THEN
  111. XMIN = 0.D0
  112. ENDIF
  113. XMAX = MAX(1.D0,U*ROE*SP*500.D0)/COFNEWT
  114. RESUP = FLA016(XMAX,roext,A,B,C,DD,EE,FF)
  115. RESU = RESUP
  116. RESUM = FLA016(XMIN,roext,A,B,C,DD,EE,FF)
  117. EPS_2 = MIN(ABS(RESUM),EPS_CON)
  118. C
  119. C Recherche de la racine sur l'intervalle (XMIN,XMAX)
  120. 10 CONTINUE
  121. X = (XMAX + XMIN) / 2.D0
  122. IF (I0 .GE. MAXNEWT) GOTO 20
  123. IF ((XMAX-XMIN) .LT. EPS_CON) GOTO 30
  124. C IF (ABS(RESU) .LT. (EPS_CON*EPS_2)) GOTO 40
  125. CC IF (ABS(RESU) .LT. (RESUP*EPS_CON*EPS_2)) GOTO 40
  126. I0 = I0 + 1
  127. RESU = FLA016(x,roext,A,B,C,DD,EE,FF)
  128. IF ((RESU*RESUP) .LE. 0.D0) THEN
  129. XMIN = X
  130. ELSE
  131. XMAX = X
  132. ENDIF
  133. GOTO 10
  134. 20 CONTINUE
  135. C
  136. C I0 .GE. MAXNEWT : non convergence du Newton
  137. C Se produit à l'usage lorsque la racine et les poles sont quasi nuls.
  138. C Par suite, le débit est mis à zero.
  139. INTERR(1) = MAXNEWT
  140. CALL ERREUR(151)
  141. X = 0.D0
  142. GOTO 40
  143. 30 CONTINUE
  144. C
  145. C (XMAX-XMIN) .LT. EPS_CON : les 2 bornes de dichotomie sont égales.
  146. C Si la fonction est de même signe en XMAX et XMIN, il n'y a pas de
  147. C racine réelle au dela des poles. Par suite, le débit est nul.
  148. RESUP = FLA016(XMAX,roext,A,B,C,DD,EE,FF)
  149. RESUM = FLA016(XMIN,roext,A,B,C,DD,EE,FF)
  150. IF ((RESUM*RESUP) .GT. 0.D0) THEN
  151. X = 0.D0
  152. ENDIF
  153. 40 CONTINUE
  154. C
  155. C On maintient quoiqu'il arrive un débit résiduel dans le PAR
  156. C de vitesse débitante u
  157. MP = X * COFNEWT
  158. IF (U*ROE*SP .GT. MP) THEN
  159. MP = U*ROE*SP
  160. ENDIF
  161. FLA015 = MP
  162. C
  163. RETURN
  164. END
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  

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