A macroscopic-shadow-corrected lattice Boltzmann method for fast, time-accurate simulation of low-Reynolds-number transient flows
arXiv:2609.39474
Abstract
Explicit lattice Boltzmann simulations of slow transient flows are constrained by the acoustic time step. Dual-time stepping can remove this restriction, but slow inner convergence has limited reported wall-clock gains to roughly four- to tenfold, without a mechanism that strengthens under refinement. We present the macroscopic-shadow-corrected lattice Boltzmann method (MSC-LBM), which applies defect correction to the unsplit kinetic residual. At each Fourier wavenumber, a Stokes-like system is solved exactly for the conserved residual moments while preserving the second-order backward-differentiation formula (BDF2) fixed point and temporal accuracy. Along the two-relaxation-time axis , makes one collision eliminate all non-hydrodynamic perturbations; vanishes at , leaving hydrodynamic slow modes whose long-wave shadow is inverted by the corrector. In completed timing campaigns, MSC-LBM achieves 27.95- and 52.71-fold wall-clock speedups at , under matched-accuracy and gates. The matched speedup rises monotonically to 122.02 at , with gains persisting across --. The three-dimensional D3Q19 extension reaches 27.85 and 8.84 under the same gates at , . A bounce-back-consistent kinetic coarse solver with adaptive relinearisation extends MSC-LBM to fully enclosed cavities, yielding setup-excluding physical-march speedups of 26.86 at and 3.84 at for . An exact per-wavenumber symbol inverse establishes the attainable off-design contraction envelope. The contraction, parameter-sweep, and timing results jointly delimit the demonstrated operating regime: low-to-moderate-Reynolds-number transients for which acoustic stepping need not dictate computational cost.
arXiv comment5 figures