Convergence of finite element right-hand-side computation from finite difference data

技术摘要：有限差分数据下有限元右端项计算的收敛性研究

问题陈述
在混合计算声学（CAA）和耦合场问题中，需要在结构化有限差分（FD）网格（通常用于流场模拟）与非结构化有限元法（FEM）网格（用于声传播）之间进行精确的场传递。核心挑战在于组装有限元右端项（RHS）向量，其定义为 $b_i = \int_{\Omega} N_i(x)f(x) d\Omega$，其中 $f(x)$ 是定义在 FD 网格上的源项场。

标准方法通常依赖于简单的插值，这可能会引入显著误差，特别是当源项表现出剧烈变化或强空间振荡（例如 Lighthill 类源项）时。本文研究了旨在确保精度、能量守恒（全局和局部）以及随网格分辨率提高而收敛的积分方法，且该方法不受源场平滑性假设的限制。

方法论
本文提出并分析了两种不同的将 FD 数据传递至 FEM RHS 的积分策略：

1. 基于场插值的高阶高斯积分法：

   • 机制： 该方法通过在每个 FEM 单元内使用高斯积分来评估积分。首先利用高阶 B-样条（阶数 $p \in \{1, \dots, 5\}$）或拉格朗日多项式（阶数 $p \in \{1, 3\}$）将 FD 场值插值到高斯积分点上。
   • 误差分析： 总误差被分解为插值误差（$\varepsilon_{interp}$）、积分误差（$\varepsilon_{quad}$）和 FEM 离散误差（$\varepsilon_{fem}$）。研究表明，为了实现收敛，插值阶数和积分阶数必须保持一致地增加。如果仅增加积分阶数而不相应增加插值阶数，插值误差将成为限制因素，从而阻碍向机器精度的收敛。

2. 基于几何裁剪的切割单元（超网格）积分法：

   • 机制： 该方法通过将问题视为两个非一致网格之间的几何交集，从而避免了网格到网格层面的插值误差。将 FD 单元和 FEM 单元叠加以创建“切割多边形”（超网格）。
   • 实现： 领域被分解为交集多边形，随后将其剖分为三角形。通过在物理空间内的子三角形上进行高斯积分来执行积分。
   • 复杂度： 一个显著的计算挑战是在物理高斯点（位于切割区域内）处评估 FEM 形函数所需的逆映射。对于非线性几何映射，这要求为每个积分点求解一个非线性系统（例如通过 Newton-Raphson 迭代）。
   • 优化： 算法利用轴对齐包围盒（AABB）索引和二分查找来高效识别相交的 FD 单元，随后使用多边形裁剪（例如 Sutherland-Hodgman 算法）来定义积分区域。

关键结果

• 收敛行为：
  • 基于插值的方法在 B-样条阶数和积分阶数增加的情况下可以收敛，前提是 FD 网格分辨率保持一致。然而，它存在一个由插值重建阶数决定的系统误差底限。对于高度振荡的源项，除非使用极高阶数，否则该方法难以达到机器精度。
  • 超网格方法消除了网格到网格界面处的插值误差和近似误差。无论源场是平滑的还是强振荡的，它都能在所有网格分辨率下收敛至机器精度。
• 源项敏感性： 使用 Lighthill 类源项（空间变化剧烈、具有尖锐过渡）和扰动对流波动方程源项（较平滑）进行的对比研究表明，超网格方法具有普遍的鲁棒性。插值类方法在处理更复杂的 Lighthill 类项时表现出更高的误差。
• 性能： 虽然超网格方法涉及更高的几何复杂度以及逆映射的计算开销，但它是经过高度优化的。对于包含大量单元的问题，其性能可以超越未优化的插值方法。两种方法的积分弱扩展性均为线性。

意义与主张
本文声称，超网格积分法为需要严格保证源项能量守恒的混合 CAA 模拟提供了一种优越的替代方案。通过消除 FD 网格与 FEM 高斯点之间的插值步骤，该方法确保了来自流场的完整能量能够传递到声学网格，而不会产生减少或人工平滑。

作者强调，虽然高阶插值可以降低误差，但如果场量不是“足够平滑”，则无法保证在实践中实现充分的收敛。因此，超网格方法被呈现为一种针对平滑及强振荡场量均通用的鲁棒解决方案。这项工作对这些积分方案进行了严谨的对比分析，强调了选择何种方法取决于源项的性质以及对严格守恒的要求。

结论
研究得出结论，对于混合声学领域中的精确场传递，当精度和守恒性至关重要时，超网格积分法是首选方法，因为它消除了插值方案中固有的系统误差。插值法在处理较平滑的源项时仍然可行，但需要仔细平衡插值阶数和积分阶数，以最小化误差。