1 Introduction

The automotive industry’s efforts for improving fuel efficiency to meet new standards and regulations have resulted in an evolution of the gasoline engine. By 2018, gasoline direct injection (GDI) engines were in 51% of all new light-duty automobiles sold in the USA, while the European Union (EU) researched this same level in 2017, surpassing light-duty diesel cars for the first time [1, 2]. It has been consistently shown that the switch from port fuel injection (PFI) to GDI increased both the number and mass of particulate matter (PM) [3,4,5]. Furthermore, the adoption of lean-GDI operation for improved fuel efficiency could lead to even greater PM emissions [6,7,8].

Whereas the switch to GDI improves fuel economy [9] compared with PFI, GDI can result in increased air-fuel stratification and fuel-spray impingement on combustion-chamber surfaces, potentially promoting PM formation [10,11,12]. Similarly, recent studies have reported on the impact of engine operating conditions [13,14,15] and fuel properties [16, 17] on GDI PM.

Simultaneous with the move away from PFI, which produced little to no PM, to GDI was the implementation of more stringent vehicle emission regulations for PM globally. The US Environmental Protection Agency (EPA) began phasing in Tier 3 PM emission standards with MY2017 [18]. The Tier 3 standard drops the allowable PM mass emission rate 70% from 10 to 3 mg/mile over the urban driving schedule of the federal test procedure (FTP), with 100% fleet phase-in to be completed by MY2022. The California Air Resources Board currently has a similar 3 mg/mile PM emission regulation as part of its low emission vehicle (LEV) III standard that also began phase-in with MY2017, but also calls for a 1-mg/mile standard with a proposed phase-in to start with MY2025 [19]. The EU implemented a particle number (> 23 nm) regulation of 6 × 1011 #/km (9.7 × 1011 #/mile) for GDI vehicles under Euro 6 in addition to the PM mass regulation of 4.5 mg/km (7.2 mg/mile) [2, 20]. China introduced the same particle number and PM mass regulation as the Euro 6 in their CN5 regulations, but required a drop to 3 mg/km (4.8 mg/mile) in their C6b targets with 2023 implementation [20]. As these mass-based PM emission regulations continue to become more stringent in the USA, and with the introduction of particle number regulations in the EU and China, there is a critical need for more research into the properties of GDI PM.

While diesel studies often refer to PM as soot due to the high fraction of elemental carbon type particulate, PM can also include organic carbon and inorganic ash particles. The compositional distribution of these types of particulate in PM can be influenced by changes in combustion strategies and fuel properties [21, 22]. For accuracy, this paper will use the more inclusive PM term, unless discussing a specific component of the PM emissions.

A wall-flow GPF, similar to a diesel particulate filter (DPF), is already a part of aftertreatment strategies to control particle number and PM emissions in the EU and other countries and may become necessary in the USA to meet future regulations. The widespread implementation of DPFs in diesel vehicles has led to a better understanding of the filtration mechanism. Partially loaded DPFs have greater PM removal efficiency than clean DPFs due to layers of previously trapped PM that act as metastable filtration layers on the DPF for subsequently generated particulate. However, in GPFs, the higher temperatures of gasoline exhaust with fuel cut events result in passive regeneration during normal vehicle operation, preventing a PM layer buildup [23]. It has also been shown [24, 25] that varying the location and loading of a washcoat on a catalyzed DPF can alter the porous structure of the particulate filter, further improving filtration efficiency. Recent studies on catalyzed GPFs have investigated how GPF ash accumulation could promote a similar type of increase in filtration efficiency [26, 27]. Ash accumulation has been reported to provide significant improvement to the filtration efficiency of a GPF over the vehicle lifetime [28]. Non-ash PM can also be expected to accumulate in vehicles which are not routinely driven at the higher speeds required for passive regeneration [29]. When PM accumulation does occur in GPFs, it can be located both on and in the porous filter walls. The potential for short-term buildup of PM on the GPF and/or the PM’s interaction with catalytic active sites, in cases where the emission control catalyst and particulate filter are combined, have created an increased interest in the oxidative reactivity of PM, which would also be relevant when consumer driving behavior requires active regeneration. An understanding of the oxidative reactivity is critical to developing both system models of the powertrain–GPF interactions and operational controls for managing PM accumulation in the GPF [30].

Reactivity measurements are commonly taken from thermal gravimetric analysis (TGA) of PM samples collected on quartz fiber filters or loose particulate knocked off wall-flow particulate filters after filter loading on an engine or a vehicle. The activation energies are then calculated using the Arrhenius equation on small segment extractions from the TGA results and therefore are non-isothermal. TGA, within a given study, can highlight the effects that a set of conditions has on the activation energy of PM samples; however, variances between experimental procedures used in different reported studies can make direct quantitative comparisons less accurate. The impact of different TGA experimental procedures, such as sample mass, temperature window integration, and ramp rate, has been noted in the literature, and activation energy calculations for PM samples from the same study vary widely, ranging from 67 to 256 kJ/mol [12, 30,31,32]. Wang et al., in a non-isothermal kinetic study, showed that ramp rate and initial PM mass used in the study impacted the calculated activation energy for the same PM [31]. Luo et al. reported the PM activation energy increased relative to the engine speed–load condition at which it was generated for series of fuels and that increasing ethanol blending reduced the activation energy at all engine conditions tested [12]. Bogara et al. reported results from an isothermal kinetic study using separate PM samples to measure oxidation rate at different temperatures, ranging from 400 to 750 °C, which gave an activation energy of 132 kJ/mol for E0 PM. In addition, they found a range of activations energies (~ 140–160 kJ/mol) for the same PM in non-isothermal studies at different ramp rates and sample masses, leading them to suggest that due to the impacts of ramp rate and initial mass of the sample on activation energy, these values should be used to understand trends but absolute values should be taken with caution [30]. Furthermore, since most TGA studies are on loose PM or PM on a microfiber filter, in contrast to an actual GPF which has a wall-flow, monolithic structure, the typical TGA approach must assume a uniform PM sample with a single activation energy regardless of whether the study uses a non-isothermal ramping method on a single sample or isothermal oxidations on multiple samples.

