Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Temperature errors in two-color pyrometry simultaneously considering reflection and combustion gas radiation

Open Access Open Access

Abstract

Temperature measurements are ubiquitous in combustion systems. However, the accuracy of surface temperature measurements of critical components operating in a harsh combustion gases environment is greatly affected by reflection and combustion gas radiation. In this paper, an analytical two-color pyrometry model was used to quantitatively analyze the temperature errors caused by the combination of reflection and H2O-CO2-CO-N2 mixture radiation. As the results indicate, the most significant contributors to the measurement errors are found to be the error arising from H2O-CO2-CO-N2 mixture absorption and emission for two-color pyrometer operating at long wavebands. The errors due to reflection predominate over the measurement errors measured at short wavebands. In a combustor where reflected radiation from high-temperature surrounding and hot/cool combustion gas is present, two-color pyrometry is practically inoperative as a consequence of its unacceptably large measurement error and high measurement sensitivity. When the intervening gas is isothermal and the optical distance from surface to detector is considered optically thin, the temperature error has linear growth with both the H2O-CO2-CO-N2 mixture concentration and viewing path length increasing. This linear change provides us a method of linear extrapolation to eliminate the effect of uncertain gaseous absorption and emission. The results of this work can be used as a theoretical support for the design and application of a two-color pyrometer in a gas-fired furnace.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Accurate surface temperature measurement in industrial processes have profound scientific importance and practical value from the viewpoints of efficient energy management, low manufacturing costs, high productivity and high product quality [1,2]. According to Kerr and Ivey that precise knowledge of the surface temperatures of turbine blades is quite vital as far as their overheating protection is concerned and is required for efficiency improvements, significant fuel consumption reductions, and financial savings [3]. The principal problem in the application of radiation pyrometry in combustion chamber is the large error caused by reflected ambient radiation and absorption and emission of radiation by interfering combustion gas.

Two-color pyrometry is an advantageous method of measuring surface temperature in industrial processes where unwanted radiation attenuation caused by intervening participating media presents [4,5]. In this method, the ratio of the radiant flux received from two separate spectral bands is used to develop a signal proportional to surface temperature. As long as the factors affecting the thermal radiation received affect both wavebands equally, the temperature inferred by this method remains unaffected [6]. Hence, two-color pyrometry has widely been used for the measurement of the temperature and concentration of soot or gas present in combustion flames [7,8]. Practical surface measurements by two-color pyrometry, however, are still prone to a number of uncertainties associated with environmental effects such as reflected ambient radiation and absorption and emission of radiation by interfering hot gases. Moreover, the environmental effect during practical measurement process is always different from that of calibration process. Therefore, it is quite important to understand these effects to achieve reliable and accurate temperature measurements.

The search for methods of modeling and eliminating the influences of emissivity and reflected radiation from extraneous sources has long challenged the field of radiation temperature measurement [912]. Saunders [13] developed a reflection model from which the errors and uncertainties for two-color pyrometers were analyzed and a mean of correcting for reflection error was presented. Araújo et al. [14] used Monte Carlo technique to estimate the reflection uncertainties for two-color pyrometry. Daniel et al. [15] proposed a reflection model to correct the reflection errors in surface temperature measurement of gas turbine blade using two-color pyrometry, taking the effects of target surface emissivity, view factor and ambient temperature into account. Nonetheless, a major problem with these aforementioned models is that they did not simultaneously consider the effect of intervening gas when dealing with reflection problems. Since pyrometer views the target through a column of gas, reflected radiation must transmit over a certain distance through the intervening gas before received by pyrometer. Any interaction between the intervening gas and the reflected radiation will change the intensity of the initial reflected radiation. Thus, a comprehensive analysis simultaneously taking both reflected ambient radiation and intervening gas radiation into account is required to precisely predict reflection error in two-color pyrometry.

Few researchers have attempted to compensate for the temperature error in two-color pyrometry due to absorption and emission of radiation by the furnace atmosphere or flue gas, in contrast with considerable investigations focusing on reflection problems. A possible reason is that most commercially radiation thermometers are designed to avoid significantly absorbing bands of flue gas. Even so, there is usually some residual overlap between the spectral responsivity of the thermometer and the absorption lines of the gas species [16]. Measurements of tube surface temperature in a large reformer furnace with a commonly used radiation thermometer operating on 1.0 µm which corresponds to the atmospheric window wavebands, undertaken by Saunders [16], indicated that the errors due to gaseous absorption and emission approach 6 °C per meter. In such situation, a more sophisticated understanding the mechanism of the measurement error caused by combustion gases is essential to the proper selection and application of radiometric temperature measurement techniques. Recently, Huang et al. [4,17] put forward a novel wavelength selection method for two-color pyrometry working in the presence of participating medium and investigated the effects of high-temperature H2O-CO2-N2 mixture radiation on two-color pyrometry. Despite this work is significant achievement, similarly, the shortfall is the fact that it solely focused on the influence of gas without simultaneously considering reflected ambient radiation effect.

