Search Results

A fundamentally based quasi-dimensional NOx emission model for spark ignition direct injection (SIDI) gasoline engines was developed. The NOx model consists of a chemical mechanism and three sub-models. The classical extended Zeldovich mechanism and N₂O pathway for NOx formation mechanism were employed as the chemical mechanism in the model. A characteristic time model for the radical species H, O and OH was incorporated to account for non-equilibrium of radical species during combustion. A model of homogeneity which correlates fundamental dimensionless numbers and mixing time was developed to model the air-fuel mixing and inhomogeneity of the charge. Since temperature has a dominant effect on NOx emission, a flame temperature correlation was developed to model the flame temperature during the combustion for NOx calculation. Measured NOx emission data from a single-cylinder SIDI research engine at different operating conditions was used to validate the NOx model.

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.

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.

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 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].

Fundamental simulations using DNS type procedures were used to investigate the ignition, combustion characteristics and the lift-off trends of a spatially evolving turbulent liquid fuel jet. In particular, the spatially evolving n-Heptane spray injected in a two-dimensional rectangular domain with an engine like environment was investigated. The computational results were compared to the experimental observations from an optical engine as reported in the literature. It was found that an initial fuel rich combustion downstream of the spray tip is followed by diffusion combustion. Investigations were also made to understand the effects of injection velocity, ambient temperature and the droplet radius on the lift-off length. For each of these parameters three different values in a given range were chosen. For both injection velocity and droplet radius, an increase resulted in a near linear increase in the lift-off length.

While Low Temperature Combustion (LTC) strategies such as Reactivity Controlled Compression Ignition (RCCI) exhibit high thermal efficiency and produce low NOx and soot emissions, low load operation is still a significant challenge due to high unburnt hydrocarbon (UHC) and carbon monoxide (CO) emissions, which occur as a result of poor combustion efficiencies at these operating points. Furthermore, the exhaust gas temperatures are insufficient to light-off the Diesel Oxidation Catalyst (DOC), thereby resulting in poor UHC and CO conversion efficiencies by the aftertreatment system. To achieve exhaust gas temperature values sufficient for DOC light-off, combustion can be appropriately phased by changing the ratio of gasoline to diesel in the cylinder, or by burning additional fuel injected during the expansion stroke through post-injection.

Diesel particulate samples were collected from a light duty engine operated at a single speed-load point with a range of biodiesel and conventional fuel blends. The oxidation reactivity of the samples was characterized in a laboratory reactor, and BET surface area measurements were made at several points during oxidation of the fixed carbon component of both types of particulate. The fixed carbon component of biodiesel particulate has a significantly higher surface area for the initial stages of oxidation, but the surface areas for the two particulates become similar as fixed carbon oxidation proceeds beyond 40%. When fixed carbon oxidation rates are normalized to total surface area, it is possible to describe the oxidation rates of the fixed carbon portion of both types of particulates with a single set of Arrhenius parameters. The measured surface area evolution during particle oxidation was found to be inconsistent with shrinking sphere oxidation.

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.

Fundamental simulations in a quiescent cell under adiabatic conditions were made to understand the effect of temperature, equivalence ratio and the components of the recirculated exhaust gas, viz., CO2 and H2O, on the combustion of n-Heptane. Simulations were made in single phase in which evaporated n-Heptane was uniformly distributed in the domain. Computations were made for two different temperatures and four different EGR levels. CO2 or H2O or N2was used as EGR. It was found that the initiation of the main combustion process was primarily determined by two competing factors, i.e., the amount of initial OH concentration in the domain and the specific heat of the mixture. Further, initial OH concentration can be controlled by the manipulating the ambient temperature in the domain, and the specific heat capacity of the mixture via the mixture composition. In addition to these, the pre combustion and the subsequent post combustion can also be controlled via the equivalence ratio.

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.

Detailed knowledge of hydrocarbon fuel combustion chemistry has grown tremendously in recent years. However, the gap between detailed chemistry and computational fluid dynamics (CFD) remains, because of the high cost of solving detailed chemistry in a large number of computational cells. This paper presents the results of applying a suite of techniques aimed at closing this gap. The techniques include use of a surrogate blend optimizer and a guided mechanism reduction methodology, as well as advanced methods for efficiently and accurately coupling the pre-reduced kinetic models with the multidimensional transport equations. The advanced methods include dynamic adaptive chemistry (DAC) and dynamic cell clustering (DCC) algorithms.

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.

