Télécharger hea010.eso

Retour à la liste

Numérotation des lignes :

  1. C HEA010 SOURCE SERRE 12/07/10 21:15:18 7431
  2. SUBROUTINE HEA010(T,TC,HEAS,HEAF,DTCSDT,MP,TRAV)
  3. C---------------------------------------------------------------------
  4. C Calcul du coefficient 'ki' de RK4 (voir HEA011)
  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 HEAS : flottant : Puissance générée dans le solide (J/s)
  15. C /S HEAF : flottant : Puissance transmise au fluide (J/s)
  16. C /S DTCSDT : flottant : dTc/dt (K/s)
  17. C /S MP : flottant : Débit de mélange traversant le HEATER (kg/s)
  18. C E/ TRAV : entier : Pointeur du segment de travail
  19. C
  20. C------------------------------
  21. C Variables de TRAV utilisées :
  22. C------------------------------
  23. C
  24. C E/ PRESSION : Pression à l'entrée du HEATER (Pa)
  25. C E/ TEMPENT : Température à l'entrée du HEATER (K)
  26. C E/ XiMOY : Fraction molaire à l'entrée du HEATER
  27. C E/ EPS_MH2 : Seuil en deça duquel on suppose que la consommation
  28. C d'H2 max est mise à 0. (kg/s)
  29. C E/ EPS_DT : Ecart entre la température des plaques et celle du
  30. C gaz (Tc-T) en deça duquel on suppose qu'il n'y a pas
  31. C d'échange convectif (K)
  32. C E/ U : Vitesse minimale dans le HEATER (= 0.01 m/s)
  33. C E/ S : Surface des plaques (m2)
  34. C E/ SP : Section de passage entrée/sortie fluide du HEATER (m2)
  35. C E/ MC : Masse des plaques (kg)
  36. C E/ CPC : Chaleur spécifique des plaques (J/kg/K)
  37. C
  38. C-------------------
  39. C Autres variables :
  40. C-------------------
  41. C
  42. C MH2 : Consommation fictive d'H2 par les plaques associée à la
  43. C puissance thermique délivrée (kg/s)
  44. C
  45. C---------------------------------------------------------------------
  46. C
  47. C Langage : ESOPE + FORTRAN 77
  48. C
  49. C Mise en oeuvre : H. Paillère (1997, TTMF)
  50. C
  51. C---------------------------------------------------------------------
  52. IMPLICIT INTEGER(I-N)
  53. IMPLICIT REAL*8 (A-H,O-Z)
  54. REAL*8 XE(7)
  55. REAL*8 MH2,MP
  56. segment trav
  57. integer iKALP,iMODEL
  58. real*8 e,L,Lch,Dh,S,sp,Ck
  59. real*8 mc,Cpc
  60. real*8 g,R,deltah
  61. real*8 M(nbesp),cstmod(ncst)
  62. real*8 Cpi(nbesp),al
  63. real*8 eps_mh2,eps_dt,eps_con,u
  64. real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  65. real*8 XHEMOY,XCO2MOY,XCOMOY
  66. endsegment
  67. C
  68. C Densité du mélange
  69. PE = PRESSION
  70. TE = TEMPENT
  71. XE(1) = XN2MOY
  72. XE(2) = XO2MOY
  73. XE(3) = XH2MOY
  74. XE(4) = XH2OMOY
  75. XE(5) = XHEMOY
  76. XE(6) = XCO2MOY
  77. XE(7) = XCOMOY
  78. ROE = FLA001(PE,TE,XE,TRAV)
  79. C
  80. C Débit minimum dans le PAR thermique (valeur par défaut)
  81. MP = U*ROE*SP
  82. C
  83. C dmh2/dt équivalent à la puissance thermique du HEATER
  84. MH2 = CSTMOD(1) / deltah
  85. C
  86. C 1. Cas sans comsommation d'hydrogène (puissance thermique résiduelle)
  87. IF (MH2 .LT. EPS_MH2) THEN
  88. IF ((TC-TE) .LT. EPS_DT) THEN
  89. C 1.1 et sans echange convectif
  90. HEAS = 0.D0
  91. HEAF = 0.D0
  92. ELSE
  93. C 1.2 et avec echange convectif (bilan d'énergie des plaques)
  94. HEAS = 0.D0
  95. HEAF = -S*FLA009(PE,TE,TC,TE,XE,TRAV)*(TC-TE)
  96. ENDIF
  97. ELSE
  98. C 2. Cas avec comsommation d'hydrogène
  99. IF ((TC-TE) .LT. EPS_DT) THEN
  100. C 2.1 et sans echange convectif
  101. HEAS = DELTAH * MH2
  102. HEAF = 0.D0
  103. ELSE
  104. C 2.2 et avec echange convectif (bilan d'énergie des plaques)
  105. MP = FLA015(T,TC,MH2,TRAV)
  106. HEAS = DELTAH * MH2
  107. HEAF = -S*FLA009(PE,TE,TC,TE,XE,TRAV)*(TC-TE)
  108. ENDIF
  109. ENDIF
  110. HEAS = HEAS * cstmod(2)
  111. DTCSDT = (HEAS + HEAF) / (MC*CPC)
  112. C
  113. RETURN
  114. END
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  

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