[论文解读] Isogeometric Boundary Element Method with Hierarchical Matrices

Jürgen Zechner, Benjamin Marussig+2 more|arXiv (Cornell University)|Jun 11, 2014

Advanced Numerical Analysis Techniques参考文献 14被引用 5

一句话总结

本文提出一种基于等几何边界元法（IGA-BEM）与层级矩阵（H-矩阵）的线性弹性问题求解方法，利用NURBS几何表示与核插值技术对满秩系统矩阵进行压缩。该方法实现了近线性复杂度与高精度，2D数值实验结果表明其具有最优收敛性与高效的矩阵压缩性能。

ABSTRACT

In this work we address the complexity problem of the isogeometric Boundary Element Method by proposing a collocation scheme for practical problems in linear elasticity and the application of hierarchical matrices. For mixed boundary value problems, a block system of matrices similar to Galerkin formulations is constructed allowing an effective application of that matrix format. We introduce a strategy for the geometric bisection of surfaces based on NURBS patches. The approximation of system matrices is carried out by means of kernel interpolation. Numerical results are shown that prove the success of the formulation.

研究动机与目标

• 解决等几何边界元法因满秩系统矩阵导致的高计算成本问题。
• 将层级矩阵（H-矩阵）集成至等几何配点BEM框架中，用于求解线性弹性问题。
• 在保持NURBS几何精度的前提下，实现大规模问题的高效矩阵压缩与求解。
• 基于NURBS节点插入与凸包性质，提出一种几何二分策略，用于矩阵块聚类。
• 通过拉格朗日多项式插值与QR压缩，确保矩阵近似具有高精度。

提出的方法

• 使用NURBS片段精确表示边界几何，消除几何离散化误差。
• 采用基于Lamé-Navier算子的配点边界元法，求解弹性问题的边界积分方程。
• 构建与Galerkin格式类似的块矩阵系统，适用于混合边界值问题。
• 通过NURBS节点插入与凸包性质生成的包围盒，实现对表面的几何二分。
• 通过在允许聚类上对基本解（Kelvin解）进行核插值，近似远场矩阵块。
• 采用低秩分解（A·S·Bᵀ）与QR分解对低秩矩阵块进行压缩，降低存储与计算成本。

实验结果

研究问题

• RQ1 H-矩阵能否有效降低等几何BEM在求解线性弹性问题时的计算复杂度？
• RQ2 如何利用基于NURBS的几何表示，实现H-矩阵压缩所需的高效几何聚类？
• RQ3 所提出的H-矩阵近似在弹性问题中单层与双层势函数的精度与收敛行为如何？
• RQ4 通过核插值与QR分解实现的矩阵压缩，如何在不损失精度的前提下提升存储与求解效率？
• RQ5 该方法在实际应用中能否在实现近线性复杂度的同时保持最优收敛率？

主要发现

• 2D数值实验表明，所提出的IGA-BEM结合H-矩阵在单层与双层势函数上均实现了最优收敛率。
• 矩阵压缩比（cₕ）随问题规模接近线性增长，表明实际计算复杂度接近线性。
• 即使在高阶压缩水平下，方法仍保持高精度，误差水平由插值阶次与秩选择控制。
• 基于拉格朗日多项式与支撑点的核插值，可有效实现远场矩阵块的低秩近似。
• 对H-矩阵块进行QR压缩可进一步降低存储需求，且精度损失可忽略不计。
• 通过将共轭导数项移至形函数中，该实现成功处理了双层算子中的法向导数，实现了精确的牵引力计算。