Télécharger la plaquette Cast3M


fr en

Télécharger fla006.eso

Retour à la liste

Numérotation des lignes : 

fla006

  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
  47.       CALL FLA005(T,MUI)
  48. C
  49. C Conductivité thermique en fonction de T de chaque constituant
  50.       CALL FLA007(T,LAMI)
  51. C
  52. C Formule de Wilke
  53.       IBOU = CPI(/1)
  54.       COEF = 1.D0 / SQRT(8.D0)
  55.       DO 20 J=1,IBOU
  56.          DO 10 I=1,IBOU
  57.             FIJ(I,J) = COEF / SQRT(1.D0+M(I)/M(J))
  58.      S      * ((1.D0 + SQRT(MUI(I)/MUI(J))*(M(J)/M(I))**0.25D0)**2.D0)
  59.  10      CONTINUE
  60.  20   CONTINUE
  61.       VAL006 = 0.D0
  62.       DO 40 I=1,IBOU
  63.          COEF = 0.D0
  64.          DO 30 J=1,IBOU
  65.             COEF = COEF + X(J)*FIJ(I,J)
  66.  30      CONTINUE
  67.          VAL006 = VAL006 + X(I)*LAMI(I)/COEF
  68.  40   CONTINUE
  69. C
  70.       FLA006 = VAL006
  71. C
  72.       RETURN
  73.       END
  74.  
  75.  
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  

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