Some researches in recent years have began to focus on the individual effects of reflection and combustion gaseous absorption and emission on single-color pyrometry [18,19]. However, few investigations have reported on the coupled effects of reflection and combustion gases on surface temperature measurement using two-color pyrometry. Motivated by these investigations, the purpose of this study is to quantify the coupled effect of reflected ambient radiation and gas absorption and emission on two-color pyrometry using an analytical two-color pyrometry model. The accuracy of this model was further verified by comparing the calculated heat flux distribution with the results of other numerical solutions using two non-isothermal and inhomogeneous combustion scenarios. Subsequently, the model was utilized to characterize the measurement error caused by the coupled influence of reflected ambient radiation and combustion gaseous absorption and emission. Finally, uncertainty and sensitivity analysis of the calculated temperature error to variations in the model’s input parameters were investigated.

2. Theoretical background

According to Saunders [13] that in many cases that the ambient temperature is not uniform, it is possible to characterise the background by a single effective temperature. Following Saunders [13] the planar geometry applied to the proposed temperature measurement system is composed of two boundless, parallel and diffusely emitting and reflecting boundaries which encompass an emitting and absorbing intervening gas layer, as depicted in Fig. 1. One boundary stands for the surroundings of practical measurement system, S2, the other boundary represents target surface S1. Because the effect of emissivity on surface temperature measurement is not the main point of this study, the emissivity of target surface is assumed to be a constant less than 1.0.

 figure: Fig. 1.

Fig. 1. Physical model of radiative temperature measurement in the presence of surroundings and combustion gases.

Download Full Size | PDF

Suppose the two-color pyrometer is placed on surroundings plane S2. The radiant flux sensed by the pyrometer includes components of emitted flux from S1 (Eemitted shown in Fig. 1) and flux originating from S2 reflected by S1 (Ereflected shown in Fig. 1), both influenced by intervening combustion gas. Obviously, the value of the radiant flux is related to target surface temperature T1, target surface emissivity ɛ1, ambient temperature T2, surroundings surface emissivity ɛ2, gas layer thickness L, gas temperature Tg and concentration Xg, which can be determined by solving the governing equations of radiative transfer in intervening combustion gas and its boundary conditions. Further details regarding the solution could be found in our previous studies [20,21].

Gas radiation properties are strong spectral dependency, which makes the spectral integration of transport equation computationally expensive even in one-dimension. The exact solution can only be given by line-by-line calculation, which is very complex and time consuming. As the practical aim of this paper is to measure the target surface temperature using two-color radiation thermometer which use a narrow spectral pass-band, a statistical narrow-band (SNB) model as described following is favored as a feasible alternative to line-by-line approach in these presented spectral radiative transfer calculations.

For an isothermal and homogeneous column of path-length l containing an individual absorbing gas at a mole fraction X and total pressure p, the SNB model with an exponential-tailed-inverse line-strength distribution provides the narrow band averaged transmissivity for Lorentz line profiles given as:

$${\bar{\tau }_{\Delta \lambda }}\textrm{ = }\exp \left[ { - \frac{{2\bar{\gamma }}}{{\bar{\delta }}}\left( {\sqrt {1 + \frac{{Xpl\bar{k}\bar{\delta }}}{{\bar{\gamma }}}} - 1} \right)} \right]$$
where $\bar{\delta }$ and $\bar{k}$, $\bar{\gamma }$ is the mean narrowband parameters. Recognize that the input parameter of radiative heat transfer equation is absorption coefficient, thus the obtained narrow band averaged transmissivity should be converted into a narrowband absorption coefficient. Different from the global gray-band approximation used in our previous study, a local gray-band approximation for implementing statistical narrow-band model, originally presented by Liu et al. [22], is chosen herein to express the radiative properties of interfering gas averaged on a narrowband. According to this approximation, an equivalent path independent narrow-band averaged local absorption coefficient can be obtained using the mean path-length of the gas layer under consideration (between xi and xi+1):
$${k_{\alpha \Delta \lambda }}({x_{i + 1/2}}) ={-} \frac{{\ln {{\bar{\tau }}_{\Delta \lambda }}({l_m})}}{{{l_m}}}$$
where subscript i stands for the spatial discretization, ${\bar{\tau }_{\Delta \lambda }}({l_m})$ represents narrow-band averaged transmissivity, and the local mean path-length lm for the planar slab is taken to be 1.8(xi+1-xi). Since ${\bar{\tau }_{\Delta \lambda }}({l_m})$ and lm are local properties of a combustion gas layer, this formulation in general yields spatially varying narrowband average absorption coefficient. It is expected that the grey-band approximation along with this local absorption coefficient is able to capture characteristics of non-grey radiation encountered in non-isothermal and inhomogeneous media since the local conditions are incorporated into this formulation.

After the detected radiant flux in waveband Δλ was obtained, for two-color pyrometry, target surface temperature can be derived by using bisection method to solve the following equation if the slope of instrumental emissivity set to 1.0:

$$\frac{{{E_{\Delta {\lambda _1}}}(T)}}{{{E_{\Delta {\lambda _2}}}(T)}} = \frac{{\int\limits_{\Delta {\lambda _1}} {\frac{{2\pi hc_0^2}}{{{n^2}{\lambda _1}^5[{e^{h{c_0}/n{\lambda _1}kT}} - 1]}}d\lambda } }}{{\int\limits_{\Delta {\lambda _2}} {\frac{{2\pi hc_0^2}}{{{n^2}{\lambda _2}^5[{e^{h{c_0}/n{\lambda _2}kT}} - 1]}}d\lambda } }}$$
where h=6.626×10−34 J•s is known as Plank’s constant, c0=2.998×108 m/s is known as speed of light in vacuum, k=1.3806×10−23 J/K is known as Boltzmann’s constant, n is refractive index.

