[论文解读] 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

一句话总结

本文提出了一种用于古典吉他流固耦合的端口-哈密顿（pH）模型，通过有限元空间离散化和状态变换，将固体力学与流体动力学相结合。该研究提出了一种结构保持的模型降阶（MOR）框架，采用基于基的投影方法，表明SVD-like和辛几何基方法在降低基大小（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.

研究动机与目标

• 开发一种统一的端口-哈密顿（pH）公式，用于古典吉他的流固耦合（FSI），整合固体力学与流体动力学，并包含结构阻尼。
• 实现瞬态时域仿真，与典型的频域声学分析形成对比。
• 比较多种结构保持的模型降阶（MOR）技术，重点关注与基无关的基生成与投影方法。
• 评估轨迹数量与基大小对MOR中近似误差与计算效率的影响。
• 建立结合商业软件与自研软件的流程，实现多物理场耦合与MOR，确保能量稳定与无源性。

提出的方法

• 通过有限元空间离散化，将线性弹性（固体）与波动方程（流体）结合，推导出二阶端口-哈密顿系统。
• 应用状态变换，将二阶系统转换为一阶pH系统，保持能量流动与耗散特性。
• 采用与基无关的结构保持投影方法，确保降阶模型中具有无源性、稳定性与能量守恒。
• 使用多种技术生成降阶基：模态截断、Krylov子空间方法、 Proper Orthogonal Decomposition（POD）、 Proper Symplectic Decomposition（PSD）以及SVD-like分解。
• 利用瞬态仿真中的时间响应快照构建数据驱动基，轨迹数量在5至100之间变化，用于敏感性分析。
• 采用动量形式的pH投影，确保在模型降阶过程中保持端口-哈密顿结构。

实验结果

研究问题

• RQ1 如何利用统一的端口-哈密顿框架，准确建模古典吉他中的流固耦合系统，同时保持能量与耗散特性？
• RQ2 哪种基生成与投影方法的组合能够生成最准确且稳定的降阶模型，同时保持原始pH系统的结构保持特性？
• RQ3 时间响应轨迹数量如何影响MOR中数据驱动基生成的近似误差与计算成本？
• RQ4 在FSI系统的结构保持MOR中，降阶基大小、近似误差与计算加速之间的权衡关系如何？
• RQ5 在精度与效率方面，辛几何与SVD-like基生成方法相较于标准POD与模态方法，优势有多大？

主要发现

• SVD-like基生成方法在基大小n=120时达到最低的相对误差（0.0017），优于所有其他方法，包括POD-State与PSD。
• 当基大小为n=12时，实现了高达540倍的加速，但伴随较高的近似误差，凸显了速度与精度之间的权衡。
• 在基大小n=120时，除少数方法外，其余方法的相对误差均低于0.005，其中SVD-like与POD-State方法达到最低值，表明在中等维度下具有高精度。
• POD-Displacement基未能有效激发系统，导致响应接近零，表明其能量捕获能力差，不适用于此FSI问题。
• 当轨迹数量超过25条后，误差的改善趋于平缓，10条轨迹已能实现精度与离线计算成本之间的良好平衡。
• 结构保持的MOR方法，特别是pH保持与能量稳定的投影方法，在稳定性与精度方面显著优于非结构保持的替代方法。