To better understand the oxidative reactivity of GDI PM in a GPF-based aftertreatment control system, the PM in the present study was loaded on small uncatalyzed GPFs during engine-dynamometer operation. To accumulate enough PM to study its oxidation kinetics, it was necessary to achieve a complete PM layer in the GPF. Although this is not fully representative of what is expected in the field, it is the only reasonable way to accumulate enough PM to study its oxidation kinetics in the mini-GPF cores required for the flow reactor setup. Additionally, the impact of oxygenated fuel blends on oxidative reactivity was investigated by comparing GPFs loaded with PM generated from gasoline with and without ethanol (E0 and E30). The GPF-loaded PM layers were then studied throughout exposure to successive oxidative burnouts. Studying PM reactivity sequentially on a GPF in a simulated exhaust gas flow reactor allowed the exposed surface area and PM layer morphology to factor into the activity measurements unlike typical TGA studies. Another benefit of having a PM layer in the GPFs is that it enables a better comparison for how the GDI-generated PM behaves contrasted to DPF PM layers reported in the literature [33, 34]. Thickness measurements made during successive burnouts with neutron imaging revealed distinct differences between the two.

Oxygenated fuels have been shown to improve efficiency as well as provide pathways to increase the use of renewable fuels. The US Department of Energy has an ongoing vehicle technology initiative, Co-Optima, on fuels and engine optimization which has actively investigated both bio-derived fuel production pathways and their impact on engine operation. Only two of the most promising blendstocks found in this program for light-duty vehicles, ethanol and isobutanol, are currently allowed in market fuels [35]. The USA currently promotes the use ethanol fuel blending of 10% (E10) for on-road vehicles and allows 15% ethanol (E15), while isobutanol (iBu) is approved for blending up to 16% in market fuels. Increased alcohol blending levels are of interest due to its potential to enhance emission reduction and improve engine efficiency at lower blending levels, with E30 being one of the most promising candidates [36, 37]. To achieve the same fuel oxygen content as E30, a 48% iBu blend was also used in this study. The impact of isobutanol blended fuels on PM emission flux and particle size has been discussed in the literature [14, 38]. Prior work has shown that blending of ethanol with gasoline results in an increased reactivity of GDI PM [12, 14, 31]. The purpose of this study was to measure the fundamental oxidation kinetics for PM loaded on wall-flow GPFs in the exhaust of a GDI engine running on several different fuels containing varying levels of ethanol and isobutanol.

2 Experimental

2.1 Engine, Fuels, and GPF Core Loading

The GPF samples were collected according to a previously described process [16]. The process is briefly described here. Particulate matter was loaded on uncatalyzed, cordierite mini-GPF cores (2.5 cm outer diameter × 7.6 cm long; NGK 12 mil/300 cpsi) by placing them four at a time in a special holder that directed the exhaust from a 2.0-L GM GDI engine, downstream of the catalyst, to flow in parallel through all four filters. This engine was specifically developed for high-efficiency operation with ethanol fuel blends. During GPF loading, the engine was operated under acceleration type conditions of 2600 rpm and 8 bar brake mean effective pressure at air-to-fuel ratios slightly higher than the engine’s minimal lambda (λ) of 0.89. At stoichiometric engine operation, there is exactly enough oxygen to combust all of the fuel hydrocarbons (λ = 1), while at λ < 1, as was used in this study, there is an excess of fuel. Prior to PM loading, the cores were ramped in air to 800 °C in a tube furnace to bake off any unknown impurities in the filter. Cores were loaded (ca. 4 g of PM per liter of filter) with PM from one of the four different volumetric fuel blends: an 87 anti-knock index (AKI) gasoline (E0), a 30% ethanol blend with the 87 AKI gasoline (E30), a 48% isobutanol blend with the 87 AKI gasoline (iBu48), or a 24% isobutanol blend with the 87 AKI gasoline (iBu24). PM collection quantities were determined by the change in filter weight before and after PM loading. Mass measurements were collected after drying the cores in a ~ 120 °C laboratory oven for about 2 h and weighed directly after removal from the oven.

2.2 Flow Reactor Experiments

2.2.1 Oxidation Study: Temperature Programed Oxidation

All oxidation studies were carried out on an automated flow reactor similar to the one described by Pihl et al. [39]. The GPF was placed at the end of a Mini-Mite® tube furnace inside a quartz tube, and the feed gases were supplied upstream of the furnace and metered in using a manifold of mass flow controllers (MFCs). The water was supplied by a high-pressure liquid chromatography pump which introduced the water directly into a heated flow of nitrogen. Fast-switching valves downstream of the MFCs but upstream of the furnace were used to control pulsed flows when needed. An MKS MultiGas 2030HS Fourier transform infrared (FTIR) spectrometer downstream of the GPF was used to monitor gas composition. The CO2 was quantified with a custom low concentration calibration developed in collaboration with MKS; all other gases were quantified using standard calibrations provided by MKS. A high-pressure liquid nitrogen dewar and an ultra-high purity compressed oxygen cylinder were used to feed the MFCs.