3. Model validation

The main radiation from sufficient burning products of most natural gas and petroleum fuels consists of the vibration–rotation infrared spectrum of carbon dioxide, carbon monoxide and water vapor. Hence, a H2O-CO2-CO-N2 mixture was considered in the present study. In order to verify the accuracies of the detected heat flux by the two-color pyrometer in the proposed temperature measurement system, numerical calculations were performed using two complicated combustion scenarios which have been previously chosen as test cases by Liu et al. [22], Chu et al. [23], Howell et al. [24], Zhou et al. [25], Bordbar et al. [26] and G. Krishnamoorthy [27] using different gas radiative property models and molecular spectroscopic databases. The first scenario considered is a typical configuration of temperature and CO2 and H2O concentration distributions encountered in a counter-flow methane-air diffusion flame. Scenario 2 is representative of an oxy- fired diffusion flame with flue gas recirculation and consists of the same profiles in Scenario 1 with the spatial distribution of CO2 mole fractions increased by 0.6 relative to that of the air flame. The distribution of temperature and radiating gas concentrations is shown in Fig. 2, which was given in detail by Liu et al. [22] and Chu et al. [23]. In the two scenarios, the radiation calculations were performed in a one-dimensional parallel-plates enclosure containing a H2O-CO2-N2 mixture at 1 atm total pressure and the two plates surfaces are assumed to be black. The separation distance is 0.5 m and two bounding plates are at 300 K. The abscissa x of Fig. 2 represents the distance from one of the plates. x = 0 and x = 0.5 correspond to the locations of one plate and the other plate, respectively. The spectral range usd to calculate radiative heat flux in the present model is 50-11250 cm-1. These combustion scenarios provide testing platform to assess the effects of non-isothermal and inhomogeneous combustion gas distributions on the radiative heat flux and was chosen here so that the results obtained from this work can be compared to other published results.

 figure: Fig. 2.

Fig. 2. Temperature and CO2 and H2O mole fraction distributions for air combustion and oxy-fuel combustion scenarios [23].

Download Full Size | PDF

Figures 3 and 4 compare the net radiative heat flux distributions predicted by present model and other gas radiative property models using different spectroscopic databases for the air combustion scenario and oxy-fuel combustion scenario. Among of those distributions, the results of LBL calculation based on HITEMP 2010, indicated by red lines in Figs. 3 and 4, are used as the benchmark solution to evaluate the present model. As shown in Figs. 3 and 4, the proposed model, which is easier to formulate and therefore much quicker computationally, predicts accurate radiative heat flux distributions for the oxy-fuel combustion scenario compared to the air combustion scenario. For the air combustion scenario, the present model is in very close/excellent agreement with benchmark solutions in the middle of the enclosure and slightly over-predicts radiative heat flux by about 10% near the walls. The slight discrepancies are attributed to different angular discretization. Kim et al. [28] pointed out that the small differences would disappear, if the absorption coefficient shown in Eq. (2) is based on the actual lengths for each angle through the medium instead of on the mean beam length. The Comparison of results of the present model with benchmark solutions shows that the proposed model can predict fairly accurate radiative heat flux for one-dimensional non-isothermal and inhomogeneous media.

 figure: Fig. 3.

Fig. 3. Comparison of distributions of the net radiative flux calculated by present model and other gas radiative property models using different spectroscopic databases for the air combustion scenario with T1 = 300 K, ɛ1= 1.0, T2 = 300 K, ɛ2= 1.0, L=0.5 m.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Comparison of distributions of the net radiative flux calculated by present model and other gas radiative property models using different spectroscopic databases for the oxy-fuel combustion scenario with T1 = 300 K, ɛ1= 1.0, T2 = 300 K, ɛ2= 1.0, L=0.5 m.

Download Full Size | PDF

4. Results and discussion

4.1 Temperature errors due to reflection and combustion gases radiation

In order to figure out the temperature errors caused by reflection and combustion gases radiation, six measurement cases, depending on the relative target temperature, ambient temperature and gas temperature, were applied in this study. The two operating wavebands for two-color pyrometry model in this section are chosen as a two adjacent narrow-bands with a bandwidth of 25 cm-1. ΔT signifies the discrepancy between the calculated temperature using two-color pyrometry model and the assumed real target temperature.

Figures 59 indicate the ΔT for two-color pyrometer working on the wavenumbers in the spectrum of 50-11250 cm-1 with different ambient temperature and gas temperature. The part input parameters used to carry out these calculations are given in Table 1. The blue scatter points in Figs. 59 represent the total temperature errors caused by the combined influence of intervening gas H2O-CO2-CO-N2 and reflection, while the black one shows the temperature errors due to individual reflection with XH2O=XCO2=XCO=XN2=0, and the red one stands for the temperature errors produced by gases mixture alone with T2= T1. The values of ΔT above 1000 K over all studied wavebands are not covered in Figs. 59.

 figure: Fig. 5.

