Research Papers: Conduction

Simultaneous Estimation of Principal Thermal Conductivities of an Anisotropic Composite Medium: An Inverse Analysis

[+] Author and Article Information
C. Balaji

e-mail: balaji@iitm.ac.in

S. P. Venkateshan

Heat Transfer and Thermal Power Laboratory,
Department of Mechanical Engineering,
Indian Institute of Technology Madras,
Chennai 600 036, India

V. Ramakrishnan

Thermal Properties Measurement Section, TSG,
ISAC, Indian Space Research Organization,
Bangalore 560 017, India

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received March 23, 2012; final manuscript received July 27, 2012; published online January 4, 2013. Assoc. Editor: Frank Cunha.

J. Heat Transfer 135(2), 021301 (Jan 04, 2013) (8 pages) Paper No: HT-12-1129; doi: 10.1115/1.4007422 History: Received March 23, 2012; Revised July 27, 2012

This paper reports the results of an experimental study to determine the principal thermal conductivities (kx,ky, and kz) of an anisotropic composite medium using an inverse heat transfer analysis. The direct problem consists of solving the three dimensional heat conduction equation in an orthotropic composite medium with the finite difference method to generate the required temperature distribution for known thermal conductivities. The measurement technique involves dissipating a known heat flux at the central region of a square sample and allowing it to conductively transfer the heat to an aluminium cold plate sink via a square copper ring. At steady state, temperatures at 28 (19 are used for retrievals due to symmetry) discrete locations are logged and used for parameter estimation. The entire measurement process is conducted in a vacuum environment. The inverse heat conduction problem (IHCP) for retrieving the orthotropic thermal conductivity tensor(parameter estimation) is then solved using a two layer feed forward back propagation artificial neural network (ANN) trained using the Levenberg–Marquardt algorithm (LMA), with temperatures as input and thermal conductivity values kx,ky, and kz as the output. The method is first validated against a stainless steel(SS-304) sample of known thermal properties followed by the determination of the orthotropic conductivities of the honeycomb composite material.

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


Beck, J., St. Clair, C., and Blackwell, B., 1985, Inverse Heat Conduction, John Wiley and Sons Inc., New York.
Sawaf, B., Ozisik, M., and Jarny, Y., 1995, “An Inverse Analysis to Estimate Linearly Temperature Dependent Thermal Conductivity Components and Heat Capacity of an Orthotropic Medium,” Int. J. Heat Mass Transfer, 38(16), pp. 3005–3010. [CrossRef]
Sawaf, B., and Ozisik, M., 1995, “Determining the Constant Thermal Conductivities of Orthotropic Materials by Inverse Analysis,” Int. Commun. Heat Mass Transfer, 22(2), pp. 201–211. [CrossRef]
Thomas, M., Boyard, N., Lefèvre, N., Jarny, Y., and Delaunay, D., 2010, “An Experimental Device for the Simultaneous Estimation of the Thermal Conductivity 3-D Tensor and the Specific Heat of Orthotropic Composite Materials,” Int. J. Heat Mass Transfer, 53(23–24), pp. 5487–5498. [CrossRef]
Deiveegan, M., Balaji, C., and Venkateshan, S., 2006, “Comparison of Various Methods for Simultaneous Retrieval of Surface Emissivities and Gas Properties in Gray Participating Media,” ASME J. Heat Transfer, 128, pp. 829–837. [CrossRef]
Deng, S., and Hwang, Y., 2006, “Applying Neural Networks to the Solution of Forward and Inverse Heat Conduction Problems,” Int. J. Heat Mass Transfer, 49(25–26), pp. 4732–4750. [CrossRef]
Krejsa, J., Woodbury, K., Ratliff, J., and Raudensky, M., 1999, “Assessment of Strategies and Potential for Neural Networks in the Inverse Heat Conduction Problem,” Inverse Probl. Sci. Eng., 7(3), pp. 197–213. [CrossRef]
Mera, N., Elliott, L., Ingham, D., and Lesnic, D., 2001, “Use of the Boundary Element Method to Determine the Thermal Conductivity Tensor of an Anisotropic Medium,” Int. J. Heat Mass Transfer, 44(21), pp. 4157–4167. [CrossRef]
Dong, C., Sun, F., and Meng, B., 2007, “A Method of Fundamental Solutions for Inverse Heat Conduction Problems in an Anisotropic Medium,” Eng. Anal. Boundary Elem., 31(1), pp. 75–82. [CrossRef]
Taktak, R., Beck, J., and Scott, E., 1993, “Optimal Experimental Design for Estimating Thermal Properties of Composite Materials,” Int. J. Heat Mass Transfer, 36(12), pp. 2977–2986. [CrossRef]
ISRO, “Thermal Properties of Spacecraft Materials—Data Bank,” TPMS, Thermal Systems Group, ISAC, ISRO, Bangalore, India.
Venkateshan, S. P., 2008, Mechanical Measurements, Ane Books, New Delhi, India.
Samarjeet, C., Yenni, G. R., Ambirajan, A., Balaji, C., and Venkateshan, S. P., 2011, “An Inverse Analysis for Estimation of Thermal Conductivity of Orthotropic Composite Medium Using Artificial Neural Network,” Proceedings of the 21st National and 10th ISHMT-ASME Heat and Mass Transfer Conference, December 27–30, IIT Madras, India.
Blackwell, B., Gill, W., Dowding, K., and Voth, T., 2000, “Determination of Thermal Conductivity of 304 Stainless Steel Using Parameter Estimation Techniques,” 34th National Heat Transfer Conference–NHTC, Pittsburgh, PA, Aug. 20–22.
Nishino, K., Yamashita, S., and Torii, K., 1995, “Thermal Contact Conductance Under Low Applied Load in a Vacuum Environment,” Exp. Therm. Fluid Sci., 10(2), pp. 258–271. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic of the problem geometry

Grahic Jump Location
Fig. 2

Simulated temperature distribution on a line along z-axis through sensor no. 3 for q = 5 W, kx=5 W/mK,ky=4 W/mK and varying kz,h=1400 W/m2K and sink temperature of −60 °C

Grahic Jump Location
Fig. 3

Schematic view of the complete experimental setup

Grahic Jump Location
Fig. 4

Exploded view of the test sample assembly

Grahic Jump Location
Fig. 5

Thermocouple sensor locations in the test sample. Numbers in parentheses show sensors at the bottom of the sample.

Grahic Jump Location
Fig. 6

Sensitivity of simulated temperatures to contact conductance (h) in one quarter of test sample, for honeycomb material, at q = 2.5 W

Grahic Jump Location
Fig. 7

Parity plot showing adequacy of the ANN in accomplishing the retrievals for honeycomb sample, at q = 2.5 W



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