Search Results

We propose a novel dual-fuel combustion model for simulating heavy-duty engines with the G-Equation. Dual-Fuel combustion strategies in such engines features direct injection of a high-reactivity fuel into a lean, premixed chamber which has a high resistance to autoignition. Distinct combustion modes are present: the DI fuel auto-ignites following chemical ignition delay after spray vaporization and mixing; a reactive front is formed on its surroundings; it develops into a well-structured turbulent flame, which propagates within the premixed charge. Either direct chemistry or the flame-propagation approach (G- Equation), taken alone, do not produce accurate results. The proposed Dual-Fuel model decides what regions of the combustion chamber should be simulated with either approach, according to the local flame state; and acts as a “kernel” model for the G- Equation model. Direct chemistry is run in the regions where a premixed front is not present.

As the global economy grows, so does the demand for heavy-duty commercial vehicles, both on-road and off-road. Currently, these vehicles are powered almost entirely by diesel engines. There is an imminent need to reduce the greenhouse gases (GHG) from this growing sector, but alternatives to the internal combustion engine face many challenges and can increase GHG emissions. For example, through simple analysis, this work will show that a Class 8 long haul on-highway truck powered entirely by battery electrics and charged from the average US electrical grid, yields significantly higher CO2 emissions per ton-mile as compared to an engine using alternative fuels. Thus, the most pragmatic and impactful way to reduce GHG emissions in commercial vehicles is using low carbon alternative fuels, such as ethanol made from renewable sources.

At recent American Society for Testing and Materials (ASTM) Subcommittee D02.01 meetings, committee members and attendees from the petroleum industry have reiterated a longstanding desire to see precision improvements to the ASTM D613 Standard Test Method for Cetane Number of Diesel Fuel Oil. The existing ASTM D613 precision limits were calculated using ASTM National Exchange Group (NEG) monthly test data from the mid-1970s through the early 1990s. Over the past few decades, many detailed studies were performed to identify and better understand the shortcomings of the cetane method (both engine equipment and instrumentation). Many of these studies concluded that inconsistent combustion is the main contributing factor behind the lack of precision in the cetane number method, followed by shortcomings in the instrumentation used to measure ignition delay.

This work explores the effectiveness of common rail fuel injectors equipped with Grouped Hole Nozzles (GHNs) in aiding the mixing process and reducing particulate matter (PM) emissions of Conventional Diesel Combustion (CDC) engines, while maintaining manageable Oxides of Nitrogen (NOx) levels. Parallel (pGHN), converging (cGHN) and diverging (dGHN) - hole GHNs were studied and the results were compared to a conventional, single hole nozzle (SHN) with the same flow area. The study was conducted on a single cylinder medium-duty engine to isolate the effects of the combustion from multi-cylinder effects and the conditions were chosen to be representative of a typical mid-load operating point for an on-road diesel engine. The effects of injection pressure and the Start of Injection (SOI) timing were explored and the tradeoffs between these boundary conditions are examined by using a response surface fitting technique, to identify an optimum operating condition.

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.

Dual-fuel reactivity controlled compression ignition (RCCI) combustion is a promising method to achieve high efficiency with near-zero NOx and soot emissions; however, the requirement to carry two fuels on board limits practical application. Advancements in catalytic reforming have demonstrated the ability to generate syngas (a mixture of CO and hydrogen) from a single hydrocarbon stream. This syngas mixture can then be used as the low reactivity fuel stream to enable single fuel RCCI combustion. The present effort uses a combination of engine experiments and system level modeling to investigate reformed fuel RCCI combustion. The impact of reformer composition is investigated by varying the syngas composition from 10% H2 to approximately 80% H2. The results of the investigation show that reformed fuel RCCI combustion is possible over a wide range of H2/CO ratios.

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.

Engine experiments were carried out on a heavy-duty single-cylinder engine to investigate the effects of Gasoline Compression Ignition on emissions and performance of a heavy-duty engine operating at a high load condition. Comparisons between gasoline fueled operation and diesel fueled operation are presented using a single, near top dead center injection. Although the fuel’s cetane numbers are very different, the combustion characteristics of the two fuels at high load are similar, with the gasoline-fueled case showing less than two crank angle degree longer ignition delay. Gasoline operation showed lower soot production at similar levels of NOx, initiating study of the impact of exhaust gas recirculation which spanned a range of NOx levels covering the range from minimal urea dosing to high urea dosing. A conventional soot-NOx tradeoff was found to exist with gasoline as exists with diesel.

Gasoline compression ignition (GCI) combustion is a promising solution to address increasingly stringent efficiency and emissions regulations imposed on the internal combustion engine. However, the high resistance to auto-ignition of modern market gasoline makes low load compression ignition (CI) operation difficult. Accordingly, a method that enables the variation of the fuel reactivity on demand is an ideal solution to address low load stability issues. Metal engine experiments conducted on a single cylinder medium-duty research engine allowed for the investigation of this strategy. The fuels used for this study were 87 octane gasoline (primary fuel stream) and diesel fuel (reactivity enhancer). Initial tests demonstrated load extension down to idle conditions with only 20% diesel by mass, which reduced to 0% at loads above 3 bar IMEPg.

Advancements in catalytic reforming have demonstrated the ability to generate syngas (a mixture of CO and hydrogen) from a single hydrocarbon stream. This syngas mixture can then be used to replace diesel fuel and enable dual-fuel combustion strategies. The role of port-fuel injected syngas, comprised of equal parts hydrogen and carbon monoxide by volume was investigated experimentally for soot reduction benefits under a transient load change at constant speed. The syngas used for the experiments was presumed to be formed via a partial oxidation on-board fuel reforming process and delivered through gaseous injectors using a custom gas rail supplied with bottle gas, mounted in the swirl runner of the intake manifold. Time-based ramping of the direct-injected fuel with constant syngas fuel mass delivery from 2 to 8 bar brake mean effective pressure was performed on a multi-cylinder, turbocharged, light-duty engine to determine the effects of syngas on transient soot emissions.

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.

