Télécharger la plaquette Cast3M


fr en

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