Semianalytical Solution of Transient Heat Transfer for Laminated Structures under Time-Varying Boundary Conditions

Li, Xi-Xia; Wang, Li; Li, Chang-Yu; Wan, Yi-Yu; Qian, Yu-Chen

doi:https://doi.org/10.1155/2021/5545566

Mathematical Problems in Engineering

On this page

Abstract Introduction Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5545566 | https://doi.org/10.1155/2021/5545566

Semianalytical Solution of Transient Heat Transfer for Laminated Structures under Time-Varying Boundary Conditions

Xi-Xia Li,¹Li Wang,¹Chang-Yu Li,¹Yi-Yu Wan,¹and Yu-Chen Qian¹

Academic Editor: Luis M. López-Ochoa

Received23 Feb 2021

Revised15 May 2021

Accepted15 Sept 2021

Published06 Oct 2021

Abstract

In order to solve the transient heat transfer problem of the laminated structure, a semianalytical method based on calculus is adopted. First, the time domain is divided into tiny time segments; the analytical solution of transient heat transfer of laminated structures in the segments is derived by using the method of separation of variables. Then, the semianalytical solution of transient heat transfer in the whole time domain is obtained by circulation. The transient heat transfer of the three-layer structure is analyzed by the semianalytical solution. Three time-varying boundary conditions (a: square wave, b: triangular wave, and c: sinusoidal wave) are applied to the surface of the laminated structure. The influence of some key parameters on the temperature field of the laminated structure is analyzed. It is found that the surface temperature of the laminated structure increases fastest when heated by square wave, and the maximum temperature can reach at 377°C, the temperature rises the most slowly when heated by the triangular wave, and the maximum temperature is 347°C. The novelty of this work is that the analytical method is used to analyze the nonlinear heat transfer problem, which is different from the general numerical method, and this method can be applied to solve the heat transfer problem of general laminated structures.

1. Introduction

In practical engineering applications, there are many problems of transient heat transfer in laminated structures [1, 2]. For example, transient heat transfer of human skin during hyperthermia, transient heat transfer of multilayer insulation materials, transient heat transfer of thermal barrier coating, or new thermoelectric materials [3–7]. More and more attention has been paid to transient heat transfer in laminated structures. The relevant research studies can be summarized as follows.

Xu et al. established a three-layer heat transfer model aiming at the skin tissue heat transfer problem by means of the finite difference method [8]. Liu et al. solved the temperature field of the three-layer heat transfer structure using the inverse Fourier transform method [9]. Campo and Arıcı exploited a radically different approach applying the Method of Lines (MOL) forthwith to the one-dimensional, unsteady, and heat conduction equation with prescribed uniform heat flux [10]. Li and Lin used the method of separation of variables to obtain the semianalytical solution of heat transfer in a laminated structure [11]. Yu et al. studied the heat transfer of multilayer structures by using the smooth particle manifold (SPH) method [12]. Salehi et al. investigated flow and heat transfer of magnetohydrodynamic squeezing nanofluid flow between two infinity parallel plates by Akbari-Ganji’s method [13]. Fu and Cheng solved the semianalytical solution of temperature field analysis of the three-dimensional structure of the laminated cylinder by using the fine integration method [14]. Ramadan solved the transient heat transfer problem of the multilayer structure by the semianalytical method and analyzed the non-Fourier characteristics of the heat transfer process [15]. Li et al. studied transient heat transfer in laminated structures by using a backward differential method [16]. Liu and Shi used the lattice Boltzmann method to solve the numerical solution of heat transfer in multilayer structures [17]. Shan et al. studied the heat transfer of a multilayer structure using a modified scale-boundary finite element method [18]. Guo et al. simulated the transient heat transfer problem of the three-layer structure by using the improved time-domain discontinuous Galerkin method and improved the numerical fluctuation problem of the traditional method [19].

