Search Results

An endoscope-based image acquisition-and-processing camera system was used for diagnostics of pilot injection combustion in a single-cylinder heavy duty diesel engine. A study of the pilot injection or light load is of interest because the spray breakup, mixing and vaporization processes are less influenced by heat feedback from the flame than in full injection cases. This allows the spray process to be decoupled from the combustion process. The experimental cases were modeled using a version of the KIVA-II code that includes improvements in the turbulence, wall heat transfer, spray, ignition and combustion models. Pilot injections of three different amounts (10, 15 and 20% of the fuel injected at 75% load and 1600 RPM) at different start-of-injection timings were studied. The imaging system included an endoscope, an intensified CID camera, a frame grabber and the control circuitry.

Two-Color imaging optics were developed and used to observe soot emission processes in a modern heavy-duty diesel engine. The engine was equipped with a common rail, electronically-controlled, high-pressure fuel injection system that is capable of up to four injection pulses per engine cycle. The engine was instrumented with an endoscope system for optical access for the combustion visualization. Multidimensional combustion and soot modeling results were used for comparisons to enhance the understanding and interpretation of the experimental data. Good agreement between computed and measured cylinder pressures, heat release and soot and NOx emissions was achieved. In addition, good qualitative agreement was found between in-cylinder soot concentration (KL) and temperature fields obtained from the endoscope images and those obtained from the multidimensional modeling.

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.

Simultaneous two-component measurements of gas velocity and multi-dimensional numerical simulation are employed to characterize the evolution of the in-cylinder turbulent flow structure in a re-entrant bowl-in-piston engine under motored operation. The evolution of the mean flow field, turbulence energy, turbulent length scales, and the various terms contributing to the production of the turbulence energy are correlated and compared, with the objectives of clarifying the physical mechanisms and flow structures that dominate the turbulence production and of identifying the source of discrepancies between the measured and simulated turbulence fields. Additionally, the applicability of the linear turbulent stress modeling hypothesis employed in the k-ε model is assessed using the experimental mean flow gradients, turbulence energy, and length scales.

In-cylinder flow and combustion processes simulated with the standard k-ε turbulence model and with an alternative model-employing a non-linear, quadratic equation for the turbulent stresses-are contrasted for both motored and fired engine operation at two loads. For motored operation, the differences observed in the predictions of mean flow development are small and do not emerge until expansion. Larger differences are found in the spatial distribution and magnitude of turbulent kinetic energy. The non-linear model generally predicts lower energy levels and larger turbulent time scales. With fuel injection and combustion, significant differences in flow structure and in the spatial distribution of soot are predicted by the two models. The models also predict considerably different combustion efficiencies and NOx emissions.

In-cylinder fluid velocity is measured in an optically accessible, fired HSDI engine at idle. The velocity field is also calculated, including the full induction stroke, using multi-dimensional fluid dynamics and combustion simulation models. A detailed comparison between the measured and calculated velocities is performed to validate the computed results and to gain a physical understanding of the flow evolution. Motored measurements are also presented, to clarify the effects of the fuel injection process and combustion on the velocity field evolution. The calculated mean in-cylinder angular momentum (swirl ratio) and mean flow structures prior to injection agree well with the measurements. Modification of the mean flow by fuel injection and combustion is also well captured.

A study was performed that takes into account both turbulence and chemical kinetic effects in the numerical simulation of diesel engine combustion in order to better understand the importance of their respective roles at changing operating conditions. An approach was developed which combines the simplicity and low computational and storage requests of the laminar-and-turbulent characteristic-time model with a detailed combustion chemistry model based on well-known simplified mechanisms. Assuming appropriate simplifications such as steady state or equilibrium for most of the radicals and intermediate species, the kinetics of hydrocarbons can be described by means of three overall steps. This approach was integrated in the KIVA-II code. The concept was validated and applied to a single-cylinder, heavy-duty engine. The simulation covers a wide range of operating conditions.

Engine design is a time consuming process in which many costly experimental tests are usually conducted. With increasing prediction ability of engine simulation tools, engine design aided by CFD software is being given more attention by both industry and academia. It is also of much interest to be able to use design information gained from an existing engine design of one size in the design of engines of other sizes to reduce design time and costs. Therefore it is important to study size-scaling relationships for engines over wide range of operating conditions. This paper presents CFD studies conducted for two production diesel engines - a light-duty GM-Fiat engine (0.5L displacement) and a heavy-duty Caterpillar engine (2.5L displacement). Previously developed scaling arguments, including an equal spray penetration scaling model and an extended, equal flame lift-off length scaling model were employed to explore the parametric scaling connections between the two engines.

