Search Results

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

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.

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 paper reports the application of Proper Orthogonal Decomposition (POD) to LES calculations for the analysis of combustion and knock tendency in a highly downsized turbocharged GDI engine that is currently under production. In order to qualitatively match the cyclic variability of the combustion process, Large-Eddy Simulation (LES) of the closed-valve portion of the cycle is used with cycle-dependent initial conditions from a previous multi-cycle analysis [1, 2, 3]. Detailed chemical modelling of fuel's auto-ignition quality is considered through an ad-hoc implemented look-up table approach, as a trade-off between the need for a reasonable representation of the chemistry and that of limiting the computational cost of the LES simulations. Experimental tests were conducted operating the engine at knock-limited spark advance (KLSA) and the proposed knock model was previously validated for such engine setup [3].

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.

Fundamental simulations using DNS type procedures were used to investigate the ignition, combustion characteristics and the lift-off trends of a spatially evolving turbulent liquid fuel jet. In particular, the spatially evolving n-Heptane spray injected in a two-dimensional rectangular domain with an engine like environment was investigated. The computational results were compared to the experimental observations from an optical engine as reported in the literature. It was found that an initial fuel rich combustion downstream of the spray tip is followed by diffusion combustion. Investigations were also made to understand the effects of injection velocity, ambient temperature and the droplet radius on the lift-off length. For each of these parameters three different values in a given range were chosen. For both injection velocity and droplet radius, an increase resulted in a near linear increase in the lift-off length.

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.

While Low Temperature Combustion (LTC) strategies such as Reactivity Controlled Compression Ignition (RCCI) exhibit high thermal efficiency and produce low NOx and soot emissions, low load operation is still a significant challenge due to high unburnt hydrocarbon (UHC) and carbon monoxide (CO) emissions, which occur as a result of poor combustion efficiencies at these operating points. Furthermore, the exhaust gas temperatures are insufficient to light-off the Diesel Oxidation Catalyst (DOC), thereby resulting in poor UHC and CO conversion efficiencies by the aftertreatment system. To achieve exhaust gas temperature values sufficient for DOC light-off, combustion can be appropriately phased by changing the ratio of gasoline to diesel in the cylinder, or by burning additional fuel injected during the expansion stroke through post-injection.

Diesel particulate samples were collected from a light duty engine operated at a single speed-load point with a range of biodiesel and conventional fuel blends. The oxidation reactivity of the samples was characterized in a laboratory reactor, and BET surface area measurements were made at several points during oxidation of the fixed carbon component of both types of particulate. The fixed carbon component of biodiesel particulate has a significantly higher surface area for the initial stages of oxidation, but the surface areas for the two particulates become similar as fixed carbon oxidation proceeds beyond 40%. When fixed carbon oxidation rates are normalized to total surface area, it is possible to describe the oxidation rates of the fixed carbon portion of both types of particulates with a single set of Arrhenius parameters. The measured surface area evolution during particle oxidation was found to be inconsistent with shrinking sphere oxidation.

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.

Fundamental simulations in a quiescent cell under adiabatic conditions were made to understand the effect of temperature, equivalence ratio and the components of the recirculated exhaust gas, viz., CO2 and H2O, on the combustion of n-Heptane. Simulations were made in single phase in which evaporated n-Heptane was uniformly distributed in the domain. Computations were made for two different temperatures and four different EGR levels. CO2 or H2O or N2was used as EGR. It was found that the initiation of the main combustion process was primarily determined by two competing factors, i.e., the amount of initial OH concentration in the domain and the specific heat of the mixture. Further, initial OH concentration can be controlled by the manipulating the ambient temperature in the domain, and the specific heat capacity of the mixture via the mixture composition. In addition to these, the pre combustion and the subsequent post combustion can also be controlled via the equivalence ratio.

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.

