Search Results

The steam reforming of CH4 plays a crucial role in the high-temperature activity of natural gas three-way catalysts. Despite existing reports on sulfur inhibition in CH4 steam reforming, there is a limited understanding of sulfur storage and removal dynamics under various lambda conditions. In this study, we utilize a 4-Mode sulfur testing approach to elucidate the dynamics of sulfur storage and removal and their impact on three-way catalyst performance. We also investigate the influence of sulfur on CH4 steam reforming by analyzing CH4 conversions under dithering, rich, and lean reactor conditions. In the 4-Mode sulfur test, saturating the TWC with sulfur at low temperatures emerges as the primary cause of significant three-way catalyst performance degradation. After undergoing a deSOx treatment at 600 °C, NOx conversions were fully restored, while CH4 conversions did not fully recover.

The prediction accuracy of a three-way catalyst (TWC) model is highly associated with the ability of the model to incorporate the reaction kinetics of the emission process as a lambda function. In this study, we investigated the O2 and H2 concentration profiles of TWC reactions and used them as critical inputs for the development of a global TWC model. We presented the experimental data and global kinetic model showing the impact of thermal degradation on the performance of the TWC. The performance metrics investigated in this study included CH4, NOx, and CO conversions under lean, rich, and dithering light-off conditions to determine the kinetics of oxidation reactions and reduction/reforming/water-gas shift reactions as a function of thermal aging. The O2 and H2 concentrations were measured using mass spectrometry to track the change in the oxidation state of the catalyst and to determine the mechanism of the reactions under these light-off conditions.

Stoichiometric natural gas engines with three-way catalysts emit less NOx and CH4 due to their higher efficiency compared to lean-burn natural gas engines. Although hydrothermal aging of three-way catalysts has been extensively studied, a deeper understanding beyond hydrothermal aging is needed to explain real-world performance, especially for natural gas engines with near-zero NOx emissions. In this investigation, field-aged three-way catalysts were characterized to identify the contribution of chemical aging to their overall performance. It was found that the sulfur species on the field-aged TWCs were entirely distributed along the catalyst length, showing a decreasing trend, whereas phosphorous contamination was mainly observed at the inlet section of the three-way catalysts, and the phosphorous concentration declined sharply along the axial length.

G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device.

In this work we studied the effects of piston bowl design on combustion in a small-bore direct-injection diesel engine. Two bowl designs were compared: a conventional, omega-shaped bowl and a stepped-lip piston bowl. Experiments were carried out in the Sandia single-cylinder optical engine facility, with a medium-load, mild-boosted operating condition featuring a pilot+main injection strategy. CFD simulations were carried out with the FRESCO platform featuring full-geometric body-fitted mesh modeling of the engine and were validated against measured in-cylinder performance as well as soot natural luminosity images. Differences in combustion development were studied using the simulation results, and sensitivities to in-cylinder flow field (swirl ratio) and injection rate parameters were also analyzed.

Sector mesh modeling is the dominant computational approach for combustion system design optimization. The aim of this work is to quantify the errors descending from the sector mesh approach through three geometric modeling approaches to an optical diesel engine. A full engine geometry mesh is created, including valves and intake and exhaust ports and runners, and a full-cycle flow simulation is performed until fired TDC. Next, an axisymmetric sector cylinder mesh is initialized with homogeneous bulk in-cylinder initial conditions initialized from the full-cycle simulation. Finally, a 360-degree azimuthal mesh of the cylinder is initialized with flow and thermodynamics fields at IVC mapped from the full engine geometry using a conservative interpolation approach. A study of the in-cylinder flow features until TDC showed that the geometric features on the cylinder head (valve tilt and protrusion into the combustion chamber, valve recesses) have a large impact on flow complexity.

Three-way catalysts have been used in a variety of stoichiometric natural gas engines for emission control. During real-world operation, these catalysts have experienced a large number of temporary and permanent deactivations including thermal aging and chemical contamination. Thermal aging is typically induced either by high engine-out exhaust temperatures or the reaction exotherm generated on the catalysts. Chemical contamination originates from various inorganic species such as Phosphorous (P) and Sulfur (S) that contain in engine fluids, which can poison and/or mask the catalyst active components. Such deactivations are quite difficult to simulate under laboratory conditions, due to the fact that multiple deactivation modes may occur at the same time in the real-world operations. In this work, a set of field-aged TWCs has been analyzed through detailed laboratory research in order to identify and quantify the real-world aging mechanisms.

