† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant No. 11272338) and the Fund from the Science and Technology on Scramjet Key Laboratory, China (Grant No. STSKFKT2013004).
The aim of the present work is to quantitatively measure the hydroxyl radical concentration by using LIF (laser-induced fluorescence) in flame. The detailed physical models of spectral absorption lineshape broadening, collisional transition and quenching at elevated pressure are built. The fine energy level structure of the OH molecule is illustrated to understand the process with laser-induced fluorescence emission and others in the case without radiation, which include collisional quenching, rotational energy transfer (RET), and vibrational energy transfer (VET). Based on these, some numerical results are achieved by simulations in order to evaluate the fluorescence yield at elevated pressure. These results are useful for understanding the real physical processes in OH-LIF technique and finding a way to calibrate the signal for quantitative measurement of OH concentration in a practical combustor.
Laser-diagnostic techniques have played an increasing role in combustion research over the last 20–30 years. These techniques are important for detecting temperatures, flow fields and some special species or radical concentration in combustion flow field, because these techniques are nonintrusive and robust. However, despite the rapid development of various techniques based on laser and spectroscopy, the quantitative measurement in the practical combustor is still difficult.
Laser-induced fluorescence (LIF) is probably the most common and important method to measure the densities of some minor species and radicals such as OH, CH, NO, and CO.[1–5] In particular, planar laser-induced fluorescence (PLIF) is more popular due to the high sensitivity, high spatiotemporal resolution. Laser-induced fluorescence relies on the ability to excite molecules or radicals from a ground state to a higher energy state. However, in LIF measurement, the induced fluorescence signal density captured by the intensified-CCD or PMT not only depends on the absolute concentration of the species excited by the laser pulse, but also is related to other processes such as electronic quenching, vibrational energy transfer (VET) and rotational energy transfer (RET).[6–8] These processes are mainly due to the molecular collisional partner, pressure, temperature and quantum state. Then, LIF measurement for quantitative species concentration becomes complex and uncertain in some quasi-practical combustors, such as gas turbine, engine, and other industry facilities. Therefore, it is important and necessary to find the relationship between the fluorescence signal and the real absolute concentrations of species in practical combustion fields.
It should be pointed out that over the years attempts have been made for quantitative measurement of minor species such as OH, NO, and CH. Some collision-insensitive techniques such as saturated LIF[9,10] and the picosecond LIF[11] have been developed to eliminate the influences of quenching, and other techniques combined with LIF and absorption spectroscopy have been used to measure the absolute density directly.[12,13] According to this way a bidirectional laser beam configuration has been demonstrated to be collision independent and self-calibrated for absolute OH number density measurements.[13] But most of the studies using these techniques mentioned above have been carried out in a laboratory-scale burner, which are excellent calibration methods. Therefore to quantitatively measure the OH concentration by LIF, the modeling and analysis on the complex LIF process in high pressure and turbulent flames are still needed. The detailed simulations of OH excitation dynamics have been conducted and continuously developed in recent years, which is conducive to understanding the phenomena in the LIF mechanism.[12,14–17]
In the present work, we analyze the physical processes of OH-LIF diagnostics in flame and make clear the influence of parameters according to pressure. Our goal is to present the real process and evaluate the fluorescence yield in high pressure combustion field, which is important and significant for performing laser diagnostics and quantitative measurement in a practical combustor and facility. In our approach, the phenomena involved in the experimentally measured signal of OH concentration are taken into account as much as possible.
Describing the various processes in LIF experiments accurately is important and critical to achieving the quantitative measurement of radical concentration in high pressure combustion. Another aspect is to obtain a fluorescence signal with high quality and signal-to-noise ratio by some optical strategies and optimized designs, which is not discussed in the current study. We prefer to discuss elementary processes: absorption line shape broadening, collisional transitions and quenching, which are all enhanced when the pressure is augmented.
Based on the typical LIF experiments and setups,[1] the fluorescent photons per unit of time captured by the photodetectors such as PMT and ICCD, Sf, can be expressed as follows:
Laser-induced fluorescence originates from an interaction process between laser and matter, simply described as a spontaneous light emission that takes place between two quantum states of some species via absorption of a photon with suitable energy. As illustrated in the energy level structure of OH shown in Fig.
The fluorescence signal can then be simulated by solving the population rate equations of the states involved.
Following the OH-LIF principle,[18] a laser pulse at wavelength λ according to a given transition between a lower rotational level J″ in vibrational level ν″ of electronic state X2Π and (J′, ν′) of A2Σ+ induces the fluorescence, which is emitted at different wavelengths. In the linear regime of LIF, if a broadband detection scheme including both A2Σ+–X2Π (0, 0) and (1, 1) bands are considered, the relationship is expressed as follows:
The electronic transition (A2Σ+–X2Π) is the most commonly used in OH-LIF diagnostics. Figure
In several publications[7,8,12,15,18,19], LIF simulation models have been developed, including these collisional processes describing the detailed dynamics of the levels that are involved. Following Ref. [18], the time-dependent population of a collision-populated rotational level j is described by
The quenching, VET and RET-rate are all induced by the collision partner species, and we define collision coefficients as
The rates are determined by the species and the number densities, therefore the quenching and other collisional effects become stronger at the same mole fraction distribution at higher pressure. In many studies, a simplified two-level system has been obtained. The effective spontaneous emission coefficient Aeff and the effective quenching rate Qeff are adopted and the quenching models and the related collision coefficients are clear.
For electronic transitions in combustion environments at higher pressure than atmospheric, one can assume that the Einstein emission coefficient is much smaller than the collisional quench rate. The fluorescence yield is then approximately equal to
However, in the present work, the simulation considering the detailed collisional effects is very complicated notably because it requires accurate knowledge of VET and RET, and this calculation cannot be dealt with without the help of a systematic model and software such as the LASKIN code.[20] Then we use this software to simulate the change of number densities in different rotational and vibrational levels.
The structure and the number density distributions of the molecular energy levels are the key to quantitative measurement. It should be specially noted that the ground state of the OH radical has two spin components which give rise to two rotational energy ladders due to spin–orbit interaction. There is an idealized case where a specific coupling between angular momenta is said to be dominant, which is called Hund’s case.[21] We denote the case with subscript 1 or 2 that has the following meaning for the total angular momentum quantum number J:
It is now important to quantitatively discuss and estimate the aspects that will influence the fluorescence signal and its interpretation based on the laser-induced fluorescence modeling. The fluorescence signal per unit of absorbing species mole fraction depends on the absorption strength (depends on the temperature through the Boltzmann population fraction), the overlap integral (depends on the line broadening), and the fluorescence yield (quenching, VET and RET). These three parameters are examined in this section.
The part ∫ Φlaser(ν)Φabs(ν)dν in function (
Based on the complex fine structure of energy levels, the broadening mechanism of absorption line is supported by two main phenomena named collisional broadening and Doppler broadening, which can be expressed as Lorentzian-type and Gaussian-type functions respectively.[23,24] Then the broadening is given by the pressure and translational temperature of the system, and another line shape Voigt is more suitable and widely accepted at high pressure.
The Voigt profile results from the convolutions of the two types of line shapes mentioned before. Theoretically, Doppler broadening depends only on temperature, and the broadening induced by collision is influenced by not only temperature but pressure and composition species. Therefore the Doppler broadening can be calculated in a straightforward manner, but evaluating the collision-induced effects is more complicated and difficult. In terms of the Voigt function, the shape function, (υ −υ0) can be expressed as[24]
The collisional broadening studies and coefficient measurement of the OH A–X transition at high temperature and pressure environments were performed by Rea et al.,[23] Battles and Hanson[25] for (0,0) band, Kessler et al.,[26] and Atakan et al.[27] for the (1,0) band. In the present study, we do not attempt to accurately measure the coefficients for the collision broadening process. Then for hard-sphere collision theory we adopt the n = −0.66, and parameter a = 0.041.
Assume that the laser is tuned to the center of the absorption line shape. Figure
The results show that as the pressure increases, the spectral overlap decreases because of the collisional broadening of OH absorption line. Secondly, the Doppler broadening increases as the temperature increases and the collisional broadening decreases. As a consequence, the Doppler broadening dominates at low pressure. However, at higher pressure, the collisional broadening dominates, so that the OH absorption line becomes thinner and the spectral overlap increases as the temperature increases. In comparison with the results in Ref. [12], there is some deviation between the exact values of spectral overlap because of the different laser characteristics, but the tendencies in the present work are in accordance with others.
According to statistical mechanics, the population distribution of an ensemble of molecules in thermal equilibrium is given by the Boltzmann equation. In laser-induced fluorescence procession, the fluorescence transition is composed of many transitions of OH radicals from the excited vibrational state, and the Boltzmann fraction fB depends on temperature.
The rotational-level populations in the state (ν″ = 0) of X2Π at different temperatures are illustrated in Fig.
In the present work, for the initial attempt to model the influence of quenching, VET and RET on the LIF signal, we evaluate the fluorescence quantum efficiency under some simple and typical conditions. For the purpose, the burned gases are considered to be a mixture of 72% N2, 19% H2O, and 9% CO2 and the temperature is varied between 800 K and 2000 K at elevated pressure from 0.1 MPa to 1.0 MPa as indicated in the literature.[12]
To evaluate the fluorescence yields at different temperatures and elevated pressure, a lot of groups have made some efforts to build and develop various models based on experiments and simulations. Paul,[7] and Garland and Crosley[29] developed models to simulate the corresponding quenching cross sections as a function of the temperature. Based on and utilizing these published results, Höinghaus et al. built a detailed rate equation model (the LASKIN program packet)[20] for simulating the energy transfer in OH laser-induced fluorescence, which includes all possible state-to-state energy transfer processes such as electronic quenching, VET, RET, depolarization, etc. Then for our case of linear excitation without saturation, the LASKIN code is employed to illustrate the influence of combustion condition on fluorescence efficiency.
In the CH4/air flame, electronic quenching is dominated by H2O and CO2, and sometimes CO is important too. In LASKIN, the quenching cross sections for various collision partners are calculated by using a two-parameter formula which is T-dependent.[18]
The quenching rate for H2O is chosen based on the default value of LASKIN, and the cross section for CO2 is cited from the model of Paul, and the collision of CO is ignored in this case. As mentioned before, Hund’s case (b) notation is used here and the population fractions of spin-split states (F1 and F2) are calculated separately. The important step in the simulation is to calculate the time-dependent populations in the upper-state rotational levels following laser excitation.
Figures
It is clear that the population of the excited level F1(6) decreases at elevated pressure for both excitation lines mainly due to quenching. For the excitation line P1(7) (0–1), the greater population is achieved at higher temperature, which can be explained by the Boltzman fraction and population distribution of the ground level. When increasing temperature, the difference between the levels close to F1(6), such as F1(5)–F1(7) and F2(5)–F2(6)–F2(7), becomes small. At the same time, there is a big gap between these two groups. Thus this tendency is a direct result of the RET in collisions with H2O (i.e., F1–F1 and F2–F2 transitions are more likely than F1–F2 and F2–F1 transitions).[18]
Figure
The LASKIN program packet provides a powerful tool for the modelling of energy transfer in laser-induced fluorescence. But it should be noted that our simulations using LASKIN are not absolutely certain because there are still a lot of uncertainties in collision cross sections for RET, VET and electronic quenching of OH radicals. We hope to achieve better agreement between experimental measurements and theoretical calculations in the next step.
The aim of the present work is to achieve quantitative measurement of hydroxyl radical concentration by using laser-induced fluorescence in flame. Therefore, the physical process of OH-LIF diagnostics in flame is analyzed and the various factors determining the detected fluorescence signal are studied numerically. Relative OH concentration can be explained on an absolute scale by calibrating the optical collection constant through using Rayleigh scattering or comparing with calculation results in some simple cases. Consequently, an important effort is devoted to the evaluation of the influences of those parameters on the fluorescence signal. The physical models of spectral absorption line shape broadening, collisional transitions and quenching at elevated pressure are built.
For the typical LIF measurement setup, the relation between the intensity of fluorescence signal and the absolute concentration of OH radicals is described clearly. The fine energy level structure of the OH molecule is illustrated to understand the process of laser-induced fluorescence and others in the case without radiation, which include collisional quenching, rotational energy transfer (RET) and vibrational energy transfer (VET). The quenching, VET and RET-rate are all induced by the collision partner species, and the rates are determined by the species and the number densities, therefore the quenching and other collisional effects become stronger at the same mole fraction distribution under higher pressure.
Based on the detailed physical models built in the present work, some numerical results are achieved by simulations in order to evaluate the fluorescence yield. Firstly, the broadening mechanisms of absorption line are studied and line shape Voigt is accepted in our case of high pressure. The spectral overlap values in a temperature range from 800 K to 2000 K and pressure range from atmosphere to 1.0 MPa are calculated for calibration. Secondly, for a given transition, the population distribution of an ensemble of molecules, which influences the fluorescence signal, is given by the Boltzmann equation. At low temperatures, most of the molecules will be in the ground energy state, while high temperatures will lead to an increase in the probability of finding molecules at higher energy levels.
Most importantly, the simulation of OH LIF with excitation from both excitation lines P1(7) (0–1) and Q1(6) (0–0) is carried out by utilizing the LASKIN program. The transition-dependent quantum yields with different temperatures and various pressures in our certain case are estimated numerically by simulating the time-integrated population distributions. It can be observed that the quenching is a main factor to influence the fluorescence yield especially at elevated pressure significantly. It is worth mentioning that the rates of quenching, VET and RET depend on not only pressure and temperature values, but also the species of production as collision partners.
These results are useful for understanding the real physical processes in OH-LIF technique and finding a way to calibrate the signal for quantitatively measuring OH concentration. The next step of this study is to investigate the fundamental combustion case such as laminar premixed flame in high pressure, develop other quantitative measurement skills and compare experimental results with those from the modeling in order to assess the ability for kinetic mechanisms to modify the parameters of the calibration method. This topic will be studied in the future.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 | |
24 | |
25 | |
26 | |
27 | |
28 | |
29 |