Search results

This study aims to present a basic mathematical model for investigating the structure of counter-flow non-premixed laminar flames propagating through uniformly-distributed organic fuel particles considering preheat, drying, vaporization, reaction and oxidizer zones.

Lycopodium particles and air are taken as biofuel and oxidizer, respectively. Dimensionalized and non-dimensionalized forms of mass and energy conservation equations are derived for each zone taking into account proper boundary and jump conditions. Subsequently, to solve the governing equations, an asymptotic method is used. For validation purpose, results achieved from the present analysis are compared with reliable data reported in the literature under certain conditions.

With regard to the comparisons, although different complex non-homogeneous differential equations are solved in this paper, acceptable agreements are observed. Finally, the impacts of significant parameters including fuel and oxidizer Lewis numbers, equivalence ratio, mass particle concentration, fuel and oxidizer mass fractions and lycopodium initial temperature on the flame temperature, flame front position and flow strain rate are elaborately explained.

An asymptotic method for mathematical modeling of counter-flow non-premixed multi-zone laminar flames propagating through lycopodium particles.

The purpose of this paper is to investigate the development process of the fire whirl in the fixed-frame facility and focus on the impacts of the fire whirl’s vortex core on the formation and flame structure of the fire whirl.

The complex turbulent reacting flame surface is captured by the large eddy simulation turbulence closure coupled with two sub-grid scale (SGS) kinetic schemes (i.e. the chemistry equilibrium and steady diffusion flamelet). Numerical predictions are validated thoroughly against the measurements by Lei et al. (2015) with excellent agreements. A double maximum tangential velocity refinement approach is proposed to quantify the vortex cores’ instantaneous location and region, addressing the missing definition in other studies.

The numerical results show that the transition process of the fire whirl is dominated by the vortex core movement, which is related to the centripetal force. The unsteadiness of the fully developed fire whirl was found depending on the instantaneous fluctuation of heat release rate. The steady diffusion flamelet scheme is essential to capture the instantaneous fluctuation. Furthermore, the axial velocity inside the vortex core is the key to determining the state of fire whirl.

Due to intensive interactions between buoyant fires and ambient rotating flow, the on-set and formation of fire whirl still remain largely elusive. This paper focused on the transition process of fire whirl between different development stages. This paper provides insights into the transition process from the inclined flame to the fire whirls based on the centripetal force.

This paper presented and compared two SGS kinetic schemes to resolve the fire whirl development process and the unsteadiness of its vortical structures. The modelling framework addresses the shortcoming of previous numerical studies where RANS turbulence closure and simplified combustion kinetics was adopted. Numerical results also revealed the fire whirl transition process and its relationship to centripetal force.

The purpose of this paper is to comprehensively evaluate the progress in the development of trapped vortex combustors (TVCs) in the past three decades. The review aims to identify the needs, predict the scope and discuss the challenges of numerical simulations in TVCs applied to gas turbines.

TVC is an emerging combustion technology for achieving low emissions in gas turbine combustors. The overall operation of such TVCs can be on very lean mixture ratio and hence it helps in achieving high combustion efficiency and low overall emission levels. This review introduces the TVC concept and the evolution of this technology in the past three decades. Various geometries that were explored in TVC research are listed and their operating principles are explained. The review then categorically arranges the progress in computational studies applied to TVCs.

Analyzing extensive literature on TVCs the review discusses results of numerical simulations of various TVC geometries. Numerical simulations that were used to optimize TVC geometry and to enhance mixing are discussed. Reactive flow studies to comprehend flame stability and emission characteristics are then listed for different TVC geometries.

To the best of the authors’ knowledge, this review is the first of its kind to discuss extensively the computational progress in TVC development specific to gas turbine engines. Earlier review on TVC covers a wide variety of applications including land-based gas turbines, supersonic Ramjets, incinerators and hence compromise on the depth of analysis given to gas turbine engine applications. This review also comprehensively group the numerical studies based on geometry, flow and operating conditions.

The widespread usage of magnetic nanoparticles (MNPs) requires their efficient synthesis during combustion process. This study aims to present a mathematical model for the oxidation of MNPs in a counter-flow non-premixed combustion system to produce MNPs, where the key sub-processes during the oxidation reaction are involved.

