Application of Tucker Decomposition in Temperature Distribution Reconstruction

: Constrained by cost, measuring conditions and excessive calculation, it is difﬁcult to reconstruct a 3D real-time temperature ﬁeld. For the purpose of solving these problems, a three-dimensional temperature distribution reconstruction algorithm based on Tucker decomposition algorithm is proposed. The Tucker decomposition algorithm is used to reduce the dimension of the measured data, and the processed core tensor is used for the temperature ﬁeld reconstruction of sparse data. Theoretical analysis and simulations show that the proposed method is feasible; the overall optimization is realized by selecting the appropriate core tensor dimensions; and the reconstruction error is less than 3%. Results indicate that the proposed method can yield a reliable reconstruction solution and can be applied to real-time applications.


Introduction
Accurate reliable measurements of 3D flame temperature distribution are highly desirable to achieve an in-depth understanding of combustion and pollutant formation processes including flame flashback prediction, thermoacoustic instability control and combustion emission reduction. YU Guan et al. experimentally studied and verified that the periodic acoustic forcing effect [1] and strategic use of transient forcing, hysteresis and mode switching can reduce the thermoacoustic amplitude of the self-excited combustion system [2]. On this basis, the application scope of open-loop control was explored [3]. C. L. Myung et al. reviewed the particulate emission characteristics of internal combustion engine exhaust [4] and the spark ignition and diesel engine combustion process [5] and compared the effects of different vehicle certification modes on particulate emission characteristics of light vehicles through experiments [6]. C. Heeger et al. explored the mechanism of flame flashback in annular slots of burners and put forward a new hypothesis of tempering mechanism [7]. A. Nauert et al. used planar laser diagnostics at high repetition rates to quantitatively study the flame during flashback [8]. F. Seffrin et al. proposed a layered burner, verified the numerical simulation model through experiments and analyzed the flow field of the burner. The experimental apparatus and results are discussed [9]. Due to the thermal and dynamic nature of the flame, an effective means for the spatial and temporal measurement of the flame temperature remains a challenge for combustion and measurement researchers.
Significant efforts have been made in past decades to measure the flame temperature distribution. Various contact and non-contact methods have been proposed and developed. There are numerous methods for measuring temperature. The contact temperature measurement methods mainly refer to the thermocouples, widely used in the combustion process [10,11]. It is convenient and efficient, but it can only be used for single-point measurements. It is impossible to install enough equipment for detailed information. Non-contact measurement methods, such as acoustic measurement [12], laser diagnostic measurement [13] and radiation imaging methods [14], can be used to reconstruct the temperature field. These techniques have the advantages of a wide measurement range, fast dynamic response and little impact on the flame [15,16]. However, these methods can generally measure the parameters on the line or surface, and many problems remain that may limit their application; for example, noise, vibration, and ash in the industrial environment may have a great impact on their accuracy. Furthermore, sensors involved in these techniques are difficult to install, and the measurement systems are complicated and expensive. All of these factors increase the difficulty of applying these devices. Numerical methods, such as computational fluid dynamics (CFD), can also be used to obtain accurate and detailed temperature information. However, numerical computations always involve a large amount of data; they are exceptionally time-consuming and not suitable for real-time applications. Therefore, it is significant to propose and develop an effective method to use these data.
By extracting eigenvector by learning from experience dataset, principal component analysis (PCA)-also known as Karhunen-Loève (KL) expansion or proper orthogonal decomposition (POD)-can analyze and utilize the most essential characteristics of these data.
As the PCA is convenient and efficient, it is used in many different domains, including signal analysis [17], data compression [18], image processing [19], chemical reaction analysis [20] and control systems [21]. In recent years, it has also been well used in 2D temperature distribution and wind field analysis and reconstruction. Chen et al. [22] combines measurement and numerical algorithms based on PCA to obtain temperature distribution in real time by using feature vectors and measurement data. Sun et al. [23] used PCA for real-time wind field reconstruction for large areas and improving the efficiency of sensor applications.
In order to develop a 3D temperature field reconstruction method that can be used in real applications, inspired by the PCA, a three-dimensional temperature distribution reconstruction algorithm based on Tucker decomposition algorithm is proposed.
Tucker decomposition was first proposed by Tucker in 1963 [24]. Since then, Tucker decomposition has been applied in many different fields, including data compression [25], image processing [26], signal analysis [27] and so on. Zhao et al. [28] used it for infrared small-target detection, and Qin et al. [29] used it in wind velocity distribution reconstruction. The method has also been used in hyperspectral change detection [30].
In this study, a method based on Tucker decomposition algorithm is proposed for the rapid and reliable reconstruction of the 3D temperature field using sparse measurements, and it was verified by simulation data and experimental data, which provides a new solution for temperature distribution measurement.

