Saved in:
Bibliographic Details
Main Authors: Zhornyak, Lyra, Forgoston, Eric, Hsieh, M. Ani
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2605.07157
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866914542962868224
author Zhornyak, Lyra
Forgoston, Eric
Hsieh, M. Ani
author_facet Zhornyak, Lyra
Forgoston, Eric
Hsieh, M. Ani
contents We present the first method to directly use a learned continuous Lagrangian to forecast the dynamics of systems governed by partial differential equations, exploiting the inherent conservative structure to achieve stable long-range predictions. We develop an optimization-based integrator that minimizes the squared Euler--Lagrange residual via a mesh-free near-symplectic construction on local space-time patches. Different from integrators for analytical models, integrators for learned models should decouple model error (phase error) from integration error (conservation error). By relying on optimization rather than time-stepping, we bypass the global coupling inherent to fixed discretizations, which slows time- and space-stepping and complicates learning. Our method scales linearly with domain size via Jacobi iteration, and places no structural requirements on the learned network, allowing it to be coupled with existing physics-guided machine learning (ML) methods. We validate our approach on a learned representation of a double pendulum, a one-dimensional wave equation, and a two-dimensional wave equation. Our method achieves error comparable to classical symplectic methods while generalizing to spatially varying dynamics and arbitrary boundary conditions without retraining.
format Preprint
id arxiv_https___arxiv_org_abs_2605_07157
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Learned Lagrangian Models of PDEs via Euler-Lagrange Residual Minimization
Zhornyak, Lyra
Forgoston, Eric
Hsieh, M. Ani
Machine Learning
We present the first method to directly use a learned continuous Lagrangian to forecast the dynamics of systems governed by partial differential equations, exploiting the inherent conservative structure to achieve stable long-range predictions. We develop an optimization-based integrator that minimizes the squared Euler--Lagrange residual via a mesh-free near-symplectic construction on local space-time patches. Different from integrators for analytical models, integrators for learned models should decouple model error (phase error) from integration error (conservation error). By relying on optimization rather than time-stepping, we bypass the global coupling inherent to fixed discretizations, which slows time- and space-stepping and complicates learning. Our method scales linearly with domain size via Jacobi iteration, and places no structural requirements on the learned network, allowing it to be coupled with existing physics-guided machine learning (ML) methods. We validate our approach on a learned representation of a double pendulum, a one-dimensional wave equation, and a two-dimensional wave equation. Our method achieves error comparable to classical symplectic methods while generalizing to spatially varying dynamics and arbitrary boundary conditions without retraining.
title Learned Lagrangian Models of PDEs via Euler-Lagrange Residual Minimization
topic Machine Learning
url https://arxiv.org/abs/2605.07157