def alpha_ortho_h2(T,density): from helmholtz_energy_residual import helmholtz_energy_residual from ideal_helmholtz_energy import ideal_helmholtz_energy N_i=[ -6.83148,0.01,2.11505,4.38353,0.211292,-1.00939,0.142086,-0.87696,0.804927,-0.710775,0.0639688,0.0710858,-0.087654,0.647088] t_i=[0.7333,1,1.1372,0.5136,0.5638,1.6248,1.829,2.404,2.105,4.1,7.658,1.259,7.589,3.946] d_i=[1,4,1,1,2,2,3,1,3,2,1,3,1,1] p_i=[0,0,0,0,0,0,0,1,1,0/1,0/1,0/1,0/1,0/1]; phi_i=[0,0,0,0,0,0,0,0,0,-1.169,-0.894,-0.04,-2.072,-1.306] Beta_i=[0,0,0,0,0,0,0,0,0,-0.4555, -0.4046,-0.0869,-0.4415,-0.5743] gamma_i=[0,0,0,0,0,0,0,0,0,1.5444,0.6627,0.763,0.6587,1.4327] D_i=[0,0,0,0,0,0,0,0,0,0.6366,0.3876,0.9437,0.3976,0.9626] Tc=32.22 rho_c=15.445*2.01594 # mol/l -> kg/m^3 factors of 1000 liter converstion, and gram conversion cancel delta=density/rho_c tau=Tc/T R=8314.472 R_hydrogen=R/2.01594 ni=2.5 ti=0.0 vi=[2.54151,-2.3661,1.00365,1.22447] ui=[856,1444,2194,6968] ho = 6883.3960754547 so = 141.315444516 alpha_o=ideal_helmholtz_energy(ni,ti,vi,ui,R,ho,so,rho_c/2.01594,Tc,tau,delta) alpha_r=helmholtz_energy_residual(tau,delta,N_i,t_i,d_i,p_i,phi_i,Beta_i,gamma_i,D_i) absolute=R_hydrogen*T*(alpha_o+alpha_r) nomalized=alpha_o+alpha_r ideal=alpha_o residual=alpha_r return absolute,nomalized,ideal,residual
def alpha_normal_h2(T,density): from helmholtz_energy_residual import helmholtz_energy_residual from ideal_helmholtz_energy import ideal_helmholtz_energy from pylab import zeros N_i=[-6.93643,0.01,2.1101,4.52059,0.732564,-1.34086,0.130985,-0.777414,0.351944,-0.0211716,0.0226312,0.032187,-0.0231752,0.0557346]; t_i=[0.6844,1,0.989,0.489,0.803,1.1444,1.409,1.754,1.311,4.187,5.646,0.791,7.249,2.986] d_i=[1,4,1,1,2,2,3,1,3,2,1,3,1,1] p_i=[0,0,0,0,0,0,0,1,1,0/1,0/1,0/1,0/1,0/1] phi_i=[0,0,0,0,0,0,0,0,0,1.685,0.489,0.103,2.506,1.607] Beta_i=[0,0,0,0,0,0,0,0,0, 0.171,0.2245,0.1304,0.2785,0.3967] gamma_i=[0,0,0,0,0,0,0,0,0,0.7164,1.3444,1.4517,0.7204,1.5445] D_i=[0,0,0,0,0,0,0,0,0,1.506,0.156,1.736,0.67,1.662] betas=zeros(len(gamma_i)) Tc=33.145 rho_c=15.508*2.01594 # mol/l -> kg/m^3 factors of 1000 liter converstion, and gram conversion cancel delta=density/rho_c tau=Tc/T R=8314.472 R_hydrogen=R/2.01594 ni=2.5 ti=0.0 vi=[1.616,-0.4117,-0.792,0.758,1.217] ui=[531,751,1989,2484,6859] R=8.314472 ho = 7206.9069892047 so = 143.4846187346 alpha_o=ideal_helmholtz_energy(ni,ti,vi,ui,R,ho,so,rho_c/2.01594,Tc,tau,delta) alpha_r=helmholtz_energy_residual(tau,delta,N_i,t_i,d_i,p_i,phi_i,Beta_i,gamma_i,D_i) absolute=R_hydrogen*T*(alpha_o+alpha_r) nomalized=alpha_o+alpha_r ideal=alpha_o residual=alpha_r return absolute,nomalized,ideal,residual
def alpha_para_h2(T,density): from helmholtz_energy_residual import helmholtz_energy_residual from ideal_helmholtz_energy import ideal_helmholtz_energy N_i =[-7.33375,0.01,2.60375,4.66279,0.682390,-1.47078,0.135801,-1.05327,0.328239,-0.0577833,0.0449743,0.0703464,-0.0401766,0.119510] t_i=[0.6855,1,1,0.489,0.774,1.133,1.386,1.619,1.162,3.96,5.276,0.99,6.791,3.19] d_i=[1,4,1,1,2,2,3,1,3,2,1,3,1,1] p_i=[0,0,0,0,0,0,0,1,1,0.0/1.0,0.0/1.0,0.0/1.0,0.0/1.0,0.0/1.0] phi_i=[0,0,0,0,0,0,0,0,0,1.7437,-0.5516 ,-0.0634 ,-2.1341 ,-1.777] Beta_i=[0,0,0,0,0,0,0,0,0,0.194,-0.2019,-0.0301,-0.2383, -0.3253] gamma_i=[0,0,0,0,0,0,0,0,0,0.8048,1.5248,0.6648,0.6832,1.493] D_i=[0,0,0,0,0,0,0,0,0,1.5487,0.1785,1.28,0.6319,1.7104] Tc=32.938 rho_c=15.538*2.01594 # mol/l -> kg/m^3 factors of 1000 liter converstion, and gram conversion cancel delta=density/rho_c tau=Tc/T R=8314.472 R_hydrogen=R/2.01594 ni=2.5 ti=0.0 vi=[4.30256,13.0289,-47.7365,50.0013,-18.6261,0.993973,0.536078] ui=[499,826.5,970.8,1166.2,1341.4,5395,10185] ho = 8172.6404795208 so = 150.0625663134 alpha_o=ideal_helmholtz_energy(ni,ti,vi,ui,R,ho,so,rho_c/2.01594,Tc,tau,delta) alpha_r=helmholtz_energy_residual(tau,delta,N_i,t_i,d_i,p_i,phi_i,Beta_i,gamma_i,D_i) absolute=R_hydrogen*T*(alpha_o+alpha_r) nomalized=alpha_o+alpha_r ideal=alpha_o residual=alpha_r return absolute,nomalized,ideal,residual