Optimal Experimental Design for Inverse Identiﬁcation of Conductive and Radiative Properties of Participating Medium

: The conductive and radiative properties of participating medium can be estimated by solving an inverse problem that combines transient temperature measurements and a forward model to predict the coupled conductive and radiative heat transfer. The procedure, as well as the estimates of parameters, are not only affected by the measurement noise that intrinsically exists in the experiment, but are also inﬂuenced by the known model parameters that are used as necessary inputs to solve the forward problem. In the present study, a stochastic Cram é r–Rao bound (sCRB)-based error analysis method was employed for estimation of the errors of the retrieved conductive and radiative properties in an inverse identiﬁcation process. The method took into account both the uncertainties of the experimental noise and the uncertain model parameter errors. Moreover, we applied the method to design the optimal location of the temperature probe, and to predict the relative error contribution of different error sources for combined conductive and radiative inverse problems. The results show that the proposed methodology is able to determine, a priori, the errors of the retrieved parameters, and that the accuracy of the retrieved parameters can be improved by setting the temperature probe at an optimal sensor position.


Introduction
Participating medium is widely presented in many engineering fields, such as aerospace engineering, energy and power systems, and information communications. The conductive and radiative properties of participating medium can be determined from transient temperature measurements by solving an inverse problem [1][2][3]. This procedure consists of comparing the measured temperatures to the responses predicted from combined conductive and radiative heat transfer simulation [4][5][6][7]. A sequence of computations is performed, and the property values are adjusted until the predictions match well with the measurements.
For the ideal case, the experimental measurements are noise-free and the predictions perfectly reflect the reality; consequently, the conductive and radiative properties may be precisely recovered. However, neither the measurements nor the predictions are strictly accurate. The measurements are distorted by measurement noise, and the predictions, simulated via combined conduction and radiation, may exhibit variations due to the fact that: (1) the solution method may simplify the complex coupled heat transfer problem, and thus, yield inaccurate temperature responses; and (2) some of the model parameters (such as geometry parameters, density, the specific heat of the material, and the boundary conditions) used in the forward problem solution are not precisely known. As the inverse problem is always ill-posed, a small deviation in the measured or predicted temperatures may lead to a considerable deviation in the recovered properties. Therefore, there is a crucial need to investigate the uncertainties of the recovered properties in the inverse scheme, and to make efforts to improve the accuracy of the retrieved parameters.
Several studies were performed to estimate the uncertainties of the recovered properties. Lazard et al. [8] retrieved the thermal diffusivity, the Plank number, and a global radiative transfer coefficient (defined from the absorption coefficient, the scattering coefficient, and the slab thickness) of an absorbing and isotropic scattering slab from transient temperature measurement, and the authors estimated the variance of the retrieved parameters. Zhao et al. [9,10] performed transient temperature measurements of fibrous insulation, and retrieved the conductive and radiative properties by solving an inverse problem; the uncertainties of the retrieved parameters were evaluated from the standard deviation of the measured temperature response. Several similar studies, which aimed to estimate the thermophysical properties of anisotropic composite [11], the thermal conductivity and heat capacity of an orthotropic medium [12], and the conductive and radiative properties of participating medium [13,14], were also performed. The above-mentioned studies considered only the experimental noise, while the uncertainties that might have existed in the known model parameters of heat transfer models were not taken into account, i.e., the predictions were assumed to be strictly accurate. Only a few research studies considered both the experimental noise and the uncertainties of model parameters, and the uncertainties of the retrieved properties were estimated using the Cramér-Rao lower bound (CRB)-based method [15][16][17][18][19][20][21][22]. These works, relative to inverse heat transfer problems, mainly focused on retrieving the thermal conductivity, thermal resistance, and heat transfer coefficient by solving inverse heat conduction problems [15][16][17]. The other studies mainly investigated the uncertainty estimation and the selection of measurement modalities for the retrieval of the magnetic material properties of electromagnetic devices (EMD) [20][21][22]. To the best of the authors' knowledge, there is no uncertainty analysis research reported for inverse conduction-radiation problems that considers both the experimental noise and the uncertainties of the model parameters.
As for the strategies for improving the accuracy of the retrieved parameters, the first approach involves reducing the errors related to the inverse identification solution, i.e.,: (1) performing accurate experimental measurements, and thereby offering perfect measured temperature responses; (2) acquiring accurate model parameters (generally measured from other experiments) before solving the inverse problem; and (3) employing an accurate method for the solution of the forward problem. However, this strategy is usually ineffective due to the limitations of experimental equipment, such that the accuracy of the transient temperature and model parameter measurements are usually difficult to improve. Another means of improving the identification accuracy is to place the sensors at optimal positions; this involves the solution of an optimal experimental design problem. In general, the optimal sensor location is such that the sensitivity of the temperature responses at the optimal positions for the parameters to be retrieved should be as large as possible. As for most transient problems, the sensitivity is a function of time, and the location of maximal sensitivity at any time duration may not be unique, and thus, the ideal position of a single sensor is not unequivocally defined; instead, the sensor should be located at positions that give the best integrated sensitivity over the entire experimental time duration. In addition, the optimal locations should be designed to minimize the noise effects, i.e., the measured temperature responses at the optimal sensor positions should be as accurate as possible, and the corresponding predictions should be less sensitive to the uncertainties of the known model parameters when solving the forward problem. It is not obvious that the maximization of integrated sensitivity and the minimization of noise effects lead to the same sensor location; therefore, the optimal design of the sensor location is comprehensively affected by the aforementioned two factors. This paper presents a stochastic Cramér-Rao bound (sCRB)-based error analysis method for estimating the uncertainties of conductive and radiative properties retrieved Energies 2021, 14, 6593 3 of 16 from transient temperature measurements by solving an inverse conductive and radiative heat transfer problem. The measurement noise and the uncertainties of known model parameters are both taken into account in the analysis, whereas the solution error that occurs as a result of the method used to solve forward problem is neglected in the present study. Furthermore, the optimal temperature sensor positions for inverse transient conductive and radiative heat transfer problems are designed to improve the accuracy of the retrieved properties on the basis of the CRB-based error analysis method. Several examples are given to illustrate the error analysis method and to show the superiority of the designed optimal sensor positions. The remainder of the present manuscript is organized as follows: Section 2 presents a combined conduction and radiation model, an inverse identification method, and the CRB-based uncertainty analysis method. Several examples, as well as the corresponding discussions, are presented in Section 3. Conclusions are drawn at the end of this manuscript.

