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.

Reactivity Controlled Compression Ignition (RCCI) is an engine combustion strategy that utilizes in-cylinder fuel blending to produce low NOx and PM emissions while maintaining high thermal efficiency. The current study investigates RCCI and conventional diesel combustion (CDC) operation in a light-duty multi-cylinder engine over transient operating conditions using a high-bandwidth, transient capable engine test cell. Transient RCCI and CDC combustion and emissions results are compared over an up-speed change from 1,000 to 2,000 rev/min. and a down-speed change from 2,000 to 1,000 rev/min. at a constant 2.0 bar BMEP load. The engine experiments consisted of in-cylinder fuel blending with port fuel-injection (PFI) of gasoline and early-cycle, direct-injection (DI) of ultra-low sulfur diesel (ULSD) for the RCCI tests and the same ULSD for the CDC tests.

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.

The higher cylinder peak pressure and pressure rise rate of modern diesel and gasoline fueled engines tend to increase combustion noise while customers demand lower noise. The multiple degrees of freedom in engine control and calibration mean there is more scope to influence combustion noise but this must first be measured before it can be balanced with other attributes. An efficient means to realize this is to calculate combustion noise from the in-cylinder pressure measurements that are routinely acquired as part of the engine development process. This publication reviews the techniques required to ensure accurate and precise combustion noise measurements. First, the dynamic range must be maximized by using an analogue to digital converter with sufficient number of bits and selecting an appropriate range in the test equipment.

The rapid increase in the sophistication of vehicle automation demands development of evaluation protocols tuned to understanding driver-automation interaction. Driving simulators provide a safe and cost-efficient tool for studying driver-automation interaction, and this paper outlines general considerations for simulator-based evaluation protocols. Several challenges confront automation evaluation, including the limited utility of standard measures of driver performance (e.g., standard deviation of lane position), and the need to quantify underlying mental processes associated with situation awareness and trust. Implicitly or explicitly vehicle automation encourages drivers to disengage from driving and engage in other activities. Thus secondary tasks play an important role in both creating representative situations for automation use and misuse, as well as providing embedded measures of driver engagement.

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.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

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.

In this paper, we studied the accuracy of computational modeling of the ignition of a pilot injectionin the Sandia National Laboratories (SNL) light-duty optical engine facility, using the physical properties of a cetane/iso-cetane Diesel Primary Reference Fuel (DPRF) mixture and the reaction kinetics of a well-validated mechanism for primary reference fuels. Local fuel-air equivalence ratio measurements from fuel tracer based planar laser-induced fluorescence (PLIF) experiments were used to compare the mixture formation predictions with KIVA-ERC-based simulations. The effects of variations in injection mass from 1 mg to 4 mg, in-cylinder swirl ratio, and near-TDC temperatures on non-combusting mixture preparation were analyzed, to assess the accuracy of the model in capturing average jet behavior, despite its inability to model the non-negligible jet-by-jet variations seen in the experiments.

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.

The University of Wisconsin-Madison Snowmobile Team has designed and constructed a clean, quiet, high performance snowmobile for entry in the 2008 Society of Automotive Engineers' Clean Snowmobile Challenge. Built on a 2003 cross-country touring chassis, this machine features a 750 cc fuel-injected four-stroke engine equipped with a fuel sensor which allows operation ranging from regular gasoline to an 85% blend of ethanol and gasoline (E85). The engine has been customized with a Mototron control system which allows for full engine optimization using a range of fuels from E00 to E85. Utilizing a heated oxygen sensor and a 3-way catalyst customized for this engine by W.C. Heraeus-GmbH, this sled reduces NOx, HC and CO emissions by up to 89% to an average specific mass of 0.484, 0.154, 4.94 g/kW-hr respectively. Finally, the Mototron system also allowed Wisconsin to extract another 4 kW from the Weber 750cc engine; producing 45 kW and 65 Nm of torque.

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.

Large-eddy simulation (LES) is a useful approach for the simulation of turbulent flow and combustion processes in internal combustion engines. This study employs the ANSYS Forte CFD package and explores several key and fundamental components of LES, namely, the subgrid-scale (SGS) turbulence models, the numerical schemes used to discretize the transport equations, and the computational mesh. The SGS turbulence models considered include the classic Smagorinsky model and a dynamic structure model. Two numerical schemes for momentum convection, quasi-second-order upwind (QSOU) and central difference (CD), were evaluated. The effects of different computational mesh sizes controlled by both fixed mesh refinement and a solution-adaptive mesh-refinement approach were studied and compared. The LES models are evaluated and validated against several flow configurations that are critical to engine flows, in particular, to fuel injection processes.

In this paper, large eddy simulation (LES) coupled with two uncertainty quantification (UQ) methods, namely latin-hypercube sampling (LHS) and polynomial chaos expansion (PCE), have been used to quantify the effects of model parameters and spray boundary conditions on diesel and gasoline spray simulations. Evaporating, non-reacting spray data was used to compare penetration, mixture fraction and spray probability contour. Two different sets of four uncertain variables were used for diesel and gasoline sprays, respectively. UQ results showed good agreement between experiments and predictions. UQ statistics indicated that discharge coefficient has stronger impact on gasoline than diesel sprays, and spray cone angle is important for vapor penetration of both types of sprays. Additionally, examination of the gasoline spray characteristics showed that plume-to-plume interaction and nozzle dribble are important phenomena that need to be considered in high-fidelity gasoline spray simulations.

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.

The results of two derivations relating to the fuel economy of hybrid-electric vehicles (vehicles which employ both a heat engine and electric drive system) are presented and their use is illustrated through the examples of the University of Wisconsin and TRW Systems Group hybrid-electric vehicles. The method of mileage estimation employs a specific fuel-consumption versus torque-speed map for the heat engine under study and knowledge of the hybrid-vehicle dynamics and road-load power. The method is easily extended to estimation of emission reductions through use of specific-emission-production versus torque-speed maps and is applicable to hybrid vehicles with other than electrical energy-storage systems.