Search Results

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 modification of the computational fluid dynamics code KIVA-II is presented that allows computations to be made in complex engine geometries. An example application is given in which three versions of KIVA-II are run simultaneously. Each version considers a separate block of the computational domain, and the blocks exchange boundary condition information with each other at their common interfaces. The use of separate blocks permits the connectedness of the overall computational domain to change with time. The scavenging flow in the cylinder, transfer pipes (ports), and exhaust pipe of a ported two-stroke engine with a moving piston was modeled in this way. Results are presented for three engine designs that differ only in the angle of their boost ports. The calculated flow fields and the resulting fuel distributions are shown to be markedly different with the different geometries.

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.

The combustion mechanism in a Compression Ignition Homogeneous Charge (CIHC) engine was studied. Previous experiments done on a four-stroke CIHC engine were modeled using the KIVA-II code with modifications to the combustion, heat transfer, and crevice flow submodels. A laminar and turbulence characteristic time combustion model that has been used for spark-ignited engine studies was extended to allow predictions of ignition. The rate of conversion from one chemical species to another is modeled using a characteristic time which is the sum of a laminar (high temperature) chemistry time, an ignition (low temperature) chemistry time, and a turbulence mixing time. The ignition characteristic time was modeled using data from elementary initiation reactions and has the Arrhenius form. It was found to be possible to match all engine test cases reasonably well with one set of combustion model constants.

Diesel engine in-cylinder combustion processes have been studied using computational models with particular attention to spray development, vaporization, fuel/air mixture formation and combustion. A thermodynamic zero-dimensional cycle analysis program was used to determine initial conditions for the multidimensional calculations. A modified version of the time-dependent, three-dimensional computational fluid dynamics code KIVA-II was used for the computations, with a detailed treatment for the spray calculations and a simplified model for combustion. The calculations were used to obtain an understanding of the potential predictive capabilities of the models. It was found that there is a strong sensitivity of the results to numerical grid resolution. With proper grid resolution, the calculations were found to reproduce experimental data for non- vaporizing and vaporizing sprays. However, for vaporizing sprays with combustion, extremely fine grids are needed.

Combustion and flow were calculated in a spark-ignited two-stroke crankcase-scavenged engine using a laminar and turbulent characteristic-time combustion submodel in the three-dimensional KIVA code. Both premixed-charge and fuel-injected cases were examined. A multi-cylinder engine simulation program was used to specify initial and boundary conditions for the computation of the scavenging process. A sensitivity study was conducted using the premixed-charge engine data. The influence of different port boundary conditions on the scavenging process was examined. At high delivery ratios, the results were insensitive to variations in the scavenging flow or residual fraction details. In this case, good agreement was obtained with the experimental data using an existing combustion submodel, previously validated in a four-stroke engine study.

High-pressure diesel sprays were simulated with an Eulerian-Lagrangian Spray and Atomization (ELSA) model, based on a multidimensional engine computational fluid dynamics (CFD) code KIVA-3V. The atomization of the dense liquid core in the near-nozzle region was modeled with turbulent mixing of the diesel fuel with the ambient gas. Under the continuum assumption of a fuel-air mixture in this region, two transport equations were solved for the liquid mass fraction and liquid surface area density. At a certain downstream location where the spray became dilute, a switch from the Eulerian to the Lagrangian approach was made to benefit from the advantages of the conventional Lagrangian droplet models, such as droplet collision and turbulent dispersion modeling. The droplet size and velocity to be initialized at this switch were determined by the local CFD cell properties.

Engine experiments and multi-dimensional modeling were used to explore Reactivity Controlled Compression Ignition (RCCI) to realize highly-efficient combustion with near zero levels of NOx and PM. In-cylinder fuel blending using port-fuel-injection of a low reactivity fuel and optimized direct-injection of higher reactivity fuels was used to control combustion phasing and duration. In addition to injection and operating parameters, the study explored the effect of fuel properties by considering both gasoline-diesel dual-fuel operation, ethanol (E85)-diesel dual fuel operation, and a single fuel gasoline-gasoline+DTBP (di-tert butyl peroxide cetane improver). Remarkably, high gross indicated thermal efficiencies were achieved, reaching 59%, 56%, and 57% for E85-diesel, gasoline-diesel, and gasoline-gasoline+DTBP respectively.

