An analytical solution of the hyperbolic heat conduction equation for the case of a finite medium symmetrically heated on both sides
M. Lewandowskaa, L. Malinowskib*
aInstitute of Physics, Szczecin University of Technology, 70-311 Szczecin, Al. Piastów 48, Poland
bFaculty of Maritime Technology, Szczecin University of Technology, 71-065 Szczecin, Al. Piastów 41, Poland
*Corresponding author E-mail address: (L. Malinowski); tel.: +48 91 4494827; fax: +48 91 4494488
Abstract
This paper presents an analytical solution of the hyperbolic heat conduction equation for the case of a thin slab symmetrically heated on both sides. In the mathematical model adopted, the heating is treated as an internal heat source with capacity dependent on coordinate and time, while walls of the slab are assumed to be insulated. The solution is obtained by Laplace transforms method taking advantage of the method of superposition. The analytical solution is validated by comparison with the results from a numerical model.
Nomenclature
athermal diffusivity
cpspecific heat at constant pressure
gcapacity of internal heat source
Ilaser incident intensity
Irarbitrary reference laser intensity
I0modified Bessel function, 0th order
kthermal conductivity
lthickness of the slab
Ldimensionless thickness of the slab
L-1inverse Laplace operator
Rsurface reflectance
sLaplace variable
ttime
tiduration of laser pulse
tkrelaxation time of heat flux
Ttemperature
Tm, T0arbitrary reference temperatures
unity step function
wspeed of heat propagation
xCartesian coordinate
Xdimensionless Cartesian coordinate
dimensionless absorption coefficient
Dirac delta function
dimensionless rate of energy absorbed in the medium
absorption coefficient
dimensionless temperature
density
dimensionless time
dimensionless duration of laser pulse
dimensionless capacity of internal heat source
constant coefficients related to the dimensionless capacity of internal heat
source
frequency of a periodic heat source
Superscript
transformed variable
Keywords: Hyperbolic equation of heat conduction; Finite medium; Analytical solution; Laplace transforms method; Superposition method
1. Introduction
In highly unsteady situations, the parabolic heat conduction equation based on the Fourier law fails, so the need for more adequate model of heat conduction, which permits the finite speed of heat flux, has arisen. There have been numerous attempts to formulate a new model in the literature [1-4], but it seems that, at present, the most frequently used is the hyperbolic model of heat conduction introduced by Cattaneo [5]. This model owns its popularity to simplicity and effectiveness. Various cases of hyperbolic heat conduction in a finite medium were studied analytically [6-12] and numerically [10, 13-15]. Recently, Torii et al. [16] solved numerically the case of a thin film subjected to a symmetrical heating on both sides. In this paper, we solve the same problem analytically by the method of Laplace transforms.
2. Model
We consider a thin slab of thickness L, initially at temperature , with constant thermophysical properties and insulated walls. At time , laser heat generation starts at both walls of the slab, giving rise to two thermal waves travelling in opposite directions. The temperature field in the slab can be described by the following hyperbolic equation of heat conduction
(1)
where tk is the relaxation time which represents a delay of the heat flow after a temperature gradient has been imposed. The relaxation time is related to the speed of propagation of thermal wave in the medium, w,by
(2)
The heat source term in Eq. (1) is modelled as
(3a)
where
(3b)
(3c)
and are the capacities of the internal heat sources acting at the left-hand side wall and at the right-hand side wall of the slab, respectively. Eq. (3b), used by Blackwell [17] and Zubair et al. [18], describes internal absorption of laser radiation. For convenience of subsequent analysis, we introduce the following dimensionless quantities
(4a)
(4b)
(4c)
(4d)
Eq. (1) is expressed in terms of the dimensionless variables (4a)-(4d) as
(5)
The dimensionless forms of Eqs. (3a)-(3c) are
(6a)
(6b)
(6c)
where
(6d)
(6e)
(6f)
The dimensionless initial conditions for the present problem are
(7a)
(7b)
The dimensionless boundary conditions are
(8a)
(8b)
Eq. (7b) is derived from the energy conservation equation on the assumption that there is no heat flow in the body at the initial moment.
3. Analytical solution
The boundary value problem of Eqs. (5), (6a) - (6f), (7a), (7b), (8a), and (8b) is solved by the method of Laplace transforms. At first, we solve Eq. (5) for
(9a)
to obtain the solution . Next, we solve Eq. (5) for
(9b)
to obtain the solution . Finally, as the problem is linear, we superimpose the two solutions
(10)
We substitute Eq. (6b) for in Eqs. (5) and (7b) to obtain, respectively
(11)
(12)
Taking the Laplace transform of Eq. (11) and using the initial conditions given by Eqs. (7a) and (12), yields
(13)
Transforming the boundary conditions given by Eqs. (8a) and (8b) gives
(14a)
(14b)
The solution of Eq. (13) satisfying boundary conditions (14a) and (14b) is the function
(15a)
where
(15b)
(15c)
(15d)
(15e)
To invert Eq. (15a) and find the time solution, we expand the terms and in binominal series
(16a)
(16b)
Moreover, we use the following pair of transforms
(17)
The inverse Laplace transform of solution (15a) is
(18a)
where
(18b)
(18c)
(18d)
(18e)
(18f)
4. Solutions for special cases of heat source capacity
The RHS of Eq. (18c) can be considerably simplified for some particular forms of [20]. Below, there are presented expressions for for the cases examined in this paper.
4.1. Instantaneous source,
(19)
4.2. Source of time independent strength,
(20)
4.3. Rectangular pulse source,
(21)
4.4. Periodic source,
(22a)
where
(22b)
(22c)
(22d)
(22e)
5. Validation of the solution
The analytical solution given by Eqs. (10), (18a) - (22e) was validated by a numerical solution of the problem given by Eqs. (5), (6a) – (6c), (7a) - (8b). The MacCormack algorithm [19] was used. Agreement between the analytical and numerical solutions was very good. We additionally confirmed the correctness of our results by using the energy balance equation for the whole slab
(23)
We also compared our results with those reported by Torii et al. [16] and found some discrepancies between them. Moreover, we found some errors in the paper by Torii et al. [16]. It seems that in the caption to Fig. 2, the equation should read . In the caption to Fig. 3, the equation should read . In the captions to Figs. 4 - 9, it should be , not . In Figs. 4a and 4b, should be equal to 0.01, not 0.1. In the caption to Fig. 6, it should be , , not , .
6. Sample calculations and discussion
Using our analytical solution, we performed sample calculations of temperature profiles in the slab for all considered types of heat source and chosen values of: , and other parameters concerning a particular type of heat source. The results of calculation are presented in Figs. 1(a) - 3(b). The heat waves excited at the vicinity of both side surfaces of the slab travel in the opposite directions, superimpose, and reflect back and forth between the surfaces. It is worth noting that the heat wave travels the distance of x in time of . Taking into account Eqs. (2), (4a), and (4b) we obtain . So, the heat wave covers the distance of L in time equal to L.
Figs. 1(a) and 1(b) displays the time-dependent temperature distribution in the slab for the instantaneous heat sources for which (par. 4.1). It is difficult to model numerically, with good accuracy, such a source, but in the analytical solution, the function takes quite a simple form. Torii et al. did not consider the instantaneous heat sources in their work [16]. It is seen in. Fig. 1(a) that the fronts of both waves meet in the middle of the slab for . They meet for all times calculated from the equation for .
Temperature profiles for the sources of time independent strength, (par. 4.2), are presented in Figs. 2(a) - (2c). For and (Fig. 2(a)), the heat produced by both sources is nearly evenly distributed through the slab which explains the flat temperature profiles for particular times. For large values of and/or L, the strength of heat generation varies a lot across the slab thickness that is manifested by substantial variation of temperature profiles (Figs. 2(b) and 2(c)). There are significant quantitative differences, growing with time, between the majority of results presented in Fig. 2(a) - 2(c) and the results for the same data sets reported by Torii et al [16]. The temperature profiles calculated from our model lie higher. For and , our results overlap with those presented by Torii et al.
Figs. 3(a) and 3(b) shows the temperature profiles in the slab for the rectangular pulse heat sources, , (par. 4.3). In the limiting cases of very small and very large values of , the time characteristic of the considered source becomes close to the time characteristic of the instantaneous source, or of the source of time independent strength, respectively. In Fig. 3(a) there is seen an overshooting of temperature for and , specific for the hyperbolic model of heat conduction.
For the periodic heat sources, , (par. 4.4), our results are the same as those reported by Tori et al [16]. For the case of pulsed heat sources considered by Torii et al. [16], we obtained higher values of temperatures; only for duration of the pulse equal to 5 the results are similar. In general, we observe better agreement between our and the Torii et al. results for large values of L and .
6. Conclusions
The problem of thermal response of a thin slab to a symmetric rapid heating in the vicinity of its both side surfaces is solved analytically. The results are validated by numerical calculations. The agreement is very good. The results are also compared with those obtained from a numerical model by Torii at al. [16]. In great part (especially for smaller values of L and ) our results are higher than those reported by Torii at al., probably because of inaccuracy in the Torii at al. calculations.
References
[1]D.D. Joseph, L. Preziosi, Heat waves, Rev. Mod. Phys. 61 (1) (1989) 41-73.
[2]D.D. Joseph, L. Preziosi, Addendum to the paper "Heat waves", Rev. Mod. Phys. 62 (2) (1990) 375-391.
[3]M.N. Ozisik, D.Y. Tzou, On the wave theory in heat conduction, J. Heat Transfer, 116 (1994) 526-535.
[4]D.Y. Tzou, Macro- to Microscale Heat Transfer, Taylor & Francis, Washington, 1997.
[5]C. Cattaneo, Sulla conduzione de calore, Atti. Semin. Mat. Fis. Univ. Modena 3 (1948) 3-21.
[6]S.H. Chan, M.J.D. Low, W.K. Mueller, Hyperbolic heat conduction in catalytic supported crystallites, AIChE Journal 17 (6) (1971) 1499-1501.
[7]Y. Taitel, On the parabolic, hyperbolic and discrete formulation of the heat conduction equation, Int. J. Heat Mass Transfer 15 (1972) 369-371.
[8]M.N. Ozisik, B. Vick, Propagation and reflection of thermal waves in a finite medium, Int. J. Heat Mass Transfer 27 (10) (1984) 1845-1854.
[9]J. Gembarovic, V. Majernik, Non-Fourier propagation of heat pulses in finite medium, Int. J. Heat Mass Transfer 31 (5) (1988) 1073-1080.
[10]A. Kar, C.L. Chan, J. Mazumder, Comparative studies on nonlinear hyperbolic and parabolic heat conduction for various boundary conditions: analytic and numerical solutions, J. Heat Transfer 114 (1992) 14-20.
[11]L. G. Hector, W-S Kim, M.N. Ozisik, Propagation and reflection of thermal waves in a finite medium due to axisymmetric surface sources, Int. J. Heat Mass Transfer 35 (4) (1992) 897-912.
[12]D.W. Tang, N. Araki, Analytical solution of non-Fourier temperature response in a finite medium under laser-pulse heating, Heat Mass Transfer 31 (1996) 359-363.
[13]D.C. Wiggert, Analysis of early-time transient heat conduction by method of characteristics, J. Heat Transfer 99 (1977) 35-40.
[14]D.E. Glass, M.N. Ozisik, D.S. McRae, B. Vick, On the numerical solution of hyperbolic heat conduction, Numer. Heat Transfer 8 (1985) 497-504.
[15]H-T Chen, J-Y Lin, Numerical analysis for hyperbolic heat conduction, Int. J. Heat Mass Transfer 36 (11) (1993) 2891-2898.
[16]S. Torii, W-J Yang, Heat transfer mechanisms in thin film with laser heat source, Int. J. Heat Mass Transfer 48 (2005) 537-544.
[17]B. F. Blackwell, Temperature profile in semi-infinite body with exponential source and convective boundary conditions, J. Heat Transfer 112 (1990) 567-571.
[18]S.M. Zubair, M.A. Chaudhry, Heat conduction in a semi-infinite solid due to time-dependent laser source, Int. J. Heat Mass Transfer 39 (14) (1996) 3067-3074.
[19]D.A. Anderson, J.C. Tannehill, R.H. Pletcher, Computational Fluid Mechanics and Heat Transfer, New York, Hemisphere, 1984.
[20]M. Lewandowska, Hyperbolic heat conduction in the semi-infinite body with a time-dependent laser heat source, Heat Mass Transfer 37 (2001) 333-342.
Figure captions
Fig. 1. Temperature profiles in the slab for the instantaneous heat sources, . , . (a) , (b) .
Fig. 2. Temperature profiles in the slab for the heat sources of time independent strength, . . (a) , , (b) , , (c) , .
Fig. 3. Temperature profiles in the slab for the rectangular pulse heat sources, . , . (a) , , (b) , .
1