Search Results

Carbon-neutral (CN) fuels will be part of the solution to reducing global warming effects of the transportation sector, along with electrification. CN fuels such as hydrogen, ammonia, biofuels, and e-fuels can play a primary role in some segments (aviation, shipping, heavy-duty road vehicles) and a secondary role in others (light-duty road vehicles). The composition and properties of these fuels vary substantially from existing fossil fuels. Fuel effects on performance and emissions are complex, especially when these fuels are blended with fossil fuels. Predictively modeling the combustion of these fuels in engine and combustor CFD simulations requires accurate representation of the fuel blends. We discuss a methodology for matching the targeted fuel properties of specific CN fuels, using a blend of surrogate fuel components, to form a fuel model that can accurately capture fuel effects in an engine simulation.

Spray-wall interactions in diesel engines have a strong influence on turbulent flow evolution and mixing, which influences the engine’s thermal efficiency and pollutant-emissions behavior. Previous optical experiments and numerical investigations of a stepped-lip diesel piston bowl focused on how spray-wall interactions influence the formation of squish-region vortices and their sensitivity to injection timing. Such vortices are stronger and longer-lived at retarded injection timings and are correlated with faster late-cycle heat release and soot reductions, but are weaker and shorter-lived as injection timing is advanced. Computational fluid dynamics (CFD) simulations predict that piston bowls with more space in the squish region can enhance the strength of these vortices at near-TDC injection timings, which is hypothesized to further improve peak thermal efficiency and reduce emissions. The dimpled stepped-lip (DSL) piston is such a design.

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.

A one-dimensional, non-equilibrium, compressible law of the wall model is proposed to increase the accuracy of heat transfer predictions from computational fluid dynamics (CFD) simulations of internal combustion engine flows on engineering grids. Our 1D model solves the transient turbulent Navier-Stokes equations for mass, momentum, energy and turbulence under the thin-layer assumption, using a finite-difference spatial scheme and a high-order implicit time integration method. A new algebraic eddy-viscosity closure, derived from the Han-Reitz equilibrium law of the wall, with enhanced Prandtl number sensitivity and compressibility effects, was developed for optimal performance. Several eddy viscosity sub-models were tested for turbulence closure, including the two-equation k-epsilon and k-omega, which gave insufficient performance.

Accurate modeling of the characteristics of diesel-engine combustion leads to more efficient design. Accurate modeling in turn depends on correctly capturing spray dynamics, turbulence, and fuel chemistry. This work presents a computational fluid dynamics (CFD) investigation of a well characterized small-bore direct injection diesel engine at Sandia National Laboratories’ Combustion Research Facility. The engine has been studied for two piston-bowls geometries and various injection timings. Simulation of these conditions test the predictive capabilities of our approach to diesel engine modeling using Ansys Forte. An experimental database covering a wide range of operating conditions is provided by the Engine Combustion Network for this engine, which is used to validate our modeling approach. Automatic and solution-adaptive meshing is used, and the recommended settings are discussed.

Accurately predicting the evolution of soot mass and soot particle numbers under engine conditions is critical to advanced engine design. A detailed soot-chemistry model that can capture soot under gasoline and diesel conditions without tuning is necessary for such predictions. Building confidence in the predictive usage of the chemistry in engine simulations requires validating the soot kinetics over a wide range of operating conditions and fuels, using data from different experimental techniques, and using sources from laboratory flames to engines. This validation study focuses on a soot-chemistry model that considers multiple nucleation, growth, and oxidation reaction pathways. It involves 14 gas-phase precursors and considers the effect of different soot-particle surface sites.

Appropriate spray modeling in multidimensional simulations of diesel engines is well known to affect the overall accuracy of the results. More and more accurate models are being developed to deal with drop dynamics, breakup, collisions, and vaporization/multiphase processes; the latter ones being the most computationally demanding. In fact, in parallel calculations, the droplets occupy a physical region of the in-cylinder domain, which is generally very different than the topology-driven finite-volume mesh decomposition. This makes the CPU decomposition of the spray cloud severely uneven when many CPUs are employed, yielding poor parallel performance of the spray computation. Furthermore, mesh-independent models such as collision calculations require checking of each possible droplet pair, which leads to a practically intractable O(np2/2) computational cost, np being the total number of droplets in the spray cloud, and additional overhead for parallel communications.

Flash boiling occurs in sprays when the ambient gas pressure is lower than the saturation pressure of the injected fuel. In the present work, a numerical study was conducted to investigate solid-cone spray behaviors under various flash-boiling conditions. A new spray cone angle correlation that is a function of injection parameters was developed and used for spray initialization at the nozzle exit to capture plume interactions and the global spray shape. The spray-breakup regime control was adjusted to enable catastrophic droplet breakup, characterized by Rayleigh-Taylor (RT) breakup, near the nozzle exit. The model was validated against experimental spray data from five different injectors, including both multi-hole and single-hole injectors, with injection pressure varying from 100 to 200 bar.

