Next Article in Journal
Study on Porosity Defect Detection in Narrow Gap Laser Welding Based on Spectral Diagnosis
Next Article in Special Issue
A Buckling Analysis of Thermoelastic Micro/Nano-Beams Considering the Size-Dependent Effect and Non-Uniform Temperature Distribution
Previous Article in Journal
Computational Materials Design for Ceramic Nuclear Waste Forms Using Machine Learning, First-Principles Calculations, and Kinetics Rate Theory
Previous Article in Special Issue
Homogeneous Zero-Index Thermal Metadevice for Thermal Camouflaging and Super-Expanding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Fourier Heat Conduction of Nano-Films under Ultra-Fast Laser

School of Thermal Engineering, Shandong Jianzhu University, Jinan 250101, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(14), 4988; https://doi.org/10.3390/ma16144988
Submission received: 23 June 2023 / Revised: 8 July 2023 / Accepted: 12 July 2023 / Published: 13 July 2023

Abstract

:
The ultra-fast laser heating process of nano-films is characterized by an ultra-short duration and ultra-small space size, in which the classical Fourier law based on the hypothesis of local equilibrium is no longer applicable. Based on the Cattaneo–Vernotte (CV) model and the dual-phase-lag (DPL) model, the two-dimensional analytical solutions of heat conduction in nano-films under ultra-fast laser are obtained using the integral transformation method. The results show that there is a thermal wave phenomenon inside the film, which becomes increasingly evident as the elapse of the lag time of the temperature gradient. Moreover, the wave amplitude in the vertical direction is much larger than that in the horizontal direction of the nano-film. By comparing the numerical result of the two models, it is found that the temperature distribution inside the nano-film based on the DPL model is gentler than that of the CV model. Additionally, the temperature distribution in the two-dimensional solution is lower than that in the one-dimensional solution under the same Knudsen number. In the comparison results of the CV model, the maximum peak difference in the thermal wave reaches 75.08 K when the Knudsen number is 1.0. This demonstrates that the horizontal energy carried by the laser source significantly impacts the temperature distribution within the film.

1. Introduction

Micro/nano-thin films are widely used in the fields of micro-electronics, energy storage and biological medicine [1,2,3,4,5,6], among which the preparation of nano-thin films is an important technical support. The pulse duration of an ultrafast laser is from nanoseconds to femtoseconds, with the unique advantages of a short action time and a high energy density [7]. The interaction of an ultra-fast laser with a nano-thin film can solve the thermal instability during the preparation of nano-thin films and improve the preparation efficiency. Therefore, the investigation of the heat conduction in nano-thin films under the action of ultra-fast laser is of great practical significance.
The Fourier’s law, which was proposed by French scientist Fourier in 1822, established the theoretical foundation of heat conduction. It has been proved that the numerical results based on the classical Fourier law have been very appropriate to solve most heat conduction problems with low heat flux, a long thermal action time and a macro scale [8,9,10]. Chen et al. [11] investigated the repeated long-pulse laser heating of solid materials based on Fourier’s law. The research results provide a theoretical foundation and guidance for the heat conduction mechanism of repeated long-pulse lasers. Turkyilmazoglu [12] proposed a Fourier theoretical model of a laser pulse heating process for two-dimensional cylindrical materials. The surface temperature increases and attenuation inside the material is detected by irradiating the pulsed laser. The results show that the temperature will rise with the increase in the initial temperature of the material or the pulse width of the laser.
However, Fourier’s law is based on the assumption that the velocity of heat conduction is infinitely great [13]. In some special conditions, such as extremely small characteristic sizes and extremely fast process times, the speed of heat conduction is finite and Fourier’s law may no longer apply [14]. To solve the problem, the Cattaneo–Vernotte (CV) model was proposed based on the hypothesis that the speed of heat conduction is limited [15,16,17]. The model first introduced the thermal wave phenomenon and the heat flux hysteresis time τq was introduced. Different from Fourier’s law, a hysteresis effect of temperature gradient is produced at the time t in a certain area of the system, where a heat flow vector appears at time t + τq instead of time t [18]. Through this improvement, the CV model has been widely used in the fields of biology medicine [19], ultra-fast lasers [20,21] and microelectronics [22,23]. Zhang et al. [24] used the CV model to investigate the range of heat-affected under step pulse actions in a semi-infinite space. The relationship between the heat distribution and time is compared using the Fourier model and the non-Fourier model. Noroozi et al. [25] applied the CV model to investigate the heat interaction between ultra-fast lasers and thin films. They concluded that the non-linearity analysis is very significant in non-Fourier heat conduction. The observation that there is a significant difference between the Fourier and non-Fourier results emphasizes the significance of a non-Fourier model for similar problems. But the CV model only considers the heat flux density vector phase retardation, so the problem of insufficient physical basis still exists [26,27,28]. Later, Tzou [29,30] proposed the DPL model to reveal the microscopic heat conduction inside the medium from a macroscopic perspective, and to explain the hysteresis time of the heat flux vector and temperature gradient. The DPL model eliminates the assumptions in the CV model, and allows the temperature gradient ahead of the heat flux vector or the heat flux vector ahead of the temperature gradient in the process of the transient state [31]. Based on the DPL model, Majchrzak [32] calculated the internal thermal processes of cylindrical microdomains (Cr, Au) under the action of an ultra-short laser pulse and introduced the numerical model of the melting and re-solidification of the materials. Xu et al. [33] proposed an effective numerical method to solve the complexity of DPL equations. The 1-D, 2-D and 3-D models in single-phase media under laser pulse were numerically simulated. The results showed that the method was efficient in predicting the thermal behavior of laser pulse heating and other non-Fourier heat conduction problems. A.J. van der Merwe [34] used modal analysis to compare the DPL, CV and Fourier models; the results showed that the difference between the CV and the DPL is related to the values of τT and τq. In addition, numerical simulation methods such as the Lattice Boltzmann method (LBM), molecular dynamics (MD) and the Monte Carlo (MC) simulation have also been applied to non-Fourier heat transport problems [35,36,37,38,39]. Mao and Xu [40] investigated the one-dimensional heat conduction of ultrafast laser-heated nano-films using the LBM model and the improved CV model, respectively. The results showed that there was an obvious thermal wave phenomenon in the film, which indicates that both the analytical and numerical methods can show the thermal wave phenomenon inside the film. Although many methods are used in non-Fourier heat conduction problems, the analysis mainly focuses on 1-D. To further investigate the mechanism of non-Fourier heat conduction, the 2-D model is built using the CV and DPL models. And, the effect of relaxation time on non-Fourier heat conduction is analyzed.
This paper consists of the following parts: In Section 2, the formulas and basic models used in this paper are introduced. In Section 3, the mathematical model selected in this paper is first verified. Then, the calculation results based on the CV and DPL models are obtained, and the non-Fourier heat conduction process under different characteristic sizes of nano-film is described in the two-dimensional case. Next, the influence of hysteresis time on the heat conduction process is analyzed in the DPL model. Finally, the results of the 1-D and 2-D models are compared. In Section 4 and Section 5, the discussion and conclusion are given, respectively.

