Issue 
Sci. Tech. Energ. Transition
Volume 77, 2022
Selected Papers from 7th International Symposium on Hydrogen Energy, Renewable Energy and Materials (HEREM), 2021



Article Number  16  
Number of page(s)  16  
DOI  https://doi.org/10.2516/stet/2022013  
Published online  13 July 2022 
Regular Article
Inverse estimation of wall temperature measurement under dispersion medium shielding based on CMAES algorithm
^{1}
China Special Equipment Inspection and Research Institute, Beijing 100013, PR China
^{2}
School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, PR China
^{3}
School of Aeronautical Science and Engineering, Beihang University, Beijing 100191, PR China
^{*} Corresponding authors: qianlinfeng@csei.org.cn; wudi@csei.org.cn
Received:
5
November
2021
Accepted:
21
April
2022
Special equipment is often accompanied by the phenomenon of high temperature and high pressure in operation, and it needs to be measured. With the deepening of research engineering application, some new problems appear, such as how to accurately measure and control the combustion temperature field in the hydrocarbon energy system, carry out combustion diagnosis, control the combustion temperature, minimize the generation of pollutants, put forward higher requirements for measurement technology. High temperature has a core impact on equipment safety, which is a key index that must be considered in design and material selection. Near the surface of these hightemperature objects, there are usually hightemperature translucent media and particles produced by combustion, which have a great impact on radiation temperature measurement. At the same time, with the change of temperature and concentration, these translucent media often have gradient refractive index change, which has a certain impact on radiation temperature measurement. These shielding objects interfere with radiation and affect the accuracy of measurement. This topic mainly focuses on the theoretical and experimental research of this kind of problem. The inversion algorithm of wall temperature reconstruction under the shielding condition of dispersed medium is proposed. Under the condition of introducing different measurement errors, the algorithm is verified, and the corresponding numerical simulation is carried out to verify the feasibility of the algorithm model.
Key words: Inverse radiation problem / Gradient refractive index / Hightemperature measurement / Special equipment
© The Author(s), published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Nomenclature
χ: The derivative of the refractive index
μ: Number of individuals selected
c_{1}: Rankone update coefficient
y: The covariance matrix update
μ_{eff}: The variance effective selection mass for the mean
F_{fitness}: Objective function
I_{exp}: The directional radiation intensity calculated by simulation experiments
F_{obj}: The objective function
a_{1,} b_{1,} a_{2,} b_{2}: The coefficient of linear distribution
φ: The ratio of wall temperature to medium temperature
w: Positive weight coefficients
c_{μ}: The rankμ update coefficient
N: A multivariate normal distribution
I_{est}: The directional radiation intensity calculated according to the solution vector
I_{*}: The outgoing radiation or radiation response signal
ς: The standard normally distributed random number
X: The value of the parameters to be inverted
x: The normalized spatial position
κ_{a}: Absorption coefficient, m^{−1}
σ: Stefan–Boltzmann constant 5.67 × 10^{−8} W/(m^{2} · K^{4})
1 Introduction
A large amount of heat generated in the combustion process of power plant boilers is mainly transmitted to the medium in the form of radiation, convection, and heat conduction [1]. Combustion is a violent chemical reaction process and the primary process of energy conversion. For furnace flame, hightemperature flue gas and metal temperature of boiler or heat exchanger are some of the important factors to measure the stability of combustion and operation condition. On fire or hightemperature objects in radiation temperature measurement technology involves internal radiation transmission and many problems [2]. The calculation of radiative transfer in participatory media usually assumes that the refractive index of a flame uniform field, namely the spread of the light always along a straightline trajectory through the research object. So although to some extent simplifies the calculation process and improves the computational efficiency. However, the refractive index tends to be unevenly distributed due to the uneven spatial distribution of the temperature, pressure, or composition of the flame [3]. If the refraction index distribution of semitransparent hightemperature mediums such as flame and hightemperature flue gas is regarded as uniform for temperature field reconstruction, the measurement accuracy and measurement range will inevitably be reduced. For metal temperature receiving radiant heat or being scoured by hightemperature flue gas, high precision temperature measurement is also required in operation. Due to the presence of flame or hightemperature smoke and other participating media on its surface, the influence of gradient refractive index formed by the change of temperature and concentration on the temperature measurement results cannot be ignored [4–6]. Therefore, it is of great scientific and practical value to study the refractive index gradient distribution and its influence on radiation transfer in flame or hightemperature semitransparent medium in actual production.
Among many temperature measurement technologies, based on features of flame radiation on temperature field inversion temperature measurement technology is the new development direction with a relatively simple device, the measurement is relatively direct, the advantage of high reliability. Such as lightfield imaging technology, its core is still the image method used to obtain the boundary of the flame radiation intensity distribution, the temperature, and optical parameters inversion. At present, light field imaging technology has great development space in the field of combustion diagnosis [7]. Essentially, it can be concluded as the reconstruction of multiple physical quantity fields in a hightemperature medium based on space radiation light field information. This is a typical radiation inverse problem with obvious discomfort [8], so it is quite difficult to solve it.
There is usually a translucent medium between the measured object and the probe in the research of furnace flame or hightemperature object temperature measurement technology. Precise gradient refractive index medium of these areas is involved in the information transmission process of field which can give people more indepth understanding of the radiation energy transfer mechanism and the radiation law of information transmission and master the gradient refractive index medium radiation transmission characteristics of field related optimization design method, manufacture technology and measurements [8–12]. Therefore, there is important theoretical significance engineering application value in studying the radiation characteristics and temperature measurement technology under the condition of gradient refractive index medium shielding.
Radiative transfer in graded refractive index media is much more complex than that in uniform refractive index media. Light in a uniform refractive index medium travels along a straight line. But in a gradient refractive index medium, the light will travel along a curved path, and sometimes even total reflection will occur inside the medium [13]. Therefore, the traditional radiative transfer method for uniform refractive index media cannot be directly applied to solve the radiative transfer in graded refractive index media. In recent years, the radiative transfer in graded refractive index media has attracted more and more attention from scholars all over the world.
Methods for calculating radiation transfer in graded refractive index media include Discrete Ordinates Method (DOM) [14–16], Finite Volume Method (FVM) [3, 17, 18], Finite Element Method (FEM) [19, 20], Spectral Element Method (SEM) [21, 22], etc. The global discretization method based on differential form radiative transfer equation can be used to deal with multidimensional complex geometric area problems conveniently and economically and complex medium and boundary characteristics problems flexibly. It can also be easily used to solve the coupled heat transfer problems of radiation, heat conduction, and convection. Especially for the finite volume method and discrete coordinate method because of their strong adaptability.
Wei et al. studied the inverse estimation of refractive index related to temperature/phase transition in translucent media. The refractive index is inverted from the radiative response of the boundary by a Stochastic Particle Swarm Optimization (SPSO) algorithm. Through the sensitivity analysis, it is proved that this method can accurately deduce the refractive index from the radiation response on the two boundaries [23]. Zhao et al. studied linear and nonlinear gradient refractive index distributions based on discrete coordinate equations. Temperature and radiative heat flow distribution in isotropic and heterotropic media [24], transient radiative transfer equation [25, 26], and coupled radiative heat transfer in multidimensional homogeneous media with arbitrarily complex shapes were solved by spectral element method. Wei et al. found that the refractive index and absorption coefficient differ in sensitivity to boundary conditions, making it difficult to estimate spatially correlated refractive index and absorption coefficient simultaneously. A twostep inversion model with different laser irradiations is proposed for this problem [27].
Although the combustion heat transfer problem has been solved by combining theory with experience. It has been widely used in large power station boilers and industrial combustion devices. With the deepening of research engineering application, some new problems appear, such as how to accurately measure and control the combustion temperature field in the hydrocarbon energy system, carry out combustion diagnosis, control the combustion temperature, minimize the generation of pollutants, put forward higher requirements for measurement technology.
This paper mainly focuses on the influence of passive radiation on the temperature measurement of flame and boiler components under the condition of translucent medium with a gradient refractive index. The reconstruction program of gradient absorption coefficient and refractive index of functional distribution is established based on Evolution Strategy based on Covariance Matrix Adaption (CMAES) algorithm. The temperature ratio was defined and the sensitivity of distribution function coefficients was analyzed under different temperature ratios, outgoing radiation, and different surface normal angles to obtain a better combination of measurement information. The temperature ratio and the included angle are selected according to the graph line of sensitivity coefficient. The form of the distribution function of the gradient absorption coefficient and refractive index of the discrete distribution in the process of reconstruction is unknown to consider the more general distribution of photothermal parameters. The sensitivity of absorption coefficient and refractive index is analyzed for the twolayer medium model. The sensitivity graph of absorption coefficient and refractive index of the threelayer medium model is analyzed and explained the rationality of the image drawn. An inversion algorithm for wall temperature reconstruction with variable refractive index medium is proposed.
2 Direct problem
Consider an absorbing, emitting, and scattering 1D parallelplane medium with variable refractive index and azimuthally symmetric radiation bounded by two diffusely reflecting plates. As shown in Figure 1, both sides of the boundary are assumed to be the diffuse gray opaque wall with the same emissivity. The left and the right boundary walls are subjected to the first kind of boundary condition with constant temperature T_{0} and T_{L} respectively. The absorption coefficient and scattering coefficient are supposed to be constant. The refractive index n(x) varies with axis x.
Fig. 1 The physical model of radiative heat transfer in 1D participating media with GRI. 
The radiative transfer equation in a 1D participating media with gradedindex in Cartesian coordinate can be expressed as,(1)
Boundary conditions for greydiffuse walls can be expressed as,(2) (3)where I(x, μ) is radiation intensity at position x and direction cosine μ, I_{b} is the blackbody radiation intensity at the temperature of the media T(x), Φ(μ′, μ) is the scattering phase function of energy transfer from an incoming direction cosine μ′ to an outgoing direction cosine μ, σ is Stefan–Boltzmann constant, κ_{a} is the absorption coefficient, κ_{s} is the scattering coefficient and χ means the derivative of the refractive index after natural logarithm to the coordinate x,(4)
3 Inverse problem
3.1 CMAES algorithm
The CMAES algorithm has been developed to be a wellestablished algorithm system that contains many variants for different types of problems with the efforts of Hansen, Ostermeier, and their colleagues [28–30]. In this study, the integral equations of light scattering signals need to be deconvoluted to reconstruct the size distribution and morphological parameters. Mathematically, it is a typical illconditioned multiple dimensional nonseparable inverse problem. Given the good performance of CMAES in solving the illconditional and nonseparable problems [29–31], this promising algorithm is suitable to solve the inverse problem of light scattering. What follows is a brief introduction of the theoretical fundamentals of the CMAES algorithm, heavily based on a tutorial written by Hansen [29] which the reader can refer to for a full and detailed derivation.
In the CMAES, if the number of objective parameters to be inverted is n, then the search point is an ndimensional vector, which is also called the objective vector. The basic equation of CMAES is to generate the new search points by sampling a multivariate distribution [29]:(5)where is the kth offspring individual in the population of λ new search points at generation g + 1, ~ denotes the same distribution on the left and right side, m^{(g)} ∈ R^{n} is the mean value of the search distribution at generation g, σ^{(g)} ∈ R^{n} is the stepsize at generation g, denotes a multivariate normal distribution with zero mean and covariance matrix C^{(g)}, and m^{(g)} + σ^{(g)}N(0, C^{(g)}), C^{(g)} ∈ R^{n × n} is the covariance matrix at generation g, λ ≥ 2 is the population size or the number of the offspring individuals. The kernel of CMAES is to update distribution parameters m^{(g)}, σ^{(g)} and C^{(g)}. Update of the mean value m^{(g)} is achieved by selection and recombination, which means the new mean m^{(g + 1)} of the search distribution is a weighted average of μ selected points from samples , ,…, [29]:(6)where μ ≤ λ is the parent population size, i.e. the number of selected points, are positive weight coefficients for recombination such that w_{1} ≥ w_{2} ≥ … ≥ w_{μ} > 0 and is the ith best individual out of , ,…, and the index i. λ denotes the index of the ith ranked individual, the “best” and “ranked” individuals are both defined by their objective function value. The smaller the individual’s objective function value is, the better the individual is, and the higher its ranking is where f is the objective function to be minimized.
The update of the covariance matrix is realized by the covariance matrix adaption based on rankμ update and rankone update. The rankμ update aims to give the estimation of covariance matrix reliable for the small population. The rankone update means that there is only one search point per generation for updating the covariance matrix in limited cases. Finally, rankμ update and rankone update are combined to generate a new covariance matrix [29]:(7)where c_{μ} ≤ 1 − c_{1} is the learning rate for the rankμ update of the covariance matrix update, c_{1} ≤ 1 − c_{μ} is the learning rate for the rankone update of the covariance matrix update, , h_{σ} is the Heaviside function, δ(h_{σ}) = (1 − h_{σ})c_{c} (2 − c_{c}) ≤ 1, is the evolution path at generation g + 1. The evolution path can be expressed by a sum of consecutive steps which is referred to as cumulation [29]:(8)where c_{c} ≤ 1 is the learning rate for cumulation for the rankone update of the covariance matrix, μ_{eff} denotes the variance effective selection mass for the mean. The factor is a normalization constant for p_{c}.
The evolution path is also utilized to control the step size. Unlike equation (8), applying the same technique in equation (9) is to build a conjugate evolution path as [29]:(9)where is the conjugate evolution path at generation g, c_{σ} < 1 is the learning rate for the cumulation for the stepsize control, is a normalization constant, , where is an eigen decomposition of C^{(g)}.
To update the stepsize σ^{(g)}, the matrix norm of conjugate evolution path is compared with its expected length E N(0, I), that is [29]:(10)where d_{σ} ≈ 1 is the damping parameter for the stepsize update, I is the unity matrix. In the present study, the strategy parameters settings of CMAES are available in [29].
To verify the performance of the CMAES algorithm in the inversion of medium photothermal parameters based on the inverse radiation problem, the forward radiation problem model of a onedimensional translucent plate considering selfradiation is selected to construct the radiation inverse problem. The CMAES algorithm is tested from three aspects: the calculation time, the calculation results, and the value of the bodyfitting function. The direct problem model and the inversion experiment flow are introduced below. The schematic diagram of the model is shown in Figure 2.
Fig. 2 Onedimensional translucent plate model. 
Only the absorption and emission characteristics of the medium are considered, L is the thickness of the medium, which is set as 1 m. κ is the absorption coefficient, which is set as 1.5 m^{−1}. Both sides of the medium are free boundaries. The radiative transfer equation is shown in equation (11):(11)where ω_{λ} is the extinction coefficient. The second and third terms on the right side of the equation contain the scattering caused by emission and incident in all directions of space, which is collectively known as the source function. The source function is,(12)
The radiation transfer equation in source function form can be written as,(13)
The flow chart of the inversion experiment is shown in Figure 3. The test was carried out at the same wavelength. It was assumed that the radiation source function was independent of the extinction coefficient and was only a quadratic polynomial of the optical thickness. The testing process mainly includes a direct problem model and optimization algorithm. During the test, the input parameter values were assumed first to be used to substitute the direct problem model to calculate the simulated measurement results. A total of 10 measurement signals were taken as the measurement results, including 5 backward radiation signals on the left boundary and 5 forward radiation signals on the right boundary. Then, the simulation measurement results and the direct problem model are brought into the optimization algorithm, and the optimization result is obtained through iterative calculation. The objective function of the optimization algorithm is,(14)where the subscript “est” and “exp” represent the directional radiation intensity calculated according to the solution vector and the directional radiation intensity calculated by simulation experiments, respectively. Finally, the error of the optimization result is calculated according to the assumed values of the input parameters.
Fig. 3 Flow chart of inversion experiment. 
This paper uses assumed parameters to test the CMAES algorithm, including no measurement error, adding 1%, 5%, 10%, 7% measurement error, and conditional 3% measurement error. During the process, the optimization algorithm was independently run 10 times, the average calculation time, success rate, average calculation result, and average calculation error were statistically analyzed. The search interval of the algorithm is [−20, 20].
Table 1 shows the average calculation results and errors of the source function coefficient inversion program of the radiation transfer equation based on the CMAES algorithm running independently 10 times. From the standard deviation of the calculated results, the results of the 10 independent measurements are relatively close which indicates that the algorithm has converged at the end of the operation. When no measurement error is added, the program can calculate more accurate measurement results. With the increase of addition error, the inversion error increases gradually and the increased amplitude is relatively slow. When 7% measurement error is added, that is the maximum measurement error adopted in this paper. The maximum error of inversion results is only 0.82%. It can be seen that the CMAES algorithm has high robustness in radiation inverse problems and the inversion results are less affected by measurement error.
Average calculation results and errors of running tests.
The average fitter function value statistics and average computing time statistics of 10 independent runs of the inversion program are analyzed in Figure 4. As can be seen from the left figure, although the fitness function value obtained by the optimization algorithm increases exponentially with the increase of measurement error, the fitness function value corresponding to adding 9% measurement error is only 0.112. It can be seen from the picture on the right, after adding the error of measurement, calculation time is around 0.5 s, namely under different measurement errors, the number of iterations of the same algorithm running. This algorithm can be attributed to the unique parameter selection mechanism, namely the CMAES algorithm by the number of iterations and by the inversion parameters calculated, but this does not affect the convergence of the algorithm itself.
Fig. 4 Fitness function values and calculation time statistics of running tests. (a) Fitter function value, (b) computation time. 
In summary, the CMAES algorithm has high robustness, strong optimization ability, and high computational efficiency when applied to the radiation inverse problem field. It is verified that the algorithm used in this paper is correct to invert the photothermal characteristic parameters of the inhomogeneous dispersion medium.
3.2 Inverse process
In general, the photothermal characteristic parameters of the dispersion medium are inhomogeneous. This section studies the inversion of the radiation characteristic parameters such as the absorption coefficient and refractive index of the dispersion medium. The onedimensional radiative heat transfer model of inhomogeneous dispersion medium considering unilateral wall radiation is shown in Figure 5. The left side is the opaque boundary. In this study, it is assumed that the wall emissivity is taken, and the right side is the free boundary. The absorption coefficient and refractive index of the dispersion medium vary with the spatial position.
Fig. 5 A physical model of radiation heat transfer in onedimensional inhomogeneous disseminated media considering radiation from one sidewall. 
Based on the evolutionary strategy algorithm of the adaptive covariance matrix tested above, this paper constructed a reconstruction algorithm for the photothermal characteristic parameters (including absorption coefficient and refractive index) of the heterogeneous dispersion medium. The inversion model of the photothermal parameters of the dispersion medium is shown in Figure 6. In the figure, the three colors respectively represent three important components of the nonmodel, namely the direct problem model, CMAES algorithm, and parameter inversion experiment. The direct problem model part refers to the previous content. In the algorithm part, as the input of the algorithm, the radiation response signals in different directions on the right boundary at different wall temperatures are selected as the measurement information. Based on the CMAES algorithm, the gradient refractive index and absorption coefficient distribution of semitransparent media were reconstructed. The selected objective function is as follows:(15)where F_{obj} is the objective function. (T_{i}, θ_{j}) is the outgoing radiation or radiation response signal on the right side in the θ_{i} direction when the left wall temperature is T_{i} in different cases. The subscript “est” and “exp” represent the calculated radiation response signal and the measured radiation response signal according to the current solution position in the search space, respectively. T_{i} is the sampling temperature, T_{i} ∈ [300, 1000]. θ_{i} is the Sampling angle, θ_{i} ∈ [0°, 90°).
Fig. 6 Inversion model of photothermal parameters of dispersed media. 
Considering the interference in the actual measurement, random noise was added to the measurement results of the simulation experiment. On the one hand, the real measurement results were simulated. On the other hand, the robustness of the CMAES algorithm in the radiation inverse problem of the photothermal parameter reconstruction of the medium was tested. The Gaussian error model used in this paper is as follows:(16)where I_{exact} is the accurate measurement result. I_{exp} is the actual measurement result. ς is the standard normally distributed random number with expectation 0 and variance 1. σ is the standard deviation, the calculation formula of σ is as follows:(17)where X_{ext} and X_{exact} respectively represent the estimated value and the true value of the parameters to be inverted. N_{var} is the dimension of the search space of the optimization algorithm, that is the number of parameters to be inverted.
The distribution of absorption coefficient and refractive index in the dispersion medium still has the property of gradient, which is shown as follows:(18)where a_{1}, b_{1}, a_{2}, b_{2} represent the coefficient of linear distribution, respectively. L represents the total thickness of the dispersion medium. This paper assumes that the thickness of the two sub media is the same. At the same time, they are the assumed parameter of the parameter inversion experiment and the parameter to be inverted of the reconstruction algorithm. x is the normalized spatial position of the thickness of the medium as the characteristic length.
3.3 Sensitivity analysis of absorption coefficient and refractive index of multilayer dispersive media
The outgoing radiation response signals at multiple angles at the right boundary under different wall temperatures are required as known information in the reconstruction algorithm of spectral characteristic parameters of dispersive media based on the CMAES algorithm. The ratio of wall temperature to medium temperature is φ:(19)
The output radiation response signal varies with the temperature ratio φ and the exit angle θ. To quantitatively analyze the relationship between the inversion parameters and the outgoing radiation response signal, the concept of sensitivity is introduced in this paper. The sensitivity describes the sensitivity between two sets of variables: when the sensitivity coefficient is small, the measured signal is less affected by the parameters to be inverted. Even without measurement error, it is difficult to reconstruct the inversion parameter distribution accurately according to the current measured signal. When the sensitivity coefficient is large, the measured signal varies greatly with the change of parameters to be inverted. When there is no measurement error, the precise inversion of parameters can be realized simply. However, when there is a small error in the measured signal, there will be a large deviation between the inversion result and the actual result.
When the sensitivity coefficient is large or small, the sensitivity of the inversion parameters can be adjusted by constructing an appropriate objective function to improve the computational efficiency of the reconstruction algorithm and the accuracy of the results. But it does not solve the problem from the root. In this paper, another idea is adopted to select the temperature ratio φ and exit angle θ corresponding to the output radiation response signal through the sensitivity coefficient. To improve the calculation efficiency, the accuracy of the results, and the robustness of the algorithm for the reconstruction of photothermal parameters of dispersive media.
The sensitivity coefficient can be calculated by taking the derivative of one variable to another variable. The sensitivity coefficient adopted in this paper can be calculated according to the following formula:(20)where is the sensitivity coefficient of the variable x_{i} to the measured signal. represents the radiation response signal in the φ direction along the angle between the right boundary and the normal direction of the surface when the temperature ratio is θ, x_{i} is the i parameter to be inversed, and its true value is x_{0}. δ is a small rate of change, in this paper δ is 0.5%.
In the twosided medium model, the sensitivity coefficient relative to the radiation response signal is shown in Figure 7 under the premise of different temperature ratios and different exit radiation angles of the absorption coefficient and refractive index of each submedium. Figures 7a and 7b correspond to the absorption coefficients of the twolayer media. The data show that the sensitivity of the absorption coefficient of the twolayer medium is greatly affected by the temperature ratio and the sensitivity coefficient decreases with the increase of the angle θ between the exit direction and the plane normal. The sensitivity of the refractive index of the twolayer medium is less affected by the temperature ratio. With the increase θ, the sensitivity coefficient of the refractive index of the first layer does not change at first and then decreases rapidly. The sensitivity coefficient of the second layer increases with the increase of θ. The trends of the two are not the same.
Fig. 7 Sensitivity analysis of parameters to be inverted. (a) a_{1}, (b) b_{1}, (c) a_{2}, (d) b_{2}. 
The essential reason for this difference is that the changes a_{2} b_{2} have different effects on the gradient refractive index distribution. It is assumed that the distribution of the gradient refractive index is linear, and a_{2} b_{2} correspond to the starting point and endpoint of the distribution function, respectively. Where a_{2} increases, corresponding to the increase of the starting point of the distribution function curve, the slope of the function (the highest coefficient of the function polynomial) decreases, that is the change of a_{2} is opposite to that of the function slope, so the changing trend of the sensitivity coefficients of the two is also opposite. The sensitivity coefficient of the linear distribution slope increases with the increase of the included angle θ when the temperature ratio is fixed, which corresponds to the decrease of the sensitivity coefficient a_{2} with the increase of θ in Figure 7c. The same analysis is carried out b_{2}.
The above analysis is extended to the highdimensional case, and the distribution of the corresponding gradient absorption coefficient and refractive index can be expressed by considering the relatively simple but representative threelayer medium model.
As in the twodimensional case, the sensitivity coefficient of the absorption coefficient of the threelayer medium has the same trend as the angle. The sensitivity coefficient decreases gradually at the same temperature ratio from left to right. There is a large gap between the first layer and the second layer because the first layer absorbs more radiation energy and plays a leading role. The changing trend of the sensitivity coefficient of the absorption coefficient of the threelayer medium is different. The first two layers decrease with the increase of θ, the third layer increases with the increase of θ. The reason is that the increase of the refractive index of the first two layers will decrease the quadratic coefficient corresponding to the distribution of the quadratic function of the refractive index and the sensitivity of the quadratic coefficient increases with the increase θ. Therefore, the refractive index of the first two layers shows the changing trend as shown in Figures 8b and 8d. By the same token, the variation trend of the sensitivity coefficient of the refractive index of the third layer medium in Figure 8f can also be explained.
Fig. 8 Sensitivity analysis of parameters to be inverted for threelayer media. (a) a_{1}, (b) a_{2}, (c) b_{1}, (d) b_{2}, (e) c_{1}, (f) c_{2}. 
According to the calculation results shown in Figure 8, the temperature ratio φ selected in this study is set as 2.3, and the value of the angle is the same as that when the gradient refractive index is distributed in a functional form.
As can be seen in Figure 8, the sensitivity coefficient of the refractive index of the medium increases sequentially from left to right and the sensitivity coefficient of the second layer is about that of the third layer, which leads to the application of the intelligent algorithm to retrieve the gradient refractive index of the medium. The first two layers may not be able to get more accurate results even when there is no measurement error, so the inversion is very difficult. In theory, the directional radiation response signal on the left can be added as the measurement information. Considering that in the direct problem model selected in this chapter, the left side is the opaque wall. Therefore, this paper only uses the directional radiation response signal on the right side as the measurement information in the gradient refractive index inversion experiment of threelayer media and inverts the distribution of three groups of absorption coefficient and refractive index without adding measurement error. Five times of independent inversion in each case. The selection of temperature ratio is the same as that in the inversion of photothermal characteristic parameters of twolayer media in inversion.
3.4 Joint inversion of absorption coefficient and refractive index for multilayer dispersive media
Firstly, the inversion of absorption coefficient and refractive index of doublelayer media under the condition of adding different measurement errors is analyzed. The calculated results are shown in Table 2. The variance of all inversion results is less than 10^{−4}, indicating that a total of 50 inversions included in this inversion experiment have achieved convergence.
Inversion results of absorption coefficient and refractive index of doublelayer media model.
The data show that the average value of the inversion results is close to the truth value of the parameters to be inverted and the error is small. It indicates that the reconstruction algorithm of the photothermal characteristics parameters of the dispersion medium constructed in this section for the bilayer medium model can calculate the true value of the inversion under the condition of the existence of errors.
The RMSE statistics of the results can obtain better measurement results even with large measurement errors. However, there are fluctuations in the RMSE data of the refractive index inversion results when the SNR of the added error is equal to 30 dB and 40 dB, which do not increase with the increase of the SNR. Moreover, the comparison of the data shows that the RMSE values of the two cases are very close. This paper believes that the reason for this phenomenon is the randomness of the inversion algorithm, so the results of this calculation are contingent.
The aptamer function and RMSE vary with SNR as shown in Figure 9. The value of aptamer function in Figure 9a increases with the decrease of signaltonoise ratio, the maximum value is less than 8 × 10^{−4}. The fluctuation described above can be seen in the red line in Figure 9b. Corresponding to the axis on the right, it can be seen that the maximum value of the axis is only 0.01. The data in the figure further verify the reliability of the inversion algorithm.
Fig. 9 The fitting function and RMSE change with SNR. (a) Bodyfitting function, (b) RMSE. 
The inversion of absorption coefficient and refractive index of threelayer media is analyzed below. The inversion results are shown in Table 3. The data show that good measurement results have been obtained from the three sets of inversion. Comparing the measurement results of the absorption coefficient and the refractive index, it is found that the standard deviation of the threelayer absorption coefficient reconstruction results is of the same or similar magnitude. But the standard deviation of the refractive index reconstruction results decreases gradually from left to right. In the sensitivity analysis of the threelayer medium model, the sensitivity of the refractive index increases gradually with the order of the medium from left to right. There is a large gap between the sensitivity coefficients of the two adjacent layers. This difference in sensitivity will lead to poor accuracy of the inversion results when there is a large error. Based on sacrificing the calculation time, the accuracy of reconstruction can be improved by increasing the directional radiation response signal under different temperature ratios.
Inversion results of absorption coefficient and refractive index of threelayer media model.
3.5 Inversion analysis of temperature of hightemperature wall obscured by dispersive media
Based on the previous reconstruction algorithm of photothermal characteristic parameter distribution of dispersion medium, this paper constructs the reconstruction algorithm of wall temperature under the condition of dispersion medium shadowing. The direct problem model refers to Figure 5. In this paper, it is assumed that the distribution of photothermal parameters of the medium is in the discrete form shown in equation (18). Taking the directional radiation response signal at the right wall as the measurement data input based on knowing the emissivity of the left wall. The wall temperature is inversed to realize the wall temperature measurement under the shadowing condition.
The following inversion experiments are carried out on the proposed wall temperature inversion algorithm. The flow chart of the experiment is shown in Figure 10. In the figure, ① indicates that the inversion results of photothermal parameters of the dispersion medium are brought back to the direct problem model. ② Indicates that the error calculation is carried out for the inversion results of the wall temperature.
Fig. 10 Experimental flow chart of wall temperature inversion. 
In the inversion experiment, the calculation result corresponding to SNR = 30 dB in Table 2 is used as the calculation result of the first step, the direct problem model is used to invert the wall temperature. The inversion was carried out under two conditions: no measurement error and added measurement error. The algorithm was independently run 10 times in each case. We statistically analyzed the average, standard deviation, and average measurement error of the calculated results.
The calculation results are shown in Table 4. The data show that in the wall temperature inversion method proposed in this section, the calculation error tolerance for the first step is large. Even if the first step of the direct problem is based on the measured signal with a high signaltonoise ratio, the inversion of temperature in the second step can still obtain better results.
Experimental calculation results of wall temperature inversion.
4 Results and discussion
4.1 Measuring result
The system consists of five parts: black body furnace, Fourier transform infrared spectrometer, electric control rotary table, control, and data acquisition system and test components. Using the experimental platform to measure the radiation signal of the translucent medium. Figure 11 shows the schematic diagram of the measurement system.
Fig. 11 Measurement system for radiation signal of a translucent body. 
In the direct problem model, the radiation signal I is a function of the wall temperature T_{w}, the refractive index n_{1} of the first sample, the extinction coefficient κ_{a1} of the first sample, the refractive index n_{2} of the second sample, and the extinction coefficient κ_{a2} of the second sample. The size and material of the two plates are the same and the same amount of smoke is injected into the two dishes at the same time during the experiment, so it is considered that the optical properties of the two samples are completely similar.
First of all, the n_{1}, n_{2}, κ_{a1} and κ_{a2} of the twolayer samples are inversed, and different angles are measured experimentally (the measuring angle refers to the angle between the normal of the light surface of the sample and the incident light source, the angle of this experiment is from 0°to 60°, every 5° is a measuring point) and different wall temperatures (the wall temperature is the temperature of the blackbody furnace, the temperature range of this experiment is from 600 K to 1700 K, every 100 K is a measuring point). The spectral radiation signal intensity at the wave number 399.1927–4001.5688 cm^{−1}.
Planck integral is introduced to process the measured radiation signals at different wavelengths and the radiation characteristic parameters of the colorimeter in this paper. The formula is as equations (21) and (22). After deducting background noise and Planck integral processing by the doubletemperature method, the radiation intensity of different wall temperatures and different sample rotation angles was obtained, as shown in Figure 12:(21) (22)
Fig. 12 Spectral signal intensity under different wall temperatures and sample rotation angles. 
The initial search range κ_{a1} κ_{a2} is set to 0 to 1 when the algorithm is used for inversion. The refractive index of smoke in the sample is similar to that of air, so the initial search range n_{1} and n_{2} is set to 1.000 to 1.005. The inversion result is: n_{1} = 1, n_{2} = 1, κ_{a1} = 0.192, κ_{a2} = 0.187. There is little difference between κ_{a1} and κ_{a2}, n_{1} and n_{2}, which confirms the assumption that the optical properties of the two samples are similar.
Then, n_{1}, n_{2}, κ_{a1} and κ_{a2} obtained in the first step are used as known conditions to invert the wall temperature T_{w} into the direct problem model. In this inversion, the radiation signal measured when the temperature T_{w} of the blackbody furnace is 1400 K, 1500 K, 1600 K, and 1700 K is inverted by a random algorithm. The initial search temperature range set by the random algorithm is 100–2000 K. The wall temperature T_{w0} obtained by inversion is shown in Table 5.
Inversion results of T_{w}.
The results show that the maximum error is less than 1% when the temperature is higher. As the true value of the wall temperature gradually decreases, the inversion accuracy gradually increases. The maximum is less than 3%. It proves the effectiveness and reliability of the wall temperature inversion algorithm based on the intelligent algorithm proposed in this paper under the condition of dispersion medium shielding.
The spectral radiation intensity of the sample at 5 μm, 10.53 μm and 16 μm at different rotation angles and different wall temperatures was measured on the experimental bench to verify the accuracy of the data processed by the Planck integral method. The inversion of k_{1} and k_{2} was carried out by the algorithm. Then, k_{1} and k_{2} obtained by inversion are used as known quantities. The spectral radiation intensities measured at different wavelengths at T_{1} = 1400 K, T_{2} = 1500 K, T_{3} = 1600 K and T_{4} = 1700 K are combined to invert the wall temperature of the blackbody furnace by the algorithm. Table 6 shows the inversion results of experimental data obtained by the Planck integral method and the inversion results without Planck integral.
Inversion results using different spectral radiation intensities as known information.
It can be seen from the above table that there is little difference between the inversion results obtained by the Planck integral method and those without Planck integral, which proves the correctness of using Planck integral to calculate the measurement information as the input of the inversion algorithm.
The tubular furnace is widely used for heating equipment in petrochemical plants. A tube furnace involves a variety of hydrocarbon chemical reactions, as a hydrocarbon energy system. The effective measures to improve the thermal efficiency of the tubular furnace are to use refractory materials with low thermal conductivity and to reduce the temperature of the external wall of the furnace to reduce heat loss. However, the improvement of the thermal efficiency of a tubular furnace is restricted by many factors. Combustion temperature measurement technology provides a direction for improving the thermal efficiency of the tubular furnace.
The calculated results show that the designed temperature reconstruction algorithm can obtain more accurate results when the temperature is less than 1700 K. Take the tubular furnace produced by a certain enterprise as an example. The highest temperature in the furnace can reach 1200 °C, less than 1700 K. Therefore, the algorithm can be used to measure the wall temperature of the tubular furnace and provide direction for its improvement.
4.2 Uncertainty analysis
The errors of κ_{a1} and κ_{a2} inversion in this experiment mainly include the measurement error of the Fourier transforms infrared spectrometer, the temperature error of the blackbody furnace and the angle error of the electric displacement table. The uncertainty of each factor is analyzed below.
 1.
