Search Results

The paper presents the development of a novel approach to the solution of detailed chemistry in internal combustion engine simulations, which relies on the analytical computation of the ordinary differential equations (ODE) system Jacobian matrix in sparse form. Arbitrary reaction behaviors in either Arrhenius, third-body or fall-off formulations can be considered, and thermodynamic gas-phase mixture properties are evaluated according to the well-established 7-coefficient JANAF polynomial form. The current work presents a full validation of the new chemistry solver when coupled to the KIVA-4 code, through modeling of a single cylinder Caterpillar 3401 heavy-duty engine, running in two-stage combustion mode.

An experimental study of luminous combustion in a modern diesel engine was performed to investigate the effect of injection parameters on NOX and soot formation via flame temperature and soot KL factor measurements. The two-color technique was applied to 2-D soot luminosity images and area-averaged soot radiation signals to obtain spatially and temporally resolved flame temperature and soot KL factor. The imaging system used for this study was based on a wide-angle endoscope that was mounted in the cylinder head and allowed different views of the combustion chamber. The experiments were carried out on a single-cylinder 2.4 liter D.I. diesel engine equipped with an electronically controlled common-rail injection system. Operating conditions were 1600 rpm and 75% load. The two-color results confirm that retarding the injection timing causes lower flame temperatures and NOX emissions but increased soot formation, independent of injection strategy.

The development of analytic models of diesel engine flow, combustion and subprocesses is described. The models are intended for use as design tools by industry for the prediction of engine performance and emissions to help reduce engine development time and costs. Part of the research program includes performing engine experiments to provide validation data for the models. The experiments are performed on a single-cylinder version of the Caterpillar 3406 engine that is equipped with state-of-the-art high pressure electronic fuel injection and emissions instrumentation. In-cylinder gas velocity and gas temperature measurements have also been made to characterize the flows in the engine.

An integrated numerical model has been developed for diesel engine computations based on the KIVA-II code. The model incorporates a modified RNG k-ε, turbulence model, a ‘wave’ breakup spray model, the Shell ignition model, the laminar-and-turbulent characteristic-time combustion model, a crevice flow model, a spray/wall impingement model that includes rebounding and breaking-up drops, and other improved submodels in the KIVA code. The model was validated and applied to model successfully different types of diesel engines under various operating conditions. These engines include a Caterpillar engine with different injection pressures at different injection timings, a small Tacom engine at different loads, and a Cummins engine modified by Sandia for optical experiments. Good levels of agreement in cylinder pressures and heat release rate data were obtained using the same computer model for all engine cases.

Due to concerns regarding pollutant and CO2 emissions, advanced combustion modes that can simultaneously reduce exhaust emissions and improve thermal efficiency have been widely investigated. The main characteristic of the new combustion strategies, such as HCCI and LTC, is that the formation of a homogenous mixture or a controllable stratified mixture is required prior to ignition. The major issue with these approaches is the lack of a direct method for the control of ignition timing and combustion rate, which can be only indirectly controlled using high EGR rates and/or lean mixtures. Homogeneous Charge Progressive Combustion (HCPC) is based on the split-cycle principle. Intake and compression phases are performed in a reciprocating external compressor, which drives the air into the combustor cylinder during the combustion process, through a transfer duct. A transfer valve is positioned between the compressor cylinder and the transfer duct.

