Search Results

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.

In the present paper optimization tools are used to recommend low-emission engine combustion chamber designs, spray targeting and swirl ratio levels for a heavy-duty diesel engine operated at high-load. The study identifies aspects of the combustion and pollution formation that are affected by mixing processes, and offers guidance for better matching of the piston geometry with the spray plume geometry for enhanced mixing. By coupling a GA (genetic algorithm) with the KIVA-CFD code, and also by utilizing an automated grid generation technique, multi-objective optimizations with goals of low emissions and fuel economy were achieved. Three different multi-objective genetic algorithms including a Micro-Genetic Algorithm (μGA), a Nondominated Sorting Genetic Algorithm II (NSGA II) and an Adaptive Range Multi-Objective Genetic Algorithm (ARMOGA) were compared for conducting the optimization under the same conditions.

The KIVA 3V code has been applied to predict combustion chamber heat flux in an air-cooled utility engine. The KIVA heat flux predictions were compared with experimentally measured data in the same engine over a wide range of operating conditions. The measured data were found to be approximately two times larger than the predicted results, which is attributed to the omission of chemical heat release in the near-wall region for the heat transfer model applied. Modifying the model with a simple scaling factor provided a good comparison with the measured data for the full range of engine load, heat flux sensor location, air-fuel ratio and spark timings tested. The detailed spatially resolved results of the KIVA predictions were then used to develop a simplified model of the combustion chamber temporally integrated heat flux using an artificial neural network (ANN).

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.

A nonintrusive diagnostic technique has been developed by which dynamic axial piston-position and tilt-angle measurements have been made in a single-cylinder research engine. A laser beam, introduced into the combustion chamber through an optical port in the cylinder head, was reflected by a polished surface on the piston crown. Motion of the reflected beam, carrying with it information on piston position and piston tilt, was monitored by a set of receiving optics. Piston motion was studied as a function of both engine speed and cylinder pressure (i.e., piston loading.) Measured axial piston-position was found to deviate from the theoretical position calculated from the measured crank-shaft position owing to the effects of tilt and piston loading. Furthermore, evidence of piston veer (tilt of the piston in a plane parallel to the axis of the wrist pin) was observed, which had an effect on the accuracy of the axial piston-position measurement.

In-cylinder heat flux and temperature measurements were obtained in an air-cooled four-stroke utility engine for a range of air-fuel ratios. For these measurements, the magnitude of the integrated heat flux peaked at the stoichiometric air-fuel ratio, with an approximately linear decrease on either side of stoichiometric. Advancing the spark generally increased the magnitude of the integrated heat flux. Evaluation of the Brake Specific Integrated Heat Flux (BSIHF) mitigated these trends, and, the effects of changes in timing were eliminated for some operating conditions Examination of the BSIHF from the compression and expansion stroke showed behavior mimicking the full cycle BSIHF. However, the fraction of the total flux contributed by this portion of the cycle varied greatly from approximately 98% of the total to approximately 75% of the total.

A single cylinder CFR research engine has been run in HCCI combustion mode at the rich and the lean limits of the homogeneous charge operating range. To achieve a variation of the degree of charge stratification, two GDI injectors were installed: one was used for generating a homogeneous mixture in the intake system, and the other was mounted directly into the side of the combustion chamber. At the lean limit of the operating range, stratification showed a tremendous improvement in IMEP and emissions. At the rich limit, however, the stratification was limited by the high-pressure rise rate and high CO and NOx emissions. In this experiment the location of the DI injector was in such a position that the operating range that could be investigated was limited due to liquid fuel impingement onto the piston and liner.

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.

An ultrahigh injection pressure, common rail fuel injection system was designed, fabricated, and evaluated. The result was a system suitable for high-power density diesel engine applications. The main advantages of the concept are a very short injection duration capability, high injection pressure independent of engine speed, a simplified electronic control valve, and good packaging flexibility. Two prototype injectors were developed. Tests were performed on an injector flow bench and in a single cylinder research engine. The first prototype delivered 320 mm3 within 2.5 milliseconds with a 2600 bar peak injection pressure. A conventional minisac nozzle was used. The second prototype employed a specially designed pintle nozzle producing a near-zero cone angle liquid jet impinging on a 9-mm cylindrical target centered on the piston bowl crown (OSKA-S system). The second prototype had the capability to deliver 316mm3 in 0.97ms.

