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.

Time-resolved measurements were made of the gas composition at the exhaust port of a direct-injection two-stroke engine operating at 2000 rpm and an air-fuel ratio of 30:1. A high-speed sampling valve capable of 1.0 ms (12 CAD) time resolution was used to collect samples 1 cm downstream of the exhaust port of the engine. The time-resolved NOx, CO2 and CO concentrations decreased continuously during the scavenging process due to the dilution by short-circuited air. The hydrocarbon emissions, however, behaved significantly differently from the other species. At the time of exhaust port opening the concentration was low, it reached a maximum value by BDC, then decreased slightly in the latter part of the scavenging event. The dilution rates calculated for the hydrocarbon data gave negative values, indicating that there was a significant production of hydrocarbons during the gas exchange period.

A single cylinder, Cummins NH, direct-injection, diesel engine has been operated in order to evaluate the effects of aromatic content and aromatic structure on diesel engine particulates. Results from three fuels are shown. The first fuel, a low sulfur Chevron diesel fuel was used as a base fuel for comparison. The other fuels consisted of the base fuel and 10% by volume of 1-2-3-4 tetrahydronaphthalene (tetralin) a single-ring aromatic and naphthalene, a double-ring aromatic. The fuels were chosen to vary aromatic content and structure while minimizing differences in boiling points and cetane number. Measurements included exhaust particulates using a mini-dilution tunnel, exhaust emissions including THC, CO2, NO/NOx, O2, injection timing, two-color radiation, soluble organic fraction, and cylinder pressure. Particulate measurements were found to be sensitive to temperature and flow conditions in the mini-dilution tunnel and exhaust system.

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.

The effect of the ring pack storage mechanism on the hydrocarbon (HC) emissions from an air-cooled utility engine has been studied using a simplified ring pack model. Tests were performed for a range of engine load, two engine speeds, varied air-fuel ratio and with a fixed ignition timing using a homogeneous, pre-vaporized fuel mixture system. The integrated mass of HC leaving the crevices from the end of combustion (the crank angle that the cumulative burn fraction reached 90%) to exhaust valve closing was taken to represent the potential contribution of the ring pack to the overall HC emissions; post-oxidation in the cylinder will consume some of this mass. Time-resolved exhaust HC concentration measurements were also performed, and the instantaneous exhaust HC mass flow rate was determined using the measured exhaust and cylinder pressure.

The focus of this paper is to discuss the modeling and control of intake plenum pressure on the Powertrain Control Research Laboratory's (PCRL) Single-Cylinder Engine (SCE) transient test system using a patented device known as the Intake Air Simulator (IAS), which dynamically controls the intake plenum pressure, and, subsequently, the instantaneous airflow into the cylinder. The IAS exists as just one of many devices that the PCRL uses to control the dynamic boundary conditions of its SCE transient test system to make it “think” and operate as though it were part of a Multi-Cylinder Engine (MCE) test system. The model described in this paper will be used to design a second generation of this device that utilizes both continuously and discretely actuating valves working in parallel.

Regenerative testing methods can be used to allow the testing of hydraulic pumps and motors at significantly higher power and flow levels than that of the power supply used. This method can also increase the accuracy of system efficiency measurements. The increase in accuracy is realized because only the power added to compensate for the system losses needs to be measured with great accuracy. Typically, for the operation points of interest this will be a much smaller quantity than the overall power of the system. Knowing that the error in flow measurements is a function of the quantity measured, the benefit of measuring the losses becomes clear. An additional, distinct advantage of regenerative testing is that no dynamometer or load is needed. This results in a much simpler test setup. This paper documents the development of such a test program at the University of Wisconsin-Madison.

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.

Appropriate spray modeling in multidimensional simulations of diesel engines is well known to affect the overall accuracy of the results. More and more accurate models are being developed to deal with drop dynamics, breakup, collisions, and vaporization/multiphase processes; the latter ones being the most computationally demanding. In fact, in parallel calculations, the droplets occupy a physical region of the in-cylinder domain, which is generally very different than the topology-driven finite-volume mesh decomposition. This makes the CPU decomposition of the spray cloud severely uneven when many CPUs are employed, yielding poor parallel performance of the spray computation. Furthermore, mesh-independent models such as collision calculations require checking of each possible droplet pair, which leads to a practically intractable O(np2/2) computational cost, np being the total number of droplets in the spray cloud, and additional overhead for parallel communications.

Large-scale flows in internal combustion engines directly affect combustion duration and emissions production. These benefits are significant given increasingly stringent emissions and fuel economy requirements. Recent efforts by engine manufacturers to improve in-cylinder flows have focused on the design of specially shaped intake ports. Utility engine manufacturers are limited to simple intake port geometries to reduce the complexity of casting and cost of manufacturing. These constraints create unique flow physics in the engine cylinder in comparison to automotive engines. An experimental study of intake-generated flows was conducted in a four-stroke spark-ignition utility engine. Steady flow and in-cylinder flow measurements were made using three simple intake port geometries at three port orientations. Steady flow measurements were performed to characterize the swirl and tumble-generating capability of the intake ports.

Nonlinear observer theories are applied to the engine estimation problem in order to reconstruct engine states based on the measured engine variables, and dynamic mean torque production and cylinder-by-cylinder engine models. Engine cylinder and intake manifold pressures are two important factors in engine control and diagnostics. This paper discusses how to design nonlinear engine cylinder pressure and intake manifold pressure observers that have good robustness and estimation accuracy. Sliding mode theory in Variable Structure Systems (VSS) have shown good performance and been successfully applied to many nonlinear systems. Accordingly, sliding observers are selected for this study.

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.