Research Papers: Conduction

Thermal Conductivity of Turbostratic Carbon Nanofiber Networks

[+] Author and Article Information
Matthew L. Bauer, Christopher B. Saltonstall, Patrick E. Hopkins, Pamela M. Norris

Department of Mechanical and
Aerospace Engineering,
University of Virginia,
Charlottesville, VA 22904-4746

Zayd C. Leseman

Department of Mechanical and Industrial Engineering,
University of New Mexico,
Albuquerque, NM 87131

Thomas E. Beechem

Sandia National Laboratories,
Albuquerque, NM 87185

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received April 8, 2015; final manuscript received November 27, 2015; published online March 15, 2016. Assoc. Editor: Alan McGaughey.The United States Government retains, and by accepting the article for publication, the publisher acknowledges that the United States Government retains, a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for United States government purposes.

J. Heat Transfer 138(6), 061302 (Mar 15, 2016) (9 pages) Paper No: HT-15-1272; doi: 10.1115/1.4032610 History: Received April 08, 2015; Revised November 27, 2015

Composite material systems composed of a matrix of nanomaterials can achieve combinations of mechanical and thermophysical properties outside the range of traditional systems. The microstructure of the system dictates the rate, in which heat moves through the material. In this work, air/carbon nanofiber networks are studied to elucidate the system parameters influencing thermal transport. Thermal properties are measured with varying initial carbon fiber fill fraction, environment pressure, loading pressure, and heat treatment temperature (HTT) through a bidirectional modification of the 3ω technique. The nanostructure of the individual fibers is characterized with small angle X-ray scattering and Raman spectroscopy providing insight to individual fiber thermal conductivity. Measured thermal conductivity of the carbon nanofiber networks varied from 0.010 W/(m K) to 0.070 W/(m K). An understanding of the intrinsic properties of the individual fibers and the interactions of the two-phase composite is used to reconcile low measured thermal conductivities with predictive modeling. Accounting for fiber-to-fiber interactions and the nuanced changes in the composite as pressure is applied is necessary to successfully model thermal transport in system.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Peebles, L. , Jr., 1995, Carbon Fibers: Formation, Structure, and Properties, CRC Press, Boca Raton, FL.
Heremans, J. , Rahim, I. , and Dresselhaus, M. S. , 1985, “ Thermal Conductivity and Raman Spectra of Carbon Fibers,” Phys. Rev. B, 32(10), pp. 6742–6747. [CrossRef]
Mayhew, E. , and Prakash, V. , 2013, “ Thermal Conductivity of Individual Carbon Nanofibers,” Carbon, 62, pp. 493–500. [CrossRef]
Qiu, L. , Zheng, X. , Zhu, J. , Su, G. , and Tang, D. , 2013, “ The Effect of Grain Size on the Lattice Thermal Conductivity of an Individual Polyacrylonitrile-Based Carbon Fiber,” Carbon, 51, pp. 265–273. [CrossRef]
Piraux, L. , Nysten, B. , Haquenne, A. , Issi, J.-P. , Dresselhaus, M. , and Endo, M. , 1984, “ The Temperature Variation of the Thermal Conductivity of Benzene-Derived Carbon Fibers,” Solid State Commun., 50(8), pp. 697–700. [CrossRef]
Chung, D. , 2001, “ Materials for Thermal Conduction,” Appl. Therm. Eng., 21(16), pp. 1593–1605. [CrossRef]
Atwater, M. A. , Mousavi, A. K. , Leseman, Z. C. , and Phillips, J. , 2013, “ Direct Synthesis and Characterization of a Nonwoven Structure Comprised of Carbon Nanofibers,” Carbon, 57, pp. 363–370. [CrossRef]
Atwater, M. A. , Phillips, J. , Doorn, S. K. , Luhrs, C. C. , Fernández, Y. , Menéndez, J. , and Leseman, Z. C. , 2009, “ The Production of Carbon Nanofibers and Thin Films on Palladium Catalysts From Ethylene–Oxygen Mixtures,” Carbon, 47(9), pp. 2269–2280. [CrossRef]
Tuinstra, F. , and Koenig, J. L. , 1970, “ Raman Spectrum of Graphite,” J. Chem. Phys., 53(3), pp. 1126–1130. [CrossRef]
Ferrari, A. C. , and Robertson, J. , 2000, “ Interpretation of Raman Spectra of Disordered and Amorphous Carbon,” Phys. Rev. B, 61(20), pp. 14095–14107. [CrossRef]
Fourdeux, A. , Herinckx, C. , Perret, R. , and Ruland, W. , 1970, “ La structure des fibres de carbone,” Compt. Rend. Acad. Sci., 25(1969), pp. 1507–1600.
Gupta, A. , Harrison, I. R. , and Lahijani, J. , 1994, “ Small-Angle X-Ray Scattering in Carbon Fibers,” J. Appl. Crystallogr., 27(4), pp. 627–636. [CrossRef]
Shioya, M. , and Takaku, A. , 1985, “ Characterization of Microvoids in Carbon Fibers by Absolute Small-Angle X-Ray Measurements on a Fiber Bundle,” J. Appl. Phys., 58(11), pp. 4074–4082. [CrossRef]
Kaganer, M. G. , 1969, Thermal Insulation in Cryogenic Engineering, Israel Program for Scientific Translations, Jerusalem, Israel.
Bauer, M. L. , and Norris, P. M. , 2014, “ General Bidirectional Thermal Characterization Via the 3ω Technique,” Rev. Sci. Instrum., 85(6), p. 064903.
Hu, X. , Padilla, A. , Xu, J. , Fisher, T. , and Goodson, K. , 2006, “ 3-Omega Measurements of Vertically Oriented Carbon Nanotubes on Silicon,” ASME J. Heat Transfer, 128(11), pp. 1109–1113. [CrossRef]
Qiu, L. , Tang, D. W. , Zheng, X. H. , and Su, G. P. , 2011, “ The Freestanding Sensor-Based 3ω Technique for Measuring Thermal Conductivity of Solids: Principle and Examination,” Rev. Sci. Instrum., 82(4), p. 045106. [CrossRef] [PubMed]
Lee, S.-M. , and Cahill, D. , 1997, “ Heat Transport in Thin Dielectric Films,” J. Appl. Phys., 81(6), pp. 2590–2595. [CrossRef]
Kim, J. H. , Feldman, A. , and Novotny, D. , 1999, “ Application of the Three Omega Thermal Conductivity Measurement Method to a Film on a Substrate of Finite Thickness,” J. Appl. Phys., 86(7), pp. 3959–3963. [CrossRef]
Sun, W.-C. , Huang, M.-J. , Chien, H.-C. , Chang, T.-Y. , and Yao, D.-J. , 2010, “ A Novel Method for Measuring Thick Film Thermal Conductivity,” 5th IEEE International Conference on Nano/Micro-Engineered and Molecular Systems (NEMS), Xiamen, China, Jan. 20–23, pp. 1052–1056.
Feldman, A. , 1999, “ Algorithm for Solutions of the Thermal Diffusion Equation in a Stratified Medium With a Modulated Heating Source,” High Temp. High Pressures, 31(3), pp. 293–298. [CrossRef]
Cahill, D. , 1990, “ Thermal Conductivity Measurement From 30 to 750 K: The 3ω Method,” Rev. Sci. Instrum., 61(2), pp. 802–808. [CrossRef]
Schiffres, S. N. , Kim, K. H. , Hu, L. , McGaughey, A. J. H. , Islam, M. F. , and Malen, J. A. , 2012, “ Gas Diffusion, Energy Transport, and Thermal Accommodation in Single-Walled Carbon Nanotube Aerogels,” Adv. Funct. Mater., 22(24), pp. 5251–5258. [CrossRef]
Prasher, R. S. , Hu, X. J. , Chalopin, Y. , Mingo, N. , Lofgreen, K. , Volz, S. , Cleri, F. , and Keblinski, P. , 2009, “ Turning Carbon Nanotubes From Exceptional Heat Conductors Into Insulators,” Phys. Rev. Lett., 102(10), p. 105901. [CrossRef] [PubMed]
Bhattacharyya, R. K. , 1980, Heat-Transfer Model for Fibrous Insulations, Thermal Insulation Performance, American Society for Testing and Materials, Philadelphia, PA, pp. 272–286.
General Electric Company, Research and Development Center, 1984, Heat Transfer and Fluid Flow Data Books, Genium Publishing, Schenectady NY.
Raed, K. , and Gross, U. , 2009, “ Modeling of Influence of Gas Atmosphere and Pore-Size Distribution on the Effective Thermal Conductivity of Knudsen and Non-Knudsen Porous Materials,” Int. J. Thermophys., 30(4), pp. 1343–1356. [CrossRef]
Yasumoto, I. , 1987, “ Accommodation Coefficients of Helium, Neon, Argon, Hydrogen, and Deuterium on Graphitized Carbon,” J. Phys. Chem., 91(16), pp. 4298–4301. [CrossRef]
Hu, L. , and McGaughey, A. J. H. , 2013, “ Energy Accommodation Between Noble Gases and Carbon Nanotubes,” J. Phys. Chem. C, 117(37), pp. 18804–18808. [CrossRef]
Amdur, I. , and Guildner, L. A. , 1957, “ Thermal Accommodation Coefficients on Gas-Covered Tungsten, Nickel and Platinum,” J. Am. Chem. Soc., 79(2), pp. 311–315. [CrossRef]
Rader, D. J. , Grasser, T. W. , Castaneda, J. N. , Trott, W. M. , and Torczynski, J. R. , 2005, “ Measurements of Thermal Accommodation Coefficients,” U.S. Department of Energy, Springfield, VA, Report No. SAND2005-6084.
Fletcher, L. , and White, F. , 1984, Heat Transfer and Fluid Flow Data Books, Genium Publishing, Schnectady, NY.
Fricke, H. , 1924, “ A Mathematical Treatment of the Electric Conductivity and Capacity of Disperse Systems—I: The Electric Conductivity of a Suspension of Homogeneous Spheroids,” Phys. Rev., 24(5), pp. 575–587. [CrossRef]
Stark, C. , and Fricke, J. , 1993, “ Improved Heat-Transfer Models for Fibrous Insulations,” Int. J. Heat Mass Transfer, 36(3), pp. 617–625. [CrossRef]
Evans, W. J. , Shen, M. , and Keblinski, P. , 2012, “ Inter-Tube Thermal Conductance in Carbon Nanotubes Arrays and Bundles: Effects of Contact Area and Pressure,” Appl. Phys. Lett., 100(26), p. 261908. [CrossRef]
Yang, J. , Waltermire, S. , Chen, Y. , Zinn, A. A. , Xu, T. T. , and Li, D. , 2010, “ Contact Thermal Resistance Between Individual Multiwall Carbon Nanotubes,” Appl. Phys. Lett., 96(2), p. 023109.


