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.

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.

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.

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.

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.

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.

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.

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.

In an attempt to increase efficiency and lower critical and highly regulated emissions (i.e., NOx, PM and CO2) many advanced combustion strategies have been investigated. Most of the current strategies fall into the category of low temperature combustion (LTC), which allow emissions mandates to be met in-cylinder along with anticipated reduction in cost and complexity. These strategies, such as homogeneous charge compression ignition (HCCI), premixed charge compression ignition (PCCI), partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI), use early injection timings, resulting in a highly lean charge with increased specific heat ratios to improve thermal efficiency and reduce PM emissions. Lower combustion temperatures also avoid the activation of NOx formation reactions.

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.

Fuel-specific differences in exhaust gas recirculation (EGR) dilution tolerance are studied in a modern, direct-injection single-cylinder research engine. A total of 6 model fuel blends are examined at a constant research octane number (RON) of 95 using n-heptane, isooctane, toluene, and ethanol. Laminar flame speeds for these mixtures, which are calculated using two different methods (an energy fraction mixing rule and a detailed kinetic simulation), span a range of about 6 cm/s. A nominal load of 350 kPa IMEPg at 2000 rpm is maintained with constant fueling and varying CA50 from 8-20 CAD aTDCf. EGR is increased until a COV of IMEP of 5% is reached. The results illustrate that flame speed affects EGR dilution tolerance; fuels with increased flame speeds have increased EGR tolerance. Specifically, flame speed correlates most closely to the initial flame kernel growth, measured as the time of ignition to 5% mass fraction burned.

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.

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 present experimental engine efficiency study explores the effects of intake pressure and temperature, and premixed and global equivalence ratios on gross thermal efficiency (GTE) using the reactivity controlled compression ignition (RCCI) combustion strategy. Experiments were conducted in a heavy-duty single-cylinder engine at constant net load (IMEPn) of 8.45 bar, 1300 rev/min engine speed, with 0% EGR, and a 50% mass fraction burned combustion phasing (CA50) of 0.5°CA ATDC. The engine was port fueled with E85 for the low reactivity fuel and direct injected with 3.5% 2-ethylhexyl nitrate (EHN) doped into 91 anti-knock index (AKI) gasoline for the high-reactivity fuel. The resulting reactivity of the enhanced fuel corresponds to an AKI of approximately 56 and a cetane number of approximately 28. The engine was operated with a wide range of intake pressures and temperatures, and the ratio of low- to high-reactivity fuel was adjusted to maintain a fixed speed-phasing-load condition.

The present study experimentally investigates spark-ignited combustion with 87 AKI E0 gasoline in its neat form and in mid-level alcohol-gasoline blends with 24% vol./vol. iso-butanol-gasoline (IB24) and 30% vol./vol. ethanol-gasoline (E30). A single-cylinder research engine is used with a low and high compression ratio of 9.2:1 and 11.85:1 respectively. The engine is equipped with hydraulically actuated valves, laboratory intake air, and is capable of external exhaust gas recirculation (EGR). All fuels are operated to full-load conditions with λ=1, using both 0% and 15% external cooled EGR. The results demonstrate that higher octane number bio-fuels better utilize higher compression ratios with high stoichiometric torque capability. Specifically, the unique properties of ethanol enabled a doubling of the stoichiometric torque capability with the 11.85:1 compression ratio using E30 as compared to 87 AKI, up to 20 bar IMEPg at λ=1 (with 15% EGR, 18.5 bar with 0% EGR).

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

The paper reports the application of Proper Orthogonal Decomposition (POD) to LES calculations for the analysis of combustion and knock tendency in a highly downsized turbocharged GDI engine that is currently under production. In order to qualitatively match the cyclic variability of the combustion process, Large-Eddy Simulation (LES) of the closed-valve portion of the cycle is used with cycle-dependent initial conditions from a previous multi-cycle analysis [1, 2, 3]. Detailed chemical modelling of fuel's auto-ignition quality is considered through an ad-hoc implemented look-up table approach, as a trade-off between the need for a reasonable representation of the chemistry and that of limiting the computational cost of the LES simulations. Experimental tests were conducted operating the engine at knock-limited spark advance (KLSA) and the proposed knock model was previously validated for such engine setup [3].