A vaporization model for realistic multi-component fuel sprays is described. The equilibrium at the interface between liquid droplets and the surrounding gas is obtained based on the UNIFAC method, which considers non-ideal molecular interactions that can greatly enhance or suppress the vaporization of the components in the system compared to predictions from ideal mixing using Raoult's Law, especially for polar fuels. The present results using the UNIFAC method are shown to be able to capture the azeotropic behaviors of polar molecule blends, such as mixtures of benzene and ethanol, benzene and iso-propanol, and ethanol and water [1]. Predicted distillation curves of mixtures of ethanol and multi-component gasoline surrogates are compared to those from experiments, and the model gives good improvements on predictions of the distillation curves for initial ethanol volume fractions ranging from 0% to 100%.

In this work, linear, non-linear and a generalized renormalization group (RNG) two-equation RANS turbulence models of the k-epsilon form were compared for the prediction of turbulent compressible flows in diesel engines. The object-oriented, multidimensional parallel code FRESCO, developed at the University of Wisconsin, was used to test the alternative models versus the standard k-epsilon model. Test cases featured the academic backward facing step and the impinging gas jet in a quiescent chamber. Diesel engine flows featured high-pressure spray injection in a constant volume vessel from the Engine Combustion Network (ECN), as well as intake flows in a high-swirl diesel engine. For the engine intake flows, a model of the Sandia National Laboratories 1.9L light-duty single cylinder optical engine was used.

In the current, numerical study RCCI combustion and emission characteristics using various fuel strategies are investigated, including methanol, ethanol, n-butanol and gasoline as the low reactivity fuel, and diesel fuel as the high reactivity fuel. A reduced Primary Reference Fuel (PRF)-alcohol chemical kinetic mechanism was coupled with a computational fluid dynamic (CFD) code to predict RCCI combustion under various operating conditions. The results show that a higher quantity of diesel was required to maintain the same combustion phasing with alcohol-diesel fuel blends, and the combustion durations and pressure rise rates of methanol-diesel (MD) and ethanol-diesel (ED) cases were much shorter and higher than those of gasoline-diesel (GD) and n-butanol-diesel (nBD) cases. The simulations also investigated the sensitivities of the direct injection strategies, intake temperature and premixed fuel ratio on RCCI combustion phasing control.

Premixed compression ignition (PCI) strategies offer the potential for simultaneously low NOx and soot emissions and diesel-like efficiency. However, these strategies are generally confined to low loads due to difficulties controlling the combustion phasing and heat release rate. Recent experiments have demonstrated that dual-fuel reactivity-controlled compression ignition (RCCI) combustion can improve PCI combustion control and expand the PCI load range. Previous studies have explored RCCI operation using port-fuel injection (PFI) of gasoline and direct-injection (DI) of diesel fuel. In this study, experiments are performed using a light-duty, single-cylinder research engine to investigate RCCI combustion using a single fuel with the addition of a cetane improver 2-ethylhexyl nitrate (EHN). The fuel delivery strategy consists of port-fuel injection of E10 (i.e., 10% ethanol in gasoline) and direct-injection of E10 mixed with 3% EHN.

The role of the fluid motion in a diesel engine on mixing and combustion was investigated using the CFD code Kiva-3v. The study considered pre-mixed charge compression ignition (PCCI) combustion that is a hybrid combustion system characterized by early injection timings and high amounts of EGR dilution to delay the start and lower the temperature of combustion. The fuel oxidizer mixture is not homogeneous at the start of combustion and therefore requires further mixing for complete combustion. PCCI combustion systems are characterized by relatively high CO and UHC emissions. This work investigates attenuating CO emissions by enhancing mixing processes through non-uniform flowfield motions. The fluid motion was characterized by the amount of average angular rotation about the cylindrical axis (swirl ratio) and the amount of non-uniform motion imparted by the relative amounts of mass inducted through tangential and helical intake ports in a 0.5L HSDI diesel engine.

