Anne Kösters

A comprehensive numerical study of Diesel fuel spray formation with OpenFOAM

The accuracy and robustness of spray models and their implementation in current commercial CFD codes vary substantially. However, common features are that the resulting spray penetration and levels of spray-generated turbulence - two factors that strongly influence the rate of heat released during combustion - are to a great extent grid size-dependent. In the work presented here a new kind of spray model has been implemented and thoroughly tested, under various ambient conditions, in the open source code OpenFOAM. In addition, since the turbulence model applied in simulations is known to strongly affect spray penetration rates, results obtained using both the standard k-e and RNG k-e models have been compared. In the new spray model, designated VSB2, the traditional Lagrangian parcel has been replaced by a so-called stochastic blob containing droplets with a distribution of sizes, rather than a number of uniform-sized droplets. These blobs do not interact with the grid directly, but through bubbles of locally determined size. One advantage of the VSB2 spray model is its less grid size dependence due to the interaction with the gas phase in the bubbles instead of the entire grid cell it currently occupy. Another benefit of the model is its robustness, resulting from rigorous calculation of equilibrium values for momentum and thermodynamic (saturation and temperature) parameters affecting transfer rates between the phases ensuring a bounded solution. Results obtained using the code are compared here with experimentally acquired data regarding spray penetration under various ambient conditions. Qualitative comparisons of the evolving spray shapes are also presented.

RANS predictions of turbulent diffusion flames: comparison of a reactor and a flamelet combustion model to the well stirred approach

The flame stabilisation process in turbulent non-premixed flames is not fully understood and several models have been developed to describe the turbulence–chemistry interaction. This work compares the performance of the multiple Representative Interactive Flamelet (mRIF) model, the Volume Reactor Fraction Model (VRFM), and the Well-Stirred reactor (WS) model in describing such flames. The predicted ignition delay and flame lift-off length of n-heptane sprays are compared to experimental results published within the Engine Combustion Network (ECN). All of the models predict the trend of ignition delay reasonably well. At a low gas pressure (42 bar) the ignition delay is overpredicted compared to the experimental data, but the difference between the models is not significant. However, the predicted lift-off lengths differ. At high pressure (87 bar) the difference between the models is small. All models slightly underpredict the lift-off length compared to the experimental data. At low gas pressure (42 bar) the mRIF model gives the best results. The VRFM and WS models predict excessively short lift-off lengths, but the VRFM model gives better results than the WS model. The flame structures of the models are also compared. The WS model and the VRFM model yield a well defined flame stabilisation point whereas the mRIF model does not. The flame of the mRIF model is more diffuse and the model is not able to predict flame propagation. All models were able to predict the experimental trends in lift-off and ignition delay, but certain differences between them are demonstrated.

A Numerical Study of the Effect of EGR on Flame Lift-off in n-Heptane Sprays Using a Novel PaSR Model Implemented in OpenFOAM

The effect of exhaust gas recirculation (EGR) on flame lift-off in non-stationary n-heptane sprays was studied under Diesel engine-like conditions using numerical simulation involving complex chemistry and a novel partially stirred reactor (PaSR) model of subgrid turbulence-chemistry interaction. The flame-stabilization mechanism is a result of complex physical and chemical interactions and cannot be described by a simplified theory. To leading order it is determined by the chemical reaction time at the leading edge, the turbulent diffusivity, and the flow velocity; so that there exists a balance between the local convection velocity and the triple-flame propagation speed. In this study of ignition and flame formation and stabilization processes, the VSB2 stochastic blob-and-bubble spray model was used in combination with the volume reactor fraction model (VRFM) implemented in OpenFOAM. The reacting volume fraction in the VRFM was determined by solving for mixture fraction, progress variable, and their variances in order to estimate the non-uniformities of the fluid cell; rather than simply taking the ratio of the mixing and chemistry time-scales. The chemistry is described by a reduced n-heptane mechanism with 36 species involved in 81 reactions. The simulated lift-off trends are compared to available experimental data from the Engine Combustion Network, Sandia National Laboratories [1].