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
