754.1,754.1,754.1,754.1,754.1,754.1 ];
p_l = [ 40.0,50.0,60.0,70.0,80.0,90.0,100.0,110.0,120.0,130.0,140.0,150.0,160.0,
170.0,180.0,190.0,200.0,210.0,220.0,230.0,240.0,250.0 ];
p_l = np.array(p_l)*1.0E3;
rho_l = [ 782.57,782.58,782.59,782.60,782.61,782.62,782.63,782.64,782.65,782.66,
782.67,782.68,782.69,782.70,782.71,782.72,782.73,782.74,782.75,782.76,
782.77,782.78 ];
u_rho_l = [ 0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,
0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62 ];
P_b = [ 0.0000010543,10.,20.,30.,40.,50.,60.,70.,80.,90.,100.,110.,120.,130.,
140.,150.,160.,170.,180.,190.,200.,210.,220.,230.,240.,250. ];
P_b = np.array(P_b)*1.0E3;
T_b = [ 140.0,274.0,287.7,296.5,303.2,308.7,313.3,317.4,321.0,324.3,327.3,330.0,
332.6,335.0,337.3,339.5,341.5,343.5,345.4,347.1,348.9,350.5,352.1,353.7,355.2,
356.6 ];
T_b_curve = ahm.curve(P_b,T_b);
M = 58.0791/1000.0;
p_l_250 = [ 40.0,50.0,60.0,70.0,80.0,90.0,100.0,110.0,120.0,130.0,140.0,150.0,
160.0,170.0,180.0,190.0,200.0,210.0,220.0,230.0,240.0,250.0 ];
p_l_250 = np.array(p_l_250)*1.0E3;
rho_l_250 = [ 837.0,837.0,837.0,837.0,837.0,837.0,837.0,837.0,837.0,837.0,837.0,
837.0,837.1,837.1,837.1,837.1,837.1,837.1,837.1,837.1,837.1,837.1];
V_l_250 = np.divide(M,rho_l_250);
################################### Tait #######################################
R = 8.3145; # Common definitions in m^3 Pa / mol K
T = 275.0; # K
P = np.linspace(np.min(p_l),np.max(p_l),50);
P_s = 33260.0; # Pa saturated pressure at 300 K
88,87,85,84,82,81,80,79,78,77,76,75,74,73,73,72,71,71,70,70,70,69,69,68,68,68, 68,67,67,67,67,67,67,66,66,66,66,66,66,66,66,65,65,65,65,65,65,65,64,64,64,64, 64,63,63,63,63,62,62,62,61,61,61,60,60,60,59,59,58,58,57,57,56,56,55,55,54,54, 53,52,52,51,51,50,49,49,48,47,46,46,45,44,43,43,42,41,40,39,39,38,37,36,35,34, 33,33,32,31,30,29,28,27,26,25,24,24,23,22,21,20,19,18,17,16,16,15,14,13,12,12, 11,10,9.6,9.1,8.7,8.4,8.2,8.1,8.1,8.3,8.5,8.9,9.3,9.8,10,11,12,12,13,14,15,15, 16,17,18,19,19,20,21,22,23,23,24,25,26,27,27,28,29,30,30,31,32,32,33,34,34,35, 36,36,37,38,38,39,40,40,41,41,42,42,43,43,44,44,45,45,46,46,47,47,48,48,49,49, 49,50,50,51,51,51,52,52,53,53,53,54,54,55,55,56,56,56,57,57,58,58,59,59,60,60, 61,61,62,63,63,64,65,65,66,67,68,68,69,70,71,72,73,74,75,76,77,78,79,80,82,83, 84,85,87,88,89,91,92,94,95,97,98,100,100,100,100,110,110,110,110,110,110,120, 120,120,120,120,120,120,130,130,130,130,130,130,130,140,140,140,140,150,150, 160,250 ]; u_rho = [float(i) for i in u_rho]; rho = ahm.curve(T_val,rho_val); M = 252.055032/1000.0; # [ kg / mol ] T = 367.057; # [ K ] V = M/rho.at(T); # [ kg / mol ] / [ kg / m^3 ] [=] [ m^3 / mol ] ################################### Soave ###################################### omega = 0.5205; T_r = np.divide(T,T_c); print T_r; # Soave uses different definitions for kappa and alpha kappa = 0.480 + 1.574 * omega - 0.176 * np.power(omega,2.0); alpha = np.power((1.0 + kappa * (1.0-np.sqrt(T_r))),2.0); a = (0.42747 * (np.power(R,2.0) * np.power(T_c,2.0) / P_c)) * alpha;