A presumed joint pdf model for turbulent combustion with varying equivalence ratio

Modeling of the joint probability density function of the mixture fraction and progress variable with a given covariance value is studied. This modeling is validated using experimental and direct numerical simulation (DNS) data. A very good agreement with experimental data of turbulent stratified flames and DNS data of a lifted hydrogen jet flame is obtained. The effect of using this joint pdf modeling to calculate the mean reaction rate with a flamelet closure in Reynolds averaged Navier-Stokes (RANS) calculation of stratified flames is studied. The covariance effect is observed to be large within the flame brush. The results obtained from RANS calculations using this modeling for stratified jet- and rod-stabilized V-flames are discussed and compared to the measurements as a posteriori validation for the joint probability density function model with the flamelet closure. The agreement between the computed and measured values of flame and turbulence quantities is found to be good. © 2012 Copyright Taylor and Francis Group, LLC.

Effects of turbulent reynolds number on the performance of algebraic flame surface density models for large eddy simulation in the thin reaction zones regime: A direct numerical simulation analysis

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Ret has been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Ret and Ka before reaching an asymptotic value for large values of Ret and Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel'dovich flame thickness ηi/δz does not exhibit any significant dependence on Ret for the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width. © 2012 Mohit Katragadda et al.

Direct simulation of turbulent entrainment due to a plume impinging on a density interface

The law for turbulent entrainment due to plumes and jets impinging on a density interface is subject to significant uncertainty, with reported differences in entrainment rates up to a factor of 10. We report preliminary results obtained by Direct Numerical Simulation which are part of a PRACE project on turbulent entrainment carried out on JUGENE at Jülich, Germany. Various interface tracking methods are discussed and the entrainment coefficient is determined.

Classical plume theory: 1937-2010 and beyond

Developing a theoretical description of turbulent plumes, the likes of which may be seen rising above industrial chimneys, is a daunting thought. Plumes are ubiquitous on a wide range of scales in both the natural and the man-made environments. Examples that immediately come to mind are the vapour plumes above industrial smoke stacks or the ash plumes forming particle-laden clouds above an erupting volcano. However, plumes also occur where they are less visually apparent, such as the rising stream of warmair above a domestic radiator, of oil from a subsea blowout or, at a larger scale, of air above the so-called urban heat island. In many instances, not only the plume itself is of interest but also its influence on the environment as a whole through the process of entrainment. Zeldovich (1937, The asymptotic laws of freely-ascending convective flows. Zh. Eksp. Teor. Fiz., 7, 1463-1465 (in Russian)), Batchelor (1954, Heat convection and buoyancy effects in fluids. Q. J. R. Meteor. Soc., 80, 339-358) and Morton et al. (1956, Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A, 234, 1-23) laid the foundations for classical plume theory, a theoretical description that is elegant in its simplicity and yet encapsulates the complex turbulent engulfment of ambient fluid into the plume. Testament to the insight and approach developed in these early models of plumes is that the essential theory remains unchanged and is widely applied today. We describe the foundations of plume theory and link the theoretical developments with the measurements made in experiments necessary to close these models before discussing some recent developments in plume theory, including an approach which generalizes results obtained separately for the Boussinesq and the non-Boussinesq plume cases. The theory presented - despite its simplicity - has been very successful at describing and explaining the behaviour of plumes across the wide range of scales they are observed. We present solutions to the coupled set of ordinary differential equations (the plume conservation equations) that Morton et al. (1956) derived from the Navier-Stokes equations which govern fluid motion. In order to describe and contrast the bulk behaviour of rising plumes from general area sources, we present closed-form solutions to the plume conservation equations that were achieved by solving for the variation with height of Morton's non-dimensional flux parameter Γ - this single flux parameter gives a unique representation of the behaviour of steady plumes and enables a characterization of the different types of plume. We discuss advantages of solutions in this form before describing extensions to plume theory and suggesting directions for new research. © 2010 The Author. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Density stratified environments: The double-tank method

