[论文解读] Test Function Spaces for Geometric Finite Elements

Oliver Sander|arXiv (Cornell University)|Jul 25, 2016

Advanced Numerical Analysis Techniques参考文献 11被引用 4

一句话总结

本文提出了几何有限元（GFE）的测试函数空间，使非线性流形值问题中的弱形式化和最优性条件成为可能，超越了最小化方法。该方法通过利用局部流形结构对GFE函数进行变分构造测试函数，推广了经典微分几何中的Jacobi场，实现了代数形式与几何形式的一阶最优性条件之间的等价性。

ABSTRACT

We construct test function spaces for geometric finite elements. Geometric finite elements (GFE) are generalizations of Lagrangian finite elements to situations where the unknown function maps into a nonlinear space. Test functions for such spaces arise as variations of GFE functions wherever the GFE function space has a local manifold structure. For any given GFE function $u_h$, the test functions form a linear space that depends on $u_h$. They generalize Jacobi fields in the same way that the GFE interpolation functions generalize geodesic curves. Having test function spaces allows to extend the GFE method to boundary value problems that do not have a minimization formulation.

研究动机与目标

• 通过构造适当的测试函数空间，将几何有限元（GFE）从基于最小化的公式推广至更广范围。
• 将测试函数定义为GFE函数的变分，推广经典微分几何中的Jacobi场。
• 在GFE函数空间中直接建立一阶最优性条件，而非仅在代数设定下。
• 建立通过代数最小化导出的最优性条件与通过几何测试函数导出的最优性条件之间的等价性。
• 为评估测试函数提供计算上可行的框架，尤其适用于测地线和基于投影的GFE。

提出的方法

• 测试函数被定义为给定GFE函数 $ u_h $ 的变分，形成依赖于 $ u_h $ 的线性空间 $ T_{u_h}V_h^M $。
• 通过在拉格朗日节点处使用每节点 $ \text{dim}(M) $ 维向量的节点基，将空间 $ T_{u_h}V_h^M $ 识别为切向量。
• 对于测地线有限元，计算测试函数在某一点的值仅需求解一个大小为 $ \dim M $ 的线性系统，从而实现高效计算。
• 对于基于投影的GFE，计算依赖于投影机制；对于球面，存在闭式表达式。
• 该构造对测地线和基于投影的GFE统一适用，仅在插值细节上有所不同。
• 通过涉及求值映射与切映射的交换图的可交换性，证明了代数与几何最优性条件之间的等价性。

实验结果

研究问题

• RQ1 如何为几何有限元定义测试函数空间，以实现在非线性流形值问题中的弱形式化？
• RQ2 离散GFE函数空间的切空间与代数配置流形 $ M^n $ 的切空间之间存在何种几何与代数关系？
• RQ3 能否在GFE函数空间中直接表述最小化问题的一阶最优性条件，而非仅在代数设定下？
• RQ4 通过代数最小化导出的最优性条件是否与通过几何测试函数导出的最优性条件等价？
• RQ5 在GFE中评估测试函数的计算成本是多少？能否实现高效的实际应用？

主要发现

• 几何有限元的测试函数被构造为离散函数 $ u_h $ 的变分，形成依赖于 $ u_h $ 的线性空间 $ T_{u_h}V_h^M $，推广了Jacobi场。
• 空间 $ T_{u_h}V_h^M $ 同构于拉格朗日节点处切空间的乘积，其节点基通过在每个节点指定切向量来定义。
• 对于测地线有限元，计算测试函数在某一点的值仅需求解一个大小为 $ \dim M $ 的线性系统，计算效率高。
• 对于球面 $ M = S^{N-1} $，测试函数存在闭式表达式，可实现快速计算。
• 一阶最优性条件 $ d\mathcal{J}[v_h](\eta_h) = 0 $ 对所有 $ \eta_h \in T_{v_h}V_h^M $ 成立，等价于代数条件 $ dJ[\bar{v}](\bar{\eta}) = 0 $，证明了几何与代数公式的自洽性。
• 所提出的测试函数构造方法为无最小化结构的问题提供了理论依据，扩展了GFE方法的应用范围。