Computation of 3D Turbulent Not-Premixed Reacting Flows Using an Implicit Unstructured Solver

2000 ◽  
Author(s):  
P. Adami ◽  
F. Martelli
Keyword(s):  
Bluff Body ◽  
Cfd Simulation ◽  
Turbulent Flame ◽  
Steady State Solution ◽  
Reacting Flows ◽  
Tvd Scheme ◽  
Reacting Flow ◽  
Riemann Solver ◽  
Turbulent Reactive Flows ◽  
Volume Method

A 3D CFD simulation of turbulent reactive flows is discussed. The original compressible version of the solver HybFlow designed for turbine rows investigation is here applied for low speed burning flow. A conserved scalar approach is considered to simulate the turbulent reacting flow field of non-premixed flames. The spatial discretization is based on an upwind finite volume method for unstructured grids using the Roe’s Riemann solver with a non-linear TVD scheme. The steady state solution is computed by means of an implicit relaxed Newton method. The linear solver is GMRES coupled with an ILU(0) preconditioning scheme. The turbulence chemistry interaction is described using a presumed β-PDF Flamelet approach. Two test applications are here presented to verify the methodology characteristics for a pilot-jet turbulent flame and a bluff-body stabilized flame both using CH4. A model combustor supplied with propane is also briefly shown as an example of application to a more realistic configuration.

Comparative Validation Study on Identification of Premixed Flame Transfer Function

2011 ◽  
Author(s):  
Luis Tay-Wo-Chong ◽  
Sebastian Bomberg ◽  
Ahtsham Ulhaq ◽  
Thomas Komarek ◽  
Wolfgang Polifke
Keyword(s):  
Experimental Data ◽  
Transfer Function ◽  
Cfd Simulation ◽  
Turbulent Flame ◽  
Maximum Amplitude ◽  
Quantitative Agreement ◽  
Combustion Model ◽  
Reacting Flow ◽  
Phase Lag ◽  
The Impact

The flame transfer function (FTF) of a premixed swirl burner was identified from time series generated with CFD simulation of compressible, turbulent, reacting flow at non-adiabatic conditions. Results were validated against experimental data. For large eddy simulation (LES), the Dynamically Thickened Flame combustion model with one step kinetics was used. For unsteady simulation in a Reynolds-averaged Navier-Stokes framework (URANS), the Turbulent Flame Closure model was employed. The FTF identified from LES shows quantitative agreement with experiment for amplitude and phase, especially for frequencies below 200 Hz. At higher frequencies, the gain of the FTF is underpredicted. URANS results show good qualitative agreement, capturing the main features of the flame response. However, the maximum amplitude and the phase lag of the FTF are underpredicted. Using a low-order network model of the test rig, the impact of the discrepancies in predicted FTFs on frequencies and growth rates of the lowest order eigenmodes were assessed. Small differences in predicted FTFs were found to have a significant impact on stability limits. Stability behavior in agreement with experimental data was achieved only with the LES-based flame transfer function.

Multi Eulerian PDF Transport Modelling of Turbulent Swirling Flame

2012 ◽  
Author(s):  
Abhinav Kapoor ◽  
Ashoke De ◽  
Rakesh Yadav
Keyword(s):  
Transport Equation ◽  
Transport Model ◽  
Turbulent Flame ◽  
Turbulence Models ◽  
Reacting Flows ◽  
Substantial Improvement ◽  
Reacting Flow ◽  
Transport Modelling ◽  
Solution Approach ◽  
Lagrangian Solution

The paper presents numerical investigation using Multi environmental Eulerian PDF (MEPDF) transport model for turbulence-chemistry interaction. A turbulent flame (SM1) from Sydney swirling burner database is simulated along with two isothermal cases (N29S054, N16S159) of different swirl numbers. MEPDF methodology, a probability density function (PDF) transport modeling, exploits the advantages of the PDF transport equation and is also computationally less expensive compared to popularly used Lagrangian solution approach of PDF transport equation. In the MEPDF approach, the PDF transport equation is represented by direct quadrature method of moments with presumed shape PDF and the closure of micro-mixing is achieved by interaction by exchange with mean (IEM) model. In the current work, the reacting flow results using MEPDF are reported for SM1 flame, which is a part of the database of turbulent reacting flows and widely considered as benchmark test cases for validating turbulent-chemistry interaction models. Initially, the non-reacting flows are simulated to properly choose the boundary conditions, turbulence models as well as the grid; followed by reacting flow calculations. SKE and RKE predictions show good agreement with each other while the other turbulence model exhibit substantially different behavior, especially for non-reacting case. However, RKE model exhibits substantial improvement in the case of reacting flows.

