LES of the Cambridge Stratified Swirl Burner using a Sub-grid pdf Approach

The sub-grid scale probability density function equation is rearranged in order to separate the resolved and sub-grid-scale (sgs) contributions to the sgs mixing term. This allows modelling that is consistent with the limiting case of negligible sub-grid scale variations, a property required for applications to laboratory premixed flames. The new method is applied to the Cambridge Stratified Swirl Burner for 6 operating conditions, 2 isothermal and 4 burning, with varying degrees of swirl and mixture stratification. The simulations are performed with the Large Eddy Simulation (LES) code BOFFIN in which the modelled pdf transport equation is solved using the Eulerian stochastic field method. Eight stochastic fields are used to account for the influence of the sub-grid fluctuations and the chemistry is modelled with a reduced version of the GRI 3.0 mechanism for methane involving 19 species and 15 reaction steps. The simulated velocities for both the isothermal and burning cases show good agreement with the experimental data. The measured temperature and major species profiles are also reproduced to a good accuracy.     

LES/CMC Simulations of Swirl-Stabilised Ethanol Spray Flames Approaching Blow-Off

Large Eddy Simulations (LES) with the Conditional Moment Closure (CMC) combustion model of swirling ethanol spray flames have been performed in conditions close to blow-off for which a wide database of experimental measurements is available for both flame and spray characterization. The solution of CMC equations exploits a three-dimensional unstructured code with a first order closure for chemical source terms. It is shown that LES/CMC is able to properly capture the flame structure at different conditions and agrees reasonably well with the measurements both in terms of mean flame shape and dynamic behaviour of the flame evaluated in terms of local extinctions and statistics of the lift-off height. Experimental measurements of the overall (liquid plus gaseous) mixture fraction, performed using the Laser-Induced Breakdown Spectroscopy technique, are also included allowing further assessment and validation of the numerical method. The sensitivity of the simulation results to the various boundary conditions is discussed.     

LES of Jets in Cross Flow and Its Application to a Gas Turbine Burner

LES computations of jets in cross flow (JICF) were performed. Experimental investigations reported in literature are reproduced with good agreement concerning the momentum field and the mixing of a passive scalar. The results validate the ability of the present LES approach to compute fuel injection of the type JICF. LES computations of fuel injection in an industrial gas turbine burner are presented. This revised version was published online in July 2006 with corrections to the Cover Date.     

Prediction of Global Extinction Conditions and Dynamics in Swirling Non-premixed Flames Using LES/CMC Modelling

The Large Eddy Simulation (LES)/three-dimensional Conditional Moment Closure (CMC) model with detailed chemistry is applied to predict the operating condition and dynamics of complete extinction (blow-off) in swirling non-premixed methane flames. Using model constants previously selected to provide relatively accurate predictions of the degree of local extinction in the piloted jet flames Sandia D ?F, the error in the blow-off air velocity predicted by LES/3D-CMC in short, recirculating flames with strong swirl for a range of fuel flow rates is within 25 % of the experimental value, which is considered a new and promising result for combustion LES that has not been applied before for the prediction of the whole blow-off curve in complex geometries. The results also show that during the blow-off transient, the total heat release gradually decreases over a duration that agrees well with experiment. The evolution of localized extinction, reactive scalars and scalar dissipation rate is analyzed. It has been observed that a consistent symptom for flames approaching blow-off is the appearance of high-frequency and high-magnitude fluctuations of the conditionally filtered stoichiometric scalar dissipation rate, resulting in an increased fraction of local extinction over the stoichiometric mixture fraction iso-surfaces. It is also shown that the blow-off time changes with the different blow-off conditions.     

Simulation of Hydrogen Auto-Ignition in a Turbulent Co-flow of Heated Air with LES and CMC Approach

Large-Eddy Simulations (LES) with the first order Conditional Moment Closure (CMC) approach of a nitrogen-diluted hydrogen jet, igniting in a turbulent co-flowing hot air stream, are discussed. A detailed mechanism (nine species, 19 reactions) is used to represent the chemistry. Our study covers the following aspects: CFD mesh resolution; CMC mesh resolution; inlet boundary conditions and conditional scalar dissipation rate modelling. The Amplitude Mapping Closure for the conditional scalar dissipation rate produces acceptable results. We also compare different options to calculate conditional quantities in CMC resolution. The trends in the experimental observations are in general well reproduced. The auto-ignition length decreases with an increase in co-flow temperature and increases with increase in co-flow velocity. The phenomena are not purely chemically controlled: the turbulence and mixing play also affect the location of auto-ignition. In order to explore the effect of turbulence, two options were applied: random noise and turbulence generator based on digital filter. It was found that stronger turbulence promotes ignition.     

Feedback-controlled forcing in hybrid LES/RANS