A genetic algorithm (GA) optimization methodology is applied to the design of the combustion system of a heavy-duty (HD) Diesel engine fueled with dimethyl ether (DME). The study has two objectives, the optimization of a conventional diffusion-controlled combustion system aiming to achieve US2010 targets and the optimization of a stoichiometric combustion system coupled with a three way catalyst (TWC) to further control NOx emissions and achieve US2030 emission standards. These optimizations include the key combustion system related hardware, bowl geometry and injection nozzle design as input factors, together with the most relevant air management and injection settings. The GA was linked to the KIVA CFD code and an automated grid generation tool to perform a single-objective optimization. The target of the optimizations is to improve net indicated efficiency (NIE) while keeping NOx emissions, peak pressure and pressure rise rate under their corresponding target levels.

The design of modern diesel-powered vehicles involves optimization and balancing of trade-offs for fuel efficiency, emissions, and noise. To meet increasingly stringent emission regulations, diesel powertrains employ aftertreatment devices to control nitrogen oxides, hydrocarbons, carbon monoxide, and particulate matter emissions and use active exhaust warm-up strategies to ensure those devices are active as quickly as possible. A typical strategy for exhaust warm-up is to operate with retarded combustion phasing, limited by combustion stability and HC emissions. The amount of exhaust enthalpy available for catalyst light-off is limited by the extent to which combustion phasing can be retarded. Diesel cetane number (CN), a measure of fuel ignition quality, has an influence on combustion stability at retarded combustion phasing. Diesel fuel in the United States tends to have a lower CN (both minimum required and average in market) than other countries.

Large-eddy simulation (LES) is a useful approach for the simulation of turbulent flow and combustion processes in internal combustion engines. This study employs the ANSYS Forte CFD package and explores several key and fundamental components of LES, namely, the subgrid-scale (SGS) turbulence models, the numerical schemes used to discretize the transport equations, and the computational mesh. The SGS turbulence models considered include the classic Smagorinsky model and a dynamic structure model. Two numerical schemes for momentum convection, quasi-second-order upwind (QSOU) and central difference (CD), were evaluated. The effects of different computational mesh sizes controlled by both fixed mesh refinement and a solution-adaptive mesh-refinement approach were studied and compared. The LES models are evaluated and validated against several flow configurations that are critical to engine flows, in particular, to fuel injection processes.

In this paper, large eddy simulation (LES) coupled with two uncertainty quantification (UQ) methods, namely latin-hypercube sampling (LHS) and polynomial chaos expansion (PCE), have been used to quantify the effects of model parameters and spray boundary conditions on diesel and gasoline spray simulations. Evaporating, non-reacting spray data was used to compare penetration, mixture fraction and spray probability contour. Two different sets of four uncertain variables were used for diesel and gasoline sprays, respectively. UQ results showed good agreement between experiments and predictions. UQ statistics indicated that discharge coefficient has stronger impact on gasoline than diesel sprays, and spray cone angle is important for vapor penetration of both types of sprays. Additionally, examination of the gasoline spray characteristics showed that plume-to-plume interaction and nozzle dribble are important phenomena that need to be considered in high-fidelity gasoline spray simulations.

A correlation was developed to predict the ignition delay of PRF blends at a wide range of engine-relevant operating conditions. Constant volume simulations were performed using Cantera coupled with a reduced reaction mechanism at a range of initial temperatures from 570-1860K, initial pressures from 10-100atm, oxygen mole percent from 12.6% to 21%, equivalence ratios from 0.30-1.5, and PRF blends from PRF0 to PRF100. In total, 6,480 independent ignition delay simulations were performed. The correlation utilizes the traditional Arrhenius formulation; with equivalence ratio (φ), pressure (p), and oxygen mole percentage (xo2) dependencies. The exponents α, β, and γ were fitted to a third order polynomial with respect to temperature with an exponential roll-off to a constant value at low temperatures to capture the behavior expressed by the reaction mechanism. The location and rate of the roll-off functions were modified by linear functions of PRF.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

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.

This paper is part of a larger body of experimental and computational work devoted to studying the role of close-coupled post injections on soot reduction in a heavy-duty optical engine. It is a continuation of an earlier computational paper. The goals of the current work are to develop new CFD analysis tools and methods and apply them to gain a more in depth understanding of the different in-cylinder environments into which fuel from main- and post-injections are injected and to study how the in-cylinder flow, thermal and chemical fields are transformed between start of injection timings. The engine represented in this computational study is a single-cylinder, direct-injection, heavy-duty, low-swirl engine with optical components. It is based on the Cummins N14, has a cylindrical shaped piston bowl and an eight-hole injector that are both centered on the cylinder axis. The fuel used was n-heptane and the engine operating condition was light load at 1200 RPM.

Due to the new challenge of meeting number-based regulations for particulate matter (PM), a numerical and experimental study has been conducted to better understand particulate formation in engines fuelled with compressed natural gas. The study has been conducted on a Heavy-Duty, Euro VI, 4-cylinder, spark ignited engine, with multipoint sequential phased injection and stoichiometric combustion. For the experimental measurements two different instruments were used: a condensation particle counter (CPC) and a fast-response particle size spectrometer (DMS) the latter able also to provide a particle size distribution of the measured particles in the range from 5 to 1000 nm. Experimental measurements in both stationary and transient conditions were carried out. The data using the World Harmonized Transient Cycle (WHTC) were useful to detect which operating conditions lead to high numbers of particles. Then a further transient test was used for a more detailed and deeper analysis.