An emissions and performance study was conducted to explore the effects of EGR and multiple injections on particulate, NOx, and BSFC. EGR is known to be effective at reducing NOx, but at high loads there is usually a large increase in particulate. Recent work has shown that multiple injections are effective at reducing particulate. Thus, it was of interest to examine the possibility of simultaneously reducing particulate and NOx with the combined use of EGR and multiple injections. The tests were conducted on a fully instrumented single cylinder version of the Caterpillar 3406 heavy duty truck engine. Tests were done at high load (75% of peak torque at 1600 RPM where EGR has been shown to produce unacceptable increases in particulate emissions. The fuel system used was an electronically controlled, common rail injector and supporting hardware. The fuel system was capable of up to four independent injections per cycle.

An experimental study has been completed which evaluated the effectiveness of using double, triple and rate shaped injections to simultaneously reduce particulate and NOx emissions. The experiments were done using a single cylinder version of a Caterpillar 3406 heavy duty D.I. diesel engine. The fuel system used was a common rail, electronically controlled injector that allowed flexibility in both the number and duration of injections per cycle. Injection timing was varied for each injection scheme to evaluate the particulate vs. NOx tradeoff and fuel consumption. Tests were done at 1600 rpm using engine load conditions of 25% and 75% of maximum torque. The results indicate that a double injection with a significantly long delay between injections reduced particulate by as much as a factor of three over a single injection at 75% load with no increase in NOx. Double injections with a smaller dwell gave less improvement in particulate and NOx at 75% load.

The three-dimensional computer code, KIVA, is being modified to include state-of-the-art submodels for diesel engine flow and combustion. Improved and/or new submodels which have already been implemented are: wall heat transfer with unsteadiness and compressibility, laminar-turbulent characteristic time combustion with unburned HC and Zeldo'vich NOx, and spray/wall impingement with rebounding and sliding drops. Progress on the implementation of improved spray drop drag and drop breakup models, the formulation and testing of a multistep kinetics ignition model and preliminary soot modeling results are described. In addition, the use of a block structured version of KIVA to model the intake flow process is described. A grid generation scheme has been developed for modeling realistic (complex) engine geometries, and initial computations have been made of intake flow in the manifold and combustion chamber of a two-intake-valve engine.

Progress on the development and validation of a CFD model for diesel engine combustion and flow is described. A modified version of the KIVA code is used for the computations, with improved submodels for liquid breakup, drop distortion and drag, spray/wall impingement with rebounding, sliding and breaking-up drops, wall heat transfer with unsteadiness and compressibility, multistep kinetics ignition and laminar-turbulent characteristic time combustion models, Zeldovich NOx formation, and soot formation with Nagle Strickland-Constable oxidation. The code also considers piston-cylinder-liner crevice flows and allows computations of the intake flow process in the realistic engine geometry with two moving intake valves. Significant progress has been made using a modified RNG k-ε turbulence model, and a multicomponent fuel vaporization model and a flamelet combustion model have been implemented.

It is well known that the accuracy of simulations of combustion processes in diesel and spark ignited (SI) engines depends on the initial conditions within the cylinder at intake valve closure (IVC). Residual gas affects the engine combustion processes through its influence on charge mass, temperature and dilution. In SI engines, there is little oxygen in the residual gas, and thus the dilution effect on flame propagation is more significant than in compression ignited (CI) engines. However, in CI engines, the ignition delay depends strongly on the in-cylinder gas temperature, which is proportional to the gas temperature at IVC. Furthermore, ignition delay is significantly affected by how much oxygen is present, which is also partly determined by the residual gas fraction. Therefore, it is of extreme importance to determine residual gas concentrations accurately.

A correction for the turbulence dissipation, based on non-equilibrium turbulence considerations from rapid distortion theory, has been derived and implemented in combination with the RNG k - ε model in a KIVA-based code. This model correction has been tested and compared with the standard RNG k - ε model for the compression and the combustion phase of two heavy duty DI diesel engines. The turbulence behavior in the compression phase shows clear improvements over the standard RNG k - ε model computations. In particular, the macro length scale is consistent with the corresponding time scale and with the turbulent kinetic energy over the entire compression phase. The combustion computations have been performed with the characteristic time combustion model. With this dissipation correction no additional adjustments of the turbulent characteristic time model constant were necessary in order to match experimental cylinder pressures and heat release rates of the two engines.

A two-dimensional transient Heat Conduction in Components code (HCC) was successfully set up and extensively used to calculate the temperature field existing in real engine combustion chambers. The Saul'yev method, an explicit, unconditionally stable finite difference method, was used in the code. Consideration of the gasket between the cylinder wall and head, and the air gap between the piston and liner were included in the code. The realistic piston bowl shape was modeled with a grid transformation and piston movement was considered. The HCC code was used to calculate the wall temperature of an Isuzu ceramic engine and a Caterpillar heavy-duty diesel engine. The code was combined with the KIVA-II code in an iterative loop, in which the KIVA-II code provided the instantaneous local heat flux on the combustion chamber surfaces, and the HCC code computed the time-averaged wall temperature distribution on the surfaces.

The three-dimensional KIVA code has been used to study the effects of injection pressure and split injections on diesel engine performance and soot and NOx emissions. The code has been updated with state-of-the-art submodels including: a wave breakup atomization model, drop drag with drop distortion, spray/wall interaction with sliding, rebounding, and breaking-up drops, multistep kinetics ignition and laminar-turbulent characteristic time combustion, wall heat transfer with unsteadiness and compressibility, Zeldovich NOx formation, and soot formation with Nagle Strickland-Constable oxidation. The computational results are compared with experimental data from a single-cylinder Caterpillar research engine equipped with a high-pressure, electronically-controlled fuel injection system, a full-dilution tunnel for soot measurements, and gaseous emissions instrumentation.

In previous work the KIVA-II code has been modified to model modem DI diesel engines and their emissions of particulate soot and oxides of nitrogen (NOx). This work presents results from a program to further validate the NOx emissions models against engine experiments with a well characterized modern engine. To facilitate a simplified comparison with experiments, a single cylinder research version of the Caterpillar 3406 heavy duty DI diesel engine was retrofitted to run as a naturally-aspirated, propane-fueled, spark-ignited engine. The retrofit includes installing a low compression ratio piston with bowl, adding a gas mixer, replacing the fuel injector assembly with a spark plug assembly and adding spark and fuel stoichiometry control hardware. Cylinder pressure and engine-out NOx emissions were measured for a range of speeds, exhaust gas residual (EGR) fractions, and spark timing settings.

In this paper, knock in a Ford single cylinder direct-injection spark-ignition (DISI) engine was modeled and investigated using the KIVA-3V code with a G-equation combustion model coupled with detailed chemical kinetics. The deflagrative turbulent flame propagation was described by the G-equation combustion model. A 22-species, 42-reaction iso-octane (iC8H18) mechanism was adopted to model the auto-ignition process of the gasoline/air/residual-gas mixture ahead of the flame front. The iso-octane mechanism was originally validated by ignition delay tests in a rapid compression machine. In this study, the mechanism was tested by comparing the simulated ignition delay time in a constant volume mesh with the values measured in a shock tube under different initial temperature, pressure and equivalence ratio conditions, and acceptable agreements were obtained.

An image acquisition-and-processing camera system was developed for in-cylinder diagnostics of a single-cylinder heavy duty diesel engine. The engine was equipped with an electronically-controlled common-rail fuel injection system that allowed both single and split (multiple) injections to be studied. The imaging system uses an endoscope to acquire luminous flame images from the combustion chamber and ensures minimum modification to the engine geometry. The system also includes an optical linkage, an image intensifier, a CID camera, a frame grabber, control circuitry and a computer. Experiments include both single and split injection cases at 90 MPa and 45 MPa injection pressures at 3/4 load and 1600 rev/min with simulated turbocharging. For the single injection at high injection pressure (90 MPa) the results show that the first luminous emissions from the ignition zone occur very close to the injector exit followed by rapid luminous flame spreading.

A three-dimensional computer code (KIVA) is being modified to include state-of-the-art submodels for diesel engine flow and combustion: spray atomization, drop breakup/coalescence, multi-component fuel vaporization, spray/wall interaction, ignition and combustion, wall heat transfer, unburned HC and NOx formation, soot and radiation and the intake flow process. Improved and/or new submodels which have been completed are: wall heat transfer with unsteadiness and compressibility, laminar-turbulent characteristic time combustion with unburned HC and Zeldo'vich NOx, and spray/wall impingement with rebounding and sliding drops.

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.

Multidimensional modeling is used to study air-fuel mixing in a direct-injection spark-ignition engine. Emphasis is placed on the effects of the start of fuel injection on gas/spray interactions, wall wetting, fuel vaporization rate and air-fuel ratio distributions in this paper. It was found that the in-cylinder gas/spray interactions vary with fuel injection timing which directly impacts spray characteristics such as tip penetration and spray/wall impingement and air-fuel mixing. It was also found that, compared with a non-spray case, the mixture temperature at the end of the compression stroke decreases substantially in spray cases due to in-cylinder fuel vaporization. The computed trapped-mass and total heat-gain from the cylinder walls during the induction and compression processes were also shown to be increased in spray cases.

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.