Search Results

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

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 multidimensional modeling, fuels have been represented predominantly by single components, such as octane for gasoline. Several bicomponent studies have been performed, but these are still limited in their ability to represent real fuels, which are blends of as many as 300 components. This study outlines a method by which the fuel composition is represented by a distribution function of the fuel molecular weight. This allows a much wider range of compositions to be modeled, and only requires including two additional “species” besides the fuel, namely the mean and second moment of the distribution. This approach has been previously presented but is applied here to multidimensional calculations. Results are presented for single component droplet vaporization for comparison with single component fuel predictions, as well as results for a multicomponent gasoline and a diesel droplet.

An integrated system model containing sub-models for diesel engine, emissions, and aftertreatment devices has been developed. The objective is to study engine-device and device-device interactions. The emissions sub-models used are for NOx and PM (particulate matter) prediction. The aftertreatment sub-models used include a diesel oxidation catalyst (DOC) and a diesel particulate filter (DPF). Controllers have also been developed to allow for transient simulations, active DPF regeneration, and prevention/control of runaway DPF regenerations. The integrated system-level model has been used to simulate DPF regeneration via exhaust fuel injection ahead of the DOC. In addition, the controller model can use intake throttling to assist in active DPF regeneration if needed. Regeneration studies have been done for both steady engine load and with load transients. High to low engine load transients are of particular interest because they can lead to runaway DPF regeneration.

A three-dimensional homogeneous equilibrium model (HEM) has been developed and implemented into an engine computational fluid dynamics (CFD) code KIVA-3V. The model was applied to simulate cavitating flow within injector nozzle passages. The effects of nozzle passage geometry and injection conditions on the development of cavitation zones and the nozzle discharge coefficient were investigated. Specifically, the effects of nozzle length (L/D ratio), nozzle inlet radius (R/D ratio) and K or KS factor (nozzle passage convergence) were simulated, and the effects of injection and chamber pressures, and time-varying injection pressure were also investigated. These effects are well captured by the nozzle flow model, and the predicted trends are consistent with those from experimental observations and theoretical analyses.

A general combustion model, in the context of large eddy simulations, was developed to simulate the full range of combustion in conventional diesel-type and HCCI-type diesels. The combustion model consisted of a Chemkin sub-model and an Extended Flamelet Time Scale (EFTS) sub-model. Specifically, Chemkin was used to simulate auto-ignition process. In the post-ignition phase, the combustion model was switched to EFTS. In the EFTS sub-model, combustion was assumed to be a combination of two elementary combustion modes: homogeneous combustion and flamelet combustion. The combustion index acted as a weighting factor blending the contributions from these two modes. The Chemkin sub-model neglected the subgrid scale turbulence-chemistry interactions whereas the EFTS model took them into account through a presumed PDF approach. The model was used to simulate an early injection mode of a Cummins DI diesel engine and a mode of a Caterpillar DI diesel engine.

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.

An integrated system model containing sub-models for a multi-cylinder diesel engine, NOx and soot(PM) emissions, diesel oxidation catalyst (DOC) and diesel particulate filter (DPF) has been developed to simulate the engine and aftertreatment systems at transient engine operating conditions. The objective of this work is two-fold; ensure correct implementation of the integrated system level model and apply the integrated model to understand the fuel economy trade-off for various DPF regeneration strategies. The current study focuses on a 1.9L turbocharged diesel engine and its exhaust system. The engine model was built in GT-Power and validated against experimental data at full-load conditions. The DPF model is calibrated for the current engine application by matching the clean DPF pressure drop for different mass flow rates. Load, boost pressure, speed and EGR controllers are tuned and linked with the current engine model.

Two advanced combustion models have been validated with the KIVA-3V Release 2 code in the context of two-stage combustion in a heavy duty diesel engine. The first model uses CHEMKIN to directly integrate chemistry in each computational cell. The second model accounts for flame propagation with the G-equation, and CHEMKIN predicts autoignition and handles chemistry ahead of and behind the flame front. A Damköhler number criterion was used in flame containing cells to characterize the local mixing status and determine whether heat release and species change should be a result of flame propagation or volumetric heat release. The purpose of this criterion is to make use of physical and chemical time scales to determine the most appropriate chemistry model, depending on the mixture composition and thermodynamic properties of the gas in each computational cell.

