Next Article in Journal
3D Sugar Printing of Networks Mimicking the Vasculature
Next Article in Special Issue
A CMOS-Thyristor Based Temperature Sensor with +0.37 °C/−0.32 °C Inaccuracy
Previous Article in Journal
Design Optimization of a Compact Double-Ended-Tuning-Fork-Based Resonant Accelerometer for Smart Spindle Applications
Previous Article in Special Issue
Adaptive Nonlinearity Compensation System for Integrated Temperature and Moisture Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Error Correction of Filament Thermocouples with Different Structures of Junction based on Inverse Filtering Method

1
School of Instrument and Electronics, North University of China, Taiyuan 030051, China
2
Key Laboratory of Instrumentation Science & Dynamic Measurement (North University of China), Ministry of Education, Taiyuan 030051, China
*
Author to whom correspondence should be addressed.
Submission received: 5 November 2019 / Revised: 23 December 2019 / Accepted: 25 December 2019 / Published: 30 December 2019

Abstract

:
Since filament thermocouple is limited by its junction structure and dynamic characteristics, the actual heat conduction process cannot be reproduced during the transient thermal shock. In order to solve this problem, we established a thermocouple dynamic calibration system with laser pulse as excitation source to transform the problem of the restoring excitation source acting on the surface temperature of thermocouple junction into the problem of solving the one-dimensional (1D) inverse heat conduction process, proposed a two-layer domain filtering kernel regularization method for double conductors of thermocouple, analyzed the factors causing unstable two-layer domain solution, and solved the regular solution of two-layer domain by the filtering kernel regularization strategy. By laser narrow pulse calibration experiment, we obtained experimental samples of filament thermocouples with two kinds of junction structures, butt-welded and ball-welded; established error estimation criterion; and obtained the optimal filtering kernel parameters by the proposed regularization strategy, respectively. The regular solutions solved for different thermocouples were very close to the exact solution under the optimal strategy, indicating that the proposed regularization method can effectively approach the actual surface temperature of the thermocouple junction.

1. Introduction