Combined Conductive and Radiative Heat Transfer in Participating Medium
Transient coupled conductive and radiative heat transfer, in an absorbing and isotropic scattering gray solid slab with a thickness of L, were considered. The physical model of the slab, as well as the associated coordinate system, are shown in Figure 1. As the geometry considered was a solid slab, convection was not considered in the present study. In addition, the geometry can be three-dimensional but only one direction is relevant; thus, only 1-D combined conductive and radiative heat transfer was investigated. The boundaries of the slab were assumed to be diffuse and gray opaque, with an emissivity of ε 0 for x = 0, and ε L for x = L, and the temperatures of the two walls were fixed at T L and T H , respectively. The extinction coefficient β, the scattering albedo ω, the thermal conductivity k c , the density ρ, and the specific heat c p of the slab were assumed to be constant in the present study. heat transfer problem. The measurement noise and the uncertainties of known model parameters are both taken into account in the analysis, whereas the solution error that occurs as a result of the method used to solve forward problem is neglected in the present study. Furthermore, the optimal temperature sensor positions for inverse transient conductive and radiative heat transfer problems are designed to improve the accuracy of the retrieved properties on the basis of the CRB-based error analysis method. Several examples are given to illustrate the error analysis method and to show the superiority of the designed optimal sensor positions. The remainder of the present manuscript is organized as follows: Section 2 presents a combined conduction and radiation model, an inverse identification method, and the CRB-based uncertainty analysis method. Several examples, as well as the corresponding discussions, are presented in Section 3. Conclusions are drawn at the end of this manuscript.