A computer model that is able to predict the occurrence of knock in spark ignition engines has been developed and implemented into the KIVA-3V code. Three major sub-models were used to simulate the overall process, namely the spark ignition model, combustion model, and end-gas auto-ignition models. The spark ignition and early flame development is modeled by a particle marker technique to locate the flame kernel. The characteristic-time combustion model is applied to simulate the propagation of the regular flame. The autoignition chemistry in the end-gas was modeled by a reduced chemical kinetics mechanism that is based on the Shell model. The present model was validated by simulating the experimental data in three different engines. The spark ignition and the combustion models were first validated by simulating a premixed Caterpillar engine that was converted to run on propane. Computed cylinder pressure agrees well with the experimental data.

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.

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.

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.

A spray model for pressure-swirl atomizers that is based on a linearized instability analysis of liquid sheets has been combined with an ignition and combustion model for stratified charge spark ignition engines. The ignition model has been advanced, such that the presence of dual spark plugs can now be accounted for. Independent validation of the spray model is achieved by investigating a pressure-swirl injector inside a pressure bomb containing air at ambient temperature. In a second step, the complete model is used to estimate the performance of a small marine DISI Two-Stroke engine operating in stratified charge mode. Simulation results and experimental data are compared for several different injection timings and the agreement is generally good such that there is confidence in the predictive quality of the combustion model. Finally the model is applied in a conceptual study to investigate possible benefits of split injections.

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 phenomenological nozzle flow model has been developed and implemented in both the FIRE and KIVA-II codes to simulate the effects of the nozzle geometry on fuel injection and spray processes. The model takes account of the nozzle passage inlet configuration, flow losses and cavitation, the injection pressure and combustion chamber conditions and provides initial conditions for multidimensional spray modeling. The discharge coefficient of the injector, the effective injection velocity and the initial drop or injected liquid ‘blob’ sizes are calculated dynamically during the entire injection event. The model was coupled with the wave breakup model to simulate experiments of non-vaporizing sprays under diesel conditions. Good agreement was obtained in liquid penetration, spray angle and drop size (Sauter Mean Diameter). The integrated model was also used to model combustion in a Cummins single-cylinder optical engine with good agreement.

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.

Numerical simulations are performed to investigate the fuel/air mixing preparation in a gasoline direct injection (GDI) engine. A two-valve OHV engine with wedge combustion chamber is investigated since automobiles equipped with this type of engine are readily available in the U.S. market. Modifying and retrofitting these engines for GDI operation could become a viable scenario for some engine manufactures. A pressure-swirl injector and wide spacing injection layout are adapted to enhance mixture preparation. The primary interest is on preparing the mixture with adequate equivalence ratio at the spark plug under a wide range of engine operating conditions. Two different engine operating conditions are investigated with respect to engine speed and load. A modified version of the KIVA-3V multi-dimensional CFD code is used. The modified code includes the Linearized Instability Sheet Atomization (LISA) model to simulate the development of the hollow cone spray.

A multi-objective genetic algorithm coupled with the KIVA3V release 2 code was used to optimize the piston bowl geometry, spray targeting, and swirl ratio levels of a high speed direct injected (HSDI) diesel engine for passenger cars. Three modes, which represent full-, mid-, and low-loads, were optimized separately. A non-dominated sorting genetic algorithm II (NSGA II) was used for the optimization. High throughput computing was conducted using the CONDOR software. An automated grid generator was used for efficient mesh generation with variable geometry parameters, including open and reentrant bowl designs. A series of new spray models featuring reduced mesh dependency were also integrated into the code. A characteristic-time combustion (CTC) model was used for the initial optimization for time savings. Model validation was performed by comparison with experiments for the baseline engine at full-, mid-, and low-load operating conditions.

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.