Surface film formation is an important phenomenon during spray impingement in a combustion chamber. The film that forms on the chamber walls and piston bowl produces soot post-combustion. While some droplets stick to the wall surface, others splash and interact with the gas present inside the combustion chamber. Accurate prediction of both the film thickness and splashed mass is crucial for surface film model development since it leads to a precise estimation of the amount of soot and other exhaust gases formed. This information could guide future studies aimed at a comprehensive understanding of the combustion process and might enable development of engines with reduced emissions. Dynamic structure Large Eddy Simulation (LES) turbulence model implemented for in-cylinder sprays [1] has shown to predict the flow structure of a spray more accurately than the Reynolds-averaged Navier-Stokes turbulence model.

In most large-eddy simulation (LES) applications to two-phase engine flows, the liquid-air interactions need to be accounted for as source terms in the respective governing equations. Accurate calculation of these source terms requires the relative velocity “seen” by liquid droplets as they move across the flow, which generally needs to be estimated using a turbulent dispersion model. Turbulent dispersion modeling in LES is very scarce in the literature. In most studies on engine spray flows, sub-grid scale (SGS) models for the turbulent dispersion still follow the same stochastic approach originally proposed for Reynolds-averaged Navier-Stokes (RANS). In this study, an SGS dispersion model is formulated in which the instantaneous gas velocity is decomposed into a deterministic part and a stochastic part. The deterministic part is reconstructed using the approximate deconvolution method (ADM), in which the large-scale flow can be readily calculated.

Diesel piston bowl geometry can affect turbulent mixing and therefore it impacts heat-release rates, thermal efficiency, and soot emissions. The focus of this work is on the effects of bowl geometry and injection timing on turbulent flow structure. This computational study compares engine behavior with two pistons representing competing approaches to combustion chamber design: a conventional, re-entrant piston bowl and a stepped-lip piston bowl. Three-dimensional computational fluid dynamics (CFD) simulations are performed for a part-load, conventional diesel combustion operating point with a pilot-main injection strategy under non-combusting conditions. Two injection timings are simulated based on experimental findings: an injection timing for which the stepped-lip piston enables significant efficiency and emissions benefits, and an injection timing with diminished benefits compared to the conventional, re-entrant piston.

To meet stringent emissions regulations, filtration devices are often used in engine exhaust systems to reduce particulate mass (PM) and particulate number (PN). Diesel particulate filters (DPFs) are a well-established means of reducing PM from diesel engines to meet emissions regulations. New emissions regulations will most likely require a similar technology on gasoline engines with direct injection, gasoline particulate filters (GPFs). Due to differences in the exhaust and particulate characteristics, the design and operation of GPFs and DPFs differ. In a DPF filtration is dominated by the buildup of a soot cake. Whereas in a GPF, much of the soot is trapped inside the porous substrate, or filter wall, where deep-bed filtration is dominant. Thus, an accurate model describing the porous filtration properties of GPF substrates is desired. The pore filtration model (PFM) was developed to more accurately model the deep-bed filtration process that occurs in a GPF.

A swept-volume method of calculating the volume swept by the flame during each time step is developed and used to improve the calculation of fuel reaction rates. The improved reaction rates have been applied to the ignition model and coupled with the level set G-equation combustion model. In the ignition model, a single initial kernel is formed after which the kernel is convected by the gas flow and its growth rate is determined by the flame speed and thermal expansion due to the energy transfer from the electrical circuit. The predicted ignition kernel size was compared with the available experimental data and good agreements were achieved. Once the ignition kernel reaches a size when the fully turbulent flame is developed, the G-equation model is switched on to track the mean turbulent flame front propagation.

The aging impact on oxygen storage capacity (OSC) and catalyst performance was investigated on one degreened and one aged (hydrothermally aged at 955 °C for 50 h) commercial three-way catalyst (TWC) by experiments and modeling. The difference of OSC between the degreened and aged TWCs was dependent on catalyst temperature. The largest difference was found at 600 °C, at which the amount of OSC decreased by 45.5%. Catalyst performance was evaluated through lightoff tests at two simulated engine exhaust conditions (lean and rich) on a micro-reactor. The aging impact on the catalyst performance was different under lean and rich environments and investigated separately. At the lean condition, oxidation of CO and C3H6 was significantly suppressed while oxidation of C3H8 was relatively less degraded. At the rich condition, the inhibition effect was more pronounced on the aged TWC and inhibiting hydrocarbon species from C3H6 partial oxidation can survive at temperatures up to 450 °C.