Tucker Decomposition Algorithm Principle
By Tucker decomposition, the N-order high-dimensional tensor can be transformed into an n-order low-dimensional core tensor and n decomposition factor matrices (N > n). The Tucker decomposition of N-order tensors χ ∈ I×J×K can be written as [31] where U (n) , n = 1, . . . , N is the decomposition factor, and ε is N-order error tensor. Therefore, the factor matrix can be obtained by the following optimization problems: Considering that the Tucker decomposition of the N-order tensor is a higher-order singular value decomposition (HOSVD), In higher-order SVD, the factor matrix U (n) ∈ I n ×J n satisfies the semi orthogonal condition U (n)T U (n) = I J : Then, to solve equation Then, both sides multiply U (n)T left and U (n)T ⊗ right, respectively, and the following can be obtained from Equation (5): Therefore, the mathematical model of the n-order tensor χ ∈ I 1 ×I 2 ×I 3 ×...I N decomposition process is expressed as follows: where ζ-tensor χ Tucker decomposition core tensor, ζ ∈ J 1 ×J 2 ×J 3 ×...×J N , J N < I N ; × Ntensor χ N-mode product and U (N) -tensor χ decomposition factor matrix in the N-module direction, U (N) ∈ I n ×J n .

Reconstruction Algorithm Based on Tucker Decomposition
In the reconstruction of three-dimensional temperature distribution, the a priori dataset is composed of n 1 three-dimensional temperature field obtained by simulation calculation or measurement. In order to facilitate the subsequent calculation and analysis, the prior dataset is recorded as a fourth-order tensor TR = {T 1 , T 2 , T 3 , . . . , T i , . . . , T n1 }, T i ∈ n 2 ×n 3 ×n 4 , TR ∈ n 1 ×n 2 ×n 3 ×n 4 . The prior dataset is decomposed into the form of core tensor and decomposition factor mode product through Tucker decomposition, as follows: where C t -the core tensor of prior dataset tensor TR Tucker decomposition, C t ∈ m 1 ×m 2 ×m 3 ×m 4 and F t(i) -decomposition factor matrix in the i-module direction of tensor TR, F t(i) ∈ n i ×m i .
where P-the decomposition factor matrix in the 1-mode direction and Q-the decomposition factor matrix in the 2-mode direction.
According to Formula (10), Formula (9) is transformed into the following form: where f 1i -the ith factor vector in the decomposition factor matrix in the direction of 1 modulus and f 1i ∈ 1×m 1 . It can be seen from Formula (11): where A-reconstruction coefficient tensor, (4) and A ∈ m 1 ×n 2 ×n 3 ×n 4 .
In the process of solving, for the three-dimensional temperature distribution to be reconstructed under the same combustion environment, the Equation (13) should also be satisfied: where f 1x -the three-dimensional temperature distribution to be solved to reconstruct the 1-mode decomposition factor vector.
According to the definition of model product, the Formula (13) is transformed into the form of vector and matrix multiplication: where t x -first-order representation of three-dimensional temperature distribution; A 1unfoldthe expansion of the coefficient matrix A along the direction of 1 modulus and A 1unfold ∈ (n 2 ×n 3 ×n 4 )×m 1 .
In the temperature measurement process, when the k measuring point positions are determined, the corresponding measurement matrix can be constructed, and the temperature measurement process can be expressed as a mathematical model: where t M -a vector of small number of temperature measurements, t M ∈ R k×1 ; M-the measurement matrix, which characterizes the location of the measurement points and M ∈ R k×n . Combining Formula (14) and Formula (15), the relationship between the measured temperature and the coefficient vector to be solved can be established: In Formula (16), t M is obtained by measurement, M is a known condition after the position of the measuring point is fixed, and A 1unfold can be calculated by Tucker decomposition of the prior data set. Set the objective function of minimizing reconstruction error based on the measured data: where J f T 1x -the objective function.
Through Formula (17), the temperature distribution reconstruction is transformed into a problem of obtaining the optimal solution of the reconstruction coefficient vector: where (18), and then calculate the Formula (14) to obtain the threedimensional reconstructed temperature distribution. The specific process of temperature field reconstruction by Tucker decomposition is shown in Figure 1 below.