To accurately describe structure of flame and determine distributions of temperature and mass fractions of both reactants and products, equations of energy and mass conservations were solved based on the prevailing assumptions that three regions, i.e. preheating, reaction and oxidizer zones exist.

The numerical simulation was first validated against experimental data and characteristics of the combustion process are discussed. Eventually, the influences of crucial parameters such as reactant Lewis numbers, strain rate ratio, particle size, inert gas and thermophoretic force on structure of flame and combustion behavior were examined. The results show that maximum flame temperature can achieve 2,205 K. Replacing nitrogen with argon and helium as carrier gases can increase flame temperature by about 27% and 34%, respectively. Additionally, maximum absolute thermophoretic force was found at approximately 9.6 × 10–8 N.

To the best of authors’ knowledge, this is the first time to numerically model the preparation of MNPs in a counter-flow non-premixed combustion configuration, which can guide large-scale experimental work in a more effective way.

Safety improvement and pollutant reduction in many practical combustion systems and especially in aero-gas turbine engines require an adequate understanding of flame ignition and stabilization mechanisms. Improved software and hardware have opened up greater possibilities for translating basic knowledge and the results of experiments into better designs. The present study deals with the large eddy simulation (LES) of an ignition sequence in a conical shaped bluff-body stabilized burner involving a turbulent non-premixed flame. The purpose of this paper is to investigate the impact of spark location on ignition success. Particular attention is paid to the ease of handling of the numerical tool, the computational cost and the accuracy of the results.

The discrete particle ignition kernel (DPIK) model is used to capture the ignition kernel dynamics in its early stage of growth after the breakdown period. The ignition model is coupled with two combustion models based on the mixture fraction-progress variable formulation. An infinitely fast chemistry assumption is first done, and the turbulent fluctuations of the progress variable are captured with a bimodal probability density function (PDF) in the line of the Bray–Moss–Libby (BML) model. Thereafter, a finite rate chemistry assumption is considered through the flamelet-generated manifold (FGM) method. In these two assumptions, the classical beta-PDF is used to model the temporal fluctuations of the mixture fraction in the turbulent flow. To model subgrid scale stresses and residual scalars fluxes, the wall-adapting local eddy (WALE) and the eddy diffusivity models are, respectively, used under the low-Mach number assumption.

Numerical results of velocity and mixing fields, as well as the ignition sequences, are validated through a comparison with their experimental counterparts. It is found that by coupling the DPIK model with each of the two combustion models implemented in a LES-based solver, the ignition event is reasonably predicted with further improvements provided by the finite rate chemistry assumption. Finally, the spark locations most likely to lead to a complete ignition of the burner are found to be around the shear layer delimiting the central recirculation zone, owing to the presence of a mixture within flammability limits.

Some discrepancies are found in the radial profiles of the radial velocity and consequently in those of the mixture fraction, owing to a mismatch of the radial velocity at the inlet section of the computational domain. Also, unlike FGM methods, the BML model predicts the overall ignition earlier than suggested by the experiment; this may be related to the overestimation of the reaction rate, especially in the zones such as flame holder wakes which feature high strain rate due to fuel-air mixing.

This work is adding a contribution for ignition modeling, which is a crucial issue in various combustion systems and especially in aircraft engines. The exclusive use of a commercial computational fluid dynamics (CFD) code widely used by combustion system manufacturers allows a direct application of this simulation approach to other configurations while keeping computing costs at an affordable level.

This study provides a robust and simple way to address some ignition issues in various spark ignition-based engines, namely, the optimization of engines ignition with affordable computational costs. Based on the promising results obtained in the current work, it would be relevant to extend this simulation approach to spray combustion that is required for aircraft engines because of storage volume constraints. From this standpoint, the simulation approach formulated in the present work is useful to engineers interested in optimizing the engines ignition at the design stage.

The purpose of the paper is to summarize recent achievements and suggest further research directions in numerical studies of round free jets with particular attention on the influence of the inlet parameters (mean velocity, turbulence intensity, length and time scales) on the jet dynamics.

The large eddy simulation (LES) and direct numerical simulation (DNS) are regarded as accurate tools which can support expensive and requiring sophisticated measurements techniques experimental studies. In the paper, the authors present challenges and recent findings related to the LES and DNS of jet type flows in isothermal, heated, excited and reactive conditions.

