A staggered pressure correction numerical scheme to compute a travelling reactive interface in a partially premixed mixture
From MaRDI portal
Publication:5020153
Abstract: We address in this paper a model for the simulation of turbulent deflagrations in industrial applications. The flow is governed by the Euler equations for a variable composition mixture and the combustion modelling is based on a phenomenological approach: the flame propagation is represented by the transport of the characteristic function of the burnt zone, where the chemical reaction is complete; outside this zone, the atmosphere remains in its fresh state. Numerically, we approximate this problem by a penalization-like approach, i.e. using a finite conversion rate with a characteristic time tending to zero with the space and time steps. The numerical scheme works on staggered, possibly unstructured, meshes. The time-marching algorithm is of segregated type, and consists in solving in a first step the chemical species mass balances and then, in a second step, mass, momentum and energy balances. For this latter stage of the algorithm, we use a pressure correction technique, and solve a balance equation for the so-called sensible enthalpy instead of the total energy balance, with corrective terms for consistency. The scheme is shown to satisfy the same stability properties as the continuous problem: the chemical species mass fractions are kept in the [0, 1] interval, the density and the sensible internal energy stay positive and the integral over the computational domain of a discrete total energy is conserved. In addition, we show that the scheme is in fact conservative, i.e. that its solution satisfy a conservative discrete total energy balance equation, with space and time discretizations which are unusual but consistent in the Lax-Wendroff sense. Finally, we observe numerically that the penalization procedure converges, i.e. that making the chemical time scale tend to zero allows to converge to the solution of the target (infinitely fast chemistry) continuous problem. Tests also evidence that the scheme accuracy dramatically depends on the discretization of the convection operator in the chemical species mass balances. October 14, 2020.
Recommendations
- Flexible Flame Structure Modelling in a Flame Front Tracking Scheme
- A conservative pressure-correction scheme for transient simulations of reacting flows
- A Conservative Adaptive Method for Flame Propagation
- A robust numerical model for premixed flames with high density ratios based on new pressure correction and IMEX schemes
- scientific article; zbMATH DE number 4197633
Cites work
- scientific article; zbMATH DE number 996053 (Why is no real title available?)
- A formally second-order cell centred scheme for convection-diffusion equations on general grids
- A numerical fluid dynamics calculation method for all flow speeds
- An unconditionally stable staggered pressure correction scheme for the compressible Navier-Stokes equations
- Conforming and nonconforming finite element methods for solving the stationary Stokes equations I
- Contact discontinuity capturing schemes for linear advection and compressible gas dynamics
- How to preserve the mass fractions positivity when computing compressible multi-component flows
- Kinetic energy control in the MAC discretization of compressible Navier-Stokes equations
- Numerical Calculation of Time-Dependent Viscous Incompressible Flow of Fluid with Free Surface
- Numerical approximation of hyperbolic systems of conservation laws
- On some implicit and semi-implicit staggered schemes for the shallow water and Euler equations
- Simple nonconforming quadrilateral Stokes element
- Staggered discretizations, pressure correction schemes and all speed barotropic flows
- The reactive Riemann problem for thermally perfect gases at all combustion regimes
- Turbulent Combustion
This page was built for publication: A staggered pressure correction numerical scheme to compute a travelling reactive interface in a partially premixed mixture
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5020153)