Combined Conductive and Radiative Heat Transfer in Participating Medium
Transient coupled conductive and radiative heat transfer, in an absorbing and isotropic scattering gray solid slab with a thickness of L, were considered. The physical model of the slab, as well as the associated coordinate system, are shown in Figure 1. As the geometry considered was a solid slab, convection was not considered in the present study. In addition, the geometry can be three-dimensional but only one direction is relevant; thus, only 1-D combined conductive and radiative heat transfer was investigated. The boundaries of the slab were assumed to be diffuse and gray opaque, with an emissivity of ε0 for x = 0, and εL for x = L, and the temperatures of the two walls were fixed at TL and TH, respectively. The extinction coefficient β, the scattering albedo ω, the thermal conductivity kc, the density ρ, and the specific heat cp of the slab were assumed to be constant in the present study. The energy conservation equation for the slab can be written as [23,24] where T(x,t) denotes the temperature at position x and time t, and , r q x t is the radiative heat flux. The gradient of the radiative heat flux can be calculated from [23,24]. The energy conservation equation for the slab can be written as [23,24] ρc p ∂T ∂t where T(x,t) denotes the temperature at position x and time t, and q r (x, t) is the radiative heat flux. The gradient of the radiative heat flux can be calculated from [23,24].
is the blackbody radiation intensity, and G is the fluence rate defined as G(r) = 4π I(ŝ,r)dΩ. The intensity I(ŝ,r) at locationr in directionŝ is governed by the radiative transfer equation (RTE), written as [23,24] where κ is the absorption coefficient and n is the refractive index, while σ s and Φ(ŝ ,ŝ) are the scattering coefficient and phase function of scattering, respectively. Corresponding to the physical model shown in Figure 1, the boundary and initial conditions can be written as where T H and T L are the temperatures at the walls of x = 0 and x = L, respectively, while T 0 is the initial temperature. The transient temperature T(x s , t) at sensor position x s can be predicted by solving Equations (1)-(4).

Inverse Method
In order to maintain the generality of the method, we assumed that the unknown conductive and radiative properties to be retrieved were labeled as u ∈ R N P ×1 , where N p is the number of unknown parameters. The model parameters were assumed to be b ∈ R N q ×1 , where N q is the number of model parameters. The predicted transient temperature at location x s was T u, b ∈ R N S ×N t , where u is the retrieved value of u, and b is the measured value of parameter b, N t is the number of sampling points, and N S is the number of sensor positions. The transient temperature history measured in the 'experiment' was expressed as W ∈ R N S ×N t .
The inverse problem was defined as an optimization problem of finding the parameter vector u, for which the transient temperature history T u, b at location x s predicted from combined conduction and radiation was closest to the experimental data W; thus, the parameter vector u could be determined by minimizing an objective function defined as Thus, The genetic algorithm (GA), which is widely used for complex, ill-posed problems [25][26][27], was employed to solve the inverse identification problem. Figure 2 shows the block diagram of the GA-based inverse method used to determine parameter vector u.

Uncertainty Estimation and Design of Experiment
In the present study, a mathematical technique based on the stochastic Cramér-Rao lower bound (sCRB) is presented; the method aimed to take the measurement noise and the model parameter uncertainties of combined conduction and radiation into account for an a priori uncertainty estimation of the retrieved properties. Assuming that u* is the actual value of parameter vector u, the actual measurement W can be expressed as is the total error vector, and the total error vector at time tk can be expressed as where exp e and pred e are the error vectors due to measurement noise and modeling uncertainties, respectively. The measurement error exp e is composed of systematic and random components, as the state-of-the-art techniques and devices used for temperature measurement provide a rather low level of systematic error, and the reproducible nature of the systematic error makes it possible to estimate the bias on the measured data by means of a calibration procedure; this manuscript restricts discussions that the measurements contain only the random component of uncertainties, and the random error is assumed to be Gaussian while distributed with a mean of zero and a variance of 2 exp,k σ . The modeling error, pred e , can also be divided into two parts: the modeling error due to the use of inaccurate model parameter vector b, and the modeling error due to the use of inaccurate physical models