LES of the isothermal jet allowed to identify the new jet instability mechanism leading to the self-sustained oscillations and to determine conditions required to trigger this phenomenon. Numerical simulation on the low-density round jet captured the phenomenon of absolute instability with a very good agreement with the experimental findings. LES/DNS of excited jet exhibited bifurcating and blooming jet and showed that the jet can be directly controlled by excitation frequency what is crucial issue also for flame shape control.

The paper shows complexity of seemingly simple jet type flow and proves that despite a huge interest in these flows and relatively deep knowledge on the jet dynamics there are still some open issues requiring further studies.

This study aims to investigate the insights of soot formation such as rate of soot coagulation, rate of soot nucleation, rate of soot surface growth and soot surface oxidation in ethylene/hydrogen/nitrogen diffusion jet flame at standard atmospheric conditions, which is very challenging to capture even with highly sophisticated measuring systems such as Laser Induced Incandescence and Planar laser-induced fluorescence. The study also aims to investigate the volume of soot in the flame using soot volume fraction and to understand the global correlation effect in the formation of soot in ethylene/hydrogen/nitrogen diffusion jet flame.

A large eddy simulation (LES) was performed using box filtered subgrid-scale tensor. A filtered and residual component of the governing equations such as continuity, momentum, energy and species are resolved and modeled, respectively. All the filtered and residual components are numerically solved using the ILU method by considering PISO pressure–velocity solver. All the hyperbolic flux uses the QUICK algorithm, and an elliptic flux uses SOU to evaluate face values. In all the cases, Courant–Friedrichs–Lewy (CFL) conditions are maintained unity.

The findings are as follows: soot volume fraction (SVF) as a function of a flame-normalized length for three different Reynolds number configurations (Re = 15,000, Re = 8,000 and Re = 5,000) using LES; soot gas phase and particulate phase insights such as rate of soot nucleation, rate of soot coagulation, rate of soot surface growth and soot surface oxidation for three different Reynolds number configurations (Re = 15,000, Re = 8,000 and Re = 5,000); and soot global correction using total soot volume in the flame volume as a function of Reynolds number and Froude number.

The originality of this study includes the following: coupling LES turbulent model with chemical equilibrium diffusion combustion conjunction with semi-empirical Brookes Moss Hall (BMH) soot model by choosing C6H6 as a soot precursor kinetic pathway; insights of soot formations such as rate of soot nucleation, soot coagulation rate, soot surface growth rate and soot oxidation rate for ethylene/hydrogen/nitrogen co-flow flame; and SVF and its insights study for three inlet fuel port configurations having the three different Reynolds number (Re = 15,000, Re = 8,000 and Re = 5,000).

The purpose of this paper is to numerically investigate the three-dimensional (3D) reacting turbulent two-phase flow field of a scaled swirl-stabilized gas turbine combustor using the commercial computational fluid dynamic (CFD) software ANSYS FLUENT. The first scope of the study aims to explicitly compare the predictive capabilities of two turbulence models namely Unsteady Reynolds Averaged Navier-Stokes and Scale Adaptive Simulation for a reasonable trade-off between accuracy of results and global computational cost when applied to simulate swirl-stabilized spray combustion. The second scope of the study is to couple chemical reactions to the turbulent flow using a realistic chemistry model and also to model the local chemical non-equilibrium(NEQ) effects caused by turbulent strain such as flame stretching.

Standard Eulerian and Lagrangian formulations are used to describe both gaseous and liquid phases, respectively. The computing method includes a two-way coupling in which phase properties and spray source terms are interchanging between the two phases within each coupling time step. The fuel used is liquid jet-A1 which is injected in the form of a polydisperse spray and the droplet evaporation rate is calculated using the infinite conductivity model. One-component (n-decane) and two-component fuels (n-decane+toluene) are used as jet-A1 surrogates. The combustion model is based on the mean mixture fraction and its variance, and a presumed-probability density function is used to model turbulent-chemistry interactions. The instantaneous thermochemical state necessary for the chemistry tabulation is determined by using initially the equilibrium (EQ) assumption and thereafter, detailed NEQ calculations through the steady flamelets concept. The combustion chemistry of these surrogates is represented through a reduced chemical kinetic mechanism (CKM) comprising 1,045 reactions among 139 species, derived from the detailed jet-A1 surrogate model, JetSurf 2.0 using a sensitivity based method, Alternate Species Elimination.

