Search Results

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.

Recent research on thermal reciprocating engines has focused on the influence of lubricant oil on the combustion process, which can lead to highly undesired super-knock events. Low-Speed Pre-Ignition (LSPI) events severely limit the further development of Direct Injection Spark Ignition Engines (DISI), preventing high efficiencies from being achieved. However, there is still a lack of knowledge about the fundamental mechanisms leading to LSPI, due to the complex phenomena involved and the interaction between lubricant oil and fuel. Understanding how the presence of lubricant oil traces affects gasoline chemical reactivity is an essential step for performing successful numerical simulations aimed at predicting the onset of LSPI phenomena. Reaction mechanisms able to predict oil-fuel interaction have been proposed, but they are computationally demanding.

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.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

Engine experiments have revealed the importance of fuel condensation on the emission characteristics of low temperature combustion. However, direct in-cylinder experimental evidence has not been reported in the literature. In this paper, the in-cylinder condensation processes observed in optically accessible engine experiments are first illustrated. The observed condensation processes are then simulated using state-of-the-art multidimensional engine CFD simulations with a phase transition model that incorporates a well-validated phase equilibrium numerical solver, in which a thermodynamically consistent phase equilibrium analysis is applied to determine when mixtures become unstable and a new phase is formed. The model utilizes fundamental thermodynamics principles to judge the occurrence of phase separation or combination by minimizing the system Gibbs free energy.

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.

This study presents a computational framework to predict the outcome of combustion process based on a given RANS initial condition by performing statistical analysis of Sankaran number, Sa, and ignition regime theory proposed by Im et al. [1]. A criterion to predict strong auto-ignition/detonation a priori is used in this study, which is based on Sankaran-Zeldovich criterion. In the context of detonation, Sa is normalized by a sound speed, and is spatially calculated for the bulk mixture with temperature and equivalence ratio stratifications. The initial conditions from previous pre-ignition simulations were used to compute the spatial Sa distribution followed by the statistics of Sa including the mean Sa, the probability density function (PDF) of Sa, and the detonation probability, PD. Sa is found to be decreased and detonation probability increased significantly with increase of temperature.

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.

Super-knock that occurs in spark ignition (SI) engines is investigated using two-dimensional (2D) numerical simulations. The temperature, pressure, velocity, and mixture distributions are obtained and mapped from a top dead center (TDC) slice of full-cycle three-dimensional (3D) engine simulations. Ignition is triggered at one end of the cylinder and a hot spot of known temperature was used to initiate a pre-ignition front to study super-knock. The computational fluid dynamics code CONVERGE was used for the simulations. A minimum grid size of 25 μm was employed to capture the shock wave and detonation inside the domain. The Reynolds-averaged Navier-Stokes (RANS) method was employed to represent the turbulent flow and gas-phase combustion chemistry was represented using a reduced chemical kinetic mechanism for primary reference fuels. A multi-zone model, based on a well-stirred reactor assumption, was used to solve the reaction terms.

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.

The experimental in-cylinder combustion process was compared with the numerical simualtion for naphtha fuel under conventional compression ignition (CI) and partially premixed combustion (PPC) conditions. The start of injection timing (SOI) with the single injection strategy was changed from late of −10 CAD aTDC to early of −40 CAD aTDC. The three-dimensional full cycle engine combustion simulation was performed coupling with gas phase chemical kinetics by the CFD code CONVERGE™. The flame index was used for evaluating the combustion evolution of premixed flame and diffusion flame. The results show that the flame index could be used as an indicator for in-cylinder homogeneity evaluation. Hydroperoxyl shows a similar distribution with the premixed combustion. Formaldehyde could be used as an indicator for low temperature combustion.

Gasoline compression ignition (GCI) engines are becoming more popular alternative for conventional spark engines to harvest the advantage of high volatility. Recent experimental study demonstrated that high reactivity gasoline fuel can be operated in a conventional mixing controlled combustion mode producing lower soot emissions than that of diesel fuel under similar efficiency and NOx level [1]. Therefore, there is much interest in using gasoline-like fuels in compression ignition engines. In order to improve the fidelity of simulation-based GCI combustion system development, it is mandatory to enhance the prediction of spray combustion of gasoline-like fuels. The purpose of this study is to model the spray characteristics of high reactivity gasoline fuels and validate the models with experimental results obtained through an optically accessible constant volume vessel under vaporizing [2] and reactive conditions [3].

One way to develop an understanding of soot formation and oxidation processes that occur during direct injection and combustion in an internal combustion engine is to image the natural luminosity from soot over time. Imaging is possible when there is optical access to the combustion chamber. After the images are acquired, the next challenge is to properly interpret the luminous distributions that have been captured on the images. A major focus of this paper is to provide guidance on interpretation of experimental images of soot luminosity by explaining how radiation from soot is predicted to change as it is transmitted through the combustion chamber and to the imaging. The interpretations are only limited by the scope of the models that have been developed for this purpose. The end-goal of imaging radiation from soot is to estimate the amount of soot that is present.

Full cycle simulations of KAUST optical diesel engine were conducted in order to provide insights into the details of fuel spray, mixing, and combustion characteristics at different start of injection (SOI) conditions. Although optical diagnostics provide valuable information, the high fidelity simulations with matched parametric conditions improve fundamental understanding of relevant physical and chemical processes by accessing additional observables such as the local mixture distribution, intermediate species concentrations, and detailed chemical reaction rates. Commercial software, CONVERGE™, was used as the main simulation tool, with the Reynolds averaged Navier-Stokes (RANS) turbulence model and the multi-zone (SAGE) combustion model to compute the chemical reaction terms. SOI is varied from late compression ignition (CI) to early partially premixed combustion (PPC) conditions.

In this work, linear, non-linear and a generalized renormalization group (RNG) two-equation RANS turbulence models of the k-epsilon form were compared for the prediction of turbulent compressible flows in diesel engines. The object-oriented, multidimensional parallel code FRESCO, developed at the University of Wisconsin, was used to test the alternative models versus the standard k-epsilon model. Test cases featured the academic backward facing step and the impinging gas jet in a quiescent chamber. Diesel engine flows featured high-pressure spray injection in a constant volume vessel from the Engine Combustion Network (ECN), as well as intake flows in a high-swirl diesel engine. For the engine intake flows, a model of the Sandia National Laboratories 1.9L light-duty single cylinder optical engine was used.

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.