[論文レビュー] Structure-Preserving Time Discretization of Port-Hamiltonian Systems via Discrete Gradient Pairs

Philipp Schulze|arXiv (Cornell University)|Nov 1, 2023

Numerical methods for differential equations被引用数 4

ひとこと要約

本稿では、状態に依存する質量行列を有する非線形ポート・ハミルトニアン系に対する構造保存型時間離散化手法として、離散勾配ペアを新たに導入する。対称的かつ正定値の質量行列に対しては、中点離散勾配ペアを構築することで、時間離散化レベルにおいて正確なパワー・バランスを保証し、低次元モデルにおける主要な物理的性質を保存する。

ABSTRACT

We discuss structure-preserving time discretization for nonlinear port-Hamiltonian systems with state-dependent mass matrix. Such systems occur, for instance, in the context of structure-preserving nonlinear model order reduction for port-Hamiltonian systems and, in this context, structure-preserving time discretization is crucial for preserving some of the properties of the time-continuous reduced-order model. For this purpose, we introduce a new class of time discretization schemes which is based on so-called discrete gradient pairs and leads to an exact power balance on the time-discrete level. Moreover, for the special case of a pointwise symmetric and positive definite mass matrix, we present an explicit construction of a discrete gradient pair. Finally, we illustrate the theoretical findings by means of a numerical example, where the time-continuous system is a nonlinear reduced-order model for an advection-diffusion problem.

研究の動機と目的

• 状態に依存する質量行列を有する非線形ポート・ハミルトニアン系の時間離散化において、物理的構造（特にパワー・バランス）を保存する課題に取り組む。
• 標準的な離散勾配法がハミルトニアン勾配が系方程式に明示的に現れる必要があるという制限を克服するため、離散勾配ペアを導入する。
• 局所的に対称的かつ正定値の質量行列を有する場合の、離散勾配ペアの明示的構成を提供する。
• 非線形拡散・対流問題の非線形低次元モデルに対して、本手法の有効性を示し、正確な時間離散パワー・バランスと2次収束性を確認する。
• 連続系と同一の時間離散的性質を保証することで、構造保存型モデル次元削減を可能にする。

提案手法

• ハミルトニアン勾配が力学方程式に明示的に現れない系のための離散勾配の一般化として、離散勾配ペアの概念を導入する。
• 離散平均値性質を満たす関数のペアとして離散勾配ペアを定義し、時間離散化レベルで正確なパワー・バランスを実現可能にする。
• 対称的かつ正定値の質量行列に対しては、ハミルトニアンの勾配に中点則を適用することで、中点離散勾配ペアを構築する。
• 離散勾配ペアに基づく時間離散化スキームを導出し、ポート・ハミルトニアン系の代数的構造を保存する。
• 時間離散系が正確なパワー・バランス方程式を満たすことを保証する：$ \frac{\mathcal{H}(x^{n+1}) - \mathcal{H}(x^n)}{\Delta t} = -y^n \cdot u^n $、非線形反復誤差を除き。
• 数値実験において非線形反復誤差を最小限に抑えるために、厳密な収束基準を設定したニュートン型ソルバーを用いてスキームを実装する。

実験結果

リサーチクエスチョン

• RQ1 離散勾配ペアを用いて、状態に依存する質量行列を有する非線形ポート・ハミルトニアン系の構造保存型時間積分スキームを構築可能か？
• RQ2 ハミルトニアンが二次形式でない場合でも、本手法が時間離散化レベルで正確なパワー・バランスを保証するか？
• RQ3 対称的かつ正定値の質量行列を有する場合に、明示的な離散勾配ペアを構成可能か？
• RQ4 離散勾配ペア法の収束特性は、インライン・ミッドポイント則などの古典的スキームと比較してどのように異なるか？
• RQ5 本手法は、偏微分方程式の非線形低次元モデルにおいて、物理的性質をどの程度正確に保存するか？

主な発見

• 提案された離散勾配ペア法は、定理4.3で示されるように、非二次ハミルトニアンに対しても時間離散化レベルで正確なパワー・バランスを達成する。
• 対称的かつ正定値の質量行列に対しては、定理4.2で示されるように、中点離散勾配ペアが明示的かつ計算可能な構成である。
• 対流・拡散問題の非線形低次元モデルに対する数値実験では、離散勾配ペア法がインライン・ミッドポイント則と比較して、パワー・バランス誤差が数個桁小さいことが確認された。
• 図2の収束プロットから確認されるように、中点離散勾配ペア法の実験的収束次数は約2であり、インライン・ミッドポイント則と同等の性能を示す。
• 本手法は、ハミルトニアン勾配が系方程式に明示的に現れない場合でも構造保存性を維持でき、従来の離散勾配手法の主な制限を克服する。
• 数値実験のためのコードは、DOI 10.5281/zenodo.10059715 を通じて公開されており、再現性とさらなる研究を支援する。