The computational cost of large eddy simulation (LES) increases rapidly with the Reynolds number when applied to attached boundary layers. This problem can be avoided by use of a Reynolds-averaged Navier–Stokes (RANS) model in the inner part of the boundary layer, which reduces the computational cost drastically. Such hybrid LES/RANS methods yield accurate results in general, but suffer from an artificial buffer layer and a shift in the velocity profile around the modeling interface. This velocity shift can be removed by use of additional forcing, but the results are very sensitive to the forcing amplitude.

The present paper proposes a feedback algorithm which efficiently finds the appropriate amplitude and thus yields accurate flow statistics. The feedback algorithm is relatively robust, both in that it is insensitive to the values of the parameters involved and that it yields accurate results with different forcing fields and for different Reynolds numbers. It is argued that the feedback algorithm is consistent with the underlying assumptions of hybrid LES/RANS and that it does not introduce additional empiricism into the method.     

Swirl measurements in a motor cylinder

 The present article analyzes different definitions of the swirl number and their application to control the mixing process in the cylinder of direct injection diesel engines. We first of all discuss the origins of the existing methods used to control swirl, which are usually based on the presence of a stirring device, or rotometer, in the cylinder head. A new possibility is then proposed for determining the swirl number, based on equivalent angular velocity, and a comparison made of results obtained using the different techniques available for measuring swirl in motor cylinders. Received: 16 December 1997/Accepted: 18 May 1998     

Modeling of Liquid Fuel Injection, Evaporation and Mixing in a Gas Turbine Burner Using Large Eddy Simulations

The interaction of turbulence, temperature fluctuation, liquid fuel transport, mixing and evaporation is studied by using Large Eddy Simulations (LES). To assess the accuracy of the different components of the methods we consider first isothermal, single phase flow in a straight duct. The results using different numerical methods incorporating dynamic Sub-Grid-Scale (SGS) models are compared with DNS and experimental data. The effects of the interactions among turbulence, temperature fluctuation, spray transport, evaporation and mixing of the gaseous fuel are studied by using different assumptions on the temperature field. It has been found that there are strong non-linear interactions among temperature-fluctuation, evaporation and turbulent mixing which require additional modeling if not full LES is used. The developed models and methods have been applied to a gas turbine burner into which liquid fuel is injected. The dispersion of the droplets in the burner is described. This revised version was published online in July 2006 with corrections to the Cover Date.     

Properties of the hybrid RANS/LES filter

One of the most important and challenging topics in the Large Eddy Simulation of turbulent flows is the connection of the LES technique to the well known and largely used RANS approach where the Navier–Stokes equations are Reynolds averaged. The hybridation of LES and RANS is not only important for its possible practical use, (a rational use of the computational means in different zones), but also from a theoretical point of view, and one possible procedure consists of blending RANS and LES models in the transition zone. In this paper a new filtering technique based on blending filters which transitions smoothly between LES and RANS is proposed and the associated universal model for the subgrid scale stresses is derived. PACS 47.27.Eq     

Hybrid URANS/LES Turbulence Simulation of Vortex Shedding Behind a Triangular Flameholder

In this study a detached eddy simulation (DES) model, which belongs to the group of hybrid URANS/LES turbulence models, is used for the simulation of vortex shedding behind a triangular obstacle. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS shear-stress transport (SST) model is used. In the other regions with higher grid resolution a LES model, which uses a transport equation for the turbulent subgrid energy, is applied. The DES model is first investigated for two standard test cases, namely decaying homogeneous isotropic turbulence and the backward facing step, respectively. For the decaying homogeneous isotropic turbulence test case the evolution of the energy spectra in wavenumber space for different times are studied for both the DES and a Smagorinsky type LES model. Different grid resolutions are analyzed with a special emphasis on the modeling constant connecting the filter length scale to the grid size. The results are compared to experimental data. The backward facing step test case is used to study the model behavior for a case with a transition region between a RANS modeling approach close to the wall and LES based modeling in the intense shear flow region. The final application is the simulation of the vortex shedding behind a triangular obstacle. First, the influence of the inlet condition formulation is studied in detail as they can have a significant influence especially for LES based models. Detailed comparisons between simulation and experiment for the flow structure past the obstacle and statistical quantities such as the shedding frequency are shown. Finally the additional temporal and spatial information provided by the DES model is used to show the predicted anisotropy of turbulence.     

Investigation of a Dynamic Hybrid RANS/LES Modelling Methodology for Finite-Volume CFD Simulations

This paper investigates a recently proposed dynamic hybrid RANS-LES framework using a general-purpose finite-volume flow solver. The new method is highly generalized, allowing coupling of any selected RANS model with any selected LES model and containing no explicit grid dependence in its formulation. Selected results are presented for three test cases: two-dimensional channel flow, backward facing step, and a nozzle flow relevant to biomedical applications. Comparison with experimental and DNS data, and with other hybrid RANS-LES approaches, highlights the advantages of the new method and suggests that further investigation is warranted.     

