Search Results

The incorporation of detailed chemistry models in internal combustion engine simulations is becoming mandatory as local, globally lean, low-temperature combustion strategies are setting the path towards a more efficient and environmentally sustainable use of energy resources in transportation. In this paper, we assessed the computational efficiency of a recently developed sparse analytical Jacobian chemistry solver, namely ‘SpeedCHEM’, that features both direct and Krylov-subspace solution methods for maximum efficiency for both small and large mechanism sizes. The code was coupled with a high-dimensional clustering algorithm for grouping homogeneous reactors into clusters with similar states and reactivities, to speed-up the chemical kinetics solution in multi-dimensional combustion simulations.

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.

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.

The simulation of combustion chemistry in internal combustion engines is challenging due to the need to include detailed reaction mechanisms to describe the engine physics. Computational times needed for coupling full chemistry to CFD simulations are still too computationally demanding, even when distributed computer systems are exploited. For these reasons the present paper proposes a time scale separation approach for the integration of the chemistry differential equations and applies it in an engine CFD code. The time scale separation is achieved through the estimation of a characteristic time for each of the species and the introduction of a sampling timestep, wherein the chemistry is subcycled during the overall integration. This allows explicit integration of the system to be carried out, and the step size is governed by tolerance requirements.

In a recent study, quantitative measurements were presented of in-cylinder spatial distributions of mixture equivalence ratio in a single-cylinder light-duty optical diesel engine, operated with a non-reactive mixture at conditions similar to an early injection low-temperature combustion mode. In the experiments a planar laser-induced fluorescence (PLIF) methodology was used to obtain local mixture equivalence ratio values based on a diesel fuel surrogate (75% n-heptane, 25% iso-octane), with a small fraction of toluene as fluorescing tracer (0.5% by mass). Significant changes in the mixture's structure and composition at the walls were observed due to increased charge motion at high swirl and injection pressure levels. This suggested a non-negligible impact on wall heat transfer and, ultimately, on efficiency and engine-out 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.

The paper presents the development of a novel approach to the solution of detailed chemistry in internal combustion engine simulations, which relies on the analytical computation of the ordinary differential equations (ODE) system Jacobian matrix in sparse form. Arbitrary reaction behaviors in either Arrhenius, third-body or fall-off formulations can be considered, and thermodynamic gas-phase mixture properties are evaluated according to the well-established 7-coefficient JANAF polynomial form. The current work presents a full validation of the new chemistry solver when coupled to the KIVA-4 code, through modeling of a single cylinder Caterpillar 3401 heavy-duty engine, running in two-stage combustion mode.

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.

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.

In this work computational and experimental approaches are combined to characterize in-cylinder flow structures and local flow field properties during operation of the Sandia 1.9L light-duty optical Diesel engine. A full computational model of the single-cylinder research engine was used that considers the complete intake and exhaust runners and plenums, as well as the adjustable throttling devices used in the experiments to obtain different swirl ratios. The in-cylinder flow predictions were validated against an extensive set of planar PIV measurements at different vertical locations in the combustion chamber for different swirl ratio configurations. Principal Component Analysis was used to characterize precession, tilting and eccentricity, and regional averages of the in-cylinder turbulence properties in the squish region and the piston bowl.

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.