Filament thermocouple has as characteristics simple structure, fast dynamic response, wide range of temperature measurement, ability to measure the ‘single point’ temperature, and its electrical signal is easy to access to the test system, so it has been widely applied in industrial production and scientific research, such as wind tunnel heat flow experiment and explosive thermal shock experiment. Therefore, it is very important to accurately obtain the temperature of the heat sources during the actual tests [1,2,3,4]. Whether there is a wind tunnel environment or an explosion field one, the heat source is usually a directional narrow pulse signal. Its rising time can reach the level of microseconds, and pulse width is in milliseconds. The response of thermocouple to it is impulse response and limited by the dynamic characteristics of the thermocouple itself, such as the thermal inertia of temperature sensing element. Its dynamic response process cannot accurately recover the thermal shock process, which generates great dynamic error [5]. Therefore, it is necessary to analyze whether thermocouple’s dynamic characteristics are competent for the transient temperature test.
The temperature measurement accuracy of the thermocouple is directly affected by the junction welding quality [6]. For filament thermocouple to reduce the heat conduction error and dynamic response error, the structures of junction welded are different. The frequently-used structures are butt-welded and ball-welded as shown in Figure 1a,b.
We can establish a dynamic calibration system with laser pulse as excitation source [7,8,9], excite the thermocouple with a narrow pulse signal, and establish an inverse mathematical model of thermocouple by means of the inter-conversion for the input/output relationship of the thermocouple. Many scholars describe this model as a transfer function, and nonlinear system identification is often used to build a mathematical model to correct the dynamic error with the output of this model and recover the temperature signal at the thermocouple junction. The problem lies in neglecting the physical process of thermocouple heat conduction; thus, the models built are unpractical for engineering applications with large temperature range [10,11,12]. Therefore, from the perspective of inverse filtering correction, the problem of inverse mathematical model identification for the thermocouple can be transformed into the problem of the solution of the process of inverse heat conduction [13].
In inverse heat conduction, our purpose is to correct the impulse response error of thermocouple to approach the actual surface temperature of thermocouple junction. The thermocouple is composed of two metal conductors made of different materials; thus, the two conductors produce two thermal properties during heat conduction. According to the working principle of thermocouple, the thermoelectric potential of the thermocouple is determined by the difference in the temperature distribution function between the two ends of the thermocouple. Therefore, we define the temperature distribution functions at both ends of the thermocouple. The measured output thermoelectric potential difference of thermocouple during thermal shock can be used to solve the surface temperature of thermocouple junction. Therefore, this process can be called the non-characteristic Cauchy problem during heat conduction [14,15]. Any tiny measurement error in the non-characteristic Cauchy problem of heat equation may cause a large error of solution. The continuous dependence of Cauchy numerical solutions were obtained under additional conditions. This is known as conditional stability [16]. In addition, the stable solution of ill-posed problem can be theoretically reconstructed under priori conditions. Due to the severity of this problem, no satisfactory stable approach was not obtained for the Cauchy problem of most classical heat conduction equations, then some regularization strategies are needed [17]. For the 1D inverse heat conduction in monolayer domain, theoretical study and computational implementation have been well studied [18,19], but for thermocouple structure, the problem is transformed into complex inverse heat conduction in two-layer domain. We can decompose multilayer inverse heat conduction into multiple monolayer problems and solve each monolayer problem layer by layer by using the relationship between the layers.
The filament thermocouple was used as a two-layer volume; the junction was defined as the interface between two layers; the heat flow direction is perpendicular to the thermocouple junction; the two metal conductors were defined as of equal length, and one of them is as the free end that its boundary temperature is 0 °C. See Figure 2.
Here, l is the length of the two conductors of the thermocouple (the free end is the first layer; then, the initial value of the working end is l ); R 1 and R 2 are the regions of the free end and the working end, respectively.
Let us say that k 1 , k 2 > 0 are the thermal conductivities of the thermocouple’s two conductors, and α 1 , α 2 > 0 are the thermal diffusion coefficients of the first layer and the second layer, respectively. The temperature distributions in the first and second layers are expressed by u 1 ( x , t ) and u 2 ( x , t ) , t is the time of the heat transfer process. The temperature distributions in R 1 and R 2 are expressed by the following partial differential equations:
u 1 t α 1 2 u 1 x 2 = 0 , 0 < x < l ,   t > 0
u 2 t α 2 2 u 2 x 2 = 0 , l < x < 2 l ,   t > 0
The additional initial value and the boundary conditions are as follows:
u 1 ( x , 0 ) = u 2 ( x , 0 ) = 0 ,   0 < x < 2 l
u 2 ( 2 l , t ) = g ( t ) ,   t > 0
u 2 x ( 2 l , t ) = 0 ,   t > 0
u 1 ( l , t ) = u 2 ( l , t ) ,   t > 0
k 1 u 1 x ( l , t ) = k 2 u 2 x ( l , t ) ,   t > 0
g ( t ) in Equation (4) is boundary temperature of the working end. When boundary temperature of the free end is 0 °C, we can use g ( t ) as junction temperature of the thermocouple.
For similar inverse heat conduction, Shcheglov analyzed the convergence of this problem by hyperbolic perturbation method [20]. For 1D inverse heat conduction in Cartesian systems with two-layer domains, some numerical methods have been proposed [21] but have not been verified in engineering application examples. The purpose of this paper is to obtain the two-layer domain analytic solution for the 1D inverse heat conduction based on the thermocouple during the laser pulse by Fourier transform (FT) [22]. For the serious morbidity of this problem, a regularization strategy was proposed to reconstruct the solution of this problem stably. In order to explain this, we proposed a two-layer domain regularization method based on the analytic solution for the inverse heat conduction. In this method, we used a filtering kernel method to reconstruct stable monolayer solution and provided the selection strategy for optimal regularization parameters. This method is simple but effective.
This paper is composed of the following parts. In the second part, we analyzed the principle and ill-posed inverse heat conduction based on the working principle of thermocouple. In the third part, we used the filtering kernel method to reconstruct the stable solution for the inverse problem under error estimation and provided the convergence estimate. In the fourth part, we reproduced the laser pulse calibration experiment and verified the effectiveness of the proposed method by using experimental samples.

2. Principle Analysis of Inverse Problem

