Broadband supercontinuum laser absorption

Letter

Vol. 40, No. 13 / July 1 2015 / Optics Letters

3141

Broadband supercontinuum laser absorption spectrometer for multiparameter gas phase combustion diagnostics NIELS GÖRAN BLUME

AND

STEVEN WAGNER*

High Temperature Process Diagnostics Group, Technische Universität Darmstadt, 64289 Darmstadt, Germany *Corresponding author: [email protected]‑darmstadt.de Received 2 April 2015; revised 2 June 2015; accepted 12 June 2015; posted 12 June 2015 (Doc. ID 237395); published 29 June 2015

We report on the development and application of a broadband absorption spectrometer utilizing a pulsed supercontinuum laser light source and dispersion compensating fiber with a single-pass absorption path to obtain absolute methane mole fractions in a laminar nonpremixed CH4 ∕air flame supported on a Wolfhard–Parker burner. The basic principle of supercontinuum broadband absorption spectroscopy (SCLAS) provides advantageous means of combustion diagnostics since the broad spectral coverage allows for use in high-pressure high-temperature environments. Furthermore, a previously validated tunable diode laser absorption spectroscopy fitting algorithm was applied to the recorded spectra and found to be applicable to SCLAS measurements as well, by comparison of fitted methane gas concentrations to reference measurements on the Wolfhard–Parker burner. The spectrometer reached spectral resolutions of up to 0.152 cm−1 , while providing a spectral coverage of over 110 cm−1 , with an absorption path length of only 41 mm. First measurements of absolute CH4 mole fractions showed the suitability of SCL-based spectroscopy for combustion diagnostics with short absorption path lengths in the nIR spectral region. Here, we achieved in-flame methane mole fraction resolutions of 3 %Vol. (1210 ppm·m) and optical resolutions of up to 1.1 × 10−2 . Based on this first validation, this method can now be extended to other species and combustion parameters such as temperature and pressure. © 2015 Optical Society of America OCIS codes: (120.1740) Combustion diagnostics; (260.2030) Dispersion; (120.6200) instrumentation.

Spectrometers

and

spectroscopic

http://dx.doi.org/10.1364/OL.40.003141

Supercontinuum broadband laser absorption spectroscopy (SCLAS) is a fairly recent development in the field of broadband spectroscopy. Based on a theoretical proposal by Orofino and Unterleitner in 1976 [1], only the combination of recent 0146-9592/15/133141-04$15/0$15.00 © 2015 Optical Society of America

developments in the field of high-speed high-bandwidth optical detection and supercontinuum light sources have allowed for application of this technique in general and specifically in the field of gas phase combustion research [2–6]. Due to its broad spectral coverage and high temporal resolution, SCLAS is a promising technique for diagnostics in the context of newly developed combinatory combustion technologies, such as homogeneous charge compression ignition or spark ignition direct injection. All these approaches to pollutant reduction share the common need for extensive simultaneous measurements of multiple combustion parameters with high temporal resolution in high-pressure and high-temperature environments, which SCLAS is able to provide. SCLAS employs recently available supercontinuum laser light sources (SCL), which emit pulsed laser light, covering the spectrum from 400 to 2400 nm at repetition rates of up to 80 MHz. The emitted pulses of SCL sources are fairly short pulses (∼100 ps), which are afterward dispersed in time to allow for frequency-resolved intensity detection. The recorded absorption spectra can be fitted against spectral databases for determination of species concentrations and other combustion parameters such as temperature. So far, SCLAS in the combustion context has been applied almost exclusively in multipass or long-pass setups. Sanders [2] used an 1840 mm absorption cell with pressures up to 10 bar to resolve spectral feature of H2 O, CO2 , and C2 H2 with a resolution of up to 1.5 cm−1 . Werblinski et al. [3] used SCLAS in a McKenna-type flat flame burner with a multipass setup and an effective absorption path length of at least 600 mm. Previous experiments of the same group [4] were performed using a pressurized flow cell with an absorption length of 300 mm and temperatures up to 900 K. For the measurements in Kaminski et al. [5], an absorption path length of 2800 mm was used to measure ambient vapor H2 O species concentrations. In Hult et al. [6], a multipass setup of 105 mm was used for an unburned mixture of 20% CH4 and air at room temperature. Previous research has been so far focused on basic studies, including very long absorption paths and unburned gas mixtures at ambient temperature and slightly elevated pressures, but is still lacking validation for

3142


