Skip to content
DraconiticResearch
Search Sources API
203,234 studies31 sourcesUpdated 49 min ago
PreprintCC BY 4.0arXiv

Warm-starting PDE solvers with any-dimensional machine learning

Abstract

Any-dimensional machine learning models, such as graph neural networks (GNNs), can be naturally trained and evaluated on inputs of different sizes and dimensions. Inspired by the GNN transferability literature, we show mathematical conditions under which a partial differential equation (PDE) learning-based solver can be trained in small dimensions and directly applied to solve a higher dimensional PDE in a zero-shot fashion. These conditions are based on symmetries in both the partial differential equation and the initial data. When the equations satisfy the symmetries but the data does not, which is the case for many PDEs arising from physics, we show that our theory gives a principled way of warm-starting low-dimensional PDE solvers for higher dimensional PDEs. We apply this method on the heat equation, Burgers' equation, and the compressible Navier--Stokes equations, improving the performance in both zero-shot and typical training regimes on high dimensional data. For example, we train a surrogate model on 2D Navier--Stokes data and achieve better results on 3D test data than a baseline surrogate model trained on 3D data, while only using 12% of the flops and 20% of the total data size.

Related

  1. Transformer-Based Warm-Starting for Feasible and Optimal Terminal Approach to Tumbling Objects with Space Manipulators
    Yuji Takubo, Maximilian Adang, Mac Schwager +1
    arXivPreprintarXiv licenceFull textOptimisation & controlarXiv
  2. PTED: A multi-dimensional two-sample test for scientific inference and generative machine learning
    Connor Stone
    arXivPreprintCC BY 4.0Full textInstrumentation & methodsStatisticsarXiv
  3. Deep-Learning Solvers and Surrogates for Infinity and p-Laplace Problems
    Tak Shing Au Yeung, Ka Chun Cheung, Hannah Potgieter +2
    arXivPreprintCC BY 4.0Full textNumerical analysis & scientific computingarXiv
  4. Towards learning-enhanced solvers for multi-parametric optimization
    Dieter Teichrib, Fabian König, Moritz Schulze Darup
    arXivPreprintCC BY-NC-ND 4.0Full textOptimisation & controlarXiv
  5. BandPC: Learning Residual-Band Preconditioner Combinations for Flexible Conjugate Gradient Solvers
    D. M. Li, Shang-Tian Yang, Xin Qiu
    arXivPreprintarXiv licenceFull textNumerical analysis & scientific computingarXiv
BibTeX
@misc{gregory2026warmstarting,
  title = {{Warm-starting PDE solvers with any-dimensional machine learning}},
  author = {Wilson G. Gregory and George A. Kevrekidis and Ben Blum-Smith and Soledad Villar},
  howpublished = {arXiv},
  year = {2026},
  month = sep,
  doi = {10.48550/arxiv.2609.38916},
  eprint = {2609.38916},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2609.38916},
  note = {CC BY 4.0},
}
RIS
TY  - UNPB
TI  - Warm-starting PDE solvers with any-dimensional machine learning
AU  - Wilson G. Gregory
AU  - George A. Kevrekidis
AU  - Ben Blum-Smith
AU  - Soledad Villar
PY  - 2026
DA  - 2026/09/30
PB  - arXiv
DO  - 10.48550/arxiv.2609.38916
UR  - https://arxiv.org/abs/2609.38916
AB  - Any-dimensional machine learning models, such as graph neural networks (GNNs), can be naturally trained and evaluated on inputs of different sizes and dimensions. Inspired by the GNN transferability literature, we show mathematical conditions under which a partial differential equation (PDE) learning-based solver can be trained in small dimensions and directly applied to solve a higher dimensional PDE in a zero-shot fashion. These conditions are based on symmetries in both the partial differential equation and the initial data. When the equations satisfy the symmetries but the data does not, which is the case for many PDEs arising from physics, we show that our theory gives a principled way of warm-starting low-dimensional PDE solvers for higher dimensional PDEs. We apply this method on the heat equation, Burgers' equation, and the compressible Navier--Stokes equations, improving the performance in both zero-shot and typical training regimes on high dimensional data. For example, we train a surrogate model on 2D Navier--Stokes data and achieve better results on 3D test data than a baseline surrogate model trained on 3D data, while only using 12% of the flops and 20% of the total data size.
ER  -