According to the working principle of the thermocouple, we hold that, after the thermocouple’s cold end is compensated, the boundary temperature at the working end g ( t ) in Equation (3), Section 1, is obtained by the actual output thermoelectric potential difference of the thermocouple, and that it is the output sample of the thermocouple in response to the laser pulse.
Assume that g   ϵ   L 2 ( 0 , ) , u 1 ( x , t ) and u 2 ( x , t ) are the solutions of the positive problem; then, in the given x   ϵ   [ 0 , 2 l ] , u 1 ( x , · ) , u 2 ( x , · ) , and their partial differentials u 1 ( x , · ) / x , u 2 ( x , · ) / x also belong to space L 2 ( 0 , ) . The solution of the inverse problem is to solve u 1 ( x , t )   ( 0 x l ) in L 2 ( 0 , ) by using these two conditions: the known g ( t ) and adiabatic boundary x = 2 l .
In a dynamic calibration experiment, there is some error in the test results; then, the actual measured value is defined as g δ ( t )   ϵ   L 2 ( 0 , ) and meets.
g δ g L 2 δ
δ represents measurement error bounds here. Since Cauchy problem is a typical ill-posed problem, any tiny error in experimental sample g δ ( t ) will cause the blow-up of the solution u 1 ( x , t ) within the interval 0 x l . The root of this problem lies in the kernel function in the high frequency part. The ill-posed inverse problem will be analyzed below:
· is expressed as the norm of L 2 ( 0 , ) , we have
f = ( 0 | f ( t ) | 2 d t ) 1 2
Then, we take the FT of f ( t )   ϵ   L 2 ( 0 , ) as
f ( ξ ) = 1 2 π 0 f ( t ) e i ξ t d t ,   ξ   ϵ   L 2 ( 0 , )
Its inverse transformation is
f ( t ) = 1 2 π 0 f ( ξ ) e i ξ t d ξ
We first conduct FT on t of the Equations (2)–(6) in the inverse problem in R 2 and can obtain the following equation:
u 2 ( x , ξ ) = c o s h ( ( 2 l x ) i ξ α 2 ) g ( ξ ) ,   l x 2 l
where
i ξ α 2 = | ξ / ( 2 α 2 ) | 1 2 ( 1 + i s g n ( ξ ) )
where s g n is a sign function.
Then, the first order partial derivative with respect to x is conducted on Equation (12):
u 2 ( x , ξ ) x = i ξ α 2 s i n h ( i ξ α 2 ( 2 l x ) ) g ( ξ )
FT is conducted with respect to ξ in Equations (12) and (14):
u 2 ( x , t ) = 1 2 π 0 g ( ξ ) c o s h ( ( 2 l x ) i ξ α 2 ) e i ξ t d ξ
u 2 ( x , ξ ) x = 1 2 π 0 g ( ξ ) i ξ α 2 × s i n h ( i ξ α 2 ( 2 l x ) ) e i ξ t d ξ
Next, FT is conducted with respect to t in Equations (1), (3), (6) and (7) in R 1 :
u 1 ( x , ξ ) = c o s h ( i ξ α 1 ( l x ) ) c o s h ( i ξ α 2 l ) g ( ξ ) + k 2 α 1 k 1 α 2 s i n h ( i ξ α 1 ( l x ) ) × s i n h ( i ξ α 2 l ) g ( ξ ) , 0 x l
Then, the first order partial derivatives with respect to x is conducted on Equation (17):
u 1 ( x , ξ ) x = i ξ α 1 s i n h ( i ξ α 1 ( l x ) ) c o s h ( i ξ α 2 l ) g ( ξ ) k 2 α 1 k 1 α 2 i ξ α 1 c o s h ( i ξ α 1 ( l x ) ) × s i n h ( i ξ α 2 l ) g ( ξ ) ,   0 x l
Then, FT is conducted with respect to ξ in Equations (17) and (18) to obtain u 1 ( x , t ) . That is, the numerical solution of the inverse problem is obtained. However, it can be seen from Equations (13) and (18) in R 1 that the reason for its ill-posedness is the existence of unbounded kernel c o s h ( ( l x ) i ξ α 1 ) , which will be amplified infinitely with noise as a factor, thus causing the blow-up of the solution. Therefore, we use the idea of a corrected kernel [23] and provide a bounded kernel to approach unbounded kernel in a new form of filtering kernel. This problem can be solved in this way.

3. Filtering Kernel Method