Uncertainty Estimation and Design of Experiment
In the present study, a mathematical technique based on the stochastic Cramér-Rao lower bound (sCRB) is presented; the method aimed to take the measurement noise and the model parameter uncertainties of combined conduction and radiation into account for an a priori uncertainty estimation of the retrieved properties. Assuming that u* is the actual value of parameter vector u, the actual measurement W can be expressed as where e tot ∈ R N S ×N t is the total error vector, and the total error vector at time t k can be expressed as The total error vector, e tot , contains two components, i.e., e tot = e exp + e pred , where e exp and e pred are the error vectors due to measurement noise and modeling uncertainties, respectively. The measurement error e exp is composed of systematic and random components, as the state-of-the-art techniques and devices used for temperature measurement provide a rather low level of systematic error, and the reproducible nature of the systematic error makes it possible to estimate the bias on the measured data by means of a calibration procedure; this manuscript restricts discussions that the measurements contain only the random component of uncertainties, and the random error is assumed to be Gaussian while distributed with a mean of zero and a variance of σ 2 exp,k . The modeling error, e pred , can also be divided into two parts: the modeling error due to the use of inaccurate model parameter vector b, and the modeling error due to the use of inaccurate physical models (such as simplification of the physical models, or the use of inaccurate numerical methods). In this study, we assumed that the physical model was perfect; thus, the modeling error was affected only by the inaccurate model parameters.
The Cramér-Rao inequality theorem states that the covariance matrix of the deviation between the true and the estimated parameters is bounded from below by the inverse of the Fisher information matrix M [15][16][17] where, the Fisher information matrix can be calculated from where M is a matrix with N p × N p dimensions, and ln L( W| u) is the log-likelihood of W given the parameter vector u; the likelihood of the data is normally distributed and is given by [15][16][17] where V k is the total error variance, expressed as The total error includes contributions of both the measurement noise and the modeling error k is the contribution of the uncertain model parameters, where G ∈ R N q ×N q is the covariance matrix of the uncertain modal parameter vector b, and Θ k is the sensitivity matrix of the temperature prediction T with respect to the uncertain parameter vector b; this can be expressed as Equation (11) can be rewritten as The first term of the right side is constant, thus The effect of the trace term is very small and can be neglected [17]; thus, the Fisher information matrix can be approximated by The lower bound for the variances of the parameters to be retrieved can be estimated as The σ 2 u i ,LB values could be used to qualitatively evaluate the retrieved results as well as the inverse identification models, and thus, could be employed in the method used to design the experiment. For inverse problems with only one parameter to be retrieved, the Fisher information matrix M can be reduced to a scalar M, σ 2 u,LB = 1/M. The algorithm for determining the optimal sensor position for inverse conductive and radiative heat transfer is shown in Figure 3 as follows: Step 1: Identify the mean value b of b and the corresponding covariance matrix G; Step 2: Identify possible sensor positions, and chose an initial sensor position; Step 3: Solve the forward problem, predict T u, b and the corresponding sensitivity Θ, then estimate the experimental error σ 2 exp ; Step 4: Estimate σ 2 u,LB for the retrieved parameter u; Step 5: Update the sensor position and go to step 3, then estimate σ 2 u,LB for all sensor positions; Step 6: Evaluate the different sensor positions and find the optimal sensor position.

Results and Discussion
We 'simulated' the measurements by using the output of the forward model with the actual values of the unknown parameters to be retrieved, and the measurements were corrupted by Gaussian noise with a mean and standard deviation of zero. In this way, we were able to perform numerical experiments to illustrate the uncertainty of the method of analysis, and to show the designed optimal sensor position.

Results and Discussion
We 'simulated' the measurements by using the output of the forward model with the actual values of the unknown parameters to be retrieved, and the measurements were corrupted by Gaussian noise with a mean and standard deviation of zero. In this way, we were able to perform numerical experiments to illustrate the uncertainty of the method of analysis, and to show the designed optimal sensor position.

