Search Results

A numerical investigation of a six-stroke direct injection compression ignition engine operation in a low temperature combustion (LTC) regime is presented. The fuel employed is a gasoline-like oxygenated fuel consisting of 90% isobutanol and 10% diethyl ether (DEE) by volume to match the reactivity of conventional gasoline with octane number 87. The computational simulations of the in-cylinder processes were performed using a high-fidelity multidimensional in-house 3D CFD code (MTU-MRNT) with improved spray-sub models and CHEMKIN library. The combustion chemistry was described using a two-component (isobutanol and DEE) fuel model whose oxidation pathways were given by a reaction mechanism with 177 species and 796 reactions.

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.

The aim of this paper is to computationally investigate the combustion behavior and energy recovery processes of a six-stroke gasoline compression ignition (6S-GCI) engine that employs a continuously variable valve duration (CVVD) technique, under highly diluted, low-temperature combustion (LTC) conditions. The effects of variation of parameters concerning injection spray targeting (number of fuel injector holes. injector nozzle size and spray included angle) and combustion chamber geometry (piston bowl design) are analyzed using an in-house 3D CFD code coupled with high-fidelity physical sub-models with the Chemkin library in conjunction with a skeletal chemical kinetics mechanism for a 14-component gasoline surrogate fuel.

In order to extend the operability limit of the gasoline compression ignition (GCI) engine, as an avenue for low temperature combustion (LTC) regime, the effects of parametric variations of engine operating conditions on the performance of six-stroke GCI (6S-GCI) engine cycle are numerically investigated, using an in-house 3D CFD code coupled with high-fidelity physical sub-models along with the Chemkin library. The combustion and emissions were calculated using a skeletal chemical kinetics mechanism for a 14-component gasoline surrogate fuel. Authors’ previous study highlighted the effects of the variation of injection timing and split ratio on the overall performance of 6S-GCI engine and the unique mixing-controlled burning mode of the charge mixtures during the two additional strokes. As a continuing effort, the present study details the parametric studies of initial gas temperature, boost pressure, fuel injection pressure, compression ratio, and EGR ratio.

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.

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.

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.

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.

Three time integration schemes and four finite volume interpolation schemes for the convection term in momentum equation were tested under turbulent planar gas jet and Sandia non-reacting vaporizing Spray-H cases. The three time integration schemes are the first-order Euler implicit scheme, the second-order backward scheme, and the second-order Crank-Nicolson scheme. The four spatial interpolation schemes are cubic central, linear central, upwind, and vanLeer schemes. Velocity magnitude contour, centerline and radial mean velocity and Reynolds stress profiles for the planar turbulent gas jet case, and fuel vapor contour and liquid and vapor penetrations for the Diesel spray case predicted by the different numerical schemes were compared. The sensitivity of the numerical schemes to mesh resolution was also investigated. The non-viscosity based dynamic structure subgrid model was used. The numerical tool used in this study was OpenFOAM.

A novel 2-zone combustion chamber concept (patent pending) was developed using multi-dimensional modeling. At minimum volume, an axial projection in the piston divides the volume into distinct zones joined by a communication channel. The projection provides a means to control the mixture formation and combustion phasing within each zone. The novel combustion system was applied to reactivity controlled compression ignition (RCCI) combustion in both light-duty and heavy-duty diesel engines. Results from the study of an 8.8 bar BMEP, 2600 RPM operating condition are presented for the light-duty engine. The results from the heavy-duty engine are at an 18.1 bar BMEP, 1200 RPM operating condition. The effect of several major design features were investigated including the volume split between the inner and outer combustion chamber volumes, the clearance (squish) height, and the top ring land (crevice) volume.

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

The paper reports the application of Proper Orthogonal Decomposition (POD) to LES calculations for the analysis of combustion and knock tendency in a highly downsized turbocharged GDI engine that is currently under production. In order to qualitatively match the cyclic variability of the combustion process, Large-Eddy Simulation (LES) of the closed-valve portion of the cycle is used with cycle-dependent initial conditions from a previous multi-cycle analysis [1, 2, 3]. Detailed chemical modelling of fuel's auto-ignition quality is considered through an ad-hoc implemented look-up table approach, as a trade-off between the need for a reasonable representation of the chemistry and that of limiting the computational cost of the LES simulations. Experimental tests were conducted operating the engine at knock-limited spark advance (KLSA) and the proposed knock model was previously validated for such engine setup [3].

A Three-Way Catalyst (TWC) model with global TWC kinetics for lean burn DISI engines were developed and validated. The model incorporates kinetics of hydrocarbons and carbon monoxide oxidations, NOx reduction, water-gas and steam reforming and oxygen storage. Ammonia (NH₃) and new nitrous oxide (N₂O) kinetics were added into the model to study NH₃ and N₂O formation in TWC systems. The model was validated over a wide range of engine operating conditions using various types of experimental data from a lean burn automotive SIDI engine. First, well-controlled time-resolved steady state data were used for calibration and initial model tests. In these steady state operations, the engine was switched between lean and rich conditions for NOx emission control. Then, the model was further validated using a large set of time-averaged steady state data. Temperature dependencies of NH₃ and N₂O kinetics in the TWC model were examined and well captured by the model.

An investigation of high speed direct injection (DI) compression ignition (CI) engine combustion fueled with gasoline injected using a triple-pulse strategy in the low temperature combustion (LTC) regime is presented. This work aims to extend the operation ranges for a light-duty diesel engine, operating on gasoline, that have been identified in previous work via extended controllability of the injection process. The single-cylinder engine (SCE) was operated at full load (16 bar IMEP, 2500 rev/min) and computational simulations of the in-cylinder processes were performed using a multi-dimensional CFD code, KIVA-ERC-Chemkin, that features improved sub-models and the Chemkin library. The oxidation chemistry of the fuel was calculated using a reduced mechanism for primary reference fuel combustion chosen to match ignition characteristics of the gasoline fuel used for the SCE experiments.

Using non-viscosity dynamic structure Large Eddy Simulations (LES) turbulence model, spray=induced turbulence is investigated on a number of different Computational Fluid Dynamics (CFD) grids of varying mesh sizes (from 0.5 to 2 mm mesh). Turbulent flow is induced inside a quiescent chamber by liquid fuel spray and then left to decay after end of injection by virtue of its molecular viscosity and turbulent dissipation. Coherent structures (CS) of this turbulent flow are constructed and visualized using λ2 definition. Using CS, analysis is performed on the turbulent flow around the liquid spray jet. These CS from LES are then compared against the results from RANS calculations as well. The visualization of CS helps to explain the mechanism of fuel-air mixing obtained from LES results and its difference with RANS calculations.