Uncertainty of Fourier transform infrared spectrometer
The uncertainty of the radiation intensity measured by the spectrometer can be divided into two types: The first type is the measurement error caused by the background noise in the measurement environment and the noise in the measurement process. It has been eliminated by the doubletemperature method above. The second type is the measurement error caused by the nonlinear spectral response of the instrument, which is the main source of the measurement error of the Fourier transform infrared spectrometer. The uncertainty of the nonlinear spectral response of the FTIR6100 spectrometer in this experiment is 0.32%.
 2.
The uncertainty of blackbody furnace temperature
The uncertainty of the temperature of the blackbody furnace is mainly caused by the fluctuation of the heating circuit under PID control. The temperature of the resistance wire will fluctuate up and down near the target temperature with the increase of time and the amplitude of fluctuation gradually decreases with the increase of time in the process of reaching the target temperature. The hightemperature blackbody furnace used in this experiment is the RT1500 type blackbody furnace produced by LAND Company in the UK. The effective emissivity of the blackbody furnace is 0.99 ± 0.001. The error range of its temperature heating is ±0.1 °C. The uncertainty is only 4% when the temperature is 25 °C (at room temperature). The uncertainty decreases gradually with the increase of temperature.
 3.
The uncertainty of the electric displacement table
The angle error of the electric displacement table refers to the difference between the actual rotation angle of the turntable and the set angle. The precision of the Y100SC06 electronic control rotary table is 1 ± 0.005°, and the uncertainty is 0.5%.
The above three uncertainties are normally distributed and do not influence each other. Therefore, the overall uncertainty synthesis formula is as follows:(23)
The formula, U_{A} represents the uncertainty of the Fourier transform infrared spectrometer, U_{B} represents the uncertainty of blackbody furnace temperature, and U_{C} represents the uncertainty of rotation angle of the electronically controlled rotary table. The uncertainty of the system calculated by the above formula is 4.04%. The value is less than 5%, which meets the requirements of extinction coefficient measurement accuracy in the field of most numbers.
5 Conclusion
In this paper, the distribution reconstruction of nonuniform photothermal characteristic parameters of dispersive media based on adaptive covariance matrix evolutionary strategy algorithm is studied. The specific research contents are as follows:

Based on the CMAES algorithm, the reconstruction program of gradient absorption coefficient and refractive index of functional distribution is established. To obtain a better combination of measurement information, the temperature ratio was defined. The sensitivity of distribution function coefficients was analyzed under different temperature ratios, outgoing radiation, and different surface normal angles. According to the graph line of sensitivity coefficient, the temperature ratio and the included Angle are selected. The reconstruction algorithm is verified and the results show that the algorithm can obtain relatively stable measurement results even with large measurement errors.

The more general distribution of photothermal characteristic parameters is considered. The form of the discrete distribution gradient absorption coefficient and refractive index distribution function is unknown during reconstruction. For the twolayer medium model, the sensitivity of absorption coefficient and refractive index is analyzed and the rationality of the image drawn in this paper is explained. The sensitivity graph of absorption coefficient and refractive index of the threelayer medium model is analyzed. The proposed reconstruction algorithm is verified by a twolayer media model and threelayer media model respectively, and the results are good.

An inversion algorithm for wall temperature reconstruction under the obscuring condition of dispersive media is proposed. Under the condition of introducing different measurement errors, the effectiveness of the algorithm is verified by experiments. The results show that the performance of the proposed inversion algorithm is good. The proposed algorithm is robust even at higher temperatures. Referring to the internal temperature data of tubular furnaces in the market, the proposed algorithm can be used to provide direction for improving the thermal efficiency of tubular furnaces.
Funding
This research was supported by National Key R&D Program of China: 2019YFF0302202.
References
 Yang S.M., Tao W.Q. (2006) Heat transfer, 4th edn, Higher Education Press, Beijing. [Google Scholar]
 Wei L.Y. (2020) Photothermal information simulation and parameter identification in gradedindex media, Harbin Institute of Technology, Harbin. https://doi.org/10.27061/d.cnki.ghgdu.2020.001579. [Google Scholar]
 Liu L.H. (2006) Finite volume method for radiation heat transfer in graded index medium, J. Thermophys. Heat Transf. 20, 1, 59–66. [CrossRef] [MathSciNet] [Google Scholar]
 Huang Y., Xia X.L., Tan H.P. (2002) Temperature field of radiative equilibrium in a semitransparent slab with a linear refractive index and gray walls, J. Quant. Spectrosc. Radiat. Transfer 74, 2, 249–261. [CrossRef] [Google Scholar]
 Huang Y., Xia X.L., Tan H.P. (2003) Comparison of two methods for solving radiative heat transfer in a gradient index semitransparent slab, Numer. Heat Transf. Fundam. 44, 1, 83–99. [CrossRef] [Google Scholar]
 Yong H., Liang X.G., Xia X.L. (2005) Monte Carlo simulation of radiative transfer in scattering, emitting, absorbing slab with gradient index, J. Quant. Spectrosc. Radiat. Transfer 92, 1, 111–120. [CrossRef] [Google Scholar]
 Ren N., Levoy M., Bredif M., Duval G., Hanrahan P. (2005) Light field photography with a handheld plenoptic camera. Stanford University CTSR. [Google Scholar]
 Zhao Z.Y., Nie J.S. (2010) Research on thermal lens effect of grooved laser crystal, ElectroOptic Technol. Appl. 05, 43–45. [Google Scholar]
 Xia X.L., Tan H.P., Ren D.P. (2005) Research status of gradedindexthermal radiation and its application prospect in spacecraft technology, Spacecr. Eng. 3, 26–33. [Google Scholar]
 Zhang X.P. (2014) Optical analysis of mirage phenomenon, Friends Phys. 30, 10, 18–19. [Google Scholar]
 Zhu J.R., Chen YW, Zhang J.Y., Ni H. (2014) The effect of light deflection on interference pattern of flame temperature field with measured by laser, Lasernal 2, 52–53. [Google Scholar]
 Wang Z.C. (2013) A study on radiative heat transfer within 1D gradedindex medium by the Dresor method, Huazhong University of Science and Technology, Wuhan. [Google Scholar]
 Huang Y. (2002) Study on thermal radiation transfer in gradedindex translucent medium, PhD Thesis in Engineering, Harbin Institute of Technology, Harbin, 12 p. [Google Scholar]
 Chandrasekhar S. (1950) Radiative transfer, Oxford University Press, New York. [Google Scholar]
 Love T.J., Grosh R.J. (1965) Radiative heat transfer in absorbing, emitting, and scattering media, J. Heat Transf. Am. Soc. Mech. Eng. 87, 2, 161–166. [Google Scholar]
 Lemonnier D., Le Dez V. (2002) Discrete ordinates solution of radiative transfer across a slab with variable refractive index, J. Quant. Spectrosc. Radiat. Transf. 73, 2, 195–204. Elsevier. [CrossRef] [Google Scholar]
 Raithby G.D., Chui E.H. (1990) A finitevolume method for predicting a radiant heat transfer in enclosures with participating media, J. Heat Transf. Am. Soc. Mech. Eng. 112, 2, 415–423. [Google Scholar]
 Chui E.H., Raithby G.D. (1992) Implicit solution scheme to improve convergence rate in radiative transfer problems, Numer. Heat Transf. Part B Fundam. 22, 3, 251–272. Taylor & Francis. [CrossRef] [Google Scholar]
 Liu L.H. (2004) Finite element simulation of radiative heat transfer in absorbing and scattering media, J. Thermophys. Heat Transf. 18, 4, 555–557. [CrossRef] [Google Scholar]
 Zhang L., Zhao J.M., Liu L.H. (2009) Finite element method for modeling radiative transfer in semitransparent graded index cylindrical medium, J. Quant. Spectrosc. Radiat. Transf. 110, 13, 1085–1096. Elsevier. [CrossRef] [Google Scholar]
 Kim A.D., Ishimaru A. (1999) A Chebyshev spectral method for radiative transfer equations applied to electromagnetic wave propagation and scattering in a discrete random medium, J. Comput. Phys. 152, 1, 264–280. Elsevier. [Google Scholar]
 Kim A.D., Moscoso M. (2002) Chebyshev spectral methods for radiative transfer, SIAM J. Sci. Comput. 23, 6, 2074–2094. SIAM. [CrossRef] [MathSciNet] [Google Scholar]
 Wei L.Y., Qi H., Niu Z.T., Wen S., Ren Y.T. (2021) Inverse heat transfer analysis to determine the temperature or phase changedependent refractive index of semitransparent materials, Inverse Probl. Sci. Eng. 29, 4, 586–608. [CrossRef] [MathSciNet] [Google Scholar]
 Zhao J.M., Liu L.H. (2007) Solution of radiative heat transfer in graded index media by least square spectral element method, Int. J. Heat Mass Transf. 50, 13, 2634–2642. Elsevier. [CrossRef] [Google Scholar]
 Zhao J.M., Liu L.H. (2008) Discontinuous spectral element approach for solving transient radiative transfer equation, J. Thermophys. Heat Transf. 22, 1, 20–28. [CrossRef] [Google Scholar]
 Zhao J.M., Liu L.H. (2008) Spectral element method with adaptive artificial diffusion for solving the radiative transfer equation, Numer. Heat Transf. Part B Fundam. 53, 6, 536–554. Taylor & Francis. [CrossRef] [Google Scholar]
 Wei L.Y., Sadaf A.I., Wen S. (2018) Simultaneous measurement of spacedependent refractive index and absorption coefficient based on a laser irradiation technique, Measur. Sci. Technol. 29, 10, 104006 (13pp). [Google Scholar]
 Hansen N. (2006) The CMA evolution strategy: A comparing review, Stud Fuzz. Soft Comput. 192, 75–102. [CrossRef] [Google Scholar]
 Hansen N. (2005) The CMA evolution strategy: A tutorial. https://lanl.arxiv.org/abs/1604.00772. [Google Scholar]
 Hansen N., Ostermeier A. (2001) Completely derandomized selfadaptation in evolution strategies, Evol. Comput. 9, 2, 159–195. [CrossRef] [PubMed] [Google Scholar]
 Ostermeier A., Gawelczyk A., Hansen N. (1994) A derandomized approach to selfadaptation of evolution strategies, Evol. Comput. 2, 4, 369–380. [CrossRef] [Google Scholar]