In R 1 , if ω = i ξ / α 1 , for the unbounded kernel c o s h ( ( l x ) i ξ / α 1 ) , we provide the filtering form as follows:
u 1 β δ ( x , ξ ) = c o s h ( ω ( l x ) ) 1 + β | cosh ( ω ) | m g δ ( ξ )
If β is a minimal value, c o s h ( ω ( l x ) ) 1 + β | cosh ( ω ) | m c o s h ( ( l x ) i ξ α 1 ) , if β > 0 , then c o s h ( ω ( l x ) ) / ( 1 + β | cosh ( ω ) | m ) is bounded.
We provide an important lemma:
Lemma 1.
If a 0 , b 0 , x 0 , ξ   ϵ   ( 0 , ) , then
| c o s h ( a + i b s g n ( ξ ) ) | 1 2 e π 2 2 e a
| cosh ( x ( a + i b s g n ( ξ ) ) ) | e x a
The following two theorems indicated that the regular solution is a good approach to the exact solution.
Theorem 1.
Assume that u 1 ( x , t ) is the regular solution of the inverse problem for the exact solution g ( t ) and u 1 β δ ( x , t ) is that for disturbance data g δ ( t ) ; the disturbance data g δ ( t ) meets g δ g L 2 ( R ) δ , and the boundary temperature u 1 ( l , · ) meets
u 1 ( l , · ) E
Here, E > 0 is a given constant; if
β = ( δ E ) m
then the following convergence estimation within the length interval 0 < x < l can be obtained
u 1 β δ ( x , · ) u 1 ( x , · ) K δ l x E x
where K is a positive constant.
It is worth noting that in Theorem 1, the convergence estimation is only limited within the interval 0 < x < l . However, any useful information on the continuity of the solution cannot be obtained at the boundary x = l . Then, in order to obtain the continuous dependence of the solution at x = l , a stronger transcendental assumption is introduced:
u 1 ( l , · ) P E
where · P represents a norm in Sobolev space H P ( R ) .
Theorem 2.
Assume that u 1 ( x , t ) is the regular solution of the inverse problem for the exact solution g ( t ) , and u 1 β δ ( x , t ) is that for disturbance data g δ ( t ) ; the disturbance data g δ ( t ) meets g δ g L 2 ( R ) δ , and the boundary temperature u 1 ( l , · ) meets Equation (25); if
β = ( δ E ) m 2
then the following convergence estimation can be obtained at x = l
u 1 β δ ( x , · ) u 1 ( x , · ) c 2 δ 1 2 E 1 2 + K 1 E m a x { ( 1 6 m ln E δ ) P , ( δ E ) 3 m 1 6 }
where K 1 is a positive constant.
Then, Theorem 2 is verified. It is divided into two parts. We can first obtain the following equation with Parseval’s equation and triangle inequality:
u 1 β ( x , · ) u 1 β δ ( x , · ) = u 1 β ( x , · ) u 1 β δ ( x , · ) = c o s h ( ω ) 1 + β | cosh ( ω ) | m ( g δ g )     δ ξ R s u p A
If
a = b = i ξ α 1 ,   σ = s g n ( ξ )
it can be seen that
A = c o s h ( ω ) 1 + β | cosh ( ω ) | m
It can be concluded from Lemma 1 that
A e a 1 + β ( 1 2 e π 2 4 ) m 2 e m a
If f ( a ) = e x a 1 + β ( 1 2 e π 2 4 ) m 2 e m a and c 1 = ( 1 2 e π 2 4 ) m 2 , it can be obtained by derivation that
f ( a ) = f ( a ) l x β c 1 e m a ( m l + x ) 1 + β c 1 e m a
It is easy to find that
a 1 = 1 m ln l x β c 1 ( m l + x )
meets f ( a ) = 0 , besides, when a > a 1 , f ( a ) < 0 ; when a < a 1 , f ( a ) > 0 . Therefore, maximum point of f ( a ) is a 1 , and
f ( a ) f ( a 1 ) = c 2 β l x m
where c 2 = ( l x ) l x m m c 1 l x m ( m l + x ) 1 l x m , thus,
A 1 e ( l x ) a 1 + β c 1 e m a c 2 β x l m
It can be verified in combination with (28) and (32) that
u 1 β ( x , · ) u 1 β δ ( x , · ) c 2 β x l m
Besides,
u 1 β ( x , · ) u 1 ( x , · ) = u 1 β ( x , · ) u 1 ( x , · ) = ( c o s h ( ω ) 1 + β | cosh ( ω ) | m c o s h ( ω ) ) g β | c o s h ( ω ) | m 1 + β | c o s h ( ω ) | m ( 1 + | ξ α 1 | 2 ) p 2 c o s h ( ω ) g E ξ R s u p B
where B = β | cosh ( ω ) | m 1 + β | cosh ( ω ) | m ( 1 + | ξ α 1 | 2 ) p 2
It can be concluded from Lemma 1 that
B 1 + β ( 1 2 e π 2 4 ) m 2 e m a ( 1 + | ξ α 1 | 2 ) p 2
If c 1 = ( 1 2 e π 2 4 ) m 2 , the estimation is conducted under two conditions:
(1) If a a 2 = 1 m ln ( β 1 3 m ) ,
β e m a 1 + β c 1 e m a ( 1 + | ξ α 1 | 2 ) p 2 β e m a 2 β 1 1 3 m
(2) If a > a 2 ,
B β e m a 1 + β c 1 e m a ( 1 + | ξ α 1 | 2 ) p 2   1 c 1 ( 1 + | ξ α 1 | 2 ) p 2 1 c 1 a 2 p = 1 c 1 ( 1 m ln ( β 1 3 m ) ) p
Therefore, it can be concluded in combination with (34) and (35) that
u 1 β δ ( x , · ) u 1 ( x , · ) u 1 β δ ( x , · ) u 1 β ( x , · ) + u 1 β ( x , · ) u 1 ( x , · ) c 2 β 1 m δ + K 1 E m a x { ( ( 1 m ln ( β 1 3 m ) ) p , β 1 1 3 m   ) } = c 2 δ 1 2 E 1 2 + K 1 E m a x { ( ( 1 6 m ln E δ ) p , ( δ E ) 3 m 1 6   ) }
where K 1 is a positive constant.
We theoretically proved that on the premise of appropriate selection rules for regularization parameters, filtering kernel method is effective. The regular solution for the inverse problem will be obtained by combining experimental samples with the proposed regularization strategy of inverse problem.

