Search Results

The Rate-Controlled Constrained-Equilibrium (RCCE) method is applied to the numerical solution of methane-air combustion. The RCCE method offers a reduction in computation time for complex chemically reacting systems because the rate equations for a small number of slowly evolving constraints need to be solved. The current work focuses on presenting both the principles of the RCCE method and its application to methane-air Homogeneous Charge Compression Ignition (HCCI) combustion. This work takes into consideration some of the previously unexplored numerical issues associated with solving the RCCE equation set. Application of the RCCE method is first demonstrated in constant and variable volume adiabatic environments and compared to the integration of the full set of kinetic rate equations for each species. Results presented here show a reduction in computational time.

A comprehensive system level modeling approach is used to understand the effects of the various physical actuators during diesel HCCI transients. Control concepts during transient operations are simulated using a set of actuators suitable for combustion control in diesel HCCI engines (intake valve actuation, injection timing, cooled EGR, intake boost pressure and droplet size). The impact of these actuating techniques on the overall engine performance is quantified by investigating the amount of actuation required, timing of actuation and the use of a combination of actuators. Combined actuation improved actuation space that can be used to phase combustion timing better and in extending the operating range. The results from transient simulations indicate that diesel HCCI operation would benefit from the combined actuation of intake valve closure, injection timing, boost and cooled EGR.

A procedure has been developed to build system level predictive models that incorporate physical laws as well as information derived from experimental data. In particular a soot model was developed, trained and tested using experimental data. It was seen that the model could fit available experimental data given sufficient training time. Future accuracy on data points not encountered during training was estimated and seen to be good. The approach relies on the physical phenomena predicted by an existing system level phenomenological soot model coupled with ‘weights’ which use experimental data to adjust the predicted physical sub-model parameters to fit the data. This approach has developed from attempts at incorporating physical phenomena into neural networks for predicting emissions. Model training uses neural network training concepts.

A fundamentally based quasi-dimensional NOx emission model for spark ignition direct injection (SIDI) gasoline engines was developed. The NOx model consists of a chemical mechanism and three sub-models. The classical extended Zeldovich mechanism and N₂O pathway for NOx formation mechanism were employed as the chemical mechanism in the model. A characteristic time model for the radical species H, O and OH was incorporated to account for non-equilibrium of radical species during combustion. A model of homogeneity which correlates fundamental dimensionless numbers and mixing time was developed to model the air-fuel mixing and inhomogeneity of the charge. Since temperature has a dominant effect on NOx emission, a flame temperature correlation was developed to model the flame temperature during the combustion for NOx calculation. Measured NOx emission data from a single-cylinder SIDI research engine at different operating conditions was used to validate the NOx model.

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.

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.

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

A kinetic carbon monoxide (CO) emission model is developed to simulate engine out CO emissions for conventional diesel combustion. The model also incorporates physics governing CO emissions for low temperature combustion (LTC). The emission model will be used in an integrated system level model to simulate the operation and interaction of conventional and low temperature diesel combustion with aftertreatment devices. The Integrated System Model consists of component models for the diesel engine, engine-out emissions (such as NOx and Particulate Matter), and aftertreatment devices (such as DOC and DPF). The addition of CO emissions model will enhance the capability of the Integrated System Model to predict major emission species, especially for low temperature combustion. In this work a CO emission model is developed based on a two-step global kinetic mechanism [8].

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.

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.

An integrated system based modeling approach has been developed to understand early Direct Injection (DI) Diesel Homogeneous Charge Compression Ignition (HCCI) process. GT-Power, a commercial one-dimensional (1-D) engine cycle code has been coupled with an external cylinder model which incorporates sub-models for fuel injection, vaporization, detailed chemistry calculations (Chemkin), heat transfer, energy conservation and species conservation. In order to improve the modeling accuracy, a multi-zone model has been implemented to account for temperature and fuel stratifications in the cylinder charge. The predictions from the coupled simulation have been compared with experimental data from a single cylinder Caterpillar truck engine modified for Diesel HCCI operation. A parametric study is conducted to examine the effect of combustion timing on four major control parameters. Overall the results show good agreement of the trends between the experiments and model predictions.

A two dimensional CFD model has been developed to study the mechanism of soot deposition on the porous wall surface in a Diesel Particulate Filter (DPF). The goal is to improve understanding of the soot deposition and its interaction with the hydrodynamic behaviour of the device. The KIVA3V CFD code and pre-processor were modified to simulate a single channel in a DPF. The code was extended to solve the conservation equations in porous media materials, to account for the sticking of particles on a porous surface and to evaluate the increasing resistance to the flow as the soot inside the trap accumulates. The code is already configured to track Lagrangian particles and these were modified to represent the soot particles in the flow. The code pre-processor was modified to allow definition of a double-symmetric geometry and to specify porous cells.

A Soot-NOx Trap (SNT) is a combinatorial aftertreatment device intended to decrease both particulate and NOx emissions simultaneously. A system-level Soot-NOx Trap model was developed by adding Lean NOx Trap kinetics to a 1D Diesel Particulate Filter model. The hybrid model was validated against each parent model for the limiting cases, then exercised to investigate the interacting redox behavior. Modulations in temperature and exhaust air-fuel ratio were investigated for their ability to facilitate particulate oxidation and NOx reduction in the trap.

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.

Two different types of mesh used for diesel combustion with the KIVA-4 code are compared. One is a well established conventional KIVA-3 type polar mesh. The other is a non-polar mesh with uniform size throughout the piston bowl so as to reduce the number of cells and to improve the quality of the cell shapes around the cylinder axis which can contain many fuel droplets that affect prediction accuracy and the computational time. This mesh is specialized for the KIVA-4 code which employs an unstructured mesh. To prevent dramatic changes in spray penetration caused by the difference in cell size between the two types of mesh, a recently developed spray model which reduces mesh dependency of the droplet behavior has been implemented. For the ignition and combustion models, the Shell model and characteristic time combustion (CTC) model are employed.

In this work the modeling aspects of fuel vaporization are studied. To start with, the effects of vaporization model on engine simulations are studied. This is done by using two different fuel surrogates. Next a set of non-reacting spray simulations were performed under different ambient and operating conditions and for two different fuels. This was done for spray model validation and to look at the effect of vaporization model on liquid penetration length. Following an observed discrepancy in one of the spray cases, effect of ambient temperature on liquid length, two sensitivity analyses were performed. These analyses take into account the effects of each spray-sub model on vaporization and effects of spray breakup constants on liquid penetration. Using the results from the sensitivity analyses and linearized stability theory an empirical correction factor was developed to correct the spray behavior at low ambient temperatures.

The effects of in-cylinder water injection on a direct injection (DI) Diesel engine were studied using a computational fluid dynamics (CFD) program based on the Kiva-3v code. The spray model is validated against experimental bomb data with good agreement for vapor penetration as a function of time. It was found that liquid penetration increased approximately 35% with 23% of the fuel volume replaced by water, due mostly to the increase in latent heat of vaporization. Engine calculations were compared to experimental results and showed very good agreement with pressure, ignition delay and fuel consumption. Trends for emissions were accurately predicted for both 44% and 86% load conditions. Engine simulations showed that the vaporization of liquid water as well as a local increase in specific heat of the gas around the flame resulted in lower Nitrogen Oxide emissions (NOx) and soot formation rates.

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.