THERMAL DIFFUSIVlTY MAPPING OF 4D CARBON-CARBON

Similar documents
MASTER ADVANCED THERMAL IMAGING OF COMPOSITES INTRODUCTION ABSTRACT. MS 6064, Building H. Wang and R. B. Dinwiddie Oak Ridge National Laboratory

C. Deemer, J. G. Sun, and W. A. Ellingson. Energy Technology Division Argonne National Laboratory Argonne, IL April.2000

PROJECT PROGRESS REPORT (03/lfi?lfibr-~/15/1998):

DE '! N0V ?

Three-Dimensional Silicon Photonic Crystals

Excitations of the transversely polarized spin density. waves in chromium. E3-r 1s

THE APPLICATION OF THERMAL DIFFUSIVITY IMAGING TO SIC-FIBER

CQNl_" RESPONSE TO 100% INTERNAL QUANTUM EFFICIENCY SILICON PHOTODIODES TO LOW ENERGY ELECTRONS AND IONS

Plasma Response Control Using Advanced Feedback Techniques

Magnetic Measurements of the Elliptical Multipole Wiggler Prototype

GA A22722 CENTRAL THOMSON SCATTERING UPGRADE ON DIII D

Data Comparisons Y-12 West Tower Data

Determine the Inside Wall Temperature of DSTs using an Infrared Temperature Sensor

Magnetic Resonance Imaging of Polymer Electrolytes and Insertion Electrodes*

ust/ aphysics Division, Argonne National Laboratory, Argonne, Illinois 60439, USA

RESONANCE INTERFERENCE AND ABSOLUTE CROSS SECTIONS IN NEAR-THRESHOLD ELECTRON-IMPACT EXCITATION OF MULTICHARGED IONS

Paper to be submitted to. Thermal Management Technology for Hydrogen Storage: Fullerene Option?

Optimization of NSLS-II Blade X-ray Beam Position Monitors: from Photoemission type to Diamond Detector. P. Ilinski

: Y/LB-16,056 OAK RIDGE Y-12 PLANT

Modeling Laser and e-beam Generated Plasma-Plume Experiments Using LASNEX

RESEARCH OPPORTUNITIES AT THE CNS (Y-12 AND PANTEX) NUCLEAR DETECTION AND SENSOR TESTING CENTERS (NDSTC)

Summary Report: Working Group 5 on Electron Beam-Driven Plasma and Structure Based Acceleration Concepts

Tell uric prof i 1 es across the Darrough Known Geothermal Resource Area, Nevada. Harold Kaufniann. Open-file Report No.

MAGNETIC RESONANCE IMAGING OF SOLVENT TRANSPORT IN POLYMER NETWORKS

(4) How do you develop an optimal signal detection technique from the knowledge of

12/16/95-3/15/96 PERIOD MULTI-PARAMETER ON-LINE COAL BULK ANALYSIS. 2, 1. Thermal Neutron Flux in Coal: New Coal Container Geometry

OAK RIDGE NATIONAL LABORATORY

UNIT OPERATIONS AND CHEMICAL PROCESSES. Alan M. Krichinsky

m w n? r i OF,THISDOCUMENT IS UNLIMITED

PRECISION STUDIES OF NUCLEI

August 3,1999. Stiffness and Strength Properties for Basic Sandwich Material Core Types UCRL-JC B. Kim, R.M. Christensen.

Determine the Inside Wall Temperature of DSTs using an Infrared Temperature Sensor

Inside Wall Temperature Measurements of DSTs Using an Infrared Temperature Sensor

Turbulent Scaling in Fluids. Robert Ecke, MST-10 Ning Li, MST-10 Shi-Yi Li, T-13 Yuanming Liu, MST-10

PROJECT PROGRESS REPORT (06/16/1998-9/15/1998):

C o f l F - 9m307-- SPREADSHEET APPLICATION TO CLASSIFY RADIOACTIVE MATERIAL FOR SHIPMENTa. A. N. Brown

sample-specific X-ray speckle contrast variation at absorption edges $ & ~ 0

Parametric Instabilities in Laser/Matter Interaction: From Noise Levels to Relativistic Regimes

UPGRADED CALIBRATIONS OF THE THOMSON SYSTEM AT DIII D

