Research Papers: Conduction

Efficient Numerical Evaluation of Exact Solutions for One-Dimensional and Two-Dimensional Infinite Cylindrical Heat Conduction Problems

[+] Author and Article Information
Te Pi

Mechanical & Materials Engineering Department,
University of Nebraska Lincoln,
W342C NH,
Lincoln, NE 68588-0526
e-mail: alitem0829@gmail.com

Kevin Cole

Mechanical & Materials Engineering Department,
University of Nebraska Lincoln,
W342C NH,
Lincoln, NE 68588-0526
e-mail: kcole1@unl.edu

James Beck

Professor Emeritus
College of Engineering,
Michigan State University,
Engineering Building, 428 S. Shaw Lane,
East Lansing, MI 48824-1226
e-mail: jameserebeck@gmail.com

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received September 14, 2016; final manuscript received May 9, 2017; published online July 25, 2017. Assoc. Editor: Alan McGaughey.

J. Heat Transfer 139(12), 121301 (Jul 25, 2017) (10 pages) Paper No: HT-16-1582; doi: 10.1115/1.4037081 History: Received September 14, 2016; Revised May 09, 2017

Estimation of thermal properties or diffusion properties from transient data requires that a model is available that is physically meaningful and suitably precise. The model must also produce numerical values rapidly enough to accommodate iterative regression, inverse methods, or other estimation procedures during which the model is evaluated again and again. Bodies of infinite extent are a particular challenge from this perspective. Even for exact analytical solutions, because the solution often has the form of an improper integral that must be evaluated numerically, lengthy computer-evaluation time is a challenge. The subject of this paper is improving the computer evaluation time for exact solutions for infinite and semi-infinite bodies in the cylindrical coordinate system. The motivating applications for the present work include the line-source method for obtaining thermal properties, the estimation of thermal properties by the laser-flash method, and the estimation of aquifer properties or petroleum-field properties from well-test measurements. In this paper, the computer evaluation time is improved by replacing the integral-containing solution by a suitable finite-body series solution. The precision of the series solution may be controlled to a high level and the required computer time may be minimized, by a suitable choice of the extent of the finite body. The key finding of this paper is that the resulting series may be accurately evaluated with a fixed number of terms at any value of time, which removes a long-standing difficulty with series solution in general. The method is demonstrated for the one-dimensional case of a large body with a cylindrical hole and is extended to two-dimensional geometries of practical interest. The computer-evaluation time for the finite-body solutions are shown to be hundreds or thousands of time faster than the infinite-body solutions, depending on the geometry.

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


