Télécharger dhvdt0.eso

Retour à la liste

Numérotation des lignes :

  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