IMAGING OF HETEROGENEOUS MATERIALS BY. P. Staples, T. H. Prettyman, and J. Lestone

GA A22677 THERMAL ANALYSIS AND TESTING FOR DIII D OHMIC HEATING COIL

AC dipole based optics measurement and correction at RHIC

PROCEEDINGS THIRD WORKSHOP GEOTHERMAL RESERVOIR ENGINEERING. December 14-15,1977

Safety Considerations for Laser Power on Metals in Contact with High Explosives-Experimental and Calculational Results

in the pinch. This paper describes the computer modeling behind the shielding design of a

ASTER. Dose Rate Visualization of Radioisotope Thermoelectric. Generators. RECElVED BEC OF THIS DOCUMEMT IS UNLlMIi L, Hanford Company

INITIAL RESULTS OF THE COMMISSIONING OF THE HRIBF RECOIL MASS SPECTROMETER

GA A24166 SUPER-INTENSE QUASI-NEUTRAL PROTON BEAMS INTERACTING WITH PLASMA: A NUMERICAL INVESTIGATION

~ _- IDOCLRXKT DMSP SATELLITE DETECTIONS OF GAMMA-RAY BURSTS. J. Terrell, NIS-2 P. Lee, NIS-2 R W.

Use of Gamma Rays from the Decay of 13.8-sec "Be to Calibrate a Germanium Gamma Ray Detector for Measurements up to 8 MeV

HEAT TRANSFER AND THERMAL STRESS ANALYSES OF A GLASS BEAM DUMP

RECEIVED SEQ

TARGET PLATE CONDITIONS DURING STOCHASTIC BOUNDARY OPERATION ON DIII D

APPLICATION SINGLE ION ACTIVITY COEFFICIENTS TO DETERMINE SOLVENT EXTRACTION MECHANISM FOR COMPONENTS OF NUCLEAR WASTE

Solid Phase Microextraction Analysis of B83 SLTs and Core B Compatibility Test Units

Survey, Alignment and Beam Stability at the Advanced Light Source

John Falconer Richard Noble

Scaling between K+ and proton production in nucleus-nucleus collisions *

Monte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions

SPHERICAL COMPRESSION OF A MAGNETIC FIELD. C. M. Fowler DISCLAIMER

Modeling of MHD Equilibria and Current Profile Evolution during the ERS Mode in TFTR

Spontaneous Lateral Composition Modulation in InAlAs and InGaAs Short-Period Superlattices

National Accelerator Laboratory

LQSAlamos National Laboratory

BWXT Y-12 Y-12. A BWXT/Bechtel Enterprise SMALL, PORTABLE, LIGHTWEIGHT DT NEUTRON GENERATOR FOR USE WITH NMIS

Plutonium 239 Equivalency Calculations

Manifestations of the axial anomaly in finite temperature QCD