The temperature dependence for oxidation of each type of PM produced was studied by temperature programed oxidation (TPO). The oxidative reactivities of the PM produced from the different fuel blends (E0, E30, iBu24, iBu48) were measured with the PM still loaded on the GPF mini-cores. The TPO experiments were run at a gas hourly space velocity (GHSV) of 20,000 h−1 in the presence of 1% O2 and 10% H2O with a N2 balance at a 10 °C/min temperature ramp from 150 to 650 °C (E0 and E30) or 150 to 700 °C (iBu24 and iBu48). For calculation purposes, all PM mass on the GPFs were assumed to be carbon, and its removal was calculated from a mass balance based on the moles of CO and CO2 measured by the FTIR downstream of the flow reactor furnace.

2.2.2 Oxidation Study: Kinetics

Kinetic studies were only performed on PM from the E0 and E30-loaded GPF cores. The same automated flow reactor described in Section 2.1 was also used for these experiments. The kinetic studies were conducted in a stepwise fashion, summarized in Table 1, on a single core sample and encompassed four different types of oxidation treatments run in the following order:

  1. (1)

    Burnout oxidation: Step 1, to remove the required moles of carbon to advance the PM layer to the next C/C0 level to be studied (C0 is the initial moles of PM loaded on the filter, and C is the cumulative moles of PM removed by oxidation);

  2. (2)

    Pulsed oxidation–temperature sweep: Steps 2–6, used for activation energy calculation;

  3. (3)

    Pulsed oxidation–O2sweep: Steps 6–10, to probe the O2 dependence of the oxidation rate;

  4. (4)

    Temperature programed desorption (TPD) oxidation, Steps 11–12, for active surface area calculations.

Table 1 Experimental steps involved in the 4 types of oxidation treatment used to study each burnout level of a GPF

The four-part kinetic study was repeated, unless otherwise noted, at five different levels of GPF burnout, nominally, C/C0 = 0, 0.2, 0.4, 0.6, and 0.8 PM burnout, where C/C0 = 0 refers to the first tests on the fully loaded filter. The term burnout will be used to refer to the total cumulative fraction of the initial PM that has been removed by any previous oxidation steps or treatments (C/C0). The conditions for these studies are summarized in Table 1 and described in the preceding sections. Multiple GPF cores loaded with PM from either E0 or E30 were used for repeating the kinetic studies.

Burnout Oxidation

Each GPF was prepared for the next burnout level to be studied by oxidizing the remaining moles of PM on the GPF to reach the next C/C0 level. This step (step 1 of Table 1) is referred to as the burnout oxidation treatment (Fig. 1). The CO and CO2 evolved from the PM oxidation were integrated in real time by the simulated exhaust flow reactor control software, and feedback control was used to shut off the O2 flow when the total moles of carbon burned reached the target threshold.

Fig. 1
figure 1

Experimental conditions used for the burnout oxidation treatment. “–” target temperature profile, “----” %O2, and “••” %H2O in the feed gas

Pulsed Oxidation–Temperature Sweep

The pulsed oxidation–temperature sweep studies, steps 2–6 of Table 1, were carried out at each burnout level. The pulse conditions and temperatures were chosen to achieve reproducible oxidation rates without consuming significant quantities (< 10%) of the PM. The conditions were chosen based on the TPO results of E30 to reduce the potential for excessive burnout at a single temperature but were also used for E0 so that a direct comparison could be made.

Similar to a pulsed-oxidation study done on DPF mini-cores [39], the reactivity was averaged over the second half (last 15 s) of the last 3 pulses at each temperature in an attempt to capture stable activity. Reactivity was measured by this method at all 5 temperatures for each target burnout levels for 2–3 loaded GPFs and averaged for each temperature-C/C0 condition studied. The PM activation energies at each burnout level were calculated from Arrhenius plots of the average isothermal reactivity measurements.

Pulsed Oxidation–O2 Sweep

A second set of pulsed oxidations, pulsed oxidationO2sweep, was performed at 440 °C according to conditions listed for steps 6–10 of Table 1. The resulting reaction rates were used to determine the O2 dependence of the PM oxidation kinetics for each burnout level.

Temperature Programed Desorption Oxidation

Following the pulsed oxidation–O2sweep, a TPD oxidation experiment was run to estimate the active specific surface area [39] of the PM. The TPD oxidation protocol is shown in Fig. 2 and represented in steps 11–12 of Table 1. Integration of the CO and CO2 evolved during the oxygen-free portion of the TPD experiment was used to calculate the PM oxygen uptake that occurred during the 200 °C O2 exposure. The PM O2 uptake could then be used to calculate the active specific surface area per mole of PM remaining (Sm; m2 mol−1) at each C/C0 level; calculations assumed a cross-sectional area of 6800 m2 mol−1 [39] for the adsorbed atomic oxygen.

Fig. 2
figure 2

Experimental conditions used for TPD oxidation treatment. “–” temperature profile, “----” %O2, and “••” %H2O in the feed gas

2.3 Neutron Imaging of GPFs

