Back to search
2505.18810

Discrete gradient methods for port-Hamiltonian differential-algebraic equations

Philipp L. Kinon, Riccardo Morandin, Philipp Schulze

correctmedium confidence
Category
Not specified
Journal tier
Strong Field
Processed
Sep 28, 2025, 12:56 AM

Audit review

The paper rigorously proves the local canonical form (Theorem 5.16) and its pHDAE consequence (Corollary 5.17) using a two–step construction: (i) a constant-rank, analytic reduction U⊤EΠ = [Ir E12; 0 0] (Corollary F.2) and (ii) a diffeomorphism built via a first-order PDE Dx2ψ1 = (Dx1ψ1)E12 that zeroes the off-pivot block, yielding Ẽ = [Ip 0 0; 0 0 0; 0 E32 I] and a specified Hamiltonian H̃1; see Theorem 5.16 and the Appendix F proof sketch that implements these steps and then derives (41a)–(41b) . The candidate’s outline omits the essential PDE-based construction and contains a false claim that if ∇H(x⋆)=0 then H∘φ is locally constant; it also does not justify how left-multiplication alone enforces the required zero structure without the right-side Jacobian change. The paper’s argument is therefore correct and complete, while the model’s is flawed/incomplete (e.g., it relies on Step 1 and Step 3 that do not hold as stated) .

Referee report (LaTeX)

\textbf{Recommendation:} minor revisions

\textbf{Journal Tier:} strong field

\textbf{Justification:}

The paper presents structure-preserving discretizations for nonlinear pHDAEs and a precise structural decomposition under analytic, constant-rank assumptions on E. The derivation of a local canonical form (and its implications for semi-explicit formulations) addresses an important modeling and numerical-analysis gap and is done carefully. Some intricate steps are relegated to the appendix; moving a brief summary of these steps (especially the PDE-based coordinate construction) into the main text would improve accessibility. Overall, the contribution is solid, technically sound, and relevant.