FlowReconBench: Accuracy, calibration and cost of learned and classical flow reconstruction from sparse sensors
Organizations: Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
Abstract
Reconstructing full flow fields from sparse sensors requires prior assumptions about unmeasured regions and physical variables. Yet it remains unclear when learned models outperform classical estimators and which aspects of their predictions can be trusted. FlowReconBench evaluates ten learned and five classical methods on three regimes: a laminar cylinder wake across Reynolds numbers, two-dimensional Kolmogorov turbulence and a large-scale three-dimensional (3D) isotropic turbulence. Methods share leakage-controlled splits and sensor draws and are scored on per-channel error, ensemble spread, calibration, spectra and cost. The leading method changes for each regime. Classical methods perform strongly in low-order and dense-sensing regimes. On the cylinder wake, linear stochastic estimation and gappy proper orthogonal decomposition lead, including for unobserved pressure. In 3D, inverse-distance weighting matches the best learned methods on observed channels. Data-driven methods become superior at decreasing sensor density or poor sensor observability. On Kolmogorov flow, continuous ranked probability score and error yield different method rankings: generative models offer broader but underdispersed ensembles, while deterministic models prioritize pointwise error. Generative models balance this tradeoff through their sampler budget, with a larger budget improving ensemble quality at higher computational cost. In 3D, scaling tests demonstrate that queryable methods are flat in memory, while full-grid construction methods run out of memory. Observing only 2 of 4 fields leaves the problem kinematically underconstrained, yet a latent model can reproduce structural similarities to the ground truth.
Figures & tables
| Regime | Dataset | Points | Channels | Split | Primary challenge |
|---|---|---|---|---|---|
| 2D laminar periodic | Cylinder wake (new) | cross-Re | Low-rank observability | ||
| 2D chaotic | Kolmogorov [ 65 ] | cross-traj. | Multiscale reconstruction | ||
| 3D chaotic | JHTDB isotropic [ 35 ] | spatial | Unobserved channels |
| Family | Cyl. | Kolm. | JHTDB |
|---|---|---|---|
| Classical (mean, NN, IDW, | ✓ | ✓ | ✓ |
| gappy POD, LSE) | |||
| Point-attention det. (Senseiver) | ✓ | ✓ | ✓ |
| Point-kernel det. (MLP-RBF) | ✓ | ✓ | n.r. |
| Fourier operator det. (Geo-FNO; FNO) | ✓ | ✓ | ✓ |
| Set-DeepONet and Set-DeepONet++ | n.r. | n.r. | ✓ |
| Rel. | ||||||||||
| Method | agg. | CRPS | Spr/Err | Cov. 90 | s/field | |||||
| SiT-point ( ) | 0.073 | 0.038 | 0.108 | 0.057 | 0.105 | 0.057 | 0.034 | 0.51 | 0.40 | 0.50 |
| Senseiver | 0.134 | 0.118 | 0.150 | 0.115 | 0.192 | 0.095 | 0.082 | — | — | 0.06 |
| MLP-RBF | 0.318 | 0.260 | 0.376 | 0.220 | 0.536 | 0.197 | 0.165 | — | — | 0.07 |
| Geo-FNO | 0.079 | 0.058 | 0.100 | 0.063 | 0.109 | 0.064 | 0.045 | — | — | 0.26 |
| LSE | 0.043 | 0.036 | 0.050 | 0.032 | 0.058 | 0.039 | 0.022 | — | — | — |
| Rel. at taps/field | ||||||
| Method | 32 | 64 | 128 | 360 | CRPS | Spr/Err |
| SiT-point | 0.103 | 0.092 | 0.091 | 0.091 | 0.040 | 0.25 |
| Senseiver | 0.130 | 0.106 | 0.101 | 0.099 | 0.063 | — |
| Geo-FNO | 0.236 | 0.182 | 0.173 | 0.157 | 0.086 | — |
| MLP-RBF | 0.370 | 0.322 | 0.303 | 0.290 | 0.156 | — |
| Ridge gappy POD | 0.270 | 0.265 | 0.262 | 0.260 | 0.124 | — |
| Rel. | ||||||||
| Method | mean | sample | CRPS | Spr/Err | Cov. 90 | s/field | GB | |
| SiT-point ( ) | 0.293 | 0.310 | 0.146 | 0.37 | 0.42 | 0.27 | 5.00 | 0.22 |
| Latent FM ( ) | 0.463 | 0.523 | 0.240 | 0.57 | 0.55 | 0.30 | 0.11 | 0.16 |
| S3GM ( ) | 0.273 | 0.363 | 0.110 | 0.95 | 0.75 | 0.93 | 6.80 | 7.30 |
| Senseiver | 0.465 | 0.465 | 0.329 | — | — | 0.023 | 0.07 | 0.81 |
| MLP-RBF | 0.500 | 0.500 | 0.337 | — | — | 0.13 | 0.11 | 1.36 |
| Method | agg. | CRPS | Spr/Err | Cov. 90 | |||||
|---|---|---|---|---|---|---|---|---|---|
| Latent flow matching | 0.249 | 0.710 | 0.260 | 0.643 | 0.465 | 0.239 | 0.69 | 0.57 | 0.056 |
| SiT-point ( ) | 0.175 | 0.957 | 0.182 | 0.904 | 0.555 | 0.323 | 0.46 | 0.40 | 11.7 |
| FNO3D-FM ( ) | 0.190 | 0.921 | 0.193 | 1.054 | 0.589 | 0.272 | 0.98 | 0.63 | 0.012 |
| CoNFiLD (single snapshot, ) a) | 0.500 | 1.024 | 0.543 | 1.034 | 0.775 | 0.378 | 0.77 | 0.62 | — |
| S3GM (normalized guidance) b) | 0.257 | 0.995 | 0.252 | 0.993 | 0.624 | 0.347 | 0.52 | 0.46 | — |
| Gen4Turb c) | 0.212 | 0.916 | 0.325 | 0.876 | 0.582 | 0.307 | 0.70 | 0.54 | — |
| Method | Params | GPU-h | Peak GB | s/field |
|---|---|---|---|---|
| Latent FM | 6.68M | 20.1 | 8.4 | 0.019 |
| SiT-point ( ) | 6.57M | 19.4 | 26.8 | 149 |
| FNO3D-FM ( ) | 6.73M | 27.4 | 48.6 | 0.097 |
| CoNFiLD a) | 118.9M | 19.2 | 10.4 | 36 (d) |
| Gen4Turb | 6.16M | 16.6 | 47.4 | 1.65 |
| S3GM | 6.35M | 18.5 | 26.3 | 73.8 |
| Local (interpolation) | Global (gappy POD) | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| Regime | Sensors | NN | IDW | Floor | Gappy POD | ||||||
| Cylinder, volume | 238 per field | 0.95–1.3 | 0.55 | 0.569 | 0.568 | 20 | 0.997 | 0.048 | 3.2 | 0.057 | |
| Cylinder, surface | 128 taps per field | 13–19 | 2.10 | 2.146 | 2.131 | 5 | 0.929 | 0.163 | 0.594 | ||
| 20 | 0.997 | 0.048 | 0.763 | ||||||||
| Kolmogorov | 655 | 0.40 | 0.58 | 0.580 | 0.530 | 160 | 0.832 | 0.513 | 3.7 | 0.607 | |
| JHTDB, and | 19,531 per channel | 0.19 | 0.216 | 0.179 | 80 | 0.955 a) | 0.648 | 11.9 | 0.953 | ||
Computational Fluid Dynamics
Computational fluid dynamics numerically models fluid motion, turbulence, heat and mass transport, and fluid–structure interactions.
See all 291 papers