Identification of Conductive Thermal Conductivity: The Optimal Experimental Design
For the first example, we considered an absorbing and isotropic scattering gray slab, as shown in Figure 1. The slab had volumetric heat capacity c v = 4×10 5 J/(m 3 ·K), refractive index n = 1.0, extinction coefficient β = 2000 m −1 , and scattering albedo ω = 0.8; furthermore, the thickness of the slab was L = 0.02 m, and the boundary emissivities for x = 0 and for x = L were ε 0 = ε L = 0.8. The initial temperature of the slab was T 0 = 300 K; the temperature of the two slab walls were fixed at T H = 1000 K for x = 0, and T L = 300 K for x = L, respectively. The conductive thermal conductivity k c was assumed to be unknown and needed to be retrieved, and the actual value of k c was 0.02 W/(m·K). The 'measured temperature' was simulated by adding Gaussian noise to the solution of combined conduction and radiation, and various noise levels, γ TS % , were considered in the present study. The standard deviation, σ exp,k , of the 'measured temperature', at time t k , could be calculated from where ξ is a normal distribution number with a mean of zero and a standard deviation of unity. We also considered that the boundary temperature T H at x = L was not accurate, and its value was distributed about the mean value of T H = 1000 K with a normally distributed noise γ TH % ; similarly, the standard deviation σ TH of T H could be calculated from σ TH = γ TH % · T H · ξ. Here, the considered time duration of the 'experiment' was t S = 200 s, and the sampling increment of time was ∆t = 2 s. Figure 4 shows the estimated lower bound, σ k c ,LB , for the standard deviation of the conductive thermal conductivity, k c , with respect to various measurement noise γ TS % and boundary temperature error γ TH % values; the values considered for γ TS % and γ TH % ranged from 1% to 9%, with an increment of 2%, and the temperature sensor was located at x s = L/2.

Identification of Conductive Thermal Conductivity: The Optimal Experimental Design
For the first example, we considered an absorbing and isotropic scattering gray slab, as shown in Figure 1. The slab had volumetric heat capacity cv = 4×10 5 J/(m 3 •K), refractive index n = 1.0, extinction coefficient β = 2000 m −1 , and scattering albedo ω = 0.8; furthermore, the thickness of the slab was L = 0.02 m, and the boundary emissivities for x = 0 and for x = L were ε0 = εL = 0.8. The initial temperature of the slab was T0 = 300 K; the temperature of the two slab walls were fixed at TH = 1000 K for x = 0, and TL = 300 K for x = L, respectively. The conductive thermal conductivity kc was assumed to be unknown and needed to be retrieved, and the actual value of kc was 0.02 W/(m•K). The 'measured temperature' was simulated by adding Gaussian noise to the solution of combined conduction and radiation, and various noise levels, S % T γ , were considered in the present study. The standard deviation, exp,k σ , of the 'measured temperature', at time tk, could be calculated from where ξ is a normal distribution number with a mean of zero and a standard deviation of unity. We also considered that the boundary temperature TH at x = L was not accurate, and its value was distributed about the mean value of TH = 1000 K with a normally distributed noise   First, it was obvious that the σ k c ,LB values increased with the increasing measurement noise γ TS % , as well as the increasing γ TH % values, which indicated that the accuracy of the retrieved parameters could be improved by performing an accurate experiment, and by using accurate model parameters. It is also interesting to note that the σ k c ,LB value for γ TS % = 3% and γ TH % = 1% was smaller than that for γ TS % = 1% and γ TH % = 9%; the decrease in measurement noise γ TS % (from 3% to 1%) did not result in a lower value of σ k c ,LB . This indicated that the accuracy of the retrieved parameter was affected comprehensively by the measurement noise and the model parameter uncertainty, and neither of the two factors could be neglected. Therefore, for some inverse identification problems, only trying to perform accurate experiments may not be sufficient in terms of improving the accuracy of the retrieved parameters. Figure 5 presents the estimated Fisher information matrix M (M is reduced to a scalar M as M has only one element), and the lower bound σ k c ,LB , for the standard deviation of the conductive thermal conductivity k c as a function of the dimensionless sensor position x s /L. The measurement noise was assumed to be γ TS % = 1%, and two different boundary temperature noise values, γ TH % = 0% and γ TH % = 5%, were investigated separately.
value of ,LB c k σ . This indicated that the accuracy of the retrieved parameter was affected comprehensively by the measurement noise and the model parameter uncertainty, and neither of the two factors could be neglected. Therefore, for some inverse identification problems, only trying to perform accurate experiments may not be sufficient in terms of improving the accuracy of the retrieved parameters.  Table 1. It can be seen from Figure 5 that the σ k c ,LB values first decreased, and then presented an increasing tendency with the increasing of the dimensionless sensor position x s /L. For γ TH % = 0%, the optimal sensor position was in the vicinity of x s /L = 0.5, and the minimum value of σ k c ,LB was about 5.5 × 10 −5 W/(m·K). Compared with the results for γ TH % = 0%, the minimum value of σ k c ,LB for γ TH % = 5% was increased to about 2.4 × 10 −4 W/(m·K); furthermore, the optimal sensor position moved from x s /L = 0.5 to a position in the vicinity of x s /L = 0.6, due to the fact that the boundary temperature error γ TH % affected the solution of the forward problem, especially for positions that were close to the boundary x = 0. Therefore, the sensor should be placed far away from the boundary to reduce its error effect.
The time-consuming Monte Carlo (MC) technique was employed to validate the designed sensor positions. We assumed that the three potential positions, x s /L = 0.5, 0.6, and 0.9, were available to place the temperature sensor for both γ TH % = 0% and γ TH % = 5%, respectively. For each sensor position and boundary temperature error γ TH % , 1000 independent inverse identifications were performed to retrieve k c ; thus, the standard deviations of the retrieved k c were calculated and compared with the σ k c ,LB value estimated via the CRB-based error analysis method. The results are presented in Table 1. It can be seen that a large discrepancy between the values estimated from the two methods was observed. This was due to the fact that the CRB-based method gave the lower bound of the uncertainty of the retrieved k c ; however, the aim of the present study was not to prove the correct quantitative error values. According to the MC simulation results, the best sensor position was x s /L = 0.5 and x s /L = 0.6 for γ TH % = 0% and γ TH % = 5%, respectively, while the worst position was x s /L = 0.9 for both γ TH % = 0% and γ TH % = 5%; this is consistent with the positions estimated using the CRB method. It indicates that the CRB method can be used to estimate the optimal experimental design for identification problems related to thermal properties.

