Télécharger la plaquette Cast3M


fr en

Télécharger fla009.eso

Retour à la liste

Numérotation des lignes : 

fla009

  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