Research Papers: Conduction

Analytical Solution for Three-Dimensional, Unsteady Heat Conduction in a Multilayer Sphere

[+] Author and Article Information
Suneet Singh

Department of Energy Science and Engineering,
Indian Institute of Technology (IIT) Bombay,
Powai, Mumbai 400076, India
e-mail: suneet.singh@iitb.ac.in

Prashant K. Jain

Reactor and Nuclear Systems Division,
Oak Ridge National Laboratory,
Oak Ridge, TN 37831


Department of Nuclear, Plasma
and Radiological Engineering,
University of Illinois at Urbana—Champaign,
216 Talbot Lab,
104 S. Wright Street,
Urbana, IL 61801

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received July 21, 2015; final manuscript received April 25, 2016; published online June 7, 2016. Assoc. Editor: William P. Klinzing.This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The U.S. government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for U.S. government purposes.

J. Heat Transfer 138(10), 101301 (Jun 07, 2016) (11 pages) Paper No: HT-15-1488; doi: 10.1115/1.4033536 History: Received July 21, 2015; Revised April 25, 2016

An analytical solution has been obtained for the transient problem of three-dimensional multilayer heat conduction in a sphere with layers in the radial direction. The solution procedure can be applied to a hollow sphere or a solid sphere composed of several layers of various materials. In general, the separation of variables applied to 3D spherical coordinates has unique characteristics due to the presence of associated Legendre functions as the eigenfunctions. Moreover, an eigenvalue problem in the azimuthal direction also requires solution; again, its properties are unique owing to periodicity in the azimuthal direction. Therefore, extending existing solutions in 2D spherical coordinates to 3D spherical coordinates is not straightforward. In a spherical coordinate system, one can solve a 3D transient multilayer heat conduction problem without the presence of imaginary eigenvalues. A 2D cylindrical polar coordinate system is the only other case in which such multidimensional problems can be solved without the use of imaginary eigenvalues. The absence of imaginary eigenvalues renders the solution methodology significantly more useful for practical applications. The methodology described can be used for all the three types of boundary conditions in the outer and inner surfaces of the sphere. The solution procedure is demonstrated on an illustrative problem for which results are obtained.

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


