Search Results

Recent research on thermal reciprocating engines has focused on the influence of lubricant oil on the combustion process, which can lead to highly undesired super-knock events. Low-Speed Pre-Ignition (LSPI) events severely limit the further development of Direct Injection Spark Ignition Engines (DISI), preventing high efficiencies from being achieved. However, there is still a lack of knowledge about the fundamental mechanisms leading to LSPI, due to the complex phenomena involved and the interaction between lubricant oil and fuel. Understanding how the presence of lubricant oil traces affects gasoline chemical reactivity is an essential step for performing successful numerical simulations aimed at predicting the onset of LSPI phenomena. Reaction mechanisms able to predict oil-fuel interaction have been proposed, but they are computationally demanding.

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.

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

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.

Although turbocharging can extend the high load limit of low temperature combustion (LTC) strategies such as reactivity controlled compression ignition (RCCI), the low exhaust enthalpy prevalent in these strategies necessitates the use of high exhaust pressures for improving turbocharger efficiency, causing high pumping losses and poor fuel economy. To mitigate these pumping losses, the divided exhaust period (DEP) concept is proposed. In this concept, the exhaust gas is directed to two separate manifolds: the blowdown manifold which is connected to the turbocharger and the scavenging manifold that bypasses the turbocharger. By separately actuating the exhaust valves using variable valve actuation, the exhaust flow is split between two manifolds, thereby reducing the overall engine backpressure and lowering pumping losses. In this paper, results from zero-dimensional and one-dimensional simulations of a multicylinder RCCI light-duty engine equipped with DEP are presented.

Part I provides a review of present designs, advanced concepts and technologies needed for the next generation of high power density diesels for midrange truck application. Part II emphasizes the use of advanced computational tools to determine the most promising design parameters for high power density diesels which are suitable for truck application. Three types of analyses are reviewed and deemed necessary for design optimization at high power density conditions, these include: I. High BMEP cycle analysis based on zero dimensional thermodynamic analysis and engine similarity rules. II. Matching of advanced turbocharging systems for high torque rise engines. III. Multi-dimensional in-cylinder computational combustion fluid dynamics (CCFD). This analytical exercise indicates that computational tools are capable of rationalizing experimental results and therefore predicting engine performance in the conceptual stage.

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.

In the current, numerical study RCCI combustion and emission characteristics using various fuel strategies are investigated, including methanol, ethanol, n-butanol and gasoline as the low reactivity fuel, and diesel fuel as the high reactivity fuel. A reduced Primary Reference Fuel (PRF)-alcohol chemical kinetic mechanism was coupled with a computational fluid dynamic (CFD) code to predict RCCI combustion under various operating conditions. The results show that a higher quantity of diesel was required to maintain the same combustion phasing with alcohol-diesel fuel blends, and the combustion durations and pressure rise rates of methanol-diesel (MD) and ethanol-diesel (ED) cases were much shorter and higher than those of gasoline-diesel (GD) and n-butanol-diesel (nBD) cases. The simulations also investigated the sensitivities of the direct injection strategies, intake temperature and premixed fuel ratio on RCCI combustion phasing control.

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.

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.

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.

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.

Gas-phase in-cylinder mixing was examined by two different methods. The first method for observing mixing involved planar Mie scattering measurements of the instantaneous number density of silicon oil droplets which were introduced to the in-cylinder flow. The local value of the number density was assumed to be representative of the local gas concentration. Because the objective was to observe the rate in which gas concentration gradients change, to provide gradients in number density, droplets were admitted into the engine through only one of the two intake ports. Air only flowed through the other port. Three different techniques were used in analyzing the droplet images to determine the spatially dependent particle number density. Direct counting, a filtering technique, and autocorrelation were used and compared. Further, numerical experiments were performed with the autocorrelation method to check its effectiveness for determination of particle number density.

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.

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.

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.

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.