in-flame measurements. Since, especially in high-pressure and high-temperature environments, the broad spectral width of SCLAS is beneficial for the fitting process and the quality of the results [7–9], extending the SCLAS diagnostic principle toward these applications is most promising for in situ diagnostics in enclosed combustion. A necessary step in this way is the validation of the SCLAS for quantitative in-flame measurements with single-pass absorption paths and path lengths on the order of cylinder diameters (∼80 mm) and shorter. In this Letter, we demonstrate the applicability of SCLAS in the context of single-pass spatially resolved absolute mole fraction measurements of combustion species within a flame, which are validated against reference measurements by multiple other diagnostic techniques. To achieve this, we utilized an atmospheric, laminar, nonpremixed CH4 ∕air flame consisting of two flame sheets on a Wolfhard–Parker burner (WHP) [10]. The WHP burner enables the quantitative validation of new measurement methods for combustion diagnostics, as it has been extensively characterized with spatially resolved measurements of more than 17 individual species, temperatures, and flow-field with a high number of independent measurement techniques. Measurements, which include, for example, the species CH4 , CO, and CO2 as well as temperatures and flow velocities, can be found in [10]. We recorded broadband methane spectra ranging over 110 cm−1 with an absorption path of 41 mm, which were, consequently, fitted against spectral database HITRAN 2012 [11] and evaluated against reference measurements for this burner [10,12,13]. For the experiments, the setup shown in Fig. 1 was used. As a laser light source, a SuperK EXTREME EXW-12 by NKT Photonics A/S was utilized. This source uses a pump wavelength of 1064 nm and generates the supercontinuum laser light in a photonic crystal fiber with a spectrum of 400– 2400 nm. To select a specific region of this spectrum for the experiment, spectral filtering was performed using the SuperK SELECT by NKT Photonics A/S, which provides acoustooptical tunable filters (AOTFs) to select up to eight spectral windows (FWHM 10–18 nm). The emitted light is coupled into a single-mode fiber and chromatically dispersed in a dispersion compensating module (OFS Fitel Denmark ApS,

Fig. 1. Experimental setup for SCLAS measurements on an atmospheric, laminar, nonpremixed CH4 ∕air flame consisting of two flame sheets stabilized using a wire screen chimney (Wolfhard–Parker burner [10]). Supercontinuum laser pulses, emitted by a supercontinuum laser light source (SCL), are spectrally filtered and dispersed chromatically in a dispersion compensating module (DCM). After passing through the flame, the light is coupled into a multimode fiber (MMF, 50 μm) and detected by a photodiode (PD) and an oscilloscope (DAQ).

Letter LLMicroDK:1500, −1500 ps∕nm at 1550 nm). The dispersed pulses are then modulated by the absorption features of the gases present in the flame of the WHP burner. After the transmission across the flame, the light is coupled into a multimode fiber (step-index fiber, 50 μm core diameter) and detected by a fiber-coupled high-speed high-bandwidth InGaAs photodiode (Newport Inc., New Focus 1544-B-50, 12 GHz bandwidth, −1000 V∕A transimpedance amplifier gain). The electrical output of the photodiode is recorded with an oscilloscope (LeCroy Corporation, SDA 820Zi-A, 20 GHz bandwidth, 40 GS∕s). The complete spectroscopic system has been characterized and verified regarding its spectral behavior with a grating spectrometer (AOTF, FWHM, etc.). Furthermore, it is necessary for the determination of the optical resolution to measure the dispersion of the DCM, since the dispersion is directly related to the spectral resolution of the system by the following equation [2]: Δν τ · ν_ :

(1)

In Eq. (1), Δν (nm) is the spectral resolution, τ the FWHM of the shortest optical pulse (ps), resolvable by the photodiode, and ν_ the scan rate (nm/ps) as a parameter of the dispersion fiber used in the setup. Based on the rise-time of 38 ps of the currently used photodiode, τ can be estimated to be 76 ps. The dispersion of the DCM was measured using a novel approach using the available AOTF module (described in detail in a future publication). Based on the measured delay (ps) of different wavelengths (nm), the dispersion can be determined as the first derivative of that delay, given that the dispersion (ps/nm) is the slope of the delay for a given length of fiber [14]. For different lengths of dispersion compensating fiber, the dispersion must be scaled according to length. For the DCM used in these experiments (9.382 km of fiber), the measured dispersion curve is shown in Fig. 2. Based on the measurements of the dispersion, the optical resolution of the system can be estimated according to Eq. (1) for the relevant wavelength range of 1650–1700 nm: Δν τ · ν_ 76 ps · 0.002 cm−1 ∕ps 0.152 cm−1 :

(2)

This optical resolution is equivalent to ∼0.04 nm at 1675 nm. To allow for the recording of broadband spectra of methane in the WHP burner, this resolution is more than sufficient and could even be decreased for high-pressure

Fig. 2. Measurement of the dispersion of the DCM used in the experiments. The delay for multiple wavelengths was recorded and, subsequently, fitted with a fourth-order polynomial function, which was thereupon used for calculation of the dispersion as the first derivative of the fourth-order polynomial function.

