Research

Paper

TESTING March 02, 2026

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

Authors

Matheus P. Severino, Fernando F. Fachini, Elmer M. Gennaro, Daniel Rodríguez, Leandro F. Souza

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.

Metadata

arXiv ID: 2603.01727
Provider: ARXIV
Primary Category: physics.flu-dyn
Published: 2026-03-02
Fetched: 2026-03-03 04:34

Related papers

Raw Data (Debug)
{
  "raw_xml": "<entry>\n    <id>http://arxiv.org/abs/2603.01727v1</id>\n    <title>Numerical method for strongly variable-density flows at low Mach number: flame-sheet regularisation and a mass-flux immersed boundary method</title>\n    <updated>2026-03-02T10:52:04Z</updated>\n    <link href='https://arxiv.org/abs/2603.01727v1' rel='alternate' type='text/html'/>\n    <link href='https://arxiv.org/pdf/2603.01727v1' rel='related' title='pdf' type='application/pdf'/>\n    <summary>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.</summary>\n    <category scheme='http://arxiv.org/schemas/atom' term='physics.flu-dyn'/>\n    <category scheme='http://arxiv.org/schemas/atom' term='physics.comp-ph'/>\n    <published>2026-03-02T10:52:04Z</published>\n    <arxiv:comment>49 pages, 14 figures, 8 tables. This work is part of the first author's Ph.D. thesis at the University of São Paulo (USP)</arxiv:comment>\n    <arxiv:primary_category term='physics.flu-dyn'/>\n    <author>\n      <name>Matheus P. Severino</name>\n    </author>\n    <author>\n      <name>Fernando F. Fachini</name>\n    </author>\n    <author>\n      <name>Elmer M. Gennaro</name>\n    </author>\n    <author>\n      <name>Daniel Rodríguez</name>\n    </author>\n    <author>\n      <name>Leandro F. Souza</name>\n    </author>\n  </entry>"
}