[論文レビュー] Port-Hamiltonian Fluid-Structure Interaction Modeling and Structure-Preserving Model Order Reduction of a Classical Guitar

Johannes Rettberg, Dominik Wittwar+5 more|arXiv (Cornell University)|Mar 18, 2022

Model Reduction and Neural Networks被引用数 1

ひとこと要約

本稿では、固体および流体の動力学を有限要素離散化と状態変換を用いて統合することで、クラシックギターにおける流体構造挙動相互作用（FSI）を記述するポート・ハミルトニアン（pH）モデルを提示する。構造保存型モデル順序低減（MOR）フレームワークを提案し、基底に依存しない射影法を用いる。特筆すべきは、特異値分解（SVD）に類似した手法およびシンプレクティックに基づく基底が、基底サイズを小さく（n=120）した状態で、最大540倍の高速化を達成するとともに、高い精度を発揮することである。非構造保存型手法に比べて優れた性能を示す。

ABSTRACT

A fluid-structure interaction model in a port-Hamiltonian representation is derived for a classical guitar. We combine the laws of continuum mechanics for solids and fluids within a unified port-Hamiltonian (pH) modeling approach by adapting the discretized equations on second-order level in order to obtain a damped multi-physics model. The high-dimensionality of the resulting system is reduced by model order reduction. The article focuses on pH-systems in different state transformations, a variety of basis generation techniques as well as structure-preserving model order reduction approaches that are independent from the projection basis. As main contribution a thorough comparison of these method combinations is conducted. In contrast to typical frequency-based simulations in acoustics, transient time simulations of the system are presented. The approach is embedded into a straightforward workflow of sophisticated commercial software modeling and flexible in-house software for multi-physics coupling and model order reduction.

研究の動機と目的

• クラシックギターにおける流体構造相互作用（FSI）を統一されたポート・ハミルトニアン（pH）形式で定式化し、固体および流体の動力学に加え、構造的減衰を統合すること。
• 通常の周波数領域の音響解析とは対照的に、時間領域における過渡応答シミュレーションを可能にすること。
• 構造保存型モデル順序低減（MOR）手法の比較を実施し、基底生成法および基底に依存しない射影法に焦点を当てる。
• MORにおけるトレーリング数および基底サイズが近似誤差および計算効率に与える影響を評価すること。
• 商用ソフトウェアと自作ソフトウェアを統合したワークフローを確立し、マルチフィジックス結合およびMORを実現。エネルギー安定性およびパassivityを保証すること。

提案手法

• 線形弾性（固体）と波動方程式（流体）を有限要素空間離散化により統合し、2階のポート・ハミルトニアン系を導出する。
• 状態変換を適用して2階系を1階系のpH系に変換し、エネルギーの流れと散逸を保存する。
• 基底に依存しない構造保存型射影法を用い、低次元モデルにおけるパassivity、安定性、エネルギー保存を保証する。
• 複数の手法を用いて低次元基底を生成：モード截断、クリロフ部分空間法、固有直交分解（POD）、固有シンプレクティック分解（PSD）、SVDに類似した分解法。
• 過渡応答シミュレーションからのスナップショットを用いてデータ駆動型基底を構築し、トレーリング数を5〜100に変化させて感度分析を実施。
• 運動量形式を用いたpH射影を実装し、モデル順序低減の過程でポート・ハミルトニアン構造を維持する。

実験結果

リサーチクエスチョン

• RQ1 クラシックギターにおける流体構造相互作用系を、エネルギーと散逸を保存する統一されたポート・ハミルトニアンフレームワークで正確にモデル化するにはどうすればよいか？
• RQ2 基底生成法と射影法の組み合わせの中で、元のpH系の構造保存性を維持しつつ、最も正確で安定した低次元モデルを生成するのはどれか？
• RQ3 時間応答トレーリング数が、MORにおけるデータ駆動型基底生成の近似誤差および計算コストに与える影響はいかほどか？
• RQ4 構造保存型MORにおいて、低次元基底サイズ、近似誤差、計算高速化の間にはどのようなトレードオフがあるか？
• RQ5 SVDに類似した基底生成法およびシンプレクティック基底生成法は、標準的なPODおよびモードアプローチに比べ、精度および効率の面でどれほど優れているか？

主な発見

• SVDに類似した基底生成法は、基底サイズn=120で相対誤差0.0017を達成し、POD-StateやPSDを含む他のすべての手法を上回る性能を示した。
• 基底サイズn=12で最大540倍の高速化を達成したが、これは高い近似誤差を伴うことを示しており、速度と精度のトレードオフが顕著に現れている。
• 基底サイズn=120では、ほとんどの手法で相対誤差が0.005未満に低下し、SVDに類似した手法およびPOD-Stateが最小値を記録した。これは中程度の次元で高い精度を達成していることを示している。
• POD-変位基底は系を十分に励起できず、応答はほぼゼロに近づき、エネルギーの捕捉が不十分であり、本FSI問題には不適切であることが判明した。
• トレーリング数が25を超えると誤差への影響が段階的に小さくなる傾向を示し、10トレーリングで精度とオフライン計算コストのバランスが良好に取れていることがわかった。
• 構造保存型MOR手法、特にpH保存型およびエネルギー安定型射影法は、非構造保存型手法に比べ、安定性および精度の両面で顕著に優れていた。