Esope source files

   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 :