In a previous study on successive regeneration of DPFs [33], neutron imaging was demonstrated to be a powerful, nondestructive technique capable of detecting thin layers of PM and ash on the filter walls with high spatial resolution throughout the volume of the filter. In the present study, a similar neutron-imaging procedure was conducted on E0 and E30 GPF cores loaded with PM like those used for the oxidation studies. Each GPF was first imaged with ~ 4 gPM/LGPF loading (C/C0 = 0). The cores were then individually placed in the simulated exhaust reactor, and a partial burnout of the PM was performed to take the core to the next PM burnout level. Both the E0 and E30 GPF cores were imaged at 6 nominal burnout levels (C/C0 = 0, 0.2, 0.4, 0.6, 0.8, and 1.0). The burnout procedure was the same as that used in the kinetic study and can be found as step 1 in Table 1.

2.3.1 Neutron Tomography

The High Flux Isotope Reactor at Oak Ridge National Laboratory generates a neutron flux of 2.5 × 1015 n cm−2 s−1 at 85 MW thermal power, making it one of the most powerful neutron sources for scientific research in the world. Some of these neutrons are diverted into over a dozen thermal or cold beamlines, and the Cold Guide 1D beamline (CG-1D) is devoted to imaging. The flux at the CG-1D aperture is typically ~ 106–107 n cm−2 s−1, and the neutron beam diverges at 2.7° traveling through a 450-cm long He-filled flight tube up to 565 cm to the imaging detection plane; by varying sample locations and/or aperture sizes, L/D (post-aperture flight path to aperture diameter, for a pinhole geometry) ratios ~ 400–800 may be obtained, resulting in different levels of resolution. Imaging is achieved by measuring the relative interaction, via absorption and scattering, of materials and neutrons between the aperture and the detector [34, 40].

At the detection plane is a 50 or 100-μm thick 6LiF/ZnS scintillator plate, orthogonal to the neutron beam, which emits photons from neutron collisions. These photons are measured with a 16-bit 2048 × 2048 pixel DW936 IkonL ANDOR™ CCD camera with a 7.4 × 7.4 cm field of view on the scintillator, for a pixel size of ≈ 0.0040 cm px−1. The resolution of this system is typically 0.0070–0.0200 cm depending on a variety of factors. The camera can transmit images at 1 frame/s, but typically 15–60 s of integrated exposure is required for sufficient contrast. Sample images are normalized between the open-beam image, which contains inhomogeneities in the beam caused by the neutron path guide structures, and a dark-field image, which represents internal camera current with no exciting photons.

The tomographic imaging procedure was conducted as in our previous study [33] on progressive regeneration in diesel particulate filters. Both GPFs were placed side by side in a special holder placed upon a rotational stage with 0.05° rotational angle precision. The GPF flow channels were oriented vertically, or perpendicular to the stage table surface. The tomographic radiographs were obtained by rotating the stage through a 183° range with 0.1° steps, and at every step, a radiograph was obtained, for a total of 1831 radiographs. The radiographs were normalized, as described above, and filtered using a median filter to remove small, localized spots in the images where incident gamma radiation saturated the detector before the tomographic reconstruction.

2.3.2 Tomographic Reconstruction

An iterative reconstruction algorithm was employed, a method which generally yields greater image contrast than an analytic algorithm such as filtered back projection. Image reconstruction was modeled as a regularized weighted least squares problem, which was solved using a SIRT-like (Simultaneous Iterative Reconstruction Technique) preconditioned gradient descent (PGD) algorithm, more specific details of which have been previously reported [41, 42].

Some computational techniques were implemented for computational and data memory efficiency. First, ordered subsets were used to accelerate the PGD algorithm, a multi-threaded approach that took advantage of the multi-core architecture of modern-day computers. The number of subsets was set to one-fourth of the number of projections during the first several iterations and was then gradually tapered off to one for the last few iterations. This approach accelerated the computation as much as possible early on and prevented ending up in a limit cycle toward the end. Second, to reduce the data memory requirements, the voxel size in the axial direction (along the GPF channel) was set to ≈ 0.03 cm, retaining the 0.004 cm voxel width across the channels to measure PM thickness.

2.3.3 Data Analysis

The data analysis procedure was conducted as in our previous study on progressive regeneration in diesel particulate filters [33]. Because the contrast between PM and the filter walls was too low to discriminate reliably and the PM layer was usually thinner than the detector resolution, the methodology for determining PM thickness outside of the filter walls relied on the differences in flow area between the outlet and the inlet channels. These differences are apparent in Fig. 3, which shows neutron tomographic sections, or slices, at the axial midpoint for the fully loaded and fully regenerated E0 GPF. In the fully loaded slice, there is a visibly apparent bias, with alternating smaller-area channels (the inlet channels with a PM layer) and larger-area channels (the outlet channels without PM); this pattern is not apparent in the fully regenerated slice. The following methodology was used to quantify these differences and estimate PM layer thickness: A set intensity threshold was used to differentiate between the open channel and the wall, and the channel pixel area was the set of contiguous pixels with sub-threshold intensities; both inlet and outlet channel areas were calculated in this manner, and the areas were computed for each channel of interest for each tomographic slice along the GPF axis. A schematic of the area computation process for a given inlet channel is shown in Fig. 4, in which the pixel area Ai for the inlet channel and pixel areas Ao,n for n (up to 4) neighboring outlet channels are calculated using image analysis; the Ao,n was averaged to obtain a single Ao. On average, the PM layer thickness, e, will be:

$$ e=\frac{b-a}{2}=\frac{\sqrt{A_{\mathrm{o}}}-\sqrt{A_{\mathrm{i}}}}{2}. $$
Fig. 3
figure 3

