1. Introduction
Platinum–rhodium alloy material has high temperature resistance, corrosion resistance and oxidation resistance, which is used to process glass fiber in engineering manufacturing. There are a lot of small openings on the platinum–rhodium glass fiber leaky plate (Pt–Rh GFLP). Under the action of high temperature molten liquid in the low stress oxidation environment for a long time, the temperature rises rapidly, which inevitably causes the oxidation volatilization of Pt–Rh GFLP. The stress concentration and temperature field concentration appear near the opening. The bearing capacity decreases, and the service life shortens. Therefore, the research on the temperature field in the vulnerable area of the Pt–Rh GFLP under the extremely fast heating environment is in urgent need of theoretical research as an engineering reference.
Since Fourier established the mathematical model of heat conduction in the 19th century, Fourier law has been widely used in various fields of heat conduction problem analysis. For the steady-state heat transfer process with long thermal action time, the Fourier law can be used to describe the relationship between heat flow density and temperature gradient, which can meet the accuracy requirements [
1,
2,
3]. However, Fourier’s law implies the assumption that the propagation speed of thermal disturbance is large and does not involve the time term of heat transfer, so it is no longer applicable to the problem of heat conduction in the extremely fast heating environment of Pt–Rh GFLP [
4,
5,
6,
7,
8]. In order to overcome the limitation of Fourier law, Cattaneo and Vernote independently proposed a non-Fourier heat conduction model with heat flow delay phase, taking into account the influence of heat flow change rate on heat conduction [
4]. Sarkar and Haji Sheikh [
5] studied the hyperbolic heat conduction problem of finite plates made of dielectric materials by Laplace transformation. Barletta and Zanchini [
7] analyzed the hyperbolic heat conduction problem of large cylinder with internal heat source and heat convection with external fluid. Tang and Araki [
8] solved the non-Fourier heat conduction problem of a finite medium under periodic surface heating. Moosaie [
9] studied the non-Fourier heat conduction problem of solid-state lasers by numerical simulation. Atefi and Talaee [
10] solved the non-Fourier temperature field of a large cylinder with boundary conditions that did not change with time by using the separation variable method. Mishra and Sahai [
11] use lattice Boltzmann method to study the non-Fourier heat conduction problem of one-dimensional cylinder and sphere. Ma et al. [
12] got a general solution of the scattered field of thermal waves from the cylindrical subsurface inclusion in the large plate. Zhang et al. [
13] studied a thermal shock problem of an elastic strip made of functionally graded materials based on the fractional heat conduction theory. Niwa et al. [
14] studied the high-temperature bending fatigue properties of oxide dispersion-strengthened platinum–rhodium alloy under high axial stress. Sellitto et al. [
15] presented generalized heat-conduction laws from a mesoscopic perspective and studied mesoscopic theories of heat transport in nanosystems. Machrafi [
16] provided for a systematic presentation of extended non-equilibrium thermodynamics in nanosystems with a high degree of applicability. Jou et al. [
17] studied the formulation of extended irreversible thermodynamics within the framework of the rational thermodynamics. Cimmelli [
18] present an overview of the modern approaches to continuum nonequilibrium thermodynamics from the perspective of their connection with the problem of heat conduction with finite speed.
The thermal wave theory (Cattaneo–Vernotte equation) and Boltzmann method were investigated in the above literature [
5,
6,
7,
8,
9,
10,
11]. According to these studies we have the following views on the Cattaneo–Vernotte model (CV model). One question is the compatibility of the CV equation with the second law of thermodynamics. For the CV model allows heat flow from low temperature to high temperature, it is easy to be considered that the CV model violates the second law of thermodynamics. However, in fact the case of deviation from equilibrium is dealt with in the CV model. The entropy of the system depends not only on the properties of equilibrium states such as temperature, but also on some other parameters such as the heat flux. This is one of the basic understandings of the extended thermodynamics. On this basis, it can be proved that the entropy produced by the thermal wave theory such as the CV model is non-negative [
19]. Therefore, the thermal wave model such as the CV model is compatible with the basic principles of the second law of thermodynamics. The essence of this question is that the CV model is mistakenly put under the conceptual framework of the equilibrium continuum model. In fact, the CV model can only be built on the conceptual framework of extended thermodynamics. Besides, the Boltzmann transfer equation (BTE) is compared with the thermal wave theory in reference [
19]. In the thermal wave theory, the nonlocality in time is generally considered and the nonlocality in space is generally ignored. BTE considers both temporal and spatial nonlocality. Considering a gas system with large thermal relaxation time [
19,
20], if only time nonlocality is retained, the system will lose the space transport completely, and it will be only related with time, which is totally inconsistent with the actual situation. This shows that the spatial nonlocality cannot be ignored, and the thermal wave theory cannot be applied in the gas system with large thermal relaxation time. In this paper, the research object is metallic material, and its thermal relaxation time is finite, so the CV model is generally applicable.
In engineering practice, the Pt–Rh GFLP is heated by ultrashort laser pulse. By controlling the frequency which refers to the periodic repetition of laser pulses and amplitude of the ultrashort laser pulse, the temperature field near the opening of the Pt–Rh GFLP can be measured, and the subsurface microstructure information of the platinum rhodium alloy material can be obtained, so as to realize the optimal design of the opening structure of the Pt–Rh GFLP. In this paper, the problem of thermal wave scattering and surface temperature distribution in a Pt–Rh GFLP structure based on the wave equation of non-Fourier heat conduction is studied, and the effects of various physical parameters on temperature distribution is analyzed. The analysis in this paper is of great significance for solving the problem of concentrated temperature distribution near openings or defects in various metal plates at extremely high temperatures.
2. Wave Equation of Heat Conduction and Its General Solution
The research object is a Pt–Rh GFLP with a circular opening and a large length in the
x direction as shown in
Figure 1. A thermal wave process can be formed inside the material as a laser pulse beam with a modulation frequency of
f irradiated on the outer surface of the leaky plate. The boundary conditions of the Pt–Rh GFLP are that the upper, lower, front, and rear boundary plane are all constant temperature surfaces equal to the ambient temperature. The boundary conditions of the left and right ends of the Pt–Rh GFLP are not discussed here for the longitudinal length of the plate is large. Because of the small area of the opening wall, it is very difficult to exchange heat with the outside world, so it is treated as an adiabatic surface in the paper. At the same time, suppose the temperature field distribution in the
z-direction is uniform, which means the problem of two-dimensional temperature field in the coordinate (
x,y) is studied.
The temperature governing equation in the solid medium without internal heat source can be obtained according to the non-Fourier heat conduction equation [
21].
where
is Laplace operator,
;
are respectively the thermal conductivity, specific heat capacity at constant pressure and density;
D =
λ/
ρ cp is the thermal diffusivity;
is the thermal wave propagation speed;
τ is the thermal relaxation time;
T is the temperature in the Pt–Rh GFLP. The density of platinum–rhodium 5% alloy leaky plate and platinum–rhodium 10% alloy leaky plate are taken
ρ = 21 g/cm
3 and
ρ = 20.5 g/cm
3 respectively; The thermal conductivity of platinum–rhodium 5% alloy leaky plate and platinum–rhodium 10% alloy leaky plate measured by Wuxi Yintai Metal Products Co. Ltd are respectively
λ = 69.9 W/(m·K) and
λ = 70.05 W/(m·K). The difference of thermal conductivity between platinum–rhodium 5% and platinum–rhodium 10% is very small, and have little effect on the numerical calculation results, therefore, they are treated as equal in their thermal conductivity in the paper, taken
λ = 70 W/(m·K) here. The specific heat capacity at constant pressure
cp = 0.133 KJ/(Kg·K) The thermal relaxation time
τ = 10
−12 s.
For a wave problem, the analytical solution can be divided into two types: time domain solution and frequency domain solution. And the two types can be converted to each other through Fourier transform method. For the internal flaw detection and detection of materials, it is necessary to use a periodic heat load to implement thermal incentives to make use of the phase and amplitude of the thermal echo, as well as the wave number and spatial attenuation coefficient, to determine the position, orientation and size of the openings in the material.
According to the Fourier decomposition theorem, the periodic heat conduction process can be regarded as the superposition of several simple harmonic motions. Therefore, the periodic non-steady-state heat conduction of the problem can be studied according to the normative analysis format of the wave theory.
Study the periodic non-steady-state heat conduction solution of Equation (1), assuming the expression of the temperature field is
Among them, the temperature amplitude should satisfy the Helmholtz equation of the following form
where Re means taking the real part;
is the average temperature of the environment; Excess temperature amplitude
;
is the Laplace operator,
;
ω is the circular frequency;
is the complex wave number,
;
is the imaginary unit;
α, β are respectively the thermal wave number and absorption coefficient, taking
α > 0,
β > 0, after normalization, that is
where
k is the thermal wave number,
k =
ω/
c when there is no diffusion effect,
µ is the thermal diffusion length.
Using complex variable function method, introducing complex variables
Equation (3) can be transformed into the following form
The general solution of the thermal wave scattering field of an opening in a Pt–Rh GFLP determined by Equation (3) is
where
An is the mode coefficient of the scattered wave generated by the circular opening in the leaky plate, which is determined by the boundary conditions;
is the Hankel function of the first kind.
The conformal mapping method can be used to meet the boundary conditions of the circular openings when solving the wave scattering plane problem of openings of any shape in solids. The mapping function of the non-circular opening boundary on the
-plane to the unit circle boundary on the
-plane can be taken as [
10,
11]
Giving proper mapping function, Equation (10) can be applied with openings of any shape, such as openings of triangular or rectangular cross section, etc. At this time, there may be
on the
-plane. The coordinates of any point can be expressed as
when taking Polar coordinates on the
-plane. In this way, the general solution of the non-circular opening thermal wave scattering field in the
-plane can be described as
For an elliptical opening with a long axis and a half axis of a
r1 and
r2, the conformal mapping function can be taken as
where
.
4. Numerical Examples
The thermal wave scattering and surface temperature distribution of elliptical opening in Pt–Rh GFLP are calculated in order to show the effectiveness of the analysis and calculation method in this paper. The radius of the circular opening is . The amplitude of the pulse mentioned below is 1270 °C.
Figure 2 correspond to the position of the circular opening along the transversal direction of the leaky plate, that is, in the
y-direction: distance from the upper boundary is taken
c1 = 200 mm, distance from the lower boundary is taken
c2 = 200 mm; The laser pulse frequency is taken
f = 80 Hz,
f = 160 Hz,
f = 240 Hz respectively; The mass fraction of rhodium element in Pt–Rh GFLP is 5%.
Figure 2a shows the temperature distribution in the section
y ∈ (−30 mm, 30 mm) with the center of the circular opening as the origin, and
Figure 2b shows the temperature distribution in the section
y ∈ (−1.5 mm, 1.5 mm) with the center of the circular opening as the origin. From the consistent trend of the three temperature distribution curves, it can be found that the temperature peaks at
y = 0 position, and the temperature begins to increase suddenly in the section near
y ∈ (−0.5 mm, 0.5 mm). The temperature gradually decreases as the distance from position increases, and the temperature distribution is symmetrical along the line
y = 0. From the comparison of the differences in the three temperature distribution curves, as the laser pulse frequency increases, the peak temperature gradually decreases, and the length of section where the temperature suddenly increases gradually decreases; In the section other than
y ∈ (−1.5 mm, 1.5 mm), the laser pulse frequency has no obvious influence on the distribution of temperature.
Figure 3 correspond to the position of the circular opening along the transversal direction of the leaky plate, that is, the
y-direction: distance from the upper boundary is taken
c1 = 200mm, distance from the lower boundary is taken
c2 = 200mm; The laser pulse frequency is taken
f = 80 Hz; The mass fraction of rhodium element in Pt–Rh GFLP is taken 5% and 10% respectively.
Figure 3a shows the distribution of temperature in the section
y ∈ (−30 mm, 30 mm) with the center of the circular opening as the origin, and
Figure 3b shows the distribution of temperature in the section
y ∈ (−1.5 mm, 1.5 mm) with the center of the circular opening as the origin. From the consistent change trend of the two temperature distribution curves, it can be found that the temperature peaks at
y = 0 position, and the temperature starts to increase suddenly in the section of
y ∈ (−0.7 mm, 0.7 mm), the temperature gradually decreases with the increase of the distance from
y = 0 position, and the temperature distribution is symmetrical along the line
y = 0. From the comparison of the differences in the three temperature distribution curves, as the laser pulse frequency increases, the peak temperature gradually decreases, and the length of section where the temperature suddenly increases gradually decreases; In the section other than
y ∈ (−1.5 mm, 1.5 mm), the laser pulse frequency has no obvious influence on the distribution of temperature.
Figure 4 correspond to the position of the circular opening along the transversal direction of the leaky plate, that is, the
y-direction: distance from the upper boundary is taken
c1 = 100 mm, distance from the lower boundary is taken
c2 = 100 mm,
c2 = 150 mm,
c2 = 200 mm,
c2 = 250 mm respectively; The laser pulse frequency is taken
f = 80 Hz; The mass fraction of rhodium element in Pt–Rh GFLP is 5%.
Figure 4a shows the distribution of temperature in the section
y ∈ (−2 mm, 2 mm) with the center of the circular opening as the origin, and
Figure 4b shows the distribution of temperature in the section
y ∈ (−
c1,
c2) with the center of the circular opening as the origin. It can be found from
Figure 4b that when
c2 ≥
c1, the larger the
c2, the lower the peak temperature. Combining curve 3, it can be found that the temperature in the section
y ∈ (−100 mm, −50 mm) is increasing, the temperature in the section
y ∈ (−50 mm, −1.5 mm) is decreasing, and the temperature in the section
y ∈ (1.5 mm, 100 mm) is decreasing. This is different from the case where the temperatures of the first two groups are symmetrically distributed along the line
y = 0. This is caused by the different distance of the circular opening from the upper and lower boundaries of the leaky plate. In the state without openings, the line
y = −50mm is the center line of the transversal direction of the leaky plate, that is, the
y-direction, when the leaky plate is heated by the laser pulse, the center line is the farthest from the boundary of the leaky plate, and has the least heat exchange with the outside environment, so its temperature is high; in the state with a circular opening, because the line
y = 0 contains the center of the circular opening, the temperature suddenly increases when heated by the laser pulse, so there is a prominent part of the temperature in the figures. It can be found from
Figure 4a that there is no significant difference in the temperature at the center line of the
y-direction when
c2 is different. According to the rules found in
Figure 4a,b, the opening center can be set reasonably outside the center line of the
y-direction on the Pt–Rh GFLP to avoid overheating at the circular opening of the leaky plate due to the superposition of factors such as difficulty in heat dissipation at the center of the leaky plate and sudden temperature increase at the circular opening of the leaky plate.
Figure 5 corresponds to the position of the circular opening along the transversal direction of the leaky plate, that is, the
y-direction: distance from the upper boundary is taken
c1 = 100 mm,
c1 = 200 mm respectively, distance from the lower boundary is taken
c2 = 100 mm,
c2 = 200 mm respectively; The laser pulse frequency is taken
f = 80 Hz; The mass fraction of rhodium element in Pt–Rh GFLP is from 0 to 30%. According to the data measured by Wuxi Yintai Metal Products Co. Ltd, during the process of the mass fraction of rhodium in the Pt–Rh GFLP changing from 0% to 30%, the changes of specific heat capacity and thermal conductivity are very small, and have little effect on the numerical calculation results, so they are treated as equal in the paper. The density of platinum–rhodium
n% is 21.45 g/cm
3 × (100 −
n)% + 12.41 g/cm
3 ×
n%.
Figure 5 shows the temperature at the center of the circular opening, that is, the position of
y = 0. It can be found from the
Figure 5 that as the mass fraction of rhodium element in Pt–Rh GFLP increases, the temperature at the center of the circular opening gradually increases. The temperature at the circular opening is higher when
c1 =
c2, and in the case of
c1 =
c2, the temperature is always equal. The temperature is lower when
c1 ≠
c2.
Figure 6 corresponds to the position of the circular opening along the transversal direction of the leaky plate, that is, the
y-direction: distance from the upper boundary is taken
c1 = 100 mm,
c1 = 200 mm respectively, distance from the lower boundary is taken
c2 = 100 mm,
c2 = 200 mm respectively; The mass fraction of rhodium element in Pt–Rh GFLP in
Figure 5 is 5%; The pulse frequency is taken from
f = 0.1 Hz to
f = 300 Hz.
Figure 6 shows the temperature at the center of the circular opening, that is, the position of
y = 0. It can be found from
Figure 6 that when the frequency is in the section (0 Hz, 5 Hz), the temperature at the center of the circular opening increases with increasing frequency, and in the section (5 Hz, 300 Hz), the temperature at the center of the circular opening gradually decreases and tends to be stable with the increase of frequency. The temperature at the circular opening is higher when
c1 =
c2, and in the case of
c1 =
c2, the temperature is always equal. The temperature is lower when
c1 ≠
c2. We can see that in
Figure 6 with a peak temperature value at about 5 Hz excitation. A reasonable explanation is that thermal wave behavior is similar to elastic waves in solid mechanics or acoustic waves in physics. The wave equation or vibration equation in dynamics is a hyperbolic equation. When the external excitation is close to the natural frequency of the system, dynamic stress concentration or resonance will occur with well-pronounced amplitude at a certain frequency. The heat conduction equation in this paper is also a hyperbolic equation, so it is reasonable to produce well-pronounced temperature excess amplitude within a certain thermal excitation frequency range. It provides an important theoretical reference for us to control the maximum temperature of the temperature field by adjusting the thermal excitation frequency. When heating with laser pulses, we can avoid the section around 5 Hz to avoid excessive temperature.