Matthew J. McNenly

An accelerated multi-zone model for engine cycle simulation of homogeneous charge compression ignition combustion

We have developed an accelerated multi-zone model for engine cycle simulation (AMECS) of homogeneous charge compression ignition (HCCI) combustion. This model incorporates chemical kinetics and is intended for use in system-level simulation software. A novel methodology to capture thermal stratification in the multi-zone model is proposed. The methodology calculates thermal stratification inside the cylinder based on a single computational fluid dynamics (CFD) calculation for motored conditions. CFD results are used for tuning zone heat loss multipliers that characterize wall heat loss from each individual engine zone based on the assumption that these heat loss multipliers can then be used at operating conditions different from those used in the single CFD run because the functional form of thermal stratification is more dependent on engine geometry than on operating conditions. The model is benchmarked against detailed CFD calculations and fully coupled HCCI CFD chemical kinetics calculations. The results indicate that the heat loss multiplier approach accurately predicts thermal stratification during the compression stroke and (therefore) HCCI combustion. The AMECS model with the thermal stratification methodology and reduced gasoline chemical kinetics shows good agreement with boosted gasoline HCCI experiments over a range of operating conditions, in terms of in-cylinder pressure and heat release rate predictions. The computational advantage of this method derives from the need for only a single motoring CFD run for a given engine, which makes the method very well suited for rapid HCCI calculations in system-level codes such as GT-Power, where it is often desirable to evaluate consecutive engine cycles.

Detailed Chemical Kinetic Modeling of Iso-octane SI-HCCI Transition

We describe a CHEMKIN-based multi-zone model that simulates the expected combustion variations in a single-cylinder engine fueled with iso-octane as the engine transitions from spark-ignited (SI) combustion to homogenous charge compression ignition (HCCI) combustion. The model includes a 63-species reaction mechanism and mass and energy balances for the cylinder and the exhaust flow. For this study we assumed that the SI-to-HCCI transition is implemented by means of increasing the internal exhaust gas recirculation (EGR) at constant engine speed. This transition scenario is consistent with that implemented in previously reported experimental measurements on an experimental engine equipped with variable valve actuation. We find that the model captures many of the important experimental trends, including stable SI combustion at low EGR (-0.10), a transition to highly unstable combustion at intermediate EGR, and finally stable HCCI combustion at very high EGR (-0.75). Remaining differences between the predicted and experimental instability patterns indicate that there is further room for model improvement.

Integration Strategies for Efficient Multizone Chemical Kinetics Models

Three integration strategies are developed and tested for the stiff, ordinary differential equation (ODE) integrators used to solve the fully coupled multizone chemical kinetics model. Two of the strategies tested are found to provide more than an order of magnitude of improvement over the original, basic level of usage for the stiff ODE solver. One of the faster strategies uses a decoupled, or segregated, multizone model to generate an approximate Jacobian. This approach yields a 35-fold reduction in the computational cost for a 20 zone model. Using the same approximate Jacobian as a preconditioner for an iterative Krylov-type linear system solver, the second improved strategy achieves a 75-fold reduction in the computational cost for a 20 zone model. The faster strategies achieve their cost savings with no significant loss of accuracy. The pressure, temperature and major species mass fractions agree with the solution from the original integration approach to within six significant digits; and the radical mass fractions agree with the original solution to within four significant digits. The faster strategies effectively change the cost scaling of the multizone model from cubic to quadratic, with respect to the number of zones. As a consequence of the improved scaling, the 40 zone model offers more than a 250-fold cost savings over the basic calculation.