Vodicka, V. , 1955, “ Eindimensionale Warmeleitung in Geschichteten Koorpern,” Math. Nachr., 14(1), pp. 47–55. [CrossRef]
Tittle, C. W. , 1965, “ Boundary Value Problems in Composite Media: Quasi-Orthogonal Functions,” J. Appl. Phys., 36(4), pp. 1486–1488. [CrossRef]
Bulavin, P. E. , and Kascheev, V. M. , 1965, “ Solution of the Non-Homogeneous Heat Conduction Equation for Multilayered Bodies,” Int. Chem. Eng., 5(1), pp. 112–115.
Mulholland, G. P. , and Cobble, M. H. , 1972, “ Diffusion Through Composite Media,” Int. J. Heat Mass Transfer, 15(1), pp. 147–160. [CrossRef]
Mikhailov, M. D. , Ozisik, M. N. , and Vulchanov, N. L. , 1983, “ Diffusion in Composite Layers With Automatic Solution of the Eigenvalue Problem,” Int. J. Heat Mass Transfer, 26(8), pp. 1131–1141. [CrossRef]
Huang, S. C. , and Chang, Y. P. , 1980, “ Heat Conduction in Unsteady, Periodic and Steady States in Laminated Composites,” ASME J. Heat Transfer, 102(4), pp. 742–748. [CrossRef]
Carslaw, H. S. , and Jaeger, J. C. , 1959, Conduction of Heat in Solids, 2nd ed., Oxford University Press, London.
Haji-Sheikh, A. , and Beck, J. V. , 1990, “ Green's Function Partitioning in Galerkin-Base Integral Solution of the Diffusion Equation,” ASME J. Heat Transfer, 112(1), pp. 28–34. [CrossRef]
Salt, H. , 1983, “ Transient Heat Conduction in a Two-Dimensional Composite Slab—I: Theoretical Development of Temperatures Modes,” Int. J. Heat Mass Transfer, 26(11), pp. 1611–1616. [CrossRef]
Salt, H. , 1983, “ Transient Heat Conduction in a Two-Dimensional Composite Slab—II: Physical Interpretation of Temperatures Modes,” Int. J. Heat Mass Transfer, 26(11), pp. 1617–1623. [CrossRef]
Mikhailov, M. D. , and Ozisik, M. N. , 1986, “ Transient Conduction in a Three-Dimensional Composite Slab,” Int. J. Heat Mass Transfer, 29(2), pp. 340–342. [CrossRef]
Yener, Y. , and Ozisik, M. N. , 1974, “ On the Solution of Unsteady Heat Conduction in Multi-Region Finite Media With Time-Dependent Heat Transfer Coefficient,” Fifth International Heat Transfer Conference, JSME, Tokyo, Vol. 1, pp. 188–192.
Wolfram Research, 2005, “ Mathematica, Version 5.2,” Wolfram Research, Champaign, IL.
de Monte, F. , 2000, “ Transient Heat Conduction in One-Dimensional Composite Slab. A ‘Natural’ Analytic Approach,” Int. J. Heat Mass Transfer, 43(19), pp. 3607–3619. [CrossRef]
de Monte, F. , 2002, “ An Analytic Approach to the Unsteady Heat Conduction Processes in One-Dimensional Composite Media,” Int. J. Heat Mass Transfer, 45(6), pp. 1333–1343. [CrossRef]
de Monte, F. , 2003, “ Unsteady Heat Conduction in Two-Dimensional Two Slab-Shaped Regions. Exact Closed-Form Solution and Results,” Int. J. Heat Mass Transfer, 46(8), pp. 1455–1469. [CrossRef]
de Monte, F. , 2004, “ Transverse Eigenproblem of Steady-State Heat Conduction for Multi-Dimensional Two-Layered Slabs With Automatic Computation of Eigenvalues,” Int. J. Heat Mass Transfer, 47(2), pp. 191–201. [CrossRef]
Haji-Sheikh, A. , and Beck, J. V. , 2002, “ Temperature Solution in Multi-Dimensional Multi-Layer Bodies,” Int. J. Heat Mass Transfer, 45(9), pp. 1865–1877. [CrossRef]
Lu, X. , Tervola, P. , and Viljanen, M. , 2005, “ A New Analytical Method to Solve Heat Equation for Multi-Dimensional Composite Slab,” J. Phys. A: Math. Gen., 38(13), pp. 2873–2890. [CrossRef]
Lu, X. , Tervola, P. , and Viljanen, M. , 2006, “ Transient Analytical Solution to Heat Conduction in Multi-Dimensional Composite Cylinder Slab,” Int. J. Heat Mass Transfer, 49(5–6), pp. 1107–1114. [CrossRef]
Singh, S. , Jain, P. K. , and Rizwan-Uddin. , 2008, “ Analytical Solution to Transient Heat Conduction in Polar Coordinates With Multiple Layers in Radial Direction,” Int. J. Therm. Sci., 47(3), pp. 261–273. [CrossRef]
Jain, P. K. , Singh, S. , and Rizwan-Uddin , 2009, “ Transient Analytical Solution to Asymmetric Heat Conduction in a Multilayer Annulus,” ASME J. Heat Transfer, 131(1), p. 011304. [CrossRef]
Singh, S. , Jain, P. K. , and Rizwan-Uddin , 2011, “ Finite Integral Transform Technique to Solve Asymmetric Heat Conduction in a Multilayer Annulus With Time Dependent Boundary Conditions,” Nucl. Eng. Des., 241(1), pp. 144–154. [CrossRef]
Jain, P. K. , Singh, S. , and Rizwan-Uddin , 2006, “ An Exact Analytical Solution for Two-Dimensional, Unsteady, Multilayer Heat Conduction in Spherical Coordinates,” Int. J. Heat Mass Transfer, 53(9–10), pp. 2133–2142.
Ozisik, M. N. , 1993, Heat Conduction, Wiley, New York.
Byerly, W. E. , 1893, An Elementary Treatise in Fourier's Series and Spherical Cylindrical and Ellipsoidal Harmonics With Applications to Problems in Mathematical Physics, Ginn & Company, New York.
Antonopoulos, K. A. , and Tzivanidis, C. , 1996, “ Analytical Solution of Boundary Value Problems of Heat Conduction in Composite Regions With Arbitrary Convection Boundary Conditions,” Acta Mech., 118(1), pp. 65–78. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Coordinate system for the problem and (b) 2D cross-sectional view of n concentric spherical layers at θ  = π/2 plane

Grahic Jump Location
Fig. 2

The plot of heat flux applied at the surface of the sphere for the illustrative example given by Eq. (45)

Grahic Jump Location
Fig. 3

Plot of eigencondition f[λ] against λ for (a) m = 1 and (b) m = 2. Roots of the transcendental eigencondition are obtained by using Mathematica [27] and are graphically verified by plotting f[λ].

Grahic Jump Location
Fig. 4

Transient temperature distribution in the radial direction along θ = 0. The change of azimuthal coordinate will not have any change in the values at θ = 0.

Grahic Jump Location
Fig. 5

Transient temperature distribution in the radial direction along θ=π/4 for different azimuthal (ϕ) locations: (a) ϕ=0, (b) ϕ=π/2, and (c) ϕ=3π/2

Grahic Jump Location
Fig. 6

Transient temperature distribution in the radial direction along θ=π/2 for different azimuthal (ϕ) locations: (a) ϕ=0, (b) ϕ=π/2, and (c) ϕ=3π/2

Grahic Jump Location
Fig. 7

Temperature contours in the θ=π/2 midplane of the three-layer sphere at different times. Coordinate system is same as in Fig. 1(b).



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