Search Results

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.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

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.

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.

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).

Engine experiments have revealed the importance of fuel condensation on the emission characteristics of low temperature combustion. However, direct in-cylinder experimental evidence has not been reported in the literature. In this paper, the in-cylinder condensation processes observed in optically accessible engine experiments are first illustrated. The observed condensation processes are then simulated using state-of-the-art multidimensional engine CFD simulations with a phase transition model that incorporates a well-validated phase equilibrium numerical solver, in which a thermodynamically consistent phase equilibrium analysis is applied to determine when mixtures become unstable and a new phase is formed. The model utilizes fundamental thermodynamics principles to judge the occurrence of phase separation or combination by minimizing the system Gibbs free energy.

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.

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.

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.

Engine experiments were conducted on a heavy-duty single-cylinder engine to explore the effects of charge preparation, fuel stratification, and premixed fuel chemistry on the performance and emissions of Reactivity Controlled Compression Ignition (RCCI) combustion. The experiments were conducted at a fixed total fuel energy and engine speed, and charge preparation was varied by adjusting the global equivalence ratio between 0.28 and 0.35 at intake temperatures of 40°C and 60°C. With a premixed injection of isooctane (PRF100), and a single direct-injection of n-heptane (PRF0), fuel stratification was varied with start of injection (SOI) timing. Combustion phasing advanced as SOI was retarded between -140° and -35°, then retarded as injection timing was further retarded, indicating a potential shift in combustion regime. Peak gross efficiency was achieved between -60° and -45° SOI, and NOx emissions increased as SOI was retarded beyond -40°, peaking around -25° SOI.