Search Results

The incorporation of detailed chemistry models in internal combustion engine simulations is becoming mandatory as local, globally lean, low-temperature combustion strategies are setting the path towards a more efficient and environmentally sustainable use of energy resources in transportation. In this paper, we assessed the computational efficiency of a recently developed sparse analytical Jacobian chemistry solver, namely ‘SpeedCHEM’, that features both direct and Krylov-subspace solution methods for maximum efficiency for both small and large mechanism sizes. The code was coupled with a high-dimensional clustering algorithm for grouping homogeneous reactors into clusters with similar states and reactivities, to speed-up the chemical kinetics solution in multi-dimensional combustion simulations.

The effects of imaging system resolution and laser sheet thickness on the measurement of the Batchelor scale were investigated in a single-cylinder optical engine. The Batchelor scale was determined by fitting a model spectrum to the dissipation spectrum that was obtained from fuel tracer planar laser-induced fluorescence (PLIF) images of the in-cylinder scalar field. The imaging system resolution was quantified by measuring the step-response function; the scanning knife edge technique was used to measure the 10-90% clip width of the laser sheet. In these experiments, the spatial resolution varied from a native resolution of 32.0 μm to 137.4 μm, and the laser sheet thickness ranged from 108 μm to 707 μm. Thus, the overall resolution of the imaging system was made to vary by approximately a factor of four in the in-plane dimension and a factor of six in the out-of-plane dimension.

The tradeoff between NOx emissions and combustion instability in an engine operating in the dual-fuel Reactivity Controlled Compression Ignition (RCCI) combustion mode was investigated using a combination of engine experiments and detailed CFD modeling. Experiments were performed on a single cylinder version of a General Motors/Fiat JTD 1.9L four-cylinder diesel engine. Gasoline was injected far upstream of the intake valve using an air assisted injector and fuel vaporization system and diesel was injected directly into the cylinder using a common rail injector. The timing of the diesel injection was swept from −70° ATDC to −20° ATDC while the gasoline percentage was adjusted to hold the average combustion phasing (CA50) and load (IMEPg) constant at 0.5° ATDC and 7 bar, respectively. At each operating point the variation in IMEP, peak PRR, and CA50 was calculated from the measured cylinder pressure trace and NOx, CO, soot and UHC were recorded.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

Reactivity Controlled Compression Ignition (RCCI) is an engine combustion strategy that produces low NO and PM emissions with high thermal efficiency. Previous RCCI research has been investigated in single-cylinder heavy-duty engines. The current study investigates RCCI operation in a light-duty multi-cylinder engine at 3 operating points. These operating points were chosen to cover a range of conditions seen in the US EPA light-duty FTP test. The operating points were chosen by the Ad Hoc working group to simulate operation in the FTP test. The fueling strategy for the engine experiments consisted of in-cylinder fuel blending using port fuel-injection (PFI) of gasoline and early-cycle, direct-injection (DI) of diesel fuel. At these 3 points, the stock engine configuration is compared to operation with both the original equipment manufacturer (OEM) and custom-machined pistons designed for RCCI operation.

A method for automatically generating hex-dominant meshes for Computational Fluid Dynamics (CFD) applications is presented in this article. Two important regions of the mesh for any CFD simulation are the interior mesh and the boundary layer mesh. The interior mesh needs to be fine in the critical flow regions to ensure accurate solutions. The proposed method uses Bubble Mesh algorithm which packs bubbles inside the geometry to generate the mesh nodes. Algorithm was tested for sample flow problems and improvements were made to interior and boundary layer mesh generation methods. The interior mesh is generated using directionality and sizing control functions specified on the points of a 3D grid generated over the entire geometry. This offers a flexible control over mesh sizing and local mesh refinement. Boundary layer mesh is important to accurately model the physics of boundary layer near the geometry walls.

Prior HCCI optical engine experiments utilizing laser-induced fluorescence (LIF) measurements of stratified fuel-air mixtures have demonstrated the utility of probability density function (PDF) statistics for correlating mixture preparation with combustion. However, PDF statistics neglect all spatial details of in-cylinder fuel distribution. The current computational paper examines the effects of spatial fuel distribution on combustion using a novel combination of a 3-D CFD model with a 1-D linear-eddy model of turbulent mixing. In the simulations, the spatial coarseness of initial fuel distribution prior to the start of heat release is varied while keeping PDF statistics constant. Several cases are run, and as the initial mixture is made coarser, combustion phasing monotonically advances due to high local equivalence ratios that persist longer. The effect of turbulent mixing is more complex.

A commercial CFD package was used to develop a three-dimensional, fully turbulent model of the compressible flow across a complex-geometry venturi, such as those typically found in small engine carburetors. The results of the CFD simulations were used to understand the effect of the different obstacles in the flow on the overall discharge coefficient and the static pressure at the tip of the fuel tube. It was found that the obstacles located at the converging nozzle of the venturi do not cause significant pressure losses, while those obstacles that create wakes in the flow, such as the fuel tube and throttle plate, are responsible for most of the pressure losses. This result indicated that an overall discharge coefficient can be used to correct the mass flow rate, while a localized correction factor can be determined from three-dimensional CFD simulations in order to calculate the static pressure at locations of interest within the venturi.