Grahic Jump Location
Fig. 1

Average Raman spectra of each sample. Each spectrum was first normalized to the G peak intensity and then averaged. (a) The averaged spectra of the unannealed samples. (b) The averaged spectra of the annealed S1 series.

Grahic Jump Location
Fig. 2

Guinier plot of SAXS intensity decay. The lack of linear region over this q regime indicates an absence of monodisperse voids observed in other carbon fiber samples.

Grahic Jump Location
Fig. 3

Scanning electron microscope image of a carbon fiber sample

Grahic Jump Location
Fig. 4

Raw data and modeled fit from sample S1 with an applied pressure of 220 kPa. The experimental uncertainty in each data point is slightly smaller than the data markers. The measured values match the inputs used in Fig. 5.

Grahic Jump Location
Fig. 5

Sensitivity of the 3ω experiment to composite samples'effective thermal conductivity and volumetric heat capacity as well as the contact resistance between the mount and sample. This plot was generated with ke = 0.024 W/(m K), Ce = 3 × 105 J/(m3 K), and Rp = 6 × 10−5 K m2/W, similar to the fourth data point in Fig. 6. The change in direction of the contact resistance sensitivity above 10,000 Hz only exists due to the presence of the SiOx layer closer to the heater. At greater applied pressures, the signal's sensitivity to the carbon fibers' thermal conductivity and heat capacity is increased. At low applied pressures, the experiment can lose sensitivity to these properties [15].

