Abstract
An analytic solution method, without integral transformation, is developed to find the exact solutions for transient heat conduction in functionally graded (FG) circular hollow cylinders with time-dependent boundary conditions. By introducing suitable shifting functions, the governing second-order regular singular differential equation with variable coefficients and time-dependent boundary conditions is transformed into a differential equation with homogenous boundary conditions. The exact solution of the system with thermal conductivity and specific heat in power functions with different orders is developed. Finally, limiting studies and numerical analyses are given to illustrate the efficiency and the accuracy of the analysis.
1. Introduction
The applications of heat conduction in functionally graded (FG) circular hollow cylinders with time-dependent boundary conditions can be widely found in many engineering fields, such as cannon barrels, heat exchanger tubes, time varied heating on walls of circular structure, and heat treatment on hollow cylinders. As such, an accurate solution method is very helpful for relevant developments.
The problem of heat conduction with time-dependent boundary conditions cannot be solved directly by the method of separation of variables. In most of the analyses, an integral transform has been used to remove the time-dependent term. For the problem of heat conduction in circular hollow uniform cylinders with time-dependent boundary conditions, the associated governing differential equation is a second-order Bessel differential equation with constant coefficients. After conducting a Laplace transformation, the analytical solution can be obtained and found in the book by Őzisik [1].
When the structure is an FG circular hollow cylinder, the associated governing differential equation is a second-order regular singular differential equation with variable coefficients. For problems with time-independent boundary conditions, various numerical methods have been developed, such as the perturbation method [2], the finite difference method [3], and the finite element method [4]. Jabbari et al. [5, 6] derived analytical solutions for thermal stresses of FG hollow cylinders whose material properties vary with power law distribution through the thickness due to radially symmetric loads and nonaxisymmetric loads. By using the Laplace transformation and a series expansion of Bessel functions, Ootao and Tanigawa [7] analyzed one-dimensional transient thermoelastic problem with the material properties varying with the powerlaw form of the radial coordinate variable. Zhao et al. [8] analyzed the temperature change when the thermal and thermoelastic properties are assumed to vary exponentially in the radial direction. Hosseini et al. [9] considered the material properties to be nonlinear with a power law distribution through the thickness, while the temperature distribution was derived analytically using the Bessel functions.
In the study of heat conduction in FG circular hollow cylinders with time-dependent boundary conditions, only limited studies can be found. Ootao et al. [10] studied the transient temperature and thermal stress distribution in an infinitely long nonhomogeneous hollow cylinder due to a moving heat source in the axial direction from the inner and/or outer surfaces using the layerwise theory in conjunction with the method of Fourier cosine and the Laplace transformations. Shao and Ma [11] employed the Laplace transform techniques and the series solving method to study thermomechanical stresses in FG circular hollow cylinders with linearly increasing boundary temperature. Jabbari et al. [12] developed an analytical solution for the one-dimensional temperature distribution, mechanical and thermal stresses in an infinitely long FG hollow cylinder under a moving heat source, which moves across the thickness of the cylinder. Asgari and Akhlaghi [13] employed the finite element method to study the transient thermal stresses in a thick hollow cylinder with finite length mad of 2D-FGM that its material properties are varied in the radial and axial directions with a power law function. The thermal boundary conditions at the inner and outer radiuses are time dependent. Singh et al. [14] applied the finite integral transform method and the separation of variables method to solve time-dependent heat conduction problem in a multilayer annulus. Malekzadeh and Heydarpour [15] used the differential quadrature method (DQM) in conjunction with the finite element method (FEM) to study the response of FG cylindrical shells under moving thermomechanical loads. Recently, Wang and Liu [16] have employed the method of separation of variables to develop the analytical solution of transient temperature fields for two-dimensional transient heat conduction in a fiber-reinforced multilayer cylindrical composite.
In the study of thermal elastic response of FG cylinders without mechanical loading, the heat conduction problem is not incorporated with the elastic field and can be studied independently. However, the thermal field will be coupled with the temperature field. In this paper, one considers the heat conduction problem of FG cylinders only. A new solution method, which is a modification on the method developed by Lee and Lin [17] and Chen et al. [18], is used to develop the analytical solution for transient heat conduction in FG circular hollow cylinders with time-dependent boundary conditions. By introducing suitable shifting functions, the governing second-order differential equation with variable coefficients and time-dependent boundary conditions is transformed into a differential equation with homogenous boundary conditions. The analytic solution of the system with thermal conductivity and specific heat in power functions with different orders is developed. Finally, limiting studies and numerical analysis are given to illustrate the efficiency and accuracy of the solution method.
2. Mathematical Modeling
Consider the transient heat conduction in an FG circular hollow cylinder with time-dependent boundary condition at the inner and outer surfaces, as shown in Figure 1. The governing differential equation of the system is The boundary conditions are and the initial condition is Here, is the space variable and is the time variable, is the specific heat, is the thermal conductivity, is the temperature, is the mass density, and is the heat source inside the circular hollow cylinder. and are the time-dependent temperatures at the inner and outer surfaces, respectively.
In terms of the following dimensionless parameters where is the reference temperature, the boundary value problem of heat conduction becomes
3. Solution Method
3.1. Change of Variables
To find the solution for the second-order differential equation with nonhomogeneous boundary conditions, the shifting variable method developed by Lee and Lin [17] and Chen et al. [18] was extended by taking where , are shifting functions to be specified and is the transformed function. Substituting (6) into (5) yields the following partial differential equation: the associated boundary conditions: and the associated initial condition
3.2. Shifting Functions
To simplify the analysis, the shifting functions are specifically chosen such that they satisfy the following differential equations and the boundary conditions These two shifting functions can be easily determined as
3.3. Reduced Homogenous Problem
With these two shifting functions, (11), the governing differential equation (7) becomes where The two nonhomogenous boundary conditions, (8), for the transformed function are reduced to homogenous ones: The transformed initial condition is
3.4. Solution of Transformed Variable
3.4.1. Characteristic Solution
To find the solution , we use the eigenfunction expansion method and assume the solution to be in the form The separation equation for the dimensionless time variable is and the dimensionless space variable satisfies the following self-adjoin operator: Let , , be the two linearly independent fundamental solutions of the system; then, the solution of the differential equation (18) can be expressed as where and are constants to be determined from the boundary conditions, (19)-(20).
After substituting solutions, (21), into the boundary conditions, (19)-(20), we obtain the following characteristic equation Consequently, the eigenvalues , can be determined. The associated th eigenfunction is determined as where and are, respectively, the fundamental solutions and associated with eigenvalues , . They are defined as and . The eigenfunctions constitute an orthogonal set in the interval , with respect to a weighting function : where
In terms of eigenfunctions, the transformed variable can be expressed as Substituting (26) into (12), multiplying it by , and integrating from to , we obtain where The general solution of (27) is is determined from the initial condition (15), and After substituting the solution of the transformed variable (26) and the shifting functions (11) back to (6), the exact solution for the general system is obtained.
4. Verification and Examples
To illustrate the previous analysis, the following examples and limiting cases are given.
Example 1. Consider the heat conduction in an uniform circular hollow cylinder with time-dependent boundary conditions. The boundary value problem of the heat conduction in dimension-less form is
In this case, , . The two shifting functions are
The two linearly independent fundamental solutions are
This leads to
where the characteristic equation is
The associated th eigenfunction is determined as
The eigenvalues and the associated eigenfunctions are obtained from (35) and (36). The two coefficients in (28)-(29) are derived as
where
As a result, the analytic closed solution of the system in dimensionless form is
when
The analytic closed solution, in dimensionless form, is reduced to
The solution is exactly the same as the one given by Őzisik [1].
Example 2. Consider the heat conduction in an FG circular hollow cylinder with time-dependent boundary conditions. The coefficients of thermal conductivity and the specific heat are and , respectively. The boundary value problem of the heat conduction in dimensionless form is
In this case, two shifting functions are
The route to two independent fundamental solutions of of (42) lies in the use of the Frobenius method which can be represented in terms of the Bessel functions:
where
Now,
where
The characteristic equation is
The associated th eigenfunction is determined as
The eigenvalues and the associated eigenfunctions are obtained from (51) and (52) by numerical analysis. Two coefficients in (28)–(30) are derived as
where
Consequently, the analytic closed solution for the system can be derived as
In Figure 2, the temperature variation of FG circular hollow cylinders with various parameters of , , and at is shown. It can be found that when is a positive constant, the temperature parameter of the mediums at increases then reaches the associated constant temperatures over time. The temperature increase rate for the system with a higher value of is greater than that of one with a lower value of . Figures 3 and 4 show the temperature distribution of FG circular hollow cylinders with various parameters of , , and at . From these figures, it can be observed that with constant parameter , the temperature of the mediums increases as parameter is increased. With constant parameter , the temperature of the mediums increases as parameter is increased.
Example 3. Consider the same physical system as discussed in Example 2. In this case, the time dependent boundary condition at , (44), is changed to the form
In this case, the eigenvalues and eigenfunctions are the same as those given in Example 2.
Now,
Following the same solution procedures as shown, the exact solution for the general system can be derived as
where
Figure 5 shows the harmonic temperature variation of FG circular hollow cylinders with constant value of , various parameters of at , , and . It can be observed that with constant value of , the amplitude of temperature oscillation for the system with a higher value of will be more than that of one with a lower value of .
In Table 1, the temperature variations of FG circular hollow cylinders with constant value of and various values of and at are given. It can be observed that with constant parameter , the temperature of the mediums will decrease as parameter is decreased.
5. Conclusions
The problem of heat conduction with general time-dependent boundary conditions cannot be solved directly by the method of separation of variables. In most of the analyses, an integral transform was used to remove the time-dependent term. In this paper, a new analytic solution method is developed to find the analytic closed solutions for the transient heat conduction in FG circular hollow cylinders with general time-dependent boundary conditions. The developed solution method is free of any kind of integral transformation.
By introducing suitable shifting functions, the governing second-order regular singular differential equation with variable coefficients and time-dependent boundary conditions is transformed into a differential equation with homogenous boundary conditions. The analytic solution of the system with thermal conductivity and specific heat in power functions with different orders is developed. Finally, limiting studies and numerical analyses are given to illustrate the efficiency and the accuracy of the analysis. The proposed solution method can also be extended to the problems with various kinds of FG materials and time-dependent boundary conditions.
Acknowledgment
It is gratefully acknowledged that this research was supported by the National Science Council of Taiwan, Taiwan, under Grant NSC 99-2221-E-006-021.