Search Results

We propose a novel dual-fuel combustion model for simulating heavy-duty engines with the G-Equation. Dual-Fuel combustion strategies in such engines features direct injection of a high-reactivity fuel into a lean, premixed chamber which has a high resistance to autoignition. Distinct combustion modes are present: the DI fuel auto-ignites following chemical ignition delay after spray vaporization and mixing; a reactive front is formed on its surroundings; it develops into a well-structured turbulent flame, which propagates within the premixed charge. Either direct chemistry or the flame-propagation approach (G- Equation), taken alone, do not produce accurate results. The proposed Dual-Fuel model decides what regions of the combustion chamber should be simulated with either approach, according to the local flame state; and acts as a “kernel” model for the G- Equation model. Direct chemistry is run in the regions where a premixed front is not present.

A new analytical method to aid in the understanding of the organic carbon (OC) phase of particulate matter (PM) from advanced compression ignition (ACI) operating modes, is presented. The presence of NO2 and unburned fuel aromatics in ACI emissions, and the low exhaust temperatures that result from this low temperature combustion strategy, provide the right conditions for the formation of carboxylic acids and nitroaromatic compounds. These polar compounds contribute to OC in the PM and are not typically measured using nonpolar solvent extraction methods such as the soluble organic fraction (SOF) method. The new extraction and detection method employs capillary electrophoresis with electrospray ionization mass spectrometry (CE-ESI MS) and was specifically developed to determine polar organic compounds in the ACI PM emissions. The new method identified both nitrophenols and aromatic carboxylic acids in the ACI PM.

Surface film formation is an important phenomenon during spray impingement in a combustion chamber. The film that forms on the chamber walls and piston bowl produces soot post-combustion. While some droplets stick to the wall surface, others splash and interact with the gas present inside the combustion chamber. Accurate prediction of both the film thickness and splashed mass is crucial for surface film model development since it leads to a precise estimation of the amount of soot and other exhaust gases formed. This information could guide future studies aimed at a comprehensive understanding of the combustion process and might enable development of engines with reduced emissions. Dynamic structure Large Eddy Simulation (LES) turbulence model implemented for in-cylinder sprays [1] has shown to predict the flow structure of a spray more accurately than the Reynolds-averaged Navier-Stokes turbulence model.

In most large-eddy simulation (LES) applications to two-phase engine flows, the liquid-air interactions need to be accounted for as source terms in the respective governing equations. Accurate calculation of these source terms requires the relative velocity “seen” by liquid droplets as they move across the flow, which generally needs to be estimated using a turbulent dispersion model. Turbulent dispersion modeling in LES is very scarce in the literature. In most studies on engine spray flows, sub-grid scale (SGS) models for the turbulent dispersion still follow the same stochastic approach originally proposed for Reynolds-averaged Navier-Stokes (RANS). In this study, an SGS dispersion model is formulated in which the instantaneous gas velocity is decomposed into a deterministic part and a stochastic part. The deterministic part is reconstructed using the approximate deconvolution method (ADM), in which the large-scale flow can be readily calculated.

To meet stringent emissions regulations, filtration devices are often used in engine exhaust systems to reduce particulate mass (PM) and particulate number (PN). Diesel particulate filters (DPFs) are a well-established means of reducing PM from diesel engines to meet emissions regulations. New emissions regulations will most likely require a similar technology on gasoline engines with direct injection, gasoline particulate filters (GPFs). Due to differences in the exhaust and particulate characteristics, the design and operation of GPFs and DPFs differ. In a DPF filtration is dominated by the buildup of a soot cake. Whereas in a GPF, much of the soot is trapped inside the porous substrate, or filter wall, where deep-bed filtration is dominant. Thus, an accurate model describing the porous filtration properties of GPF substrates is desired. The pore filtration model (PFM) was developed to more accurately model the deep-bed filtration process that occurs in a GPF.

A swept-volume method of calculating the volume swept by the flame during each time step is developed and used to improve the calculation of fuel reaction rates. The improved reaction rates have been applied to the ignition model and coupled with the level set G-equation combustion model. In the ignition model, a single initial kernel is formed after which the kernel is convected by the gas flow and its growth rate is determined by the flame speed and thermal expansion due to the energy transfer from the electrical circuit. The predicted ignition kernel size was compared with the available experimental data and good agreements were achieved. Once the ignition kernel reaches a size when the fully turbulent flame is developed, the G-equation model is switched on to track the mean turbulent flame front propagation.

Modeling ignition kernel development in spark ignition engines is crucial to capturing the sources of cyclic variability, both with RANS and LES simulations. Appropriate kernel modeling must ensure that energy transfer from the electrodes to the gas phase has the correct timing, rate and locations, until the flame surface is large enough to be represented on the mesh by the G-Equation level-set method. However, in most kernel models, geometric details driving kernel growth are missing: either because it is described as Lagrangian particles, or because its development is simplified, i.e., down to multiple spherical flames. This paper covers the geometric aspects of kernel development, which makes up the core of a Triangulated Lagrangian Ignition Kernel model. One (or multiple, if it restrikes) spark channel is initialized as a one-dimensional Lagrangian particle thread.

Low temperature combustion (LTC) strategies present a means of reducing soot and oxides of nitrogen (NOx) emissions while simultaneously increasing efficiency relative to conventional combustion modes. By sufficiently premixing fuel and air before combustion, LTC strategies avoid high fuel-to-air equivalence ratios that lead to soot production. Dilution of the mixture lowers the combustion temperatures to reduce NOx production and offers thermodynamic advantages for improved efficiency. However, issues such as high heat release rates (HRRs), incomplete combustion, and difficulty in controlling the timing of combustion arise with low equivalence ratios and combustion temperatures. Ignition delay (the time until the start of combustion) is a way to quantify the time available for fuel and air to mix inside the cylinder before combustion. Previous studies have used ignition delay to explain trends seen in LTC such as combustion stability and HRRs.

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.

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.