Numerical results of the gas velocity, the gas temperature and the species molar fractions are compared with their experimental counterparts obtained from a steady state flame available in the literature. It is observed that, SAS coupled to the tabulated flamelet-based chemistry, predicts reasonably the main flame trends, while URANS even provided with the same combustion model and computing resources, leads to a poor prediction of the global flame trends, emphasizing the asset of a proper resolution when simulating spray flames.

The steady flamelet model even coupled with a robust turbulence model does not reproduce accurately the trend of species with slow oxidation kinetics such as CO and H2, because of the restrictiveness of the solutions space of flamelet equations and the assumption of unity Lewis for all species.

This work is adding a contribution for spray flame modeling and can be seen as an extension to the significant efforts for the modeling of gaseous flames using robust turbulence models coupled with the tabulated flamelet-based chemistry approach to considerably reduce computing cost. The exclusive use of a commercial CFD code widely used in the industry allows a direct application of this simulation approach to industrial configurations while keeping computing cost reasonable.

This study is useful to engineers interested in designing combustors of gas turbines and others combustion systems fed with liquid fuels.

The purpose of this paper is to review the applications of the chemical reactor network (CRN) approach for modeling the combustion in gas turbine combustors and classify the CRN construction methods that have been frequently used by researchers.

This paper initiates with introducing the CRN approach as a practical tool for precisely predicting the species concentrations in the combustion process with lower computational costs. The structure of the CRN and its elements as the ideal reactors are reviewed in recent studies. Flow field modeling has been identified as the most important input for constructing the CRNs; thus, the flow field modeling methods have been extensively reviewed in previous studies. Network approach, component modeling approach and computational fluid dynamics (CFD), as the main flow field modeling methods, are investigated with a focus on the CRN applications. Then, the CRN construction approaches are reviewed and categorized based on extracting the flow field required data. Finally, the most used kinetics and CRN solvers are reviewed and reported in this paper.

It is concluded that the CRN approach can be a useful tool in the entire process of combustion chamber design. One-dimensional and quasi-dimensional methods of flow field modeling are used in the construction of the simple CRNs without detailed geometry data. This approach requires fewer requirements and is used in the initial combustor designing process. In recent years, using the CFD approach in the construction of CRNs has been increased. The flow field results of the CFD codes processed to create the homogeneous regions based on construction criteria. Over the past years, several practical algorithms have been proposed to automatically extract reactor networks from CFD results. These algorithms have been developed to identify homogeneous regions with a high resolution based on the splitting criteria.

This paper reviews the various flow modeling methods used in the construction of the CRNs, along with an overview of the studies carried out in this field. Also, the usual approaches for creating a CRN and the most significant achievements in this field are addressed in detail.

This paper aims to assess the sensitivity of numerical simulation results of turbulent reactive flow to the formulation of inlet boundary conditions. The analysis concerns the profiles of the mean velocity the turbulence kinetic energy k and its dissipation rate ϵ. It is intended to provide guidance to the determination of inlet conditions when only global flow data are available. This situation can be met both in simple laboratory experiments and in industrial full-scale applications, when measurements are either incomplete or infeasible, resulting in lack of detailed inlet data.

Two turbulence–chemistry interaction models were studied: eddy dissipation concept and partially stirred reactor. Three different velocity profiles and related turbulence statistics were applied to present feasible scenarios and their consequences. Simulations with the most appropriate inlet data were accompanied with profiles of turbulent quantities obtained with a proposed method. This method was contrasted to other approaches popular in the literature: the pre-inlet pipe and the separate cold flow simulations of a burner. The methodology was validated on two laboratory-scale jet flames: Delft Jet-in-Hot-Coflow and Sandia CHN B. The simulations were carried out with open source code OpenFOAM.

The proposed relations for turbulence kinetic energy and its dissipation rate at the inlet are found to provide results comparable to those obtained with the use of experimental data as inlet boundary conditions. Moreover, from a certain location downstream the jet, weakly dependent on the Reynolds number, the influence of inlet conditions on flow statistics was found to be negligible.

This work reveals the consequences of the use of rather crude assumptions made for inlet boundary conditions. Proposed formulas for the profiles for k and epsilon are attractive alternatives to other approaches aiming to determine the inlet boundary conditions for turbulent jet flows.