Télécharger fla009.eso

Retour à la liste

Numérotation des lignes :

  1. C FLA009 SOURCE SERRE 12/07/10 21:15:14 7431
  2. FUNCTION FLA009(P,T,TC,TPROP,X,TRAV)
  3. C---------------------------------------------------------------------
  4. C Calcul du coefficient d'échange thermique entre les plaques du PAR
  5. C et le gaz via la corrrélation d'Elenbaas (en W/m2/K)
  6. C---------------------------------------------------------------------
  7. C
  8. C---------------------------
  9. C Parametres Entree/Sortie :
  10. C---------------------------
  11. C
  12. C E/ P : flottant : Pression de référence (Pa)
  13. C E/ T : flottant : Température du gaz à l'entrée du PAR (K)
  14. C E/ TC : flottant : Température des plaques (K)
  15. C E/ TPROP : flottant : Température de référence (K)
  16. C E/ X(7) : flottant : Fraction molaire de N2, O2, H2 et H2O
  17. C /S FLA009 : flottant : Coefficient d'échange thermique (W/m2/K)
  18. C
  19. C------------------------------
  20. C Variables de TRAV utilisées :
  21. C------------------------------
  22. C
  23. C E/ E : Distance entre deux plaques
  24. C E/ L : Hauteur des plaques
  25. C E/ G : Gravité (=9.81 m/s2)
  26. C E/ R : Constante des gaz (= 8.314 J/mol/K)
  27. C E/ IKALP : Flag permettant d'activer la correction RALOC
  28. C E/ M(7) : Masse molaire des constituants du mélange (kg/mol)
  29. C E/ IMODEL : Type de model
  30. C E/ CSTMOD(3) : Coeff. correctif du coefficient d'échange (HEATER)
  31. C
  32. C---------------------------------------------------------------------
  33. C
  34. C Langage : ESOPE + FORTRAN 77
  35. C
  36. C Mise en oeuvre : H. Paillère (1997, TTMF)
  37. C
  38. C---------------------------------------------------------------------
  39. IMPLICIT INTEGER(I-N)
  40. IMPLICIT REAL*8 (A-H,O-Z)
  41. REAL*8 FLA009
  42. REAL*8 X(7)
  43. REAL*8 LAMBDA,NU
  44. segment trav
  45. integer iKALP,iMODEL
  46. real*8 e,L,Lch,Dh,S,sp,Ck
  47. real*8 mc,Cpc
  48. real*8 g,R,deltah
  49. real*8 M(nbesp),cstmod(ncst)
  50. real*8 Cpi(nbesp),al
  51. real*8 eps_mh2,eps_dt,eps_con,u
  52. real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  53. real*8 XHEMOY,XCO2MOY,XCOMOY
  54. endsegment
  55. C
  56. C Conductivité thermique du mélange (en W/m/K)
  57. LAMBDA = FLA006(TPROP,X,TRAV)
  58. C
  59. C Calcul du Rayleigh modifié
  60. RA = G*(TC-T)*E*E*E*FLA001(P,TPROP,X,TRAV)*P*FLA002(X,TRAV)
  61. RA = RA / (R*TPROP*TPROP*FLA004(TPROP,X,TRAV)*LAMBDA)
  62. RAB = RA * E / L
  63. C
  64. C Corrélation de Elenbaas
  65. IF (RAB .LT. 1D-3) THEN
  66. NU = 0.D0
  67. ELSE
  68. NU = RAB/24.D0 * (1.D0 -EXP(-35.D0/RAB)) ** 0.75D0
  69. ENDIF
  70. VAL009 = LAMBDA * NU / E
  71. C
  72. C Correction RALOC éventuelle
  73. IF (IKALP .EQ. 1) THEN
  74. DT = TC - T
  75. IF (DT .GE. 75.D0) THEN
  76. VAL009 = VAL009 + 40.D0
  77. ELSEIF (DT .GE. 25.D0) THEN
  78. VAL009 = VAL009 + 30.D0
  79. ENDIF
  80. ENDIF
  81. C
  82. IF (IMODEL .NE. 1) THEN
  83. val009 = val009 * cstmod(3)
  84. ENDIF
  85. C
  86. FLA009 = VAL009
  87. C
  88. RETURN
  89. END
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  

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