Three time integration schemes and four finite volume interpolation schemes for the convection term in momentum equation were tested under turbulent planar gas jet and Sandia non-reacting vaporizing Spray-H cases. The three time integration schemes are the first-order Euler implicit scheme, the second-order backward scheme, and the second-order Crank-Nicolson scheme. The four spatial interpolation schemes are cubic central, linear central, upwind, and vanLeer schemes. Velocity magnitude contour, centerline and radial mean velocity and Reynolds stress profiles for the planar turbulent gas jet case, and fuel vapor contour and liquid and vapor penetrations for the Diesel spray case predicted by the different numerical schemes were compared. The sensitivity of the numerical schemes to mesh resolution was also investigated. The non-viscosity based dynamic structure subgrid model was used. The numerical tool used in this study was OpenFOAM.

Rapid vehicle powertrain development has become a technological breakthrough for the design and implementation of vehicles that meet and exceed the fuel efficiency, cost, and performance targets expected by today’s consumer. Recently, advances in large scale additive manufacturing have provided the means to bridge hardware-in-the-loop with preproduction mule chassis testing. This paper details a case study from Oak Ridge National Laboratory bridging the powertrain-in-the-loop development process with vehicle systems implementation using big area additive manufacturing (BAAM). For this case study, the use of a component-in-the-loop laboratory with math-based models is detailed for the design of a battery electric powertrain to be implemented in a printed prototype mule. The ability for BAAM to accelerate the mule development process via the concept of computer-aided design to part is explored.

This paper is part of a larger body of experimental and computational work devoted to studying the role of close-coupled post injections on soot reduction in a heavy-duty optical engine. It is a continuation of an earlier computational paper. The goals of the current work are to develop new CFD analysis tools and methods and apply them to gain a more in depth understanding of the different in-cylinder environments into which fuel from main- and post-injections are injected and to study how the in-cylinder flow, thermal and chemical fields are transformed between start of injection timings. The engine represented in this computational study is a single-cylinder, direct-injection, heavy-duty, low-swirl engine with optical components. It is based on the Cummins N14, has a cylindrical shaped piston bowl and an eight-hole injector that are both centered on the cylinder axis. The fuel used was n-heptane and the engine operating condition was light load at 1200 RPM.

Due to the new challenge of meeting number-based regulations for particulate matter (PM), a numerical and experimental study has been conducted to better understand particulate formation in engines fuelled with compressed natural gas. The study has been conducted on a Heavy-Duty, Euro VI, 4-cylinder, spark ignited engine, with multipoint sequential phased injection and stoichiometric combustion. For the experimental measurements two different instruments were used: a condensation particle counter (CPC) and a fast-response particle size spectrometer (DMS) the latter able also to provide a particle size distribution of the measured particles in the range from 5 to 1000 nm. Experimental measurements in both stationary and transient conditions were carried out. The data using the World Harmonized Transient Cycle (WHTC) were useful to detect which operating conditions lead to high numbers of particles. Then a further transient test was used for a more detailed and deeper analysis.

Reactivity Controlled Compression Ignition (RCCI) has been shown to be an attractive concept to achieve clean and high efficiency combustion. RCCI can be realized by applying two fuels with different reactivities, e.g., diesel and gasoline. This motivates the idea of using a single low reactivity fuel and direct injection (DI) of the same fuel blended with a small amount of cetane improver to achieve RCCI combustion. In the current study, numerical investigation was conducted to simulate RCCI and HCCI combustion and emissions with various fuels, including gasoline/diesel, iso-butanol/diesel and iso-butanol/iso-butanol+di-tert-butyl peroxide (DTBP) cetane improver. A reduced Primary Reference Fuel (PRF)-iso-butanol-DTBP mechanism was formulated and coupled with the KIVA computational fluid dynamic (CFD) code to predict the combustion and emissions of these fuels under different operating conditions in a heavy duty diesel engine.

A novel 2-zone combustion chamber concept (patent pending) was developed using multi-dimensional modeling. At minimum volume, an axial projection in the piston divides the volume into distinct zones joined by a communication channel. The projection provides a means to control the mixture formation and combustion phasing within each zone. The novel combustion system was applied to reactivity controlled compression ignition (RCCI) combustion in both light-duty and heavy-duty diesel engines. Results from the study of an 8.8 bar BMEP, 2600 RPM operating condition are presented for the light-duty engine. The results from the heavy-duty engine are at an 18.1 bar BMEP, 1200 RPM operating condition. The effect of several major design features were investigated including the volume split between the inner and outer combustion chamber volumes, the clearance (squish) height, and the top ring land (crevice) volume.

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.

A novel 2-zone combustion system was examined at medium load operation consistent with loads in the light duty vehicle drive cycle (7.6 bar BMEP and 2600 rev/min). Pressure rise rate and noise can limit the part of the engine map where pre-mixed combustion strategies such as HCCI or RCCI can be used. The present 2-zone pistons have an axial projection that divides the near TDC volume into two regions (inner and outer) joined by a narrow communication channel defined by the squish height. Dividing the near TDC volume provides a means to prepare two fuel-air mixtures with different ignition characteristics. Depending on the fuel injection timing, the reactivity of the inner or outer volume can be raised to provide an ignition source for the fuel-air mixture in the other, less reactive volume. Multi-dimensional CFD modeling was used to design the 2-zone piston geometry examined in this study.