The objective of this study is to understand the effects of fuel volatility on soot emissions. This effect is investigated in two experimental configurations: a jet flame and a model gas turbine combustor. The jet flame provides information about the effects of fuel on the spatial development of aromatics and soot in an axisymmetric, co-flow, laminar flame. The data from the model gas turbine combustor illustrate the effect of fuel volatility on net soot production under conditions similar to an actual engine at cruise. Two fuels with different boiling points are investigated: n-heptane/n-dodecane mixture and n-hexadecane/n-dodecane mixture. The jet flames are nonpremixed and rich premixed flames in order to have fuel conditions similar to those in the primary zone of an aircraft engine combustor. The results from the jet flames indicate that the peak soot volume fraction produced in the n-hexadecane fuel is slightly higher as compared to the n-heptane fuel for both nonpremixed and premixed flames. Comparison of aromatics and soot volume fraction in nonpremixed and premixed flames shows significant differences in the spatial development of aromatics and soot along the downstream direction. The results from the model combustor indicate that, within experiment uncertainty, the net soot production is similar in both n-heptane and n-hexadecane fuel mixtures. Finally, we draw conclusions about important processes for soot formation in gas turbine combustor and what can be learned from laboratory-scale flames.

Introduction

The increasing demand for aviation fuels, unstable fuel prices, stricter environmental rules, and need for energy security have led to an increased interest in alternative jet fuels. Extensive research in alternative jet fuels has led to approval of Fisher-Tropsch fuel (Annex 1 ASTM D7566) and hydroprocessed renewable jet (Annex 2 ASTM D7566) to be used as blended fuels with current petroleum-derived jet fuels [1]. The components of the new alternative fuels may have very different molecular weights and structures than those of the currently approved fuels [2]. This variation in the fuel molecular-structure distribution generates new interest in understanding the effects of fuel composition on operability of current and future gas turbine engines.

The research efforts in the current project aim to aid in developing fundamental understanding of the effects of fuel molecular structure on emissions. To achieve this goal over a wide range of experimental conditions, a series of experiments and companion simulations were planned that progress in their level of complexity from simple laboratory flames to sector rigs [3]. The study at Penn State involves understanding the chemical effect of fuel molecular structure on aromatics and soot in jet flames, and testing the effect of fuel volatility on emissions in a model gas turbine combustor. The chemical effects of fuel molecular structure on soot have been discussed in Refs. [4] and [5].

The primary aim of the work presented here is to understand the effect of fuel volatility on emissions in an actual engine-like environment. The hypothesis of the study is that fuels with varying volatility will affect fuel droplet vaporization processes and hence, on fuel–air mixing and emissions. A fuel with a higher boiling point, i.e., a heavier hydrocarbon fuel, is expected to produce larger initial droplet sizes than a lighter hydrocarbon [6,7]; if the droplet vaporization time is longer than the convention time of the droplet from the injector exit to the flame front, the droplets will burn in a diffusion flame at a stoichiometric condition. As a result of this potentially incomplete vaporization, a heavier hydrocarbon fuel will likely have higher soot emissions as compared to a fuel with a lower boiling point.

In addition to testing the fuels in the model gas turbine combustor, the chemical effect of the two fuel mixtures on the aromatics and soot is investigated in a laminar, co-flow, axisymmetric burner under nonpremixed and rich premixed condition at atmospheric pressure. The aromatics are not only soot precursors, but are also classified as hazardous air pollutants themselves [8]. Nonpremixed flames play an important role in practical combustion systems; these flames display all major soot formation processes, which include soot nucleation, growth, and burnout processes, similar to those in real engines. The rich premixed conditions are selected in order to mimic conditions similar to those in the primary zone of an aircraft engine combustor. The n-alkane fuel mixtures studied here include n-heptane/n-dodecane (C7) and n-hexadecane/n-dodecane (C16) fuels. The fuel selection and experimental conditions are discussed in the Experimental Setup section. The C7 results are also discussed in more detail in Refs. [4] and [5]; we use them here as a means to understand fuel molecular weight and volatility effects, which are not discussed in Refs. [4] and [5].

Previous studies have investigated the soot development of higher-hydrocarbon fuels [911]. The n-hexadecane fuel, studied here, is representative of the molecular weight of components of jet fuel [12] and diesel fuel [13]. For example, Eddings et al. [14] studied a six component surrogate (Hex-11 and Hex-12) as a possible surrogate to match the volatility (the boiling point distribution) and sooting tendency of JP-8. The Hex-11 and Hex-12 had 5 mol % and 12 mol % of n-hexadecane as one of the surrogate components, respectively. The sooting limit and polycyclic aromatic hydrocarbons (PAH) formation of n-hexadecane has been studied in a microflow reactor at equivalence ratios of 1.5–4.5, and 4, respectively [15]. Pastor et al. [10] studied mixtures of n-decane and n-hexadecane fuels in a diesel engine. They found an increase in the n-hexadecane concentration resulted in an increase in the amount of soot formed in the flame. Douce et al. [16] studied soot formation from n-hexadecane in a shock tube. However, to the best of our knowledge, no study has reported aromatics and soot formation in both nonpremixed and rich-premixed jet flames involving n-hexadecane as a fuel. The results from the jet flame studies provide important validation-quality data for soot models.

The goal of this work is to understand the impact of fuel molecular weight on soot, and to compare laboratory flame results and the behavior of these fuels in actual gas turbines. The remainder of the paper presents the experimental configurations used in this study, as well as the measurement techniques. We provide details of the soot formation and oxidation processes in the laboratory flame, with quantitative measurements of soot volume fraction for the two fuels. Finally, we compare these soot measurements from jet flames to the model gas turbine combustor results. The results presented here on soot volume fraction in jet flames can be used to test chemical mechanisms and soot models for n-hexadecane in a well-characterized configuration. In addition, the results obtained from the model gas turbine combustor provide validation data for soot emissions in an environment similar to a gas turbine engine.

Experimental Setup and Techniques

Laminar Flame Experiment.

