function y = kscalc(T,cmpnm)

% calculates the thermal conductivity of solids
% kscalc (temperature, number of compound)
% 1 silver
% 2 copper
% 3 iron
% 4 lead
% 5 gold
% 6 aluminum
% 7 zinc
% 8 tin
% 9 ceramic
% 10 glass
% k is in units watts/m/K

if cmpnm==1
A=-7E-10;
B=4E-7;
C=-.0008;
D=4.506;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==2
A=-9E-10;
B=2E-6;
C=-.0019;
D=4.4396;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==3
A=0;
B=4E-7;
C=-.0008;
D=4.506;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==2
A=-9E-10;
B=2E-6;
C=-.0019;
D=4.4396;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==3
A=0;
B=4E-7;
C=-.0012;
D=1.0236;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==4
A=-3E-9;
B=4E-6;
C=-.0014;
D=.5017;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==5
A=-8E-10;
B=2E-6;
C=-.0016;
D=3.5548;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==6
A=-8E-10;
B=6E-7;
C=-.0003;
D=2.4437;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==7
A=4E-10;
B=-4E-7;
C=-.0004;
D=1.359;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==8
A=-4E-9;
B=5E-6;
C=-.0023;
D=1.0351;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==9
A=-2E-11;
B=7E-8;
C=-7E-5;
D=.0563;
y=100*(A*T^3 + B*T^2 + C*T +D);
end

if cmpnm==10
A=6E-11;
B=-9E-8;
C=6E-5;
D=.0006;
y=100*(A*T^3 + B*T^2 + C*T +D);
end
%Mollie Finley and Ryan Dunnavant 4/2004

HomeGasesLiquidsSolidsSources
© 2003 Finley and Dunnavant