[论文解读] Space-Time Least-Squares Isogeometric Method for Parabolic Problems

Monica Montardini, Matteo Negri+2 more|arXiv (Cornell University)|Sep 26, 2018

Advanced Numerical Analysis Techniques被引用 1

一句话总结

该论文提出了一种基于高阶光滑样条的时空最小二乘等几何方法，用于抛物型问题，利用张量积结构设计了一种鲁棒的预条件子，通过快速对角化求解类似西尔维斯特方程的方程。该方法确保计算成本几乎与自由度成正比，且与多项式次数无关，对于高阶等几何离散化具有高度高效性和可并行性。

ABSTRACT

In this paper, we propose a space-time least-squares isogeometric method to solve parabolic evolution problems, well suited for high-degree smooth splines in the space-time domain. We focus on the linear solver and its computational efficiency: thanks to the proposed formulation and to the tensor-product construction of space-time splines, we can design a preconditioner whose application requires the solution of a Sylvester-like equation, which is performed efficiently by the fast diagonalization method. The preconditioner is robust w.r.t. spline degree and mesh size. The computational time required for its application, for a serial execution, is almost proportional to the number of degrees-of-freedom and independent of the polynomial degree. The proposed approach is also well-suited for parallelization.

研究动机与目标

• 解决高阶等几何方法在抛物型问题中因线性系统病态而导致的计算低效问题。
• 设计一种时空弱形式，利用空间和时间上光滑的高阶样条实现高精度解。
• 设计一种鲁棒且高效的预条件子，确保在不同多项式次数和网格尺寸下均保持性能。
• 确保预条件子应用成本几乎与自由度成线性关系，且独立于样条次数。
• 通过张量积结构实现解算过程的高效并行化，适用于大规模问题。

提出的方法

• 在时空域中使用高阶B样条的等几何分析方法对抛物型问题进行公式化。
• 对时空弱形式应用最小二乘变分公式，确保系统矩阵的对称性和正定性。
• 利用时空样条的张量积结构，将系统分解为类似西尔维斯特的矩阵方程。
• 基于张量结构设计预条件子，将线性系统求解简化为求解西尔维斯特方程。
• 应用快速对角化方法高效求解西尔维斯特方程，利用空间和时间算子的特征值分解。
• 确保预条件子在多项式次数或网格加密增加时仍保持有效性和计算效率。

实验结果

研究问题

• RQ1 能否通过时空最小二乘等几何方法在使用高阶光滑样条时，对抛物型问题实现鲁棒性能？
• RQ2 如何在保持高阶精度的同时最小化求解所得线性系统的计算成本？
• RQ3 所提出的预条件子在不同多项式次数和网格尺寸下是否仍保持有效性和高效性？
• RQ4 能否利用张量积结构高效并行化求解过程？
• RQ5 预条件子应用的可扩展性和计算复杂度在自由度和多项式次数方面如何？

主要发现

• 所提方法在空间和时间上均使用高阶光滑样条，实现了高阶精度，具备优异的逼近特性。
• 预条件子应用成本几乎与自由度成线性关系，且独立于多项式次数。
• 快速对角化方法高效求解了由张量积结构产生的类似西尔维斯特方程。
• 预条件子对网格加密和样条次数增加均具有鲁棒性，在广泛离散化范围内保持高效。
• 由于时空系统的可分结构，该方法非常适合并行执行。
• 计算实验表明，求解器性能接近最优，且随着多项式次数增加，性能退化极小。