The effect of fuel molecular weight on the spatial development of aromatics and soot is investigated in an axisymmetric, co-flow, laminar flame configuration at atmospheric pressure. The burner design is similar to the co-annular burner used by Santoro et al. [17]. The burner design is modified to operate it using prevaporized higher hydrocarbon fuels. The fuel tube is 4 mm above the exit plane of the burner. The detailed design of the burner and boundary conditions is discussed in Refs. [4] and [5]. The schematic of the burner and the burner system is shown in Fig. 1. The fuel is delivered to the two vaporization units using two different pumps: an Isco syringe pump supplies the n-dodecane fuel, and a Nexus 6000D syringe pump supplies the “second” fuel of the binary mixture. The fuel vaporization unit consists of a Swagelok tee filled with fiberglass.

Fig. 1
Schematic of burner and the burner system
Fig. 1
Schematic of burner and the burner system
Close modal

The vaporized fuel is carried by nitrogen gas from the vaporization unit to the burner exit. The energy for fuel vaporization is provided by several heating tapes wrapped around the outside surface of the tubing system. The air inlet, which is only used to study premixed flames, is situated after the two vaporization units. The temperature of the system is maintained at 300 ± 10 °C to prevent fuel condensation in the system. The gas chromatography analysis of the recondensed fuel from the tube exit indicates absence of fuel pyrolysis in the tubing system. The further details of the burner system are described in Ref. [4].

The co-flow air is set at a constant flow rate of 4 scfm. The flame is shielded from the surroundings using a brass chimney (400 mm long) with screens on the top. Four slots are machined at symmetrical locations on the chimney to provide access for the measurement techniques. A stepper motor (Velmex PK266) and a Velmex bislide, with step resolution of 6.50 μm and 5.00 μm, are used to move the burner in the horizontal and vertical directions, respectively.

The two-dimensional (2D) soot volume fraction in the jet flames is obtained using laser-induced incandescence (LII) [18]. The LII is calibrated using laser extinction (LE). The 2D aromatics profile is obtained using laser-induced fluorescence (LIF). The optical setup, data collection, and analysis of measurement techniques used in this study are identical to that described in Refs. [4] and [5]. A brief description of the technique is provided here. The LII is obtained using a Nd:YAG laser (Spectra Physics GCR 270-10, 10 Hz, 7 ns, 532 nm). The laser beam is expanded into a laser sheet and focused on the flame center using a combination of a plano-convex cylindrical lens and a convex lens. The laser fluence is maintained at 0.3 J/cm2 to be in the “plateau” region of LII [18] for the LII experiment. In the LIF experiment, a 266 nm, 585 μJ laser pulse from the Nd:YAG laser is used for exciting the aromatic molecules.

The LII and LIF signal is recorded using a Princeton Instruments ICCD camera with a sensor size of 512 × 512 pixels using an f/2.8 105 mm UV lens. A spatial resolution of 0.1 mm/pixel of camera is achieved using the optical setup. The LII signal is collected using a narrow band interference filter centered at 430 nm with 10 nm FWHM [19]. The LII signal collection duration is 100 ns, after a delay of 20 ns following the laser pulse. For LIF signal collection, camera gate-width is 170 ns; the camera gate is opened 20 ns before the laser pulse. The study by Zizak et al. [20] found that the fluorescence wavelength shifts to longer wavelength as the aromatic molecule size increases. This property is implemented in the current study to identify two different classes of aromatics. The signal for the small (one and two ring) aromatics is obtained using two filters (WG320 and UG11), which have transmission in the wavelength band 320–380 nm. The signal for large (three–five ring) aromatics is obtained using a filter (500FS20-50), which has transmission in the wavelength band 488–514 nm. The small and large aromatics are further discussed in Ref. [4]. The LII and LIF signals are corrected for background noise, pixel-to-pixel responsivity of the detector, and variation in the laser intensity across the laser sheet.

A calibration factor of 3350 LII counts equals 1 ppm of soot is obtained using LE. A value of m̃=1.570.56i is used as the refractive index of soot to be consistent with previous published work [19,21]. The estimated uncertainty in soot volume fraction measurement through laser extinction and calibration of LII data is ±15%. The measurement uncertainty are discussed in detail in Refs. [4] and [5]. The LIF measurements are presented as raw counts, because the fluorescence in each wavelength band is collected from a number of aromatic species; the fluorescence quantum yield and temperature dependency of absorption cross section of each aromatic species are not well characterized. The temperature measurement in the jet flame is obtained using a 125 μm, R-type thermocouple and is corrected for radiative heat transfer loss [4].

Model Gas Turbine Combustor.

In practical devices like gas turbine engines, the physical properties of the fuel can have an important impact on the engine's performance. Experiments in a model gas turbine combustor can help to understand the fuel volatility effect on emissions in a realistic combustion device. The schematic of the model gas turbine combustor is shown in Fig. 2. It consists of two sections: the inner section where the combustion occurs and the outer co-flow section for regenerative cooling.

Fig. 2
Schematic of model gas turbine combustor [7]
Fig. 2
Schematic of model gas turbine combustor [7]
Close modal

The compressed air from the air delivery system enters the outer co-flow section of the combustor. The passage through this section preheats the air, thereby cooling the quartz window on the combustor chamber and removing heat from the exhaust stream. The preheated air then enters an electric heater, where it can be heated to a maximum temperature of 800 K. The heated air enters the combustion chamber after passing through a 7.39 mm inlet venture and a 45 deg flat vane swirler. The schematic of the air swirler is shown in Fig. 3. The geometric swirl number based on Ref. [22] is 0.77.

Fig. 3
Schematic diagram of the air swirler [7]
Fig. 3
Schematic diagram of the air swirler [7]
Close modal

The liquid fuel injector used in this study is the same as used in Ref. [21], shown in Fig. 4. It consists of two concentric stainless steel tubes. The liquid fuel under pressure from the fuel delivery system enters the inner tube of the injector and is forced out through the 0.005 in diameter, 12 radial injector holes. The injection point is 36.8 mm upstream of the dump plane. A part of the primary air is diverted to flow through the annular region surrounding the inside tube of the injector to prevent the fuel from vaporizing inside the inner tube. The annulus is tapered at the end to create a shearing force to break the liquid fuel into small droplets. A stainless steel bluff body is attached at the end of the injector for flame anchoring.

Fig. 4
Schematic diagram of the fuel injector [7]
Fig. 4
Schematic diagram of the fuel injector [7]
Close modal