Grahic Jump Location
Fig. 6

Best fit results of data taken on S1 (without heat treatment) with measured density of 90 kg/m3. Measurements were taken at six unique compressive pressures, with the sample thermal conductivity, sample volumetric heat capacity, and the sample to mount contact resistance simultaneously fit. This data are representative of data taken on each sample. The fitting error propagated from sources of uncertainty in the system is estimated at 9% for the sample's thermal conductivity, 6% for the volumetric heat capacity, and 4% for the contact resistance between the sample and the passivation layer.

Grahic Jump Location
Fig. 7

Measured thermal conductivity values from sample S5 at both vacuum and atmospheric pressure. Carbon fill fraction is calculated from the best fit volumetric heat capacity using Eq. (2). Also, plotted is the SFM accounting for varying sample properties as pressure is applied.

Grahic Jump Location
Fig. 8

Normalized modeled response of 200 kg/m3 dense carbon nanofiber composite (S5) as pressure is applied. The total thermal conductivity response, ke, only increases due to the role of thermal conductance through fiber contacts. Alternatively, increasing carbon fiber fill fraction does not increase effective thermal transport while excluding additional fiber-to-fiber contacts. The contribution of gas conductivity decreases at greater pressures due to decreasing pore sizes. The bending over of fibers at greater applied pressures also causes a decrease in the thermal conductivity.

Grahic Jump Location
Fig. 9

Best fit 3ω thermal conductivity results for nonheat-treated samples S1, S2, S3, S4, and S5 (densities of 90, 110, 130, 160, and 200 kg/m3, respectively). All the data are acquired at room temperature and atmospheric pressure. The thermal conductivity predictions from the SFM (dashed lines) are also plotted for initial fill fractions ranging from 4% to 12% in steps of 2% (ϕ0). As pressure is applied, the fill fraction increases to a maximum value of ϕ=40%.

Grahic Jump Location
Fig. 10

Best fit data from the heat-treated sample series (S1). The measured thermal conductivity increases significantly from the as-grown sample compared to the heat-treated samples, with the sample which underwent a 1500 °C HTT showing an additional increase compared to the more moderately heated samples. The change in effective thermal conductivity with respect to carbon fiber fill fraction has increased with heat treatment. Over the data, the SFM is plotted with carbon fiber thermal conductivities of 40, 80, 120, and 160 W/(m K), respectively.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In