clf;

% 1. General Electric ZrC, pressed and sintered, density 5.94 g/cm^3
kappa1 = [
  528.2   67.8
  530.4   70.7
  531.5   71.4
  535.4   66.3
  818.2   30.3
  829.3   30.4
  833.2   33.5
  833.7   34.6
  845.4   37.6
  846.5   41.2
  ];
plot (kappa1(:,1), kappa1(:,2), 'o');
hold on;
% 2. General Electric ZrC, pressed and sintered, density 5.94 g/cm^3
kappa2 = [
 580.9   30.3
 595.9   26.5
 597.1   25.8
 785.4   44.6
 787.1   44.3
 793.2   42.7
 793.2   45.0
 1097.1  15.3
 1101.5  15.1
 1105.4  16.2
 1107.1  14.3
 1108.2  15.0
 1110.9  16.3
  ];
plot (kappa2(:,1), kappa2(:,2), 'o');
% 3. General Electric ZrC, pressed and sintered, density 5.94 g/cm^3
kappa3 = [
1308.2   38.4
1308.2   34.9
1313.7   41.8
1702.6   42.5
1702.6   43.0
1705.4   45.9
1705.4   45.0
1972.1   38.1
1983.2   51.3
2191.5   39.2
2199.8   40.4
2211.0   51.9
2436.0   44.7
2436.0   43.6
2444.3   48.5
2647.1   52.5
2649.8   50.2
  ];
plot (kappa3(:,1), kappa3(:,2), 'o');
% 4. General Electric ZrC, pressed and sintered, density 5.94 g/cm^3
kappa4 = [
1449.8   22.1
1449.8   19.6
1933.2   46.0
1938.7   39.8
2408.2   47.9
2424.8   47.2
2424.8   43.0
2424.8   43.4
  ];
plot (kappa4(:,1), kappa4(:,2), 'o');


% 5. Wah Chang Co. ZrC (11.7C), lattice parameter 4.699A, density = 6.06 g/cm^3
kappa5 = [
1965    29.3
2120    29.9
2225    29.9
2265    36.0
2470    33.1
  ];
plot (kappa5(:,1), kappa5(:,2), 'r*');
% 6. Wah Chang Co. ZrC (11.8C), lattice parameter 4.699A, density = 6.06 g/cm^3
kappa6 = [
  1490   26.2
  1560   26.6
  1655   29.3
  1795   30.3
  1905   30.5
  2055   30.8
  2170   33.8
  2310   34.6
  2400   34.1
  ];
plot (kappa6(:,1), kappa6(:,2), 'r*');


% 7. Carborundum Co. ZrC (11.0C), density = 6.13 g/cm^3
kappa7 = [
  813.2   30.8
  873.2   32.4
  898.2   34.7
  938.2   35.4
  958.2   33.7
  983.2   35.4
  998.2   35.1
  1008.2  34.1
  1008.2  32.6
  1073.2  32.0
  1083.2  35.6
  1098.2  35.1
  1118.2  32.4
  1138.2  32.0
  1183.2  33.5
  1208.2  35.6
  1293.2  35.1
  1523.2  36.8
  1583.1  38.5
  1613.2  41.8
  1728.2  39.5
  1793.2  40.2
  1843.2  37.7
  2053.2  40.4
  2208.2  43.1
  ];
plot (kappa7(:,1), kappa7(:,2), 'g^');
% 8. Carborundum Co. ZrC (12.1C), density = 6.17 g/cm^3
kappa8 = [
  813.2   30.1
  963.2   32.6
  993.2   32.6
  1038.2  32.0
  1048.2  35.1
  1053.2  34.7
  1163.2  32.6
  1163.2  33.3
  1238.2  39.3
  2423.2  42.3
  ];
