Télécharger la plaquette Cast3M


fr en

Télécharger dhvdt0.eso

Retour à la liste

Numérotation des lignes : 

dhvdt0

  1. C DHVDT0    SOURCE    SERRE     12/07/10    21:15:06     7431
  2.       FUNCTION DHVDT0(P,T)
  3.  
  4. C--------------------------------------------------------------------
  5. C Derivee partielle de enthalpie specifique vapeur par rapport a la
  6. C Temperature T
  7. C Issue de la fonction enthalpie specifique vapeur HVS0(P,T)
  8. C--------------------------------------------------------------------
  9. C  E/  P    : Pression partielle de vapeur (en Pa)
  10. C  E/  T    : Température (en K)
  11. C   /S DHVDT0 : Derivee de enthalpie specifique (en J/kg/K)
  12. C             d/dT (Hvap(Pvap,T))
  13. C--------------------------------------------------------------------
  14.       IMPLICIT INTEGER(I-N)
  15.       IMPLICIT REAL*8 (A-H,O-Z)
  16. C
  17. C     COEFFICIENTS
  18.  
  19.       A1=1.748383D3
  20.       A2=0.277768D0
  21.       A3=0.275D-3
  22.  
  23.  
  24.       B0 = -4512.70756D0
  25.       B1 = 4245585.9D0
  26.       B2 = -1584624342.D0
  27.       R  = 461.523D0
  28.  
  29. C     premier terme de la formule
  30.       VAL0  = A1+(A3*T+A2)*T
  31.  
  32. C     ECRITURE DE EPSILON(Pv,T)
  33. C     *************************************************
  34.       EPSILON= P*((B0*T+B1)*T+B2)
  35.       EPSILON= (EPSILON/(T*T*T*T))+((R*R)/4)
  36.       EPSILON = SQRT(EPSILON)
  37.  
  38. C     ECRITURE DE VAL1
  39. C     *************************************************
  40.       VAL1=(-R/2)+EPSILON
  41.  
  42. C     ECRITURE DE VAL2
  43. C     *************************************************
  44.       S1=2*B0*T*T*T
  45.       S2=3*B1*T*T
  46.       S3=4*B2*T
  47.       S4=S1+S2+S3
  48.       S5= (B0*T+B1)*T+B2
  49.       VAL2=S4/S5
  50.  
  51. C     ECRITURE DE BETA1
  52. C     *************************************************
  53.       X1=(6*B0*T*T)+(6*B1*T)+4*B2
  54.       X2= (B0*T*T)+(B1*T)+B2
  55.       X3=(2*B0*T*T*T)+(3*B1*T*T)+(4*B2*T)
  56.       X4=(2*B0*T)+B1
  57.       BETA1=((X1*X2)-(X3*X4))/(X2*X2)
  58.  
  59. C     ECRITURE DE BETA2
  60. C     *************************************************
  61.       Y1=(((2*B0*T+3*B1)*T)+4*B2)*(-P)
  62.       Y2=2*EPSILON*(T*T*T*T*T)
  63.       BETA2=Y1/Y2
  64.  
  65. C     ECRITURE DE LA FORMULE BRUTE
  66. C     *************************************************
  67.       DHVDT0 =VAL0 + VAL1*BETA1+VAL2*BETA2
  68.       RETURN
  69.       END
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  

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