Esope source files

   1 : C FLA006    SOURCE    SERRE     12/07/10    21:15:12     7431           
   2 :       FUNCTION FLA006(T,X,TRAV)
   3 : C---------------------------------------------------------------------
   4 : C Calcul de la conductivité thermique du mélange (en W/m/K) en fonction  
   5 : C de la composition du mélange pour une température donnée par la
   6 : C formule de Wilke
   7 : C---------------------------------------------------------------------
   8 : C
   9 : C---------------------------
  10 : C Parametres Entree/Sortie :
  11 : C---------------------------
  12 : C
  13 : C E/  T      : flottant : Température de référence
  14 : C E/  X(7)   : flottant : Fraction molaire de N2, O2, H2 et H2O
  15 : C  /S FLA006 : flottant : Conductivité thermique du mélange (W/m/K)
  16 : C
  17 : C------------------------------
  18 : C Variables de TRAV utilisées :
  19 : C------------------------------
  20 : C
  21 : C E/  M(7)   : Masse molaire des constituants du mélange (J/kg/K)
  22 : C
  23 : C---------------------------------------------------------------------
  24 : C
  25 : C Langage : ESOPE + FORTRAN 77
  26 : C
  27 : C Mise en oeuvre : H. Paillère (1997, TTMF)
  28 : C
  29 : C---------------------------------------------------------------------
  30 :       IMPLICIT INTEGER(I-N)
  31 :       IMPLICIT REAL*8 (A-H,O-Z)
  32 :       REAL*8 FLA006
  33 :       REAL*8 X(7),MUI(7),LAMI(7),FIJ(7,7)
  34 :       segment trav
  35 :            integer iKALP,iMODEL
  36 :            real*8 e,L,Lch,Dh,S,sp,Ck
  37 :            real*8 mc,Cpc
  38 :            real*8 g,R,deltah
  39 :            real*8 M(nbesp),cstmod(ncst)
  40 :            real*8 Cpi(nbesp),al
  41 :            real*8 eps_mh2,eps_dt,eps_con,u
  42 :            real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  43 :            real*8 XHEMOY,XCO2MOY,XCOMOY
  44 :       endsegment
  45 : C
  46 : C Viscosité dynamique en fonction de T de chaque constituant
  46 :       CALL FLA005(T,MUI)
  47 : C
  48 : C Conductivité thermique en fonction de T de chaque constituant
  48 :       CALL FLA007(T,LAMI)
  49 : C
  50 : C Formule de Wilke
  51 :       IBOU = CPI(/1)
  52 :       COEF = 1.D0 / SQRT(8.D0)
  53 :       DO 20 J=1,IBOU
  54 :          DO 10 I=1,IBOU
  55 :             FIJ(I,J) = COEF / SQRT(1.D0+M(I)/M(J))
  56 :      S      * ((1.D0 + SQRT(MUI(I)/MUI(J))*(M(J)/M(I))**0.25D0)**2.D0)
  57 :  10      CONTINUE
  58 :  20   CONTINUE
  59 :       VAL006 = 0.D0
  60 :       DO 40 I=1,IBOU
  61 :          COEF = 0.D0
  62 :          DO 30 J=1,IBOU
  63 :             COEF = COEF + X(J)*FIJ(I,J)
  64 :  30      CONTINUE
  65 :          VAL006 = VAL006 + X(I)*LAMI(I)/COEF
  66 :  40   CONTINUE
  67 : C
  68 :       FLA006 = VAL006
  69 : C
  70 :       RETURN
  71 :       END
  72 : 
  73 : 
  74 :  
  75 :  
  76 :  
  77 :  
  78 :  
  79 :  
  80 :