Fig. 5. The predicted temperature errors for two-color pyrometer with Tg= 800 K, T2 = 1600 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. The predicted temperature errors for two-color pyrometer with Tg= 800 K, T2 = 800 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. The predicted temperature errors for two-color pyrometer with Tg= 1600 K, T2 = 1600 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. The predicted temperature errors for two-color pyrometer with Tg= 1600 K, T2 = 800 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. The predicted total errors for two-color pyrometer for different ambient temperature with T1=1200 K and Tg=1200 K.

Download Full Size | PDF

Tables Icon

Table 1. Model input parameters

4.1.1 Case 1 (Tg <T1)

The ΔT produced by low-temperature intervening gas and high-temperature environment (Tg=400 K, T2=1600 K) and by low-temperature intervening gas and environment (T2=800 K, Tg=800 K) are shown in Figs. 5 and 6, respectively. Different from single-color thermometry [21], the values of ΔT for two-color pyrometer caused by intervening gas may larger or smaller than zero, depending on the operating wavelengths of two-color pyrometer. Roughly speaking, as Figs. 5 and 6 indicate, for two-color pyrometer operating at short wavebands, the temperature error caused by reflection predominates over the total errors, while the temperature error due to the gaseous absorption and emission takes an overwhelming component over the total error measured at long wavebands. The similar phenomenon occurs in single-color pyrometer [21]. The results suggest that a short-waveband thermometer operator should pay more attention on the effect of reflection compared to that of gas mixture, particularly when ambient temperature is high.

In a furnace where reflected radiation from high-temperature surrounding and low-temperature gas emission are always present, a two-color pyrometer is virtually useless because of the unacceptably large errors in the temperature measurement as shown in Fig. 5(a). As indicated in Fig. 6(a), the calculated temperature error due to reflection of low-temperature surrounding is insignificant, relative to that caused by low-temperature gaseous effects. As a result, the total error measured at a low ambient temperature and a low combustion gas temperature is dominated by the contribution due to the gas mixture absorption and emission effect.

Low-temperature gas absorption result in the underestimate of the target radiance, the opposite situation appears with the reflection of high-temperature environment, because it can cause an increase in the apparent target radiance. Hence, it is quite natural to imagine that the temperature error due to hot gas mixture can be properly compensated by adopting a low-temperature surrounding. Note that this method is suitable for single-color thermometry, but it is not necessarily suitable for two-color thermometry. Because two-color thermometers infer temperature from the ratio of the measured radiance in two different wavebands, the reduction or addition in signal of each waveband do not necessarily cause the increase or decrease in this ratio.

4.1.2 Case 2 (Tg >T1)

The ΔT caused by high-temperature intervening gas and environment (Tg=1600 K, T2=1600 K) and by high-temperature intervening gas and low-temperature environment (Tg=1600 K, T2=800 K) are depicted in Figs. 7 and 8, respectively. The emission from gas will significantly increase not only the emission of the target but also the reflection from surrounding, resulting in an immense error in the radiometric temperature measurement. For measurement systems consisted of high-temperature combustion gas and high-temperature surrounding, large uncertainties error of two-color pyrometer making them unsuitable for engineering applications.

When the high-temperature gas emission occurs in combination with the low-temperature surrounding, the total error is dominated by the contribution due to the gas mixture absorption and emission. Therefore, like the aforementioned treatment on the total error measured at a low ambient temperature and a low combustion gas temperature, the measured error in this case also can be minimized by employing atmospheric window wavebands with the condition that the high-temperature H2O-CO2-CO-N2 mixture present in low concentration.

4.1.3 Case 3 (Tg =T1)

Figure 9 shows the ΔT as a function of operating wavenumber with T1=1200 K, Tg=1200 K and ambient temperature T2 set as several values around target temperature. As the results indicate, the gaseous absorption and emission effects still greatly affects two-color radiation thermometer readings especially for long wavebands, even though Tg=T1. This is because the reflected ambient radiation absorbed by intervening gas is not equal to the gaseous emission terms when ambient temperature T2 is not equal to gas temperature. Accordingly, not accounting for gaseous absorption and emission effects will lead to the misjudgement of reflection errors.

4.2 Measurement uncertainty analysis

From a theoretical standpoint, the present model can be directly applied to correct for the coupled effects of reflection and combustion gas absorption and emission. From a practical perspective, however, the input parameters of the model, including ambient temperature T2, target surface emissivity ɛ1, intervening gas concentration Xg and temperature Tg, are often nearly impossible to measure accurately. Analyzing and summarizing the influence of measurement uncertainty of these parameters on the accuracy of radiometric temperature measurement is conducive to the recognition, understanding, and avoidance of measuring temperature errors.

4.2.1 Ambient temperature T2

The ΔT for two-color thermometry at various ambient temperatures T2 were calculated with T1=1200 K, XH2O=0.2, XCO2=0.2, XCO=0.2, XN2=0.4, L=100 cm, ɛ1=0.85, ɛ2=0.9, p=1 atm. Herein, the operating wavebands for two-color pyrometry are chosen as 0.9 um, 1.0 um, 1.25 um, and 3.9 µm, because these bands are the atmospheric window wavebands and are often used to design commercially radiation thermometers for special applications like measuring surface temperature in combustion chamber. The wavenumber combination of the wavebands mentioned above is shown in Table 2. Given the gas temperature has a vital effect on the measurement errors depending on whether the gas temperature is higher or lower than the target temperature, four different gas temperature are considered as shown in Fig. 10.

 figure: Fig. 10.