Detailed knowledge of hydrocarbon fuel combustion chemistry has grown tremendously in recent years. However, the gap between detailed chemistry and computational fluid dynamics (CFD) remains, because of the high cost of solving detailed chemistry in a large number of computational cells. This paper presents the results of applying a suite of techniques aimed at closing this gap. The techniques include use of a surrogate blend optimizer and a guided mechanism reduction methodology, as well as advanced methods for efficiently and accurately coupling the pre-reduced kinetic models with the multidimensional transport equations. The advanced methods include dynamic adaptive chemistry (DAC) and dynamic cell clustering (DCC) algorithms.

Surface film formation is an important phenomenon during spray impingement in a combustion chamber. The film that forms on the chamber walls and piston bowl produces soot post-combustion. While some droplets stick to the wall surface, others splash and interact with the gas present inside the combustion chamber. Accurate prediction of both the film thickness and splashed mass is crucial for surface film model development since it leads to a precise estimation of the amount of soot and other exhaust gases formed. This information could guide future studies aimed at a comprehensive understanding of the combustion process and might enable development of engines with reduced emissions. Dynamic structure Large Eddy Simulation (LES) turbulence model implemented for in-cylinder sprays [1] has shown to predict the flow structure of a spray more accurately than the Reynolds-averaged Navier-Stokes turbulence model.

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.

Improved wall models for LES are presented in this paper. The classical Werner-Wengle (WW) wall shear stress model was used along with the eddy viscosity near walls. A sub-grid scale turbulent kinetic energy was employed in a model for the eddy viscosity. To provide heat flux results, a modified classical variable-density wall heat transfer model, which includes the variation of the turbulent Prandtl number in the boundary layer, was also employed. The fully turbulent developed flow in a square duct with constant wall temperature was used to validate the model and the friction factor and Nusselt number predictions are in good agreement with experimental results. The resulting time and spatially averaged velocity and temperature wall functions from the new wall models match well with the law-of-the-wall experimental data. Additionally, the model was validated using experimental data from a Caterpillar engine operated with conventional diesel combustion.

A modified version of the three dimensional CFD code KIVA-3, which accommodates moving valves and a moving piston crown, has been applied to a heavy-duty, four cycle, dual intake valve, direct injection, diesel engine The fluid domain encompasses an intake runner, two valved ports and a cylinder with a Mexican hat bowl-in-piston configuration In the first part of the study, the modified KIVA-3 code was used to simulate the flow through the port and cylinder of the engine without a piston The intake valves were cycled (300 rpm at the camshaft) Two turbulence models were compared in this part of the study, standard k - ε and the RNG modified k - ε as discussed in [11] The results were compared to particle image velocimetry (PIV) images Large scale flow features of the computer simulations agreed moderately well with the ensemble averaged flow bench results There was very little difference between the results from the two turbulence models Motored engine simulations including a piston were conducted to characterize the in-cylinder gas flow and mixing during the intake and compression strokes Characterization methods were developed which yield insight into the in-cylinder gas motion and mixing during both strokes Intake and residual gases are tracked separately, both large scale convection and turbulent mixing are investigated, flow critical points are examined to provide information about flow topology and turbulence production is correlated with the evolution of flow structures Results show that mixing of the intake and residual gases is very non-uniform Many complex flow structures develop during intake and are destroyed during compression However, several structures survive through compression and contribute to enhanced mixing near top dead center These significant structures have been identified and tracked back to intake The flow field near top dead center exhibits spatial inhomogeneities in temperature and small scale mixing parameters such as turbulence kinetic energy and its dissipation rate

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.

The paper reports the application of a look-up table approach within a LES combustion modelling framework for the prediction of knock limit in a highly downsized turbocharged DISI engine. During experimental investigations at the engine test bed, high cycle-to-cycle variability was detected even for relatively stable peak power / full load operations of the engine, where knock onset severely limited the overall engine performance. In order to overcome the excessive computational cost of a direct chemical solution within a LES framework, the use of look-up tables for auto-ignition modelling perfectly fits with the strict mesh requirements of a LES simulation, with an acceptable approximation of the actual chemical kinetics. The model here presented is a totally stand-alone tool for autoignition analysis integrated with look-up table reading from detailed chemical kinetic schemes for gasoline.