The present study attempted to model experimental results obtained on an optical engine at the Sandia National Laboratory. Measurements of in-cylinder unburned hydrocarbon (UHC) distributions were provided using advanced optical diagnostics on a near production type piston. Previous multidimensional modeling provided accurate pressure profiles and heat release rate (HRR) predictions. However, the experimental UHC distribution was not matched, and the model predicted UHC extending from the bowl into the squish region in the expansion stroke. To explore the causes of this discrepancy a parametric study was performed using a variety of initial conditions, boundary conditions and model constants to explore their effects on the UHC distribution. Of the initial conditions, the swirl ratio was found to have the biggest impact on the UHC distribution.

The present work proposes a methodology for diesel engine development using scaling laws and computational optimization with multi-dimensional CFD tools. A previously optimized 450cc HSDI diesel engine was down-scaled to 400cc size using recently developed scaling laws. The scaling laws were validated by comparing the performance of these two engines, including pressure, HRR, peak and averaged temperature, and pollutant emissions. A novel optimization methodology, which is able to simultaneously optimize multiple operating conditions, was proposed. The method is based on multi-objective genetic algorithms, and was coupled with the KIVA3V Release 2 code to further optimize the down-scaled diesel engine. An adaptive multi-grid chemistry model was used in the KIVA3V code to reduce the computational cost of the optimization. The computations were conducted using high-throughput computing with the CONDOR system.

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.

This paper concerns an innovative concept to control HCCI combustion in diesel-fuelled engines. It was named Homogenous Charge Progressive Combustion (HCPC) and operates on the split-cycle principle. In previous papers the feasibility of this combustion concept was shown for light-duty diesel engines. This paper illustrates a CFD study concerning a heavy-duty version of the HCPC engine. The engine displaces 13 liters and develops 700 kW indicated power at 2200 rpm with 49% maximum indicated efficiency and clean combustion.

Dual-fuel reactivity controlled compression ignition (RCCI) engine experiments were conducted with port fuel injection of isooctane and direct injection of n-heptane. The experiments were conducted at a nominal load of 4.75 bar IMEPg, with low isooctane equivalence ratios. Two sets of experiments explored the effects of direct injection timing with single and double injections, and multi-dimensional CFD modeling was used to explore mixture preparation and timing effects. The findings were that if fuel-liner impingement is to be avoided, double injections provide a 40% reduction in CO and HC emissions, resulting in a 1% increase in thermal efficiency. The second engine experiment showed that there is a linear relationship between reactivity (PRF number) and intake temperature. It was also found that if the premixed fuel fraction is above a certain limit, the high-temperature heat release (HTHR) can be manipulated by changing the global PRF number of the in-cylinder fuel blend.

The simulation of combustion chemistry in internal combustion engines is challenging due to the need to include detailed reaction mechanisms to describe the engine physics. Computational times needed for coupling full chemistry to CFD simulations are still too computationally demanding, even when distributed computer systems are exploited. For these reasons the present paper proposes a time scale separation approach for the integration of the chemistry differential equations and applies it in an engine CFD code. The time scale separation is achieved through the estimation of a characteristic time for each of the species and the introduction of a sampling timestep, wherein the chemistry is subcycled during the overall integration. This allows explicit integration of the system to be carried out, and the step size is governed by tolerance requirements.

In an attempt to increase efficiency and lower critical and highly regulated emissions (i.e., NOx, PM and CO2) many advanced combustion strategies have been investigated. Most of the current strategies fall into the category of low temperature combustion (LTC), which allow emissions mandates to be met in-cylinder along with anticipated reduction in cost and complexity. These strategies, such as homogeneous charge compression ignition (HCCI), premixed charge compression ignition (PCCI), partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI), use early injection timings, resulting in a highly lean charge with increased specific heat ratios to improve thermal efficiency and reduce PM emissions. Lower combustion temperatures also avoid the activation of NOx formation reactions.