Search Results

Two different types of mesh used for diesel combustion with the KIVA-4 code are compared. One is a well established conventional KIVA-3 type polar mesh. The other is a non-polar mesh with uniform size throughout the piston bowl so as to reduce the number of cells and to improve the quality of the cell shapes around the cylinder axis which can contain many fuel droplets that affect prediction accuracy and the computational time. This mesh is specialized for the KIVA-4 code which employs an unstructured mesh. To prevent dramatic changes in spray penetration caused by the difference in cell size between the two types of mesh, a recently developed spray model which reduces mesh dependency of the droplet behavior has been implemented. For the ignition and combustion models, the Shell model and characteristic time combustion (CTC) model are employed.

Reactivity-controlled compression ignition (RCCI) has been shown to be capable of providing improved engine efficiencies coupled with the benefit of low emissions via in-cylinder fuel blending. Much of the previous body of work has studied the benefits of RCCI operation using high injection pressures (e.g., 500 bar or greater) with common rail injection (CRI) hardware. However, low-pressure fueling technology is capable of providing significant cost savings. Due to the broad market adoption of gasoline direct injection (GDI) fueling systems, a market-type prototype GDI injector was selected for this study. Single-cylinder light-duty engine experiments were undertaken to examine the performance and emissions characteristics of the RCCI combustion strategy with low-pressure GDI technology and compared against high injection pressure RCCI operation. Gasoline and diesel were used as the low-reactivity and high-reactivity fuels, respectively.

This paper investigates flow and combustion in a full-cycle simulation of a four-stroke, three-valve HCCI engine by visualizing the flow with pathlines. Pathlines trace massless particles in a transient flow field. In addition to visualization, pathlines are used here to trace the history, or evolution, of flow fields and species. In this study evolution is followed from the intake port through combustion. Pathline analysis follows packets of intake charge in time and space from induction through combustion. The local scalar fields traversed by the individual packets in terms of velocity magnitude, turbulence, species concentration and temperatures are extracted from the simulation results. The results show how the intake event establishes local chemical and thermal environments in-cylinder and how the species respond (chemically react) to the local field.

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.

A transport equation residual model incorporating refined G-equation and detailed chemical kinetics combustion models has been developed and implemented in the ERC KIVA-3V release2 code for Gasoline Direct Injection (GDI) engine simulations for better predictions of flame propagation. In the transport equation residual model a fictitious species concept is introduced to account for the residual gases in the cylinder, which have a great effect on the laminar flame speed. The residual gases include CO2, H2O and N2 remaining from the previous engine cycle or introduced using EGR. This pseudo species is described by a transport equation. The transport equation residual model differentiates between CO2 and H2O from the previous engine cycle or EGR and that which is from the combustion products of the current engine cycle.

The ability to make fully resolved turbulent scalar field measurements has been demonstrated in an internal combustion engine using one-dimensional fluorobenzene fluorescence measurements. Data were acquired during the intake stroke in a motored engine that had been modified such that each intake valve was fed independently, and one of the two intake streams was seeded with the fluorescent tracer. The scalar energy spectra displayed a significant inertial subrange that had a −5/3 wavenumber power dependence. The scalar dissipation spectra were found to extend in the high-wavenumber regime, to where the magnitude was more than two decades below the peak value, which indicates that for all practical purposes the measurements faithfully represent all of the scalar dissipation in the flow.

An investigation of high speed direct injection (DI) compression ignition (CI) engine combustion fueled with gasoline injected using a triple-pulse strategy in the low temperature combustion (LTC) regime is presented. This work aims to extend the operation ranges for a light-duty diesel engine, operating on gasoline, that have been identified in previous work via extended controllability of the injection process. The single-cylinder engine (SCE) was operated at full load (16 bar IMEP, 2500 rev/min) and computational simulations of the in-cylinder processes were performed using a multi-dimensional CFD code, KIVA-ERC-Chemkin, that features improved sub-models and the Chemkin library. The oxidation chemistry of the fuel was calculated using a reduced mechanism for primary reference fuel combustion chosen to match ignition characteristics of the gasoline fuel used for the SCE experiments.