4. Calibration System and Experimental Analysis

Although the widely used excitation signal for the dynamic calibration of temperature sensor is temperature step signal, it is difficult to generate ideal step temperature signal, which has brought errors to the dynamic calibration of temperature sensor with step excitation signal, but semiconductor laser can generate very narrow pulse excitation signal. Narrow pulse high-temperature excitation signal is easy to produce and has broader frequency components than step signal; then, it can fully cover the spectrum of thermocouple sensor and excite all of its dynamic characteristics [24]. The RFL-A500D high-power semiconductor laser manufactured by Wuhan Raycus in China is used as the excitation source. Its central wavelength is 915 nm, and its output power is 515 W and is adjustable. The laser pulse calibration system is constructed in this way. See Figure 3 for the composition of the calibration system.
Figure 4 shows our experimental environment; the thermocouple calibrated is a K -type filament thermocouple with a junction diameter of 1 mm; its measuring range is 0–1200 °C. The conductors are made of nickel and chromium, (nickel: thermal conductivity is 90.0 W/m·k, thermal capacity is 0.46 kJ/kg·°C and density is 8.9 g/cm3; chromium: thermal conductivity is 93.7 W/m·k, thermal capacity is 0.45 kJ/kg·°C and density is 7.19 g/cm3). The filament thermocouple is fixed on the bracket. Uniform laser beams [25,26,27] are formed after laser pulses converge in the laser flux uniformity system to heat the surface of thermocouple junction. The fiber optic probe of high speed infrared thermometer (the fast response time is 6 (μs)) is aligned to the surface of thermocouple junction. The measured temperature signal is used as a response signal g ( t ) that acts as a laser pulse on the junction surface. The thermocouple potential difference signal is amplified by an amplifier; and the thermocouple laser pulse response signal g ( t ) and infrared thermometer output signal are synchronously collected by the data acquisition and processing system. Besides, the digital cold end compensation processing is conducted on the free end of the thermocouple, with the output signal of the infrared thermometer as the exact solution u 1 ( l , t ) .
The output power of the laser is adjusted by percentage; thus, we set the output power of the laser to 90%. The theoretical maximum temperature at 90% power should not be higher than 1000 °C. The pulse width is set to 5 ms. The sampling frequency of data acquisition and processing system is 1 MHz. Both butt-welded and ball-welded thermocouples are excited at an ambient temperature of 21 °C. The pulse response data and junction surface temperature data of the thermocouples at 90% powers are obtained. See Figure 5 (the following operation results and charts are obtained in MATLAB).
It can be seen that peak value of infrared temperature is 694.8 °C and the rise time is 12.5 ms. Different response times of the two junction structures are obtained and both are slower than infrared signal, and there is large deviation in the peak value of temperature.
The following formula is used to generate error data.
g δ ( t ) = g ( t ) + σ · r a n d n ( s i z e ( g ( t ) ) )
where t is sample time, determined by sample size and sampling rate, t i = ( i 1 ) Δ t , Δ t = 1 / ( n 1 ) , i = 1 , 2 , , n . According to Equation (8) in Section 2, the value of δ at the junction is:
δ : = g δ ( t ) g ( t ) L 2 = 1 n i = 0 n ( g i δ g i ) 2
On the basis of Equations (26), (27), (39) and (40) in Section 3, β = ( δ E ) m 2 is used to take different m values under different σ conditions to obtain a series of regular solutions u 1 β δ ; the value of E is the maximum range of thermocouple; σ is respectively 0.01, 0.001, 0.0001, and m is integer from 0 to 20. We first observed the changes in δ as shown in Figure 6a,b.
It can be seen from Figure 6 that the two structures have the same tendency, and the smaller the σ value, the more obvious the attenuation. However, when σ values are different, the smaller the σ value, the faster the convergence. When m 8 , σ = 0.0001 , the error attenuation tends to be stable first, which indicates that the regularization parameters selected tend to be optimal.
When σ = 0.0001 , m is 2, 5 and 8, respectively; the regular solution and exact solution under 90% power output are as shown in Figure 7a–d.
The calculation results under different conditions are acquired and listed in Table 1 and Table 2.
It can be seen from Figure 7 and Table 1 and Table 2 that when the value of m is small, although the regular solutions of two junction structures have some distortion, both of them can still approach the exact solution. Furthermore, when the regularization parameters reach the optimal value, the regular solution and the exact solution are very close, we calculated the relative error of butt-welded signal between the regular solution and the exact solution in the pulse region when σ = 0.0001 and m = 8 . See Figure 8.
The maximum relative error is 7.64%; the relative error of peak value is even less; thus, the proposed regularization method is effective for recovering the junction surface temperature of filament thermocouple.