ASME, 2009, “ Standard for Verification and Validation in Computational Fluid Dynamics and Heat Transfer,” American Society of Mechanical Engineers, New York, Standard No. ASME V V 20-2009.
McMasters, R. L. , Dowding, K. J. , Beck, J. V. , and Yen, D. , 2002, “ Methodology to Generate Accurate Solutions for Verification in Transient Three-Dimensional Heat Conduction,” J. Numer. Heat Transfer, Part B, 41(6), pp. 521–541. [CrossRef]
Beck, J. V. , McMasters, R. L. , Dowding, K. J. , and Amos, D. E. , 2006, “ Intrinsic Verification Methods in Linear Heat Conduction,” Int. J. Heat Mass Transfer, 49(17–18), pp. 2984–2994. [CrossRef]
Tian, T. , and Cole, K. D. , 2012, “ Anisotropic Thermal Conductivity Measurements of Carbon-Fiber/Epoxy Composites,” Int. J. Heat Mass Transfer, 55(23–24), pp. 6530–6537. [CrossRef]
Guo, J. , Wang, X. , and Wang, T. , 2007, “ Thermal Characterization of Microscale Conductive and Nonconductive Wires Using Transient Electrothermal Technique,” J. Appl. Phys., 101(6), p. 063537. [CrossRef]
Woodbury, K. A. , Beck, J. V. , and Najafi, H. , 2014, “ Filter Solution of Inverse Heat Conduction Problem Using Measured Temperature History as Remote Boundary Condition,” Int. J. Heat Mass Transfer, 72, pp. 139–147. [CrossRef]
Beck, J. V. , Blackwell, B. , and Clair, C. R. S., Jr. , 1985, Inverse Heat Conduction—III: Posed Problems, Wiley, New York.
Busswell, G. S. , Banerjee, R. , Thambynayagam, R. M. , and Spath, J. B. , 2006, “ Generalized Analytical Solution for Reservoir Problems With Multiple Wells and Boundary Conditions,” Intelligent Energy Conference and Exhibition, Amsterdam, The Netherlands, Apr. 11–13, SPE Paper No. SPE-99288.
Cheng, W.-L. , Huang, Y.-H. , Lu, D.-T. , and Yin, H.-R. , 2011, “ A Novel Analytical Transient Heat-Conduction Time Function for Heat Transfer in Steam Injection Wells Considering the Wellbore Heat Capacity,” Energy, 36(7), pp. 4080–4088. [CrossRef]
Yeh, H.-D. , Yang, S.-Y. , and Peng, H.-Y. , 2003, “ A New Closed-Form Solution for a Radial Two-Layer Drawdown Equation for Groundwater Under Constant-Flux Pumping in a Finite-Radius Well,” Adv. Water Resour., 26(7), pp. 747–757. [CrossRef]
Clow, G. D. , 2015, “ A Green’s Function Approach for Assessing the Thermal Disturbance Caused by Drilling Deep Boreholes in Rock or Ice,” Geophys. J. Int., 203(3), pp. 1877–1895. [CrossRef]
McMasters, R. L. , Dinwiddie, R. B. , and Haji-Sheikh, A. , 2007, “ Estimating the Thermal Conductivity of a Film on a Known Substrate,” J. Thermophys. Heat Transfer, 21(4), pp. 681–687. [CrossRef]
Zeng, H. Y. , Diao, N. R. , and Fang, Z. H. , 2002, “ A Finite Line-Source Model for Boreholes in Geothermal Heat Exchangers,” Heat Transfer Asian Res., 31(7), pp. 558–567. [CrossRef]
Lamarche, L. , and Beauchamp, B. , 2007, “ A New Contribution to the Finite Line-Source Model for Geothermal Boreholes,” Energy Build., 39(2), pp. 188–198. [CrossRef]
Healy, J. J. , de Groot, J. J. , and Kestin, J. , 1976, “ The Theory of the Transient Hot-Wire Method for Measuring Thermal Conductivity,” Physica B+C, 82(2), pp. 392–408. [CrossRef]
Beck, J. V. , and Litkouhi, B. , 1988, “ Heat Conduction Number System,” Int. J. Heat Mass Transfer, 31(3), pp. 505–515. [CrossRef]
Cole, K. D. , Woodbury, K. , Amos, D. E. , Beck, J. V. , Crittenden, P. E. , De Monte, F. , and Haji-Sheikh, A. , Guimaraes, G. , McMasters, R. , and Roberty, N. , 2017, “ EXACT Analytical Conduction Toolbox,” University of Nebraska–Lincoln, Lincoln, NE, accessed Oct. 4, 2015, http://exact.unl.edu
Cole, K. D. , Beck, J. V. , Haji-Sheikh, A. , and Litkouhi, B. , 2011, Heat Conduction Using Green’s Functions, 2nd ed., CRC Press, Boca Raton, FL.
Press, W. H. , Teukolsky, S. A. , Vetterling, W. T. , and Flannery, B. P. , 2007, Numerical Recipes: The Art of Scientific Computing, 3rd ed., Cambridge University Press, Cambridge, UK.
de Monte, F. , Beck, J. V. , and Amos, D. E. , 2012, “ Solving Two-Dimensional Cartesian Unsteady Heat Conduction Problems for Small Values of the Time,” Int. J. Therm. Sci., 60, pp. 106–113. [CrossRef]
Cole, K. D. , 2004, “ Fast-Converging Series for Heat Conduction in the Circular Cylinder,” J. Eng. Math., 49(3), pp. 217–232. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Geometry for R20 with R1<r<∞ and (b) geometry for R22 with R1<r<R2

Grahic Jump Location
Fig. 2

Plot integrand of R20’s solution F(β,t̃) with t̃=0.5,1,5,10

Grahic Jump Location
Fig. 3

Plot of integrand in improper integral from R20’s solution F¯(z,t̃) for 1/b = 10 with t̃=0.5,1,5,10

Grahic Jump Location
Fig. 4

(a) R20Z22 geometry model with radial domain (R1<r<∞) and (b) R21Z22 geometry model with radial domain (R1<r<R2)

Grahic Jump Location
Fig. 5

(a) Dimensionless temperature of R20Z22 case versus z̃ at t̃=0.0001,0.0003,0.0005,0.001 with accuracy A = 3; (b) dimensionless temperature of R20Z22 case versus z̃ at t̃=0.0001,0.0003,0.0005,0.001 with accuracy A = 4; and (c) dimensionless temperature of R20Z22 case versus z̃ at t̃=0.0001,0.0003,0.0005,0.001 with accuracy A = 5

Grahic Jump Location
Fig. 6

(a) R00Z22 geometry model with infinite radial domain (0<r<∞) and (b) R02Z22 geometry model with finite radial domain (0<r<R2)



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