arxiv: 2605.14082 · v1 · submitted 2026-05-13 · 🧮 math.NA · cs.NA

Recognition: 2 theorem links

· Lean Theorem

Goal-Oriented Time Adaptivity for Linear Port-Hamiltonian Differential-Algebraic Equations of Index~1

Aashutosh Sharma , Andreas Bartel , Manuel Schaller

Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:02 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords port-Hamiltonian systemsdifferential-algebraic equationstime adaptivitydual weighted residualenergy balancea posteriori error estimationelectrical circuit simulation

0 comments

The pith

Goal-oriented time adaptivity controls energy balance violations in linear port-Hamiltonian DAEs of index 1.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a method to adapt time steps in simulations of port-Hamiltonian differential-algebraic equations so that the discrete scheme respects the continuous energy balance more closely. It does this by applying the dual weighted residual technique to derive goal-oriented a posteriori error indicators that drive local grid refinement. The port-Hamiltonian structure is used to replace expensive linear solves in the error estimator with a cheaper dissipativity-exploiting block-Jacobi approximation. The resulting adaptive algorithm is tested on linear electrical circuit models, where it reduces energy drift without a large increase in total steps. A reader interested in structure-preserving numerics would care because energy balance is a defining physical invariant that standard integrators often violate.

Core claim

We propose an approach that controls the energy balance violation for port-Hamiltonian differential algebraic equations via time adaptivity using a posteriori grid refinement techniques based on the dual weighted residual method. In particular, we show how one may leverage the port-Hamiltonian structure to efficiently compute the error estimators using a dissipativity-exploiting block-Jacobi approximation. We illustrate the efficacy of the method by means of simulations of electrical circuit models.

What carries the argument

Dissipativity-exploiting block-Jacobi approximation inside a dual-weighted-residual error estimator that drives goal-oriented time-step refinement for port-Hamiltonian DAEs.

If this is right

Adaptive time grids can be generated automatically that keep the discrete energy balance within a prescribed tolerance for linear index-1 pH-DAEs.
The computational cost of the error estimator stays linear in the number of degrees of freedom because the block-Jacobi approximation avoids full coupled solves.
The same indicator can be reused to control other goal functionals that depend on the energy flow, not only the balance itself.
The method extends the range of problems for which structure-preserving time integration is feasible without redesigning the integrator.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same structure-exploiting estimator idea could be combined with existing discrete-gradient schemes to obtain both local accuracy and global conservation.
If the block-Jacobi approximation degrades for strongly dissipative or stiff circuits, the method would naturally fall back to more expensive but exact error indicators.
Because the refinement criterion is goal-oriented, the same framework could be applied to output functionals other than energy, such as voltage at a specific node.

Load-bearing premise

The block-Jacobi approximation that exploits dissipativity remains sufficiently accurate for the chosen goal quantity so that the resulting error indicators still drive useful mesh refinement.

What would settle it

A sequence of refined grids on a known linear circuit where the measured energy-balance residual fails to decrease monotonically once the block-Jacobi estimator is replaced by the exact dual-weighted residual.

Figures

Figures reproduced from arXiv: 2605.14082 by Aashutosh Sharma, Andreas Bartel, Manuel Schaller.

**Figure 1.** Figure 1: Spectral radius ρ(Γi) of the amplification matrix (65), tested on (70). Left: distribution over [0, T] on a uniform (N = 200) and DWR-adaptive (N = 302) grid; the adaptive grid places short intervals near t = 0 where the input is active, resulting in ρ(Γi) close to one in that region. Right: ρ(Γi) versus step size ki ; all values confirm the theoretical decay of (1 + ki µmin) −1 . 6 Numerical experiments T… view at source ↗

