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.

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

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.

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.

In this work the modeling aspects of fuel vaporization are studied. To start with, the effects of vaporization model on engine simulations are studied. This is done by using two different fuel surrogates. Next a set of non-reacting spray simulations were performed under different ambient and operating conditions and for two different fuels. This was done for spray model validation and to look at the effect of vaporization model on liquid penetration length. Following an observed discrepancy in one of the spray cases, effect of ambient temperature on liquid length, two sensitivity analyses were performed. These analyses take into account the effects of each spray-sub model on vaporization and effects of spray breakup constants on liquid penetration. Using the results from the sensitivity analyses and linearized stability theory an empirical correction factor was developed to correct the spray behavior at low ambient temperatures.

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.

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.

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.

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.

Knock in a Rotax-914 engine was modeled and investigated using an improved version of the KIVA-3V code with a G-equation combustion model, together with a reduced chemical kinetics model. The ERC-PRF mechanism with 47 species and 132 reactions [1] was adopted to model the end gas auto-ignition in front of the flame front. The model was validated by a Caterpillar SI engine and a Rotax-914 engine in different operating conditions. The simulation results agree well with available experimental results. A new engineering quantified knock criterion based on chemical mechanism was then proposed. Hydroperoxyl radical (HO2) shows obvious accumulation before auto-ignition and a sudden decrease after auto-ignition. These properties are considered to be a good capability for HO2 to investigate engine knock problems.

Numerical investigation is carried out to explore various strategies of combustion mode switching in a diesel engine operating at high power. Numerical results are compared with high power single cylinder (CAT 3401E) experiments for combustion phasing and emission characteristics. In this study CFD calculations are carried out using the KIVA CFD code with Large Eddy Simulation turbulence model and Direct Chemistry Solver sub-models. The advanced turbulence and combustion sub-models enabled more realistic visualization of the effects of single-cycle mode switching on in-cylinder flow structures, fuel-air mixing behavior and combustion phasing. Two circumstances of mode switch are presented in this study. Mode switches are performed from traditional High Temperature Combustion to early injection PCCI combustion and vice versa. In this investigation several aspects of combustion control are investigated.

Numerical investigation is carried out in order to explore diesel combustion using advanced turbulence and combustion models. Turbulence is modeled by one-equation non-viscosity dynamic structure Large Eddy Simulation (LES) model. Sub-grid fuel-air mixing is calculated using a dynamic scale similarity sub-grid scalar dissipation model to represent the local state of combustion. Fuel-air mixing time scale is used in order to determine the local in-homogeneity and rate of mixing of fuel and air. Diesel combustion is studied and compared with experimental results for high power diesel engine setup at different conditions representing both low temperature combustion and traditional high temperature combustion regimes. Further studies are carried out in diesel engine to investigate in-cylinder fuel air mixing and the onset of ignition.