2. Method

2.1. Physical Model

As shown in Figure 1, the ultra-fast laser is assumed to be irradiated vertically from the center of the film. The influence radius of laser heat source and longitudinal penetration depth are r0 and δ, respectively. These two parameters determine the range inside the film into which the ultra-fast laser energy can be radiated. The vertical size and horizontal size of the film are H and Lr, respectively. Where Lr = 750 nm, H = l/Kn. Kn is the Knudsen number, a ratio of mean free path of the energy carriers l to the vertical size H. To simplify the analysis, the film is assumed to be an adiabatic boundary.

2.2. Mathematical Model

The CV model and the DPL model are denoted by Equations (1) and (2), respectively [15,16,17,29,30].
q r , t + τ q = k T r , t ,
q r , t + τ q = k T r , t + τ T .
where, k is the thermal conductivity of medium, T is the temperature gradient which is a vector function of the time variable t and the position vector r, τT is the phase lag time of T and τq is the phase lag time of the heat flux vector q.
The derivation process of the more complex DPL model is shown below. Because the derivation process of the CV model is similar to the DPL model, it is not explicitly outlined here. Combining Equation (2) with the energy conservation equation, the governing equation based on the DPL model can be expressed as Equation (3):
k 2 T r , z , t r 2 + 1 r T r , z , t r + 2 T r , z , t z 2 + τ T 3 T r , z , t r 2 t + 1 r 2 T r , z , t r t + 3 T r , z , t z 2 t + S r , z , t + τ q S r , z , t t = C v T r , z , t t + C v τ q 2 T r , z , t t 2 .
The definite conditions are expressed as follows:
T r , z , t z / z = 0 = 0 , T r , z , t z / z = H = 0 , T r , z , t r / r = L r = 0 ,
T r , z , t / t = 0 = T 0 , T r , z , t t / t = 0 = S r , z , t C v .
where Cv is volumetric heat capacity and S(r,z,t) is the energy absorption rate of the ultra-fast laser which can be expressed as follows:
S r , z , t = 0.94 J 1 R t p δ exp z δ r 2 r 0 2 1.992 t t p .
where J is the ultra-fast laser energy density, tp is the laser pulse duration and R is the surface reflectivity. To solve the Equations (3)–(5), the dimensionless parameters are introduced as Equation (7):
z * = z H , r * = r L r , T * = T T 0 T 0 , t * = t t p ,
Then the dimensionless governing equations are obtained:
a 2 T * r * , z * , t * t * 2 + T * r * , z * , t * t * = k t p C v 1 L r 2 2 T r * , z * , t * r * 2 + 1 r * T r * , z * , t * r * + 1 H 2 2 T r * , z * , t * z * 2 + a b 1 L r 2 t p 2 T r * , z * , t * r * 2 + 1 r * T r * , z * , t * r * + 1 H 2 t p 3 T r * , z * , t * z * 2 t * + 1 1.992 a S 0 Q * r * , z * , t * ,
Q * r * , z * , t * = exp ζ z * ε * 2 r * 2 1.992 t * ,
S v = 0.94 J 1 R T 0 C v δ , μ = H δ , ε = L r r , a = τ q t p , b = τ T τ q .
where Sv, μ , ξ and ε are the transition parameters with no special meaning; a is the ratio of τq to tp, b is the ratio of τT to τq. Then, the definite solution conditions can be expressed as follows:
T * r * , z * , t * z * / z * = 0 = 0 , T * r * , z * , t * z * / z * = 1 = 0 , T * r * , z * , t * r * / r * = 1 = 0 ,
T * r * , z * , t * / t * = 0 = 0 ; T * r * , z * , t * t * / t * = 0 = S 0 Q * r * , z * , t * .
Equations (8)–(10) are solved using the integral transformation method. In order to distinguish, one overline is added above the expression symbol ( T ¯ ) to symbolize the integral transformation of one variable, and two overlines symbolize the integral transformation of two variables ( T * ¯ ¯ ). The forward integral transformation and inverse integral transformation of r* are shown as:
T * r * , z * , t * = m = 0 J 0 β m r * N β m T ¯ β m , z * , t * ,
T ¯ β m , z * , t * = 0 1 r * J 0 β m r * T r * , z * , t * d r * .
where J 0 β m r * is the Bessel functions of order zero of the first class and β m is the root of J 1 β m r * = 0 , m = 1 , 2 , 3 , N β m = J 0 2 β m / 2 , β m 0 . When m = 0 and β 0 = 0 , N β 0 = 1 / 2 . Then the integral transformation of z* are shown as follows:
T * ¯ β m , z * , t * = n = 1 Z η n z * N η n T * ¯ ¯ β m , η n , t * ,
T * ¯ ¯ β m , η n , t * = 0 1 Z * η n z * T ¯ * β m , z * , t * .
where n = 0 , 1 , 2 , η n = n π , N η n = 1 / 2 , η n 0 1 , η n = 0 .
By substitution and simplification, the dimensionless governing Equations (7)–(10) can be expressed as follows:
a 2 T * ¯ ¯ β m , η n , t * t * 2 + 1 + a b γ m n T * ¯ ¯ β m , η n , t * t * + γ m n T * ¯ ¯ β m , η n , t * = 1 1.992 a S 0 Q * ¯ ¯ β m , η n , t * ,
γ m n = β m 2 L r 2 + η n 2 H 2 k t p C v ,
Q * ¯ ¯ β m , η n , t * = 0 1 r * J 0 β m r * exp ε * 2 r * 2 d r * 0 1 cos η n z * exp ζ z * 1.992 t * d z * ,
T * ¯ ¯ β m , η n , t * / t * = 0 = 0 ; T * ¯ ¯ β m , η n , t * t * / t * = 0 = S 0 Q * ¯ ¯ β m , η n , t * .
Equation (17) is a second-order nonhomogeneous linear ordinary differential equation; then, the expression of T * ¯ ¯ r * , z * , t * can be obtained as:
T * ¯ ¯ r * , z * , t * = m = 0 n = 0 J 0 β m r * cos η n z * N β m N η n C 1 F 1 β m , η n , t * + C 2 F 2 β m , η n , t * + C exp 1.992 t * .
Let Δ = (1 + abγmn)2 − 4mn, the function and parameter in the Equation (21) can be expressed as:
C = 1 0.498 a S 0 0 1 r * J 0 β m r * exp ζ z * d r * 0 1 cos η n z * exp ε * 2 r * 2 d z * 1.992 2 a 1.992 1 + a b γ m n + γ m n .
(1) When Δ > 0,
F 1 β m , η n , t * = e 1 + a b γ m n + Δ 2 a t * , F 2 β m , η n , t * = e 1 + a b γ m n Δ 2 a t * ,
(2) When Δ = 0,
F 1 β m , η n , t * = e 1 + a b γ m n 2 a t * , F 2 β m , η n , t * = t * e 1 + a b γ m n 2 a t * ,
(3) When Δ < 0,
F 1 β m , η n , t * = e 1 + a b γ m n 2 a t * cos Δ 2 a , F 2 β m , η n , t * = e 1 + a b γ m n 2 a t * sin Δ 2 a .
C 1 = C F 2 β m , η n , t * t * 1.992 C + S 0 0 1 cos η n z * exp ζ z * d r * 0 1 r * J 0 β m r * exp ε * 2 r * 2 d z * F 2 β m , η n , t * F 1 β m , η n , t * F 2 β m , η n , t * t * F 2 β m , η n , t * F 1 β m , η n , t * t * C 2 = C F 1 β m , η n , t * t * + 1.992 C + S 0 0 1 cos η n z * exp ζ z * d r * 0 1 r * J 0 β m r * exp ε * 2 r * 2 d z * F 1 β m , η n , t * F 1 β m , η n , t * F 2 β m , η n , t * t * F 2 β m , η n , t * F 1 β m , η n , t * t * .
Substituting Equation (21) into Equations (13) and (15), T * r * , z * , t * can be obtained as follows:
T * r * , z * , t * = m = 0 n = 0 J 0 β m r * cos η n z * N β m N η n C 1 F 1 β m , η n , t * + C 2 F 2 β m , η n , t * + C exp 1.992 t * .

