Search Results

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.

The effects of residual gas mixing were studied in a single-cylinder, air-cooled utility engine using both external exhaust gas recirculation (EGR) and internal residual retention. EGR was introduced far upstream of the throttle to ensure proper mixing. Internal residual was changed by varying the length of the valve overlap period. EGR was measured in the intake system; the total in-cylinder diluent was directly measured using a skip-fire, cylinder dumping technique. A sweep of diluent fraction was performed at different engine speeds, engine loads, fuel mixture preparation systems, and ignition timings. An optimum level of diluent, where the combined hydrocarbon and NOx emissions were minimal, was found to exist for each operating condition. Higher levels of diluent, either through internal retention or external recirculation, caused the combined emissions to increase.

This computational study compares predictions and experimental results for the use of direct injected iso-octane fuel during the negative valve overlap (NVO) period to achieve HCCI combustion. The total fuel injection was altered in two ways. First the pre-DI percent, (the ratio of direct injected fuel during the NVO period “pre-DI” to the secondary fuel supplied at the intake manifold “PI”), was varied at a fixed pre-DI injection timing, Secondly the timing of the pre-DI injection was varied while all of the fuel was supplied during the NVO period. A multi-zone, two-dimensional CFD simulation with chemistry was performed using KIVA-3V release 2 implemented with the CHEMKIN solver. The simulations were performed during the NVO period only.

A single cylinder Yamaha research engine was operated with gasoline HCCI combustion using negative valve overlap (NVO). The injection strategy for this study involved using fuel injected directly into the cylinder during the NVO period (pre-DI) along with a secondary injection either in the intake port (PI) or directly into the cylinder (DI). The effects of timing of the pre-DI injection along with the percent of fuel injected during the pre-DI injection were studied in two sets of experiments using secondary PI and DI injections in separate experiments. Results have shown that by varying the pre-DI timing and pre-DI percent the main HCCI combustion timing can be influenced as a result of varied heat release during the negative valve overlap period along with hypothesized varied degrees of reformation of the pre-DI injected fuel. In addition to varying the main combustion timing the ISFC, emissions and combustion stability are all influenced by changes in pre-DI timing and percent.

Hydraulically intensified medium pressure common rail (MPCR) electronic fuel injection systems are an attractive concept for heavy-duty diesel engine applications. They offer excellent packaging flexibility and thorough engine management system integration. Two different concepts were evaluated in this study. They are different in how the pressure generation and injection events are related. One used a direct principle, where the high-pressure generation and injection events occur simultaneously producing a near square injection rate profile. Another concept was based on an indirect principle, where potential energy (pressure) is first stored inside a hydraulic accumulator, and then released during injection, as a subsequent event. A falling rate shape is typically produced in this case. A unit pump, where the hydraulic intensifier is separated from the injector by a high-pressure line, and a unit injector design are considered for both concepts.

A multi-objective genetic algorithm methodology was applied to a heavy-duty diesel engine at three different operating conditions of interest. Separate optimizations were performed over various fuel injection nozzle parameters, piston bowl geometries and swirl ratios (SR). Different beginning of injection (BOI) timings were considered in all optimizations. The objective of the optimizations was to find the best possible fuel economy, NOx, and soot emissions tradeoffs. The input parameter ranges were determined using design of experiment methodology. A non-dominated sorting genetic algorithm II (NSGA II) was used for the optimization. For the optimization of piston bowl geometry, an automated grid generator was used for efficient mesh generation with variable geometry parameters. The KIVA3V release 2 code with improved ERC sub-models was used. The characteristic time combustion (CTC) model was employed to improve computational efficiency.

The combustion and emissions of a diesel/natural gas dual-fuel engine are studied. Available engine experimental data demonstrates that the dual-fuel configuration provides a potential alternative to diesel engine operation for reducing emissions. The experiments are compared to multi-dimensional model results. The computer code used is based on the KIVA-3V code and consists of updated sub-models to simulate more accurately the fuel spray atomization, auto-ignition, combustion and emissions processes. The model results show that dual-fuel engine combustion and emissions are well predicted by the present multi-dimensional model. Significant reduction in NOx emissions is observed in both the experiments and simulations when natural gas is substituted for diesel fuel. The HC emissions are under predicted by numerical model as the natural gas substitution is increased.

