USing EES to calculate the dryness of the refrigeratant R600a flowing in a capillary

[electricmouse]

moring to all!
I've recently used EES to simulate the flowing of R600a in a capillary.And I find that using the correlation: A=0.719*B^0.403*C^-0.590*D^0.157 can causing a quite big diviation up to 35% from the experimental data. However the correlation is conclude from the experimental data and the author declared that the bias was within 10%. So I want to know if anyone have similar experience? The article is <An experimental analysis of adiabatic capillary tubes>. The after is my EES code:

module getm_ch(L_cap,D_cap,P_in,DELTAT_sc,m_h)
T_f=temperature(Isobutane,p=p_in,x=0)
T_sc=T_f-DELTAT_sc
mu_f=viscosity(Isobutane,P=P_in,T=T_sc)
v_f=volume(Isobutane,P=P_in,T=T_sc)
c_pf=specheat(Isobutane,P=P_in,T=T_sc)
A=m/(D_cap*mu_f)
B=D_cap^2*p_in/(v_f*mu_f^2)
C=L_cap/D_cap
D=D_cap^2*c_pf*DELTAT_sc/(v_f^2*mu_f^2)

A=0.719*B^0.403*C^(-0.590)*D^0.157

m_h=m*convert(kg/s,kg/h)
"结果输出"
END

call getm_ch(2.926,0.00077,721000,4.2,m)
call getm_ch(2.926,0.00077,923000,3.7,m1)
call getm_ch(2.926,0.00077,1119000,15.7,m2)

[electricmouse]

To addition, the correlation and code are all to calculate the choked mass flow-rate in the capillary under certain conditions.