An analysis of the soot formation and oxidation process in a conventional direct-injection (DI) diesel flame was conducted using numerical simulations. An improved multi-step phenomenological soot model that includes particle inception, particle coagulation, surface growth and oxidation was used to describe the soot formation and oxidation process. The soot model has been implemented into the KIVA-3V code. Other model Improvements include a piston-ring crevice model, a KH/RT spray breakup model, a droplet wall impingement model, a wall-temperature heat transfer model, and the RNG k-ε turbulence model. The Shell model was used to simulate the ignition process, and a laminar-and-turbulent characteristic time combustion model was used for the post-ignition combustion process. Experimental data from a heavy-duty, Cummins N14, research DI diesel engine operated with conventional injection under low-load conditions were selected as a benchmark.

Combustion simulations utilizing the characteristic time combustion model have been performed for four DI diesel engines ranging in size from heavy-duty to large-bore designs. It has been found that the pre-factor to the turbulent characteristic time acts as a scaling parameter between the engines. This phenomenon is explained in terms of the non-equilibrium behavior of the turbulent time and length scales, as is encountered in the rapidly distorting, spray-induced flows of DI diesel engines. In fact, the equilibrium assumption between turbulence production and dissipation, which forms the basis for the employed k-ε-type turbulence models, does not hold in these situations. For such flows, the real turbulent dissipation time scale is locally proportional to the turbulent characteristic time scale which is determined by a typical eddy turnover time.

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.

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.

In this work computational and experimental approaches are combined to characterize in-cylinder flow structures and local flow field properties during operation of the Sandia 1.9L light-duty optical Diesel engine. A full computational model of the single-cylinder research engine was used that considers the complete intake and exhaust runners and plenums, as well as the adjustable throttling devices used in the experiments to obtain different swirl ratios. The in-cylinder flow predictions were validated against an extensive set of planar PIV measurements at different vertical locations in the combustion chamber for different swirl ratio configurations. Principal Component Analysis was used to characterize precession, tilting and eccentricity, and regional averages of the in-cylinder turbulence properties in the squish region and the piston bowl.

To model sprays from pressure-swirl atomizers, the connection between the injector and the downstream spray must be considered. A new model for pressure-swirl atomizers is presented which assumes little knowledge of the internal details of the injector, but instead uses available observations of external spray characteristics. First, a correlation for the exit velocity at the injector exit is used to define the liquid film thickness. Next, the film must be modeled as it becomes a thin, liquid sheet and breaks up, forming ligaments and droplets. A linearized instability analysis of the breakup of a viscous, liquid sheet is used as part of the spray boundary condition. The spray angle is estimated from spray photographs and patternator data. A mass averaged spray angle is calculated from the patternator data and used in some of the calculations.

Particle Image Velocimetry (PIV) was used to make gas velocity and turbulence measurements in a motored diesel engine. The experiments were conducted using a single-cylinder version of the Caterpillar 3406 production engine. One of the exhaust valves and the fuel injector port were used to provide optical access to the combustion chamber so that modifications to the engine geometry were minimal, and the results are representative of the actual engine. Measurements of gas velocity were made in a plane in the piston bowl using TiO2 seed particles. The light sheet necessary for PIV was formed by passing the beam from a Nd:YAG laser through the injector port and reflecting the beam off a conical mirror at the center of the piston. PIV data was difficult to obtain due to significant out-of-plane velocities. However, data was acquired at 25° and 15° before top dead center of compression at 750 rev/min.

In the Atomization regime, liquid jets breakup either within the nozzle or immediately upon entering the chamber gas and drops much smaller than the jet diameter are formed. The mechanism of Atomization, which is presently unknown, was investigated by the simultaneous use of two photographic techniques. The initial transient was observed with a 106 frames/s camera and the steady state by a technique similar to spark photography. The experiment range was: liquid pressure 500 to 2500 psia; five mixtures of water and glycerol to vary the liquid viscosity; air, nitrogen, helium, and xenon at up to 600 psia as chamber gases to separate gas pressure from gas density effects; and 14 nozzle designs. Not changed were the temperature (room value), the nozzle diameter (340 μ), and the surface tension (70 dyne/cm).

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

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.