Search Results

This work for the Coordinating Research Council (CRC) explores dependencies on the opportunity for fuel to impinge on internal engine surfaces (i.e., fuel–wall impingement) as a function of fuel properties and engine operating conditions and correlates these data with measurements of stochastic preignition (SPI) propensity. SPI rates are directly coupled with laser–induced florescence measurements of dye-doped fuel dilution measurements of the engine lubricant, which provides a surrogate for fuel–wall impingement. Literature suggests that SPI may have several dependencies, one being fuel–wall impingement. However, it remains unknown if fuel-wall impingement is a fundamental predictor and source of SPI or is simply a causational factor of SPI. In this study, these relationships on SPI and fuel-wall impingement are explored using 4 fuels at 8 operating conditions per fuel, for 32 total test points.

The possibility to operate current diesel engines in dual-fuel mode with the addition of an alternative fuel is fundamental to accelerate the energy transition to achieve carbon neutrality. The simulation of the dual- fuel combustion process with 0D/1D combustion models is fundamental for the performance prediction, but still particularly challenging, due to chemical interactions of the mixture. The authors defined a novel data-driven workflow for the development of combustion reaction mechanisms and used it to generate a dual-fuel mechanism for Ammonia and Diesel Primary Reference Fuels (DPRF) suitable for efficient combustion simulations in heavy duty engines, with variable cetane number Diesel fuels. A baseline reaction mechanism was created by merging the detailed ammonia mechanism by Glarborg et al. with reaction pathways for n- hexadecane and 2,2,4,4,6,8,8-heptamethylnonane from a well-established multi-component fuel mechanism.

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.

G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device.

This work explores the effect of the particle matter index (PMI) and aromatic content on fuel wall impingement associated with stochastic pre-ignition (SPI). Statically significant measurements of SPI rates are directly coupled with laser induced florescence (LIF) measurements of fuel dilution from spray-linear impingement. Literature suggests that PMI is could be correlated with the number of SPI events, but the root cause(s) of PMI and SPI are directly causational or are a predicator of SPI. Three fuels have been used in this study with 3 different PMI and two different aromatic contents. The fuels are direct injected at two different injection timings, an earlier injection timing which initially targets the piston crown, 310°CA bTDC, and a later injection timing that the liner, 220°CA bTDC start of injection timings (SOI) respectively. The earlier 310 SOI injection increases soot, whereas the later 220°CA SOI targets the liner and increases wall-wetting.

This work employs a novel laser induced fluorescence (LIF) diagnostic to measure fuel dilution in a running single cylinder research engine operated at stochastic pre ignition (SPI) and non-SPI prone conditions. Measurements of LIF based fuel dilution are quantified over a range of engine loads and fuel injection timings for two fuels. The in situ LIF measurements of fuel/lubricant interactions illustrate regions of increased fuel dilution from fuel-wall interactions and is believed to be a fundamental underpinning to generating top ring zone liquid conditions conducive to SPI. A novel level of dye doped in the fuel, between 50 to 500 ppm was used as the fluorescence source, at engine operating speed of 2000r/min from 5 to 18 bar of IMEPg injection timings was swept for two fuels of varying volatility.

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability.

This work explores pathways to achieve diesel-like, high-efficiency combustion with stoichiometric 3-way catalyst compatible spark ignition (SI). A high stroke-to-bore engine design (1.5:1) with cooled exhaust gas recirculation (EGR) and high compression ratio (rc) was used to improve engine efficiency by up to 30% compared with a production turbocharged gasoline direct injection spark ignition engine. To achieve efficiency improvements, engine experiments were coupled with computational fluid dynamics simulations to guide and explain experimental trends between the original engine and the high stroke-to-bore ratio design (1.5:1). The effects of EGR and late intake valve closing (IVC) and fuel characteristics are investigated through their effects on knock mitigation. Direct injection of 91 RON E10 gasoline, 99 RON E0 gasoline, and liquified petroleum gas (i.e., propane/autogas) were evaluated with geometric rc ranging from 13.3:1 to 16.8:1.

This work explores the impact of the interaction of lubricant and fuel properties on the propensity for stochastic pre-ignition (SPI). Findings are based on statistically significant changes in SPI tendency and magnitude, as determined by measurements of cylinder pressure. Specifically, lubricant detergents, lubricant volatility, fuel volatility, fuel chemical composition, fuel-wall impingement, and engine load were varied to study the physical and chemical effects of fuel-lubricant interactions on SPI tendency. The work illustrates that at low loads, with fuels susceptible to SPI events, lubricant detergent package effects on SPI were non-significant. However, with changes to fuel distillation, fuel-wall impingement, and most importantly engine load, lubricant detergent effects could be observed even at reduced loads This suggests that there is a thermal effect associated with the higher load operation.

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.

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.

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.

Numerous studies have demonstrated that exhaust gas recirculation (EGR) can attenuate knock propensity in spark ignition (SI) engines at naturally aspirated or lightly boosted conditions [1]. In this study, we investigate the role of cooled EGR under higher load conditions with multiple fuel compositions, where highly retarded combustion phasing typical of modern SI engines was used. It was found that under these conditions, EGR attenuation of knock is greatly reduced, where EGR doesn’t allow significant combustion phasing advance as it does under lighter load conditions. Detailed combustion analysis shows that when EGR is added, the polytropic coefficient increases causing the compressive pressure and temperature to increase. At sufficiently highly boosted conditions, the increase in polytropic coefficient and additional trapped mass from EGR can sufficiently reduce fuel ignition delay to overcome knock attenuation effects.

Natural gas is a promising alternative fuel for internal combustion engines due to its rich reserves and low price, as well as good physical and chemical properties. Its low carbon structure and high octane number are beneficial for CO2 reduction and knock mitigation, respectively. Diesel and natural gas dual fuel combustion is a viable pathway to utilize natural gas in diesel engines. To achieve high efficiency and low emission combustion in a practical diesel engine over a wide range of operating conditions, understanding the performance responses to engine system parameter variations is needed. The controllability of two combustion strategies, diesel pilot ignition (DPI) and single injection reactivity controlled compression ignition (RCCI), were evaluated using the multi-dimension CFD simulation in this paper.

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.

This work explores the dependence of fuel ignition delay on stochastic pre-ignition (SPI). Findings are based on bulk gas thermodynamic state, where the effects of kinetically controlled bulk gas pre-spark heat release (PSHR) are correlated to SPI tendency and magnitude. Specifically, residual gas and low temperature PSHR chemistry effects and observations are explored, which are found to be indicative of bulk gas conditions required for strong SPI events. Analyzed events range from non-knocking SPI to knocking SPI and even detonation SPI events in excess of 325 bar peak cylinder pressure. The work illustrates that singular SPI event count and magnitude are found to be proportional to PSHR of the bulk gas mixture and residual gas fraction. Cycle-to-cycle variability in trapped residual mass and temperature are found to impose variability in singular SPI event count and magnitude.

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.