plot (kappa8(:,1), kappa8(:,2), 'g^');
% 9. Carborundum Co. ZrC (11.0C), density = 6.18 g/cm^3
kappa9 = [
  1368.2   35.8
  1493.2   41.6
  1608.2   38.9
  1633.2   42.7
  1673.2   35.6
  1703.2   38.5
  1708.2   40.8
  1773.2   39.7
  1783.2   37.0
  1838.2   38.7
  1883.2   42.0
  1923.2   40.8
  1973.2   38.9
  2088.2   42.0
  2338.2   43.9
  ];
plot (kappa9(:,1), kappa9(:,2), 'g^');


% 10. pressed and sintered ZrC (10.606C), electrical resistivity = 100.0E-6 ohm.cm
kappa10 = [298.2 11.6];
plot (kappa10(:,1), kappa10(:,2), 'k>'); 
% 11. pressed and sintered ZrC (9.725C), electrical resistivity = 123.6E-6 ohm.cm
kappa11 = [298.2 10.3];
plot (kappa11(:,1), kappa11(:,2), 'k>'); 
% 12. pressed and sintered ZrC (9.035C), electrical resistivity = 130.2E-6 ohm.cm
kappa12 = [298.2 11.0];
plot (kappa12(:,1), kappa12(:,2), 'k>'); 
% 13. pressed and sintered ZrC (8.639C), electrical resistivity = 130.0E-6 ohm.cm
kappa13 = [298.2 9.4];
plot (kappa13(:,1), kappa13(:,2), 'k>'); 
% 14. pressed and sintered ZrC (7.644C), electrical resistivity = 166.2E-6 ohm.cm
kappa14 = [298.2 8.2];
plot (kappa14(:,1), kappa14(:,2), 'k>'); 


% 15. specific heat = 0.125 cal/g/K
kappa15 = [
1573  38.5
1673  40.6
1773  42.3
1873  43.9
  ];
plot (kappa15(:,1), kappa15(:,2), 'm+');


% Wah Chang Co.; hot pressed; density 6.26 g/cm^3 (theoretical %
% density 6.59 g/cm^3); lattice constant 4.6987 A; 87.1 Zr, %
% 11.2 C (0.6 free carbon), carbon/metal ratio 0.924; and 0.6 O2 - %
% The Thermal Diffusivity and Conductivity of Zirconium Carbide and %
% Niobium Carbide from 100 to 2500K, by B.H. Morrison and L.L. Sturgess, %
% Rev. Int. Hautes Temper. et Refract., 7, pp 351-358, 1970 %

kappa_Morrison = [
    100 23.4
    150 20.9
    200 19.4
    250 19.4
    300 19.6
    400 20.2
    500 20.5
    600 20.7
    700 20.9
    800 21.8
    900 22.2
    1000 23.3
    1200 25.1
    1400 27.1
    1600 28.9
    1800 30.7
    2000 32.3
    2200 34.5
    2400 36.1
    2600 38.3 ];
plot (kappa_Morrison(:,1), kappa_Morrison(:,2), 'ks');

% CRC Materials Science and Engineering Handbook 2nd ed., %
% edited by J.F. Shackelford, W. Alexander and J.S. Park (CRC Press, %
% 1994), pp.279-280, for Zirconium Monocarbide (ZrC) %

kappa_CRC = [
    20 0.049
%     50 0.098 %
%     150 0.069 %
    188 0.065
%     288 0.061 %
    600 0.08
    800 0.083
    1000 0.086
    1200 0.089
    1400 0.092
    1600 0.096
    1800 0.099
    2000 0.103
    2200 0.105 ];
kappa_CRC(:,1) = kappa_CRC(:,1) + 273.16;
kappa_CRC(:,2) = kappa_CRC(:,2) * 4.186800 * 100;
plot (kappa_CRC(:,1), kappa_CRC(:,2), 'rd');
plot (kappa_CRC(:,1), kappa_CRC(:,2), 'r-');

xlabel ('T [K]');
ylabel ('\kappa [W/m/K]');
title ('Experimental ZrC Thermal Conductivity');
set(text(240,24,'CRC values'), 'Color', 'red');

axis([0 2700 8 53]);
set(gca,'XTick',0:300:2700)
grid on;

print -depsc oldie.eps