The effect of injection rate shape on spray evolution and emission characteristics is investigated and a methodology for active control of fuel injection is proposed. Extensive validation of advanced vaporization and primary jet breakup models was performed with experimental data before studying the effects of systematic changes of injection rate shape. Excellent agreement with the experiments was obtained for liquid and vapor penetration lengths, over a broad range of gas densities and temperatures. Also the predicted flame lift-off lengths of reacting diesel fuel sprays were in good agreement with the experiments. After the validation of the models, well-defined rate shapes were used to study the effect of injection rate shape on liquid and vapor penetration, flame lift-off lengths and emission characteristics.

The recently developed KIVA-GA computer code was used in the current study to optimize the combustion chamber geometry of a heavy -duty diesel truck engine and a high-speed direct-injection (HSDI) small-bore diesel engine. KIVA-GA performs engine simulations within the framework of a genetic algorithm (GA) global optimization code. Design fitness was determined using a modified version of the KIVA-3V code, which calculates the spray, combustion, and emissions formation processes. The measure of design fitness includes NOx, unburned HC, and soot emissions, as well as fuel consumption. The simultaneous minimization of these factors was the ultimate goal. The KIVA-GA methodology was used to optimize the engine performance using nine input variables simultaneously. Three chamber geometry related variables were used along with six other variables, which were thought to have significant interaction with the chamber geometry.

Detailed chemical kinetics was implemented in the KIVA-3V multidimensional CFD code to study the combustion process in Homogeneous Charge Compression Ignition (HCCI) engines. The CHEMKIN code was implemented such that the chemistry and flow solutions were coupled. Detailed reaction mechanisms were used to simulate the fuel chemistry of ignition and combustion. Effects of turbulent mixing on the reaction rates were also considered. The model was validated using the experimental data from two modified heavy-duty diesel engines, including a Volvo engine and a Caterpillar engine operated at the HCCI mode. The results show that good levels of agreement were obtained using the present KIVA/CHEMKIN model for a wide range of engine conditions, including various fuels, injection systems, engine speeds, and EGR levels. Ignition timings were predicted well without the need to adjust any kinetic constants.

An improved discrete particle ignition kernel (DPIK) model and the G-equation combustion model have been developed and implemented in KIVA-3V. In the ignition model, the spark ignition kernel growth is tracked by Lagrangian markers and the spark discharge energy and flow turbulence effects on the ignition kernel growth are considered. The predicted ignition kernel size was compared with the available measurements and good agreement was obtained. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the level set method (G-equation) is used to track the mean turbulent flame propagation. It is shown that, by ignoring the detailed turbulent flame brush structure, fine numerical resolution is not needed, thus making the models suitable for use in multidimensional modeling of SI engine combustion.

The present study uses a numerical model to investigate the effects of flow turbulence on premixed iso-octane HCCI engine combustion. Different levels of in-cylinder turbulence are generated by using different piston geometries, namely a disc-shape versus a square-shape bowl. The numerical model is based on the KIVA code which is modified to use CHEMKIN as the chemistry solver. A detailed reaction mechanism is used to simulate the fuel chemistry. It is found that turbulence has significant effects on HCCI combustion. In the current engine setup, the main effect of turbulence is to affect the wall heat transfer, and hence to change the mixture temperature which, in turn, influences the ignition timing and combustion duration. The model also predicts that the combustion duration in the square bowl case is longer than that in the disc piston case which agrees with the measurements.

A new primary breakup model was developed and applied to simulate the diesel fuel spray and atomization process. The continuous liquid fuel jet was simulated by a discrete Lagrangian particle method, and the primary breakup of the jet was calculated using a new 1-D Eulerian method that provides the jet breakup time and drop size distribution. A set of correlations of the breakup characteristics, including the breakup time and drop size, were developed for a range of operating conditions. The correlations were then used in the KIVA code to predict the jet primary breakup. For drop secondary breakups, the Kelvin-Helmholtz/Rayleigh-Taylor hybrid model was employed. The new primary breakup model was first validated by comparison to experimental breakup length and jet liquid tip penetration lengths. Predictions of the new breakup model were also compared with experimental data and predictions of the standard breakup model.

