* NOTE *

THIS VERSION OF KLCALC IS A SLIGHT MODIFICATION OF MOLLY FINLEY AND RYAN DUNNAVANT’S CODE FROM 2004.

function y = klcalc(T,index,mw,cnms,lrho)

gamma=1;

found=0;

kappa=1.3807E-23; %J/K

Nn=6.023e23; %1/gmol

if cnms(index,[1:5])=='water'

A=.00015;

B=-.1512;

C=58.912;

D=-6001.7;

a=0;

b=-.0239;

c=62.405;

rho=(a*(T)^2+b*(T)+c)*mw(index);

vs=A.*(T)^3 +B.*(T)^2+C.*(T)+D;

found=1;

Vn=mw(index)/rho/1000; % if rho in kg/m^3 in m^3/gmol

y=2.8*(Nn/Vn)^(2/3)*kappa*vs; %in units kg m /s^3/K

elseif cnms(index,[1:6])=='nitrog'

A=0;

B=0;

C=-10.145;

D=-1134.8;

a=-.0158;

b=-8.528;

c=-1122.8;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='oxygen'

A=.0037;

B=2.1299;

C=403.26;

D=25887;

a=-.0003;

b=-.2517;

c=-1.4332;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='fluori'

A=-1e-7;

B=-5e-5;

C=-.0079;

D=-.4171;

a=-.0003;

b=-.3049;

c=-6.0883;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='methan'

A=0;

B=0;

C=-9.7087;

D=-228.79;

a=-.0002;

b=-.1596;

c=6.0081;

rho=(a*(T-273.15)^2+b*(T-273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='ethane'

A=0;

B=0;

C=-7.371;

D=664.62;

a=-3e-5;

b=-.0461;

c=14.254;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='pentan'%not in 301 database

A=-2e-5;

B=.0012;

C=-4.8315;

D=1132;

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='ethene'%not in 301 database

A=0;

B=0;

C=-7.0581;

D=579.4;

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='butane'%not in 301 database

A=2e-7;

B=.0021;

C=-5.3588;

D=1028.5;

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='isobut'

A=-9e-6;

B=-.0006;

C=-5.4483;

D=933.76;

a=-1e-5;

b=-0.0188;

c=10.009;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='heptan'%not in 301 database

A=0;

B=0;

C=-4.3199;

D=1244.1;

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='hexane'%not in 301 database

A=0;

B=.0045;

C=-4.4999;

D=1164;

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='ammoni'

A=0;

B=0;

C=-7.7709;

D=1503.4;

a=-.0001;

b=-.0812;

c=37.473;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:6])=='propan'

A=0;

B=0;

C=-6.714;

D=860.28;

a=-8e-6;

b=-.0253;

c=12.156;

rho=(a*(T+273.15)^2+b*(T+273.15)+c)*mw(index);

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

elseif cnms(index,[1:7])=='propene' %not in 301 database

A=0;

B=0;

C=-6.6531;

D=836.7;

found=1;

vs=A.*(T-273.15)^3 +B.*(T-273.15)^2+C.*(T-273.15)+D;

else

disp('you suck')

vs=-99; % dummy variable

rho=-99999;

end

if vs<0

disp('This is not in the liquid range')

end

if found==0

cnms(1,:);

disp('Is not in the database for liquids')

end

rho;

Vn=mw(index)/rho/1000; % if rho in kg/m^3 in m^3/gmol

y=2.8*(Nn/Vn)^(2/3)*kappa*vs; %in units kg m /s^3/K