Numerical Simulation and Analysis
On the basis of the previous research, in order to reduce the computational cost, numerical simulation is used to quantitatively demonstrate the feasibility of the proposed algorithm based on the Tucker decomposition and determines the effectiveness of the method.
A space of 4 m × 4 m × 4 m is constructed as a 3D temperature field reconstruction area. Set three different 3D temperature distribution numerical models as the research objects of the temperature reconstruction in this section, which are used to calculate the prior dataset and the test dataset. Apply the proposed 3D temperature distribution reconstruction algorithm to reconstruct the temperature profile under different test conditions. The effectiveness and feasibility of the reconstruction algorithm are analyzed. Restricted by the existing calculation conditions, the three-dimensional temperature distribution data is characterized by selecting 2D temperature field sections at different heights (0.5 m, 1.5 m, 2.5 m and 3.5 m).

Single-Peak 3D Temperature Distribution Reconstruction
Set the single-peak 3D temperature distribution model as follows: where T 1 (x, y)-temperature at the points (x, y, z) of the single-peak 3D temperature distribution (K); x-X coordinate(m); y-Y coordinate(m); z-Z coordinate(m) and a-initial boundary condition.
In the reconstruction area, use Equation (19) to calculate the temperature value. In order to reduce the amount of calculation, the X-axis and Y-axis directions take 0.1 m as the step length, and the Z-axis direction takes 1 m as the step length to start sampling from the origin. To form 3D temperature distribution data, 4 × 40 × 40 temperature values are collected. Calculate the temperature field under different boundary conditions, construct a priori dataset with dimensions of 30 × 4 × 40 × 40 and test datasets with dimensions of 5 × 4 × 40 × 40. Set the prior dataset and the test dataset of the single-peak 3D temperature distribution model and the boundary conditions as shown in Appendix A  Tables A1 and A2. According to the 3D temperature distribution reconstruction calculation process, the prior dataset is decomposed by the Tucker decomposition algorithm to obtain the reconstruction coefficient matrix. Select 16 coordinate points as temperature measuring points to reconstruct the 3D temperature profile.
The temperature distribution of test condition 5 calculated is shown in Figure 2. The temperature field result of the test condition 5 obtained by the Tucker algorithm is as follows.
The boundary condition of test condition 5 is outside the constraint conditions of the priori dataset. As shown in Figure 3, the proposed algorithm still realizes the effective reconstruction. To analyze the adaptability of the 3D temperature distribution reconstruction algorithm under different working conditions, the temperature profile data of the five test working conditions are reconstructed, and the reconstruction error is shown in Figure 4 below. As shown in Figure 4, the temperature field reconstruction results under five test conditions are analyzed. It can be seen that in the construction of single-peak 3D temperature distribution, the three-dimensional temperature field reconstruction algorithm based on the Tucker decomposition can effectively reconstruct the temperature profile data under different conditions.