5. Conclusions

Since the thermocouple’s dynamic characteristics are affected by its junction structure and cannot faithfully recover the actual temperature of the measured medium, we constructed a laser pulse calibration system to quantitatively analyze its dynamic characteristics and describe the problem as the solution for 1D inverse heat conduction from the point of view of inverse filtering. For thermocouple of a two-layer conductor structure, this problem cannot be solved by most regularization methods in monolayer domain. We can use inverse heat conduction in stable monolayer domain to extend and solve the two-layer domain problem. We find the kernel function, which causes the instability of the numerical solution in monolayer domain, and use the stable kernel function to filter it to provide the error estimation from the two-layer domain filtering kernel method. With the experimental data from laser pulse calibration, we obtain the optimal regularization parameter, verify that the proposed regularization method can effectively recover the surface temperature of the thermocouple, and make up the deficiency for the dynamic characteristics of the thermocouple methodologically.

Author Contributions

Conceptualization, C.Z. and Z.Z.; methodology, C.Z. and Z.Z.; software, C.Z.; validation, C.Z. and Z.Z.; formal analysis, C.Z. and Z.Z.; investigation, C.Z. and Z.Z.; resources, Z.Z.; data curation, C.Z. and Z.Z.; writing—original draft preparation, C.Z. and Z.Z.; writing—review and editing, C.Z.; visualization, C.Z. and Z.Z.; supervision, Z.Z.; project administration, Z.Z.; funding acquisition, Z.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

The paper is sponsored by the Fund for Shanxi ‘1331 Project’ Key Subject Construction of China and the National Natural Science Foundation of China (Approved No. 51575499).

Acknowledgments

We would like to express our gratitude to all those Fuwen Lai and Wuliang Yin who have offered us professional instructions in this study and manuscript writing, and we gratefully acknowledge Yanfeng Li who have done a lot of experimental work early.

Conflicts of Interest