Letter environments. Since the optical resolution is directly related to the applied dispersion, a reduction in resolution is achieved by a reduction in dispersion. Such a reduction, on the other hand, would allow for an increase in repetition rate, since the time duration of a dispersed pulse is reduced as well, and no overlap of pulses is present at higher repetition rates. For the experiments, a modified WHP burner design was used employing nitrogen purge slots to suppress end-flame effects [13]. The modified WHP burner consists of a center slot of methane (length: 41 mm) with two adjacent slots of air. The two diffusion flames are stabilized against a wire screen chimney, and the complete setup is mounted on a two-axis motion stage to allow for precise positioning of the flame. The WHP burner was operated with an exit velocity for methane of 11 cm∕s and for both air slots of 22 cm∕s as well as for nitrogen of 22 cm∕s. Using the motion stage, measurements were taken at different transversal positions as well as at different heights above the burner surface. Based on the HITRAN 2012 spectral database [11], multiple multiline Voigt shapes (over 1000 absorption lines) were fitted against the measured spectra using a previously validated Levenberg–Marquardt algorithm [13] to minimize the residual. A higher-order polynomial function was fitted to the background as a reference intensity. The obtained direct absorption spectra were evaluated using the extended Lambert–Beer law, as shown in Eq. (3): Iλ I 0 λ · exp−ST · gλ − λ0 · N · L · Trt · Et: (3) In Eq. (3), I λ describes the light intensity after absorption, I 0 λ the light intensity before absorption, ST the line strength, gλ − λ0 the line form function, N the number density, L the absorption path length, Trt the broadband transmission fluctuation (i.e., particles), and Et the broadband background emission. Given the near static conditions of this experiment, Trt is almost identical to one. Using Eq. (3) and previously performed flame temperature measurements as well as pressure measurements in combination with the tabulated line strength values and the ideal gas law, the concentration of methane was determined and compared to reference measurements found in [10]. Figure 3 shows a raw measurement of methane in the WHP burner. For this measurement, 9870 single pulses (2 MHz repetition rate) of the SCL were averaged, resulting in a spectral sampling rate approximately 200 Hz (5 ms). After application of time-wavelength correlation (nonlinear dispersion of DCM), application of background correction and fitting of appropriate spectral features, an exemplary result spectrum is shown in Fig. 4. The estimated fractional absorption from the local residual of 1.133 × 10−2 gives a signal to noise ratio of 24 at 489 K (3 %Vol. or 1210 ppm·m, [15]). In advanced processing algorithms, advanced combustion parameters such as temperature and species concentration of additional combustion educts and products (i.e., H2 O) will be derived from single SCLAS broadband spectra as opposed to the need to combine multiple narrowly tuned laser sources for multiparameter measurements. SCLAS might suffer from slightly less optical resolution but can overcome this drawback by providing robust and versatile broadband measurements


3143

Fig. 3. Raw measurement using the SCLAS setup of methane gas concentration in WHP burner at transversal position 0.5 mm and height 7 mm above burner surface.

that allow for a unified approach to multiparameter combustion measurements. By combining the recorded spectra of multiple positions, a spatially resolved methane profile across the flame sheets can be measured. In Fig. 5, the recorded spectra at the lateral positions 0.5 mm (489 K), 1.0 mm (571 K), 1.5 mm (685 K), and 2.0 mm (822 K) are shown at their corresponding temperature. To evaluate the performance of SCLAS, the measured spatial concentrations for methane according to [10] were compared against our SCLAS measurements, as shown in Fig. 6. It can be derived from the comparison that the species concentrations of methane obtained with SCLAS are consistent with the concentrations of Norton et al. [10], thus validating SCLAS as a combustion diagnostics tool for absolute species concentration measurements in flames. Overall, in this Letter we were able to validate SCLAS as a versatile and feasible tool for quantitative in situ combustion diagnostics. It was shown that SCLAS is able to correctly determine species concentrations inside a laminar flame.

Fig. 4. Exemplary methane spectrum as recorded with the SCLAS setup at lateral WHP position of 0.5 mm with a temperature of 489 K. For the fit 1001 methane lines from HITRAN 2012 [11] (Voigt line-shape functions) were fitted simultaneously with a higherorder polynomial reference intensity I 0, as described in [13].

3144


Fig. 5. Methane mole fraction measurements at 7 mm above burner surface and different transversal positions. Different lateral positions are plotted at their corresponding temperature, which increases from the center outward from 489 to 822 K.

