Search Results

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.

In the present work, five surrogate components (n-Hexadecane, n-Tetradecane, Heptamethylnonane, Decalin, 1-Methylnaphthalene) are proposed to represent liquid phase of diesel fuel, and another different five surrogate components (n-Decane, n-Heptane, iso-Octane, MCH (methylcyclohexane), Toluene) are proposed to represent vapor phase of diesel fuel. For the vapor phase, a 5-component surrogate chemical kinetic mechanism has been developed and validated. In the mechanism, a recently updated H2/O2/CO/C1 detailed sub-mechanism is adopted for accurately predicting the laminar flame speeds over a wide range of operating conditions, also a recently updated C2-C3 detailed sub-mechanism is used due to its potential benefit on accurate flame propagation simulation. For each of the five diesel vapor surrogate components, a skeletal sub-mechanism, which determines the simulation of ignition delay times, is constructed for species C4-Cn.

This study investigated the potential of generating reactive chemical species (including radicals) through pilot fuel injection in negative valve overlap for improving the combustion and emissions performances of spark ignition gasoline engines under low load and low speed operating conditions. Several Ford sub-models were used for simulating the physics and chemistry processes of injecting a small amount of fuel in NVO (negative valve overlap). Effects of different NVO degrees and different pilot injection timings, factors for fuel conversion were simulated and investigated. CO and H2 conversions during NVO, CO and H2 amounts before spark timing were used for comparing different schemes.

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.

The variability in gasoline energy content, though most frequently not a consumer concern, is an issue of concern for vehicle manufacturers in demonstrating compliance with regulatory requirements. Advancements in both vehicle technology, test methodology, and fuel formulations have increased the level of visibility and concern with regard to the energy content of fuels used for regulatory testing. The R factor was introduced into fuel economy calculations for vehicle certification in the late 1980s as a means of addressing batch-to-batch variations in the heating value of certification fuels and the resulting variations in fuel economy results. Although previous studies have investigated values of the R factor for modern vehicles through experimentation, subsequent engine studies have made clear that it is difficult to distinguish between the confounding factors that influence engine efficiency when R is being studied experimentally.

Six blendstocks identified by the Co-Optimization of Fuels & Engines Program were used to prepare fuel blends using a fixed blendstock for oxygenate blending and a target RON of 97. The blendstocks included ethanol, n-propanol, isopropanol, isobutanol, diisobutylene, and a bioreformate surrogate. The blends were analyzed and used to establish interaction factors for a non-linear molar blending model that was used to predict RON and MON of volumetric blends of the blendstocks up to 35 vol%. Projections of efficiency increase, volumetric fuel economy increase, and tailpipe CO2 emissions decrease were produced using two different estimation techniques to evaluate the potential benefits of the blendstocks. Ethanol was projected to provide the greatest benefits in efficiency and tailpipe CO2 emissions, but at intermediate levels of volumetric fuel economy increase over a smaller range of blends than other blendstocks.

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.

This work presents an application of two sub-models relative to chemical-kinetics-based turbulent pre-mixed combustion modeling approach on the simulation of burn rate and emissions of spark ignition engines. In present paper, the justification of turbulent pre-mixed combustion modeling directly based on chemical kinetics plus a turbulence model is given briefly. Two sub-models relative to this kind of pre-mixed combustion modeling approach are described generally, including a practical PRF (primary reference fuel) chemical kinetic mechanism which can correctly capture the laminar flame speed under a wide range of Ford SI (spark ignition) engines/operating conditions, and an advanced spark plug ignition model which has been developed by Ford recently.

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.

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.

The latent heat-of-vaporization (HoV) of blends of biofuel and hydrocarbon components into gasolines has recently experienced expanded interest because of the potential for increased HoV to increase fuel knock resistance in direct-injection (DI) engines. Several studies have been conducted, with some studies identifying an additional anti-knock benefit from HoV and others failing to arrive at the same conclusion. Consideration of these studies holistically shows that they can be grouped according to the level of fuel octane sensitivity variation within their fuel matrices. When comparing fuels of different octane sensitivity significant additional anti-knock benefits associated with HoV are sometimes observed. Studies that fix the octane sensitivity find that HoV does not produce additional anti-knock benefit. New studies were performed at ORNL and NREL to further investigate the relationship between HoV and octane sensitivity.

New regulations requiring increases in vehicle fuel economy are challenging automotive manufacturers to identify fuel-efficient engines for future vehicles. Lean gasoline direct injection (GDI) engines offer significant increases in fuel efficiency over the more common stoichiometric GDI engines already in the marketplace. However, particulate matter (PM) emissions from lean GDI engines, particularly during stratified combustion modes, are problematic for lean GDI technology to meet U.S. Environmental Protection Agency Tier 3 and other future emission regulations. As such, the control of lean GDI PM with wall-flow filters, referred to as gasoline particulate filter (GPF) technology, is of interest. Since lean GDI PM chemistry and morphology differ from diesel PM (where more filtration experience exists), the functionality of GPFs needs to be studied to determine the operating conditions suitable for efficient PM removal.

This work presents turbulent premixed combustion modeling in spark ignition engines using G-equation based turbulent combustion model. In present study, a turbulent flame speed expression proposed and validated in recent years by two co-authors of this paper is applied to the combustion simulation of spark ignition engines. This turbulent flame speed expression has no adjustable parameters and its constants are closely tied to the physics of scalar mixing at small scales. Based on this flame speed expression, a minor modification is introduced in this paper considering the fact that the turbulent flame speed changes to laminar flame speed if there is no turbulence. This modified turbulent flame speed expression is implemented into Ford in-house CFD code MESIM (multi-dimensional engine simulation), and is validated extensively.

The recirculation of gases from the crankcase and valvetrain can potentially lead to the entrainment of lubricant in the form of aerosols or mists. As boost pressures increase, the blow-by flow through both the crankcase and the valve cover increases. The resulting lubricant can then become part of the intake charge, potentially leading to fouling of intake components such as the intercooler and the turbocharger. The entrained aerosol which can contain the lubricant and soot may or may not have the same composition as the bulk lubricant. The complex aerodynamic processes that lead to entrainment can strip out heavy components or volatilize light components. Similarly, the physical size and numbers of aerosol particles can be dependent upon the lubricant formulation and engine speed and load. For instance, high rpm and load may increase not only the flow of gases but the amount of lubricant aerosol.