Search Results

Manufacturers of heavy-duty diesel engines are facing increasingly stringent, emission standards. These standards have motivated new research efforts towards improving the performance of diesel engines. The objective of the present program is to develop a comprehensive analytical model of the diesel combustion process that can be used to explore the influence of design changes. This will enable industry to predict the effect of these changes on engine performance and emissions. A major benefit of the successful implementation of such models is that engine development time and costs would be reduced through their use. The computer model is based on the three-dimensional KIVA-II code, with state-of-the-art submodels for spray atomization, drop breakup / coalescence, multi-component fuel vaporization, spray/wall interaction, ignition and combustion, wall heat transfer, unburned HC and NOx formation, and soot and radiation.

A comprehensive biodiesel combustion model is presented for use in multi-dimensional engine simulations. The model incorporates realistic physical properties in a vaporization model developed for multi-component fuel sprays and applies an improved mechanism for biodiesel combustion chemistry. Previously, a detailed mechanism for methyl decanoate and methyl-9-decenoate was reduced from 3299 species to 85 species to represent the components of biodiesel fuel. In this work, a second reduction was performed to further reduce the mechanism to 69 species. Steady and unsteady spray simulations confirmed that the model adequately reproduced liquid penetration observed in biodiesel spray experiments. Additionally, the new model was able to capture expected fuel composition effects with low-volatility components and fuel blend sprays penetrating further.

The objective of this investigation is to optimize a light-duty diesel engine in order to minimize soot, NOx, carbon monoxide (CO), unburned hydrocarbon (UHC) emissions and peak pressure rise rate (PPRR) while improving fuel economy in a low oxygen environment. Variables considered are the injection timings, fractional amount of fuel per injection, half included spray angle, swirl, and stepped-bowl piston geometry. The KIVA-CHEMKIN code, a multi-dimensional computational fluid dynamics (CFD) program with detailed chemistry is used and is coupled to a multi-objective genetic algorithm (MOGA) along with an automated grid generator. The stepped-piston bowl allows more options for spray targeting and improved charge preparation. Results show that optimal combinations of the above variables exist to simultaneously reduce emissions and fuel consumption. Details of the spray targeting were found to have a major impact on the combustion process.

In a recent study, quantitative measurements were presented of in-cylinder spatial distributions of mixture equivalence ratio in a single-cylinder light-duty optical diesel engine, operated with a non-reactive mixture at conditions similar to an early injection low-temperature combustion mode. In the experiments a planar laser-induced fluorescence (PLIF) methodology was used to obtain local mixture equivalence ratio values based on a diesel fuel surrogate (75% n-heptane, 25% iso-octane), with a small fraction of toluene as fluorescing tracer (0.5% by mass). Significant changes in the mixture's structure and composition at the walls were observed due to increased charge motion at high swirl and injection pressure levels. This suggested a non-negligible impact on wall heat transfer and, ultimately, on efficiency and engine-out emissions.

We propose a novel dual-fuel combustion model for simulating heavy-duty engines with the G-Equation. Dual-Fuel combustion strategies in such engines features direct injection of a high-reactivity fuel into a lean, premixed chamber which has a high resistance to autoignition. Distinct combustion modes are present: the DI fuel auto-ignites following chemical ignition delay after spray vaporization and mixing; a reactive front is formed on its surroundings; it develops into a well-structured turbulent flame, which propagates within the premixed charge. Either direct chemistry or the flame-propagation approach (G- Equation), taken alone, do not produce accurate results. The proposed Dual-Fuel model decides what regions of the combustion chamber should be simulated with either approach, according to the local flame state; and acts as a “kernel” model for the G- Equation model. Direct chemistry is run in the regions where a premixed front is not present.

The physical in-cylinder processes and ignition during cold starting have been studied using computational models, with particular attention to the influences of blowby, heat transfer during the compression stroke, spray development, vaporization and fuel/air mixture formation and ignition. Two different modeling approaches were used. A thermodynamic zero dimensional cycle analysis program in which the fuel injection effects were not modeled, was used to determine overall and gas exchange effects. The three-dimensional KIVA-II code was used to determine details of the closed cycle events, with modified atomization, blowby and spray/wall impingement models, and a simplified model for ignition. The calculations were used to obtain an understanding of the cold starting process and to identify practical methods for improving cold starting of direct injection diesel engines.

Real transportation fuels, such as gasoline and diesel, are mixtures of thousands of different hydrocarbons. For multidimensional engine applications, numerical simulations of combustion of real fuels with all of the hydrocarbon species included exceeds present computational capabilities. Consequently, surrogate fuel models are normally utilized. A good surrogate fuel model should approximate the essential physical and chemical properties of the real fuel. In this work, we present a novel methodology for the formulation of surrogate fuel models based on local optimization and sensitivity analysis technologies. Within the proposed approach, several important fuel properties are considered. Under the physical properties, we focus on volatility, density, lower heating value (LHV), and viscosity, while the chemical properties relate to the chemical composition, hydrogen to carbon (H/C) ratio, and ignition behavior. An error tolerance is assigned to each property for convergence checking.

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