Fig. 10. Effect of ambient temperature T2 on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, XH2O=0.2, XCO2=0.2, XCO=0.2, XN2=0.4, L=100 cm, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different gas temperature: a) Tg=2000K, b) Tg=1600 K, c) Tg=1200 K, d) Tg=800 K.

Download Full Size | PDF

Tables Icon

Table 2. The wavenumber combination (η1, η2) for the atmospheric window wavebands

The results presented in Fig. 10 indicate that when ambient temperature is less than target temperature, the effect of ambient temperature on measuring temperature error is negligible except for the 3.9 um wavebands with the interfering gas at a high temperature. When ambient temperature is large than target temperature, the error increases nearly exponentially with ambient temperature. The rate of the increase is related with the operating waveband of two-color pyrometer. The shorter the waveband, the faster the rate.

4.2.2 Target emissivity ɛ1

Figure 11 represents the ΔT as a function of operating wavenumber with T1=1200 K, XH2O=0.2, XCO2=0.2, XCO=0.2, XN2=0.4, L=100 cm, ɛ2= 0.9, p=1 atm and the target surface emissivity at various values. Herein, four different combinations of ambient temperature and gas temperature are considered as shown in Fig. 11(a), 11(b), 11(c), and 11(d). It is noted that the measuring temperature error decreases with the target emissivity approaching 1.0. When the measurement system is at a low ambient temperature, the longer the waveband, the faster the rate of the increase of ΔT. On the contrary, for the high ambient temperature, the shorter the waveband, the faster the rate, which means the measuring temperature error is more sensitive to the change in target emissivity under this circumstance. As a result of the effect of high-temperature combustion gas, the measured error is a value not equal to zero when the target emissivity is equal to 1.0.

 figure: Fig. 11.

Fig. 11. Effect of target emissivity ɛ1 on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, XH2O=0.2, XCO2=0.2, XCO=0.2, L=100 cm, ɛ2= 0.9, p=1 atm and different combination of ambient temperature and gas temperature: a) Tg=1600 K, T2=1600 K, b) Tg=1600 K, T2=800 K, c) Tg=800 K, T2=1600 K, d) Tg=800 K, T2=800 K.

Download Full Size | PDF

For the measurement system composed of high-temperature interfering gas and low-temperature surrounding, the measured error due to low target emissivity is much greater than the measured error of other measurement systems. Therefore, we should pay special attention to the measurement accuracy of target emissivity in this measurement system. When the ambient temperature is larger than target temperature, the shorter the operating waveband, the larger the measured error. In contrast, the longer the operating waveband, the larger the measured error for low ambient temperature.

4.2.3 Gas temperature Tg

Figure 12 depicts the effect of gas temperature Tg on the ΔT for two-color thermometry. Herein, four different ambient temperatures are considered, the measurement uncertainty caused by the change in gas temperature less than 1200 K is insignificant except for 0.9 um wavebands with the surrounding at a high temperature. For gas temperature larger than target temperature, the value of ΔT increases nearly exponentially with gas temperature. The ΔT has the slowest growth rate at 1.0 µm waveband, followed by1.25 µm and 0.9 µm waveband, and the growth rate of ΔT at 3.9 µm band is the fastest. Consequently, two-color pyrometer working on 3.9 µm waveband will be more susceptible to the measurement uncertainty of hot gas temperature.

 figure: Fig. 12.

Fig. 12. Effect of gas temperature Tg on the ΔT for two-color pyrometer working at atmospheric window wavebands for T1=1200 K, XH2O=0.1, XCO2=0.1, XCO=0.1, XN2=0.7, L=100 cm, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different ambient temperature: a) T2=2000K, b) T2=1600 K, c) T2=1200 K, d) T2=800 K.

Download Full Size | PDF

4.2.4 Gas concentration Xg and measuring path length L

Figure 13 and Fig. 14 indicate the effect of H2O-CO2-CO-N2 concentrations characterized by XH2O:XCO2:XCO=1:1:1 and XN2=1-(XH2O+XCO2+XCO) and measuring path lengths on the ΔT for two-color pyrometer with T1=1200 K, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and four different combinations of ambient temperature and gas temperature, respectively. As shown in Fig. 13 and Fig. 14, the ΔT exhibit a linear relationship with both of gas concentration Xg and measuring path length L increasing. More important is that this linear relationship suggests a simple way to eliminate the effect of uncertain gaseous absorption and emission from the temperature measurement by two-color radiation thermometry. The effect of gaseous absorption and emission is eliminated as measuring path length or gas concentration approach zero. Thus a linear extrapolation of ΔT vs. L or X to L = 0 or X=0 yields a real target temperature immune from gaseous absorption and emission effect. Further detailed investigations on the linear extrapolation method will be summarized in our next study.

 figure: Fig. 13.

Fig. 13. Effect of measuring path length L on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, XH2O=0.1, XCO2=0.1, XCO=0.1, XN2=0.7, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different combination of ambient temperature and gas temperature: a) Tg=1600 K, T2=1600 K, b) Tg=1600 K, T2=800 K, c) Tg=800 K, T2=1600 K, d) Tg=800 K, T2=800 K.

Download Full Size | PDF

 figure: Fig. 14.

