Computational Fluid Dynamics Study of a Nonpremixed Turbulent Flame Using openfoam: Effect of Chemical Mechanisms and Turbulence Models

2021 ◽  
Vol 143 (11) ◽  
Author(s):  
Hermann Chopkap Noume ◽  
Valentin Bomba ◽  
Marcel Obounou ◽  
Henri Ekobena Fouda ◽  
Flavian Emmanuel Sapnken

Abstract This article presents a study of the influence of chemical mechanisms and turbulence models on Reynolds-averaged Navier–Stokes (RANS) simulations of the CH4/H2/N2-air turbulent diffusion flame, i.e., the so-called DLR-A flame. The first part of this study is focused on the assessment of the influence of four chemical models on predicted profiles of the DLR-A flame. The chemical mechanisms considered are as follows: (i) a C2 compact skeletal mechanism, which is derived from the GRI3.0 mechanism using an improved multistage reduction method, (ii) a C1 skeletal mechanism containing 41 elementary reactions amongst 16 species, (iii) the global mechanism by Jones and Lindstedt, (iv) and a global scheme consisting of the overall reactions of methane and dihydrogen. RANS numerical results (e.g., velocities, temperature, species, or the heat production rate profiles) obtained running the reactingFOAM solver with the four chemical mechanisms as well as the standard k − ɛ model, the partially stirred reactor (PaSR) combustion model, and the P − 1 radiation model indicate that the C2 skeletal mechanism yields the best agreement with measurements. In the second part of this study, four turbulence models, namely, the standard k − ɛ model, the renormalization group (RNG) k − ɛ model, realizable k − ɛ model, and the k − ω shear stress transport (SST) model, are considered to evaluate their effects on the DLR-A flame simulation results obtained with the C2 skeletal mechanism. Results reveal that the predictions obtained with the standard k − ɛ and the RNG k − ɛ models are in very good agreement with the experimental data. Hence, for simple jet flame with moderately high Reynolds number such as the DLR-A flame, the standard k-epsilon can model the turbulence with a very good accuracy.

2019 ◽  
Vol 142 (3) ◽  
Author(s):  
Hermann Chopkap Noume ◽  
Valentin Bomba ◽  
Marcel Obounou

Abstract The present work assesses the capabilities of a compact skeletal mechanism, derived using an in-house reduction code, to accurately model chemical processes in a turbulent CH4/H2/N2 flame. To this end, a numerical investigation of the DLR-A flame is performed using the free and open-source code openfoam with the derived mechanism. Specifically, the numerical investigation is performed using the Reynolds-averaged Navier–Stokes (RANS) approach and a compact skeletal mechanism consisting of 51 elementary reactions among 21 species. The skeletal mechanism is derived from the GRI3.0 mechanism using an improved multistage reduction method. The k − ɛ model is used as a closure for the RANS equations, while the source terms in the species and energy transport equations are closed by the partially stirred reactor (PaSR) model. The radiation term is modeled by the P-1 model. The numerical results show a good agreement with the experimental data.


Author(s):  
Usman Allauddin ◽  
Michael Pfitzner

Recently, a fractal-based algebraic flame surface density (FSD) premixed combustion model has been derived and validated in the context of large eddy simulation (LES). The fractal parameters in the model, namely the cut-off scales and the fractal dimension were derived using theoretical models, experimental and direct numerical simulation (DNS) databases. The model showed good performance in predicting the premixed turbulent flame propagation for low to high Reynold numbers (Re) in ambient as well as elevated pressure conditions. Several LES combustion models have a direct counterpart in the Reynolds-averaged Navier–Stokes (RANS) context. In this work, a RANS version of the aforementioned LES subgrid scale FSD combustion model is developed. The performance of the RANS model is compared with that of the original LES model and validated with the experimental data. It is found that the RANS version of the model shows similarly good agreement with the experimental data.


Author(s):  
N. Parameswara Rao ◽  
K. Arul Prakash

Numerical simulation of complex three-dimensional flow through the spiral casing has been studied using a finite element method. An explicit Eulerian velocity correction scheme has been employed to solve the Reynolds averaged Navier-stokes equations. The simulation has been performed to describe the flow in high Reynolds number (106) regime and two k-ε turbulence models (standard k-ε and RNG k-ε) have been used for computing the turbulent flow. A streamline upwind Petrov Galerkin technique has been used for spatial discretisation. The velocity field and the pressure distribution inside the spiral casing has been studied. It has been observed that very strong secondary flow is evolved on the cross-stream planes.


2002 ◽  
Vol 124 (4) ◽  
pp. 924-932 ◽  
Author(s):  
Scott Morton ◽  
James Forsythe ◽  
Anthony Mitchell ◽  
David Hajek

An understanding of vortical structures and vortex breakdown is essential for the development of highly maneuverable vehicles and high angle of attack flight. This is primarily due to the physical limits these phenomena impose on aircraft and missiles at extreme flight conditions. Demands for more maneuverable air vehicles have pushed the limits of current CFD methods in the high Reynolds number regime. Simulation methods must be able to accurately describe the unsteady, vortical flowfields associated with fighter aircraft at Reynolds numbers more representative of full-scale vehicles. It is the goal of this paper to demonstrate the ability of detached-eddy Simulation (DES), a hybrid Reynolds-averaged Navier-Stokes (RANS)/large-eddy Simulation (LES) method, to accurately predict vortex breakdown at Reynolds numbers above 1×106. Detailed experiments performed at Onera are used to compare simulations utilizing both RANS and DES turbulence models.


