A robust numerical model for premixed flames with high density ratios based on new pressure correction and IMEX schemes

Authors:
F. Paravento
Affiliations:
Laboratory for Aero & Hydrodynamics, Delft University of Technology, The Netherlands
Venue:
Journal of Computational Physics
Year:
2010

Citing 11
Cited 0

Numerical computation of internal & external flows: fundamentals of numerical discretization

Numerical computation of internal & external flows: fundamentals of numerical discretization
High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations

SIAM Journal on Numerical Analysis
Numerical methods for ordinary differential systems: the initial value problem

Numerical methods for ordinary differential systems: the initial value problem
A semi-implicit numerical scheme for reacting flow: I. stiff chemistry

Journal of Computational Physics
An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow

SIAM Journal on Scientific Computing
A PDE-based fast local level set method

Journal of Computational Physics
Algorithm 547: Fortran Routines for Discrete Cubic Spline Interpolation and Smoothing [E1], [E3]

ACM Transactions on Mathematical Software (TOMS)
Weighted ENO Schemes for Hamilton--Jacobi Equations

SIAM Journal on Scientific Computing
Central Differencing Based Numerical Schemes for Hyperbolic Conservation Laws with Relaxation Terms

SIAM Journal on Numerical Analysis
Level Set Evolution without Re-Initialization: A New Variational Formulation

CVPR '05 Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05) - Volume 1 - Volume 01
A solution for the odd-even decoupling problem in pressure-correction algorithms for variable density flows

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

In this study we present a model for the interaction of premixed flames with obstacles in a channel flow. Although the flow equations are solved with Direct Numerical Simulation using a low Mach number approximation, the resolution used in the computation is limited (~1mm) hence the inner structure of the flame and the chemical scales are not solved. The species equations are substituted with a source term in the energy equation that simulates a one-step global reaction. A level set method is applied to track the position of the flame and its zero level is used to activate the source term in the energy equation only at the flame front. An immersed boundary method reproduces the geometry of the obstacles. The main contribution of the paper is represented by the proposed numerical approach: an IMEX (implicit-explicit) Runge-Kutta scheme is used for the time integration of the energy equation and a new pressure correction algorithm is introduced for the time integration of the momentum equations. The approach presented here allows to calculate flames which produce high density ratios between burnt and unburnt regions. The model is verified by simulating first simple solutions for one- and two-dimensional flames. At last, the experiments performed by Masri and Ibrahim with square and rectangular bodies are calculated.