[论文解读] Legendre-Gauss-Lobatto grids and associated nested dyadic grids

Kolja Brix, Claudio Canuto+1 more|arXiv (Cornell University)|Oct 31, 2013

Advanced Numerical Methods in Computational Mathematics参考文献 9被引用 9

一句话总结

本论文提出了一类与勒让德-高斯-洛巴托（LGL）网格相关联的新型嵌套二元网格，解决了高阶谱方法和不连续伽辽金法中LGL节点非嵌套的关键问题。通过推导LGL网格的精细解析性质，并利用递归加密算法构建二元网格，作者证明了所得网格在局部上与LGL网格等价，且保持最优基数，从而为高阶PDE求解器提供了稳健的预条件化方法。

ABSTRACT

Legendre-Gauss-Lobatto (LGL) grids play a pivotal role in nodal spectral methods for the numerical solution of partial differential equations. They not only provide efficient high-order quadrature rules, but give also rise to norm equivalences that could eventually lead to efficient preconditioning techniques in high-order methods. Unfortunately, a serious obstruction to fully exploiting the potential of such concepts is the fact that LGL grids of different degree are not nested. This affects, on the one hand, the choice and analysis of suitable auxiliary spaces, when applying the auxiliary space method as a principal preconditioning paradigm, and, on the other hand, the efficient solution of the auxiliary problems. As a central remedy, we consider certain nested hierarchies of dyadic grids of locally comparable mesh size, that are in a certain sense properly associated with the LGL grids. Their actual suitability requires a subtle analysis of such grids which, in turn, relies on a number of refined properties of LGL grids. The central objective of this paper is to derive just these properties. This requires first revisiting properties of close relatives to LGL grids which are subsequently used to develop a refined analysis of LGL grids. These results allow us then to derive the relevant properties of the associated dyadic grids.

研究动机与目标

• 解决勒让德-高斯-洛巴托（LGL）网格非嵌套的根本性局限，该问题阻碍了高阶不连续伽辽金方法中高效预条件化的实现。
• 构建一组与LGL网格局部均匀等价的嵌套二元网格层级，以支持辅助空间方法中辅助问题的稳定高效求解。
• 建立LGL网格的严格解析性质，包括节点间距与区间大小的单调性，作为构建关联二元网格的基础。
• 证明所得二元网格相对于LGL网格保持最优基数，并满足强局部等价性界。
• 提供一种构造性算法，用于生成此类二元网格，确保嵌套性的同时保持网格质量和局部均匀性。

提出的方法

• 重新审视并推导LGL网格的精细性质，包括区间大小单调性与边界附近节点聚集特性，以切比雪夫-高斯-洛巴托节点为参考。
• 提出一种递归网格加密算法（算法2），通过仅在必要时分割区间，生成与LGL网格保持局部等价的二元网格。
• 通过条件确保局部均匀性：对于任意与LGL区间Δ相交的二元区间D，其比值|Δ|/|D|受依赖于α和Cg的常数有界。
• 利用LGL区间大小在加密下的单调性（定理4.2）证明粗网格中的区间在加密后不会增大，从而确保稳定性。
• 应用命题3验证二元网格满足所需的局部等价性与基数界，依赖于推导出的LGL节点性质。
• 引入参数α以控制加密阈值，确保仅添加最少数量的额外点以维持嵌套性。

实验结果

研究问题

• RQ1 能否构建一个嵌套的二元网格层级，使其与勒让德-高斯-洛巴托（LGL）网格局部均匀等价？
• RQ2 为确保关联的二元网格保持局部均匀性与最优基数，LGL网格所需的精细解析性质是什么？
• RQ3 在高阶不连续伽辽金法与辅助空间预条件化背景下，如何克服LGL网格的非嵌套性？
• RQ4 同一区域中二元区间与对应LGL区间的大小关系如何？能否实现统一有界？
• RQ5 能否证明所得二元网格中的点数与原始LGL网格的点数相当？

主要发现

• 由算法2生成的二元网格与LGL网格在局部(A,B)一致等价，其比值|Δ|/|D|的界依赖于α及LGL区间大小估计中的常数Cg。
• 二元网格的基数满足#D_N ≃ #G^LGL_N，即点数的增长渐近速率与LGL网格相同。
• 该构造确保仅添加最少数量的额外点即可实现嵌套性，如图7中的数值图示所证实。
• 证明依赖于LGL区间大小在加密下单调性的性质，该性质通过与切比雪夫-高斯-洛巴托节点比较得以建立。
• 二元网格满足辅助空间方法中的关键要求：网格嵌套且与LGL网格保持局部网格尺寸可比性。
• 结果证实了在高阶谱方法与不连续伽辽金法中构建稳健预条件器的理论基础。