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 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 focus of the present study was to characterize Reactivity Controlled Compression Ignition (RCCI) using a single-fuel approach of gasoline and gasoline mixed with a commercially available cetane improver on a multi-cylinder engine. RCCI was achieved by port-injecting a certification grade 96 research octane gasoline and direct-injecting the same gasoline mixed with various levels of a cetane improver, 2-ethylhexyl nitrate (EHN). The EHN volume percentages investigated in the direct-injected fuel were 10, 5, and 2.5%. The combustion phasing controllability and emissions of the different fueling combinations were characterized at 2300 rpm and 4.2 bar brake mean effective pressure over a variety of parametric investigations including direct injection timing, premixed gasoline percentage, and intake temperature. Comparisons were made to gasoline/diesel RCCI operation on the same engine platform at nominally the same operating condition.

An experimental investigation of non-intrusive combustion sensing was performed using a tri-axial accelerometer mounted to the engine block of a small-bore high-speed 4-cylinder compression-ignition direct-injection (CIDI) engine. This study investigates potential techniques to extract combustion features from accelerometer signals to be used for cycle-to-cycle engine control. Selection of accelerometer location and vibration axis were performed by analyzing vibration signals for three different locations along the block for all three of the accelerometer axes. A magnitude squared coherence (MSC) statistical analysis was used to select the best location and axis. Based on previous work from the literature, the vibration signal filtering was optimized, and the filtered vibration signals were analyzed. It was found that the vibration signals correlate well with the second derivative of pressure during the initial stages of combustion.

This study investigates n-dodecane split injections of “Spray A” from the Engine Combustion Network (ECN) using two different turbulence treatments (RANS and LES) in conjunction with a Conditional Moment Closure combustion model (CMC). The two modeling approaches are first assessed in terms of vapor spray penetration evolutions of non-reacting split injections showing a clearly superior performance of the LES compared to RANS: while the former successfully reproduces the experimental results for both first and second injection events, the slipstream effect in the wake of the first injection jet is not accurately captured by RANS leading to an over-predicted spray tip penetration of the second pulse. In a second step, two reactive operating conditions with the same ambient density were investigated, namely one at a diesel-like condition (900K, 60bar) and one at a lower temperature (750K, 50bar).

Three-dimensional reactive computational fluid dynamics (CFD) plays a crucial role in IC engine development tasks complementing experimental efforts by providing improved understanding of the combustion process. A widely adopted combustion model in the engine community for (partially) premixed combustion is the G-Equation where the flame front is represented by an iso-level of an arbitrary scalar G. A convective-reactive equation for this iso-surface is solved, for which the turbulent flame speed ST must be provided. In this study, the commonly used and well-established Damköhler approach is compared to a novel correlation, derived from an algebraic closure for the scalar dissipation of reaction progress as proposed by Kolla et al. [1].

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.

An existing three-stage ignition delay model which has seen successful application to Primary Reference Fuels (PRFs) has been extended to six surrogate fuels which constitute potential candidates for future Homogeneous Charge Compression Ignition (HCCI) engines. The fuels include petroleum-derived and oxygenated components and can be divided into low, intermediate and high cetane number groups. A new methodology to obtain the model parameters is presented which relies jointly on simulation and experimental data: in a first step, constant volume adiabatic reactor simulations using chemical kinetic mechanisms are performed to generate ignition delays for a very wide range of conditions, namely variations in equivalence ratio, Exhaust Gas Recirculation (EGR), pressure and temperature.

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.

The interaction of turbulent premixed methane combustion with the surrounding flow field can be studied using optically accessible test rigs such as a rapid compression expansion machine (RCEM). The high flexibility offered by such a test rig allows its operation at various thermochemical conditions at ignition. However, limitations inherent to such test rigs due to the absence of an intake stroke do not allow turbulence production as found in IC-engines. Hence, means to introduce turbulence need to be implemented and the relevant turbulence quantities have to be identified in order to enable comparability with engine relevant conditions. A dedicated high-pressure direct injection of air at the beginning of the compression phase is considered as a measure to generate adjustable turbulence intensities at spark timing and during the early flame propagation.

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.

Modern Diesel engines have become ever more complex systems with many degrees of freedom. Simultaneously, with increasing computational power, simulations of engines have become more popular, and can be used to find the optimum set up of engine operation parameters which result in the desired point in the emission-efficiency trade off. With increasing number of engine operation parameter combinations, the number of calculations increase exponentially. Therefore, adequate models for combustion and emissions with limited calculation costs are required. For obvious reasons, the accuracy of the ignition timing is a key point for the following combustion and emission model quality. Furthermore, the combination of mixing and chemical processes during the ignition delay is very challenging to model in a fast way for a wide range of operation conditions.

Large Eddy Simulations (LES) provide instantaneous values indispensable to conduct statistical studies of relevant fluctuating quantities for diesel sprays. However, numerous realizations are generally necessary for LES to derive statistically averaged quantities necessary for validation of the numerical framework by means of measurements and for conducting sensitivity studies, leading to extremely high computational efforts. In this context, the aim of this work is to explore and validate alternatives to the simulation of 20-50 single realizations at considerably lower computational costs, by taking advantage of the axisymmetric geometry and the Quasi-Steady-State (QSS) condition of the near nozzle flow at a certain time after start-of-injection (SOI).

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.

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.

Pneumatic hybridization of internal combustion engines may prove to be a viable and cost-efficient alternative to electric hybridization. This paper evaluates the effects of pneumatic hybridization of various engine concepts using the criteria of fuel efficiency, driveability, emissions, and cost efficiency. The most promising engine concept is found to be the pneumatic hybridization combined with downsizing and supercharging spark-ignited engines. With this concept, a fuel consumption reduction of over 30% compared to a standard engine with the same rated power can be achieved. The poor driveability usually associated with heavily downsized and supercharged engines is completely overcome by injecting additional air during transients. The most important design issues for this new concept are discussed and several possible solutions are presented. Following these considerations, the first fully functional hybrid pneumatic engine was realized.

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.

Plug-In Hybrid Vehicles (PHEVs) represent the middle point between Hybrid Electric Vehicles (HEVs) and Electric Vehicles (EVs), thus combining benefits of the two architectures. PHEVs can achieve very high fuel economy while preserving full functionality of hybrids - long driving range, easy refueling, lower emissions etc. These advantages come at an expense of added complexity in terms of available fuel. The PHEV battery is recharged both though regenerative braking and directly by the grid thus adding extra dimension to the control problem. Along with the minimization of the fuel consumption, the amount of electricity taken from the power grid should be also considered, therefore the electricity generation mix and price become additional parameters that should be included in the cost function.

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.