Identification of Conductive and Radiative Properties: The Optimal Experimental Design
For problems regarding identification of conductive and radiative multiple properties, we considered the same physical model that was discussed in Section 3.1. The conductive thermal conductivity k c , extinction coefficient β, and scattering albedo ω of the slab were assumed to be unknown, and thus, needed to be retrieved, and their actual values were such that k c = 0.02 W/(m·K), β = 2000 m −1 , and ω = 0.8, respectively. The time duration of the 'experiment' was t S = 1000 s, and the sampling increment of time was ∆t = 2 s. The other parameters including the geometry parameter, the boundary condition parameters, and other properties were the same as those presented in Section 3.1.
For optimal experimental design problems involving the retrieving of only one parameter, the optimal sensor position could be easily identified according to the lower bound for the standard deviation values of the parameter to be retrieved. The optimal sensor position for multiple-parameter identification problems could not be determined directly from the lower bound for the standard deviation σ 2 u i ,LB of the parameter to be retrieved, as the minimum σ 2 u i ,LB for each parameter would not necessarily lead to the same sensor location. For this reason, it was necessary to define a new parameter to evaluate the retrieved parameters; in the present study, the parameter EU% was defined where N t is the number of sampling points, T S,pred (u i,fic , x e , t k ) is the predicted temperature at time t k and location x e using the fictitious parameter value u i,fic , and in the present study, we assumed that x e = L/2. The parameter EU% measured the integrated uncertainty of the recovered transient temperature response; the lower the EU%, the better the retrieved parameters. Thus, the best sensor position was the one that featured the lowest EU%. Figure 6 presents the estimated EU% with respect to various measurement noise γ TS % and boundary temperature error γ TH % values. The values considered for γ TS % and γ TH % ranged from 1 to 5%, with an increment of 1%. The temperature sensor was located at x s /L = 0.5. As with those used for one-parameter identification problems, the accuracy of the retrieved parameters could have been improved by performing more accurate experiments, and by using accurate model parameters when solving inverse conductive and radiative heat transfer problems. and % TH γ ranged from 1 to 5%, with an increment of 1%. The temperature sensor was located at xs/L = 0.5. As with those used for one-parameter identification problems, the accuracy of the retrieved parameters could have been improved by performing more accurate experiments, and by using accurate model parameters when solving inverse conductive and radiative heat transfer problems. It is interesting to note that the positions for each parameter were not consistent; therefore, the overall optimal sensor position could not be directly determined. Figure 7b shows the EU% values with respect to the dimensionless sensor location xs/L; the position corresponding to the minimum value of EU% was xs/L = 0.66, and the position should have been the optimal sensor position for the multiple-property identification problem. In fact, as the variation of xs/L between 0.4 and 0.8 led to only slight changes in the EU% values, the temperature sensor could have been placed at any position within this range. Therefore, the main consideration for researchers or engineers should be the ease and reliability of sensor installation.  Figure 7a presents the values of σ u,LB /u fic with respect to dimensionless sensor location x s /L; the measurement noise considered here was γ TS % = 2.0%. It is obvious that smaller σ u,LB /u fic values led to better retrieved results. The dimensionless positions corresponding to the minimum value of σ u,LB /u fic for k c , β and ω were x s /L = 0.76, 0.75, and 0.19, respectively, and the positions should have been the optimal sensor position for each parameter. It is interesting to note that the positions for each parameter were not consistent; therefore, the overall optimal sensor position could not be directly determined. It is interesting to note that the positions for each parameter were not co sistent; therefore, the overall optimal sensor position could not be directly determined Figure 7b shows the EU% values with respect to the dimensionless sensor locat xs/L; the position corresponding to the minimum value of EU% was xs/L = 0.66, and position should have been the optimal sensor position for the multiple-property ident cation problem. In fact, as the variation of xs/L between 0.4 and 0.8 led to only slig changes in the EU% values, the temperature sensor could have been placed at any posit within this range. Therefore, the main consideration for researchers or engineers shou be the ease and reliability of sensor installation.  Figure 7b shows the EU% values with respect to the dimensionless sensor location x s /L; the position corresponding to the minimum value of EU% was x s /L = 0.66, and the position should have been the optimal sensor position for the multiple-property identification problem. In fact, as the variation of x s /L between 0.4 and 0.8 led to only slight changes in the EU% values, the temperature sensor could have been placed at any position within this range. Therefore, the main consideration for researchers or engineers should be the ease and reliability of sensor installation. Figure 8a,b show the values of σ u,LB /u fic with respect to the dimensionless sensor position x s /L for the measurement noise value of γ TS % = 2.0%, and the boundary temperature error of γ TH % = 2.0% and γ TH % = 4.0%, respectively. Compared with the results reported in Figure 7a, the σ u,LB /u fic increased with the increasing of the γ TH % values. The optimal sensor positions for k c , β and ω moved slightly further away from the boundary x = 0 as γ TH % increased. Figure 8c reports the EU% as a function of the dimensionless sensor location for various γ TH % values of 0.0%, 2.0%, and 4.0%, respectively. It is obvious that the EU% values increased dramatically with the increasing of γ TH %, which means that the uncertainty of the boundary temperature had an obvious effect on the retrieved results. Furthermore, the optimal sensor position displayed a tendency to move slightly further away from the boundary x = 0. As the EU% values corresponding to x s /L between 0.4 and 0.8 changed slightly within this range, the temperature sensor could have been placed at any position within the range.  Figure 8c reports the EU% as a function of the dimensionless sens location for various % TH γ values of 0.0%, 2.0%, and 4.0%, respectively. It is obvious th the EU% values increased dramatically with the increasing of % TH γ , which means th the uncertainty of the boundary temperature had an obvious effect on the retrieved sults. Furthermore, the optimal sensor position displayed a tendency to move slightly f ther away from the boundary x = 0. As the EU% values corresponding to xs/L between and 0.8 changed slightly within this range, the temperature sensor could have been plac at any position within the range. As only one sensor position was designed, the total error variance k V and the perimental error variance k S , shown in Equation (13), were reduced to scalar Vk and Sk; th if Ok = Vk − Sk, the quantity Ok could then be used to measure the error variance caused As only one sensor position was designed, the total error variance V k and the experimental error variance S k , shown in Equation (13), were reduced to scalar V k and S k ; thus, if O k = V k − S k , the quantity O k could then be used to measure the error variance caused by an inaccurate boundary temperature T H . Therefore, the S k /V k and O k /V k values could be used to measure the relative contribution of the experimental noise and boundary temperature error to the total error at time t k , and thus, the integrated error contribution of experimental noise and boundary temperature T H over the whole time duration could be expressed as Figure 9a presents S k /V k and O k /V k as a function of time t. At this point, the measurement noise was still γ TS % = 2.0%, and the boundary temperature error considered in the present study was such that γ TH % = 2.0% and γ TH % = 4.0%, respectively. The value of γ TH % affected the contribution of the boundary temperature uncertainty to the overall error, and the O k /V k values for γ TH % = 2.0% were always smaller than those for γ TH % = 4.0%. The S k /V k and O k /V k values varied with time; first, the experimental error predominated, but, as time went on, the contribution of the boundary temperature error gradually increased, and started to become predominant, indicating that the main contribution factors changed in an alternating manner during the entirety of the time duration.