Lagrangian spray modeling represents a critical boundary condition for multidimensional simulations of in-cylinder flow structure, mixture formation and combustion in internal combustion engines. Segregated models for injection, breakup, collision and vaporization are usually employed to pass appropriate momentum, mass, and energy source terms to the gas-phase solver. Careful calibration of each sub-model generally produces appropriate results. Yet, the predictiveness of this modeling approach has been questioned by recent experimental observations, which showed that at trans- and super-critical conditions relevant to diesel injection, classical atomization and vaporization behavior is replaced by a mixing-controlled phase transition process of a dense fluid. In this work, we assessed the shortcomings of classical spray modeling with respect to real-gas and phase-change behavior, employing a multicomponent phase equilibrium solver and liquid-jet theory.

In this work we studied the effects of piston bowl design on combustion in a small-bore direct-injection diesel engine. Two bowl designs were compared: a conventional, omega-shaped bowl and a stepped-lip piston bowl. Experiments were carried out in the Sandia single-cylinder optical engine facility, with a medium-load, mild-boosted operating condition featuring a pilot+main injection strategy. CFD simulations were carried out with the FRESCO platform featuring full-geometric body-fitted mesh modeling of the engine and were validated against measured in-cylinder performance as well as soot natural luminosity images. Differences in combustion development were studied using the simulation results, and sensitivities to in-cylinder flow field (swirl ratio) and injection rate parameters were also analyzed.

Sector mesh modeling is the dominant computational approach for combustion system design optimization. The aim of this work is to quantify the errors descending from the sector mesh approach through three geometric modeling approaches to an optical diesel engine. A full engine geometry mesh is created, including valves and intake and exhaust ports and runners, and a full-cycle flow simulation is performed until fired TDC. Next, an axisymmetric sector cylinder mesh is initialized with homogeneous bulk in-cylinder initial conditions initialized from the full-cycle simulation. Finally, a 360-degree azimuthal mesh of the cylinder is initialized with flow and thermodynamics fields at IVC mapped from the full engine geometry using a conservative interpolation approach. A study of the in-cylinder flow features until TDC showed that the geometric features on the cylinder head (valve tilt and protrusion into the combustion chamber, valve recesses) have a large impact on flow complexity.

Aromatics constitute a significant portion of refinery fuels. Characterizing the impact of various aromatic components on combustion and emissions facilitates formulation of surrogate fuels for engine simulations. The impact of blending aromatics in fuel surrogates is usually nonlinear for ignition characteristics responsible for knocking in spark engines and for combustion phasing in diesel engines. In this work, we have characterized the behavior of nine aromatics components under engine-relevant conditions. A self-consistent and validated detailed kinetics mechanism has been developed for gasoline and diesel surrogates that contains toluene, ethylbenzene, n-propylbenzene, n-butylbenzene, isomers of xylene, 1,2,4-trimethylbenzene, and 1-methylnaphthalene. Numerical experiments using 0-D and 1-D models have been performed to study the relative behavior of these aromatics for different reacting conditions.

This study focuses on Lagrangian spray models that are commonly used in engine CFD simulations. In modeling sprays, droplet collision is one of the physical phenomena that must be accounted for. There are two main parts of droplet collision models for sprays - detecting colliding pairs of droplets and predicting the outcomes of these collisions. For the first part, we focus on the efficiency of the algorithm. We present an implementation of the arbitrary adaptive collision mesh model of Hou and Schmidt [1], and examine its efficiency in dealing with large simulations. Through theoretical analysis and numerical tests, we show that the computational cost of this model scales pseudo-linearly with respect to the number of parcels in the sprays. Regarding the second part, we examine the variations in existing phenomenological models used for predicting binary droplet collision outcomes. A quantitative accuracy metric is used to evaluate the models with respect to the experimental data set.

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.

Diesel piston bowl geometry can affect turbulent mixing and therefore it impacts heat-release rates, thermal efficiency, and soot emissions. The focus of this work is on the effects of bowl geometry and injection timing on turbulent flow structure. This computational study compares engine behavior with two pistons representing competing approaches to combustion chamber design: a conventional, re-entrant piston bowl and a stepped-lip piston bowl. Three-dimensional computational fluid dynamics (CFD) simulations are performed for a part-load, conventional diesel combustion operating point with a pilot-main injection strategy under non-combusting conditions. Two injection timings are simulated based on experimental findings: an injection timing for which the stepped-lip piston enables significant efficiency and emissions benefits, and an injection timing with diminished benefits compared to the conventional, re-entrant piston.

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.

Advanced research in Spark-ignition (SI) engines has been focused on dilute-combustion concepts. For example, exhaust-gas recirculation is used to lower both fuel consumption and pollutant emissions while maintaining or enhancing engine performance, durability and reliability. These advancements achieve higher engine efficiency but may deteriorate combustion stability. One symptom of instability is a large cycle-to-cycle variation (CCV) in the in-cylinder flow and combustion metrics. Large-eddy simulation (LES) is a computational fluid dynamics (CFD) method that may be used to quantify CCV through numerical prediction of the turbulent flow and combustion processes in the engine over many engine cycles. In this study, we focus on evaluating the capability of LES to predict the in-cylinder flows and gas exchange processes in a motored SI engine installed with a transparent combustion chamber (TCC), comparing with recently published data.

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.