Computing statistical Euler limits of the Navier--Stokes equations in three dimensions
arXiv:2608.23786
Abstract
We develop a Monte Carlo lattice Boltzmann method to compute statistical solutions to the three-dimensional incompressible Navier--Stokes and Euler equations. Entropic space-time adaptive relaxation of the higher order kinetic moments yields stable numerical solutions with decreasing viscosity. We provide a convergence analysis that is conditional on four explicitly stated assumptions regarding the discrete dynamics. Under diffusive scaling, the laws of the discrete ensemble converge along a subsequence to a limit satisfying the Foias--Temam Liouville formulation of the Navier--Stokes equations. Consequently, provided the structure function scaling holds uniformly, the vanishing viscosity limit of these measures satisfies the multi-point statistical Euler hierarchy of Fjordholm, Mishra, and Weber. The limit measures inherit the known weak-strong uniqueness principle on the interval of existence of a strong Euler solution. Under explicit scaling assumptions, a Kuznetsov-type argument yields a fractional 1-Wasserstein convergence rate. We present three-dimensional computations of time-dependent statistical solutions along the inviscid limit of the incompressible Navier--Stokes equations together with convergence measurements in the Wasserstein metric. Numerical experiments on a randomized Taylor--Green vortex with 24-dimensional initial uncertainty recover Kolmogorov's K41 scaling for energy spectra and structure functions, exhibit the failure of pathwise strong convergence, and yield Wasserstein convergence rates of about at the onset of turbulence. Finally, error measurements with respect to spectral hyperviscosity computations indicate that the computed limit measure is independent of the numerical regularization.