A Lattice Boltzmann Method with Adaptive Relaxation Parameter and Lax--Friedrichs-Type Equilibrium for Hyperbolic Systems

Abstract: In this paper, a lattice Boltzmann method with an adaptive relaxation parameter and a Lax--Friedrichs-type equilibrium is proposed for hyperbolic systems with source terms. The equilibrium distribution recovers the conservative variables and physical fluxes while incorporating dissipation determined by characteristic-speed bounds. To balance accuracy and robustness, the relaxation parameter is selected from a local smoothness indicator based on characteristic projections. In smooth regions, the parameter approaches the low-dissipation limit, retaining second-order accuracy; near discontinuities, it is automatically reduced to introduce localized dissipation and suppress nonphysical oscillations. The stabilization acts directly through the local collision step and preserves the standard collide-and-stream structure, without a posteriori recomputation or interface-based limiting. Maxwell iteration establishes second-order consistency in smooth regions, and a weighted $L<sup>2$-stability estimate is proved for linear hyperbolic systems with periodic boundary conditions under a standard CFL condition. Numerical experiments for scalar advection, Euler, shallow-water, and reactive Euler equations demonstrate the expected accuracy for smooth solutions and robust resolution of challenging one- and two-dimensional discontinuous problems, including wet--dry fronts and reactive discontinuities. A large-scale simulation of a 2D cellular detonation further demonstrates the capability of the method to resolve long-time multidimensional shock--reaction interactions.