Simultaneous measurements of the radial and the tangential components of velocity are obtained in a high-speed, direct-injection diesel engine typical of automotive applications. Results are presented for engine operation with fuel injection, but without combustion, for three different swirl ratios and four injection pressures. With the mean and fluctuating velocities, the r-θ plane shear stress and the mean flow gradients are obtained. Longitudinal and transverse length scales are also estimated via Taylor's hypothesis. The flow is shown to be sufficiently homogeneous and stationary to obtain meaningful length scale estimates. Concurrently, the flow and injection processes are simulated with KIVA-3V employing a RNG k-ε turbulence model. The measured turbulent kinetic energy k, r-θ plane mean strain rates ( 〈Srθ〉, 〈Srr〉, and 〈Sθθ〉 ), deviatoric turbulent stresses , and the r-θ plane turbulence production terms are compared directly to the simulated results.

The possibility to operate current diesel engines in dual-fuel mode with the addition of an alternative fuel is fundamental to accelerate the energy transition to achieve carbon neutrality. The simulation of the dual- fuel combustion process with 0D/1D combustion models is fundamental for the performance prediction, but still particularly challenging, due to chemical interactions of the mixture. The authors defined a novel data-driven workflow for the development of combustion reaction mechanisms and used it to generate a dual-fuel mechanism for Ammonia and Diesel Primary Reference Fuels (DPRF) suitable for efficient combustion simulations in heavy duty engines, with variable cetane number Diesel fuels. A baseline reaction mechanism was created by merging the detailed ammonia mechanism by Glarborg et al. with reaction pathways for n- hexadecane and 2,2,4,4,6,8,8-heptamethylnonane from a well-established multi-component fuel mechanism.

RNG k-ε closure turbulence dissipation equations are evaluated employing the CFD code KIVA-3V Release 2. The numerical evaluations start by considering simple jet flows, including incompressible air jets and compressible helium jets. The results show that the RNG closure turbulence model predicts lower jet tip penetration than the "standard" k-ε model, as well as being lower than experimental data. The reason is found to be that the turbulence kinetic energy is dissipated too slowly in the downstream region near the jet nozzle exit. In this case, the over-predicted R term in RNG model becomes a sink of dissipation in the ε-equation. As a second step, the RNG turbulence closure dissipation models are further tested in complex engine flows to compare against the measured evolution of turbulence kinetic energy, and an estimate of its dissipation rate, during both the compression and expansion processes.

Two-dimensional computations of premixed-charge engine combustion were made using the KIVA-II code. The purpose of the study was to assess the influence of heat transfer and turbulence model boundary conditions on engine combustion predictions. Combustion was modeled using a laminar- and turbulent-characteristic-time model. Flow through the piston-cylinder-ring crevice was accounted for using a phenomenological crevice-flow model. The predictions were compared to existing cylinder pressure and wall heat transfer experimental data under motoring and fired conditions, at two engine speeds. Two different wall heat transfer model formulations were considered. The first is the standard wall function method. The second is based on solutions to the one-dimensional unsteady energy equation, formulated such that the standard wall function method is recovered in the quasi-steady limit. Turbulence was modeled using the standard k-ε turbulence model equations.

Biodiesel-fueled engine simulations were performed using the KIVA3v-Release 2 code coupled with Chemkin-II for detailed chemistry. The model incorporates a reduced mechanism that was created from a methyl decanoate/methyl-9-decenoate mechanism developed at the Lawrence Livermore National Laboratory. A combination of Directed Relation Graph, chemical lumping, and limited reaction rate tuning was used to reduce the detailed mechanism from 3299 species and 10806 reactions to 77 species and 209 reactions. The mechanism was validated against its detailed counterpart and predicted accurate ignition delay times over a range of relevant operating conditions. The mechanism was then combined with the ERC PRF mechanism to include n-heptane as an additional fuel component. The biodiesel mechanism was applied in KIVA using a discrete multi-component model with accurate physical properties for the five common components of real biodiesel fuel.

Fuel spray atomization and breakup processes within a direct-injection spark-ignition (DISI) engine and outside the engine were modeled using a modified KIVA-3V code with improved spray models. The structures of the predicted sprays were qualitatively compared with planar images. The considered sprays were created by a prototype pressure-swirl injector and the planar images were obtained by laser sheet imaging in an optical DISI engine. In the out-of-engine case, the spray was injected into atmospheric air, and was modeled in a two dimensional bomb. In the engine case, the injection started from 270° ATDC, and full 3-D computations in the same engine were performed. In both cases, two liquid injection pressure conditions were applied, that is, 3.40 MPa and 6.12 MPa. The model gives good prediction of the tip penetration, and external spray shape, but the internal structure prediction has relatively lower accuracy, especially near the spray axis.