Many recent studies have shown that the Reactivity Controlled Compression Ignition (RCCI) combustion strategy can achieve high efficiency with low emissions. However, it has also been revealed that RCCI combustion is difficult at high loads due to its premixed nature. To operate at moderate to high loads with gasoline/diesel dual fuel, high amounts of EGR or an ultra low compression ratio have shown to be required. Considering that both of these approaches inherently lower thermodynamic efficiency, in this study natural gas was utilized as a replacement for gasoline as the low-reactivity fuel. Due to the lower reactivity (i.e., higher octane number) of natural gas compared to gasoline, it was hypothesized to be a better fuel for RCCI combustion, in which a large reactivity gradient between the two fuels is beneficial in controlling the maximum pressure rise rate.

Although turbocharging can extend the high load limit of low temperature combustion (LTC) strategies such as reactivity controlled compression ignition (RCCI), the low exhaust enthalpy prevalent in these strategies necessitates the use of high exhaust pressures for improving turbocharger efficiency, causing high pumping losses and poor fuel economy. To mitigate these pumping losses, the divided exhaust period (DEP) concept is proposed. In this concept, the exhaust gas is directed to two separate manifolds: the blowdown manifold which is connected to the turbocharger and the scavenging manifold that bypasses the turbocharger. By separately actuating the exhaust valves using variable valve actuation, the exhaust flow is split between two manifolds, thereby reducing the overall engine backpressure and lowering pumping losses. In this paper, results from zero-dimensional and one-dimensional simulations of a multicylinder RCCI light-duty engine equipped with DEP are presented.

Fundamental simulations using DNS type procedures were used to investigate the ignition, combustion characteristics and the lift-off trends of a spatially evolving turbulent liquid fuel jet. In particular, the spatially evolving n-Heptane spray injected in a two-dimensional rectangular domain with an engine like environment was investigated. The computational results were compared to the experimental observations from an optical engine as reported in the literature. It was found that an initial fuel rich combustion downstream of the spray tip is followed by diffusion combustion. Investigations were also made to understand the effects of injection velocity, ambient temperature and the droplet radius on the lift-off length. For each of these parameters three different values in a given range were chosen. For both injection velocity and droplet radius, an increase resulted in a near linear increase in the lift-off length.

The contribution to the engine-out hydrocarbon (HC) emissions from fuel that escapes the main combustion event in piston ring crevices was estimated for an air-cooled, V-twin utility engine. The engine was run with a homogeneous pre-vaporized mixture system that avoids the presence of liquid films in the cylinder, and their resulting contribution to the HC emissions. A simplified ring pack gas flow model was used to estimate the ring pack contribution to HC emissions; the model was tested against the experimentally measured blowby. At high load conditions the model shows that the ring pack returns to the cylinder a mass of HC that exceeds that observed in the exhaust, and thus, is the dominant contributor to HC emissions. At light loads, however, the model predicts less HC mass returned from the ring pack than is observed in the exhaust. Time-resolved HC measurements were performed and used to assess the effect of combustion quality on HC emissions.

The role of mixture stratification on combustion rate has been investigated in a constant volume combustion vessel in which mixtures of different equivalence ratios can be added in a spatially and temporally controlled fashion. The experiments were performed in a regime of low fluid motion to avoid the complicating effects of turbulence generated by the injection of different masses of fluid. Different mixture combinations were investigated while maintaining a constant overall equivalence ratio and initial pressure. The results indicate that the highest combustion rate for an overall lean mixture is obtained when all of the fuel is contained in a stoichiometric mixture in the vicinity of the ignition source. This is the result of the high burning velocity of these mixtures, and the complete oxidation which releases the full chemical energy.

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.

In this paper, knock in a Ford single cylinder direct-injection spark-ignition (DISI) engine was modeled and investigated using the KIVA-3V code with a G-equation combustion model coupled with detailed chemical kinetics. The deflagrative turbulent flame propagation was described by the G-equation combustion model. A 22-species, 42-reaction iso-octane (iC8H18) mechanism was adopted to model the auto-ignition process of the gasoline/air/residual-gas mixture ahead of the flame front. The iso-octane mechanism was originally validated by ignition delay tests in a rapid compression machine. In this study, the mechanism was tested by comparing the simulated ignition delay time in a constant volume mesh with the values measured in a shock tube under different initial temperature, pressure and equivalence ratio conditions, and acceptable agreements were obtained.