The fuel–air mixture then flows into the 45 mm diameter, 307 mm long combustion chamber. The combustion products exit the combustor through a choked nozzle. The choked inlet and outlet help to isolate the impact of combustion instabilities on inlet air or fuel flow rate and inhibit pressure oscillations from the exit section to propagate to the combustion zone. The exit nozzle is maintained at a temperature of 480 K using air and water cooling. Two quartz windows of 12.5 mm in diameter are installed at the exit of the combustor to provide optical access for laser extinction technique. The quartz windows are purged with nitrogen gas to maintain the temperature of the windows and prevent soot deposition on them.

Line-of-sight soot volume fraction is measured at the exit of the combustor, 24.8 cm downstream of the dump plane. The LE optical setup used for the model gas turbine combustor in this study is identical to the one used in Ref. [21]. The output signals from both photodiodes are converted to voltages using a trans-impedance amplifier and then fed into a lock-in amplifier to reduce the noise from the surrounding environment. The data are collected by a National Instruments data acquisition system (NI-PCI-6110) at a sampling rate of 10 samples/s for 4 s.

Fuel Selection and Experimental Matrix

The fuel matrix, shown in Table 1, is comprised of two binary mixtures: n-heptane/n-dodecane and n-hexadecane/n-dodecane, a high boiling point fuel; n-hexadecane has a boiling point of 287 °C as compared to 216 °C for n-dodecane and 98 °C for n-heptane. The two binary mixtures are chosen such that the carbon fraction of n-dodecane and the added hydrocarbon is constant.

Table 1

Fuel matrix

NameComponentVolume fraction (%)H2 content (%)H/C ratio
C7n-Heptane31.315.62.20
n-Dodecane68.7
C16n-Hexadecane28.415.42.15
n-Dodecane71.6
NameComponentVolume fraction (%)H2 content (%)H/C ratio
C7n-Heptane31.315.62.20
n-Dodecane68.7
C16n-Hexadecane28.415.42.15
n-Dodecane71.6

The flames investigated are nonpremixed (jet equivalence ratio Φjet = Inf) and rich premixed (Φjet = 24 and 6). The total carbon flow rate is held constant for all the flames at 0.011 mol/min. The carbon flow rate is defined as the summation of the product of each component flow rate and the corresponding number of carbon atoms. The carbon flow rate is fixed to match carbon flow from the base flame, a 75% n-dodecane and 25% m-xylene mixture [4]. The flow rate of nitrogen, 0.2 slpm, is used to carry the vaporized fuel to obtain a stable flame.

The effect of fuel volatility on soot emissions, for the fuels listed in Table 1, is investigated in a model gas turbine combustor. This experiment replicates the complexity of bulk mixing, turbulent mixing, and spray as is present in a real combustor operating on liquid fuels. The combustor has been used in previous studies on combustion instability, soot, and emissions [21,23]. In addition to the two binary fuel mixtures, JP-8 (POSF 5699) fuel is tested to compare the results obtained in this study to previous published data [21], which involved testing of soot emission from JP-8 under the same experimental conditions. The experimental conditions for the model gas turbine combustor are outlined in Table 2. The air mass flow rate and inlet temperature are kept constant at 32 g/s and 560 ± 10 K for all the fuels tested. The inlet air temperature is measured at the location upstream of the critical venturi.

Table 2

Experimental conditions for model gas turbine combustor

Air mass flow rate32 g/s
Inlet air temperature560 ± 10 K
InjectorMP, 12 hole, 0.005 in diameter
Swirl angle45 deg
Chamber pressure0.51 MPa (5 atm)
Chamber length307 mm
Air mass flow rate32 g/s
Inlet air temperature560 ± 10 K
InjectorMP, 12 hole, 0.005 in diameter
Swirl angle45 deg
Chamber pressure0.51 MPa (5 atm)
Chamber length307 mm

The soot production can be impacted by the fuel droplet size and droplet evaporation time. The initial droplet size depends on the type of injector used [6]. For a particular size droplet, the droplet evaporation time would depend on the initial temperature of the droplet and the adiabatic flame temperature. A higher inlet temperature would increase the adiabatic flame temperature, which would result in faster droplet evaporation. The faster evaporation results in more premixed combustion, leading to less soot production. The inlet air temperature is held constant for all the fuels in order to make the comparison between different fuels possible. The adiabatic temperature of both fuels n-heptane and n-hexadecane is quite similar, close to 2265 K [24]. Hence, we expect the flame temperature to be similar for both the fuels at a particular equivalence ratio. The global equivalence ratio (Φglobal) is varied from 0.9 to 1.8 by changing the fuel flow rate. An exit nozzle of 10.2 mm diameter is used to choke the flow to create a mean chamber pressure of 0.51 MPa. The combustor chamber pressure is held steady at 5 atm, within measurement uncertainty, for Φglobal from 0.9 to 1.8. The flames are both statically and dynamically stable under these operating conditions. The fuels used had purity of >99% and 95% for the jet flame experiment and the model combustor, respectively.

Results and Discussion

Laminar Flame: Visible Flame Heights.

Figure 5 shows the visible flame images for all the flames obtained using a Canon camera at same operating settings (ISO = 80, f-stop, and exposure time). From left to right, the flames are C7Φjet = Inf flame, C16Φjet = Inf flame, C7Φjet = 24 flame, C16Φjet = 24 flame, C7Φjet = 6 flame, and C16Φjet = 6 flame.

Fig. 5
Visible flame images obtained at same camera operating condition: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Fig. 5
Visible flame images obtained at same camera operating condition: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Close modal

In this work, the flame height is defined as the distance from the tube exit to the flame tip. The flame tip is obtained from the visible flame images, shown in Fig. 5. The flame tip has a fluctuation of ±2 mm over the duration of the experiment. At a particular Φjet, within measurement uncertainty, the flame heights are similar for both the fuels. As the level of premixing is increased, the flame height decreases. The flame heights decrease by approximately 18% for both the fuels from the Φjet = Inf condition to Φjet = 6 condition. In these six buoyancy-driven flames, the similar flame height, at a particular Φjet, indicates that the residence time in the flame for soot formation and oxidation is similar for the two fuels, if the soot nucleation is at similar height above burner (HAB).

Laminar Flame: Development of Aromatics.