From the literature analysis, it can be seen that, for the heat transfer of laminated structures, numerical methods are mostly used, and the boundary conditions generally do not change with time. The transient heat transfer of multilayer materials is a nonlinear problem because the thermophysical parameters of each layer are different; in addition, the boundary conditions change with time, and it is difficult to obtain the analytical solution of the problem directly [20–23]. In this paper, the boundary conditions are differentiated in the time domain, and a special method of separating variables is used to obtain the analytical solution of the temperature field of the laminated structure in a small time period, and then, the cycle is carried out to obtain the analytical solution of the temperature field in the whole time domain. Then, the temperature field of the three-layer structure heated by time-varying heat flow is analyzed by using the method presented in this paper. The temperature field of the structure under three heating conditions (square wave, triangular wave, and sine wave) has been obtained, and the influence of surface convective heat transfer coefficient, heating frequency, and other parameters on the temperature field distribution was studied.

2. Mathematical Model

As shown in Figure 1, the number of layers of the laminated structure is n, and the heat flow of f₁(t) is heated on the left side of the laminated structure. At the same time of heating, the surface of the laminated structure is cooled by fluid. The temperature on the right side of the laminated structure is f_n(t). The thermophysical parameters of each layer of the laminated structure are consistent, but they vary from layer to layer.

A mathematical and physical model can be established based on the problem studied, and the corresponding governing equation can be expressed as [24]where is the heat generation and is the interacts with the environment by convective heat transfer from upper and lower surfaces. In view of some special conditions, there is a heat transfer phenomenon in the laminated material through a certain medium. Here, h_i is defined as the internal heat transfer coefficient of ith layer.

The boundary conditions are as follows.

When x = 0,

When ,

The initial condition is

After dimensionless, the governing equation, boundary condition, and initial conditions can be expressed as follows:

The boundary condition of dimensionless is as follows.

When ,

When , in other words, at the left end of the laminated structure, the temperature changes according to , where :

The initial conditions for dimensionless is

3. Solution Method

3.1. Differential in Time Domain

In general, for the heat transfer problem of the laminated structure under time-varying boundary conditions, the analytical solution cannot be obtained directly by the separation of variables method because the thermophysical parameters of each layer of the structure are not consistent and the boundary conditions change with time. In order to solve the transient heat transfer problem of the laminated structure with time-varying boundary conditions, it is necessary to carry out differentiation in the time domain. Suppose there is a sufficiently small time interval , according to the idea of calculus, in the time domain of [0 ], the following hypothesis is true.

When ,

When (i = 1,2,3, …, n),where is the differential approximation function of the left boundary condition of the laminated structure in [0 ]: is the function of temperature change at the junction of layer i and layer i + 1 of the laminated structure and is the slope of the function , which is an unknown constant. At this point, boundary conditions on both sides of each layer of the inner layered structure at time [0 ] are known. Therefore, the method of separation of variables can be used to solve the nonlinear problem.

3.2. The First Analytic Solution in Time [0 ]

For the first layer in the first [0 ] time domain, the governing equation is as follows:

The corresponding boundary conditions are as follows.

When ,

The initial condition is

Through analysis, it can be known that the boundary conditions of governing equation are inhomogeneous. In order to homogeneous the boundary conditions, equations (16) and (17) are assumed [25]:

Under homogeneous boundary conditions, needs to satisfy equations (18) and (19):

By substituting equation (17) into equations (18) and (19), then a, b, c, and d in equation (17) can be derived:

By substituting equation (16) into equation (12), equation (21) can be obtained:

The characteristic equation of equation (21) is

In order to separate variables, let

By substituting (23) into (22), equations (24) and (25) can be obtained:

It can be seen from equations (13), (14), (18), and (19) that

Substituting equation (23) into equations (26) and (27), we can obtain

By combining equations (24), (28), and (29), which is a second-order homogeneous ordinary differential equation, its characteristic function iswhere is the positive solution of

Equation (32) is established by calculation:

Equation (32) expresses that satisfies the orthogonality property. In order to separate variables, let

By substituting equation (33) into equation (21) and integrating within [0, ], equation (34) can be obtained:

Equation (34) is a first-order ordinary differential equation with respect to , and its solution is [26]where