We present an alternative method of producing density stratifications in the laboratory based on the 'double-tank' method proposed by Oster (Sci Am 213:70-76, 1965). We refer to Oster's method as the 'forced-drain' approach, as the volume flow rates between connecting tanks are controlled by mechanical pumps. We first determine the range of density profiles that may be established with the forced-drain approach other than the linear stratification predicted by Oster. The dimensionless density stratification is expressed analytically as a function of three ratios: the volume flow rate ratio n, the ratio of the initial liquid volumes λ and the ratio of the initial densities ψ. We then propose a method which does not require pumps to control the volume flow rates but instead allows the connecting tanks to drain freely under gravity. This is referred to as the 'free-drain' approach. We derive an expression for the density stratification produced and compare our predictions with saline stratifications established in the laboratory using the 'free-drain' extension of Oster's method. To assist in the practical application of our results we plot the region of parameter space that yield concave/convex or linear density profiles for both forced-drain and free-drain approaches. The free-drain approach allows the experimentalist to produce a broad range of density profiles by varying the initial liquid depths, cross-sectional and drain opening areas of the tanks. One advantage over the original Oster approach is that density profiles with an inflexion point can now be established. © 2008 Springer-Verlag.

DNS of EGR-type turbulent flame in MILD condition

Three-dimensional direct numerical simulation (DNS) of exhaust gas recirculation (EGR)-type turbulent combustion operated in moderate and intense low-oxygen dilution (MILD) condition has been carried out to study the flame structure and flame interaction. In order to achieve adequate EGR-type initial/inlet mixture fields, partially premixed mixture fields which are correlated with the turbulence are carefully preprocessed. The chemical kinetics is modelled using a skeletal mechanism for methane-air combustion. The results suggest that the flame fronts have thin flame structure and the direct link between the mean reaction rate and scalar dissipation rate remains valid in the EGR-type combustion with MILD condition. However, the commonly used canonical flamelet is not fully representative for MILD combustion. During the flame-flame interactions, the heat release rate increases higher than the maximum laminar flame value, while the gradient of progress variable becomes smaller than laminar value. It is also proposed that the reaction rate and the scalar gradient can be used as a marker for the flame interaction. © 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Laminar and turbulent radial jets

In this paper the radial free jet which is produced by a continuous discharge of fluid from the space between two identical, parallel, circular, concentric discs into an infinite region of stagnant fluid of the same density and viscosity is investigated. Both laminar and turbulent jets are considered with analytical solutions being obtained near to the origin of the jet and at large distances along the jet. These asymptotic solutions are matched using a computational technique, and the numerical predictions show very good agreement with all the available experimental data.

Long range exhaustion-a mathematical model for the axisymmetric air flow of a local exhaust ventilation hood assisted by a turbulent radial jet

Abstract-Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art Banged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air Bow to be reduced and provide a means of identifying and assessing the various parameters that control the air Bow. The mathematical model is formulated in terms of the Stokes steam function, ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained Bow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions. | Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art flanged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air flow to be reduced and provide a means of identifying and assessing the various parameters that control the air flow. The mathematical model is formulated in terms of the Stokes steam function, Ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained flow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions.

Flame interactions in turbulent premixed twin V-flames

Multiple flame-flame interactions in premixed combustion are investigated using direct numerical simulations of twin turbulent V-flames for a range of turbulence intensities and length scales. Interactions are identified using a novel automatic feature extraction (AFE) technique, based on data registration using the dual-tree complex wavelet transform. Information on the time, position, and type of interactions, and their influence on the flame area is extracted using AFE. Characteristic length and time scales for the interactions are identified. The effect of interactions on the flame brush is quantified through a global stretch rate, defined as the sum of flamelet stretch and interaction stretch contributions. The effects of each interaction type are discussed. It is found that the magnitude of the fluctuations in flamelet and interaction stretch are comparable, and a qualitative sensitivity to turbulence length scale is found for one interaction type. Implications for modeling are discussed. © 2013 Copyright Taylor and Francis Group, LLC.

Tomographic reconstruction of OH* chemiluminescence in two interacting turbulent flames

The tomographic reconstruction of OH* chemiluminescence was performed on two interacting turbulent premixed bluff-body stabilized flames under steady flow conditions and acoustic excitation. These measurements elucidate the complex three-dimensional (3D) vortex-flame interactions which have previously not been accessible. The experiment was performed using a single camera and intensifier, with multiple views acquired by repositioning the camera, permitting calculation of the mean and phase-averaged volumetric OH* distributions. The reconstructed flame structure and phase-averaged dynamics are compared with OH planar laser-induced fluorescence and flame surface density measurements for the first time. The volumetric data revealed that the large-scale vortex-flame structures formed along the shear layers of each flame collide when the two flames meet, resulting in complex 3D flame structures in between the two flames. With a fairly simple experimental setup, it is shown that the tomographic reconstruction of OH* chemiluminescence in forced flames is a powerful tool that can yield important physical insights into large-scale 3D flame dynamics that are important in combustion instability. © 2013 IOP Publishing Ltd.