Télécharger la plaquette Cast3M


fr en

Télécharger detoh.eso

Retour à la liste

Numérotation des lignes : 

detoh

  1. C DETOH     SOURCE    CHAT      05/01/12    22:44:24     5004
  2.       SUBROUTINE DETOH(T,X,AN,AH2,AH2O,AN2,AO2,TINT,H,IINT,NBC,NBT)
  3. C---------------------------------------------------------------------
  4. C Calcul de l'enthalpie spécifique du mélange à la température T
  5. C---------------------------------------------------------------------
  6. C
  7. C---------------------------
  8. C Parametres Entree/Sortie :
  9. C---------------------------
  10. C
  11. C E/  T         : Température du mélange (K)
  12. C E/  X         : Fraction molaire (-)
  13. C E/  AN(NBC)   : Masse molaire des constituants (gramme)
  14. C E/  AH2(NBT)  : Coeff du devt en T des propriétés de h2
  15. C E/  AH2O(NBT) : Coeff du devt en T des propriétés de h2o
  16. C E/  AN2(NBT)  : Coeff du devt en T des propriétés de n2
  17. C E/  AO2(NBT)  : Coeff du devt en T des propriétés de o2
  18. C E/  TINT      : Température seuil (K)
  19. C  /S H         : Enthalpie spécifique du mélange (J/kg)
  20. C E/  IINT      : Nombre de coeff nécessaire à T donnée
  21. C E/  NBC       : Nombre de constituants du mélange
  22. C E/  NBT       : Nombre de coeff total
  23. C
  24. C---------------------------------------------------------------------
  25. C
  26. C Langage : FORTRAN
  27. C
  28. C Auteurs : F.DABBENE et E.STUDER 02/99
  29. C
  30. C---------------------------------------------------------------------
  31.       IMPLICIT INTEGER(I-N)
  32.       IMPLICIT REAL*8 (A-H,O-Z)
  33.       DIMENSION X(*),AN(*),AH2(*),AH2O(*),AN2(*),AO2(*)
  34. C
  35. C- Décalage dans la table suivant la borne TINT
  36. C
  37.       IINTM2 = IINT - 2
  38.       IINTM1 = IINT - 1
  39.       II = 0
  40.       IF (T.LT.TINT) II = IINT
  41. C
  42. C- Calcul de l'enthalpie spécifique de chaque constituant en J/mole
  43. C
  44.       HH2  = 0.D0
  45.       HO2  = 0.D0
  46.       HN2  = 0.D0
  47.       HH2O = 0.D0
  48.       DO 10 I = 1,IINTM2
  49.         HH2  = HH2  + (AH2 (I+II) * (T**(I-1)) /I)
  50.         HO2  = HO2  + (AO2 (I+II) * (T**(I-1)) /I)
  51.         HN2  = HN2  + (AN2 (I+II) * (T**(I-1)) /I)
  52.         HH2O = HH2O + (AH2O(I+II) * (T**(I-1)) /I)
  53.  10   CONTINUE
  54.       HH2  = HH2  + (AH2 (IINTM1+II) / T)
  55.       HO2  = HO2  + (AO2 (IINTM1+II) / T)
  56.       HN2  = HN2  + (AN2 (IINTM1+II) / T)
  57.       HH2O = HH2O + (AH2O(IINTM1+II) / T)
  58. C
  59.       HH2  = HH2  * 8.314D0 * T
  60.       HO2  = HO2  * 8.314D0 * T
  61.       HN2  = HN2  * 8.314D0 * T
  62.       HH2O = HH2O * 8.314D0 * T
  63. C
  64. C- Calcul de l'enthalpie spécifique du mélange en J/kg
  65. C- (c'est pourquoi on divise par la masse molaire AM)
  66. C
  67.       H  = X(1)*HH2   + X(2)*HO2   + X(3)*HH2O  + X(4)*HN2
  68.       AM =(X(1)*AN(1) + X(2)*AN(2) + X(3)*AN(3) + X(4)*AN(4))*1.D-3
  69.       H  = H / AM
  70. C
  71.       RETURN
  72.       END
  73.  
  74.  
  75.  

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