1 : C DHVDP0 SOURCE TTMF3 12/07/05 21:15:10 7425 2 : FUNCTION DHVDP0(P,T) 3 : C-------------------------------------------------------------------- 4 : C derivee partielle de Enthalpie vapeur spécifique par rapport à Pvap 5 : C-------------------------------------------------------------------- 6 : C 7 : C E/ P : Pression partielle de vapeur (en Pa) 8 : C E/ T : Température (en K) 9 : C /S DHVDP0 : Derivee partielle Enthalpie spécifique de la vapeur 10 : C par rapport a la pression de vapeur Pvap (en J/kg/K) 11 : C d/dP (Hvap(Pvap,T)) 12 : C Issue de la fonction enthalpie specifique HVS0(P,T) 13 : C-------------------------------------------------------------------- 14 : IMPLICIT INTEGER(I-N) 15 : IMPLICIT REAL*8 (A-H,O-Z) 16 : 17 : C COEFFICIENTS 18 : 19 : B0 = -4512.70756D0 20 : B1 = 4245585.9D0 21 : B2 = -1584624342.D0 22 : R = 461.523D0 23 : 24 : C ECRITURE DE EPSILON (Pv,T) 25 : C ************************************************* 26 : 27 : EPSILON= P*((B0*T+B1)*T+B2) 28 : EPSILON= EPSILON/(T*T*T*T)+((R*R)/4) 29 : EPSILON = sqrt(EPSILON) 30 : 31 : VAL1 = ((2*B0*T+3*B1)*T)+4*B2 32 : VAL2 = 2*EPSILON*T*T*T 33 : 34 : C ECRITURE DE LA FORMULE BRUTE 35 : C *** 36 : DHVDP0 = VAL1/VAL2 37 : RETURN 38 : END 39 : 40 : 41 : 42 :
© Cast3M 2003 - All rights reserved.
Disclaimer