Search Results

A novel 2-zone combustion chamber concept (patent pending) was developed using multi-dimensional modeling. At minimum volume, an axial projection in the piston divides the volume into distinct zones joined by a communication channel. The projection provides a means to control the mixture formation and combustion phasing within each zone. The novel combustion system was applied to reactivity controlled compression ignition (RCCI) combustion in both light-duty and heavy-duty diesel engines. Results from the study of an 8.8 bar BMEP, 2600 RPM operating condition are presented for the light-duty engine. The results from the heavy-duty engine are at an 18.1 bar BMEP, 1200 RPM operating condition. The effect of several major design features were investigated including the volume split between the inner and outer combustion chamber volumes, the clearance (squish) height, and the top ring land (crevice) volume.

A novel 2-zone combustion system was examined at medium load operation consistent with loads in the light duty vehicle drive cycle (7.6 bar BMEP and 2600 rev/min). Pressure rise rate and noise can limit the part of the engine map where pre-mixed combustion strategies such as HCCI or RCCI can be used. The present 2-zone pistons have an axial projection that divides the near TDC volume into two regions (inner and outer) joined by a narrow communication channel defined by the squish height. Dividing the near TDC volume provides a means to prepare two fuel-air mixtures with different ignition characteristics. Depending on the fuel injection timing, the reactivity of the inner or outer volume can be raised to provide an ignition source for the fuel-air mixture in the other, less reactive volume. Multi-dimensional CFD modeling was used to design the 2-zone piston geometry examined in this study.

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

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.

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.

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.

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.

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.

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.

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.

Various single and split injection schemes are studied to provide a better understanding of fuel distribution during cold starting in DI diesel engines. Improved spray-wall interaction, fuel film and multicomponent vaporization models are used to analyze the combustion processes. Better combustion characteristics are obtained for the split injection schemes than with a single injection. An analysis of the fuel impingement processes identifies the mechanisms involved in producing the differences in vaporization and combustion of the fuel. A greater amount of splashing occurred for the split injections compared to a single injection. This behavior is attributed to the decreased film thickness (less dissipation of impingement energy), the decreased impingement area (obtained by increasing the impingement Weber number), and most importantly, the reduced frequency of drop impingement.

A novel algebraic similarity model for subgrid scalar dissipation rate has been developed as part of the Large Eddy Simulation (LES) package KIVA3V-LES for diesel engine study. The model is proposed from an a priori study using Direct Numerical Simulation (DNS) of forced isotropic turbulence. In the a posteriori test, fully resolved turbulent passive scalar field measurements are used to validate the model in actual engine flows. For reason of the length limit by SAE and the specific interest in engine applications, only a prior test and a posteriori test in engine flows are included in this paper. A posteriori tests in isotropic cube flow, turbulent round jet and flame cases will be presented in separate papers. An engine LES simulation of multi consecutive cycles was performed in this study.

Diesel engine combustion is simulated using Large Eddy Simulation (LES) with a multi-mode combustion (MMC) model. The MMC model is based on the combination of chemical kinetics, chemical equilibrium, and quasi-steady flamelet calculations in different local combustion regimes. The local combustion regime is identified by two combustion indices based on the local temperature and the extent of mixture homogeneity. The LES turbulence model uses the dynamic structure model (DSM) for sub-grid stresses. A new spray model in the LES context is used, and the Reynolds-averaged Navier-Stokes (RANS) based wall model is retained with the LES derived scales. These models are incorporated in the KIVA3V-ERC-Release 2 code for engine combustion simulations. A wide range of diesel engine operating conditions were chosen to validate the combustion model.

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.

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.

Natural gas is a promising alternative fuel for internal combustion engines due to its rich reserves and low price, as well as good physical and chemical properties. Its low carbon structure and high octane number are beneficial for CO2 reduction and knock mitigation, respectively. Diesel and natural gas dual fuel combustion is a viable pathway to utilize natural gas in diesel engines. To achieve high efficiency and low emission combustion in a practical diesel engine over a wide range of operating conditions, understanding the performance responses to engine system parameter variations is needed. The controllability of two combustion strategies, diesel pilot ignition (DPI) and single injection reactivity controlled compression ignition (RCCI), were evaluated using the multi-dimension CFD simulation in this paper.