Search results

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.

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).

This study aims to numerically investigate the flow features and mixing/combustion efficiencies in a turbulent reacting jet in cross-flow by a hybrid Eulerian-Lagrangian methodology.

A high-order hybrid solver is employed where, the velocity field is obtained by solving the Eulerian filtered compressible transport equations while the species are simulated by using the filtered mass density function (FMDF) method.

The main features of a reacting JICF flame are reproduced by the large-eddy simulation (LES)/FMDF method. The computed mean and root-mean-square values of velocity and mean temperature field are in good agreement with experimental data. Reacting JICF’s with different momentum ratios are considered. The jet penetrates deeper for higher momentum ratios. Mixing and combustion efficiency are improved by increasing the momentum ratio.

The authors investigate the flow and combustion characteristics in subsonic reacting JICFs for which very limited studies are reported in the literature.

The influence of computational aspects on simulation results is quantitatively investigated for the specific case of a turbulent piloted jet diffusion flame (Sandia Flame D). It is illustrated that, with a fixed turbulence and chemistry model, the results can heavily depend on the numerical aspects. The influence of inlet boundary conditions has already been discussed in an earlier paper. In this work, attention is focused onto the order of accuracy of the spatial discretization in the numerical scheme and onto the position of the outlet boundary. It is stressed that the purpose is not to judge the quality of the applied models, but to illustrate the possible impact of numerical influence factors. The conclusion is a warning message and a demonstration that all numerical aspects must be completely described when calculation results are presented.

The dynamics of a transitional propane jet diffusion flame with a fuel‐jet velocity of 2.2 m/s has been studied using an implicit, third‐order‐accurate, upwind numerical scheme. The large‐scale vortices outside the flame surface and the small‐scale ones inside were simulated simultaneously, and their interactions with the flame surface were investigated. Numerical experiments were conducted to gain insight into the influence of buoyancy and shear‐layer forcing on the development of the outer and inner vortices. In the presence of buoyancy forces, the outer vortices developed as part of the solution, and the vortex‐crossing frequency was approximately 15 Hz. The inner structures were manifested from a weak perturbation the vorticity that was introduced at the nozzle exit, and, at 185 Hz, these vortices were found to travel farther downstream. It was also found that the inner vortices do not play a role in the formation of the outer vortices, and vice versa. However, the growth of the inner vortices in the downstream locations is strongly influenced by the slowly moving outer vortices.

The purpose of this numerical work is to present and test a new approach for large-scale scalar advection (splicing) in large eddy simulations (LES) that use the linear eddy sub-grid mixing model (LEM) called the LES-LEM.

The new splicing strategy is based on an ordered flux of spliced LEM segments. The principle is that low-flux segments have less momentum than high-flux segments and, therefore, are displaced less than high-flux segments. This strategy affects the order of both inflowing and outflowing LEM segments of an LES cell. The new splicing approach is implemented in a pressure-based fluid solver and tested by simulation of passive scalar transport in a co-flowing turbulent rectangular jet, instead of combustion simulation, to perform an isolated investigation of splicing. Comparison of the new splicing with a previous splicing approach is also done.

The simulation results show that the velocity statistics and passive scalar mixing are correctly predicted using the new splicing approach for the LES-LEM. It is argued that modeling of large-scale advection in the LES-LEM via splicing is reasonable, and the new splicing approach potentially captures the physics better than the old approach. The standard LES sub-grid mixing models do not represent turbulent mixing in a proper way because they do not adequately represent molecular diffusion processes and counter gradient effects. Scalar mixing in turbulent flow consists of two different processes, i.e. turbulent mixing that increases the interface between unmixed species and molecular diffusion. It is crucial to model these two processes individually at their respective time scales. The LEM explicitly includes both of these processes and has been used successfully as a sub-grid scalar mixing model (McMurtry et al., 1992; Sone and Menon, 2003). Here, the turbulent mixing capabilities of the LES-LEM with a modified splicing treatment are examined.

The splicing strategy proposed for the LES-LEM is original and has not been investigated before. Also, it is the first LES-LEM implementation using unstructured grids.

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 study aims to methodologically investigate heat transfer effects on reacting flow inside a liquid-fueled, swirl-stabilized burner. Furthermore, particular attention is paid to turbulence modeling and the results of Reynolds-averaged Navier–Stokes and large eddy simulation approaches are compared in terms of velocity field and flame temperature.

Simulations consist liquid fuel distribution using Eulerian–Lagrangian approach. Flamelet-Generated Manifold combustion model, which is a mixture fraction-progress variable formulation, is used to obtain reacting flow field. Discrete ordinates method is also added for modeling radiation heat transfer effect inside the burner. As a parametric study, different thermal boundary conditions namely: adiabatic wall, constant temperature and heat transfer coefficient are applied. Because of the fact that the burner is designed for operating with different materials, the effects of burner material on heat transfer and combustion processes are investigated. Additionally, material temperatures have been calculated using 1 D method. Finally, soot particles, which are source of luminous radiation in gas turbine combustors, are calculated using Moss-Brookes model.

The results show that the flow behavior is obviously different in recirculation region for both turbulence modeling approach, and this difference causes change on flame temperature distribution, particularly in the outer recirculation zone and region close to swirler. In thermal assessment of the burner, it is predicted that material of the burner walls and the applied thermal boundary conditions have significant influence on flame temperature, wall temperature and flow field. The radiation heat transfer also makes a strong impact on combustion inside the burner; however, luminous radiation arising from soot particles is negligible for the current case.

These types of burners are widely used in research of gas turbine combustion, and it can be seen that the heat transfer effects are generally neglected or oversimplified in the literature. This parametric study provides a basic understanding and methodology of the heat transfer effects on combustion to the researchers.

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.

A three‐dimensional computer simulation of a combustion chamber used in the glass production industry is presented. A numerical solution technique is used to solve the governing time‐averaged partial differential equation and the physical modelling for turbulence, combustion and thermal radiation. A two‐equation turbulence model is employed along with a combustion model based on a fast kinetics statistical approach. A radiation model is used along with the Hottel mixed grey gas model. To solve the governing differential equations an implicit technique of finite‐difference kind is applied. The economy of the computations is very considerably enhanced by the separate calculation of the burner and bulk glass combustion chamber regions, in a manner which takes account of the differing physical nature of their flows. The burner outlet region is calculated with an axisymmetric model. Such two‐dimensional calculations allowed a good resolution of the burner outlet, and provide the inlet conditions for the three‐dimensional calculations of the glass furnace. The prediction procedure is applied to an industrial glass furnace, which operates with oxy‐fuel conditions. Measurements of mean gas temperature and concentrations were performed at different locations in the furnace. The calculated flame length, temperature field and concentrations are with satisfactory agreement with the measured ones.