Télécharger la plaquette Cast3M


fr en

Télécharger dzvdt0.eso

Retour à la liste

Numérotation des lignes : 

dzvdt0

  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