The authors declare no conflict of interest and the funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Michalski, L.; Eckersdorf, K.; McGhee, J. Temperature Measurement; Wiley: Chichester, UK, 1991. [Google Scholar]
  2. Wiśniewski, S. Temperature Measurement in Engines and Thermal Facilities; WNT: Warszawa, Poland, 1983. (In Polish) [Google Scholar]
  3. Agarwal, S.; Sahoo, N.; Singh, R.K. Experimental techniques for thermal product determination of coaxial surface junction thermocouples during short duration transient measurements. Int. J. Heat Mass Transf. 2016, 103, 327–335. [Google Scholar] [CrossRef]
  4. Prajapati, H.; Deshmukh, N.N. Design and development of thin wire sensor for transient temperature measurement. Measurement 2019, 140, 582–589. [Google Scholar] [CrossRef]
  5. Magdalena, J. Reduction of dynamic error in measurements of transient fluid temperature. Arch. Thermodyn. 2011, 32, 55–66. [Google Scholar]
  6. Li, J.P.; Chen, H.; Zhang, S.Z.; Yu, H.R. On the response of coaxial surface thermocouples for transient aerodynamic heating measurements. Exp. Therm. Fluid Sci. 2017, 86, 141–148. [Google Scholar] [CrossRef]
  7. Castellini, P.; Rossi, G.L. Dynamic characterization of temperature sensors by laser excitation. Rev. Sci. Instrum. 1996, 67, 2595–2601. [Google Scholar] [CrossRef]
  8. Li, Y.F.; Zhang, Z.J.; Hao, X.J.; Chen, H. Dynamic calibration method of temperature sensor based on Quasi-delta pulse temperature excitation. Microw. Opt. Technol. Lett. 2018, 60, 212–219. [Google Scholar] [CrossRef]
  9. Li, Y.; Zhang, Z.; Cheng, H.; Juan, W. Research on Dynamic Compensation Technology of Thermocouple Sensor. In Proceedings of the International Seminar on Measurement, Hangzhou, China, 4–5 September 2017; pp. 41–45. [Google Scholar]
  10. Yang, Z.X.; Meng, X.F. Research on the Dynamic calibration of thermocouple and temperature excitation signal generation method based on shock-tube theory. J. Eng. Gas. Turbines Power 2014, 136, 071602. [Google Scholar] [CrossRef]
  11. Yang, Z.X.; Meng, X.F. Study on transient temperature generator and dynamic compensation technology. Appl. Mech. Mater. 2014, 511, 161–164. [Google Scholar] [CrossRef]
  12. Zimmerschied, R.; Isermann, R. Nonlinear time constant estimation and dynamic compensation of temperature sensors. Control Eng. Pract. 2010, 18, 300–310. [Google Scholar] [CrossRef]
  13. Beck, J.V. Nonlinear estimation applied to the nonlinear inverse heat conduction problem. Int. J. Heat Mass Transf. 1970, 13, 703–716. [Google Scholar] [CrossRef]
  14. Eldén, L. Approximations for a Cauchy problem for the heat equation. Inverse Prob. 1987, 3, 263–273. [Google Scholar] [CrossRef]
  15. Eldén, L. Hyperbolic approximations for a Cauchy problem for the heat equation. Inverse Prob. 1988, 4, 59–71. [Google Scholar] [CrossRef]
  16. Tikhonov, A.N.; Arsenin, V.Y. Solution of Ill-Posed Problems; Winston/Wiley: New York, NY, USA, 1977. [Google Scholar]
  17. Engl, H.W.; Hanke, M.; Neubauer, A. Regularization of Inverse Problems; Kluwer Academic Publisher: Dordrecht, The Netherlands; Boston, MA, USA; London, UK, 1996. [Google Scholar]
  18. Wei, T.; Yamamoto, M. Reconstruction of a moving boundary from Cauchy data in one-dimensional heat equation. Inverse Prob. Sci. Eng. 2009, 17, 551–567. [Google Scholar] [CrossRef]
  19. Tanana, V.P.; Sidikova, A.I. Approximate solution of an inverse boundary value problem for a system of differential equations of parabolic type and estimation of the error of this solution. Trudy Inst. Mat. Mekhaniki UrO RAN 2019, 25, 247–264. [Google Scholar]
  20. Shcheglov, A.I. A method for the approximate solution of an inverse problem for the heat equation. Zhurnal Vychislitelnoi Mat. Mat. Fiz. 1992, 32, 904–916. [Google Scholar]
  21. Aviles-Ramos, C.; Haji-Sheikh, A.; Beck, J.V. Exact solution of heat conduction in composite materials and application to inverse problems. ASME Heat Transf. 1997, 351, 305–314. [Google Scholar] [CrossRef]
  22. Sidikova, A.I. The study of an inverse boundary problem for the heat conduction equation. Sib. J. Numer. Math. 2019, 22, 81–98. [Google Scholar] [CrossRef]
  23. Qian, Z.; Fu, C.L. Regularization strategies for a two-dimensional inverse heat conduction problem. Inverse Probl. 2007, 23, 1053–1068. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Zu, J.; Zhang, H.Y. Dynamic calibration method of high-pressure transducer based on quasi-δ function excitation source. Measurement 2012, 45, 1981–1988. [Google Scholar] [CrossRef]
  25. Cheng, H.; Zhang, Z.; Li, Y. Study on laser beam homogenization in temperature sensor calibration. Laser Mag. 2017, 38, 44–47. (In Chinese) [Google Scholar]
  26. Manjhi, S.K.; Kumar, R. Performance assessment of K-type, E-type and J-type coaxial thermocouples on the solar light beam for short duration transient measurements. Meas. J. Int. Meas. Confed. 2019, 146, 343–355. [Google Scholar] [CrossRef]
  27. Kim, S.J.; Park, G.C. Development of algorithm to measure temperatures of liquid/gas phases using micro-thermocouple and experiment with optical chopper. J. Mech. Sci. Technol. 2007, 21, 184–195. [Google Scholar] [CrossRef]
