Télécharger dzvdp0.eso

Retour à la liste

Numérotation des lignes :

dzvdp0
  1. C DZVDP0 SOURCE TTMF3 12/07/05 21:15:17 7425
  2. FUNCTION DZVDP0(P,T)
  3. C--------------------------------------------------------------------
  4. C derivee partielle de Z par rapport a P
  5. C Issue de la fonction compressibilite 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 DZVDP0 :Derivee de facteur compressibilite Z / Pvap (en J/kg/Pa)
  10. C d/dP (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 =(B0*T+B1)*T+B2
  24. VAL1=R*T*T*T*T
  25. VAL2=VAL0/VAL1
  26.  
  27. C ECRITURE DE EPSILON(Pv,T)
  28. C *************************************************
  29. EPSILON= P*((B0*T+B1)*T+B2)
  30. EPSILON= (EPSILON/(T*T*T*T)+((R*R)/4))
  31. EPSILON = SQRT(EPSILON)
  32.  
  33. C ECRITURE DE VAL3
  34. C *************************************************
  35. VAL3=(R/2)+EPSILON
  36.  
  37. C ECRITURE DE VAL4
  38. C *************************************************
  39. S1=(VAL0*P)/(EPSILON*T*T*T*T*2)
  40. S2=VAL3-S1
  41. S3=S2/(VAL3*VAL3)
  42.  
  43. C ECRITURE DE LA FORMULE BRUTE
  44. C *************************************************
  45. DZVDP0 =S3*VAL2
  46. RETURN
  47. END
  48.  
  49.  
  50.  
  51.  
  52.  

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