research papers
Densitometry and temperature measurement of
by X-ray Compton scatteringaDepartment of Electronics and Informatics, Gunma University, 1-5-1 Tenjin-cho, Kiryu, Gunma 376-8515, Japan, bDepartment of Mechanical and Systems Engineering, Okayama University, Tsushima-Naka 3, Okayama 700-8530, Japan, and cJapan Synchrotron Radiation Research Institute (JASRI), 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5198, Japan
*Correspondence e-mail: sakuraih@gunma-u.ac.jp
Measurement of in situ technique to probe inside a combustion reaction.
by high-energy X-ray Compton scattering is reported. The intensity of Compton-scattered X-rays has shown a position dependence across the flame of the allowing us to estimate the temperature distribution of the combustion flame. The energy spectra of Compton-scattered X-rays have revealed a significant difference across the combustion reaction zone, which enables us to detect the combustion reaction. These results demonstrate that high-energy X-ray Compton scattering can be employed as anKeywords: Compton scattering; laminar flame; imaging; temperature; combustion reaction.
1. Introduction
Recent automobile technology requires measurement of complex combustion phenomena in order to achieve cleaner exhaust and less carbon dioxide emission. The combustion occurs due to complex interactions involving heat, mass and momentum transfer. The local flow temperature is a key consequence of the complex phenomena and must be accurately measured. Non-intrusive temperature measurements offer important advantages in combusting flows since they do not introduce disturbances that could alter the flow characteristics and change the temperature distribution (Xiao et al., 2000). Holographic temperature measurements have been suggested to accurately determine the in flames in order to infer the temperature distribution (Tieng et al., 1992; Xiao et al., 2000). However, these measurements are sensitive to mechanical noise, and hence it is not easy to obtain accurate measurements in a combustion instrument.
High-energy X-ray Compton scattering can provide a non-intrusive measurement of the temperature distribution of the N, is given by
in a flame with high accuracy and robustness, because it does not have the drawbacks of the holographic measurements. The intensity of Compton-scattered X-rays for monochromatic X-ray beams, dwhere Φ0 is the incident into a target object, t1 is the incident X-ray transmittance to a probing volume in the object, t2 is the scattered X-ray transmittance from the probing volume to an X-ray detector, ρe is the average electron density over the probing volume dV, and dσKN/dΩ is the Klein–Nishina differential (Sharaf, 2001). In the case of light-element gases, t1 and t2 are treated as constants to a good approximation. Therefore the intensity of Compton-scattered X-rays, dN, probes the electron density ρe. The ratio of molecular density to electron density, α, is defined by the following equation,
Here ρm denotes the molecular density. Since the flame system can be expressed by the law (Tieng et al., 1992), we can estimate the temperature distribution from the electron density as follows,
Here, ρm = n/V denotes the molecular density, T the temperature and R the constant. The ratio of the molecular density to the electron density, α, can be obtained assuming a combustion reaction. Therefore the intensity of Compton-scattered X-rays gives a temperature T if we regard the pressure P as a constant value of the ambient pressure, for example.
Furthermore, energy spectra of Compton-scattered X-rays give an electron momentum distribution, which is known as the Compton profile, J(pz), given as
Here, p [= (px, py, pz)] is the momentum of an electron in a target material. E1 denotes an incident X-ray energy, E2 a Compton-scattered X-ray energy with a scattering angle θ, m the electron mass and c the velocity of light. Since the molecules are not oriented in a certain direction, we can regard pz = px = py. Electron momentum density is denoted by n(p). ψi(r) denotes the wavefunction of the ith electronic state under the impulse approximation and hence χi(p) is the wavefunction of ith electronic state in momentum space. Therefore, X-ray Compton scattering provides information on a wavefunction, in other words, chemical bonds and molecular species, through the measurement of electron momentum distribution (Cooper, 1985). Recently, Itou et al. (2015) have applied this method to a discharging coin cell and successfully detected the change in the local electron density due to lithium ion insertion in the positive electrode.
Despite the merit of a non-intrusive measurement due to the large penetration of X-rays used in the X-ray Compton scattering technique, there are few reports of studies on combustion or flames. This is because Compton scattering from the gas phase is weak. However, the recent availability of high-energy and brilliant synchrotron radiation X-rays have made it possible to perform accurate measurements of Compton-scattered X-rays from low-density materials such as rare gases and N2 in ambient atmosphere (Sakurai et al., 2011; Kobayashi et al., 2011). In this study, we report the development of a densitometry and temperature measurement technique for use in flames by using Compton scattering.
2. Experiment
The experiment was performed at the BL08W beamline (Sakurai, 1998) at SPring-8. The X-ray source is an elliptical multipole wiggler operating in the linear polarization mode (Maréchal et al., 1998). The synchrotron radiation is monochromated and horizontally focused by an asymmetric Johann-type Si monochromator with (620) reflecting planes (Yamaoka et al., 2000) to deliver 115.6 keV X-ray beams. The entrance slit system defines the size of the incident X-rays to be 1 mm in height and 1 mm in width. Scattered X-rays from the in the probing volume were detected by a Ge solid state detector (Ge-SSD) with a collimating silt of 1 mm at a scattering angle of 124.9°. Therefore the probing volume was about 1 mm3. The intensity of the incident X-ray beams was monitored by an ion chamber for data normalization. An intensity map of the Compton scattered X-rays was measured for 10 s with 1 mm steps. The Compton scattered map was measured within 2.5 h.
A cylindrical din = 10 mm) was used to provide a laminar mixture of propane (6 wt%, 450 ml min−1) and air (94 wt%, 8 l min−1). The burner provided a laminar propane/air jet at the exit. The air and combustible gases consisted of pure synthetic gases. The temperature distribution was confirmed by a micro-thermocouple probe consisting of a Pt wire butt-welded with Pt-13% Rh wire. The wires were coated by SiO2 to avoid catalytic effects. The diameter of the micro-thermocouple probe was 0.5 mm. The burner was mounted on a movable stage, and the intensity of the Compton scattered X-rays from the probing volume was measured by scanning the along the vertical (y) direction and horizontal (x) direction on a plane perpendicular to the incident X-ray beam direction.
(inner diameter3. Results and discussion
Fig. 1(a) shows a side-view photograph of the self-sustaining flame. The bright flame front, at which CH* and OH* radicals emit light, shows a cone shape. Fig. 1(b) shows a cross-sectional map of Compton scattered X-ray intensities for the flame shown in Fig. 1(a). The relative intensities of Compton-scattered X-rays are shown by a graduated color scale. The orange region shows strong X-ray intensities and the blue region shows weak intensities. The two sides of the orange triangle region correspond to the bright flame front in Fig. 1(a) and thus the cross-sectional map reproduces the shape of the flame. Since the intensity of the Compton-scattered X-rays has a linear relationship with the electron density averaged over the probing volume of 1 mm3 [see equation (1)], the orange region inside the flame cone has a relatively higher electron density and the blue region outside the cone has a relatively lower electron density.
Fig. 2(a) shows the relative electron densities at the position of y = 5 mm and −10 mm < x < 20 mm shown in Fig. 1(b). Here the electron densities are normalized to that of the ambient air at x = 20 mm. The electron density shows a sudden decrease at the combustion reaction zone of |x| = 4.5 mm. The electron density remains almost constant in the inner flame region, |x| < 4 mm, which is shown as the orange region in Fig. 1(b). The electron density is one-fifth of the ambient air just in the outer flame region, |x| = 5 mm, which is shown as the blue region in Fig. 1(b).
We assume the following ideal chemical reaction (Tieng et al., 1992),
Therefore, the ratio of the molecular density to the electron density for air [7.8(O2 + 3.76N2); αair] is 0.0693; for the mixture of propane and air [C3H8 + 7.8(O2 + 3.76N2); α0] is 0.0679; and the ratio for gas after the combustion (3CO2 + 4H2O + 29.3N2 + 2.8O2; α1) is 0.0697. Since the system can be expressed by the law (Tieng et al., 1992), we can estimate the temperature distribution from equation (3). Here, α is regarded as a constant value of αair = 0.0693 within an error of 2% in the chemical reaction in (8). The ambient pressure P is regarded as a constant. The ambient temperature is 298 K at x = 20 mm.
The estimated temperature is shown in Fig. 2(b). The ambient temperature of 298 K is almost maintained in the inner flame region, |x| < 4 mm. The temperature increases suddenly at the flame front, reaching 1500 K just outside the flame (|x| = 5 mm). The value of 1500 K is confirmed by the thermocouple measurement as shown in Fig. 2(b) and is consistent with the previous report (Tieng et al., 1992). This suggests the validity of the temperature estimate given by the present Compton scattering analysis and suggests that the chemical reaction is dominated by equation (8). The blue region shown in Fig. 1(b), which is the low electron density region, corresponds to a temperature of 1500 K, although no luminescence is observed in the high-temperature region of 1500 K in Fig. 1(a).
The temperature given by the thermocouple measurement at x = 2.5 mm is 1200 K although the temperature measured by the present Compton scattering analysis is the ambient temperature of 298 K as shown in Fig. 2(b). It has been shown that the temperature suddenly decreases to ambient temperature within about a 1 mm distance from the flame front in the inner flame region (Tieng et al., 1992; Fristrom & Westenberg, 1957). Therefore the temperature at the position x = 2.5 mm, which is 2 mm away from the flame front of x = 4.5 mm, should be the ambient temperature of 298 K. This error in the thermocouple measurement may come from thermal conduction from the high-temperature region of 1500 K in the present thermocouple wire. Therefore the present Compton scattering analysis gives accurate position sensitivity in the measurement of the temperature distribution of the flame.
Fig. 3(a) shows energy spectra of Compton-scattered X-rays at the position of y = 5 mm and x = 0 mm, 4 mm, 4.5 mm and 5 mm shown in Figs. 1(b) and 2. Compton-scattered X-rays have a peak at 85.27 keV, which corresponds to pz = 0 (a.u.) in equations (4)–(6). The peak broadening comes from the Doppler shift of X-ray photons scattered by moving electrons in molecules. The intensities of Compton-scattered X-rays, which are defined as the sum of the counts between 80 keV and 90 keV, depend on the position x. The intensity decreases with increasing value of x, which has been discussed with reference to Figs. 1(b) and 2(a).
The energy spectra shown in Fig. 3(a) can be converted into electron momentum distributions by using equation (7), which are called Compton profiles, as shown by equations (4)–(6). Fig. 3(b) shows Compton profiles, J(pz), obtained from Fig. 3(a). The linear background contributions in Fig. 3(a) are removed. The Compton profiles shown in Fig. 3(b) give information on chemical bonding as explained earlier. In order to highlight the effects of chemical bonding, difference Compton profiles, ΔJxmm(pz), are obtained as follows,
Here Jxmm(pz) denotes a Compton profile at x mm as shown in Fig. 3(b). A normalization factor, Nxmm, is obtained so as to satisfy equation (10). The difference Compton profile ΔJxmm(pz) reflects the difference of chemical bonding between x = 0 mm, where the gas is a mixture of air and and x mm, where the gas is under combustion reaction.
Fig. 4 shows ΔJxmm(pz) just inside the combustion reaction zone (x = 4 mm), at the combustion reaction zone (x = 4.5 mm) and just outside the combustion reaction zone (x = 5 mm). ΔJxmm(pz) just inside the combustion reaction zone (x = 4 mm) is almost zero within experimental error (see Fig. 4). This shows that the chemical reaction does not take place. This is consistent with the fact that the ambient temperature of 298 K is almost maintained just inside the combustion reaction zone, |x| < 4 mm, as shown in Fig. 2(b). It should be noted that ΔJxmm(pz) in the combustion reaction zone (x = 4.5 mm) is almost zero within experimental error as shown in Fig. 4, although bright luminescence of the combustion reaction zone is observed as shown in Fig. 1(a) and the temperature is 800 K as shown in Fig. 2(b). This suggests that the main type of chemical bonding remains but the combustion reaction, such as a radical luminescence, starts at the combustion reaction zone. ΔJxmm(pz) just outside the combustion reaction zone (x = 5 mm) deviates from zero as shown in Fig. 4. This shows that the main type of chemical bonding, or molecular species, has changed drastically. This means that the combustion reactions complete just outside the combustion reaction zone and high-temperature region of 1500 K with no luminescence extending into the outer region of the combustion reaction zone as a result of the flame gas flow.
The features of ΔJxmm(pz) just outside the combustion reaction zone (x = 5 mm) in Fig. 4, which correspond to a change of chemical bonds, are a dip structure with |pz| < 0.5 a.u., positive values with |pz| < 1 a.u. and negative values with |pz| > 1 a.u. as shown in Fig. 4. The calculated ΔJ5mm(pz) assuming the ideal chemical reaction in equation (8) by using CRYSTAL09 (Dovesi et al., 2005, 2009) is also shown in Fig. 4. Although the dip structure with |pz| < 0.5 a.u. seems to be reproduced by the calculation, positive values with |pz| < 1 a.u. and negative values with |pz| > 1 a.u. are not reproduced. Since the ideal chemical reaction in equation (8) dominates the system as explained by Fig. 2(b), positive values with |pz| < 1 a.u. and negative values with |pz| > 1 a.u. can be attributed to partially generated radical species such as OH* and atomic H (Beduneau et al., 2009).
Fig. 5 shows calculated Compton profiles using CRYSTAL09 for air (O2 + 3.76N2), CO2, H2O, C3H8 and Hartree–Fock calculation (Biggs et al., 1975) for atomic H. The Compton profiles shown in Fig. 5 are normalized using equations (9) and (10). The CO2, which has smaller J(pz) values with |pz| < 0.5 a.u. than that of air, is ascribed to the dip structure with |pz| < 0.5 a.u. The C3H8 and the atomic H have larger J(pz) values with |pz| < 1 a.u. and lower values with |pz| > 1 a.u. than that of air. Because the combustion is complete and the temperature increases to 1500 K just out of the flame front as shown in Fig. 2, there seems to be no uncombusted C3H8 (Tieng et al., 1992; Fristrom & Westenberg, 1957). Thus the positive values with |pz| < 1 a.u. and negative values with |pz| > 1 a.u. are attributed to partially generated radical species such as atomic H accompanied by OH* radical. Therefore it is suggested that partially generated radical species such as OH* and atomic H contribute to the region just outside the combustion reaction zone (x = 5 mm) in Fig. 4 although the ideal chemical reaction in equation (8) dominates.
In conclusion, the intensity of Compton-scattered X-rays from the flame of a * and atomic H by the combustion reaction. Since it uses high-energy X-rays with high-penetrating power, it will be possible that the Compton scattering technique can be developed as an in situ and operando measurement for probing the inside of a combustion reaction.
has been measured. We have evaluated the local temperature from the and the temperature distribution across the flame is obtained with accurate position sensitivity. The energy spectra of Compton-scattered X-rays have shown a small but significant variation across the combustion reaction zone, which reflects partially generated radical species such as OHAcknowledgements
These experiments were performed with the approval of the Japan Synchrotron Research Institute (JASRI) (Proposal No. 2013B1158) and Human Resource Cultivation Center (HRCC), Gunma University (Proposal 2013). This work was supported by JSPS KAKENHI 24560364 [Grant-in-aid for Scientific Research (C)].
References
Beduneau, J. L., Kawahara, N., Nakayama, T., Tomita, E. & Ikeda, Y. (2009). Combust. Flame, 156, 642–656. CrossRef CAS Google Scholar
Biggs, F., Mendelsohn, L. B. & Mann, J. B. (1975). At. Data Nucl. Data Tables, 16, 201–309. CrossRef CAS Google Scholar
Cooper, M. J. (1985). Rep. Prog. Phys. 48, 415–481. CrossRef CAS Web of Science Google Scholar
Dovesi, R., Orlando, R., Civalleri, B., Roetti, C., Saunders, V. R. & Zicovich-Wilson, C. M. (2005). Z. Kristallogr. 220, 571–573. Web of Science CrossRef CAS Google Scholar
Dovesi, R., Saunders, V. R., Roetti, C., Orlando, R., Zicovich-Wilson, C. M., Pascale, F., Civalleri, B., Doll, K., Harrison, N. M., Bush, I. J., D'Arco, P. & Llunell, M. (2009). CRYSTAL09 User's Manual. University of Torino, Torino, Italy. Google Scholar
Fristrom, R. M. & Westenberg, R. M. (1957). Combust. Flame, 1, 217–228. CrossRef CAS Google Scholar
Itou, M., Orikasa, Y., Gogyo, Y., Suzuki, K., Sakurai, H., Uchimoto, Y. & Sakurai, Y. (2015). J. Synchrotron Rad. 22, 161–164. Web of Science CrossRef CAS IUCr Journals Google Scholar
Kobayashi, K., Itou, M., Hosoya, T., Tsuji, N., Sakurai, Y. & Sakurai, H. (2011). J. Phys. B, 44, 115102. CrossRef Google Scholar
Maréchal, X.-M., Hara, T., Tanabe, T., Tanaka, T. & Kitamura, H. (1998). J. Synchrotron Rad. 5, 431–433. Web of Science CrossRef IUCr Journals Google Scholar
Sakurai, Y. (1998). J. Synchrotron Rad. 5, 208–214. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sakurai, H., Ota, H., Tsuji, N., Itou, M. & Sakurai, Y. (2011). J. Phys. B, 44, 065001. CrossRef Google Scholar
Sharaf, J. M. (2001). Appl. Radiat. Isot. 54, 801–809. CrossRef PubMed CAS Google Scholar
Tieng, S. M., Lai, W. Z. & Fujiwara, T. (1992). Meas. Sci. Technol. 3, 1179–1187. CrossRef Google Scholar
Xiao, X., Choi, C. W. & Puri, I. K. (2000). Combust. Flame, 120, 318–332. CrossRef Google Scholar
Yamaoka, H., Hiraoka, N., Ito, M., Mizumaki, M., Sakurai, Y., Kakutani, Y., Koizumi, A., Sakai, N. & Higashi, Y. (2000). J. Synchrotron Rad. 7, 69–77. Web of Science CrossRef CAS IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.