Two-Peak 3D Temperature Distribution Reconstruction
Set two-peak 3D temperature distribution as follows: where T 2 (x, y)-points (x, y, z) of two-peak 3D temperature distribution model (K); and a, b-initial boundary conditions. Set the prior dataset and test dataset of the two-peak three-dimensional temperature distribution model and the boundary conditions as shown in Appendix A Tables A3 and A4.
The data acquisition and reconstruction process of the two-peak temperature distribution model is consistent with the calculation procedure of the single-peak mentioned above. The calculated temperature distribution of test condition 5 and the reconstruction result of the temperature field of test condition 5 calculated by the algorithm are shown in Figure 5. The 3D temperature field reconstruction algorithm is used to reconstruct the temperature distribution data of the five test conditions. The error is shown in Figure 6. Analysis of Figures 5 and 6 shows that in the five test conditions of the double-peak 3D temperature distribution model, the proposed algorithm can effectively realize the temperature distribution reconstruction.

Three-Peak 3D Temperature Distribution Reconstruction
Set the three-peak 3D temperature distribution as follows: where T 3 (x, y)-points (x, y, z) of three-peak 3D temperature distribution model (K). Set the model prior dataset and test dataset boundary conditions as shown in Appendix A Tables A5 and A6.
The calculated temperature distribution of test condition 5 and the reconstruction result of the temperature field of test condition 5 calculated by the algorithm are shown in Figure 7 below. The Tucker decomposition algorithm is used to reconstruct the temperature distribution data of the five test conditions. The error is shown in Figure 8. An analysis of Figures 7 and 8 shows that in the five test conditions of the three-peak 3D temperature distribution model, the 3D temperature distribution reconstruction algorithm proposed in this paper can effectively realize the temperature distribution reconstruction.
Based on the above analysis, this section uses the calculation data of single-peak, two-peak and three-peak three-dimensional temperature distribution models to verify the effectiveness and practicability of the proposed algorithm. The reconstruction results illustrate the feasibility of applying Tucker decomposition algorithm to the 3D temperature distribution reconstruction.

Simulation Experiment
Based on the preliminary demonstration of the effectiveness of the algorithm, in order to analyze the application effect of the reconstruction algorithm in the combustion, this paper takes the Flame D in Sandia laboratory as the research object and builds a numerical simulation model of methane combustion. The practicability of the reconstruction algorithm is analyzed by the simulation calculation data, and the validity of the temperature distribution reconstruction algorithm proposed in this paper is further demonstrated.

Methane Burning Duty Flame
The duty flame is a small flame set around the high-speed fuel outlet of the burner to achieve ignition and stable combustion. At Sandia Laboratories, Flame D is a typical jet flame obtained by burning at initial conditions with a methane inlet velocity of 49.6 m/s and a hot air velocity of 11.4 m/s formed by the duty flame [32,33].
Sandia Labs has measured the flame burning area of Flame D through laser-induced fluorescence, and the measurement data has been made public. The burner outlet structure used in Flame D is shown in Figure 9. The geometric dimensions of the burner outlet are as follows: the diameter of the fuel outlet is d = 7.2 mm; the diameter of the wake heat flow outlet is 18.2 mm, which is an annular structure; the outside of the heat flow outlet is air flow. In order to facilitate the establishment of the model in the subsequent process, the diameter d of the fuel outlet is used as the unit length to characterize the size and coordinate position [33].
The fuel consists of methane and air mixed in a volume ratio of 1:3. The accompanying heat flow is formed by the combustion of mixed gases such as acetylene and hydrogen. The specific parameter data are shown in Table 1.

Numerical Simulation Model Settings
Flame D is a typical turbulent jet flame. The Reynolds number of the flame is 22,400, and the fuel inlet velocity is relatively fast, which belongs to the fully developed fast turbulent combustion process. The Realizable model is used in the simulation calculation. According to the actual parameters, the following settings are made: a cylindrical area with a diameter of 40 d and a height of 80 d is designated as the calculation area, and it is divided into 1,823,760 grids. The geometry of the Flame D combustion simulation model is shown in Figure 10.  In order to ensure that the numerical simulation model can accurately simulate the Flame D combustion process, the parameters of the model are set based on the actual operating conditions, and appropriate optimization is carried out on the basis of taking into account the calculation efficiency and accuracy.

