🤖 AI Summary
This work addresses the lack of uncertainty quantification in numerical solutions of the incompressible Navier–Stokes equations by proposing a Bayesian sequential inference framework. The approach models discretized dynamics as a state-space system, integrates the Feynman–Kac stochastic representation with spectral methods, and incorporates a particle learning mechanism to mitigate weight degeneracy. For the first time, this framework enables a deep integration of Navier–Stokes solvers with Bayesian inference, supporting non-Gaussian heavy-tailed error modeling, sequential data assimilation under partial observations, and scale-augmented latent variables. Demonstrated on two- and three-dimensional test cases, the method robustly propagates uncertainty through time, significantly enhances resilience to model error, and yields full posterior distributions over quantities of interest rather than point estimates.
📝 Abstract
We develop a Bayesian methodology for numerical solution of the incompressible Navier--Stokes equations with quantified uncertainty. The central idea is to treat discretized Navier--Stokes dynamics as a state-space model and to view numerical solution as posterior computation: priors encode physical structure and modeling error, and the solver outputs a distribution over states and quantities of interest rather than a single trajectory. In two dimensions, stochastic representations (Feynman--Kac and stochastic characteristics for linear advection--diffusion with prescribed drift) motivate Monte Carlo solvers and provide intuition for uncertainty propagation. In three dimensions, we formulate stochastic Navier--Stokes models and describe particle-based and ensemble-based Bayesian workflows for uncertainty propagation in spectral discretizations. A key computational advantage is that parameter learning can be performed stably via particle learning: marginalization and resample--propagate (one-step smoothing) constructions avoid the weight-collapse that plagues naive sequential importance sampling on static parameters. When partial observations are available, the same machinery supports sequential observational updating as an additional capability. We also discuss non-Gaussian (heavy-tailed) error models based on normal variance-mean mixtures, which yield conditionally Gaussian updates via latent scale augmentation.