2.3. Cases Descriptions

Nano silicon film is selected in the present work, the physical parameters of nano-silicon thin film and an ultra-fast laser are shown in Table 1. These parameters are used in all the following calculations.
According to the above formula derivation and specific parameters, the CV and DPL models can be applied to calculate the heat conduction process inside the nano-film under the ultra-fast laser. In Table 1, the laser pulse duration tp = 0.65 ps. Combined with Equation (7), it can be seen that the laser pulse starts to launch when t* = 0, and ends when t* = 1. In order to show the heat conduction process inside the film as completely as possible, the transient temperature distributions at different dimensionless time t* are shown in the following results. In addition, Knudsen number plays an important role in terms of the scale effect, the increase in Kn (Kn = l/H) means the decrease in the vertical size (H) when the mean free path of phonons (l) is invariable. In the present work, the cases of Kn = 0.1, 0.5, 1 and 2 are mainly calculated around the vertical size of the film equivalent to the mean free path of phonons, and the corresponding vertical sizes (H) are 410 nm, 82 nm, 41 nm and 21 nm.

3. Analysis of Results

3.1. Authentication

To validate the methodology in the present work, a simple one-dimensional silicon film heat transfer model is established and verified with Escobar [36]. To be consistent with the mathematical model, the laser heat source is removed and the dimensionless transformation changes to the same as Chapter four in Ref [36]. Where the dimensionless temperature T* is defined as (T − T2)/(T1 − T2), T is the temperature at a random point in the system, T1 is the hot boundary temperature set at 301 K and T2 is the cold boundary temperature set at 300 K. Then, the same problem is solved by the DPL model in the present work.
Figure 2a shows the transient temperature distribution at Kn = 0.1 in the diffusive state. All three results show that the diffusion behavior of heat transport is due to the large characteristic length and time scale of the system. The present work is well agreed with the Cattaneo and LBM, and the classical Fourier law also applies in this case [36]. When Kn = 1, the heat transfer changes to a transitional state as shown in Figure 2b. The results of LBM show that the temperature jumps at both ends of the film and the linear temperature distribution inside the film. The present work and Cattaneo show the phenomenon of discontinuities in the temperature distribution, which is one of the characteristics of non-Fourier heat conductions under the CV or DPL model. Compared with LBM, there are some differences in the present work. This is because the error exists in the discretization process during the implementation of the LBM numerical method, while the exact analytical solution is absolutely error-free because it is independent of the size of the domain. The above non-Fourier thermal conductivity models such as Cattaneo and DPL can be viewed as being based on the assumption that there is a time lag between heat flow and temperature gradient, so the results of the present work are more closely related to the Cattaneo. Therefore, the graph truly proves the accuracy of the current mathematical model.