Characteristics of Bluff Body Stabilized Turbulent Premixed Flames

2011 ◽  
Author(s):  
Alejandro M. Briones ◽  
Balu Sekar ◽  
Hugh Thornburg
Keyword(s):  
Bluff Body ◽  
Three Dimensional ◽  
Leading Edge ◽  
Reacting Flows ◽  
Chemical Mechanism ◽  
Reacting Flow ◽  
Bluff Bodies ◽  
Two Dimensional ◽  
Aspect Ratios ◽  
Trailing Edges

Non-reacting and reacting flows past typical flameholders are modeled with URANS and LES. The continuity, momentum, energy, species, and turbulence governing equations are solved using two- and three-dimensional configurations. Either 2-step global or 44-step reduced chemical mechanism for C3H8-air combustion, accounting for turbulence-chemistry interaction, and with temperature- and species-dependent thermodynamic and transport properties is utilized. For square and rectangular bluff bodies the flow separates at the leading edges, whereas for triangular bluff body separation occurs only at the trailing edges. These bluff bodies exhibit two shear layers at the trailing edges that shed asymmetric vortices. For rectangular bluff bodies with aspect ratios (AR) less than 2.3 there is backflow from the wake. With increasing AR from unity, backflow is gradually diminished, and the von Ka´rma´n Strouhal number (StvK) decreases. For 2.0<AR<2.3, StvK jumps to a higher value and separation again occurs at the trailing edges for AR = 2.3. Further increase in AR decreases StvK again. The simulations with URANS qualitatively and quantitatively match experimental results for StvK vs. AR. Quantitative discrepancies are, however, found for AR≥2.3. In addition, two-dimensional non-reacting flows with URANS are sufficient to predict StvK. Moreover, two-dimensional simulations of reacting flow indicate that the flame promotes static and dynamic stability for AR = 1.0 and 2.3. The flame is dynamically unstable for AR = 2.0, exhibiting a von Ka´rma´n flow pattern. Stable flames anchored at the most downstream separation location (e.g., the flame anchored at AR = 1.0 is attached to the leading edge, whereas that of AR = 2.3 is attached to the trailing edge). Realizable k-ε URANS and LES simulations for the triangular cylinder closely match the experimental StvK for both non-reacting and reacting flows. Nonetheless, LES predicts a smaller recirculation length than k-ε URANS. LES predicts a flow field in which Be´rnard/von Ka´rma´n (BvK) instability is suppressed, whereas URANS predicts a competition between the Kelvin-Helmholtz (KH) instability and BvK.

Experimental investigation of isothermal and reacting flow characteristics downstream of axisymmetric bluff body stabilizers under stratified or fully premixed inlet mixture conditions

2020 ◽  
Author(s):  
Γεώργιος Πατεράκης
Keyword(s):  
Experimental Investigation ◽  
Turbulent Flow ◽  
Bluff Body ◽  
Operating Conditions ◽  
Reacting Flow ◽  
Wide Range ◽  
Flow Features ◽  
Air Mixing ◽  
The Impact ◽  
Stratified Flames