This paper presents a theoretical analysis of the fuel and air flows within the idle circuit found in simple carburetors. The idle circuit is modeled numerically using a dynamic model that considers the resistances of the flow paths as well as the inertia of the fuel. The modeling methodology is flexible, in that the organization and techniques can be applied to any configuration and geometry. The numerical model calculates the fuel flow response of carburetor idle/transition circuits to pressure variations associated with air flow through the venturi and around the throttle plate. The model is implemented for a typical small engine carburetor and the nominal results are presented for this specific design.

Homogeneous or partially premixed charge compression ignition combustion is considered to be an attractive alternative to traditional internal combustion engine operation because of its extremely low levels of pollutant emissions. However, since it is difficult to control the start of combustion timing, direct injection of fuel into the combustion chamber is often used for combustion phasing control, as well as charge preparation. In this paper, numerical simulations of compression ignition processes using gasoline fuel directly injected using a low pressure, hollow cone injector are presented. The multi-dimensional CFD code, KIVA3V, that incorporates various advanced sub-models and is coupled with CHEMKIN for modeling detailed chemistry, was used for the study. Simulation results of the spray behavior at various injection conditions were validated with available experimental data.

The introduction of high-pressure injection systems in D.I. diesel engines has highlighted already known drawbacks of in-cylinder turbulence modeling. In particular, the well known equilibrium hypothesis is far from being valid even during the compression stroke and moreover during the spray injection and combustion processes when turbulence energy transfer between scales occurs under non-equilibrium conditions. The present paper focuses on modeling in-cylinder engine turbulent flows. Turbulence is accounted for by using the RNG k-ε model which is based on equilibrium turbulence assumptions. By using a modified version of the Kiva-3 code, different mathematically based corrections to the computed macro length scale are proposed in order to account for non-equilibrium effects. These new approaches are applied to a simulation of a recent generation HSDI Diesel engine at both full load and partial load conditions representative of the emission EUDC cycle.

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.

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.

This research investigates how the handling of mixing and heat transfer in a multi-zone kinetic solver affects the prediction of carbon monoxide and hydrocarbon emissions for simulations of HCCI engine combustion. A detailed kinetics multi-zone model is now more closely coordinated with the KIVA3V computational fluid dynamics code for simulation of the compression and expansion processes. The fluid mechanics is solved with high spatial and temporal resolution (40,000 cells). The chemistry is simulated with high temporal resolution, but low spatial resolution (20 computational zones). This paper presents comparison of simulation results using this enhanced multi-zone model to experimental data from an isooctane HCCI engine.

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.

A new search technique, called Non-Gradient Step-Controlled algorithm (NGSC), is presented. The NGSC was applied independently from pre-selected starting points and as a supplement to a Genetic Algorithm (GA) to optimize a HSDI diesel engine using split injection strategies. It is shown that the NGSC handles well the challenges of a complex response surface and factor high-dimensionality, which demonstrates its capability as an efficient and accurate tool to seek “local” convergence on complex surfaces. By directly tracking the change of a merit function, the NGSC places no requirement on response surface continuity / differentiability, and hence is more robust than gradient-dependent search techniques. The directional search mechanism takes factor interactions into consideration, and search step size control is adopted to facilitate search efficiency.

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.

The three-dimensional CFD codes KIVA-II and KIVA-3 have been used together to study the effects of intake generated in-cylinder flow structure on fuel-air mixing and combustion in a direct injected (DI) Diesel engine. In order to more accurately account for the effect of intake flow on in-cylinder processes, the KIVA-II code has been modified to allow for the use of data from other CFD codes as initial conditions. Simulation of the intake and compression strokes in a heavy-duty four-stroke DI Diesel engine has been carried out using KIVA-3. Flow quantities and thermodynamic field information were then mapped into a computational grid in KIVA-II for use in the study of mixing and combustion. A laminar and turbulent timescale combustion model, as well as advanced spray models, including wave breakup atomization, dynamic drop drag, and spray-wall interaction has been used in KIVA-II.

Fine-wire resistance thermometers were used to measure compression gas temperatures in a motoring (nonfiring) cycle CFR engine. Temperature versus crankangle curves were obtained for the compression and expansion strokes by means of tungsten wires ranging in diameter from 0.15–1.00 mils and at speeds from 600–1800 rpm. The results were compared with the infrared pyrometer; the peak temperature and peak crankangle lags were determined as a function of the wire diameter and engine speed. Attempts to evaluate the instantaneous energy balance around the wire resulted in a negative heat transfer coefficient, for which no current satisfactory explanation is available, although other observers have reported similar phenomena. The tungsten resistance thermometer is simple to build, easy to install, and requires no modification of the engine block for use during motoring. Thus, it is suitable for comparing the compression temperatures of different design engines.

A one-dimensional, non-equilibrium, compressible law of the wall model is proposed to increase the accuracy of heat transfer predictions from computational fluid dynamics (CFD) simulations of internal combustion engine flows on engineering grids. Our 1D model solves the transient turbulent Navier-Stokes equations for mass, momentum, energy and turbulence under the thin-layer assumption, using a finite-difference spatial scheme and a high-order implicit time integration method. A new algebraic eddy-viscosity closure, derived from the Han-Reitz equilibrium law of the wall, with enhanced Prandtl number sensitivity and compressibility effects, was developed for optimal performance. Several eddy viscosity sub-models were tested for turbulence closure, including the two-equation k-epsilon and k-omega, which gave insufficient performance.