Wavelet Diffusion Neural Operator
About
Simulating and controlling physical systems described by partial differential equations (PDEs) are crucial tasks across science and engineering. Recently, diffusion generative models have emerged as a competitive class of methods for these tasks due to their ability to capture long-term dependencies and model high-dimensional states. However, diffusion models typically struggle with handling system states with abrupt changes and generalizing to higher resolutions. In this work, we propose Wavelet Diffusion Neural Operator (WDNO), a novel PDE simulation and control framework that enhances the handling of these complexities. WDNO comprises two key innovations. Firstly, WDNO performs diffusion-based generative modeling in the wavelet domain for the entire trajectory to handle abrupt changes and long-term dependencies effectively. Secondly, to address the issue of poor generalization across different resolutions, which is one of the fundamental tasks in modeling physical systems, we introduce multi-resolution training. We validate WDNO on five physical systems, including 1D advection equation, three challenging physical systems with abrupt changes (1D Burgers' equation, 1D compressible Navier-Stokes equation and 2D incompressible fluid), and a real-world dataset ERA5, which demonstrates superior performance on both simulation and control tasks over state-of-the-art methods, with significant improvements in long-term and detail prediction accuracy. Remarkably, in the challenging context of the 2D high-dimensional and indirect control task aimed at reducing smoke leakage, WDNO reduces the leakage by 78% compared to the second-best baseline. The code can be found at https://github.com/AI4Science-WestlakeU/wdno.git.
Related benchmarks
| Task | Dataset | Result | Rank | |
|---|---|---|---|---|
| PDE System Rollout Forecasting | Active Matter | VRMSE (01:16)4.481 | 14 | |
| PDE System Rollout Forecasting | Rayleigh-Bénard | VRMSE (Rollout 01:32)19.346 | 14 | |
| PDE System Rollout Forecasting | Turbulent Radiative Layer 2D | VRMSE (Rollouts 01-32)1.517 | 13 | |
| Autoregressive rollout | Active Matter | CRPS (1:16)0.218 | 11 | |
| PDE Emulation | Turbulent Radiative Layer 2D | Spectral Coherence RMSE0.27 | 11 | |
| PDE Emulation | Active Matter | Spectral Coherence RMSE0.303 | 11 | |
| Spectral Coherence Estimation | Turbulent Radiative Layer 2D | Low Frequency Spectral Coherence RMSE0.199 | 11 | |
| Spectral Coherence Estimation | Active Matter | Low Frequency Spectral Coherence RMSE0.366 | 11 | |
| Autoregressive rollout | Turbulent Radiative Layer 2D | CRPS (Rollout 01:32)2.46 | 11 | |
| Autoregressive rollout | Rayleigh-Bénard | CRPS (Rollout 01:32)0.051 | 11 |