Search Results

Spray-wall interactions in diesel engines have a strong influence on turbulent flow evolution and mixing, which influences the engine’s thermal efficiency and pollutant-emissions behavior. Previous optical experiments and numerical investigations of a stepped-lip diesel piston bowl focused on how spray-wall interactions influence the formation of squish-region vortices and their sensitivity to injection timing. Such vortices are stronger and longer-lived at retarded injection timings and are correlated with faster late-cycle heat release and soot reductions, but are weaker and shorter-lived as injection timing is advanced. Computational fluid dynamics (CFD) simulations predict that piston bowls with more space in the squish region can enhance the strength of these vortices at near-TDC injection timings, which is hypothesized to further improve peak thermal efficiency and reduce emissions. The dimpled stepped-lip (DSL) piston is such a design.

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.

A one-dimensional, non-equilibrium, compressible law of the wall model is proposed to increase the accuracy of heat transfer predictions from computational fluid dynamics (CFD) simulations of internal combustion engine flows on engineering grids. Our 1D model solves the transient turbulent Navier-Stokes equations for mass, momentum, energy and turbulence under the thin-layer assumption, using a finite-difference spatial scheme and a high-order implicit time integration method. A new algebraic eddy-viscosity closure, derived from the Han-Reitz equilibrium law of the wall, with enhanced Prandtl number sensitivity and compressibility effects, was developed for optimal performance. Several eddy viscosity sub-models were tested for turbulence closure, including the two-equation k-epsilon and k-omega, which gave insufficient performance.

Appropriate spray modeling in multidimensional simulations of diesel engines is well known to affect the overall accuracy of the results. More and more accurate models are being developed to deal with drop dynamics, breakup, collisions, and vaporization/multiphase processes; the latter ones being the most computationally demanding. In fact, in parallel calculations, the droplets occupy a physical region of the in-cylinder domain, which is generally very different than the topology-driven finite-volume mesh decomposition. This makes the CPU decomposition of the spray cloud severely uneven when many CPUs are employed, yielding poor parallel performance of the spray computation. Furthermore, mesh-independent models such as collision calculations require checking of each possible droplet pair, which leads to a practically intractable O(np2/2) computational cost, np being the total number of droplets in the spray cloud, and additional overhead for parallel communications.

Lagrangian spray modeling represents a critical boundary condition for multidimensional simulations of in-cylinder flow structure, mixture formation and combustion in internal combustion engines. Segregated models for injection, breakup, collision and vaporization are usually employed to pass appropriate momentum, mass, and energy source terms to the gas-phase solver. Careful calibration of each sub-model generally produces appropriate results. Yet, the predictiveness of this modeling approach has been questioned by recent experimental observations, which showed that at trans- and super-critical conditions relevant to diesel injection, classical atomization and vaporization behavior is replaced by a mixing-controlled phase transition process of a dense fluid. In this work, we assessed the shortcomings of classical spray modeling with respect to real-gas and phase-change behavior, employing a multicomponent phase equilibrium solver and liquid-jet theory.

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.

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.

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.

A computational fluid dynamics (CFD) guided combustion system optimization was conducted for a heavy-duty compression-ignition engine with a gasoline-like fuel that has an anti-knock index (AKI) of 58. The primary goal was to design an optimized combustion system utilizing the high volatility and low sooting tendency of the fuel for improved fuel efficiency with minimal hardware modifications to the engine. The CFD model predictions were first validated against experimental results generated using the stock engine hardware. A comprehensive design of experiments (DoE) study was performed at different operating conditions on a world-leading supercomputer, MIRA at Argonne National Laboratory, to accelerate the development of an optimized fuel-efficiency focused design while maintaining the engine-out NOx and soot emissions levels of the baseline production engine.

Heavy-duty vehicles, currently the second largest source of fuel consumption and carbon emissions are projected to be fastest growing mode in transportation sector in future. There is a clear need to increase fuel efficiency and lower emissions for these engines. The Opposed-Piston Engine (OP Engine) has the potential to address this growing need. In this paper, results are presented for a 9.8L three-cylinder two-stroke OP Engine that shows the potential of achieving 55% brake thermal efficiency (BTE), while simultaneously satisfying emission targets for tail pipe emissions. The two-stroke OP Engines are inherently more cost effective due to less engine parts. The OP Engine architecture presented in this paper can meet this performance without the use of waste heat recovery systems or turbo-compounding and hence is the most cost effective technology to deliver this level of fuel efficiency.