3.2. Result of the CV Model

In this section, the CV model is applied to analyze the heat conduction process of ultra-fast laser heating nano-thin films. As can be seen from Figure 3a, after the laser irradiates the film, the temperature of the surface subjected to the laser action rises rapidly, then decreases with the elapse of dimensionless time t*, and the heat is gradually transferred to the inside of the film. With the elapse of dimensionless time, the peak temperature point is no longer on the heated surface but appears inside the film. This phenomenon means that the heat inside the nano-thin film is transferred in the manner of a thermal wave. This is because the heat in the film cannot be diffused in time as a result of the hysteresis of the heat flux in the heating process, and the accumulation of heat can produce the phenomenon of a thermal wave. The CV model indicates that the heat transfers not in the manner of diffusion but in a limited speed, which offers a clear difference between the CV model and the Fourier model. In addition, the temperature drops rapidly after reaching the peak, such as in Figure 3a at t* = 5. This is because the energy for the temperature increase in the area near the heated surface is mainly from the laser heat source. With the depth increases in the vertical, the absorption of the laser energy decreases while the heat conduction plays a dominant role. The impact of the laser source becomes less and less as the depth is further expanded, so the heat conduction inside the film is the main mechanism controlling the temperature increase in this region.
With Kn increases from 0.1 to 0.5, as shown in Figure 3b, the thermal wave inside the film becomes gentler compared with Kn = 0.1, and the temperature inside the film is close to equilibrium at t* = 100. Figure 3c shows that the temperature peak inside the film reaches the bottom at t* = 20, and then moves in an opposite direction at t* = 30 due to the adiabatic boundary conditions. The phenomenon occurs repeatedly until the temperature inside the film reaches a stable state, although, at t* = 100 in Figure 3d, there is still a thermal wave slightly. Compared with Figure 3a–d, it can be seen that the propagation of the thermal wave inside the film becomes gentler, and the time for the temperature to approach the equilibrium state is shorter. This is because the vertical size of the film decreases as Kn increases and each unit length of the film absorbs more energy at the same time.
Figure 4 shows the transient temperature distribution in the horizontal direction corresponding to Figure 3. The horizontal size is 750 nm invariably, which does not change with Kn. As can be seen from Figure 4, the heat in the horizontal direction of the film is also transferred in the manner of a thermal wave, but the amplitude is smaller than that in the vertical. It can be seen from Equation (6) that the energy transmitted using the ultra-fast laser in the horizontal is far lower than that in the vertical direction, resulting in the heat accumulation per unit length of the film in the horizontal being much less than that in the vertical. Hence, the thermal wave amplitude in the horizontal direction is small. Additionally, when there is no thermal wave inside the film, it does not mean the heat is not transferred in the manner of a thermal wave, but the energy accumulation is not enough to make the value of the next position higher than the previous one.
From Figure 4a–d, it can be found that the horizontal size affected by the heat increases gradually with the elapse of dimensionless time, but it has not been transferred to the horizontal boundary of the film. This is because the heat transferred in the horizontal direction is much smaller than in the vertical direction. In addition, the horizontal size of the film is too large, resulting in the heat arriving at equilibrium before reaching the horizontal boundary. As is shown in Figure 3c, the transient temperature profile inside the film is transferred in an opposite direction to the vertical direction at t* = 50, which will also affect the horizontal temperature. As can be seen in Figure 4c, the phenomenon of a thermal wave in the horizontal direction disappears at t* = 50, and then the temperature at t* = 50 exceeds that at t* = 30 again. Therefore, the heat transfer in the vertical and horizontal direction will influence each other after the laser irradiates the thin film.