Tomographic slices across the axial center plane for the fully loaded (left) and fully regenerated (right) GPFs. The bright arcs on the filter outer surface are hydrogenous adhesives in tape-holding aluminum foil wrapping in place

Fig. 4
figure 4

Schematic of open-channel areas measured with image analysis and geometric factors used in determining PM layer thickness

As a result, the PM layer thickness for a given inlet channel was calculated as a function of axial distance, for all 6 regeneration stages. An ensemble of axial thickness functions from several dozen inlet channels was computed at each regeneration step, and the axially resolved mean thickness profile was computed as a final result.

The measured dispersion within the ensemble of PM layer thickness profiles was pronounced, but the mean profiles for each regeneration step followed a clear trend as more PM was removed. The C/C0 = 1.0 case, which removed all of the PM, was used to determine if there was an intrinsic difference in dimensions between the inlet and outlet channels. Using the method of areas as described above, the measured PM layer thickness for C/C0 = 1.0 was ± 0.0010 cm, placing it well below the measured thickness at C/C0 = 0.8; see below for demonstration.

3 Results and Discussion

3.1 Oxidation Study: TPO

The results of TPO experiments on GPFs for each of the four fuel blends are shown in Fig. 5. The TPO indicates that the PM oxidative reactivities for the four fuel blends were impacted by both type and quantity of the alcohol blended. It was previously reported [16] that the chemical composition, particle size, and quantity of PM generated were also impacted by the type of alcohol in the fuel blend used to operate the engine. The temperatures at which 50% (T50) and 90% (T90) of the PM were oxidized during the TPO experiments (Table 2) highlight the higher oxidation temperatures required for PM from isobutanol fuel blends and dramatically lower temperatures needed for an ethanol blend compared with neat gasoline (no alcohol, E0). The oxidative reactivity of PM based on these TPO experiments follows E30 > E0 > iBu24 ≈ iBu48. The TPO profile for E30 is shifted to lower temperatures and has a broader, more symmetric profile than the other fuels. Less active low-temperature oxidation (≤ 550 °C) is seen for E0, iBu24, and iB48 than for the E30 PM, resulting in asymmetric oxidation profiles and highlighting the difference ethanol can make on the PM properties. The greater low-temperature reactivity of E30 PM may be related to previously reported oxygen incorporation into the backbone of E30 PM, analytically detected as methyl propenoate [16]. The addition of ethanol to the gasoline dropped the T90 of the resulting PM by 65 °C. In contrast, the addition of isobutanol increased the T90 of the corresponding PM more than 100 °C compared with the neat gasoline PM, for both concentrations of isobutanol fuel blends studied.

Fig. 5
figure 5

PM oxidative reactivity from TPO experiments of 4 different gasoline blended fuels: E0 (100% gasoline): E30 (30% ethanol blend): iBu24 (24% isobutanol blend): iBu48 (48% isobutanol blend)

Table 2 Results from PM TPO study: temperatures at which 50% of PM is oxidized (T50) and 90% (T90)

3.2 Oxidation Study: Kinetics

3.2.1 Activation Energy by Pulsed Oxidation

As described in the Section 2, Table 1, isothermal, pulsed-oxidation studies were used to obtain isothermal oxidation rates in a kinetically controlled range, like those shown in Fig. 6. These rates were used in the Arrhenius equation to calculate activation energies for the PM collected with the E0 and E30 fuel blends. The activation energies, at each of the five different levels of PM burnout, for all GPF samples studied are plotted in Fig. 7. Figure 8 simplifies the results of Fig. 7 showing activation energy averages with min/max error bars for the five nominal burnout levels of the GPF PM studied (C/C0 = 0, 0.2, 0.4, 0.6, and 0.8). For comparison, Fig. 8 includes values reported by Pihl et al. [39] for a similar DPF PM oxidation study.

Fig. 6
figure 6

a Arrhenius plots for oxidation of PM on an E0 GPF with linear best fits used for activation energy calculations. Oxidation is plotted in terms of rate, r (mol/s), of PM oxidation divided by the total volume (L) of the GPF sample as a function of temperature, T (K). b The five, 30 s duration, O2 pulses with a 150-s delay between pulses for a E0 GPF sample at 430 °C with a C/C0 = 0.4. The blue boxes in b indicate the last 15 s of the last three pulses used to calculate an average oxidation rate which was averaged with data from two other E0 GPF samples to calculate the ▲ data point in a outlined in black

Fig. 7
figure 7

Activation energy (EA) plot of all data as a function of fraction of moles of carbon oxidized from loaded mini-GPF cores (C/C0 oxidation): EA calculated from pulsed-oxidation experiments performed after a burnout oxidation treatment (Δ); directly after another pulsed-oxidation experiment (○); and after a TPD oxidation (●). The average EA for the 5 targeted C/C0 oxidation levels (▲, dark line) and the corresponding min and max values from the averages (no markers, lighter lines) are also shown

Fig. 8
figure 8

Average activation energy (EA) as a function of carbon oxidized on the loaded particulate filter (C/C0) for E30 PM (orange triangles), E0 PM (blue diamonds), and reported diesel PM (black x). Error bars for both activation energy (vertical) and C/C0 (horizontal) represent min and max values for each average used. All averages for E0 and E30 were from 3 GPF cores, except the E0 C/C0 oxidation values at 18% and 74% which are of only 2 GPF cores. Trendline for E30 does not include nominal 0.8 C/C0 point

The activation energies of the PM for both fuel blends in this study (E0 and E30) at each burnout level, shown in Fig. 7 and Fig. 8, suggest that as PM is removed by oxidation, the chemical stability of the remaining PM increases, resulting in higher activation energies. The activation energies reported [39] for diesel particulate, still loaded on a particulate filter, indicate no change in the chemical stability (127 kJ/mol) of the PM using comparable burnout levels (~ 0.2 to 0.7 C/C0). The PM from both E0 and E30 GDI fuels, reported here, initially has lower activation energies (ca. 60–80 kJ/mol) than commonly reported for diesel PM, but by 60% burnout, these GDI PM activation energy values exceeded typical diesel activation energies. The PM generated from E30 remains chemically less stable than that from the E0 under oxidative conditions up to 70–75% burnout, at which point the activation energy of both types of PM was seen to converge at about 170 kJ/mol; see Fig. 8.

The oxidative chemical stability of the PM from 0 to 60% burnout increased linearly for both E0 and E30. As can be seen in the average activation energy plot, Fig. 8, the chemical stability for E30 increased more rapidly with successive PM burnouts than E0, making the PM harder to oxidize. While the rate of increase in the activation energy as a function of burnout continued for E0 up to 74% burnout, a sharp increase was observed for E30 from 60 to 73% burnout. At the final burnout level studied, E0 and E30 both converged at an activation energy near 170 kJ/mol. Prior to the PM activation energies of E0 and E30 converging, E30 PM had a lower overall chemical stability than the E0 PM, which agrees with the oxidative reactivity seen in the TPO experiments discussed earlier where more of the E30 PM burned out at a lower temperature than the E0 PM.

The activation energy measurement for one of the E0 PM samples at the nominally 0.8 burnout point dropped significantly compared with the other two E0 GPF samples (see orange diamonds in Fig. 7). However, since this drastically different data point could not be repeated within the scope of the project, it was not included in the average value at this nominal burnout level.

All of the trends discussed thus far were observed when the pulsed oxidation–temperature sweep, steps 2–6 in Table 1, was run following a step 1 burnout oxidation. Two different types of deviation from the steps 1–12 sequence were also performed in the course of the study and resulted in widely different impacts on the measured PM activation energy, as can be seen with the open and filled circle data points in Fig. 7. The open circles in Fig. 7 occurred when steps 2–6 were immediately repeated, resulting in the activation energy calculated from the second set of pulsed oxidations being greater than expected from the trends for both E0 and E30 samples. The second type of deviation resulted in lower than expected activation energies, filled circles in Fig. 7, and occurred when the E0 activation energies were calculated from steps 2–6 run after a TPD oxidation, steps 11–12, without running step 1 in between. This treatment was performed twice on the same E0 PM sample core, and in both cases, the calculated activation energy was the same, 61 kJ/mol, even though the C/C0 had changed and was the lowest calculated for E0 under any set of conditions but fell close to the average seen for the fully loaded E30 PM (C/C0 = 0).

Once it was observed that the sequence of the oxidation experiments influenced the measured activation energies, a consistent sequence of experiments was performed at each level of PM oxidation studied, such that pulsed oxidation (steps 2–6) was only performed immediately following a burnout oxidation (step 1) experiment. Thus, more studies would be required to determine if these impacts on activation energy would continue at higher burnout levels. The changes in activation energy observed as a result of the PM’s previous oxidizing environment suggest that multiple types of PM exist which have different barriers to oxidation.

Assuming the experiments were operating under the targeted isothermal, kinetically controlled conditions, the activation energy changes seen at successive C/C0 burnout levels may also indicate the presence of different types of PM, with different activation energies, in varying fractional combinations. Since the PM in each of the GPFs studied was generated at the same constant engine conditions, all of the PM generated from a given fuel would be expected to have a consistent chemical makeup. This chemical makeup would correspond to the fractional distribution of the different types of PM and associated hydrocarbons at C/C0 = 0. Earlier work [16] from our lab reported that some of the chemical properties, characterized from other PM samples from this study, varied between the E0 and E30 PM. PM samples characterized in the earlier report were collected on supported Teflon membrane filters during the same engine campaign as the GPFs discussed in this manuscript. Briefly, the chemical speciation of the PM samples was performed using a gas chromatography–mass spectrometer (GC–MS) to measure the hydrocarbon compounds produced during a thermal desorption step at 325 °C under helium and a pyrolysis step at 500 °C, collectively called a thermal desorption pyrolysis (TDP). All reported observations were made after normalizing for mass of PM analyzed and fraction of gasoline in the fuel blend. Note that the thermal desorption was done below the temperatures used in this GPF kinetic oxidation study while the pyrolysis step was at a higher temperature than the GPF samples were exposed to during steps 2–6, from which the activation energies were derived. The reported speciation of hydrocarbons released in the thermal desorption step indicated that the E0 PM had a greater fraction of polyaromatic hydrocarbons (PAHs) compared with E30. The pyrolysis step at 500 °C served to break down the particulate into subcomponents; those which elute earlier in the GC–MS measurement are typically considered to reveal PM backbone structural differences. Unique to E30 PM was the presence of methyl propenoate, an oxygenate, not seen in the E0 PM, suggesting that fuel oxygen had been incorporated into the backbone of the E30 PM. During the 500 °C pyrolysis step, E30 PM was shown to contain a higher concentration of multi-ring PAHs than E0 when normalized to the fraction of gasoline in the fuel blend. Considering this previous work from our lab, the types of PM could be categorized into (1) lower boiling point hydrocarbons, (2) higher boiling point hydrocarbons, (3) carbonaceous PM, often referred to as soot, and (4) oxygenated PM. Based on the different chemical properties of these different types of PM material, the changes in their fractional distribution could be expected to alter the stability and reactivity of the PM matrix as seen in the activation energies reported here.

The isothermal DPF study reported by Pihl et al. [39], which observed a constant activation energy of 127 kJ/mol for diesel PM on a DPF at all levels of PM burnout, followed an oxidation sequence similar to steps 1–12 that were repeated in order at each increasing burnout level. However, the Pihl et al. study used a burnout oxidation temperature of 600 °C (step 1) immediately prior to its pulsed oxidationtemperature sweep (steps 2–6) while the GPF study, reported herein, used 440 °C in step 1. When the activation energy was calculated for GPF PM samples immediately after a TPD oxidation (i.e., followed step order 11–12, 2–12, 2–6) which took the temperature to 650 °C for two consecutive C/C0 PM oxidation levels (filled circles, Fig. 7), the activation energy did not change, indicating that steps 11–12 exposed similar PM material. Interestingly, when the test sequence returned to the typical steps 1–12, 1–12, etc. order, the increasing activation energies were again observed.

3.2.2 Oxygen Dependence by Pulsed Oxidation

Using the reaction rates obtained from the pulsed oxidationO2sweep experiments (steps 6–10 in Table 1), the rate of PM oxidation was found to be 0.8 order in oxygen. This approximately first-order dependence on oxygen was consistent for oxidation at all C/C0 on both E0 and E30 PM.

3.2.3 PM Active Surface Area by TPD Oxidation

If the morphology of the PM layer impacted the activation energy measurement, the PM active surface area (SA) may contribute more to the reactivity than the fraction or quantity of PM on the GPFs. To assess this, the active SAs were calculated from production of CO and CO2 during the TPD oxidation experiments, steps 11–12, shown in Fig. 9. These active surface area measurements were taken following steps 2–10 at each C/C0 burnout level. The molar specific active surface area (Sm, m2/molPM) at each C/C0 was normalized to the GPF’s initially measured molar SA (Sm,0), C/C0 = 0, to allow better correlation across the multiple GPF samples studied for each fuel. The normalized data (Sm/Sm,0) are also plotted as a function of PM burnout, C/C0, in Fig. 9. The trend lines seen in Fig. 9 suggest that a linear relationship exists between the burnout level and the normalized data for both E0 and E30 PM. A slight deviation in the linear relationship for the E0 PM occurs when the TPD oxidation followed steps 2–10 run without a burnout oxidation (i.e., followed step order 11–12, 2–12). The Sm/Sm,0 observed in these cases (filled circles, Fig. 9) deviated from the trend found when calculated from runs that followed the repeated step 1–12 order. The molar active surface area also dropped lower than the trendline at the C/C0 = 0.8 burnout level for the E0 (gray band, Fig. 9). In both cases, the PM also had an activation energy which dropped below the trend, as discussed above (Fig. 7 and Fig. 8).

Fig. 9
figure 9

Normalized molar active surface area (Sm/Sm,0) as a function of the amount of carbon oxidized from the GPF (C/C0) is plotted in the primary y axis (filled symbols). The smaller, black, open symbols corresponding active surface areas without normalization (Sm) and are plotted against the secondary y axis. Data is from 2 GPFs loaded with E0 PM (top) and from 3 GPFs loaded with E30 PM (bottom). The E0 linear best fits do not include the filled circles which represent molar surface area values obtained from TPD oxidation when there was a deviation from the order described in Table 1 or the diamond located in the gray band

The greater active surface area per mole of PM remaining on the E30 GPF vs. E0 indicates that either more of the PM on E30 is physically accessible to the oxygen, and/or that the PM from the different fuels is physically or chemically different. The lower-temperature reactivity, activation energies, and compositional differences discussed earlier all suggest that the higher active surface area of E30 can be attributed to different types of PM and/or the fractional distribution of these types compared with E0.

3.3 Neutron Imaging of PM Layer During Oxidation

A morphological study of the PM layer thickness was performed using neutron imaging of PM-loaded mini-GPF cores as successive fractions of the PM were removed by burnout oxidation (C/C0 = 0, 0.2, 0.4, 0.6, 0.8, and 1.0) according to step 1 of Table 1. The two cores imaged were loaded from either E0 or E30 engine operation at the same time as those used in the TPO and kinetic studies reported above. The E0 and E30 cores were first imaged in the neutron beam fully loaded to C/C0 = 0. Between neutron-imaging experiments, 20% of the PM loading was removed by burnout, and then the cores were imaged using the neutron beam. The burnout and neutron imaging were repeated four more times for each core with the final oxidation (C/C0 = 1.0) removing all remaining PM such that the final images captured the channel structure of the empty, unloaded GPF cores.

While the PM deposit profile, measured by neutron imaging, down the length of the channels was similar for both E0 and E30 PM fully loaded, changes in the PM thickness profile from the different fuels differed at successive PM burnout levels. The results in Fig. 10 are condensed in Fig. 11, showing the fractional change in average PM layer thickness at 7, 20, 40, and 60 mm down the length of the GPF channels with each sequential burnout. Note that 0 mm from the front of the GPF is in the channel plug region of the filter and hinders thickness analysis; 7 mm was selected as representative of the front end of the filter as it is past the channel plugs.

Fig. 10
figure 10

Thickness profiles, derived from neutron imaging, of the PM layers along the length of the GPF channel. The x axis indicates the length (mm) from the inlet of the GPF. Plots show the changes in the average E0 (left) and E30 (right) PM thickness profiles from the same GPF channels after successive fractions of the PM mass are removed by burnout oxidation

Fig. 11
figure 11

Simplified PM thickness profiles from the neutron-imaging results presented in Fig. 10 as increasing fractions of the PM were removed by burnout oxidation. Plots a and c display the same E0 results while plots b and d display the same E30 result in 2D and 3D plots, respectively

The un-regenerated E0 PM layer (C/C0 = 0, Fig. 10 and Fig. 11) had a thicker deposit at the back of the GPF. The PM layer then burned off proportionally down the middle length (20–60 mm) of the channels for C/C0 = 0.2 and 0.4, but the thickness was only reduced by ~ 10% for each burnout, suggesting that at least a fraction of the PM removed during burnout was located throughout the depth of the layer, potentially increasing its porosity. The neutron-imaging approach does not account for any PM located within the walls of the GPF that could also contribute to the smaller reduction in thickness. The thickness of the PM layer near the inlet of the GPF, which began thinner, dropped proportionally to the 20% mass reduction between each burnout. After removing another 20% of the E0 PM mass deposit (C/C0 = 0.6), no thickness change was observed at the inlet of the GPF, while the rest of the layer down the length of the channels showed a 30% drop in thickness after the burnout. The result of these different thickness changes along the length of the PM layer resulted in consistent PM thickness along the entire channel length after the 60% burnout. During the burnout to C/C0 = 0.8, the front of the layer was more readily removed while a roughly 20% reduction in thickness was seen down the channels. With only 20% of the E0 PM left in the GPF, the thickness of the PM layer still had nearly 30% of its initial, fully loaded thickness, suggesting the PM removed during the prior burnout came from throughout the depth of the layer or within the wall.

The thickness of the E30 PM layer behaved differently (Fig. 10 and Fig. 11) during the same set of successive burnout steps as the E0-loaded GPF. After the first 20% of the PM mass was removed, little to no change in the layer thickness was observed down the length of the channels. This suggests that either a significant increase in porosity within the PM deposit had to have occurred without the layer collapsing, or the majority of the PM removed was located in the filter walls. At C/C0 = 0.4, a uniform layer thickness was observed down the length of the channel, showing a 25% drop in thickness compared with the previous burnout. The layer thickness change from 0.4 to 0.6 burnout showed another ~ 25% drop in thickness except at the front of the channel where the thickness change was less than 20%. During the 0.8 burnout, most of the remaining PM located in the front of the channels was removed, but the fraction of PM removed down the rest of the channel gradually diminished, with no reduction in the thickness at 60 mm down the channels. These results are in contrast to a neutron-imaging study of PM layers in a diesel-loaded DPF for which Toops et al. reported [33] the PM layer thickness along the channels collapsed by a third after only 20% burnout, and the PM layers at the back of the channels, near the filter plug, burned out first. A study by Sappok et al. [43] on DPF regeneration had previously reported that the soot/PM cake reached complete regeneration only near the inlet. Toops et al. [33] suggested the presence of ash at the back of the DPF channels in Sappok’s study likely blocked oxygen access near the rear.

4 Conclusions

The results of this study suggest that GDI PM trapped in a GPF is not uniform, and the oxidative reactivity of GDI PM depends on fuel composition and exposure history, and evolves in a complex manner over the course of filter regeneration.

  • The oxidation reactivity of PM generated by a GDI engine varies significantly with the composition of the fuel. E30 generates PM with much higher oxidation reactivity than E0 gasoline, while iBu24 and iBu48 generate PM with much lower reactivities.

  • The oxidation kinetics of GDI PM depend on its exposure history. PM that has been exposed to 650 °C under inert conditions is much less reactive than PM that has been recently exposed to O2 and H2O at moderate temperatures.

  • The chemical reactivity and physical properties of GDI PM collected in a GPF continuously vary as the PM is oxidized. The activation energy for oxidation by O2 increases linearly as the PM is oxidized, indicating that overall reactivity of the PM is decreasing as the PM is consumed. The active surface area of the PM also increases linearly as the PM is oxidized, contrary to what might be expected based on the reactivity trends.

  • The variations in chemical and physical properties with oxidation indicate that GDI PM behaves very differently from diesel PM.

  • The thickness of the PM layer in the GPF channel does not decrease uniformly over the course of oxidation. The observed variations in PM layer thickness are consistent with multiple stages of GPF regeneration that could include preferential oxidation of PM embedded in the filter wall; increasing porosity of the PM layer; and subsequent collapse of the more porous layer. The evolution of the PM layer during GPF regeneration is different for different fuels.

While a GPF used in the aftertreatment system of a GDI vehicle may not be expected to accumulate as much of a PM layer [23] as was present in these studies, the intrinsic kinetics of GDI PM oxidation, and how they vary with fuel composition and with the progression of the regeneration process, should provide useful insights to designers of control strategies for engines equipped with GPFs. Furthermore, the decreasing reactivity of GDI PM as oxidation progresses could result in accumulation of very unreactive PM within GPFs over certain drive cycles. The findings from this study will help in the development of active regeneration strategies for such instances. Studies that compare the PM reactivities presented here to PM generated during the first seconds of a cold start and/or kinetics from the use of other oxidants would further develop the understanding of GDI PM emissions.