1997 ◽  
Vol 119 (4) ◽  
pp. 900-905 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

In this paper, computational results are presented for three-dimensional high-Reynolds number turbulent flows over a simplified submarine model. The simulation is based on the solution of Reynolds-Averaged Navier-Stokes equations and two-equation turbulence models by using a preconditioned time-stepping approach. A multiblock method, in which the block loop is placed in the inner cycle of a multi-grid algorithm, is used to obtain versatility and efficiency. It was found that the calculated body drag, lift, side force coefficients and moments at various angles of attack or angles of drift are in excellent agreement with experimental data. Fast convergence has been achieved for all the cases with large angles of attack and with modest drift angles.


2020 ◽  
Vol 14 (2) ◽  
pp. 6663-6678
Author(s):  
Akshay Sherikar ◽  
P. J. Disimile

The objective of this study is to expound on the deliverables of a steady-state RANS (Reynolds Averaged Navier Stokes) simulation in one of the simplest flows, Couette flow, at a very high Reynolds number. To that end, a process to perform better grid sensitivity testing is introduced. Three two-equation turbulence models ( , , and ) are compared against each other as well as pitted against formal literature on the subject and core flow velocities, slopes, wall-bounded velocities, shear stresses and kinetic energies are analyzed.  applied with enhanced wall functions is consistently found to be in better agreement with previous studies. Finally, plane turbulent Couette flow at  51,099, the range at which it has not been studied experimentally, numerically or analytically in former studies, is simulated. The results are found to be consistent with the trends asserted by literature and preliminary computations of this study.


Author(s):  
Kazui Fukumoto ◽  
Yoshifumi Ogami

This research aims at developing a turbulent diffusion combustion model based on the chemical equilibrium method and chemical kinetics for simplifying complex chemical mechanisms. This paper presents a combustion model based on the chemical equilibrium method and the eddy dissipation concept (CE-EDC model); the CE-EDC model is validated by simulating a H2-air turbulent diffusion flame. In this model, the reaction rate of fuels and intermediate species is estimated by using the equations of the EDC model. Further, the reacted fuels and intermediate species are assumed to be in chemical equilibrium; the amount of the other species is determined from the amount of the reacted fuels, intermediate species, and air as reactants by using the Gibbs free energy minimization method. An advantage of the CE-EDC model is that the amount of the combustion products can be determined without using detailed chemical mechanisms. The results obtained by using this model were in good agreement with the experimental and computational data obtained by using the EDC model. Using this model, the amount of combustion products can be calculated without using detailed chemical mechanisms. Further, the accuracy of this model is same as that of the EDC model.


Author(s):  
Mehran Masoumifar ◽  
Suyash Verma ◽  
Arman Hemmati

Abstract This study evaluates how Reynolds-Averaged-Navier-Stokes (RANS) models perform in simulating the characteristics of mean three-dimensional perturbed flows in pipes with targeted wall-shapes. Capturing such flow features using turbulence models is still challenging at high Reynolds numbers. The principal objective of this investigation is to evaluate which of the well-established RANS models can best predict the flow response and recovery characteristics in perturbed pipes at moderate and high Reynolds numbers (10000-158000). First, the flow profiles at various axial locations are compared between simulations and experiments. This is followed by assessing the well-known mean pipeflow scaling relations. The good agreement between our computationally predicted data using Standard k-epsilon model and those of experiments indicated that this model can accurately capture the pipeflow characteristics in response to introduced perturbation with smooth sinusoidal axial variations.


Author(s):  
Yiannis Constantinides ◽  
Owen H. Oakley

The prediction of deepwater riser Vortex Induced Vibrations (VIV) is one of the most challenging areas in the offshore industry. Numerous experimental and numerical studies have been performed in an effort to improve the understanding and prediction of cylinder VIV behavior. This paper presents the numerical simulation of rigid circular sections, both bare and fitted with strakes, using a second order accurate finite element computational fluid dynamics (CFD) method. Two turbulence models are examined: the Spalart-Allmaras Reynolds Averaged Navier Stokes (RANS) and the Detached Eddy Simulation (DES). Pragmatic high Reynolds number simulations of fixed and moving cylinders are presented and compared with laboratory experiments. Flow visualization provides insights on how strakes mitigate VIV. Comparisons between RANS and DES results are also presented and discussed.


2014 ◽  
Vol 757 ◽  
Author(s):  
Stephen B. Pope

AbstractFor inhomogeneous turbulent flows at high Reynolds number, it is shown that the redistribution term in Reynolds-stress turbulence models can be determined from the velocity–acceleration correlation. It is further shown that the drift coefficient in the generalized Langevin model (which is used in probability density function (PDF) methods) can be determined from the Reynolds stresses and the velocity–acceleration correlation. These observations are valuable, since the second moments of velocity and acceleration can be measured in experiments, in direct numerical simulations and in well-resolved large-eddy simulations (LES), and hence these turbulence-model quantities can be determined. The redistribution is closely related to the pressure–rate-of-strain, and the unknown in the PDF equation is closely related to the conditional mean pressure gradient (conditional on velocity). In contrast to the velocity–acceleration moments, these pressure statistics are much more difficult to obtain, and our knowledge of them is quite limited. It is also shown that the generalized Langevin model can be re-expressed to provide a direct connection between the drift term and the fluid acceleration. All of these results are first obtained using the constant-property Navier–Stokes equations, but it is then shown that the results are simply extended to variable-density flows.


Sign in / Sign up

Export Citation Format

Share Document