3.3. Result of the DPL Model

The results using the DPL model are analyzed below. Figure 5 shows the transient temperature distribution in the vertical direction, and Figure 6 shows it in the horizontal direction corresponding to Figure 5, where b = 0.2 and different Kns are chosen. Compared with the CV model, the temperature peak obtained based on the DPL model is relatively low under the same Kn and t*, the thermal wave is gentler and the propagation speed is relatively slow. That is, the peak position is slightly closer to the laser action surface of the DPL results than in that of the CV results. The reason is that the DPL model considers the hysteresis time of the heat flux and the temperature gradient between hot carriers.
Without considering the external heat source, the temperature distribution of the DPL model at dimensionless time t* is equal to the temperature distribution of the CV model at dimensionless time t* × (1 + b). The detailed explanation is that the heat at the peak temperature of the film is transferred to other areas when the laser irradiates the film. Because of the influence of the hysteresis of the temperature gradient of the DPL model, the heat at peak temperature is transferred more in the film than in that of the CV model. Therefore, compared with the CV model, the temperature peak of the DPL model is lower but the temperature near the peak is higher. In other words, with the dimensionless time elapse, the thermal wave from the DPL model is relatively gentle and the position of the temperature peak changes slowly in the vertical direction.

3.4. Influence of the Hysteresis Time of the Temperature Gradient

This section analyzes the influence of the hysteresis time of the temperature gradient on the heat conduction process by changing the b value, which is the ratio of τT to τq as shown in Equation (10). Figure 7 shows the results of different data b at Kn = 0.1 and 0.5, respectively. As is shown in Figure 7, the peak temperature becomes lower and lower with the increase in b value. The thermal wave phenomenon gradually disappears and the temperature distribution becomes gentler. This is because the hysteresis effect of the temperature gradient is more evident after the increase in b value. When b = 0, there is only the hysteresis of the heat flux vector inside the film, which implies that results are equivalent to the CV model. When b > 0, the temperature distribution will lag by b × t* dimensionless time because of the hysteresis time of the temperature gradient. Therefore, the temperature distribution becomes gentler following the increase in b value under the same t*. In addition, there is no thermal wave phenomenon in the horizontal direction as shown in Figure 8 due to the heat transferred by the ultra-fast laser source in the horizontal direction being lower than that in the vertical direction. This indicates that the appearance of the thermal wave is not only related to the hysteresis time of the temperature gradient but is greatly related to the amount of the energy transferred by the laser source.

3.5. Comparative between 2-D and 1-D Models

Figure 9 and Figure 10 show the comparison of temperature distribution in the vertical direction of the film between the 2-D and 1-D results of the CV model and the DPL model, respectively, where b = 0.2 in the DPL model. As can be seen from Figure 9, there is no evident difference between the 2-D and 1-D results at t* = 0.2 and 1.5 owing to the short heat conduction process; the ultra-fast laser source mainly acts on the vertical direction at the moment. When t* = 17.32, the difference between the two results become obvious, and the maximum difference is 120.42 K and 124.65 K in Figure 9a,b, respectively, which appears at the left boundary of the film. Although the peak of the thermal wave is reached at the same time (t* = 17.32) and position (z = 30 nm), the 2-D results are lower than those in the 1-D results and the peak difference of the thermal wave is 70.63 K and 75.08 K in Figure 9 (a) and (b), respectively. This is because the 2-D model takes into account the effect of the laser on the horizontal direction of the film. The conduction of heat in the horizontal direction is small at the initial time, but more and more heat in the 2-D model will transfer to the horizontal direction with the elapse of dimensionless time. In addition, as the Kn increases, the film contains more heat, and the difference between the 2-D and 1-D model also increases.
Figure 10 shows that the difference between the 2-D and 1-D of the DPL model is slightly greater than that of the CV model. This implies that the hysteresis time of the temperature gradient makes the influence more pronounced in the horizontal direction for the heat transfer process. By comparing the 2-D and 1-D results, it is found that the heat transfer of the horizontal dimension has extraordinary impact on the temperature distribution inside the film when the heat conduction process takes a long time. The accuracy and objectivity of the 2-D model are verified, and it is also illustrated that the horizontal heat transfer cannot be ignored in the investigation of heat transfer in thin films, nor can it be simplified into a 1-D solution.

4. Discussion

From the above results, we conclude that the heat is transferred in the manner of a thermal wave inside the nano-film under an ultra-fast laser based on the non-Fourier model. This conclusion is consistent with previous work on the Lattice Boltzmann method [40]. Different from the Lattice Boltzmann method, the non-Fourier model in the present work can investigate the influence of the hysteresis time of heat flux and temperature gradient on the heat transfer process. And, the heat conduction dimension is extended from 1-D to 2-D, which reveals that the heat inside the system is greatly overestimated in the 1-D model. Normally, the experiment is the most essential method to obtain the information of heat transfer inside the film, but it is almost impossible to perform a complete experiment due to its ultra-small characteristics’ size and ultra-fast physical processes. Therefore, theoretical research and simulation analysis constructed on valid mathematical models become an alternative method to exploring the mechanisms of the heat transfer process. To overcome the limitations, the next stage of our work involves conducting experiments for verification and comparison.