Numerical Simulation Calculation Results Analysis
Sandia Labs has released a measurement data set of Flame D, taking the cross-sections formed by +z and +x in the model as the research object. Compare the simulation data with the measured data to analyze the accuracy of the simulation.
The comparison between the simulation calculation data of the z-axis observation point and the temperature measurement data of the corresponding position is shown in the following figure.
It can be seen from Figure 11 that there is a certain deviation between the simulation results and the measured data. Analyzing the influencing factors, the error between the calculation result and the measured data mainly comes from the solution process and the model building process. Since the grid independence analysis has been carried out, the influence of grid division on the calculation results can be ignored, so the error in the solution process mainly comes from the discrete process of the governing equations. In the process of model establishment, the selection of turbulence model, chemical reaction model and parameter setting, as well as the simplification and assumption in the process of model building, are the main factors that cause the deviation between the calculated results and the measured data. Nevertheless, as shown in Figure 11, the temperature calculated values at different positions are basically consistent with the general trend of the actual temperature distribution, that is, the simulation calculation results are close to the actual measured values.
In addition, the z-axis coordinates are taken as 30 d, 45 d, 60 d and 75 d for the radial observation point temperature data comparison as follows.
Similar to the conclusion in Figure 11, the calculated temperature values at different radial positions shown in Figure 12 are basically consistent with the overall trend of the actual temperature distribution. Therefore, the simulation model can better simulate the temperature distribution of Flame D. In summary, although the numerical calculation method has approximations and assumptions that cannot be avoided in the process of solving and model building, it can still accurately simulate the basic properties of the simulation object and meet the actual needs of engineering applications. The methane combustion simulation model established in this paper can simulate the real temperature distribution of Flame D accurately.

Prior Data and Test Data Acquisition
As mentioned in the previous section, the established numerical simulation model can simulate the methane combustion process, and the calculated temperature distribution data is close to the actual value and can reflect the distribution of the actual combustion temperature. Therefore, this model can be applied to build prior datasets as well as the test dataset for the temperature distribution reconstruction. Considering the turbulence intensity and the range of the prior dataset, the boundary conditions for methane combustion simulation calculation are set as shown in Table A7.
In order to construct the temperature distribution reconstruction calculation samples, five test cases were set up. The boundary conditions of the test condition are different from the prior condition, and the boundary condition of the test condition 5 is set outside the range of the boundary condition of the prior condition to test the robustness of the temperature distribution reconstruction algorithm. The boundary condition settings of the test conditions are shown in Table 2.