Figure 6 shows the 2D spatial development of small aromatics for the two fuels at three different Φjet. In all the 2D LIF and soot volume fraction figures comparing the six flames, Figs. 6(a)6(f) show the 2D profiles for C7Φjet = Inf flame, C16Φjet = Inf flame, C7Φjet = 24 flame, C16Φjet = 24 flame, C7Φjet = 6 flame, and C16Φjet = 6 flame, respectively. The minimum and maximum value of the scale are shown on the top of each figures, for example, [0 27,000], [0 15,000], and [0 6] means the lowest value of the scale is 0; 27,000, 15,000, and 6 corresponds to the maximum value of the scale. At a particular Φjet, the 2D profiles are plotted on the same scale for the two fuels. However, the scale is different for the three Φjet to show the effect of premixing on the spatial development of small aromatics, large aromatics, and soot distinctly.

Fig. 6
Two-dimensional small aromatics (counts) obtained for: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Fig. 6
Two-dimensional small aromatics (counts) obtained for: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Close modal

In Fig. 6, an iso-contour of approximately 15% of the peak LIF count is shown in white on the 2D plots for the three Φjet to distinguish the spatial extent of the small aromatics for the two fuels. Figures 6(a) and 6(b) show that the spatial distribution of the small aromatics is similar for both the fuels under nonpremixed condition; however, the peak small aromatics LIF count is approximately 35% higher in the C16 fuel as compared to the C7 fuel. The first location of the detection of small aromatic LIF signal is similar for both the fuels; the small aromatics LIF signal is visible in the flame annular region at approximately 3 mm HAB and at 6 mm HAB on the flame centerline. The LIF signal from small aromatics initially increases as the one and two ring aromatic species formation processes are dominant. Further downstream, the signal peaks and then decreases as the small aromatics consumption processes through oxidation or formation of larger aromatics are dominant.

Figures 6(c) and 6(d) shows that for the Φjet = 24 condition, the spatial distribution of small aromatics is similar for the two fuels; the small aromatics LIF signal is first detected in the annular region at lower HAB at approximately 4 mm HAB as compared to 7 mm HAB on the flame centerline. Similar to the Φjet = Inf flame, the peak small aromatics LIF count is higher in the C16 flame as compared to the C7 flame by approximately 30%. The small aromatics spatial distribution for the two fuels at the Φjet = 24 condition is similar to the Φjet = Inf condition. However, for the two fuels, the peak count of small aromatics decreases at Φjet = 24 as compared to the corresponding Φjet = Inf flame.

With increased level of air addition to the fuel, at Φjet = 6 shown in Figs. 6(e) and 6(f), the spatial distributions and the peak level of small aromatics are similar in both the fuels. It could be possible that the LIF technique is unable to resolve the differences in the level of small aromatics in the two fuels at the Φjet = 6 condition because of low signal-to-noise ratio. In addition to the decrease in the peak LIF count of small aromatics in the Φjet = 6 flame as compared to the Φjet = 24 flame, substantial differences are visible in the spatial distribution of small aromatics in the Φjet = 6 as compared to the Φjet = Inf and 24 flames. First, in the Φjet = Inf and 24 flames, the peak small aromatics LIF count is higher in the annular region than on the centerline. However, in the Φjet = 6 flames, the peak LIF count of small aromatics is near the centerline at 20–25 mm HAB. Furthermore, with the increase in the premixing, the location where small aromatics LIF signal is first detected on the flame centerline shifts to higher HAB; this location is approximately 6 mm HAB for Φjet = Inf and 24 flames, and 10 mm HAB for Φjet = 6 flame. This shift in the location and reduction in small aromatics with more premixing could be due to a combination of dilution, chemical, thermal, residence time of small aromatics, and LIF signal quenching due to oxygen, as discussed in detail in Ref. [5].

The differences evident in the peak level of small aromatics of the two fuels can be due to differences in temperature in the flame and/or chemical effects of the fuel. The comparison of radial temperature measurement at 5 mm HAB, shown in Figs. 7(a)7(c), and the centerline temperature measurement until 20 mm HAB, shown in Figs. 7(d)7(f), for the two fuels and three Φjet indicate that within measurement uncertainty, the temperature is similar for the two fuels at a given Φjet. Therefore, in a similar temperature field, the difference in the peak level of small aromatics of the two fuels at a particular Φjet is likely due to the chemical effect of the fuel breakdown process and growth to the one and two ring aromatics. We used a detailed chemical mechanism developed as part of the program [3] to compare the pyrolysis of the pure n-dodecane and the n-heptane/n-dodecane fuel blend. In addition, we used the LLNL high temperature n-alkane mechanism [25] to compare the pyrolysis of the C16 and C7 fuels. Figure 8 shows the consumption of n-heptane, n-dodecane, and n-hexadecane in a pyrolysis simulation under homogeneous reactor conditions at a constant temperature and pressure of 1600 K and 1 atm. These simulations show that n-dodecane and n-hexadecane decompose more quickly than n-heptane. The faster decomposition of n-hexadecane results in higher concentrations of reaction products, including the C3/C4 species that form the first aromatic ring. For example, Fig. 9 shows the comparison of the mole fraction of C3H3 and C4H6 species for the C7 and C16 fuel mixtures in a pyrolysis (with no air in reactants) simulation under homogeneous reactor condition at constant temperature (1600 K) and pressure (1 atm) using CHEMKIN. These simulations show that the peak mole fraction of C3H3 and C4H6 is higher in C16 fuel mixture as compared to the C7 fuel mixture. The C3H3 radical can lead to benzene formation through a self-addition reaction [26]. The C4H6 can form i-C4H5 or n-C4H5, which can form benzene through pathways described in Ref. [26]. The homogeneous reactor simulation with addition of air to the fuel (Φjet = 24 and 6) under the same temperature and pressure conditions shows that the n-hexadecane decomposes faster than n-heptane even under premixed conditions.