Fig. 14. Effect of H2O-CO2-CO-N2 concentration, characterized by XH2O:XCO2:XCO=1:1:1 and XN2=1-(XH2O+XCO2+XCO) on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, L=100 cm, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different combination of ambient temperature and gas temperature: a) Tg=1600 K, T2=1600 K, b) Tg=1600 K, T2=800 K, c) Tg=800 K, T2=1600 K, d) Tg=800 K, T2=800 K.

Download Full Size | PDF

For the temperature measurement system composed of low temperature gas, viewing path length and gas mixture concentration have a little effect on the measured temperature error. When the interfering gas is at a high temperature, the ΔT increase as the viewing path length or gas mixture concentration increases except for the error at the 3.9 um wavebands with the surroundings at a high temperature. Similarly, the ΔT has the slowest growth rate at 1.0 µm waveband, which means this waveband is more suitable for radiation temperature measurement in long measuring distance and high gas concentration environment.

4.3 Measurement sensitivity analysis

According to the measurement uncertainty analysis mentioned above, we could conclude that ambient temperature T2, target emissivity ɛ1, gas temperature Tg and concentration Xg, and measuring path length L are important parameters that affect the output of the proposed model. To quantitatively understand how the impact factors affect the behavior of the detection system, a sensitivity analysis method was employed herein. Detailed description of the analysis is given in Ref. [17].

The sensitivity coefficient of each parameter, which defined as the change in the quantity of the output of the proposed model divided by the change in the quantity of input parameter while keeping other parameters set at a given central value, is displayed in Table 3. As the Table shows, comparing with the temperature of low-temperature gas and surrounding, gas concentration and measuring path length, the temperature of hot gas and hot surrounding, and target emissivity have a relatively significant impact on the measured temperature using two-color pyrometry. As a result, the accurate measurement of the temperature of hot gas and hot surrounding, and target emissivity is necessary for the good estimation of real target temperature. Instead, slightly larger measurement uncertainties on the temperature of low-temperature gas and surrounding, gas concentration and measuring path length are permitted, since they have relatively small impact on the estimation of real target temperature.

Tables Icon

Table 3. Sensitivity coefficient of each parameter on the calculated temperature for two-color pyrometer operating at atmospheric window wavebands

5. Conclusion

This paper performed a theoretical analysis which has as aim to quantify the error committed in the temperature measurement by technique two-color pyrometry. In particular, the temperature error was evaluated as a function of operating wavelength of two-color pyrometry, ambient temperature, target emissivity, H2O-CO2-CO-N2 mixture temperature and concentration, and measuring path length. The target surface and background are considered to be a gray body with constant emissivity. Sensitivity analysis of the derived temperature to variations in the model’s input parameters were investigated. The main predictions for the coupling influences of intervening gas mixture and reflected ambient radiation on two-color pyrometry reveal the following.

  • (1) For the temperature measurement system consisted of high-temperature surrounding and high-temperature or low-temperature combustion gas, the unacceptably large measurement errors and high measurement sensitivity render two-color pyrometry virtually useless.
  • (2) Generally speaking, for two-color ratio thermometer operating at short wavebands, the total errors due to the coupling influences are dominated by the error arising from reflection. The errors due to the gas mixture absorption and emission predominate over the total error measured at long wavebands.
  • (3) When the intervening gas is isothermal and the optical distance from the surface to the detector is considered optically thin, the calculated temperature error for two-color pyrometer having a spectral responsivity near 3.9, 1.25, 1.0 and 0.9 µm is the linear growth in magnitude with both gas mixture concentration and viewing path length increasing. This linear change provides us a method of linear extrapolation to eliminate the effect of uncertain gaseous absorption and emission.
  • (4) Either in terms of measurement error or measurement sensitivity, the effect of reflection from low-temperature surrounding is insignificant. The absorption and emission from low-temperature combustion gas has a fairly important effect on the temperature measurement for two-color radiation thermometry operating at long wavebands.

Funding

National Natural Science Foundation of China (51775112, 52005103, 71801046); Guangdong Provincial Applied Science and Technology Research and Development Program (2019B1515120095); Intelligent Manufacturing PHM Innovation Team Program (2018KCXTD029, TDYB2019010); Chongqing Science and Technology Commission (2019jcyj-zdxmX0013).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. D. P. N. Dewitt and G. D. Nutter, Theory and Practice of Radiation Thermometry (1988).

2. Y. Zhao, J. Lv, K. Zheng, J. Tao, and J. Liang, “Multispectral radiometric temperature measurement algorithm for turbine blades based on moving narrow-band spectral windows,” Opt. Express 29(3), 4405–4421 (2021). [CrossRef]  

3. C. Kerr and P. Ivey, “An overview of the measurement errors associated with gas turbine aeroengine pyrometer systems,” Meas. Sci. Technol. 13(6), 873–881 (2002). [CrossRef]  

4. Y. Huang, M. Long, D. Chen, H. Duan, K. Tan, L. Gui, and T. Liu, “A new wavelength selection criterion for two-color pyrometer interfered with participating media,” Infrared Phys. Technol. 93, 136–143 (2018). [CrossRef]  

5. X. Zhou, M. J. Hobbs, B. S. White, J. P. R. David, J. R. Willmott, and C. H. Tan, “An InGaAlAs-InGaAs two-color photodetector for ratio thermometry,” IEEE Trans. Electron Devices 61(3), 838–843 (2014). [CrossRef]  