Various single and split injection schemes are studied to provide a better understanding of fuel distribution during cold starting in DI diesel engines. Improved spray-wall interaction, fuel film and multicomponent vaporization models are used to analyze the combustion processes. Better combustion characteristics are obtained for the split injection schemes than with a single injection. An analysis of the fuel impingement processes identifies the mechanisms involved in producing the differences in vaporization and combustion of the fuel. A greater amount of splashing occurred for the split injections compared to a single injection. This behavior is attributed to the decreased film thickness (less dissipation of impingement energy), the decreased impingement area (obtained by increasing the impingement Weber number), and most importantly, the reduced frequency of drop impingement.

A modified version of the three dimensional CFD code KIVA-3, which accommodates moving valves and a moving piston crown, has been applied to a heavy-duty, four cycle, dual intake valve, direct injection, diesel engine The fluid domain encompasses an intake runner, two valved ports and a cylinder with a Mexican hat bowl-in-piston configuration In the first part of the study, the modified KIVA-3 code was used to simulate the flow through the port and cylinder of the engine without a piston The intake valves were cycled (300 rpm at the camshaft) Two turbulence models were compared in this part of the study, standard k - ε and the RNG modified k - ε as discussed in [11] The results were compared to particle image velocimetry (PIV) images Large scale flow features of the computer simulations agreed moderately well with the ensemble averaged flow bench results There was very little difference between the results from the two turbulence models Motored engine simulations including a piston were conducted to characterize the in-cylinder gas flow and mixing during the intake and compression strokes Characterization methods were developed which yield insight into the in-cylinder gas motion and mixing during both strokes Intake and residual gases are tracked separately, both large scale convection and turbulent mixing are investigated, flow critical points are examined to provide information about flow topology and turbulence production is correlated with the evolution of flow structures Results show that mixing of the intake and residual gases is very non-uniform Many complex flow structures develop during intake and are destroyed during compression However, several structures survive through compression and contribute to enhanced mixing near top dead center These significant structures have been identified and tracked back to intake The flow field near top dead center exhibits spatial inhomogeneities in temperature and small scale mixing parameters such as turbulence kinetic energy and its dissipation rate

The paper reports the application of a look-up table approach within a LES combustion modelling framework for the prediction of knock limit in a highly downsized turbocharged DISI engine. During experimental investigations at the engine test bed, high cycle-to-cycle variability was detected even for relatively stable peak power / full load operations of the engine, where knock onset severely limited the overall engine performance. In order to overcome the excessive computational cost of a direct chemical solution within a LES framework, the use of look-up tables for auto-ignition modelling perfectly fits with the strict mesh requirements of a LES simulation, with an acceptable approximation of the actual chemical kinetics. The model here presented is a totally stand-alone tool for autoignition analysis integrated with look-up table reading from detailed chemical kinetic schemes for gasoline.

A novel algebraic similarity model for subgrid scalar dissipation rate has been developed as part of the Large Eddy Simulation (LES) package KIVA3V-LES for diesel engine study. The model is proposed from an a priori study using Direct Numerical Simulation (DNS) of forced isotropic turbulence. In the a posteriori test, fully resolved turbulent passive scalar field measurements are used to validate the model in actual engine flows. For reason of the length limit by SAE and the specific interest in engine applications, only a prior test and a posteriori test in engine flows are included in this paper. A posteriori tests in isotropic cube flow, turbulent round jet and flame cases will be presented in separate papers. An engine LES simulation of multi consecutive cycles was performed in this study.

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.

A novel 2-zone combustion system was examined at medium load operation consistent with loads in the light duty vehicle drive cycle (7.6 bar BMEP and 2600 rev/min). Pressure rise rate and noise can limit the part of the engine map where pre-mixed combustion strategies such as HCCI or RCCI can be used. The present 2-zone pistons have an axial projection that divides the near TDC volume into two regions (inner and outer) joined by a narrow communication channel defined by the squish height. Dividing the near TDC volume provides a means to prepare two fuel-air mixtures with different ignition characteristics. Depending on the fuel injection timing, the reactivity of the inner or outer volume can be raised to provide an ignition source for the fuel-air mixture in the other, less reactive volume. Multi-dimensional CFD modeling was used to design the 2-zone piston geometry examined in this study.

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.