Letter complex combustion scenarios, including high-pressure chemical reactors or internal combustion engines. Such environments provide new challenges for the SCLAS system, since the pressure levels are significantly higher as well as optical accessibility being limited. Regarding the data evaluation, we see great potential for specifically designed SCLAS fitting algorithms that fully utilize the spectral information contained in SCLAS spectra. As mentioned earlier, the temporal resolution can be increased in high-pressure environments as well, since reduced requirements for optical resolution allow for an increase in pulse repetition rate. Furthermore, an optimization of the SCL source itself for the nIR range would be beneficial to the signal quality and the overall spectrometer performance, since the noise level inherent to the supercontinuum generation could be reduced when reducing the spectral range to the nIR [16]. Funding. Technische Darmstadt) (UNICO).

Universität

Darmstadt

(TU

REFERENCES

Fig. 6. Comparison of spatially resolved methane concentration comparing data from Norton et al. [10] and our recent SCLAS measurements.

Additionally, the setup presented here allowed for spatially resolved quantitative in-flame gas concentration measurements with an absorption path length of only 41 mm in a single-pass setup. Based on these results, SCLAS can be considered as a validated quantitative diagnostic system for in-flame combustion diagnostics. Given the fully fiber-coupled nature of SCLAS, this diagnostic technique is, furthermore, well suited for use in harsh combustion environments as a broadband diagnostic tool. In the future, we will continue to apply SCLAS to the WHP burner setup and record, as well as evaluate, complete spatial profiles to further characterize the behavior and applicability of SCLAS in combustion environments. Additionally, we currently investigate the applicability of this technique in more

1. T. A. Orofino and F. C. Unterleitner, Appl. Opt. 15, 1907 (1976). 2. S. Sanders, Appl. Phys. B 75, 799 (2002). 3. T. Werblinski, F. Mittmann, M. Altenhoff, T. Seeger, L. Zigan, and S. Will, Appl. Phys. B 118, 153 (2015). 4. T. Werblinski, S. R. Engel, R. Engelbrecht, and S. Will, Opt. Express 21, 13656 (2013). 5. C. F. Kaminski, R. Watt, A. Elder, J. Frank, and J. Hult, Appl. Phys. B 92, 367 (2008). 6. J. Hult, R. S. Watt, and C. F. Kaminski, Opt. Express 15, 11385 (2007). 7. C. S. Goldenstein, R. M. Spearrin, I. A. Schultz, J. B. Jeffries, and R. K. Hanson, Meas. Sci. Technol. 25, 55101 (2014). 8. L. A. Kranendonk, J. W. Walewski, T. Kim, and S. T. Sanders, Proc. Combust. Inst. 30, 1619 (2005). 9. J. Seitzman, R. Tamma, and B. Scully, “Broadband infrared sensor for active control of high pressure combustors,” in 36th AIAA Aerospace Sciences Meeting and Exhibit, Atlanta, Georgia, 1998. 10. T. S. Norton, K. C. Smyth, J. H. Miller, and M. D. Smooke, Combust. Sci. Technol. 90, 1 (1993). 11. L. S. Rothman, I. E. Gordon, Y. Babikov, A. Barbe, D. Chris Benner, P. F. Bernath, M. Birk, L. Bizzocchi, V. Boudon, L. R. Brown, A. Campargue, K. Chance, E. A. Cohen, L. H. Coudert, V. M. Devi, B. J. Drouin, A. Fayt, J.-M. Flaud, R. R. Gamache, J. J. Harrison, J.-M. Hartmann, C. Hill, J. T. Hodges, D. Jacquemart, A. Jolly, J. Lamouroux, R. J. Le Roy, G. Li, D. A. Long, O. M. Lyulin, C. J. Mackie, S. T. Massie, S. Mikhailenko, H. S. P. Müller, O. V. Naumenko, A. V. Nikitin, J. Orphal, V. Perevalov, A. Perrin, E. R. Polovtseva, C. Richard, M. A. H. Smith, E. Starikova, K. Sung, S. Tashkun, J. Tennyson, G. C. Toon, V. G. Tyuterev, and G. Wagner, J. Quant. Spectrosc. Radiat. Transfer 130, 4 (2013). 12. J. H. Miller, S. Elreedy, B. Ahvazi, F. Woldu, and P. Hassanzadeh, Appl. Opt. 32, 6082 (1993). 13. S. Wagner, B. T. Fisher, J. W. Fleming, and V. Ebert, Proc. Combust. Inst. 32, 839 (2009). 14. L. G. Cohen and C. Lin, Appl. Opt. 16, 3136 (1977). 15. F. Stritzke, O. Diemel, and S. Wagner, Appl. Phys. B 119, 143 (2015). 16. S. Dupont, Z. Qu, S.-S. Kiwanuka, L. E. Hooper, J. C. Knight, S. R. Keiding, and C. F. Kaminski, Laser Phys. Lett. 11, 75601 (2014).