[论文解读] Fast integral equation methods for linear and semilinear heat equations in moving domains

Jun Wang, Leslie Greengard+2 more|arXiv (Cornell University)|Oct 2, 2019

Electromagnetic Scattering and Analysis参考文献 75被引用 5

一句话总结

本论文提出了一种快速、高阶、时空自适应的积分方程方法，用于求解移动或复杂静止区域中的线性和半线性热方程。通过结合基于快速高斯变换的热势求解器与自适应网格加密及高阶积分法，该方法在无需矩阵求逆的情况下实现了稳定性和高精度，从而即使在刚性问题中也能实现显式时间推进。

ABSTRACT

We present a family of integral equation-based solvers for the linear or semilinear heat equation in complicated moving (or stationary) geometries. This approach has significant advantages over more standard finite element or finite difference methods in terms of accuracy, stability and space-time adaptivity. In order to be practical, however, a number of technical capabilites are required: fast algorithms for the evaluation of heat potentials, high-order accurate quadratures for singular and weakly integrals over space-time domains, and robust automatic mesh refinement and coarsening capabilities. We describe all of these components and illustrate the performance of the approach with numerical examples in two space dimensions.

研究动机与目标

• 解决传统有限差分/有限元方法在复杂时变几何区域中求解热方程时的局限性。
• 通过使用避免矩阵求逆的积分方程公式，克服显式格式的稳定性限制。
• 为具有任意几何形状的移动区域中的热势开发高阶、自适应求解器。
• 实现对具有复杂边界条件的二维区域中线性和半线性热方程的高效、稳定且精确模拟。
• 为将该方法扩展至扩散方程组及变系数扩散问题提供基础。

提出的方法

• 利用抛物型势理论，通过初始、体积和层热势来表述热方程的解。
• 采用快速高斯变换（FGT），结合分层聚类与自适应加密，以快速评估热势。
• 使用层级受限四叉树数据结构，实现高效的时空分区与局部加密。
• 对时空边界上的奇异积分和弱奇异积分应用高阶积分规则（例如16阶勒让德多项式）。
• 采用带有自适应时间步长和时空网格加密的预测-校正时间推进格式。
• 将快速FGT与边界积分方程相结合，高效计算狄利克雷、诺伊曼和罗宾边界条件。

实验结果

研究问题

• RQ1 快速积分方程方法是否能在移动或复杂几何形状的热方程求解中实现高阶精度与稳定性？
• RQ2 如何将大规模粒子系统中热势评估的计算成本降低至O(N)复杂度？
• RQ3 该方法在处理具有刚性反应项的半线性项时表现如何？
• RQ4 自适应时空网格策略在演化区域中对精度与效率有何影响？
• RQ5 该方法能否扩展至处理二维区域中变系数或分段常数扩散系数的问题？

主要发现

• 该方法在时间方向达到二阶精度，空间方向通过16阶勒让德多项式积分规则实现高阶精度。
• 对于包含73条闭合曲线的周期性方盒，T = 1000Δt时的相对L²误差为2.7 × 10⁻⁶，且Δt ≈ Δx²。
• 每时间步的CPU时间为17.6秒，其中10.8秒用于非层级热势，6.27秒用于层势。
• 在p=2的富江模型下，周期性边界条件下表现出有限时间爆破行为，其行为与常微分方程解1/(1/u_c - (t - t_c))高度吻合。
• 在周期性方盒情况下，解最初因扩散而衰减，之后在t_c ≈ 0.099后趋于近似恒定，与常微分方程模型一致。
• 该方法在处理具有复杂边界几何形状和高分辨率离散化的多连通静态区域时表现出鲁棒性与高效性。