The current work describes an experimental investigation of isothermal and turbulent reacting flow field characteristics downstream of axisymmetric bluff body stabilizers under a variety of inlet mixture conditions. Fully premixed and stratified flames established downstream of this double cavity premixer/burner configuration were measured and assessed under lean and ultra-lean operating conditions. The aim of this thesis was to further comprehend the impact of stratifying the inlet fuelair mixture on the reacting wake characteristics for a range of practical stabilizers under a variety of inlet fuel-air settings. In the first part of this thesis, the isothermal mean and turbulent flow features downstream of a variety of axisymmetric baffles was initially examined. The effect of different shapes, (cone or disk), blockage ratios, (0.23 and 0.48), and rim thicknesses of these baffles was assessed. The variations of the recirculation zones, back flow velocity magnitude, annular jet ejection angles, wake development, entrainment efficiency, as well as several turbulent flow features were obtained, evaluated and appraised. Next, a comparative examination of the counterpart turbulent cold fuel-air mixing performance and characteristics of stratified against fully-premixed operation was performed for a wide range of baffle geometries and inlet mixture conditions. Scalar mixing and entrainment properties were investigated at the exit plane, at the bluff body annular shear layer, at the reattachment region and along the developing wake were investigated. These isothermal studies provided the necessary background information for clarifying the combustion properties and interpreting the trends in the counterpart turbulent reacting fields. Subsequently, for selected bluff bodies, flame structures and behavior for operation with a variety of reacting conditions were demonstrated. The effect of inlet fuel-air mixture settings, fuel type and bluff body geometry on wake development, flame shape, anchoring and structure, temperatures and combustion efficiencies, over lean and close to blow-off conditions, was presented and analyzed. For the obtained measurements infrared radiation, particle image velocimetry, laser doppler velocimetry, chemiluminescence imaging set-ups, together with Fouriertransform infrared spectroscopy, thermocouples and global emission analyzer instrumentation was employed. This helped to delineate a number of factors that affectcold flow fuel-air mixing, flame anchoring topologies, wake structure development and overall burner performance. The presented data will also significantly assist the validation of computational methodologies for combusting flows and the development of turbulence-chemistry interaction models.

scholarly journals CFD Simulation of the Dispersion of Exhaust Gases in a Traffic Loaded Street of Astana, Kazakhstan

2016 ◽  
Vol 9 (2) ◽  
pp. 158-166
Author(s):  
Ardak Akhatova ◽  
Assylan Kassymov ◽  
Meruyert Kazmaganbetova ◽  
Luis Ramon Rojas-Solórzano
Keyword(s):  
Cfd Simulation ◽  
Maximum Level ◽  
Weather Conditions ◽  
Control Measures ◽  
Convergence Criterion ◽  
Emission Rates ◽  
Wall Functions ◽  
Mesh Density ◽  
Volume Method ◽  
Future Work

The aim of this paper is to consider one of the most traffic-loaded regions of Astana city (Kazakhstan) and to determine the concentration of carbon-monoxide (CO) in the air during the peak hours. CFD analysis based on the SolidWorks-EFD platform was used to simulate the dispersion of contaminants given the estimated emission rates and weather conditions at the crossroad of Bogenbay Batyr and Zhenis Avenues in Astana. Turbulence prediction was based on k-ε model with wall functions. The governing equations were discretized using the finite volume method and a 2nd order spatial scheme. The mesh verification was based on 1% convergence criterion for a 50% of mesh density increment; air pressure near the wall of a selected building was chosen as the parameter to control the convergence. Numerical results are presented for prevailing conditions during all 4 seasons of the year, demonstrating that the highest levels of CO are recorded in summer and reach the values up to 11.2 ppm which are still lower than the maximum level admitted for humans. Nevertheless, obtained results show that Astana is gradually becoming a city that is likely to reach the critical levels of pollutants in the nearest future if control measures are not taken with enough anticipation. As for a future work, it is proposed to perform in-situ validation of specific scenarios to check and support the results obtained with CFD and to develop then specific policies for tackling the problem before it becomes evident.

scholarly journals A New Model Approach for Convective Wall Heat Losses in DQMOM-IEM Simulations for Turbulent Reactive Flows

2018 ◽  
Vol 141 (5) ◽  
Author(s):  
Yeshaswini Emmi ◽  
Andreas Fiolitakis ◽  
Manfred Aigner ◽  
Franklin Genin ◽  
Khawar Syed
Keyword(s):  
Experimental Data ◽  
Heat Losses ◽  
Reacting Flow ◽  
Jet Flame ◽  
New Model ◽  
Turbulent Reactive Flows ◽  
Novel Method ◽  
The One ◽  
The Impact ◽  
Model Approach

