# Physics-constrained robust learning of open-form partial differential equations from limited and noisy data

arXiv: https://arxiv.org/abs/2309.07672 · HTML: https://mlsanity.com/paper/2309.07672/

## Key facts

- Cited by 8 papers in the map (data through 14 Oct 2026).
- Builds on 25 of its 60 references that are in the map.
- First posted on arXiv 14 Sep 2023.

## Builds on (most cited first, 25 of 25)

- [PyTorch: An Imperative Style, High-Performance Deep Learning Library](https://mlsanity.com/paper/1912.01703/) (2019)
- [Scientific Machine Learning through Physics-Informed Neural Networks: Where we are and What's next](https://mlsanity.com/paper/2201.05624/) (2022)
- [Deep Hidden Physics Models: Deep Learning of Nonlinear Partial Differential Equations](https://mlsanity.com/paper/1801.06637/) (2018)
- [PDE-Net 2.0: Learning PDEs from Data with A Numeric-Symbolic Hybrid Deep Network](https://mlsanity.com/paper/1812.04426/) (2018)
- [Physics-informed learning of governing equations from scarce data](https://mlsanity.com/paper/2005.03448/) (2020)
- [Learning Equations for Extrapolation and Control](https://mlsanity.com/paper/1806.07259/) (2018)
- [Weak SINDy For Partial Differential Equations](https://mlsanity.com/paper/2007.02848/) (2020)
- [Ensemble-SINDy: Robust sparse model discovery in the low-data, high-noise limit, with active learning and control](https://mlsanity.com/paper/2111.10992/) (2021)
- [SINDy-PI: A Robust Algorithm for Parallel Implicit Sparse Identification of Nonlinear Dynamics](https://mlsanity.com/paper/2004.02322/) (2020)
- [CAN-PINN: A Fast Physics-Informed Neural Network Based on Coupled-Automatic-Numerical Differentiation Method](https://mlsanity.com/paper/2110.15832/) (2021)
- [DeepMoD: Deep learning for Model Discovery in noisy data](https://mlsanity.com/paper/1904.09406/) (2019)
- [DLGA-PDE: Discovery of PDEs with incomplete candidate library via combination of deep learning and genetic algorithm](https://mlsanity.com/paper/2001.07305/) (2020)
- [Any equation is a forest: Symbolic genetic algorithm for discovering open-form partial differential equations (SGA-PDE)](https://mlsanity.com/paper/2106.11927/) (2021)
- [DL-PDE: Deep-learning based data-driven discovery of partial differential equations from discrete and noisy data](https://mlsanity.com/paper/1908.04463/) (2019)
- [PDE-READ: Human-readable Partial Differential Equation Discovery using Deep Learning](https://mlsanity.com/paper/2111.00998/) (2021)
- [Robust discovery of partial differential equations in complex situations](https://mlsanity.com/paper/2106.00008/) (2021)
- [DySMHO: Data-Driven Discovery of Governing Equations for Dynamical Systems via Moving Horizon Optimization](https://mlsanity.com/paper/2108.00069/) (2021)
- [WeakIdent: Weak formulation for Identifying Differential Equations using Narrow-fit and Trimming](https://mlsanity.com/paper/2211.03134/) (2022)
- [Discovery of partial differential equations from highly noisy and sparse data with physics-informed information criterion](https://mlsanity.com/paper/2208.03322/) (2022)
- [Noise-aware Physics-informed Machine Learning for Robust PDE Discovery](https://mlsanity.com/paper/2206.12901/) (2022)
- [Online Weak-form Sparse Identification of Partial Differential Equations](https://mlsanity.com/paper/2203.03979/) (2022)
- [Bayesian Deep Learning for Partial Differential Equation Parameter Discovery with Sparse and Noisy Data](https://mlsanity.com/paper/2108.04085/) (2021)
- [Physics-Guided Discovery of Highly Nonlinear Parametric Partial Differential Equations](https://mlsanity.com/paper/2106.01078/) (2021)
- [AutoKE: An automatic knowledge embedding framework for scientific machine learning](https://mlsanity.com/paper/2205.05390/) (2022)
- [Parsimony-Enhanced Sparse Bayesian Learning for Robust Discovery of Partial Differential Equations](https://mlsanity.com/paper/2107.07040/) (2021)

## Built on by (most cited first, 8 of 8)

- [LLM4ED: Large Language Models for Automatic Equation Discovery](https://mlsanity.com/paper/2405.07761/) (2024)
- [A Data-Driven Framework for Discovering Fractional Differential Equations in Complex Systems](https://mlsanity.com/paper/2412.03970/) (2024)
- [Data-driven discovery of governing differential equations across physical systems](https://mlsanity.com/paper/2606.09638/) (2026)
- [Discovering an interpretable mathematical expression for a full wind-turbine wake with artificial intelligence enhanced symbolic regression](https://mlsanity.com/paper/2406.00695/) (2024)
- [Evaluating PDE discovery methods for multiscale modeling of biological signals](https://mlsanity.com/paper/2506.20694/) (2025)
- [Symbolic identification of tensor equations in multidimensional physical fields](https://mlsanity.com/paper/2507.01466/) (2025)
- [Efficient temporal prediction of compressible flows in irregular domains using Fourier neural operators](https://mlsanity.com/paper/2601.01922/) (2026)
- [Hypothesize, Evaluate, Refine: A Scientific Agent for PDE Discovery with Unknown Spatial Coefficient Fields](https://mlsanity.com/paper/2608.27475/) (2026)

## Abstract (arXiv)

Unveiling the underlying governing equations of nonlinear dynamic systems remains a significant challenge. Insufficient prior knowledge hinders the determination of an accurate candidate library, while noisy observations lead to imprecise evaluations, which in turn result in redundant function terms or erroneous equations. This study proposes a framework to robustly uncover open-form partial differential equations (PDEs) from limited and noisy data. The framework operates through two alternating update processes: discovering and embedding. The discovering phase employs symbolic representation and a novel reinforcement learning (RL)-guided hybrid PDE generator to efficiently produce diverse open-form PDEs with tree structures. A neural network-based predictive model fits the system response and serves as the reward evaluator for the generated PDEs. PDEs with higher rewards are utilized to iteratively optimize the generator via the RL strategy and the best-performing PDE is selected by a parameter-free stability metric. The embedding phase integrates the initially identified PDE from the discovering process as a physical constraint into the predictive model for robust training. The traversal of PDE trees automates the construction of the computational graph and the embedding process without human intervention. Numerical experiments demonstrate our framework's capability to uncover governing equations from nonlinear dynamic systems with limited and highly noisy data and outperform other physics-informed neural network-based discovery methods. This work opens new potential for exploring real-world systems with limited understanding.