TH TH
The Sk/Vk and Ok/Vk values varied with time; first, the experimental error predominat but, as time went on, the contribution of the boundary temperature error gradually creased, and started to become predominant, indicating that the main contribution fact changed in an alternating manner during the entirety of the time duration. Figure 9b shows the integrated error contribution of the two error resources. First can be seen that the boundary error contribution was predominant for both % TH γ = 2. and % TH γ = 4.0%, and the error contribution was about 67% and 88%, respective Therefore, in order to improve the accuracy of the retrieved conductive and radiat properties, an effective method would involve trying to improve the accuracy of boundary temperature, TH, rather than concentrating on transient temperature measu ments. and TH.

Conclusions
In the present work, we proposed a stochastic Cramér-Rao bound (sCRB)-based n merical methodology to estimate the error of the conductive and radiative properties participating medium that was recovered from transient temperature measurements solving inverse heat transfer problems. The measurement noise and the inaccurate mo parameters were both taken into account in the analysis. The inverse identification pro lems of retrieving only one parameter and retrieving multiple parameters were illustrat separately. The proposed sCRB-based method was numerically validated by the tim consuming Monte Carlo simulations, and it was shown that the method was able to d termine, a priori, the error of the retrieved parameters. Based on the method, the optim  Figure 9b shows the integrated error contribution of the two error resources. First, it can be seen that the boundary error contribution was predominant for both γ TH % = 2.0% and γ TH % = 4.0%, and the error contribution was about 67% and 88%, respectively. Therefore, in order to improve the accuracy of the retrieved conductive and radiative properties, an effective method would involve trying to improve the accuracy of the boundary temperature, T H , rather than concentrating on transient temperature measurements.

Conclusions
In the present work, we proposed a stochastic Cramér-Rao bound (sCRB)-based numerical methodology to estimate the error of the conductive and radiative properties of participating medium that was recovered from transient temperature measurements by solving inverse heat transfer problems. The measurement noise and the inaccurate model parameters were both taken into account in the analysis. The inverse identification problems of retrieving only one parameter and retrieving multiple parameters were illustrated separately. The proposed sCRB-based method was numerically validated by the time-consuming Monte Carlo simulations, and it was shown that the method was able to determine, a priori, the error of the retrieved parameters. Based on the method, the optimal temperature sensor positions were designed to improve the accuracy of the retrieved parameters, and the relative error contributions of the error sources were also estimated.
The results show that: (1) the optimal sensor position is comprehensively determined by the factors of measurement noise as well as the uncertainties of inaccurate model parameters, and the optimal position varies with the levels of the error sources; (2) for problems regarding multiple parameter identification, the optimal position for each parameter may not be consistent, and thus, the optimal sensor position for the identification problem should be evaluated by the comprehensive parameter EU%, which is defined in Equation (21); and (3) the relative error contributions for each error source vary according to their error level, and the estimated relative error contributions can provide suggestions for improving the accuracy of the retrieved parameters.