A level set G-equation model has been developed to model the combustion process in spark ignition engines. The spark ignition process was modeled using an improved version of the Discrete Particle Ignition Kernel (DPIK) model. The two models were implemented into the KIVA-3V code to simulate SI engine combustion under both premixed and direct injection conditions. In the ignition model, the ignition kernel growth is tracked by Lagrangian markers, and spark discharge energy and flow turbulence effects on the kernel growth are considered. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the G-equation model is used to track the mean turbulent flame evolution. When combined with a characteristic time scale combustion model, the models were also used to simulate stratified combustion in DISI engines, where the triple flame structure must be considered.

An optimization study combining multidimensional CFD modeling and a global, evolutionary search technique known as the Genetic Algorithm has been carried out. The subject of this study was a 2-stroke, spark-ignited, direct-injection, single-cylinder research engine (SCRE). The goal of the study was to optimize the part load operating parameters of the engine in order to achieve the lowest possible emissions, improved fuel economy, and reduced wall heat transfer. Parameters subject to permutation in this study were the start-of-injection (SOI) timing, injection duration, spark timing, fuel injection angle, dwell between injections, and the percentage of fuel mass in the first injection pulse. The study was comprised of three cases. All simulations were for a part load, intermediate-speed condition representing a transition operating regime between stratified charge and homogeneous charge operation.

A radiation model based on the Discrete Ordinates Method (DOM) was incorporated into the KIVA3v multidimensional code to study the effects of soot and radiation on diesel engine performance at high load. A thermophoretic soot deposition model was implemented to predict soot concentrations in the near-wall region, which was found to affect radiative heat flux levels. Realistic, non-uniform combustion chamber wall surface temperature distributions were predicted using a finite-element-based heat conduction model for the engine metal components that was coupled with KIVA3v in an iterative scheme. The more accurate combustion chamber wall temperatures enhanced the accuracy of both the radiation and soot deposition models as well as the convective heat transfer model. For a basline case, (1500 rev/min, 100% load) it was found that radiation can account for as much as 30% of the total wall heat loss and that soot deposition in each cycle is less than 3% of the total in-cylinder soot.

An approach to accurately capture overall behavior in a system level model of DI Diesel HCCI engine operation is presented. The modeling methodology is an improvement over the previous effort [36], where a multi-zone model with detailed chemical kinetics was coupled with an engine cycle simulation code. This multi-zone technique was found to be inadequate in capturing the fuel spray dynamics and its impact on mixing. An improved methodology is presented in this paper that can be used to model fully and partially premixed charge compression ignition engines. A Computational Fluid Dynamics (CFD) driven model is used where the effects of fuel injection, spray evolution, evaporation, and turbulent mixing are considered. The modeling approach is based on the premise that once the initial spray dynamics are correctly captured, the overall engine predictions during the combustion process can be captured with good accuracy.

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.

The residual gas fraction was measured in an air-cooled single-cylinder utility engine by directly sampling the trapped cylinder charge during a programmed misfire. Tests were performed for a range of fuel mixture preparation systems, cam timings, ignition timings, engine speeds and engine loads. The residual fraction was found to be relatively insensitive to the fuel mixture preparation system, but was, to a moderate degree, sensitive to the ignition timing. The residual fraction was found to be strongly affected by the amount of valve overlap and engine speed. The effects of engine speed and ignition timing were, in part, due to the in-cylinder conditions at EVO, with lower temperatures favoring higher residual fractions. The data were compared to existing literature models, all of which were found to be lacking.

A numerical study of the effects of injection parameters and operating conditions for diesel-fuel HCCI operation is presented with consideration of Variable Geometry Sprays (VGS). Methods of mixture preparation are explored that overcome one of the major problems in HCCI engine operation with diesel fuel and conventional direct injection systems, i.e., fuel loss due to wall impingement and the resulting unburned fuel. Low pressure injection of hollow cone sprays into the cylinder of a production engine with the spray cone angle changing during the injection period were simulated using the multi-dimensional KIVA-3V CFD code with detailed chemistry. Variation of the starting and ending spray angles, injection timing, initial cylinder pressure and temperature, swirl intensity, and compression ratio were explored. As a simplified case of VGS, two-pulse, hollow-cone sprays were also simulated.