Reconstruction of turbulent fluctuations for hybrid RANS/LES simulations using a Synthetic-Eddy Method

A coupling methodology between an upstream Reynolds Averaged Navier–Stokes (RANS) simulation and a Large Eddy Simulation (LES) further downstream is presented. The focus of this work is on the RANS-to-LES interface inside an attached turbulent boundary layer, where an unsteady LES content has to be explicitly generated from a steady RANS solution. The performance of the Synthetic-Eddy Method (SEM), which generates realistic synthetic eddies at the inflow of the LES, is investigated on a wide variety of turbulent flows, from simple channel and square duct flows to the flow over an airfoil trailing edge. The SEM is compared to other existing methods of generation of synthetic turbulence for LES, and is shown to reduce substantially the distance required to develop realistic turbulence downstream of the inlet.     

Dissolution Dynamics of Liquid/Liquid Binary Mixtures Within a Micromodel

We report optical observations of the dissolution behaviour of glycerol/water, soybean oil/hexane, and isobutyric acid (IBA)/water binary mixtures within a glass micromodel built as a 2D regular network of capillary tubes with diameter of 0.2 mm. The micromodel is initially filled with solute and then is horizontally immersed into a thermostatic solvent-filled bath. The micromodel is open at its corners for solute dissolution to occur with no pressure gradients being applied. Our study shows that the solvent penetration into the micromodel is diffusion-dominated in completely miscible binary mixtures (glycerol/water and soybean oil/hexane). This is, however, non-Fickian diffusion with the dissolution rate, dV/dt, being proportional to $D^{1/3}t^{-0.4}$ for almost the entire duration of the experiment (V is the volume occupied by the solvent, D is the diffusion coefficient, and t is time). For the partially miscible IBA/water mixture the experiments performed at undercritical temperatures revealed that the diffusive transport was negligible despite the mixture being out of its thermodynamic equilibrium. The water phase penetrated into some of the channels, but IBA was never completely displaced/dissolved from the micromodel and numerous interfaces remained visible after very long-time periods.     

Quantitative Analysis of Flame Instabilities in a Hele-Shaw Burner

We show in this paper that a Hele-Shaw burner can be used for studying the development of premixed flame instabilities in a quasi-two dimensional configuration. It is possible to ignite a plane flame at the top of the cell, and to measure quantitatively the growth rates of the instability by image analysis. Experiments are performed with propane and methane-air mixtures. It is found that the most unstable wavelength, and the maximum linear growth rate of perturbations, directly measured in the present experiments, have the same order of magnitude as those previously measured on flames propagating freely downwards in wide tubes.     

LES of the flow around several cuboids in a row

This paper presents Large Eddy Simulations (LES) of flow around a four-vehicle platoon when one of the platoon members was forced to undergo in-line oscillations. The LES were made at the Reynolds number of 10⁵ based on the height of the vehicles. Combinations of two different frequencies corresponding to non-dimensional frequencies at the Strouhal numbers St₁ = 0.025 and St₂ = 0.013 and two oscillation amplitudes were used in this study. The methodology was validated by comparisons with data from previous experimental investigations. In order to highlight the dynamic effects, comparisons were made with steady results on a single vehicle and on a four-vehicle platoon. Large differences were found in the flow structures between quasi-steady and dynamic results. Furthermore, the behavior of the drag coefficient of the upstream neighbor of the oscillating model was investigated.     

A Hybrid RANS/LES Simulation of Turbulent Channel Flow

Hybrid models combining large eddy simulation (LES) with Reynolds-averaged Navier–Stokes (RANS) simulation are expected to be useful for wall modeling in the LES of high Reynolds number flows. Some hybrid simulations of turbulent channel flow have a common defect; the mean velocity profile has a mismatch between the RANS and LES regions due to a steep velocity gradient at the interface. This mismatch is reproduced and examined using a simple hybrid model; the Smagorinsky model is switched to a RANS model increasing the filter width. It is suggested that a rapid spatial variation in the eddy viscosity is responsible for an underestimate of the grid-scale shear stress and for the steep velocity gradient. To reduce the mean velocity mismatch a new scheme is proposed; additional filtering is introduced to define two kinds of velocity components at the interface between the two regions. The two components are used to remove inconsistency in the velocity equations due to a rapid variation in the filter width. Using the new scheme, simulations of channel flow are carried out with the simple hybrid model. It is shown that the grid-scale shear stress becomes large enough and most of the mean velocity mismatch is removed. Simulations for higher Reynolds numbers are carried out with the k–ε model and the one-equation subgrid-scale model. Although it is necessary to improve the turbulence models and the treatment of the buffer region, the new scheme is shown to be effective for reducing the mismatch and to be useful for developing better hybrid simulations. Received 5 April 2002 and accepted 8 January 2003 Published online 25 March 2003 Communicated by M.Y. Hussaini