Simulation Data Reconstruction Analysis
The feasibility of the reconstruction algorithm is verified by using the temperature distribution data calculated from the numerical simulation model.
The temperature of 1600 locations on the three planes of z = 30 d, 45 d and 60 d of the methane combustion simulation model were selected as the research objects of the three-dimensional temperature distribution reconstruction. The selected positions include 36 measuring points on the axial 30 d, 45 d, and 60 d published by Sandia Labs so as to facilitate the experimental data reconstruction and analysis. Based on the simulation calculation results of the prior and test conditions, a three-dimensional temperature distribution reconstruction prior dataset with dimensions of 24 × 3 × 40 × 40 and a test dataset with dimensions of 5 × 3 × 40 × 40 were constructed.
The 3-D temperature distribution reconstruction algorithm is used to reconstruct the test conditions' temperature distribution. The measurement data is taken from the simulation calculation data. Based on the measurement matrix, the temperature of the corresponding position in the test dataset is extracted as the temperature measurement data, which is used for the temperature distribution reconstruction.
Taking the test case 5 as an example, the simulation results of the temperature distribution in the three target sections are shown in Figure 13.  Figure 13 shows the simulation results of the three-dimensional temperature distribution calculated by using the simulation numerical software according to the geometric parameters of the established model, combustion parameters and boundary conditions of test condition 5. The white points in the Figure 13 represent the arrangement of temperature measuring points.
The temperatures of the 12 measuring points shown in Figure 13 are used as temperature measurement data, and the reconstruction calculation is performed using the core tensor with dimensions of 1 × 2 × 2 × 3.
After reconstruction calculation, it can be seen that the reconstruction error of test case 5 is 4.13%, and the temperature deviation of the reconstruction results is less than 90 K, which realizes the effective original temperature distribution reconstruction and verifies the 3-D temperature distribution reconstruction algorithm.
In order to better show the reconstruction results, the 3D temperature distribution reconstruction calculation and analysis process presents the reconstruction effect as a contour map of the absolute value of the deviation between the reconstructed temperature value and the calculated temperature value. The absolute value of the temperature deviation between the reconstruction calculation results and the simulation calculation results is shown in Figure 14.
In order to analyze the adaptability of the reconstruction algorithm under different working conditions, the three-dimensional temperature distribution reconstruction algorithm was used to reconstruct the five test conditions' data, and the reconstruction errors are shown in Figure 15.
As shown in Figures 14 and 15, the three-dimensional temperature distribution reconstruction algorithm can effectively reconstruct the five different test conditions' data, which further verifies the practicability and effectiveness.

Effect of Core Tensor Dimension on Reconstruction Accuracy
On the basis of the validity of the reconstruction algorithm, the influence of the core tensor dimension on the reconstruction accuracy during the temperature distribution reconstruction is analyzed to further improve the algorithm accuracy. The temperatures of 1600 positions were selected on the three planes of z = 30 d, 45 d and 60 d in the methane combustion simulation model as the research object of 3D temperature distribution reconstruction.
The three-dimensional temperature distribution reconstruction algorithm based on Tucker decomposition uses the key information of the prior dataset carried by the core tensor to realize the reconstruction. The core tensor of the prior dataset has four dimensions, so it is necessary to analyze the values of the four dimensions one by one.
The prior dataset of 3D temperature distribution reconstruction is a tensor of 24 × 3 × 40 × 40. The dimension in the 2-modular direction in the prior data is small, so the calculation amount of optimizing the dimension value in the 2-modular direction is also small. Therefore, the effect of the core tensor dimension on the reconstruction accuracy is explored from the 2-modular direction.
Set up the 2-module direction dimension influence analysis experiment, and the experimental steps are as follows: (1) In order to avoid too few measurement data points affecting the analysis results, 50 measurement data points were used in the reconstruction. In order to avoid the influence of the dimension of other modes, the dimension of other modes is fixed to 10; (2) Perform Tucker decomposition on the prior dataset to obtain a core tensor with a set dimension of 10 × n × 10 × 10, n = 1,2,3; (3) Use the core tensor, decomposition factor matrix and measured temperature to reconstruct the 3D temperature field, and calculate the reconstruction error. Add noise data with SNR = 40 dB to the measured temperature and compare the reconstruction results in the presence of measurement error; (4) Repeat the above steps 1000 times to obtain the average reconstruction error in order to eliminate the influence of randomness in the temperature measurement so a universal conclusion can be obtained; (5) Repeat the above steps until the dimension n has traversed all the values.
Taking the test case 5 as an example, the different 2-mode dimensions reconstruction results are shown in Table 3. As can be seen from Table 3, as the dimensionality in the 2-mode direction increases, the reconstruction error decreases whether using the original temperature measurement data or using the noise-added measurement data. Therefore, the 2-modular dimension of the core tensor is set to be three in the follow-up research work.
When the 2-module dimension is three, the effect of the 1-module dimension on the reconstruction accuracy is analyzed. The 3-module and 4-module dimension are set to 10, and the value range of the 1-module dimension n is 1 to 10. Other parameters and calculation steps are consistent with the above analysis and will not be repeated here.
Taking the test case 5 as an example, the reconstruction results of different 1-module dimensions are shown in the table below: In Table 4, it can be seen that the reconstruction error using the original temperature measurement data will gradually decrease with the dimension increase; the reconstruction error using noise-added measurement data increases with the increase of dimension. From the perspective of balancing the reconstruction accuracy and noise sensitivity, the 1-module dimension of the core tensor is set to two in the subsequent reconstruction work. Similarly, it is determined that the 3-module dimension is five, and the 4-module dimension is two. Finally, the optimized core tensor dimension is determined to be 2 × 3 × 5 × 2. The optimized core tensor dimension is used to add noise with different signal-tonoise ratios to Test Condition 5, and the comprehensive performance of the optimization algorithm under different noise conditions is investigated, as shown in Figure 16. The data of five test cases are reconstructed, and the reconstruction effect with SNR = 40 dB noise-added is analyzed. The results are shown in Figure 17.
It can be seen from Figures 16 and 17 that the optimization effect of the core tensor dimension is stable in the five test conditions. Compared with the reconstruction error, the optimization algorithm can reduce the reconstruction error under different test conditions, and the optimization effect is more obvious with the increase of noise intensity. The optimized reconstruction algorithm has a high reconstruction accuracy and good anti-noise performance.

