[논문 리뷰] Structure-Preserving Time Discretization of Port-Hamiltonian Systems via Discrete Gradient Pairs

Philipp Schulze|arXiv (Cornell University)|2023. 11. 01.

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 이 방법은 PDE의 비선형 감소된 차수 모델에서 물리적 성질을 어느 정도 유지하는가?

주요 결과

• 제안된 이산 기울기 쌍 방법은 비선형 해밀토니안이더라도 정리 4.3에 의해 시간 이산화 수준에서 정확한 전력 균형을 달성한다.
• 대칭적이고 양의 정부호인 질량 행렬의 경우, 중점 이산 기울기 쌍은 정리 4.2에 의해 명시적이고 계산 가능한 구성 방법을 제공한다.
• 열화-확산 문제의 비선형 감소된 차수 모델에 대한 수치 실험 결과, 이산 기울기 쌍 방법은 은밀한 중점 규칙보다 전력 균형 오차가 수 개의 주자수만큼 작다.
• 중점 이산 기울기 쌍 방법의 실험적 수렴 차수는 약 2이며, 그림 2의 수렴도를 통해 은밀한 중점 규칙과 동일한 것을 확인하였다.
• 이 방법은 시스템 방정식에 해밀토니안 기울기가 명시적으로 나타나지 않아도 구조 보존를 유지할 수 있으며, 이는 이전의 이산 기울기 접근법의 핵심 한계를 극복한다.
• 수치 실험을 위한 코드는 DOI 10.5281/zenodo.10059715 를 통해 공개되어 있어 재현성과 향후 연구를 지원한다.