Research Papers

Estimating Over Temperature and Its Duration in a Flat Plate With Sudden Changes in Heating and Cooling

[+] Author and Article Information
Chien-Shing Lee

School of Aeronautics and Astronautics,
Purdue University,
West Lafayette, IN 47907-2045
e-mail: cslee@purdue.edu

Tom I-P. Shih

School of Aeronautics and Astronautics,
Purdue University,
West Lafayette, IN 47907-2045
e-mail: tomshih@purdue.edu

Kenneth Mark Bryden

Department of Mechanical Engineering,
Iowa State University,
Ames, IA 50011
e-mail: kmbryden@iastate.edu

1Corresponding author.

Contributed by the International Gas Turbine Institute (IGTI) of ASME for publication in the JOURNAL OF TURBOMACHINERY. Manuscript received December 3, 2015; final manuscript received December 9, 2015; published online February 9, 2016. Editor: Kenneth C. Hall.

J. Turbomach 138(6), 061007 (Feb 09, 2016) (12 pages) Paper No: TURBO-15-1291; doi: 10.1115/1.4032306 History: Received December 03, 2015; Revised December 09, 2015

When the operating condition of a gas turbine engine changes from one steady-state to another, the cooling must ensure that the solid's temperatures never exceed the maximum allowable throughout the transient process. Exceeding the maximum allowable temperature is possible even though cooling is increased to compensate for the increase in heating because there is a time lag in how the solid responds to changes in its convective heating and cooling environments. In this paper, a closed-form solution (referred to as the 1D model) is derived to estimate the over temperature and its duration in a flat plate subjected to sudden changes in heating and cooling rates. For a given change in heating rate, the 1D model can also be used to estimate the minimum cooling needed to ensure that the new steady-state temperature will not exceed the maximum allowable. In addition, this model can estimate the temperature the material must be cooled to before imposing a sudden increase in heat load to ensure no over temperature throughout the transient process. Comparisons with the exact solutions show the 1D model to be accurate within 0.1%. This 1D model was generalized for application to problems in multidimensions. The generalized model was used to estimate the duration of over temperature in a two-dimensional problem involving variable heat transfer coefficient (HTC) on the cooled side of a flat plate and provided results that match the exact solution within 5%.

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


Lee, C.-S. , Shih, T. I.-P. , and Bryden, K. M. , 2013, “ CFD Conjugate Analysis of a Jet-Impingement Configuration With Sudden Changes in Heating and Cooling Loads,” ASME Paper No. IMECE2013-65542.
Perlmutter, M. , and Siegel, R. , 1961, “ Unsteady Laminar Flow in a Duct With Unsteady Heat Addition,” ASME J. Heat Transfer, 83(4), pp. 432–440. [CrossRef]
Sucec, J. , 1987, “ Exact Solution for Unsteady Conjugated Heat Transfer in the Thermal Entrance Region of a Duct,” ASME J. Heat Transfer, 109(2), pp. 295–299. [CrossRef]
Sucec, J. , and Weng, H.-Z. , 2002, “ Unsteady, Conjugated Duct Heat Transfer Solution With Wall Heating,” AIAA J. Thermophys. Heat Transfer, 16(1), pp. 128–134. [CrossRef]
Takahashi, T. , Watanabe, K. , and Takahashi, T. , 2001, “ Transient Analyses of Conjugate Heat Transfer of a First Stage Rotor Blade in Start-Up and Shut-Down,” ASME Paper No. 2001-GT-0171.
Rossette, A. H. , C, Z. M. , Demeulenaere, A. , and Roque Lopez Hernandez, J. A. , 2009, “ The Effect of Start-up Cycle in Ceramic Coating Used as Thermal Barrier for a Gas Turbine Bucket,” Appl. Therm. Eng., 29(14–15), pp. 3056–3065. [CrossRef]
Sun, Z. , Chew, J. W. , Hills, N. J. , Volkov, K. N. , and Barnes, C. J. , 2010, “ Efficient Finite Element Analysis/Computational Fluid Dynamics Thermal Coupling for Engineering Applications,” ASME J. Turbomach., 132(3), p. 031016. [CrossRef]
Sun, Z. , Chew, J. W. , Hills, N. J. , Lewis, L. , and Mabilat, C. , 2011, “ Coupled Aerothermomechanical Simulation for a Turbine Disk Through a Full Transient Cycle,” ASME J. Turbomach., 134(1), p. 011014. [CrossRef]
Carslaw, H. S. , and Jaeger, J. C. , 1959, Conduction of Heat in Solids, Clarendon Press, Oxford, UK.
Arpaci, V. S. , 1966, Conduction Heat Transfer, Addison-Wesley, Reading, MA.
Ozisik, M. N. , 1993, Heat Conduction, 2nd ed., Wiley, Hoboken, NJ.
Shih, T. I.-P. , 2010, “ Airbreathing Engines in Propulsion and Power,” Encyclopedia of Aerospace Engineering, Vol. 2, R. Blockley and W. Shyy , eds., Wiley, Chichester, UK, pp. 889–1074.
Shih, T. I.-P. , and Yang, V. , eds., 2014, Turbine Aerodynamics, Heat Transfer, Materials, and Mechanics (Progress in Astronautics and Aeronautics, Vol. 243), American Institute of Aeronautics and Astronautics, Reston, VA, pp. 311–421.
Zieliñska, M. , Yavorska, M. , Porêba, M. , and Sieniawski, J. , 2010, “ Thermal Properties of Cast Nickel Based Superalloys,” Arch. Mater. Sci. Eng., 44(1), pp. 35–38.
Lax, P. D. , and Richtmyer, R. D. , 1956, “ Survey of the Stability of Linear Finite Difference Equations,” Commun. Pure Appl. Math., 9(2), pp. 267–293. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic of the problem used to develop the 1D model

Grahic Jump Location
Fig. 2

Temperature profiles in the plate during the first time domain

Grahic Jump Location
Fig. 3

Temperature profiles in the plate during the second and third time domains

Grahic Jump Location
Fig. 4

Duration of transient process

Grahic Jump Location
Fig. 5

Duration of over temperature predicted by the exact solution and by the 1D model

Grahic Jump Location
Fig. 6

Evolution of the temperature in the plate as a function of time for ∅ = 0, 0.25, 0.5, and0.75

Grahic Jump Location
Fig. 7

Evolution of the maximum temperature in the plate as a function of time for ∅ = 0, 0.25, 0.5, and 0.75

Grahic Jump Location
Fig. 8

Schematic of the 2D problem

Grahic Jump Location
Fig. 9

Dimensionless surface temperature on the heated side of the plate as a function of the Biot numbers on the cooled and heated sides under steady-state conditions

Grahic Jump Location
Fig. 10

Dimensionless temperature at x = 0 as a function of Fourier number and dimensionless initial temperature

Grahic Jump Location
Fig. 11

Dimensionless temperature at x=0 as a function of the Fourier number and the Biot number, Bi1

Grahic Jump Location
Fig. 12

Dimensionless temperature at x=0 as a function of the Fourier number and the Biot number, Bi2

Grahic Jump Location
Fig. 13

Time constant as a function of thermal resistance and capacitance for time domains II and III




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