Search Results

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.

A multi-dimensional model was used to calculate interactions between spray drops and gas motions close to the nozzle in dense high-pressure sprays. The model also accounts for the phenomena of drop breakup, drop collision and coalescence, and the effect of drops, on the gas turbulence. The calculations used a new method to describe atomization (a boundary condition in current spray codes). The method assumes that atomization and drop breakup are indistinguishable processes within the dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops (‘blobs’) that have a size equal to the nozzle exit diameter. The injected ‘blobs’ breakup due to interaction with the gas as they penetrate, yielding a core region which contains relatively large drops. The computed core length agrees well with available measurements of core length in high-pressure sprays.

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.

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.

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.

A numerical simulation and optimization study was conducted for a medium speed direct injection diesel engine. The engine's operating characteristics were first matched to available experimental data to test the validity of the numerical model. The KIVA-3V ERC CFD code was then modified to allow independent spray events from two rows of nozzle holes. The angular alignment, nozzle hole size, and injection pressure of each set of nozzle holes were optimized using a micro-genetic algorithm. The design fitness criteria were based on a multi-variable merit function with inputs of emissions of soot, NOx, unburned hydrocarbons, and fuel consumption targets. Penalties to the merit function value were used to limit the maximum in-cylinder pressure and the burned gas temperature at exhaust valve opening. The optimization produced a 28.4% decrease in NOx and a 40% decrease in soot from the baseline case, while giving a 3.1% improvement in fuel economy.

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

Managing the injection rate profile is a powerful tool to control engine performance and emission levels. In particular, Common Rail (C.R.) injection systems allow an almost completely flexible fuel injection event in DI-diesel engines by permitting a free mapping of the start of injection, injection pressure, rate of injection and, in the near future, multiple injections. This research deals with the development of a network-based numerical tool for understanding operating condition limits of the Common Rail injector. The models simulate the electro-fluid-mechanical behavior of the injector accounting for cavitation in the nozzle holes. Validation against experiments has been performed. The model has been used to provide insight into the operating conditions of the injector and in order to highlight the application to injection system design.

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.

Multidimensional computations were carried out to simulate the in-cylinder fuel/air mixing process of a direct-injection spark-ignition engine using a modified version of the KIVA-3 code. A hollow cone spray was modeled using a Lagrangian stochastic approach with an empirical initial atomization treatment which is based on experimental data. Improved Spalding-type evaporation and drag models were used to calculate drop vaporization and drop dynamic drag. Spray/wall impingement hydrodynamics was accounted for by using a phenomenological model. Intake flows were computed using a simple approach in which a prescribed velocity profile is specified at the two intake valve openings. This allowed three intake flow patterns, namely, swirl, tumble and non-tumble, to be considered. It was shown that fuel vaporization was completed at the end of compression stroke with early injection timing under the chosen engine operating conditions.

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.

Fuel Injection System (FIS) research on injection pressure, timing control, and rate shaping, and studies on the modeling of injector nozzle flows and their effect on fuel spray characteristics are usually conducted separately. Only recently has the fuel injection and spraying process been studied as a complete system, i.e., including both the high-pressure fuel delivery and its effect on the nozzle flow characteristics, including nozzle cavitation. A methodology for coupling the fuel injection system and its effect on spray characteristics is presented here. The method is applied to an example case of a conventional pump-line-nozzle system. Mathematical models for characterizing the flows from the pump to the nozzle are formulated and solved using the Method of Characteristics and finite difference techniques. The nozzle internal flow is modeled using zero-dimensional flow models, in which the nozzle cavitation and its effect on the nozzle exit flow are accounted for.

Cycle-resolved analysis of velocity data obtained in the re-entrant bowl of a fired high-;speed, direct-injection diesel engine, demonstrates an unambiguous, approximately 100% increase in late-cycle turbulence levels over the levels measured during motored operation. Model predictions of the flow field, obtained employing RNG k-ε turbulence modeling in KIVA-3V, do not capture this increased turbulence. A combined experimental and computational approach is taken to identify the source of this turbulence. The results indicate that the dominant source of the increased turbulence is associated with the formation of an unstable distribution of mean angular momentum, characterized by a negative radial gradient. The importance of this source of flow turbulence has not previously been recognized for engine flows. The enhanced late-cycle turbulence is found to be very sensitive to the flow swirl level.