Fig. 7
Comparison of temperature measurement for the two fuels. Top row: radial measurement at 5 mm HAB; (a) Φjet = Inf, (b) Φjet = 24, and (c) Φjet = 6. Bottom row: axial measurement until 20 mm HAB for (d) Φjet = Inf, (e) Φjet = 24, and (f) Φjet = 6.
Fig. 7
Comparison of temperature measurement for the two fuels. Top row: radial measurement at 5 mm HAB; (a) Φjet = Inf, (b) Φjet = 24, and (c) Φjet = 6. Bottom row: axial measurement until 20 mm HAB for (d) Φjet = Inf, (e) Φjet = 24, and (f) Φjet = 6.
Close modal
Fig. 8
Comparison of the consumption of n-dodecane, n-heptane, and n-hexadecane in a pyrolysis simulation under homogeneous reactor condition at constant temperature (T = 1600 K) and pressure (1 atm) using CHEMKIN
Fig. 8
Comparison of the consumption of n-dodecane, n-heptane, and n-hexadecane in a pyrolysis simulation under homogeneous reactor condition at constant temperature (T = 1600 K) and pressure (1 atm) using CHEMKIN
Close modal
Fig. 9
Comparison of C3H3 and C4H6 mole fraction for C7 and C16 fuel mixtures in a pyrolysis simulation under homogeneous reactor condition at constant temperature (T = 1600 K) and pressure (1 atm) using CHEMKIN
Fig. 9
Comparison of C3H3 and C4H6 mole fraction for C7 and C16 fuel mixtures in a pyrolysis simulation under homogeneous reactor condition at constant temperature (T = 1600 K) and pressure (1 atm) using CHEMKIN
Close modal

To complete the discussion of the LIF results, Fig. 10 shows the 2D spatial development of large aromatics for the two fuels at three Φjet. Figures 10(e) and 10(f) includes LII signal at 30–35 mm HAB, which is comparable to the LIF signal at 20–25 mm HAB. The LII signal appears because the laser beam intensity could not be lowered further than 585 μJ, as doing so would reduce the LIF signal-to-noise ratio in other flame conditions. As a result of the low signal levels, no conclusions can be drawn about the large aromatics in the Φjet = 6 flames.

Fig. 10
Two-dimensional large aromatics (counts) obtained for: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Fig. 10
Two-dimensional large aromatics (counts) obtained for: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Close modal

Figures 10(a) and 10(b) show that the spatial distribution of large aromatics is similar for C16 and C7 nonpremixed flames; the large aromatics LIF signal is first detected at approximately 6 mm HAB and 11 mm HAB in the annular region and on the flame centerline, respectively. Similar to the small aromatics shown in Figs. 6(a) and 6(b), the peak LIF count of large aromatics is higher in the flame annular region as compared to the centerline. Quantitatively, the peak large aromatic LIF count is similar, within measurement uncertainty, for both the fuel mixtures. Figures 10(c) and 10(d) show that the spatial distribution of large aromatics for Φjet = 24 flames is similar to the Φjet = Inf flames, shown in Figs. 10(a) and 10(b). Within experimental uncertainty, the peak LIF count of large aromatics is similar for both the fuels. Similar to small aromatics discussed earlier, the peak LIF count of large aromatics decrease as compared to the nonpremixed flame with the increase in air addition to fuel jet for both C16 and C7 fuels. As expected, for each flame condition, the HAB where large aromatics LIF signal is first detected is at higher HAB as compared to the location at which small aromatics LIF signal is first detected, because large aromatics form from the small aromatics. The spatial or temporal development of small aromatics to large aromatics is discussed in detail in Refs. [4] and [5].

Laminar Flame: Soot Volume Fraction.

Figure 11 shows the 2D spatial development of soot volume fraction for the two fuels at three Φjet. An iso-contour of approximately 10% of peak soot volume fraction is shown to illustrate the spatial extent of LII signal and compare the results at three Φjet. The horizontal stripes visible in the images are an artifact from variation in the local fluence of the laser sheet. Figure 11 shows that the LII is first detected in the flame annular region, and further downstream, the soot field develops toward the flame centerline for all six flames. The spatial development of soot in this co-flow flame is consistent with the previous published literature on nonpremixed and partially premixed flames [17].

Fig. 11
Two-dimensional soot volume fraction (ppm) obtained for: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Fig. 11
Two-dimensional soot volume fraction (ppm) obtained for: (a) C7Φjet = Inf, (b) C16Φjet = Inf, (c) C7Φjet = 24, (d) C16Φjet = 24, (e) C7Φjet = 6, and (f) C16Φjet = 6
Close modal

With the increase in premixing level, the spatial extent of the LII signal reduces for both the fuels; the LII signal is detected until approximately 58 mm, 54 mm, and 45 mm HAB for the Φjet = Inf, 24, and 6 flames, respectively. Furthermore, the location at which the LII signal is first detected in the annular region shifts further downstream with increased premixing for both the fuels; the LII signal is first detected at approximately 12 mm, 15 mm, and 20 mm HAB in the flame annular region for the Φjet = Inf, 24, and 6 flames, respectively. Quantitatively, the peak soot volume fraction is slightly higher in the C16 flame as compared to the C7 flame for both nonpremixed and rich premixed flames. Even though the differences in the soot volume fraction among the two fuels at a particular Φjet is within the measurement uncertainty, the trend of slightly higher soot volume fraction in the C16 flames as compared to the corresponding C7 flames is repeatable. The trend of increase in soot volume fraction with increase in the carbon number in the n-alkane family is consistent with previous work [27,28]; the threshold sooting index, a global sooting parameter, showed a trend of increasing threshold sooting index with increase in the carbon number of the n-alkane fuels.

The peak soot volume fraction is lower for the two fuels at the Φjet = 24 and 6 conditions as compared to the nonpremixed flame. The percentage reduction in peak soot volume fraction is similar for both the fuels; the peak soot volume fraction reduces by approximately 18% and 50% for the Φjet = 24 and 6 flames, as compared to the nonpremixed flames.

In order to understand the spatial location of peak soot volume fraction, we plot the soot volume fraction on a streamline that passes through the flame center and maximum soot volume fraction along the downstream direction. Figures 12(a)12(c) show the comparison of the centerline soot volume fraction for the two fuels at Φjet = Inf, 24, and 6 flames, respectively. For all the three Φjet, the soot volume fraction increases initially due to dominant soot growth processes, reaches a peak, and then decreases at higher HAB as oxidation processes are dominant.

Fig. 12
Centerline comparison for C16 and C7 fuels: Soot volume fraction for (a) Φjet = Inf, (b) Φjet = 24, and (c) Φjet = 6. Maximum soot volume fraction for (d) Φjet = Inf, (e) Φjet = 24, and (f) Φjet = 6. Note: plots are 1 mm moving average.
Fig. 12
Centerline comparison for C16 and C7 fuels: Soot volume fraction for (a) Φjet = Inf, (b) Φjet = 24, and (c) Φjet = 6. Maximum soot volume fraction for (d) Φjet = Inf, (e) Φjet = 24, and (f) Φjet = 6. Note: plots are 1 mm moving average.
Close modal

