Eric L. Pheterson - Homework 5 - Problem 10.61

function p10_61()
fprintf('\n\nEric L. Pheterson - Homework 5 - Problem 10.61\n\n');

Ts=150+273; %K
D=100E-3; %m

%Saturated Ethylene Glycol Vapor
Tsat=470; %K
hfg=812E3; %J/kg
pv=0; %kg/m^3

%Liquid Ethylene Glycol @ T-373K
pl=1058.5; %kg/m^3
cpl=2742; %J/kg-K
ul=0.215E-2; %N-s/m^2
kl=0.263; %W/m-K

Ja = cpl*(Tsat-Ts)/hfg
hfgp=hfg*(1+0.68*Ja)

A = pi*D^2

C=0.826;
g=9.81; %m/s^2
hd=C*(g*pl*(pl-pv)*kl^3*hfgp/(ul*(Tsat-Ts)*D))^(1/4)

m=hd*A*(Tsat-Ts)/hfgp