Convergence and error analysis of PINNs
Nathan Doumèche, Gérard Biau, Claire Boyer
TL;DR
This work analyzes the theoretical behavior of physics-informed neural networks (PINNs) in two regimes: hybrid modeling and PDE solvers. It first shows that vanilla PINNs can overfit and fail risk-consistency, motivating ridge regularization to achieve risk-consistency for linear and nonlinear PDEs. It then develops a Sobolev-regularization framework that yields well-posed weak formulations and, for linear PDEs, strong convergence results, including PDE-solver and hybrid-modeling cases with quantified physics-consistency. The results provide practical guidance on regularization choices, architecture considerations, and performance guarantees in PINN-based PDE solving and physics-informed learning. Overall, the paper links functional-analytic regularization to improved convergence and physics fidelity in PINNs, supporting reliable extrapolation and robust solution reconstruction.
Abstract
Physics-informed neural networks (PINNs) are a promising approach that combines the power of neural networks with the interpretability of physical modeling. PINNs have shown good practical performance in solving partial differential equations (PDEs) and in hybrid modeling scenarios, where physical models enhance data-driven approaches. However, it is essential to establish their theoretical properties in order to fully understand their capabilities and limitations. In this study, we highlight that classical training of PINNs can suffer from systematic overfitting. This problem can be addressed by adding a ridge regularization to the empirical risk, which ensures that the resulting estimator is risk-consistent for both linear and nonlinear PDE systems. However, the strong convergence of PINNs to a solution satisfying the physical constraints requires a more involved analysis using tools from functional analysis and calculus of variations. In particular, for linear PDE systems, an implementable Sobolev-type regularization allows to reconstruct a solution that not only achieves statistical accuracy but also maintains consistency with the underlying physics.
