Télécharger la plaquette Cast3M


fr en

Télécharger fla001.eso

Retour à la liste

Numérotation des lignes : 

fla001

  1. C FLA001    SOURCE    SERRE     12/07/10    21:15:09     7431
  2.       FUNCTION FLA001(P,T,X,TRAV)
  3. C------------------------------------------------------------------------
  4. C Calcul de la densité du mélange en fonction de la pression, de la
  5. C température et de la composition du mélange (on suppose que la loi des
  6. C gaz parfaits s'applique aussi à la vapeur d'eau)
  7. C------------------------------------------------------------------------
  8. C
  9. C---------------------------
  10. C Parametres Entree/Sortie :
  11. C---------------------------
  12. C
  13. C E/  P      : flottant : Pression de référence (Pa)
  14. C E/  T      : flottant : Température de référence (K)
  15. C E/  X(7)   : flottant : Fraction molaire de N2, O2, H2 et H2O
  16. C  /S FLA001 : flottant : Densité du mélange (kg/m3)
  17. C
  18. C------------------------------
  19. C Variables de TRAV utilisées :
  20. C------------------------------
  21. C
  22. C E/  R      : Constante des gaz (= 8.314 J/mol/K)
  23. C E/  M(7)   : Masse molaire des constituants du mélange (kg/mol)
  24. C
  25. C------------------------------------------------------------------------
  26. C
  27. C Langage : ESOPE + FORTRAN 77
  28. C
  29. C Mise en oeuvre : H. Paillère (1997, TTMF)
  30. C
  31. C------------------------------------------------------------------------
  32.       IMPLICIT INTEGER(I-N)
  33.       IMPLICIT REAL*8 (A-H,O-Z)
  34.       REAL*8 FLA001
  35.       REAL*8 X(7)
  36.       segment trav
  37.            integer iKALP,iMODEL
  38.            real*8 e,L,Lch,Dh,S,sp,Ck
  39.            real*8 mc,Cpc
  40.            real*8 g,R,deltah
  41.            real*8 M(nbesp),cstmod(ncst)
  42.            real*8 Cpi(nbesp),al
  43.            real*8 eps_mh2,eps_dt,eps_con,u
  44.            real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  45.            real*8 XHEMOY,XCO2MOY,XCOMOY
  46.       endsegment
  47. C
  48.       ibou = cpi(/1)
  49.       sum  = 0.D0
  50.       do 10 i=1,ibou
  51.          sum = sum + x(i)*m(i)
  52.  10   continue
  53.       FLA001=P*sum/(R*T)
  54. C
  55.       RETURN
  56.       END
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  

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