The current trend of downsizing used in gasoline engines, while reducing fuel consumption and CO2 emissions, imposes severe thermal loads inside the combustion chamber. These critical thermodynamic conditions lead to the possible auto-ignition (AI) of fresh gases hot-spots around Top-Dead-Center (TDC). At this very moment where the surface to volume ratio is high, wall heat transfer influences the temperature field inside the combustion chamber. The use of a realistic wall temperature distribution becomes important in the case of a downsized engine where fresh gases hot spots found near high temperature walls can initiate auto-ignition. This paper presents a comprehensive numerical methodology for an accurately prediction of thermodynamic conditions inside the combustion chamber based on Conjugate Heat Transfer (CHT).

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.

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.

This work involves modeling internal and near-nozzle flows of a gasoline direct injection (GDI) nozzle. The Engine Combustion Network (ECN) Spray G condition has been considered for these simulations using the nominal geometry of the Spray G injector. First, best practices for numerical simulation of the two-phase flow evolution inside and the near-nozzle regions of the Spray G injector are presented for the peak needle lift. The mass flow rate prediction for peak needle lift was in reasonable agreement with experimental data available in the ECN database. Liquid plume targeting angle and liquid penetration estimates showed promising agreement with experimental observations. The capability to assess the influence of different thermodynamic conditions on the two-phase flow nature was established by predicting non-flashing and flashing phenomena.

The 3-Zones Extended Coherent Flame Model (ECFM3Z) and the Tabulated Kinetics for Ignition (TKI) auto-ignition model are widely used for RANS simulations of reactive flows in Diesel engines. ECFM3Z accounts for the turbulent mixing between one zone that contains compressed air and EGR and another zone that contains evaporated fuel. These zones mix to form a reactive zone where combustion occurs. In this mixing zone TKI is applied to predict the auto-ignition event, including the ignition delay time and the heat release rate. Because it is tabulated, TKI can model complex fuels over a wide range of engine thermodynamic conditions. However, the ECFM3Z/TKI combustion modeling approach requires an efficient predictive spray injection calculation. In a Diesel direct injection engine, the turbulent mixing and spray atomization are mainly driven by the liquid/gas coupling phenomenon that occurs at moving liquid/gas interfaces.

Reynolds-averaged Navier-Stokes (RANS) modeling is expected to deliver an ensemble-averaged result for the majority of turbulent flows. This could lead to the conclusion that multi-cycle internal combustion engine (ICE) simulations performed using RANS must exhibit a converging numerical solution after a certain number of consecutive cycles. However, for some engine configurations unsteady RANS simulations are not guaranteed to deliver an ensemble-averaged result. In this paper it is shown that, when using RANS modeling to simulate multiple engine cycles, the cycle-to-cycle variations (CCV) generated from different initial conditions at each cycle are not damped out even after a large number of cycles. A single-cylinder GDI research engine is simulated using RANS modeling and the numerical results for 20 consecutive engine cycles are evaluated for two specific operating conditions.

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.

Reactivity Controlled Compression Ignition (RCCI) has been shown to be an attractive concept to achieve clean and high efficiency combustion. RCCI can be realized by applying two fuels with different reactivities, e.g., diesel and gasoline. This motivates the idea of using a single low reactivity fuel and direct injection (DI) of the same fuel blended with a small amount of cetane improver to achieve RCCI combustion. In the current study, numerical investigation was conducted to simulate RCCI and HCCI combustion and emissions with various fuels, including gasoline/diesel, iso-butanol/diesel and iso-butanol/iso-butanol+di-tert-butyl peroxide (DTBP) cetane improver. A reduced Primary Reference Fuel (PRF)-iso-butanol-DTBP mechanism was formulated and coupled with the KIVA computational fluid dynamic (CFD) code to predict the combustion and emissions of these fuels under different operating conditions in a heavy duty diesel engine.