Figures 12(d)12(f) show the maximum soot volume fraction for the Φjet = Inf, 24, and 6 flames, respectively, for the two fuels. At a particular HAB, higher value of maximum soot volume fraction as compared to the centerline soot volume fraction, indicates that the peak soot volume location is in the flame annular region. The C16 and C7Φjet = Inf flames have a peak soot volume fraction of approximately 5 ppm, in both the flame annular region and on the centerline. Similar to the nonpremixed flames, the C16 and C7Φjet = 24 flames have a peak soot volume fraction of approximately 4 ppm, in both the flame annular region and on the centerline. In contrast to the Φjet = Inf and 24 flames, Figs. 12(c) and 12(f) show that the Φjet = 6 flames have peak soot volume fraction of approximately 2.7 ppm on the flame centerline for both the fuels.

The LII measurements in the jet flames using prevaporized fuels shows that the C16 fuel result in slightly higher peak soot volume fractions as compared to the C7 fuel under both nonpremixed and rich premixed conditions. With a nearly identical temperature field and soot residence time, the differences in the soot volume fraction between the two fuels are most likely due to the chemical effect of the fuel. In the laminar flame experiment, both fuels were prevaporized, whereas in a model combustor, both the differences in chemical structure and fuel volatility may impact the combustor soot emissions. In order to understand the effect of the fuel volatility of the two fuels: n-heptane and n-hexadecane, which varies by 189 °C, on soot emission, the fuel mixtures are tested in a model combustor.

Soot Emissions From the Model Gas Turbine Combustor.

Figure 13 shows the comparison of line-of-sight soot volume fraction results as a function of Φglobal for the two binary mixtures and JP-8 fuel. Φglobal is defined as the ratio of inlet fuel-to-air ratio. The plot shows result from two different experiments for each of the three fuels. Each experiment is comprised of three or four sweeps in Φglobal from 0.9 to 1.8 and then again from 1.8 to 0.9 in increments of 0.1. The uncertainty in the equivalence ratio is estimated to be±2.6% [21]. The line-of-sight soot volume fraction at the exit of the combustor provides an estimate of the net soot produced in the combustor.

Fig. 13
Line-of-sight soot volume fraction results for JP-8, C7 fuel and C16 fuel as a function of Φglobal
Fig. 13
Line-of-sight soot volume fraction results for JP-8, C7 fuel and C16 fuel as a function of Φglobal
Close modal

The soot emission increases with the increase in the Φglobal for all the three fuels; this result is reflected in the laminar flame studies, where higher Φjet results in higher peak soot volume fraction. Additionally, the increase in Φglobal would result in lower adiabatic flame temperature, thereby leading to slower droplet evaporation. The slower evaporation rate may result in relatively bigger size droplets, at Φglobal = 1.8 as compared to at 1.0, to convect to the flame front, thereby increasing soot formation due to fuel rich regions. The comparison of soot volume fraction among the three fuels shows that the net soot produced by the two binary mixtures is substantially less as compared to the JP-8 fuel in the range of Φglobal from 1.0 to 1.8. At Φglobal = 1.8, the soot emission from JP-8 is approximately 6 ppm, nearly a factor of six higher than the soot emission from the two fuel mixtures, which emit approximately 1 ppm of soot. Furthermore, for the two n-alkane fuels, the Φglobal at which measurable amount of soot using the LE technique is first detected is at a richer Φglobal of approximately 1.4 as compared to at approximately Φglobal = 1.2 for the JP-8 fuel. Additionally, the factor by which the soot volume fraction is higher in JP-8 fuel as compared to the two fuel mixtures, increases with the increase in Φglobal from 1.4 to 1.8. Within experimental uncertainty, the C16 and C7 fuels have a similar amount of net soot production for Φglobal of 1.0–1.8.

The JP-8 fuel emits more soot as compared to the two binary fuel mixtures, because JP-8 contains a range of aromatic species, 20% by mass [29]; the aromatic components in a fuel can result in a higher soot production as compared to a fuel that has no aromatic content. The higher soot production in a fuel containing aromatic species is consistent with our recent study done in jet flames [4,5]. In the nonpremixed flames, replacing 25% volume fraction of n-dodecane with m-xylene, thereby maintaining constant carbon flow rate, resulted in an increase in the peak soot volume fraction by approximately three times the amount produced in the n-dodecane flame. Furthermore, the m-xylene/n-dodecane flame emits smoke since the soot produced is not completely oxidized. This increase in soot production is because the m-xylene reacts to form radicals such as m-xylyl that react to form two-ring PAH, thereby accelerating the formation of larger aromatic species that are involved in the formation of incipient soot particles. The faster production of PAH also leads to higher concentrations of these species, which increases the amount of soot mass through condensation processes.

The hypothesis of this study was that fuels with different boiling points would affect the fuel droplet vaporization processes, and that a fuel with a higher boiling point would produce higher levels of soot. However, Fig. 13 shows that soot emission from both the C16 and C7 fuels is quite similar in the gas turbine combustor. It is possible that the current design of the injector and combustor results in the complete vaporization of the fuel droplets before the combustion zone. If the fuel droplets are vaporized before the combustion zone, the impact of fuel boiling point on soot emission may not be evident in the current experiment, just as it was not in the prevaporized laminar flame. In future, the fuel–air mixing and droplet size measurements need to be done to better interpret the results. The chemical effect of the fuel molecular weight on soot is expected to remain. However, the experiment in jet flames show that the chemistry has relatively small effect on the peak soot volume fraction, which cannot be resolved by the LE technique applied on the model gas turbine combustor.

Conclusions

The effect of fuel volatility on soot emissions is tested in a model gas turbine combustor. The change in the fuel volatility did not affect the soot emissions in the model combustor operating at Φglobal = 1.0–1.8. The absence of the impact of fuel volatility on soot emission could be likely due to unique design of the injector and combustor, and hence, it is recommended to test the fuel matrix in other types of combustors. The JP-8 fuel has significantly higher soot emission as compared to the two fuel mixtures. In addition, the LE technique first detects soot emission approximately at richer Φglobal = 1.4 for the n-alkane fuel mixtures as compared to at Φglobal = 1.2 for the JP-8 fuel.