Figure 1. Two types of filamentous thermocouple: (a) butt-welded and (b) ball-welded.
Figure 1. Two types of filamentous thermocouple: (a) butt-welded and (b) ball-welded.
Micromachines 11 00044 g001
Figure 2. Filament thermocouple.
Figure 2. Filament thermocouple.
Micromachines 11 00044 g002
Figure 3. Transient temperature calibration system.
Figure 3. Transient temperature calibration system.
Micromachines 11 00044 g003
Figure 4. Experimental environment.
Figure 4. Experimental environment.
Micromachines 11 00044 g004
Figure 5. Output signal with two junction structures at 90% power.
Figure 5. Output signal with two junction structures at 90% power.
Micromachines 11 00044 g005
Figure 6. Attenuation of δ with the gradual increase of m when σ is 0.01, 0.001 and 0.0001: (a) butt-welded signal and (b) ball-welded signal.
Figure 6. Attenuation of δ with the gradual increase of m when σ is 0.01, 0.001 and 0.0001: (a) butt-welded signal and (b) ball-welded signal.
Micromachines 11 00044 g006
Figure 7. Calculation results of two junction structures: (a) butt-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 2 ; (b) butt-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 5 ; (c) butt-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 8 ; (d) ball-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 8 .
Figure 7. Calculation results of two junction structures: (a) butt-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 2 ; (b) butt-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 5 ; (c) butt-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 8 ; (d) ball-welded signal, exact solution and regular solution when σ = 0.0001 ,   m = 8 .
Micromachines 11 00044 g007
Figure 8. Relative errors between regular solution and exact solution.
Figure 8. Relative errors between regular solution and exact solution.
Micromachines 11 00044 g008
Table 1. σ = 0.0001 and m = 2, 5, 8; the peak value and rise time of butt-welded thermocouple after regularization.
Table 1. σ = 0.0001 and m = 2, 5, 8; the peak value and rise time of butt-welded thermocouple after regularization.
Title 1OriginalM = 2M = 5M = 8
Peak value (°C)262.1621.6674.5692.5
Rise time (ms)25.913.712.712.5
Table 2. σ = 0.0001 and m = 8; the peak value and rise time of butt-welded thermocouple after regularization.
Table 2. σ = 0.0001 and m = 8; the peak value and rise time of butt-welded thermocouple after regularization.
Title 1OriginalM = 8
Peak value (°C)304.4691.2
Rise time (ms)43.212.5

Share and Cite

MDPI and ACS Style

Zhao, C.; Zhang, Z. Dynamic Error Correction of Filament Thermocouples with Different Structures of Junction based on Inverse Filtering Method. Micromachines 2020, 11, 44. https://0-doi-org.brum.beds.ac.uk/10.3390/mi11010044

AMA Style

Zhao C, Zhang Z. Dynamic Error Correction of Filament Thermocouples with Different Structures of Junction based on Inverse Filtering Method. Micromachines. 2020; 11(1):44. https://0-doi-org.brum.beds.ac.uk/10.3390/mi11010044

Chicago/Turabian Style

Zhao, Chenyang, and Zhijie Zhang. 2020. "Dynamic Error Correction of Filament Thermocouples with Different Structures of Junction based on Inverse Filtering Method" Micromachines 11, no. 1: 44. https://0-doi-org.brum.beds.ac.uk/10.3390/mi11010044

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