A numerical framework for low-speed flows with large thermal variations

• Published in: Computers & fluids Vol. 265; p. 105989
• Main Authors: Kuan, Tzuo Wei It, Szmelter, Joanna
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 30.10.2023
• Subjects: Differentially Heated Cavity; Low Mach number flows; MPDATA; Non-Oberbeck–Boussinesq effects; Non-oscillatory Forward-in-Time schemes; MPDATA; Low Mach number flows; Non-Oberbeck–Boussinesq effects; Non-oscillatory Forward-in-Time schemes; Differentially Heated Cavity
• ISSN: 0045-7930; 1879-0747; 1879-0747
• DOI: 10.1016/j.compfluid.2023.105989

In the numerical simulation of fluid dynamic problems there are situations in which low-speed flows exhibit a non-negligible density variation driven by thermal or compositional distributions. A new method for solving the Low Mach Number equations governing such flows is developed. The basis of the formulation is the semi-implicit scheme build on the Non-oscillatory Forward-in-Time integrator encompassing Multidimensional Positive Definite Advection Transport Algorithm and a robust Krylov solver. The spatial discretisation is constructed using the median dual finite volumes and the method is second-order-accurate in space and time. A collocated arrangement of all flow variables is implemented. The validity of the presented method is evaluated using the Differentially Heated Cavity flows. Three-dimensional simulations demonstrate its efficacy and ability to capture the Non-Oberbeck–Boussinesq effects. Ideally, a numerical method should be able to treat different flow regimes without strong limitations. When implemented, it is shown that in the limit of small density variations the new formulation correctly reproduces simulations of incompressible flows subject to the Oberbeck–Boussinesq approximation. Furthermore, the three-dimensional computations set to the quasi-two-dimensional cases are in excellent agreement with the established two-dimensional Differentially Heated Cavity benchmarks. •The first Low-Mach-Number (LMN) solver from the NFT-MPDATA class is developed.•3D simulations are validated with 2D and 3D Differentially Heated Cavity (DHC) flows.•Validations include DHC flows with low- and high- temperature variations.•Quantitative data for 3D laminar DHC flows exhibiting NOB effects are documented.