6. B. Müller and U. Renz, “Development of a fast fiber-optic two-color pyrometer for the temperature measurement of surfaces with varying emissivities. Rev Sci Instrum,” Rev. Sci. Instrum. 72(8), 3366–3374 (2001). [CrossRef]  

7. S. Deep, Y. Krishna, and G. Jagadeesh, “Temperature characterization of a radiating gas layer using digital-single-lens-reflex-camera-based two-color ratio pyrometry,” Appl. Opt. 56(30), 8492–8500 (2017). [CrossRef]  

8. T. Fu, X. Cheng, and Z. Yang, “Theoretical evaluation of measurement uncertainties of two-color pyrometry applied to optical diagnostics,” Appl. Opt. 47(32), 6112–6123 (2008). [CrossRef]  

9. Y. Lü, X. He, Z. H. Wei, Z. Y. Sun, and S. T. Chang, “Ambient temperature-independent dual-band mid-infrared radiation thermometry,” Appl. Opt. 55(9), 2169–2174 (2016). [CrossRef]  

10. R. R. Corwin and A. Rodenburghii, “Temperature error in radiation thermometry caused by emissivity and reflectance measurement error,” Appl. Opt. 33(10), 1950–1957 (1994). [CrossRef]  

11. J. Xing, B. Peng, Z. Ma, X. Guo, L. Dai, W. Gu, and W. Song, “Directly data processing algorithm for multi-wavelength pyrometer (MWP),” Opt. Express 25(24), 30560 (2017). [CrossRef]  

12. L. Jiafeng, S. Dai, W. Chen, N. Peng, W. Song, and J. Xing, “Generalized inverse matrix-exterior penalty function (GIM-EPF) algorithm for data processing of multi-wavelength pyrometer (MWP),” Opt. Express 26(20), 25706–25720 (2018). [CrossRef]  

13. P. Saunders, “Reflection errors and uncertainties for dual and multiwavelength pyrometers,” High Temp.-High Press. 32(2), 239–249 (2000). [CrossRef]  

14. A. Araújo and N. Martins, “Monte Carlo simulations of ambient temperature uncertainty determined by dual-band pyrometry,” Meas. Sci. Technol. 26(8), 085016 (2015). [CrossRef]  

15. D. Ketui, C. Feng, and S. Gao, “Single wavelength and ratio pyrometry reflection errors in temperature measurement of gas turbine blade,” Measurement 86, 133–140 (2016). [CrossRef]  

16. P. Saunders, Radiation Thermometry: Fundamentals and Applications in the Petrochemical Industry (SPIE press, 2007), Vol. 78.

17. Y. Huang, M. Long, H. Fan, L. Gui, D. Chen, and H. Duan, “Quantifying the Effects of Combustion Gases’ Radiation on Surface Temperature Measurements Using Two-Color Pyrometry,” Energy Fuels 33(4), 3610–3619 (2019). [CrossRef]  

18. S. Gao, L. Wang, C. Feng, and K. D. Kipngetich, “Analyzing the influence of combustion gas on a gas turbine by radiation thermometry,” Infrared Phys. Technol. 73, 184–193 (2015). [CrossRef]  

19. H. Li, G. Wang, N. Nirmalan, S. Dasgupta, and E. R. Furlong, “Passive absorption/emission spectroscopy for gas temperature measurements in gas turbine engines,” in ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition, (American Society of Mechanical Engineers, 2011), 19–28.

20. Y. Huang, M. Long, D. Chen, K. Tan, H. Duan, and P. Xu, “Effect of hot water vapor on strand surface temperature measurement in steel continuous casting,” Int. J. Therm. Sci. 138, 467–479 (2019). [CrossRef]  

21. Y. Huang, M. Long, J. Zhao, H. Fan, D. Chen, K. Tan, and H. Duan, “Coupled effects of reflection and absorptive gas mixture on surface temperature determined by single color pyrometer,” J. Quant. Spectrosc. Radiat. Transfer 228, 111–123 (2019). [CrossRef]  

22. F. Liu, Ö. Gülder, G. Smallwood, and Y. Ju, “Non-grey gas radiative transfer analyses using the statistical narrow-band model,” Int. J. Heat Mass Transfer 41(14), 2227–2236 (1998). [CrossRef]  

23. H. Chu, F. Liu, and H. Zhou, “Calculations of gas thermal radiation transfer in one-dimensional planar enclosure using LBL and SNB models,” Int. J. Heat Mass Transfer 54(21-22), 4736–4745 (2011). [CrossRef]  

24. J. Ma, B.-W. Li, and J. R. Howell, “Thermal radiation heat transfer in one-and two-dimensional enclosures using the spectral collocation method with full spectrum k-distribution model,” Int. J. Heat Mass Transfer 71, 35–43 (2014). [CrossRef]  

25. S. Zheng, Y. Yang, and H. Zhou, “The effect of different HITRAN databases on the accuracy of the SNB and SNBCK calculations,” Int. J. Heat Mass Transfer 129, 1232–1241 (2019). [CrossRef]  

26. M. H. Bordbar, G. Węcel, and T. Hyppänen, “A line by line based weighted sum of gray gases model for inhomogeneous CO2–H2O mixture in oxy-fired combustion,” Combust. Flame 161(9), 2435–2445 (2014). [CrossRef]  