The two binary fuel mixtures are also tested in co-flow flame under nonpremixed and rich premixed conditions to simulate the fuel rich conditions present in the primary zone of the combustor. This experiment aids in understanding the chemical impact of the two fuel mixtures on aromatics and soot volume fraction. The LII measurement shows that the C16 fuel produces slightly higher peak soot volume fraction as compared to the C7 fuel at all three Φjet; however, the difference in soot volume fraction between the two fuels is within the measurement uncertainty of ± 15%. Additionally, the C16 produces approximately 35% higher peak level of small aromatics as compared to the C7 fuel at Φjet = Inf and 24. For both nonpremixed and premixed flames, the peak large aromatics level is similar for the two fuel mixtures, within the measurement uncertainty. The addition of air results in the decrease in the peak level of small aromatics, large aromatics, and soot volume fraction in Φjet = 24 and 6, as compared to in Φjet = Inf, for both the fuels. The LII-LIF results presented here provide a unique dataset to test chemical mechanisms and soot models for the n-alkane fuels. Furthermore, the results from model combustor can be used to validate numerical models in an engine-like environment.

Acknowledgment

We are grateful to Dr. Tim Edwards, Dr. Scott Stouffer, and Dr. Mel Roquemore for providing fuels used in the study. The authors would also like to acknowledge Dr. Venkatesh Iyer for his help with the experiment.

Funding Data

  • Strategic Environmental Research and Development Program (Grant No. WP-2145).

Nomenclature

C7 =

n-heptane/n-dodecane mixture

C16 =

n-hexadecane/n-dodecane mixture

HAB =

height above the burner

Inf =

nonpremixed flame

LE =

laser extinction

LIF =

laser induced fluorescence

LII =

laser induced incandescence

PAH =

polycyclic aromatic hydrocarbons

2D =

two-dimensional

Φjet =

jet equivalence ratio in jet flames

Φglobal =

global equivalence ratio in model gas turbine combustor

References