Conclusions
In this paper, a fast and accurate optimization algorithm based on Tucker decomposition has been proposed to reconstruct 3D temperature fields with sparse data. The relevant mathematical derivations from Tucker decomposition have been carried out to determine a 3D temperature field reconstruction algorithm. The dataset used to generate the core tensor plays a very important role in this method. When using CFD results to establish datasets, it is essential to sample large-scale boundary conditions that can cover most of the working conditions in the flow field. Simulations were performed, and numerical results demonstrated that the Tucker decomposition algorithm combined with sparse measurement is capable of reconstructing the 3D temperature field quickly and accurately. To evaluate the performance of this algorithm, a numerical simulation model of methane combustion has been established based on the Flame D of Sandia Laboratory. Simulation experimental results showed that this method can reconstruct temperature distribution rapidly using limited measurement data with the relative of less than 5.08%, which verifies the feasibility of this algorithm. It can greatly reduce the computational complexity and measurement cost. Finally, this paper analyzes the influence of core tensor dimension on reconstruction accuracy in the process of a temperature distribution reconstruction calculation. Core tensor can be extracted from a high-dimensional experience dataset. The first few dimensions carry most of the useful information. With the use of a higher dimension core vector, too much useless information is introduced to the reconstruction process, and the reconstruction results are more sensitive to the noise. Thus, it is necessary to obtain a balance between reconstruction accuracy and noise sensitivity. The results show that the reconstruction error is less than 2% by optimizing the core tensor dimension.
Furthermore, the present study lays the groundwork for future research into some interesting fields with several objectives, including the following: 1.
In order to show more object field characteristics, this study uses a wide range of boundary conditions to establish the dataset. This helps to generate the core tensor of Tucker decomposition, but how to construct dataset more effectively remains to be solved. Further research can focus on the source of empirical data, effective dataset sampling and the impact of dataset noise.

2.
The proposed method needs little measurement data to yield an acceptable reconstruction result. The sensor number has some effect on the reconstruction result: the greater the number of sensors used, the better result. Therefore, it is necessary to explore the number of appropriate sensors needed to balance the measurement cost and reconstruction accuracy.

3.
When only a few sensors are used for reconstruction, sensor placement has a great influence on reconstruction accuracy. Testing the feasibility of finding a sensor placement optimization algorithm that can yield a much better reconstruction result.
Appendix A  Table A2. Boundary conditions of test dataset (single-peak model in a cubic area).

Working Condition Boundary Condition (a)
Test Condition 1  723  Test Condition 2  1107  Test Condition 3  1426  Test Condition 4 1846 Test Condition 5 2200 Table A3. Boundary conditions of prior dataset (Two-peak model in a cubic area).   Table A6. Boundary conditions of test dataset (Three-peak model in a cubic area).