Télécharger detoh.eso

Retour à la liste

Numérotation des lignes :

  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