All Tables
Inversion results of absorption coefficient and refractive index of doublelayer media model.
Inversion results of absorption coefficient and refractive index of threelayer media model.
Inversion results using different spectral radiation intensities as known information.
All Figures
Fig. 1 The physical model of radiative heat transfer in 1D participating media with GRI. 

In the text 
Fig. 2 Onedimensional translucent plate model. 

In the text 
Fig. 3 Flow chart of inversion experiment. 

In the text 
Fig. 4 Fitness function values and calculation time statistics of running tests. (a) Fitter function value, (b) computation time. 

In the text 
Fig. 5 A physical model of radiation heat transfer in onedimensional inhomogeneous disseminated media considering radiation from one sidewall. 

In the text 
Fig. 6 Inversion model of photothermal parameters of dispersed media. 

In the text 
Fig. 7 Sensitivity analysis of parameters to be inverted. (a) a_{1}, (b) b_{1}, (c) a_{2}, (d) b_{2}. 

In the text 
Fig. 8 Sensitivity analysis of parameters to be inverted for threelayer media. (a) a_{1}, (b) a_{2}, (c) b_{1}, (d) b_{2}, (e) c_{1}, (f) c_{2}. 

In the text 
Fig. 9 The fitting function and RMSE change with SNR. (a) Bodyfitting function, (b) RMSE. 

In the text 
Fig. 10 Experimental flow chart of wall temperature inversion. 

In the text 
Fig. 11 Measurement system for radiation signal of a translucent body. 

In the text 
Fig. 12 Spectral signal intensity under different wall temperatures and sample rotation angles. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.