Neural networks are useful tools for optimization studies since they are very fast, so that while capturing the accuracy of multi-dimensional CFD calculations or experimental data, they can be run numerous times as required by many optimization techniques. This paper describes how a set of neural networks trained on a multi-dimensional CFD code to predict pressure, temperature, heat flux, torque and emissions, have been used by a genetic algorithm in combination with a hill-climbing type algorithm to optimize operating parameters of a diesel engine over the entire speed-torque map of the engine. The optimized parameters are mass of fuel injected per cycle, shape of the injection profile for dual split injection, start of injection, EGR level and boost pressure. These have been optimized for minimum emissions. Another set of neural networks have been trained to predict the optimized parameters, based on the speed-torque point of the engine.

Neural networks have been used for engine computations in the recent past. One reason for using neural networks is to capture the accuracy of multi-dimensional CFD calculations or experimental data while saving computational time, so that system simulations can be performed within a reasonable time frame. This paper describes three methods to improve upon neural network predictions. Improvement is demonstrated for in-cylinder pressure predictions in particular. The first method incorporates a physical combustion model within the transfer function of the neural network, so that the network predictions incorporate physical relationships as well as mathematical models to fit the data. The second method shows how partitioning the data into different regimes based on different physical processes, and training different networks for different regimes, improves the accuracy of predictions.

A modeling study has been conducted in order to characterize the heat transfer in an automotive diesel exhaust system. The exhaust system model, focusing on 2 exhaust pipes, has been created using a transient 1-D engine flow network simulation program. Model results are in excellent agreement with experimental data gathered before commencement of the modeling study. Predicted pipe exit stream temperatures are generally within one percent of experimental values. Sensitivity analysis of the model was the major focus of this study. Four separate variables were chosen for the sensitivity analysis. These being the external convective heat transfer coefficient, external emissivity, mass flow rate of exhaust gases, and amplitude of incoming pressure fluctuations. These variables were independently studied to determine their contribution to changes in exhaust gas stream temperature and system heat flux. There are two primary benefits obtained from conducting this analysis.

The present study investigated HCCI combustion in a heavy-duty diesel engine both experimentally and numerically. The engine was equipped with a hollow-cone pressure-swirl injector using gasoline direct injection. Characteristics of HCCI combustion were obtained by very early injection with a heated intake charge. Experimental results showed an increase in NOx emission and a decrease in UHC as the injection timing was retarded. It was also found that optimization can be achieved by controlling the intake temperature together with the start-of-injection timing. The experiments were modeled by using an engine CFD code with detailed chemistry. The CHEMKIN code was implemented into KIVA-3V such that the chemistry and flow solutions were coupled. The model predicted ignition timing, cylinder pressure, and heat release rates reasonably well. The NOx emissions were found to increase as the injection timing was retarded, in agreement with experimental results.

Artificial neural networks (ANN) have been recognized as universal approximators for nonlinear continuous functions and actively applied in engine research in recent years [1, 2, 3, 4, 5, 6, 7 and 8]. This paper describes the methodology and results of using the ANN to model a turbocharged DI diesel engine. The engine was simulated using the CFD code (KIVA-ERC) over a wide range of operating conditions, and numerical simulation results were used to train the ANN. An efficient data collection methodology using the Design of Experiments (DOE) techniques was developed to select the most characteristic engine operating conditions and hence the most informative data to train the ANN. This approach minimizes the time and cost of collecting training data from either computational or experimental resources. The trained ANN was then used to predict engine parameters such as cylinder pressure, cylinder temperature, NOx and soot emissions, and cylinder heat transfer.

