Efficient supersonic flows through high-order guided equilibrium with lattice Boltzmann

TL;DR

This paper introduces a high-order lattice Boltzmann method with an extended equilibrium state for efficient simulation of supersonic, thermal, and compressible flows in 3D, demonstrating stability and accuracy in complex scenarios.

Contribution

The paper develops a novel DDF-LBM with 13-moment equilibrium and analytical relaxation times, enabling stable, efficient simulation of high-speed compressible flows with fewer velocities.

Findings

Successfully simulates 1D Riemann problem with discontinuities.

Accurately models flow past a NACA0012 airfoil at Mach 1.5.

Outperforms previous models in efficiency and stability.

Abstract

A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on an extended equilibrium state that is constructed to reproduce the first 13 moments of the Maxwell-Boltzmann distribution exactly. This extends the validity of the standard 5-constraint (mass, momentum and energy) approach and leads to the correct simulation of thermal, compressible flows with only 39 discrete velocities in 3D. The stability of this BGK-LBM is reinforced by relying on Knudsen-number-dependent relaxation times that are computed analytically. Hence, high-Reynolds number, supersonic flows can be simulated in an efficient and elegant manner. While the 1D Riemann problem shows the ability of the proposed approach to handle discontinuities in the zero-viscosity limit, the simulation of the flow past a NACA0012 airfoil (Mach number $\mathrm{Ma}=1.5$, Reynolds number $\mathrm{Re=10^4}$) confirms the excellent behavior of this model in a low-viscosity and supersonic regime. The proposed model is substantially more efficient than the previous 5-moment D3Q343 DDF-LBM and opens up a whole new world of compressible flow applications that can be realistically tackled with a purely LB approach.