**Figure 2.** Figure 2: Convergence of uniform and DWR adaptive dG(0). (a) Error in QoI [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗

**Figure 3.** Figure 3: Evolution of time grid. Importantly, the mesh remains relatively coarse for t > 0.6, confirming that once the system reaches a more dissipative state (cf. u in (67)), the local energy balance is already well resolved without the need for additional degrees of freedom [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗

**Figure 4.** Figure 4: Structure of the DWR error indicator. Left: local energy residual magnitude [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗

**Figure 5.** Figure 5: RCL ladder transmission line [23, 24] with building block given in dashed box. The right end [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗

**Figure 6.** Figure 6: RCL ladder transmission line. Left: node voltages at nodes 1, 6, 11, 16, 21, 26 on the final adaptive grid; the excitation propagates from left to right with increasing delay and decay of amplitude. Right: energy-balance violation P i G2 i under uniform (•) and DWR-adaptive (▲) dG(0) refinement; the dotted line gives the fitted rate O(N −3.0 ). Physical behavior. Figure 6a shows the node voltage vk(t) at s… view at source ↗

**Figure 7.** Figure 7: Effect of the Hamiltonian-weighted QoI on the transmission-line benchmark (70) ( [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗

**Figure 8.** Figure 8: Late adaptive grid (N = 143, |Mex| = 12, 14 DWR refinements). (a) Relative indicator error (solid) and adjoint error (dashed), with marked-set agreement (right axis, green). Stabilization at k ∗ = 110. (b) DWR indicators ranked by |ηi | for k ∈ {5, 20, 40}. Jacobi approximations match the exact curve within the marked region well before k ∗ . 22 [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗

read the original abstract

Port-Hamiltonian systems provide a highly-structured framework for modeling of physical systems. By definition, they encode a balance equation relating energy changes to supplied and dissipated energy. Capturing this energy balance in discrete approximations is a fundamental challenge and often has been achieved by designing particular schemes such as discrete gradient methods. In this work, we propose an approach that controls the energy balance violation for port-Hamiltonian differential algebraic equations via time adaptivity using a posteriori grid refinement techniques based on the dual weighted residual method. In particular, we show how one may leverage the port-Hamiltonian structure to efficiently compute the error estimators using a dissipativity-exploiting block-Jacobi approximation. We illustrate the efficacy of the method by means of simulations of electrical circuit models.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper offers a DWR-based time adaptivity scheme for controlling energy balance violations in linear port-Hamiltonian DAEs via a dissipativity-exploiting block-Jacobi adjoint approximation, but lacks bounds or quantitative checks on that approximation.

read the letter

The main point is a goal-oriented adaptive time-stepping method for linear index-1 port-Hamiltonian DAEs that uses dual-weighted residuals to drive refinement so the discrete energy balance stays closer to the continuous one. The key efficiency trick is approximating the adjoint system with a block-Jacobi scheme that exploits the dissipativity structure instead of solving the full coupled adjoint each time, and they demonstrate it on electrical circuit models. What the paper does reasonably well is frame a practical numerical problem—energy drift in structure-preserving discretizations—and show how the port-Hamiltonian form can be turned into a cheaper error estimator. Applying DWR specifically to the energy-balance functional and simplifying the adjoint via block-Jacobi is a concrete engineering move that could save work in circuit simulation or similar physical models. The circuit examples at least give a sense that the adaptive grids behave as intended. The soft spot is the missing support for the central approximation. The block-Jacobi scheme is load-bearing for practicality, yet the text supplies no a-priori bound on its deviation from the true adjoint solution and no side-by-side numbers comparing the approximate versus exact error indicators on the test cases. Without that, it is hard to know whether the adaptivity is reliably targeting the claimed quantity. The abstract also omits any derivation of the estimators or tabulated quantitative results, so the evidence for efficacy is mostly that the code runs and produces plausible grids. This is aimed at specialists working on numerical methods for DAEs and port-Hamiltonian systems who need implementable adaptivity ideas rather than new theory. A reader in that niche could extract usable implementation details, but anyone needing rigorous error control or proven reliability of the approximation will want more. I would send it to peer review. The core construction is worth a referee's time to check the analysis and validation that the abstract leaves out.

Referee Report

1 major / 2 minor

Summary. The manuscript proposes a goal-oriented time-adaptive method for linear port-Hamiltonian DAEs of index 1 that controls violations of the energy balance equation. It employs a posteriori grid refinement via the dual-weighted residual (DWR) method and exploits the port-Hamiltonian dissipativity structure through a block-Jacobi approximation to the adjoint system for efficient estimator computation. The approach is demonstrated on electrical-circuit examples.

Significance. If the block-Jacobi approximation to the adjoint remains sufficiently accurate, the method could offer an efficient structure-preserving route to goal-oriented adaptivity for energy-balance control in pH-DAE simulations, a practically relevant task in circuit and multibody modeling. The work correctly identifies the energy-balance violation as a natural goal functional and shows how the pH structure can be used to simplify the estimator; however, the absence of any a-priori accuracy statement or quantitative validation of the approximation prevents a firm assessment of its reliability or novelty relative to existing DWR literature.

major comments (1)

[§3–4 (adjoint approximation and numerical examples)] The central construction (abstract and §3–4) relies on a dissipativity-exploiting block-Jacobi approximation to the adjoint system for the DWR estimator. No a-priori bound is given on the difference between this approximation and the true coupled adjoint solution, nor is any quantitative comparison (e.g., relative difference in the goal functional or effect on the adapted mesh) reported for the electrical-circuit examples. This approximation is load-bearing for the claim that the resulting indicators reliably drive refinement of the energy-balance violation.

minor comments (2)

[Abstract] The abstract states that the method “controls the energy balance violation” but supplies no concrete definition of the goal functional or the precise measure of violation used in the DWR estimator.
[§2] Notation for the port-Hamiltonian DAE (E, J, R, B, …) is introduced without an explicit reference to the standard linear pH-DAE form; a short displayed equation would improve readability.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments on the significance of the work. We address the major concern regarding the block-Jacobi approximation below and will revise the manuscript to include additional validation.

read point-by-point responses

Referee: [§3–4 (adjoint approximation and numerical examples)] The central construction (abstract and §3–4) relies on a dissipativity-exploiting block-Jacobi approximation to the adjoint system for the DWR estimator. No a-priori bound is given on the difference between this approximation and the true coupled adjoint solution, nor is any quantitative comparison (e.g., relative difference in the goal functional or effect on the adapted mesh) reported for the electrical-circuit examples. This approximation is load-bearing for the claim that the resulting indicators reliably drive refinement of the energy-balance violation.

Authors: We agree that the manuscript would benefit from explicit validation of the block-Jacobi approximation. The approximation is motivated by the dissipativity structure of the port-Hamiltonian system, which allows the adjoint to be decoupled into independent forward and backward solves while preserving the essential energy-balance properties. Although no a-priori bound is derived, the numerical examples already demonstrate that the adapted meshes successfully reduce the energy-balance violation to the desired tolerance. In the revised manuscript we will add quantitative comparisons: for each circuit example we will report the relative difference between the goal functional evaluated with the approximate versus the fully coupled adjoint, together with a side-by-side comparison of the resulting adapted time grids. These additions will confirm that the approximation error remains small enough to preserve the reliability of the refinement indicators. revision: yes

standing simulated objections not resolved

A rigorous a-priori error bound on the difference between the block-Jacobi approximation and the true coupled adjoint solution.

Circularity Check

0 steps flagged

No significant circularity; derivation uses standard DWR with structure-exploiting approximation independent of inputs

full rationale

The paper proposes a goal-oriented adaptive time-stepping scheme for linear port-Hamiltonian DAEs that controls energy-balance violation via the dual-weighted residual method. The central technical step is the introduction of a dissipativity-exploiting block-Jacobi approximation to the adjoint system, which is presented as a practical efficiency device rather than a fitted or self-defined quantity. No equation is shown to reduce by construction to a prior fit, no prediction is statistically forced by the same data used for calibration, and no load-bearing uniqueness theorem or ansatz is imported solely via self-citation. The derivation chain therefore remains self-contained against external benchmarks such as the standard DWR theory and the port-Hamiltonian energy-balance identity.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only abstract available; no explicit free parameters, axioms, or invented entities can be identified. The approach implicitly relies on standard assumptions of numerical analysis for linear index-1 DAEs and port-Hamiltonian structure.

pith-pipeline@v0.9.0 · 5432 in / 1044 out tokens · 53624 ms · 2026-05-15T02:02:29.587510+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel (J uniqueness) echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

dissipativity-exploiting block-Jacobi approximation... Γᵢ is a strict contraction... α>0 eigenvalue of eS
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

energy balance violation as QoI... dual weighted residual method

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

[1]

Springer Basel AG, 2003

Wolfgang Bangerth and Rolf Rannacher.Adaptive Finite Elements Methods for Differential Equa- tions. Springer Basel AG, 2003

work page 2003
[2]

Goal-oriented time adaptivity for port-Hamiltonian systems

Andreas Bartel and Manuel Schaller. Goal-oriented time adaptivity for port-Hamiltonian systems. Journal of Computational and Applied Mathematics, 461:116450, 2025. ISSN 0377-0427

work page 2025
[3]

Port-Hamiltonian systems modelling in electrical engineering, 2023

Andreas Bartel, Markus Clemens, Michael G¨ unther, Birgit Jacob, and Timo Reis. Port-Hamiltonian systems modelling in electrical engineering, 2023. URLhttps://arxiv.org/abs/2301.02024

work page arXiv 2023
[4]

Port-Hamiltonian descriptor systems.Mathematics of Control, Signals, and Systems, 30(4), 2018

Christopher Beattie, Volker Mehrmann, Hongguo Xu, and Hans Zwart. Port-Hamiltonian descriptor systems.Mathematics of Control, Signals, and Systems, 30(4), 2018

work page 2018
[5]

An optimal control approach to a posteriori error estimation in finite element methods.Acta Numerica, 10:1–102, 2001

Rolf Becker and Rolf Rannacher. An optimal control approach to a posteriori error estimation in finite element methods.Acta Numerica, 10:1–102, 2001

work page 2001
[6]

Port-Hamiltonian of rigid multi- body systems, 2025

Thomas Berger, Ren´ e Hochdahl, Timo Reis, and Robert Seifried. Port-Hamiltonian of rigid multi- body systems, 2025. URLhttps://arxiv.org/abs/2504.17063

work page arXiv 2025
[7]

A port-Hamiltonian system perspective on electromagneto-quasistatic field formulations of Darwin-type, 2024

Markus Clemens, Marvin-Lucas Henkel, Fotios Kasolis, and Michael G¨ unther. A port-Hamiltonian system perspective on electromagneto-quasistatic field formulations of Darwin-type, 2024. URL https://arxiv.org/abs/2404.18767

work page arXiv 2024
[8]

Endtmayer, U

B. Endtmayer, U. Langer, and T. Wick. Two-side a posteriori error estimates for the dual-weighted residual method.SIAM Journal on Scientific Computing, 42(1):A371–A394, 2020. 25

work page 2020
[9]

Optimal control of port-Hamiltonian descriptor systems with minimal energy supply.SIAM Journal on Control and Optimization, 60(4):2132–2158, 2022

Timm Faulwasser, Bernhard Maschke, Friedrich Philipp, Manuel Schaller, and Karl Worthmann. Optimal control of port-Hamiltonian descriptor systems with minimal energy supply.SIAM Journal on Control and Optimization, 60(4):2132–2158, 2022

work page 2022
[10]

Roland W. Freund. The SPRIM algorithm for structure-preserving order reduction of general RCL circuits. In P. Benner, M. Hinze, and E.J.W. ter Maten, editors,Model Reduction for Circuit Simulation, volume 74 ofLect. Notes Electr. Eng., chapter 2, pages 25–52. Springer, Dordrecht, 2011

work page 2011
[11]

Time integration and discrete Hamiltonian systems.Journal of Nonlinear Science, 6:449–467, 1996

Oscar Gonzalez. Time integration and discrete Hamiltonian systems.Journal of Nonlinear Science, 6:449–467, 1996

work page 1996
[12]

Efficient model predictive control for parabolic PDEs with goal oriented error estimation.SIAM Journal on Scientific Computing, 44(1):A471–A500, 2022

Lars Gr¨ une, Manuel Schaller, and Anton Schiela. Efficient model predictive control for parabolic PDEs with goal oriented error estimation.SIAM Journal on Scientific Computing, 44(1):A471–A500, 2022

work page 2022
[13]

Springer Science & Business Media, 2012

Birgit Jacob and Hans J Zwart.Linear port-Hamiltonian systems on infinite-dimensional spaces. Springer Science & Business Media, 2012

work page 2012
[14]

PhD thesis, Humboldt-Universit¨ at zu Berlin, Mathematisch-Naturwissenschaftliche Fakult¨ at I, 2015

Lennart Jansen.A Dissection concept for DAEs. PhD thesis, Humboldt-Universit¨ at zu Berlin, Mathematisch-Naturwissenschaftliche Fakult¨ at I, 2015

work page 2015
[15]

Discrete gradient methods for port- Hamiltonian differential-algebraic equations.Preprint arXiv:2505.18810, 2025

Philipp L Kinon, Riccardo Morandin, and Philipp Schulze. Discrete gradient methods for port- Hamiltonian differential-algebraic equations.Preprint arXiv:2505.18810, 2025

work page arXiv 2025
[16]

Discrete-time port-Hamiltonian systems: A definition based on symplectic integration.Systems & Control Letters, 133:104530, 2019

Paul Kotyczka and Laurent Lefevre. Discrete-time port-Hamiltonian systems: A definition based on symplectic integration.Systems & Control Letters, 133:104530, 2019

work page 2019
[17]

Linear algebra properties of dissipative Hamiltonian descriptor systems.SIAM Journal on Matrix Analysis and Applications, 39(3):1489– 1519, 2018

Christian Mehl, Volker Mehrmann, and Michal Wojtylak. Linear algebra properties of dissipative Hamiltonian descriptor systems.SIAM Journal on Matrix Analysis and Applications, 39(3):1489– 1519, 2018

work page 2018
[18]

Matrix pencils with coefficients that have positive semidefinite Hermitian parts.SIAM Journal on Matrix Analysis and Applications, 43(3): 1186–1212, 2022

Christian Mehl, Volker Mehrmann, and Michal Wojtylak. Matrix pencils with coefficients that have positive semidefinite Hermitian parts.SIAM Journal on Matrix Analysis and Applications, 43(3): 1186–1212, 2022

work page 2022
[19]

Structure-preserving discretization for port-Hamiltonian descriptor systems

Volker Mehrmann and Riccardo Morandin. Structure-preserving discretization for port-Hamiltonian descriptor systems. In2019 IEEE 58th Conference on Decision and Control (CDC), pages 6863–6868, 2019

work page 2019
[20]

Control of port-Hamiltonian differential-algebraic systems and applications.Acta Numerica, 32:395–515, 2023

Volker Mehrmann and Benjamin Unger. Control of port-Hamiltonian differential-algebraic systems and applications.Acta Numerica, 32:395–515, 2023

work page 2023
[21]

A priori error estimates for space-time finite element discretiza- tion of parabolic optimal control problems part I: Problems without control constraints.SIAM J

Dominik Meidner and Boris Vexler. A priori error estimates for space-time finite element discretiza- tion of parabolic optimal control problems part I: Problems without control constraints.SIAM J. Control Optim., 47(3):1150–1177, 2008

work page 2008
[22]

PhD thesis, Technische Universit¨ at Berlin, 2024

Riccardo Morandin.Modeling and numerical treatment of port-Hamiltonian descriptor systems. PhD thesis, Technische Universit¨ at Berlin, 2024

work page 2024
[23]

PortHamiltonianBenchmarkSystems, 2021

Tim Moser et al. PortHamiltonianBenchmarkSystems, 2021. URLhttps://algopaul.github.io/ PortHamiltonianBenchmarkSystems

work page 2021
[24]

PhD thesis, Technical University of Munich, 2023

Tim Fabian Moser.Structure-Preserving Model Reduction of Port-Hamiltonian Descriptor Systems. PhD thesis, Technical University of Munich, 2023

work page 2023
[25]

Structure-preserving time discretization of port-Hamiltonian systems via discrete gradient pairs.Preprint arXiv:2311.00403, 2023

Philipp Schulze. Structure-preserving time discretization of port-Hamiltonian systems via discrete gradient pairs.Preprint arXiv:2311.00403, 2023

work page arXiv 2023
[26]

Springer Series in Com- putational Mathematics, Berlin, 2006

Vidar Thom´ ee.Galerkin Finite Element Methods for Parabolic Problems. Springer Series in Com- putational Mathematics, Berlin, 2006

work page 2006
[27]

Port-Hamiltonian systems theory: An introductory overview.Foundations and Trends in Systems and Control, 1(2–3):173–378, 2014

Arjan van der Schaft and Dimitri Jeltsema. Port-Hamiltonian systems theory: An introductory overview.Foundations and Trends in Systems and Control, 1(2–3):173–378, 2014. 26

work page 2014