Multidimensional computations were made of spark-ignited premixed-charge combustion in a pancake-combustion-chamber engine with a centrally located spark plug and in two pent-roof-chamber engines, one with a central spark plug and the other with dual lateral spark plugs. A global combustion submodel was used that accounts for laminar kinetics and turbulent mixing effects. The predictions were compared with available measurements in the pancake-chamber engine over a range of loads, speeds, and equivalence ratios. In all cases the computed and measured cylinder pressures agreed well in trends and magnitudes (within 8%) for the entire duration of combustion. Fair agreements were also obtained between predicted and measured values of wall heat flux and emission index of nitric oxide. In the pent-roof-chamber engines the predicted maximum cylinder pressures also agreed well with measurements (within 12%) in cases with MBT (Minimum spark advance for Best Torque) or advanced spark timing.

An analytical mechanics model was employed to predict the delamination of several thermal-barrier-coated pistons that had been previously tested in a high-output, single-cylinder diesel engine. Some of the coatings delaminated during engine operation. Results are presented for two thicknesses of the same coating material, and for two similar coatings with different levels of stiffness. All the coating thermomechanical properties such as thermal conductivity, density, volumetric heat capacity, thickness, elastic modulus, coefficient of thermal expansion, Poisson ratio and toughness, were measured prior to engine testing. Previous measurements of the piston transient heat flux, based on fast-response surface temperature data, in the same engine were used as an input to calculate the multilayer wall temperature distribution. A theoretical methodology was employed to evaluate and predict the coating durability.

In the present study a reduced chemical reaction mechanism for biodiesel surrogate fuel was developed and validated for multi-dimensional engine combustion simulations. An existing detailed methyl butanoate mechanism that contained 264 species and 1219 reactions was chosen to represent the oxygenated portion of the fuel. The reduction process included flux analysis, ignition sensitivity analysis, and optimization of reaction rate constants under constant volume conditions. The current reduced mechanism consists of 41 species and 150 reactions and gives predictions in excellent agreement with those of the comprehensive mechanism. In order to validate the mechanism under biodiesel-fueled engine conditions, it was combined with another skeletal mechanism for n-heptane oxidation. This combined reaction mechanism can be used to adjust the energy content of the fuel, and account for diesel/biodiesel blend engine simulations.

A new ignition and combustion model has been developed and tested for use in premixed spark-ignition engines. The ignition model is referred to as the Discrete Particle Ignition Kernel (DPIK) model, and it uses Lagrangian markers to track the flame-front growth. The model includes the effects of electrode heat transfer on the early flame kernel growth process, and it is used in conjunction with a characteristic-time-scale combustion model once the ignition kernel has grown to a size where the effects of turbulence on the flame must be considered. A new term which accounts for the effect of air-fuel ratio, was added to the combustion model for modeling combustion in very lean and very rich mixtures. The flame kernel size predicted by the DPIK model was compared with measurements of Maly and Vogel. Furthermore, predictions of the electrode heat transfer were compared with data of Kravchik and Heywood. In both comparisons the model predictions were in good agreement with the experiments.

The development of a new Nitric Oxide (NOx) reaction mechanism has been conducted by adding species, including hydrogen cyanide (HCN) and the CH radical to a reduced chemistry diesel combustion model. The additional chemical reactions were added to the ERC's reduced 12-step NOx mechanism, which consists of N, NO, N2O, and NO2. The new NOx mechanism was implemented into the KIVA/ERC-CHEMKIN code and was found to be able to predict the experimentally observed trend that the amount of engine-out NOx decreases as engine load is increased, which is not reproduced by the current reduced NOx mechanism. HCN and CH were found to be species that bridge CxHy products and N radicals via the reaction CH+N2→HCN+N under high equivalence ratio conditions, and Zeldovich NO formation is suppressed by the formation of HCN, a species in the Fenimore NO formation pathway. The additional species and reactions were also found to influence the prediction of engine-out soot emissions.

Several diffusion combustion scaling models were experimentally tested in two geometrically similar single-cylinder diesel engines with a bore diameter ratio of 1.7. Assuming that the engines have the same in-cylinder thermodynamic conditions and equivalence ratio, the combustion models primarily change the fuel injection pressure and engine speed in order to attain similar performance and emissions. The models tested include an extended scaling model, which scales diffusion flame lift-off length and jet spray penetration; a simple scaling model, which only scales spray penetration at equal mean piston speed; and a same speed scaling model, which holds crankshaft rotational velocity constant while also scaling spray penetration. Successfully scaling diffusion combustion proved difficult to accomplish because of apparent differences that remained in the fuel-air mixing and heat transfer processes.