Combined with the above derivation process, it can be seen that the solution of the first layer of the laminated structure in the first differential period can be expressed by

3.3. The ith Layer Analytic Solution in Time [0 ]

The difference between the mathematical and physical models of the second to ith layers and the first layer is that the left side of the first layer is the time-varying heat flow boundary condition and the right side is the time-varying temperature boundary condition. The left and right sides of layers 2 to i are time-varying temperature boundary conditions. The analytical solution of the ith (i ≥ 2) layer in the first differential period can be expressed as equation (40) by the same method as that of the first layer:

The parameters in equation (40) are shown as follows:where is the positive root of

3.4. The Solution of the Temperature Field in the Whole Time Domain

According to Sections 3.1–3.3, the analytic solutions of each layer structure in the first [0 ] time domain are derived. However, there is an unknown coefficient which included in the analytic solutions. According to the equal temperature at the junction and the equal heat flow out and in at the junction, equation (48) can be obtained:

If is substituted into equation (48), equation (49) can be obtained:

So,

can be solved by equation (49). At this time, the analytic solution of the first [0 ] time domain can be completely solved, and can be solved. With at this time as the initial value of the next differential period, the semianalytic solution of each layer of the layered structure in the entire time domain can be solved by executing the loop.

4. Case Study and Discussion

4.1. Description of the Problem

In order to verify the correctness of the semianalytical solution derived in this paper, the semianalytical method has been used to analyze the temperature field variation of a three-layer structure under time-varying boundary conditions. The dimensions of the three-layer structure and thermophysical parameters of each layer are shown in Table 1 [10]. Time-varying boundary conditions are applied to the left side of the laminated structure. The change of heating heat flow with time is shown in Figure 2. Three heating conditions are considered in this study (a: square wave heating, b: triangular wave, and c: sinusoidal wave). The period and amplitude of the three heating cases are the same. The former period of time is heating and cooling at the same time. In the later period, only cooling is carried out. When heating, the heat flow amplitude is 20 W/cm². When cooling, the fluid temperature is set as 0°C. The convective heat transfer coefficient of the fluid and the left boundary is set as h₀ = 10 W/m²K or h₀ = 300 W/m²K, the initial temperature of the structure is T_ina = 37°C, and the right boundary condition of the structure is f₃(t) = 37°C.

4.2. Analysis and Discussion

Figure 3 shows the change of surface temperature of the laminated structure with time. The data corresponding to line O are the result of the method of presented in this paper, and it is found that the results of this paper are completely consistent with those of Liu et al. [9]. Lines A, B, and C correspond to three different boundary conditions (a: square wave, b: triangular wave, and c: sinusoidal wave). When h₀ = 10 W/m²K, as can be seen from Figure 3 in the top 10 s, temperature rise by fluctuations on the surface of the laminated structure. When the square wave is heated, the amplitude of the wave is the largest, the sine wave is the second, and the triangle wave is the smallest; that is, it is caused by different effective heating areas. After 10 s heating, the surface temperature of the structure showed a downward trend due to the continued fluid cooling. Lines D, E, and F correspond to three different boundary conditions (a: square wave, b: triangular wave, and c: sinusoidal wave). When h₀ = 300 W/m²K, the surface temperature of the laminated structure also fluctuates. However, compared with h₀ = 10 W/m²K, it increases more slowly. The reason is that the convective heat transfer coefficient is larger and the cooling effect is better. Line G and line A are the results calculated by the finite element method and the semianalytical method in this paper under exactly the same conditions. It can be seen from Figure 3 that they are completely coincident, which further verifies the correctness of the method in this paper.

Figure 4 shows the temperature change over time at the junction of the first and second layers of the laminated structure, where lines A, B, and C, respectively, represent the temperature change with time under the conditions of square wave heating, triangular wave heating, and sine wave heating, when h = 10 W/m² K, and lines D, E, and F, respectively, represent the temperature change with time under the conditions of square wave heating, triangular wave heating, and sine wave heating when h = 300 W/m² K. Its variation trend is consistent with the surface trend of the laminated structure. The first 10 s heating process also has temperature fluctuation, but it is relatively gentle.

Figure 5 shows the temperature change over time at the junction of the second and third layers of the laminated structure. Here, the data of A, B, and C are calculated with h₀ = 10 W/m²K and the data of D, E, and F are calculated with h₀ = 300 W/m²K. It can be seen from Figure 5 that the temperature change is not fluctuating at this position. It is inferred that this location is far away from the heat source, and the fluctuation of heat flux of heat source cannot cause the fluctuation of temperature at this position. It can also be seen from Figure 5 that the larger the convective heat transfer coefficient, the slower the temperature rise. The temperature rise in this position is much lower than the temperature rise on the surface of the structure.

Figure 6 shows the temperature distribution along the length direction of the laminated structure at different times. The blue line is the result of the method presented in this paper, which is consistent with that of Liu et al. [9]. Thus, the correctness of the method presented can be verified. It can be seen from Figure 6 that, at 9.5 s, the temperature from the left to the right decreases successively because the fluctuating heating on the surface of the laminated structure has just finished. When h₀ = 300 W/m²K, the heat dissipation is better and the temperature rises rapidly, so the overall temperature of the structure is higher than that of h₀ = 10 W/m²K at the end of heating. It can also be seen that the temperature is highest when square wave heating is applied and lowest when triangle wave heating is applied. When t = 15 s, the heating ended and a period of cooling was carried out. It can be seen that, from the left to right side of the laminated structure, the temperature first increased and then decreased because the cooling effect was better near the surface of the structure on the left side.

Figure 7 shows the effect of heating period on the surface temperature change at the surface of the laminated structure. In Figures 7(a)–7(c), respectively, represent the temperature changes with period T_c = 0.5 s, T_c = 1 s, and T_c = 2 s. It can be seen from Figure 7 that the larger the heating period is, the larger the temperature fluctuation amplitude is. It can also be seen from Figure 7 that the temperature change is the most obvious when the square wave is heated, and the temperature change curve when the triangle wave is heated is relatively smooth. It is inferred that the larger the period is, the larger the limited heating area is, so the faster the temperature rises.

(a)

(b)

(c)

Figure 8 shows the effect of heating period on the surface temperature change at the junction of the first and second floors of the laminated structure. Figures 8(a)–8(c), respectively, represent the temperature change with period T_c = 0.5 s, T_c = 1 s, and T_c = 2 s. It can be seen from Figure 8 that the temperature change trend at the junction is consistent with the surface temperature change trend. However, because it is far away from the surface heating source, the range of temperature change is smaller than the surface. At the same time, it can be seen that the longer the period is, the smaller the frequency is and the larger the range of temperature fluctuation is.

(a)

(b)

(c)

4.3. Summary of Results

Through the above analysis and discussion, the relevant results can be summarized, as shown in Table 2.

5. Conclusion

In this paper, a general semianalytical solution for transient heat transfer of the laminated structure under time-varying boundary conditions is derived, and the transient heat transfer of a three-layered structure under time-varying heating is studied by using the semianalytical solution. The specific work and the conclusion can be summarized as follows:(1)The time-varying boundary conditions and temperature change functions at the boundary are assumed to be differential and linear.(2)In a small time period, the analytic solutions of each layer of the laminated structure are deduced by using the special separation of variable method.(3)The analytic solutions of each layer in each minute period are obtained by a loop, and then, the semianalytic solutions in the whole time domain are obtained.(4)It is found that when the boundary of the laminated structure is heated by time-varying heat flow, the difference of effective heating area will lead to different amplitude of temperature change. Square wave heating temperature variation amplitude is the largest, sine wave is the second, and triangle wave is the smallest.(5)The larger the heating frequency of the time-varying heat flux is, the smaller the amplitude of temperature fluctuation is.

Nomenclature

:	The number of layers of the laminated structure
:	Heat flux applied to the left side of the laminated structure (W/m²)
:	The temperature on the right side of the laminated structure (°C)
:	Density of ith layer of the laminated structure (kg/m³)
:	Specific heat capacity of ith layer of the laminated structure J/(kg·K)
:	Temperature of ith layer of the laminated structure (°C)
:	Position coordinates (m)
:	Time (s)
:	Thermal conductivity of ith layer of the laminated structure (W/m²·K)
:	The internal heat transfer coefficient of ith layer W/ (m²·K)
:	Base temperature (°C)
:	Heat generation (J)
:	Convective heat transfer coefficient of the surface W/(m²·K)
:	Temperature of cooling fluid (°C)
:	Thickness of the laminated structure (m)
:	Initial temperature (°C)
:	Period of heating (s)
:	The dimensionless representation of ,
:	The dimensionless representation of ,
:	The dimensionless representation of ,
:	The dimensionless representation of ,
:	The dimensionless representation of t,
:	The dimensionless representation of x,
:	The dimensionless representation of ,
:	The dimensionless representation of ,
:	The dimensionless representation of ,
:	The dimensionless representation of ,
:	The dimensionless representation of , .

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the financial supports by the Department of Education of Guangdong Province, under associated Grant nos. 2021KTSCX172 and 2017KTSCX218.

References

M. Wang and J. Yong, “Research progress of multilayer insulation materials,” Material guide, vol. S2, pp. 461–465, 2016.
View at: Google Scholar
Z. K. Wu, F. L. Li, and D. Y. Kwak, “Inverse heat source problem for one-dimensional heat conduction equation based on least-square regularization method,” Computational physics, vol. 1, pp. 49–56, 2016.
View at: Google Scholar
F. Xu, T. J. Lu, K. A. Seffen, and E. Y. K. Ng, “Mathematical modeling of skin bioheat transfer,” Applied Mechanics Reviews, vol. 62, no. 5, pp. 1201–1210, 2009.
View at: Publisher Site | Google Scholar
G. H. Liu and R. H. Xie, “Theoretical analysis and test on influence parameters of air tightness performance of refrigerated vehicles in multi-temperature zone,” Journal of Agricultural Machinery, vol. 48, no. 1, pp. 289–296, 2017.
View at: Google Scholar
Z. S. LiS. Zhu, “Continuous models of heat transfer characteristics of superlattices and layered structures and their applications in energy material design,” Journal of Physics, vol. 56, no. 11, pp. 205–213, 2016.
View at: Google Scholar
C. Y. Li, S. M. Lin, H.-Y. Dai, and D.-L. Liu, “Transient heat transfer analysis of laminated materials based on extended separation variable method,” Journal of Physics, vol. 21, pp. 245–252, 2018.
View at: Google Scholar
M. M. Peiravi and J. Alinejad, “Hybrid conduction, convection and radiation heat transfer simulation in a channel with rectangular cylinder,” Journal of Thermal Analysis and Calorimetry, vol. 140, no. 6, pp. 2733–2747, 2020.
View at: Publisher Site | Google Scholar
F. Xu, T. J. Lu, and K. A. Seffen, “Biothermomechanics of skin tissues,” Journal of the Mechanics and Physics of Solids, vol. 56, no. 5, pp. 1852–1884, 2008.
View at: Publisher Site | Google Scholar
K. C. Liu, H. T. Chen, and P. J. Cheng, “Inverse investigation of non-Fourier heat conduction in tissue,” Journal of Thermal Biology, vol. 21, no. 2, pp. 68–73, 2016.
View at: Publisher Site | Google Scholar
A. Campo and M. Arıcı, “Numerical analysis of unsteady heat conduction in regular solid bodies comprising natural convection to nearby fluids,” Heat Transfer-Asian Research, vol. 47, no. 7, pp. 887–906, 2018.
View at: Publisher Site | Google Scholar
C. Y. Li and S. M. Lin, “Semi-analytical solution of bio-heat conduction for multi-layers skin subjected to laser heating and fluid cooling,” Journal of Mechanics in Medicine and Biology, vol. 17, no. 2, Article ID 1750029, 2017.
View at: Publisher Site | Google Scholar
Z. Yu, X. Shi, S. R. Brohi, H. Qian, F. Liu, and Y. Yang, “Measurement of blood flow in an intracranial artery bypass from the internal maxillary artery by intraoperative duplex sonography,” Journal of Uitrasound in Medicine, vol. 36, no. 2, pp. 439–447, 2017.
View at: Google Scholar
S. Salehi, A. Nori, K. Hosseinzadeh, and D. D. Ganji, “Hydrothermal analysis of MHD squeezing mixture fluid suspended by hybrid nanoparticles between two parallel plates,” Case Studies in Thermal Engineering, vol. 21, Article ID 100650, 2020.
View at: Publisher Site | Google Scholar
M. H. Fu and C. Z. Chen, “Semi-analytical and fine Integral method for three-dimensional temperature field analysis of laminated cylinders,” Journal of Applied Mechanics, vol. 29, no. 1, pp. 15–20, 2012.
View at: Google Scholar
K. Ramadan, “Semi-analytical solutions for the dual phase lag heat conduction in multilayered media,” International Journal of Thermal Sciences, vol. 48, no. 1, pp. 14–25, 2009.
View at: Publisher Site | Google Scholar
J. E. Li, B. L. Wang, and D. M. Chang, “Non-Fourier heat conduction and thermal stress of laminated materials,” Journal of Solid Mechanics, vol. s1, pp. 248–253, 2011.
View at: Google Scholar
F. Liu and W. P. Shi, “Nonlinear heat conduction equations are simulated by lattice Boltzmann method,” Applied Mathematics and Mechanics, vol. 36, no. 11, pp. 1158–1166, 2015.
View at: Google Scholar
L. Shan, J. Liu, G. Lin, and P. Zhang, “Modified scaled boundary finite element analysis of 3D steady-state heat conduction in anisotropic layered media,” International Journal of Heat and Mass Transfer, vol. 108, no. 1, pp. 2462–2471, 2017.
View at: Google Scholar
P. Guo, W. H. Wu, and J. Zhao, “Numerical simulation of heat conduction process of the underlying skin tissue under thermal shock,” Journal of Computational Mechanics, vol. 6, pp. 683–689, 2017.
View at: Google Scholar
L. Li, L. Zhou, and Y. Mo, “An expanded lattice Boltzmann method for dual phase lag model,” International Journal of Heat and Mass Transfer, vol. 93, no. 2, pp. 834–838, 2016.
View at: Publisher Site | Google Scholar
C. Li, G. Y. Feng, and H. M. Yang, “Analytical expression of medium temperature and stress of liquid cooled strip,” Journal of Physics, vol. 5, pp. 123–131, 2016.
View at: Google Scholar
J. Ma, Y. Sun, and J. Yang, “Analytical solution of non-Fourier heat conduction in a square plate subjected to a moving laser pulse,” International Journal of Heat and Mass Transfer, vol. 115, no. 2, pp. 606–610, 2017.
View at: Publisher Site | Google Scholar
M. M. Peiravi, J. Alinejad, D. D. Ganji, and S. Maddah, “3D optimization of baffle arrangement in a multi-phase nanofluid natural convection based on numerical simulation,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 30, no. 5, pp. 2583–2605, 2020.
View at: Google Scholar
S. Kumar and A. Srivastava, “Thermal analysis of laser-irradiated tissue phantoms using dual phase lag model coupled with transient radiative transfer equation,” International Journal of Heat and Mass Transfer, vol. 90, pp. 466–479, 2015.
View at: Publisher Site | Google Scholar
S. M. Lin, “Analytical solutions of bio-heat conduction on skin in Fourier and non-Fourier models,” Journal of Mechanics in Medicine and Biology, vol. 13, no. 4, 2013.
View at: Publisher Site | Google Scholar
S.-M. Lin and C.-Y. Li, “Analytical solutions of non-Fourier bio-heat conductions for skin subjected to pulsed laser heating,” International Journal of Thermal Sciences, vol. 110, no. 7, pp. 146–158, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Xi-Xia Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies