Numerical method for strongly variable-density flows at low Mach number: flame-sheet regularisation and a mass-flux immersed boundary method

TL;DR

This paper presents a novel numerical approach for simulating low Mach number flows with strong temperature gradients, incorporating flame-sheet regularisation and an extended immersed boundary method for complex geometries.

Contribution

It introduces a simplified low-Mach-number flow model with regularisation and an extended immersed boundary method to handle mass flux, suitable for combustion and complex geometries.

Findings

Method accurately captures flame front discontinuities.

Numerical scheme is stable and robust in test cases.

Effective handling of mass flux across complex boundaries.

Abstract

A low-Mach-number flow, in the laminar regime, has intrinsically two characteristic spatial scales for a given time scale, or two characteristic temporal scales for a given spatial scale, and these dual scales are very different due to the disparity between the flow and acoustic speed. Therefore low-Mach-number flows impose mathematical and computational challenges in their description. Standard numerical methods for compressible flows, which are typically designed for problems with a single dominant spatial and temporal scale, require alternative approaches such as preconditioning techniques or solvers tailored for low-Mach-number equations. The present work introduces a simplified fluid dynamics model for flows at low Mach number, based on the fractional time-step method. The proposed approach is suitable for handling strong temperature gradients and thermal diffusion, as encountered in combustion systems. To address discontinuities at the flame front in reacting-flow cases, due to the hypothesis of infinitely fast chemistry, a regularisation procedure is employed. Additionally, the immersed boundary method (IBM) is extended to handle mass flux across the boundary surface, enabling simulations of fuel ejection from an arbitrary burner geometry, using a convenient Cartesian grid. The numerical method utilises a predictor-corrector scheme for time integration on a collocated grid, with flux interpolation to prevent numerical pressure oscillations (``odd-even decoupling''). Relevant test cases are used to verify the methods and their implementations, demonstrating correctness and robustness.