A new model approach is presented in this work for including convective wall heat losses in the direct quadrature method of moments (DQMoM) approach, which is used here to solve the transport equation of the one-point, one-time joint thermochemical probability density function (PDF). This is of particular interest in the context of designing industrial combustors, where wall heat losses play a crucial role. In the present work, the novel method is derived for the first time and validated against experimental data for the thermal entrance region of a pipe. The impact of varying model-specific boundary conditions is analyzed. It is then used to simulate the turbulent reacting flow of a confined methane jet flame. The simulations are carried out using the DLR in-house computational fluid dynamics code THETA. It is found that the DQMoM approach presented here agrees well with the experimental data and ratifies the use of the new convective wall heat losses model.

scholarly journals Turbulent Flame Shape Switching at Conditions Relevant for Gas Turbines

2019 ◽  
Vol 142 (1) ◽  
Author(s):  
Ivan Langella ◽  
Johannes Heinze ◽  
Thomas Behrendt ◽  
Lena Voigt ◽  
Nedunchezhian Swaminathan ◽  
...  
Keyword(s):  
Turbulent Combustion ◽  
Gas Turbines ◽  
Turbulent Flame ◽  
Interaction Model ◽  
Numerical Approach ◽  
Reacting Flow ◽  
Combustion Chambers ◽  
Lean Burn ◽  
Large Eddy ◽  
Flame Behavior

Abstract A numerical investigation is conducted to shed light on the reasons leading to different flame configurations in gas turbine (GT) combustion chambers of aeronautical interest. Large eddy simulations (LES) with a flamelet-based combustion closure are employed for this purpose to simulate the DLR-AT big optical single sector (BOSS) rig fitted with a Rolls-Royce developmental lean burn injector. The reacting flow field downstream this injector is sensitive to the intricate turbulent–combustion interaction and exhibits two different configurations: (i) a penetrating central jet leading to an M-shape lifted flame; or (ii) a diverging jet leading to a V-shaped flame. The LES results are validated using available BOSS rig measurements, and comparisons show the numerical approach used is consistent and works well. The turbulent–combustion interaction model terms and parameters are then varied systematically to assess the flame behavior. The influences observed are discussed from physical and modeling perspectives to develop physical understanding on the flame behavior in practical combustors for both scientific and design purposes.

scholarly journals Turbulent flame-vortex dynamics of bluff-body premixed flames

2021 ◽  
Vol 223 ◽  
pp. 28-41
Author(s):  
Marissa K. Geikie ◽  
Cal J. Rising ◽  
Anthony J. Morales ◽  
Kareem A. Ahmed
Keyword(s):  
Vortex Dynamics ◽  
Bluff Body ◽  
Turbulent Flame ◽  
Premixed Flames

Effect of the Soil Inclination on Natural Convection in Half-Elliptical Greenhouses

2020 ◽  
Vol 50 ◽  
pp. 70-78
Author(s):  
Djedid Taloub ◽  
Abdelkarim Bouras ◽  
Zied Driss
Keyword(s):  
Natural Convection ◽  
Cfd Simulation ◽  
Numerical Study ◽  
Boussinesq Approximation ◽  
Elliptic Cavity ◽  
Horizontal Cavity ◽  
Inclined Cavity ◽  
Volume Method ◽  
Heat Transfers ◽  
Work Done

A numerical study of the natural convection of laminar heat transfers in the stationary state in a half-elliptic inclined cavity, which represents a continuation of the work done, we studied the influence of the tilt of the cavity by varying the angle — entered 0 degrees, which corresponds to the horizontal cavity, up to 15 degrees. For each value of δ we varied the Rayleigh number from 2.13 103 to 106. The system of equations governing the problem solved numerically by the fluent calculation code based on the finite volume method. Based on the Boussinesq approximation. Both bottom and upper walls maintained at a constant temperature. The interest of this study is to see the influence of the tilt of the half-elliptic cavity on the structure of the flow and the distribution of temperature. These results can exploited in semi-elliptic agricultural greenhouses that rest on sloping soils. We chose a Prandtl number 0.71 that corresponds to the air. Keywords: Heat transfer; half-elliptical; Natural convection; Laminar flow; Multicellular; CFD simulation

Sign in / Sign up

Export Citation Format

Share Document