27. G. Krishnamoorthy, “A new weighted-sum-of-gray-gases model for oxy-combustion scenarios,” Int. J. Energy Res. 37(14), 1752–1763 (2013). [CrossRef]  

28. T. K. Kim, J. A. Menart, and H. S. Lee, “Nongray radiative gas analyses using the SN discrete ordinates method,” J. Heat Transfer. 113(4), 946–952 (1991). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (14)

Fig. 1.
Fig. 1. Physical model of radiative temperature measurement in the presence of surroundings and combustion gases.
Fig. 2.
Fig. 2. Temperature and CO2 and H2O mole fraction distributions for air combustion and oxy-fuel combustion scenarios [23].
Fig. 3.
Fig. 3. Comparison of distributions of the net radiative flux calculated by present model and other gas radiative property models using different spectroscopic databases for the air combustion scenario with T1 = 300 K, ɛ1= 1.0, T2 = 300 K, ɛ2= 1.0, L=0.5 m.
Fig. 4.
Fig. 4. Comparison of distributions of the net radiative flux calculated by present model and other gas radiative property models using different spectroscopic databases for the oxy-fuel combustion scenario with T1 = 300 K, ɛ1= 1.0, T2 = 300 K, ɛ2= 1.0, L=0.5 m.
Fig. 5.
Fig. 5. The predicted temperature errors for two-color pyrometer with Tg= 800 K, T2 = 1600 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.
Fig. 6.
Fig. 6. The predicted temperature errors for two-color pyrometer with Tg= 800 K, T2 = 800 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.
Fig. 7.
Fig. 7. The predicted temperature errors for two-color pyrometer with Tg= 1600 K, T2 = 1600 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.
Fig. 8.
Fig. 8. The predicted temperature errors for two-color pyrometer with Tg= 1600 K, T2 = 800 K. (a) Total error (b) Reflection error (c) Gaseous absorption and emission error.
Fig. 9.
Fig. 9. The predicted total errors for two-color pyrometer for different ambient temperature with T1=1200 K and Tg=1200 K.
Fig. 10.
Fig. 10. Effect of ambient temperature T2 on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, XH2O=0.2, XCO2=0.2, XCO=0.2, XN2=0.4, L=100 cm, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different gas temperature: a) Tg=2000K, b) Tg=1600 K, c) Tg=1200 K, d) Tg=800 K.
Fig. 11.
Fig. 11. Effect of target emissivity ɛ1 on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, XH2O=0.2, XCO2=0.2, XCO=0.2, L=100 cm, ɛ2= 0.9, p=1 atm and different combination of ambient temperature and gas temperature: a) Tg=1600 K, T2=1600 K, b) Tg=1600 K, T2=800 K, c) Tg=800 K, T2=1600 K, d) Tg=800 K, T2=800 K.
Fig. 12.
Fig. 12. Effect of gas temperature Tg on the ΔT for two-color pyrometer working at atmospheric window wavebands for T1=1200 K, XH2O=0.1, XCO2=0.1, XCO=0.1, XN2=0.7, L=100 cm, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different ambient temperature: a) T2=2000K, b) T2=1600 K, c) T2=1200 K, d) T2=800 K.
Fig. 13.
Fig. 13. Effect of measuring path length L on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, XH2O=0.1, XCO2=0.1, XCO=0.1, XN2=0.7, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different combination of ambient temperature and gas temperature: a) Tg=1600 K, T2=1600 K, b) Tg=1600 K, T2=800 K, c) Tg=800 K, T2=1600 K, d) Tg=800 K, T2=800 K.
Fig. 14.
Fig. 14. Effect of H2O-CO2-CO-N2 concentration, characterized by XH2O:XCO2:XCO=1:1:1 and XN2=1-(XH2O+XCO2+XCO) on the ΔT for two-color pyrometer working at atmospheric window wavebands with T1=1200 K, L=100 cm, ɛ1= 0.85, ɛ2= 0.9, p=1 atm and different combination of ambient temperature and gas temperature: a) Tg=1600 K, T2=1600 K, b) Tg=1600 K, T2=800 K, c) Tg=800 K, T2=1600 K, d) Tg=800 K, T2=800 K.

Tables (3)

Tables Icon

Table 1. Model input parameters

Tables Icon

Table 2. The wavenumber combination (η1, η2) for the atmospheric window wavebands

Tables Icon

Table 3. Sensitivity coefficient of each parameter on the calculated temperature for two-color pyrometer operating at atmospheric window wavebands

Equations (3)

Equations on this page are rendered with MathJax. Learn more.

τ ¯ Δ λ  =  exp [ 2 γ ¯ δ ¯ ( 1 + X p l k ¯ δ ¯ γ ¯ 1 ) ]
k α Δ λ ( x i + 1 / 2 ) = ln τ ¯ Δ λ ( l m ) l m
E Δ λ 1 ( T ) E Δ λ 2 ( T ) = Δ λ 1 2 π h c 0 2 n 2 λ 1 5 [ e h c 0 / n λ 1 k T 1 ] d λ Δ λ 2 2 π h c 0 2 n 2 λ 2 5 [ e h c 0 / n λ 2 k T 1 ] d λ
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.