The fuel film thickness resulting from diesel fuel spray impingement was measured in a chamber at conditions representative of early injection timings used for low temperature diesel combustion. The adhered fuel volume and the radial distribution of the film thickness are presented. Fuel was injected normal to the impingement surface at ambient temperatures of 353 K, 426 K and 500 K, with densities of 10 kg/m3 and 25 kg/m3. Two injectors, with nozzle diameters of 100 μm and 120 μm, were investigated. The results show that the fuel film volume was strongly affected by the ambient temperature, but was minimally affected by the ambient density. The peak fuel film thickness and the film radius were found to increase with decreased temperature. The fuel film was found to be circular in shape, with an inner region of nearly constant thickness. The major difference observed with temperature was a decrease in the radial extent of the film.

An ongoing goal of the Powertrain Control Research Laboratory (PCRL) at the University of Wisconsin-Madison has been to expand and improve the ability of the single cylinder internal combustion research engine to represent its multi-cylinder engine counterpart. To date, the PCRL single cylinder engine test system is able to replicate both the rotational dynamics (SAE #2004-01-0305) and intake manifold dynamics (SAE #2006-01-1074) of a multi cylinder engine using a single cylinder research engine. Another area of interest is the replication of multi-cylinder engine cold start emissions data with a single-cylinder engine test system. For this replication to occur, the single-cylinder engine must experience heat transfer to the engine coolant as if it were part of a multi-cylinder engine, in addition to the other multi-cylinder engine transient effects.

The lift-off length plays a significant role in spray combustion as it influences the air entrainment upstream of the lift-off location and hence the soot formation. Accurate prediction of lift-off length thus becomes a prerequisite for accurate soot prediction in lifted flames. In the present study, KIVA-3v coupled with CHEMKIN, as developed at the Engine Research Center (ERC), is used as the CFD model. Experimental data from the Sandia National Labs. is used for validating the model predictions of n-heptane lift-off lengths and soot formation details in a constant volume combustion chamber. It is seen that the model predictions, in terms of lift-off length and soot mass, agree well with the experimental results for low ambient density (14.8 kg/m3) cases with different EGR rates (21% O2 - 8% O2). However, for high density cases (30 kg/m3) with different EGR rates (15% O2 - 8% O2) disagreements were found.

Multi-zone CFD simulations with detailed kinetics were used to model iso-octane HCCI experiments performed on a single-cylinder research engine. The modeling goals were to validate the method (multi-zone combustion modeling) and the reaction mechanism (LLNL 857 species iso-octane) by comparing model results to detailed exhaust speciation data, which was obtained with gas chromatography. The model is compared to experiments run at 1200 RPM and 1.35 bar boost pressure over an equivalence ratio range from 0.08 to 0.28. Fuel was introduced far upstream to ensure fuel and air homogeneity prior to entering the 13.8:1 compression ratio, shallow-bowl combustion chamber of this 4-stroke engine. The CFD grid incorporated a very detailed representation of the crevices, including the top-land ring crevice and head-gasket crevice. The ring crevice is resolved all the way into the ring pocket volume. The detailed grid was required to capture regions where emission species are formed and retained.

Transient engine operations are modeled and simulated with a 1-D code (GT Power) using heat release and emission data computed by a 3-D CFD code (Kiva3). During each iteration step of a transient engine simulation, the 1-D code utilizes the 3-D data to interpolate the values for heat release and emissions. The 3-D CFD computations were performed for the compression and combustion stroke of strategically chosen engine operating points considering engine speed, torque and excess air. The 3-D inlet conditions were obtained from the 1-D code, which utilized 3-D heat release data from the previous 1-D unsteady computations. In most cases, only two different sets of 3-D input data are needed to interpolate the transient phase between two engine operating points. This keeps the computation time at a reasonable level. The results are demonstrated on the load response of a generator which is driven by a medium-speed diesel engine.

Fundamental simulations in a quiescent cell under adiabatic conditions were made to understand the effect of temperature, equivalence ratio and the components of the recirculated exhaust gas, viz., CO2 and H2O, on the combustion of n-Heptane. Simulations were made in single phase in which evaporated n-Heptane was uniformly distributed in the domain. Computations were made for two different temperatures and four different EGR levels. CO2 or H2O or N2was used as EGR. It was found that the initiation of the main combustion process was primarily determined by two competing factors, i.e., the amount of initial OH concentration in the domain and the specific heat of the mixture. Further, initial OH concentration can be controlled by the manipulating the ambient temperature in the domain, and the specific heat capacity of the mixture via the mixture composition. In addition to these, the pre combustion and the subsequent post combustion can also be controlled via the equivalence ratio.