DPF regenerations involve a trade-off between fuel economy and DPF durability. High temperature regenerations of DPFs have fewer fuel penalties but simultaneously tend to give higher substrate temperatures, which can reduce thermal reliability. In order to weaken the trade-off, the integrated system-level model [1,2,3,4] is used to conduct optimization studies and explore novel regeneration strategies for DPF aftertreatment devices. The integrated model developed in the Engine Research Center (ERC) includes sub-models for engines, emissions, aftertreatment devices and controllers. Based on the engine and regeneration fuel economy, multiple and single cycle regeneration tests are performed and analyzed. The optimal soot loadings to initiate and terminate regenerations are discussed. A pulsed regeneration strategy, which is characterized by injecting multiple pulses of fuel (upstream of a DOC) during regenerations, is investigated.

Large Eddy Simulation (LES) with a flamelet time scale combustion model is used to simulate diesel combustion. The flamelet time scale model uses a steady-state flamelet library for n-heptane indexed by mean mixture fraction, mixture fraction variance, and mean scalar dissipation rate. In the combustion model, reactions proceed towards the flamelet library solution at a time scale associated with the slowest reaction. This combination of a flamelet solution and a chemical time scale helps to account for unsteady mixing effects. The turbulent sub-grid stresses are simulated using a one-equation, non-viscosity LES model called the dynamic structure model. The model uses a tensor coefficient determined by the dynamic procedure and the subgrid kinetic energy. The model has been expanded to include scalar mixing and scalar dissipation. A new model for the conditional scalar dissipation has been developed to better predict local extinction.

To help account for fuel distribution during combustion in diesel engines, a fuel film model has been developed and implemented into the KIVA-II code [1]. Spray-wall interaction and spray-film interaction are also incorporated into the model. Modified wall functions for evaporating, wavy films are developed and tested. The model simulates thin fuel film flow on solid surfaces of arbitrary configuration. This is achieved by solving the continuity, momentum and energy equations for the two dimensional film that flows over a three dimensional surface. The major physical effects considered in the model include mass and momentum contributions to the film due to spray drop impingement, splashing effects, various shear forces, piston acceleration, dynamic pressure effects, and convective heat and mass transfer.

The three-dimensional CFD codes KIVA-II and KIVA-3 have been used together to study the effects of intake generated in-cylinder flow structure on fuel-air mixing and combustion in a direct injected (DI) Diesel engine. In order to more accurately account for the effect of intake flow on in-cylinder processes, the KIVA-II code has been modified to allow for the use of data from other CFD codes as initial conditions. Simulation of the intake and compression strokes in a heavy-duty four-stroke DI Diesel engine has been carried out using KIVA-3. Flow quantities and thermodynamic field information were then mapped into a computational grid in KIVA-II for use in the study of mixing and combustion. A laminar and turbulent timescale combustion model, as well as advanced spray models, including wave breakup atomization, dynamic drop drag, and spray-wall interaction has been used in KIVA-II.

A turbulent combustion model based on the coherent flamelet model was developed in this study and applied to diesel engines. The combustion was modeled in three distinct but overlapping phases: low temperature ignition kinetics using the Shell ignition model, high temperature premixed burn using a single step Arrhenius equation, and the flamelet based diffusion burn. Two criteria for transitions based on temperature, heat release rate, and the local Damköhler number were developed for the progression of combustion between each of these phases. The model was implemented into the computational computer code KIVA-II. Previous experiments on a Caterpillar model E 300, # 1Y0540 engine, a Tacom LABECO research engine, and a single cylinder version of a Cummins N14 production engine were used to validate the cylinder averaged predictions of the model.

The objective of the present work was to investigate a rapid method of obtaining the convection velocity of the bulk gas near the spark plug gap of a firing engine at the time of ignition. To accomplish this, a simple model was developed which utilized both the secondary current and voltage signals, from a conventional spark discharge. The model assumed the spark path was elongated in a rectangular U-shape by the flow. Based on experimentally measured electrical signals the mean convection velocity was computed. The convection velocity calculated by the model first needed calibration which was accomplished with a bench test that used a hot wire anemometer. The technique has a weak correlation at low velocities of 1-2 m/s, but correlates well at higher velocities up to 15 m/s.