1.
ASTM,
2013
, “
Standard Specification for Aviation Turbine Fuel Containing Synthesized Hydrocarbons
,” ASTM International, West Conshohocken, PA, Standard No.
D7566-17
https://www.astm.org/Standards/D7566.htm.
2.
Edwards
,
T.
,
Moses
,
C.
, and
Dryer
,
F.
,
2010
, “
Evaluation of Combustion Performance of Alternative Aviation Fuels
,”
AIAA
Paper No. 2010-7155.
3.
Roquemore
,
W.
, and
Litzinger
,
T.
,
2016
, “
The Science of Emissions From Alternative Fuels
,” Strategic Environmental Research and Development Program/Environmental Security Technology Certification Program, Alexandria, VA, Report No.
WP-2145
.https://www.serdp-estcp.org/Program-Areas/Weapons-Systems-and-Platforms/Noise-and-Emissions/Air-Emissions/WP-2145
4.
Wang
,
Y.
,
Makwana
,
A.
,
Iyer
,
S.
,
Linevsky
,
M.
,
Santoro
,
R. J.
,
Litzinger
,
T. A.
, and O'Connor, J., 2017, “
Effect of Fuel Composition on Soot and Aromatic Species Distributions in Laminar, Co-Flow Flames—Part 1: Non-Premixed Fuel
,”
Combust. Flame
, epub.
5.
Makwana
,
A.
,
Wang
,
Y.
,
Iyer
,
S.
,
Linevsky
,
M.
,
Santoro
,
R. J.
,
Litzinger
,
T. A.
, and
O'Connor
,
J.
, 2017, “
Effect of Fuel Composition on Soot and Aromatic Species Distributions in Laminar, Co-Flow Flames—Part 2: Premixed Fuel
,”
Combust. Flame
, epub.
6.
Lefebvre
,
A. H.
,
1998
,
Gas Turbine Combustion
,
CRC Press
, Boca Raton, FL.
7.
Mordaunt
,
C.
,
2005
, “
Dual Fuel Issues Related to Performance, Emissions, and Combustion Instability in Lean Premixed Gas Turbine Systems
,” Ph.D. thesis, Pennsylvania State University, University Park, PA.
8.
Nisbet
,
I. C.
, and
LaGoy
,
P. K.
,
1992
, “
Toxic Equivalency Factors (TEFs) for Polycyclic Aromatic Hydrocarbons (PAHs)
,”
Regul. Toxicol. Pharmacol.
,
16
(
3
), pp.
290
300
.
9.
Inal
,
F.
, and
Senkan
,
S. M.
,
2002
, “
Effects of Equivalence Ratio on Species and Soot Concentrations in Premixed n-Heptane Flames
,”
Combust. Flame
,
131
(
1
), pp.
16
28
.
10.
Pastor
,
J. V.
,
García-Oliver
,
J. M.
,
García
,
A.
,
Micó
,
C.
, and
Möller
,
S.
,
2016
, “
Application of Optical Diagnostics to the Quantification of Soot in n-Alkane Flames Under Diesel Conditions
,”
Combust. Flame
,
164
, pp.
212
223
.
11.
Zhang
,
T.
,
Zhao
,
L.
, and
Thomson
,
M. J.
,
2017
, “
Effects of n-Propylbenzene Addition to n-Dodecane on Soot Formation and Aggregate Structure in a Laminar Coflow Diffusion Flame
,”
Proc. Combust. Inst.
,
36
(
1
), pp.
1339
1347
.
12.
Colket
,
M.
,
Edwards
,
T.
,
Williams
,
S.
,
Cernansky
,
N. P.
,
Miller
,
D. L.
,
Egolfopoulos
,
F.
,
Lindstedt
,
P.
,
Seshadri
,
K.
,
Dryer
,
F. L.
, and
Law
,
C. K.
,
2007
, “
Development of an Experimental Database and Kinetic Models for Surrogate Jet Fuels
,”
AIAA
Paper No. 2007-770.
13.
Farrell
,
J.
,
Cernansky
,
N.
,
Dryer
,
F.
,
Law
,
C.
,
Friend
,
D.
,
Hergart
,
C.
,
McDavid
,
R.
,
Patel
,
A.
,
Mueller
,
C. J.
, and
Pitsch
,
H.
,
2007
, “
Development of an Experimental Database and Kinetic Models for Surrogate Diesel Fuels
,”
SAE
Paper No. 2007-01-0201.
14.
Eddings
,
E. G.
,
Yan
,
S.
,
Ciro
,
W.
, and
Sarofim
,
A. F.
,
2005
, “
Formulation of a Surrogate for the Simulation of Jet Fuel Pool Fires
,”
Combust. Sci. Technol.
,
177
(
4
), pp.
715
739
.
15.
Nakamura
,
H.
,
Suzuki
,
S.
,
Tezuka
,
T.
,
Hasegawa
,
S.
, and
Maruta
,
K.
,
2015
, “
Sooting Limits and PAH Formation of n-Hexadecane and 2,2,4,4,6,8,8-Heptamethylnonane in a Micro Flow Reactor With a Controlled Temperature Profile
,”
Proc. Combust. Inst.
,
35
(
3
), pp.
3397
3404
.
16.
Douce
,
F.
,
Djebaïli-Chaumeix
,
N.
,
Paillard
,
C.-E.
,
Clinard
,
C.
, and
Rouzaud
,
J.-N.
,
2000
, “
Soot Formation From Heavy Hydrocarbons Behind Reflected Shock Waves
,”
Proc. Combust. Inst.
,
28
(
2
), pp.
2523
2529
.
17.
Santoro
,
R.
,
Semerjian
,
H.
, and
Dobbins
,
R.
,
1983
, “
Soot Particle Measurements in Diffusion Flames
,”
Combust. Flame
,
51
, pp.
203
218
.
18.
Kohse-Höinghaus
,
K.
, and
Jeffries
,
J. B.
,
2002
,
Applied Combustion Diagnostics
, Taylor & Francis, New York.
19.
Mouis
,
A.
,
Menon
,
A.
,
Katta
,
V.
,
Litzinger
,
T.
,
Linevsky
,
M.
,
Santoro
,
R.
,
Zeppieri
,
S.
,
Colket
,
M.
, and
Roquemore
,
W.
,
2012
, “
Effects of m-Xylene on Aromatics and Soot in Laminar, N2-Diluted Ethylene Co-Flow Diffusion Flames From 1 to 5 atm
,”
Combust. Flame
,
159
(
10
), pp.
3168
3178
.
20.
Zizak
,
G.
,
Cignoli
,
F.
,
Montas
,
G.
,
Benecchi
,
S.
, and
Donde
,
R.
, 1996, “
Detection of Aromatic Hydrocarbons in the Exhaust Gases of a Gasoline IC Engine by Laser-Induced Fluorescence Technique
,”
Recent Res. Dev. Appl. Spectrosc.
,
1
, pp.
17
24
.https://www.researchgate.net/profile/Roberto_Donde/publication/284088644_Detection_of_aromatic_hydrocarbons_in_the_exhaust_gases_of_a_gasoline_IC_engine_by_laser-induced_fluorescence_technique/links/564e012608ae1ef9296c2802/Detection-of-aromatic-hydrocarbons-in-the-exhaust-gases-of-a-gasoline-IC-engine-by-laser-induced-fluorescence-technique.pdf
21.
Iyer
, V
. R.
,
Iyer
,
S. S.
,
Linevsky
,
M. J.
,
Litzinger
,
T. A.
,
Santoro
,
R. J.
,
Dooley
,
S.
,
Dryer
,
F. L.
, and
Mordaunt
,
C. J.
,
2014
, “
Simulating the Sooting Propensity of JP-8 With Surrogate Fuels From Hydrocarbon Fluids
,”
J. Propul. Power
,
30
(
5
), pp.
1410
1418
.
22.
Beér
,
J. M.
, and
Chigier
,
N. A.
,
1972
,
Combustion Aerodynamics
, Applied Science Publishers, Ltd., London.
23.
Santoro
,
R. J.
, and
Mordaunt
,
C.
,
2002
, “
Dual Fuel Issues Related to Performance, Emissions and Combustion Instability in Gas Turbine Systems
,” Pennsylvania State University, University Park, PA.
24.
Glassman
,
I.
,
Yetter
,
R. A.
, and
Glumac
,
N. G.
,
2014
,
Combustion
,
Academic Press
, New York.
25.
Sarathy
,
S. M.
,
Westbrook
,
C. K.
,
Mehl
,
M.
,
Pitz
,
W. J.
,
Togbe
,
C.
,
Dagaut
,
P.
,
Wang
,
H.
,
Oehlschlaeger
,
M. A.
,
Niemann
,
U.
, and
Seshadri
,
K.
,
2011
, “
Comprehensive Chemical Kinetic Modeling of the Oxidation of 2-Methylalkanes From C7 to C20
,”
Combust. Flame
,
158
(
12
), pp.
2338
2357
.
26.
McEnally
,
C. S.
,
Pfefferle
,
L. D.
,
Atakan
,
B.
, and
Kohse-Höinghaus
,
K.
,
2006
, “
Studies of Aromatic Hydrocarbon Formation Mechanisms in Flames: Progress Towards Closing the Fuel Gap
,”
Prog. Energy Combust. Sci.
,
32
(
3
), pp.
247
294
.
27.
Olson
,
D.
, and
Pickens
,
J.
,
1984
, “
The Effects of Molecular Structure on Soot Formation—I: Soot Thresholds in Premixed Flames
,”
Combust. Flame
,
57
(
2
), pp.
199
208
.
28.
Olson
,
D.
,
Pickens
,
J.
, and
Gill
,
R.
,
1985
, “
The Effects of Molecular Structure on Soot Formation—II: Diffusion Flames
,”
Combust. Flame
,
62
(
1
), pp.
43
60
.
29.
Dryer
,
F. L.
,
Jahangirian
,
S.
,
Dooley
,
S.
,
Won
,
S. H.
,
Heyne
,
J.
,
Iyer
, V
. R.
,
Litzinger
,
T. A.
, and
Santoro
,
R. J.
,
2014
, “
Emulating the Combustion Behavior of Real Jet Aviation Fuels by Surrogate Mixtures of Hydrocarbon Fluid Blends: Implications for Science and Engineering
,”
Energy Fuels
,
28
(
5
), pp.
3474
3485
.