$SB1-X5.226
� Function EBULLICION7(D,DELTAT_e,DELTAP_sat,G,x,sigma,Cpl,h_lg,rho_liq,rho_vap,mu_liq,mu_vap,k_liq,Pr)
fi=((k_liq^0.79)*(Cpl^0.45)*(rho_liq^0.49))/((sigma^0.5)*(mu_liq^0.29)*(h_lg^0.24)*(rho_liq^0.24))
X_tt=(((1-x)/x)^0.9)*((rho_liq/rho_vap)^0.5)*((mu_liq/mu_vap)^0.1)
y=1/X_tt
F=0.99597+1.7772*y-0.077198*y^2+0.0032843*y^3-0.000076809*y^4+8.8730E-07*y^5-3.9372E-09*y^6
Re_liq=(G*(1-x)*D)/mu_liq
Re_TP=(F^1.25)*Re_liq
S=0.96172-0.000010368*Re_TP+5.7136E-11*Re_TP^2-1.6500E-16*Re_TP^3+2.5107E-22*Re_TP^4-1.9052E-28*Re_TP^5+5.6741E-35*Re_TP^6
h_EN=0.00122*fi*(DELTAP_sat^0.75)*(DELTAT_e^0.24)*S
h_C=0.023*(k_liq/D)*(Re_liq^.8)*(Pr^0.4)*F
EBULLICION7=h_C+h_EN
END
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
