Télécharger dzvdt0.eso

Retour à la liste

Numérotation des lignes :

  1. C DZVDT0 SOURCE TTMF3 12/07/05 21:15:18 7425
  2. FUNCTION DZVDT0(P,T)
  3. C--------------------------------------------------------------------
  4. C Derivee partielle de Z par rapport a T
  5. C Issue de la fonction ZVAP0(P,T)
  6. C--------------------------------------------------------------------
  7. C E/ P : Pression partielle de vapeur (en Pa)
  8. C E/ T : Température (en K)
  9. C /S DZVDT0 : Derivee de facteur compressibilite Z / Pvap (en J/kg/K)
  10. C d/dT (Zvap(Pvap,T))
  11. C--------------------------------------------------------------------
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8 (A-H,O-Z)
  14. C
  15. C COEFFICIENTS
  16.  
  17. B0 = -4512.70756D0
  18. B1 = 4245585.9D0
  19. B2 = -1584624342.D0
  20. R = 461.523D0
  21.  
  22. C premier terme de la formule
  23. VAL0 = ((2*B0*T+3*B1)*T+4*B2)*P
  24. VAL1=R*T*T*T*T*T
  25. VAL2=VAL0/VAL1
  26.  
  27.  
  28. C ECRITURE DE EPSILON(Pv,T)
  29. C *************************************************
  30. EPSILON= P*((B0*T+B1)*T+B2)
  31. EPSILON= (EPSILON/(T*T*T*T))+((R*R)/4)
  32. EPSILON = sqrt(EPSILON)
  33.  
  34. C ECRITURE DE VAL3
  35. C *************************************************
  36. VAL3=(R/2)+EPSILON
  37.  
  38. C ECRITURE DE VAL4
  39. C *************************************************
  40. S1 = ((B0*T+B1)*T+B2)*P
  41. S2 = 2*T*EPSILON*T*T*T
  42. S3 = (S1/S2)/(VAL3*VAL3)
  43. S4 = 1/(VAL3)
  44. VAL4=S3-S4
  45.  
  46. C ECRITURE DE LA FORMULE BRUTE
  47. C *************************************************
  48. DZVDT0 =VAL4*VAL2
  49. RETURN
  50. END
  51.  
  52.  
  53.  
  54.  
  55.  

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