Esope source files

   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 :