5. Conclusions

Based on the DPL model and the CV model, the non-Fourier heat conduction process of nano-films under ultra-fast laser action is investigated. Both the DPL and the CV models can show that the heat inside the nano-film is transferred in thermal waves after laser irradiation, but the temperature distribution based on the DPL model is gentler than that of the CV model.
With the decrease in the size film, the propagation of the thermal wave inside the film becomes gentler and the time for the temperature to approach the equilibrium state is shorter. And, the thermal wave in the vertical direction is more obvious than that in the horizontal direction of the film.
When the hysteresis time of the heat flux is invariable, the greater the hysteresis time of the temperature gradient and the gentler the temperature distribution. This indicates that the thermal wave phenomenon will appear in the film when the hysteresis time of the temperature gradient is less than the hysteresis time of the heat flux.
The result of the 1-D model is substantially higher than that of the 2-D model after a long period of heat conduction due to the disregard of the horizontal heat transfer. Furthermore, with the decrease in the size film, the differences will be more evident. This study verifies the accuracy and objectivity of the 2-D model and also illustrates that the horizontal heat transfer cannot be neglected in the research of heat conduction in nano-films.

Author Contributions

Conceptualization, Y.M. and K.Y.; Methodology, Y.M., J.L., M.Y. and K.Y.; Software, Y.M.; Validation, Y.M. and X.L.; Investigation, S.L., J.L., M.Y. and X.L.; Writing—original draft, Y.M., S.L. and K.Y.; Writing—review and editing, Y.M., J.L., M.Y. and K.Y.; Supervision, J.L., M.Y. and K.Y.; Project administration, Y.M.; Funding acquisition, K.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Natural Science Foundation of Shandong Province (No. ZR2017BEE052), Sub-project of National Key R&D Program (No. 2016YFC0700803-01). This work was also supported by the Plan of Introduction and Cultivation for Young Innovative Talents in Colleges and Universities of Shandong Province, and the scientific and technological project of Suzhou City in 2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Catania, F.; Oliveira, H.; Lugoda, P. Thin-film electronics on active substrates: Review of materials, technologies and applications. J. Phys. D Appl. Phys. 2022, 55, 323002. [Google Scholar] [CrossRef]
  2. Abhik, S.; Sudip, M.; Kanoria, M. Effect of nonlocality in the vibration of a microscale ribbon due to laser pulse. Waves Random Complex Media 2020, 32, 2751–2775. [Google Scholar]
  3. Badgujar, A.C.; Dusane, R.O.; Dhage, S.R. Solution-processed CIGS thin film solar cell by controlled selenization process. Mater. Today Proc. 2022, 52, 829–833. [Google Scholar] [CrossRef]
  4. Abouabassi, K.; Sala, A.; Atourki, L. Electrodeposited CuSbSe2 thin films based solar cells on various substrates. J. Nanopart. Res. 2022, 24, 221. [Google Scholar] [CrossRef]
  5. Al-Qaissi, H.; Al-qarni, Y.; Al-Aniezi, A. Optical characterizations of synthetic nano PcGaCl thin film. In vitro, hemocompatibility, cytotoxicity, and PDT evaluation. Dig. J. Nanomater. Biostruct. 2022, 17, 861–869. [Google Scholar] [CrossRef]
  6. Asgar, A.; Barman, A.; Sanatan, D. Electromagnetic phenomena in cilia actuated peristaltic transport of hybrid nano-blood with Jeffrey model through an artery sustaining regnant magnetic field. Waves Random Complex Media 2022, 1–32. [Google Scholar]
  7. He, S.Q.; He, T.H. Investigation on the transient response of a porous half-space with strain and thermal relaxations due to laser pulse heating. Waves Random Complex Media 2021, 1–25. [Google Scholar] [CrossRef]
  8. Joseph, D.D.; Preziosi, L. Addendum to the Paper Thermal waves. Rev. Mod. Phys. 1990, 62, 375–391. [Google Scholar] [CrossRef]
  9. Jou, D.; Casas-Vazquez, J.; Lebon, G. Extended Irreversible Thermodynamics; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  10. Ghazanfarian, J.; Shomali, Z.; Abbassi, A. Macro-to nanoscale heat and mass transfer: The lagging behavior. Int. J. Thermophys. 2015, 36, 1416–1467. [Google Scholar] [CrossRef]
  11. Chen, G.; Wang, Y.D.; Zhang, J.J. An analytical solution for two-dimensional modeling of repetitive long pulse laser heating material. Int. J. Heat Mass Transf. 2017, 104, 503–509. [Google Scholar] [CrossRef]
  12. Turkyilmazoglu, M. An analytic parametric study of rounded laser pulse heating. Optik 2021, 240, 167001. [Google Scholar] [CrossRef]
  13. Mishra, S.C.; Sahai, H. Analyses of non-Fourier heat conduction in 1-D cylindrical and spherical geometry—An application of the lattice Boltzmannmethod. Int. J. Heat Mass Transf. 2012, 55, 7015–7023. [Google Scholar] [CrossRef]
  14. Frankel, J.I.; Vick, B.; Özisik, M.N. General formulation and analysis of hyperbolic heat conduction in composite media. Int. J. Heat Mass Transf. 1987, 30, 1293–1305. [Google Scholar] [CrossRef]
  15. Cattaneo, C. A form of heat conduction equation which eliminates the paradox of instantaneous propagation. Compte Rendus 1958, 247, 431–435. [Google Scholar]
  16. Cattaneo, C. Sulla conduzione del calore. In Atti del Seminario matematico e Fisico dell’Università di Modena; Società tipografica Modenese: Modena, Italy, 1948; Volume 3, pp. 83–101. [Google Scholar]
  17. Vernotte, P. Some possible complications in the phenomena of thermal conduction. Compte Rendus 1961, 252, 2190–2193. [Google Scholar]
  18. Liu, Y.; Li, L.; Lou, Q. A hyperbolic lattice Boltzmann method for simulating non-Fourier heat conduction. Int. J. Heat Mass Transf. 2019, 131, 772–780. [Google Scholar] [CrossRef]
  19. Lin, S.M.; Li, C.Y. Analytical solutions of non-Fourier bio-heat conductions for skin subjected to pulsed laser heating. Int. J. Therm. Sci. 2016, 110, 146–158. [Google Scholar] [CrossRef]
  20. Xu, G.Y.; Wang, J.B. Analytical solution of time fractional Cattaneo heat equation for finite slab under pulse heat flux. Appl. Math. Mech. 2018, 39, 1465–1476. [Google Scholar] [CrossRef]
  21. Lenarczyk, M.; Domanski, R. Non-Fourier effects in the thermal protection against high-power ultra-fast laser pulses. HTRSE-2018 2018, 70, 03008. [Google Scholar] [CrossRef] [Green Version]
  22. Faouzi, N.; Hichem, G.; Mohamed, F.B.A. Influence of Joule effect on thermal response of nano FinFET transistors. Superlattices Microstruct. 2021, 156, 106980. [Google Scholar]
  23. Nasri, F.; Ben Aissa, M.F.; Belmabrouk, H. Microscale thermal conduction based on Cattaneo-Vernotte model in Silicon On Insulator and Double Gate MOSFETs. Appl. Therm. Eng. 2015, 76, 206–211. [Google Scholar] [CrossRef]
  24. Zhang, X.M.; Peng, S.; Zhang, L. Variational equation of non-fourier heat conduction. Heat Transf. Res. 2018, 49, 275–285. [Google Scholar] [CrossRef]
  25. Noroozi, M.J.; Saedodin, S.; Ganji, D.D. A new approximate-analytical method to solve non-Fourier heat conduction problems. Kuwait J. Sci. 2017, 44, 87–96. [Google Scholar]
  26. Joseph, D.D.; Preziosi, L. Heat waves. Rev. Mod. Phys. 1989, 61, 41–73. [Google Scholar] [CrossRef]
  27. Rubin, M.B. Hyperbolic heat conduction and the second law. Int. J. Eng. Sci. 1992, 30, 1665–1676. [Google Scholar] [CrossRef]
  28. Yang, Z.B.; Wang, Z.K.; Tian, S.H. Analysis and Modelling of Non-Fourier Heat Behavior Using the Wavelet Finite Element Method. Materials 2019, 12, 1337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Tzou, D.Y. The generalized lagging response in small-scale and high-rate heating. Int. J. Heat. Mass Transf. 1995, 38, 3231–3234. [Google Scholar] [CrossRef]
  30. Tzou, D.Y. Macro-to Microscale Heat Transfer: The Lagging Behavior; Taylor & Francis: Washington, DC, USA, 1997. [Google Scholar]
  31. Xu, M.; Wang, L. Thermal oscillation and resonance in dual-phase-lagging heat conduction. Int. J. Heat Mass Transf. 2002, 45, 1055–1061. [Google Scholar] [CrossRef]
  32. Majchrzak, E.; Mochnacki, B. Numerical Simulation of Thermal Processes in a Domain of Thin Metal Film Subjected to an Ultrashort Laser Pulse. Materials 2018, 11, 2116. [Google Scholar] [CrossRef] [Green Version]
  33. Xu, B.B.; Gao, X.W.; Cui, M. High precision simulation and analysis of non-Fourier heat transfer during laser processing. Int. J. Heat Mass Transf. 2021, 17, 121574. [Google Scholar] [CrossRef]
  34. Van der Merwe, A.J.; Van Rensburg, N.F.J.; Sieberhagen, R.H. Comparing the dual phase lag, Cattaneo-Vernotte and Fourier heat conduction models using modal analysis. Appl. Math. Comput. 2021, 698, 125934. [Google Scholar] [CrossRef]
  35. Chen, B.Q.; Kan, A.K.; Chen, Z.F. Investigation on Effective Thermal Conductivity of Fibrous Porous Materials as Vacuum Insulation Panels’ Core Using Lattice Boltzmann Method. Energies 2023, 19, 3692. [Google Scholar] [CrossRef]
  36. Escobar, R. Lattice Boltzmann Modeling of Phonon Transport in Silicon Films. Ph.D. Thesis, Carnegie Mellon University, Pittsburgh, PA, USA, 2005. [Google Scholar]
  37. Chen, J.; Zhang, G.; Li, B.W. Molecular Dynamics Simulations of Heat Conduction in Nanostructures: Effect of Heat Bath. J. Phys. Soc. Jpn. 2010, 79, 074604. [Google Scholar] [CrossRef] [Green Version]
  38. Kumar, S.R.; Reddy, B.K.; Balaji, C.A. Application of Hybrid Monte Carlo Algorithm in Heat Transfer. J. Heat Transf. 2017, 139, 082004. [Google Scholar] [CrossRef]
  39. Palluotto, L.; Dumont, N.; Rodrigues, P. Assessment of randomized Quasi-Monte Carlo method efficiency in radiative heat transfer simulations. J. Quant. Specteosc. Radiat. Transf. 2019, 236, 106570. [Google Scholar] [CrossRef]
  40. Mao, Y.D.; Xu, M.T. Lattice Boltzmann numerical analysis of heat transfer in nano-scale silicon films induced by ultra-fast laser heating. Int. J. Therm. Sci. 2015, 89, 210–221. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of ultra-fast laser-heating nano-films.