Modeling ignition kernel development in spark ignition engines is crucial to capturing the sources of cyclic variability, both with RANS and LES simulations. Appropriate kernel modeling must ensure that energy transfer from the electrodes to the gas phase has the correct timing, rate and locations, until the flame surface is large enough to be represented on the mesh by the G-Equation level-set method. However, in most kernel models, geometric details driving kernel growth are missing: either because it is described as Lagrangian particles, or because its development is simplified, i.e., down to multiple spherical flames. This paper covers the geometric aspects of kernel development, which makes up the core of a Triangulated Lagrangian Ignition Kernel model. One (or multiple, if it restrikes) spark channel is initialized as a one-dimensional Lagrangian particle thread.

One way to develop an understanding of soot formation and oxidation processes that occur during direct injection and combustion in an internal combustion engine is to image the natural luminosity from soot over time. Imaging is possible when there is optical access to the combustion chamber. After the images are acquired, the next challenge is to properly interpret the luminous distributions that have been captured on the images. A major focus of this paper is to provide guidance on interpretation of experimental images of soot luminosity by explaining how radiation from soot is predicted to change as it is transmitted through the combustion chamber and to the imaging. The interpretations are only limited by the scope of the models that have been developed for this purpose. The end-goal of imaging radiation from soot is to estimate the amount of soot that is present.

In this work, linear, non-linear and a generalized renormalization group (RNG) two-equation RANS turbulence models of the k-epsilon form were compared for the prediction of turbulent compressible flows in diesel engines. The object-oriented, multidimensional parallel code FRESCO, developed at the University of Wisconsin, was used to test the alternative models versus the standard k-epsilon model. Test cases featured the academic backward facing step and the impinging gas jet in a quiescent chamber. Diesel engine flows featured high-pressure spray injection in a constant volume vessel from the Engine Combustion Network (ECN), as well as intake flows in a high-swirl diesel engine. For the engine intake flows, a model of the Sandia National Laboratories 1.9L light-duty single cylinder optical engine was used.

Oxygen storage capacity (OSC) is one of the most critical characteristics of a three-way catalyst (TWC) and is closely related to the catalyst aging and performance. In this study, a dynamic OSC model involving two oxygen storage sites with distinct kinetics was developed. The dual-site OSC model was validated on a bench reactor and a natural gas engine. The model was capable of predicting temperature dependence on OSC with H2, CO and CH4 as reductants. Also, the effects of oxygen concentration and space velocity on the amount of OSC were captured by the model. The validated OSC model was applied to simulate lean breakthrough phenomena with varied space velocities and oxygen concentrations. It is found that OSC during lean breakthrough is not a constant for a particular TWC catalyst and is dependent on space velocity and oxygen concentration. Specifically, breakthrough time exhibits a non-linear, inverse correlation to oxygen flux.

In an attempt to increase efficiency and lower critical and highly regulated emissions (i.e., NOx, PM and CO2) many advanced combustion strategies have been investigated. Most of the current strategies fall into the category of low temperature combustion (LTC), which allow emissions mandates to be met in-cylinder along with anticipated reduction in cost and complexity. These strategies, such as homogeneous charge compression ignition (HCCI), premixed charge compression ignition (PCCI), partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI), use early injection timings, resulting in a highly lean charge with increased specific heat ratios to improve thermal efficiency and reduce PM emissions. Lower combustion temperatures also avoid the activation of NOx formation reactions.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

Low-temperature combustion (LTC) strategies have been an active area of research due to their ability to achieve high thermal efficiency while avoiding the formation of NOx and particulate matter. One of the largest challenges with LTC is the relative lack of authority over the heat release rate profile, which, depending on the particular injection strategy, either limits the maximum attainable load, or creates a tradeoff between noise and efficiency at high load conditions. We have shown previously that control over heat release can be dramatically improved through a combination of reactivity stratification in the premixed charge and a diffusion-limited injection that occurs after the conclusion of the low-temperature heat release, in a strategy called direct dual fuel stratification (DDFS).