RICE COAL COMBUSTION: EFFECT OF PROCESS CONDITIONS ON CHAR REACTIVITY. Quarterly Technical Report Performance Period: 10/1/94 42/31/94 (Quarter #13)

Simulation of Double-Null Divertor Plasmas with the UEDGE Code

A Few-Group Delayed Neutron Model Based on a Consistent Set of Decay Constants. Joann M. Campbell Gregory D. Spriggs

J. T. Mihalczo. P. 0. Box 2008

Development of a High Intensity EBIT for Basic and Applied Science

Bulk Modulus Capacitor Load Cells

Fission-Fusion Neutron Source

Direct Determination of the Stacking Order in Gd 2 O 3 Epi-Layers on GaAs

Use of AVHRR-Derived Spectral Reflectances to Estimate Surface Albedo across the Southern Great Plains Cloud and Radiation Testbed (CART) Site

GA A22689 SLOW LINER FUSION

0STI. E. Hammerberg, XNH MD SIMULATIONS OF DUSTY PLASMA CRYSTAL FORMATION: PRELIMINARY RESULTS M. J. S. Murillo, XPA

Gary E. Giles Computational Physics and Engineering Division. Jon S. Bullock Y-17npwlopment Division

GA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIII-D

Measurement of Low Levels of Alpha in 99M0Product Solutions

8STATISTICAL ANALYSIS OF HIGH EXPLOSIVE DETONATION DATA. Beckman, Fernandez, Ramsay, and Wendelberger DRAFT 5/10/98 1.

Start-up Noise in 3-D Self-AmpMed

Laser-Modulator System

Cable Tracking System Proposal. Nick Friedman Experimenta1 Facilities Division. June 25,1993. Advanced Photon Source Argonne National Laboratory

Alex Dombos Michigan State University Nuclear and Particle Physics

Analysis of Shane Telescope Aberration and After Collimation

GA A25658 PRECISION X-RAY OPTICAL DEPTH MEASUREMENTS IN ICF SHELLS

A Distributed Radiator, Heavy Ion Driven Inertial Confinement Fusion Target with Realistic, Multibeam Illumination Geometry

Diffractive Dijet Search with Roman Pots at CDF

INTERMOLECULAR POTENTIAL FUNCTIONS AND HIGH RESOLUTION MOLECULAR SPECTROSCOPY OF WEAKLY BOUND COMPLEXES. Final Progress Report

Abstract. results that address this question for the main linacs of the NLC. We will show that the effects of alignment drifts can indeed be handled.

SMALL BIPOLARONS IN BORON CARBIDES: PAIR BREAKING IN SEMICLASSICAL HOPPING* David Emin Sandia National Laboratories Albuquerque, NM , USA

EFFECTS OF AXIAL HEAT CONDUCTION AND MATERIAL PROPERTIES ON THE PERFORMANCE CHARACTERISTICS OF A THERMAL TRANSIENT ANEMOMETER PROBE

Initial Results from the Lick Observatory Laser Guide Star Adaptive Optics System

Application of a Three-Dimensional Prognostic Model During the ETEX Real-Time Modeling Exercise: Evaluatin of Results (u)

U. S. Department of Energy

Analysis of D Pellet Injection Experiments in the W7-AS Stellarator*

Flowing Interval Spacing Parameter for Matrix Diffusion in the Saturated Zone

STUDY OF CYCLOTRON RESONANCE AND MAGNETO-PHOTOLUMINESCENCE OF N-TYPE MODULATION DOPED InGaAs QUANTUM WELL LAYERS AND THEIR CHARACTERIZATIONS

Transcription:

THERMAL DIFFUSIVlTY MAPPING OF 4D CARBON-CARBON COMPOSITES H. Wang and R. B. Dinwiddie High Temperature Materials Laboratory, Metals and Ceramics Division Oak Ridge National Laboratory Oak Ridge, TN 37831-6064 ABSTRACT High resolution, 2-D thermal diffusivity maps of carbon-carbon composites were obtained by a state-of-the-art infrared thermal imaging system. Unlike the traditional single-point IR detector used for thermal diffusivity measurements, the IR camera is capable of capturing images in its 256 x 256 pixel Focal Plane Array detector in a snap-shot mode. The camera takes up to 200 images at a rate of 120 frameshecond. The temperature resolution of the IR camera is 0.015"Cand the spatial resolution is 20 pm. Thermal diffusivity was calculated for each pixel. Four-direction carbon-carbon composites were used for the thermal diffusivity mapping study. The fiber bundles along the heat flow direction were found to have 25% higher diffusivity values than the surrounding matrix. The diffusivity map also showed detailed local variations in diffusivity which were impossible to measure using a single-point detector. Accurate diffusivity maps are very important to the design of composite materials. INTRODUCTION Composite materials, such as carbon-carbon composites, have been used as structural components for high-temperature applications. One of the important characteristics of a composite is its thermal transport property, i.e. thermal conductivity for steady state heat conduction or thermal diffusivity for transient heat conduction. Accurate measurements of thermal transport properties are crucial to the component design. Thermal diffusivity[11 is the most often used technique for determination of transport properties since thermal conductivity, K, can be calculated from the product of density, p, specific heat, C,, and thermal diffusivity, 'h "The submitted manuscript has been authored by a contractor of the U.S. Government under contract No. DE-AC05-960R22464. Accordingly, the U.S. Government retains a nonexclusive, royalty-free license to publish or reproduce the published form of this contribution, or allow others to do so, for U.S. Government purposes."

DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, make any warranty,express or implied, or assumesany legal liability or mponsibilityfor the accuracy, completeness, or usefulwss of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily comtitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.

DISCLAIMER Portions of this document may be illegible in electronic image products. Images are pfoduced from the best available original document.

a: that is, K = ape,. The flash diffusivity technique is an ASTM standard material testing procedure[2]. Data analysis algorithms developed by Clark and Taylor[3] and Koski[4] are often used to calculate thermal diffusivity. The traditional flash diffusivity technique utilizes a single-point infrared detector. It measures the average temperature rise at the back surface of a specimen after a heat pulse. This technique requires the specimen to be homogeneous and the heating of the specimen to be uniform. Most solid state materials with the exception of composites can meet these requirements. Composites can be divided into three basic groups: dispersed, layered and fiber-reinforced composites. Thermal diffusivity of dispersed composites can be measured by the flash technique as long as the specimen is thick enough to represent the bulk material. Layered composites have been studied extensively for many years[5-71. Thermal diffusivity of two to four layered materials can be calculated from various one-dimensional heat flow models. Fiber-reinforced composites have more complex structures and the thermal diffusivity measurements are less straightforward. Fibers can be continuous or chopped and the stacking of the fibers can be periodic or random. Taylor[8] showed at least four unit cells are needed to get an accurate measurement for composites having periodic structures. In the case of aligned fibers, the volume percentage of the fibers must be very low in order to get a good coupling effect between the fibers and the matrix. In continuous fiber ceramic composites, the thermal transport property of the fibers may be very different from the matrix and the fiber volume can be high. Heat conduction can change significantly depending on the relative orientation between the heat flow and the fibers. Thermal diffusivity of fiber-reinforced composites has been measured using the traditional flash technique. Some researchers have developed models to estimate the thermal diffusivity of fiber reinforced composites[9-12]. Accurate thermal transport properties of both the matrix and the fibers must be known in order to apply the rule-of-mixtures. However, thermal properties of the material may change after processing and sintering. It is very important to be able to measure thermal diffusivity of different phases in a composite. Recently, thermal diffusivity imaging of CFCC was studied by Ahuja et al[13]. They used an IR camera to generate thermal diffusivity maps and studied artificial defects of at least 4 mm x 4 mm in dimension. In the present study, a high-resolution infrared camera was used to create thermal diffusivity maps of carbon-carbon composites. Unlike the single-point IR detector, the camera has a 256 x 256 pixel InSb Focal Plane Array (FPA). Each pixel in the FPA is an individual IR detector. The camera can operate in a snapshot mode where all the pixels are exposed to IR radiation simultaneously for the same amount of time. This feature allows the camera to be used as a 2-D IR detector in the flash technique. The camera can take up to 200 images at a rate of

120Hz. Thermal diffusivity of each pixel was calculated and formed a thermal diffusivity map with 20 pm spatial resolution. Four-direction (4D)carbon-carbon composites made by Fiber Materials Inc. with Amoco Performance Products P-25 fibers were tested in this study. A preliminary study[l4] showed the possibility of generating diffusivity maps for the carbon-carbon composites. More powerful image analyzing software and better spatial resolution were utilized in the present study. EXPERDENTAL4 An infrared camera manufactured by Amber, Goleta CA, was used as the detector described in the flash technique. Figure 1 shows the experimental setup for diffusivity mapping. A typical specimen was a disk about 12.5mm in diameter. It was placed in a stainless steel sample holder and mounted in front of a high power xenon flash lamp. A low thermal conductivity clay was used to hold the specimen and minimize heat loss via the sample holder. The IR camera monitored the back surface of the specimen. A photo diode was used to detect the flash and send a trigger signal to the control workstation. The camera took 200 images at a rate of 120 Hz following the trigger. The images were stored in the computer. An Excelm Macro program was written to calculate thermal diffusivity pixel by pixel. Spatial resolution of 20 pm was used (7.5 pm resolution can be achieved using a microscope lens with a field of view of 1 mm x 1 mm). Parker's technique[l] was used to calculate the thermal diffusivity. It assumes no heat loss and measures the time for the back surface to reach half of the maximum temperature, to,5. Knowing the specimen thickness, d, thermal diffusivity is simply calculated as: a = 0.1388d2/b,5 where a is the thermal diffusivity in cm2/s. 'The results were stored in a spreadsheet and can be easily plotted by graphic software packages. A 4D coarse weaved carbon-carbon composite sample, from Fiber Materials Inc., 12.5 mm in diameter and 14 mm thick was used for thermal diffusivity mapping. There were 4O00 fibers in a bundle. The composite was produced from P-25 pitch precursor carbon fibers woven in three directions in the horizontal plane, 120 degrees apart with the fourth fiber direction perpendicular to the plane. The overall fiber volume was 48%. The carbon fiber preform was densified through a series of repetitive impregnation, pressure carbonization and graphitization cycles until a density of 2.02 g/crn3was achieved. The inpregant was a coal tar based pitch material manufactured by Allied Signal which infiltrated the spaces between the

0~~~ Xenon Flash Lamp or IR Detector ~ Specimen Workstation 00 0. Power Supply Figure 1. Flash diffusivity measurement using an IR detector or camera. fibers as a liquid, and was converted to pure carbon by heating to 65OOC. The resulting carbon reinforced carbon block was heat treated at 225OOC to graphitize the carbon. The specimen was cut parallel to the horizontal planes. The direction of heat flow was parallel to the fourth set of fiber bundles. RESULTS AND DISCUSSION The rapid heating of carbon-carbon composite specimens was captured by the IR camera. Figure 2 is an R image 32 milliseconds after the heat pulse. 'Ihe background was subtracted to reduce the noise level. The fiber bundles in the same direction as the heat flow showed higher temperatures which indicated they heated up faster than the matrix and fiber bundles in the horizontal planes. The high temperature sensitivity of the camera also showed temperature differences between the matrix materials and the horizontal fiber bundles. The horizontal fiber bundles showed slightly higher heating than the matrix. An intensity histogram of 65,536 pixels in the same image is shown in Figure 3.?he first peak at zero intensity comes from the black area representing the sample holder. The second broad peak at intensity level 16 represents the matrix and horizontal fiber bundles. A third peak

Figure 2. IR image 32ms after the heat pulse.fiber bundles are 1.5 mm in diameter. -a v).-nx 6000 0 L. 2E a 4000 2 2000 v - 0 \ right Areas yl. ~ 5 1 0 1 5 2 0 2 5 3 0 intensity Unit Figure 3. Intensity histogram of Figure 2. Three peaks can be found in the plot.

Di 1 Figure 4. Thermal diffusivity map of a 4D carbon-carbon composite. The grey matrix area has an average diffusivity of 1.20 cm% and the bright fiber bundles have average diffusivity of 1.50 cm2/s. The black edge is the sample holder. at intensity level 22 appeared as a shoulder. It represents the bright fiber bundles parallel to the heat flow. Using the whole sequence of 200 images, thermal diffusivity values of each pixel were calculated. A 256 x 256 pixel thermal diffusivity map is shown as a surface plot in Figure4. Thermal diffusivity of the fiber bundles in the same direction as the heat flow averages 1.50 & 0.17 cm2/s and the average thermal diffusivity of the fiber bundles in the horizontal plane and the matrix was 1.20 & 0.15 cm2/s. Averaging was done via the imaging software. The average thermal diffusivity of the entire specimen was 1.30 cm2/s. The same 14 mm thick specimen was also tested using a single-point InSb infrared detector. The measured half rise time was 132 5 ms and the thermal diffusivity calculated using Parker s method[l] was 1.52 2 0.01 cm2/s. This result

'. " indicated that the average thermal diffusivity measured by the single-point detector tends to be closer to the high conductivity fiber bundle value. Therefore, it did not represent the true average thermal diffusivity value of the specimen. The singlepoint IR detector would respond first to the temperature rise of the fiber bundles and subsequently to the horizontal fiber bundles and the matrix. Thus the singlepoint detector output is a summation of several rise and fall signals with the fastest and highest rise due to the fiber bundles parallel to the direction of heat flow. The high fraction of the fiber bundle in the sample, nearly 2596, may contribute to the high diffusivity value. The thermal diffusivity map provided more information on the carbon-carbon composite than a single diffusivity value. The coupling effect between the matrix and the fiber bundle showed as a gradual change in diffusivity at the interface. The center of the fiber bundles showed the highest diffusivity. The as-received P-25 fibers had low thermal conductivity, 25 W/mK, to start with. During high temperature graphtization the microstructure of the fibers changed and thermal conductivity increased dramatically. The fiber bundles in the 4D carbon-cabbon composite had thermal conductivities ranging from 170 to 200 W/mK. The carbon-carbon composite was chosen as an example to demonstrate the thermal diffusivity mapping technique. The ability to generate detailed thermal diffusivity maps will certainly be beneficial to the design of high temperature materials. SUMMARY Thermal images of heat conduction through 4D carbon-carbon composites were captured by an IR camera. Fiber bundles orientaed parallel to the heat flow were found conducting heat faster than the matrix. A thermal diffusivity map of the 4D carbon-carbon composite was generated. This technique has shown promising applications in characterizing thermal transport properties of composite materials. ACKNOWLEDGMENT This work was supported by the U.S. Department of Energy, Assistant Secretary for Energy Efficiency and Renewable Energy, Office of Transportation Technologies, as part of the High Temperature Materials Laboratory User Program under contract DE-AC05-960R22464, managed by Lockheed Martin Energy Research Corporation. This research work was also supported in part by a Postdoctoral Research Associates Program administered jointly by ORISE and ORNL.

REFERENCES [I] W. J. Parker, R. J. Jenkins, C.P. Butler and G.L. Abbott, Flash Method of Determining Thermal Diffusivity, Heat Capacity and Thermal Conductivity, J. Appl. Phy. 32(9) 1679 (1961) [2] ASTM Designation: E 1461-92, Standard Test Method for Thermal Diffusivity of Solids by the Flash Method, Annual Book of ASTM Standards, Vol. 14.02, 933-940 (1992) [3] L.M. Clark and R.E. Taylor, Radiation Loss in the Flash Method for Thermal Diffusivity, J. Appl. Phys., Vo1.46 No.2, 7 14-7 19 (1975) [4] J. A. Koski, Improved Data Reduction Methods for Laser Pulse Diffusivity Determination with the Use of Minicomputers, Proceeding of the 8th Symposium of Thermophysical Properties, Vol. II:94-103 (1981) [5] C.S. Ang, H.S. Tan and S.L. Chan, Three-layer Thermal Diffusivity Problem Applied to Measurements on Mercury, J. Appl. Phys., Vo1.44 No.2, 687-691 (1973) [6] R.F. Bulmer and R. Taylor, Measurement by the Flash Method of Thermal Diffusivity in Two-layer Composite Samples, High Temperatures-High Pressures, Vol. 6 491-497 (1974) [7] T. Ishiguro, A. M h o, N. Araki and N. Noda, Transient Temperature Response in Functionally Gradient Materials, International Journal of Thermophysics, Vol. 14 No. 1 101-121 (1993) [8] R. E. Taylor, J. Jortner and H. Groot, Thermal Diffusivity of Fiber-reinforced Composites Using Laser Flash Technique, Carbon 23,215-222 (1985) [9] K.K. Chawla, Ceramic Matrix Composites, Chapman and Hall, 216-218 (1993) [ 101 E. Behren, Thermal Conductivity of Composite Materials, J. Composite Materials, (2)2 1968 [ 111 D.P.H. Hasselman and L.F. Johnson, J. Composite Materials, 27 (508) 1987 [12] H. Bhatt, K.Y. Donaldson, D.P.H. Hasselman and R.T. Bhatt, Effect of Fiber Orientation on the Thermal Conductivity of a Uniaxial Carbon-reinforced Aluminoborosilicate Glass-matrix Composite for Various Specimen Geometries, J. Mater. Sci. (27) 6653 (1992) [13] S. Ahuja, W.A. Ellingson, J. S. Steckenrider and S. King, Thermal Diffusivity Imaging of Continuous-Fiber Ceramic Composite (CFCC) Materials and Components, Thermal Conductivity 23, Ed. by R.B. Dinwiddie et al., published by Technomic, 31 1-321 (1996) [14] H.Wang and R.B. Dinwiddie, Advanced Thermal Imaging of Composite Materials, Advances in Ceramic-Matrix Composites 111, Ceramic Transactions Vol. 74, (1996)