Figure 1. Schematic diagram of ultra-fast laser-heating nano-films.
Materials 16 04988 g001
Figure 2. Validation of present work with different discrete data points of published numerical model [36]. (a) Kn = 0.1; (b) Kn = 1.
Figure 2. Validation of present work with different discrete data points of published numerical model [36]. (a) Kn = 0.1; (b) Kn = 1.
Materials 16 04988 g002
Figure 3. The transient temperature profile inside the film at the center of the action surface in the vertical direction based on the CV model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Figure 3. The transient temperature profile inside the film at the center of the action surface in the vertical direction based on the CV model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Materials 16 04988 g003
Figure 4. The transient temperature profile inside the film at the center of the action surface in the horizontal direction based on the CV model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Figure 4. The transient temperature profile inside the film at the center of the action surface in the horizontal direction based on the CV model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Materials 16 04988 g004
Figure 5. The transient temperature profile inside the film at the center of the action surface in the vertical direction based on the DPL model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Figure 5. The transient temperature profile inside the film at the center of the action surface in the vertical direction based on the DPL model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Materials 16 04988 g005
Figure 6. The transient temperature profile inside the film at the center of the action surface in the horizontal direction based on the DPL model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Figure 6. The transient temperature profile inside the film at the center of the action surface in the horizontal direction based on the DPL model. (a) Kn = 0.1; (b) Kn = 0.5; (c) Kn = 1; (d) Kn = 2.
Materials 16 04988 g006aMaterials 16 04988 g006b
Figure 7. The transient temperature profile inside the film in the vertical direction at different b values when t* = 10. (a) Kn = 0.1; (b) Kn = 0.5.
Figure 7. The transient temperature profile inside the film in the vertical direction at different b values when t* = 10. (a) Kn = 0.1; (b) Kn = 0.5.
Materials 16 04988 g007
Figure 8. The transient temperature profile inside the film in the horizontal direction at different b values when t* = 10 at Kn = 0.1.
Figure 8. The transient temperature profile inside the film in the horizontal direction at different b values when t* = 10 at Kn = 0.1.
Materials 16 04988 g008
Figure 9. The transient temperature profile inside the film in the vertical direction of 2-D and 1-D based on the CV model. (a) Kn = 0.5; (b) Kn = 1.
Figure 9. The transient temperature profile inside the film in the vertical direction of 2-D and 1-D based on the CV model. (a) Kn = 0.5; (b) Kn = 1.
Materials 16 04988 g009
Figure 10. The transient temperature profile inside the film in the vertical direction of 2-D and 1-D based on the DPL model. (a) Kn = 0.5; (b) Kn = 1.
Figure 10. The transient temperature profile inside the film in the vertical direction of 2-D and 1-D based on the DPL model. (a) Kn = 0.5; (b) Kn = 1.
Materials 16 04988 g010
Table 1. Physical parameters of nano-silicon thin film and an ultrafast laser.
Table 1. Physical parameters of nano-silicon thin film and an ultrafast laser.
T0 = 300 Kl = 41 nmJ = 732 J/m2
K = 148 W/(m·K)τq = 6.5 ps δ = 15.3 nm
Cv = 2.3 × 106 J/(m3·K)tp = 0.65 psR = 0.93
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mao, Y.; Liu, S.; Liu, J.; Yu, M.; Li, X.; Yang, K. Non-Fourier Heat Conduction of Nano-Films under Ultra-Fast Laser. Materials 2023, 16, 4988. https://doi.org/10.3390/ma16144988

AMA Style

Mao Y, Liu S, Liu J, Yu M, Li X, Yang K. Non-Fourier Heat Conduction of Nano-Films under Ultra-Fast Laser. Materials. 2023; 16(14):4988. https://doi.org/10.3390/ma16144988

Chicago/Turabian Style

Mao, Yudong, Shouyu Liu, Jiying Liu, Mingzhi Yu, Xinwei Li, and Kaimin Yang. 2023. "Non-Fourier Heat Conduction of Nano-Films under Ultra-Fast Laser" Materials 16, no. 14: 4988. https://doi.org/10.3390/ma16144988

APA Style

Mao, Y., Liu, S., Liu, J., Yu, M., Li, X., & Yang, K. (2023). Non-Fourier Heat Conduction of Nano-Films under Ultra-Fast Laser. Materials, 16(14), 4988. https://doi.org/10.3390/ma16144988

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop