[论文解读] Mimetic framework on curvilinear quadrilaterals of arbitrary order

Jasper Kreeft, Artur Palha+1 more|arXiv (Cornell University)|Nov 18, 2011

Numerical methods for differential equations参考文献 28被引用 47

一句话总结

本文提出了一种在曲线四边形网格上实现高阶仿射谱元框架的方法，通过离散算子与其连续对应算子之间的精确交换性，精确保持微分形式的拓扑与几何结构。通过利用代数拓扑、Hodge分解和度量相容投影，该方法确保离散算子在曲边单元上也能精确模拟连续微分算子，实现最优收敛率与稳定性。

ABSTRACT

In this paper higher order mimetic discretizations are introduced which are firmly rooted in the geometry in which the variables are defined. The paper shows how basic constructs in differential geometry have a discrete counterpart in algebraic topology. Generic maps which switch between the continuous differential forms and discrete cochains will be discussed and finally a realization of these ideas in terms of mimetic spectral elements is presented, based on projections for which operations at the finite dimensional level commute with operations at the continuous level. The two types of orientation (inner- and outer-orientation) will be introduced at the continuous level, the discrete level and the preservation of orientation will be demonstrated for the new mimetic operators. The one-to-one correspondence between the continuous formulation and the discrete algebraic topological setting, provides a characterization of the oriented discrete boundary of the domain. The Hodge decomposition at the continuous, discrete and finite dimensional level will be presented. It appears to be a main ingredient of the structure in this framework.

研究动机与目标

• 开发一种高阶仿射离散化框架，以在曲线四边形网格上保持微分形式的拓扑与几何结构。
• 通过与外微分算子和Hodge星算子交换的投影算子，建立连续微分形式与离散余链之间的一一对应关系。
• 通过在连续与离散层面显式引入对偶网格与方向性（内/外），将有限体积法与有限元法统一于同一仿射框架下。
• 通过在拉回与重构映射下保持交换关系，确保在曲边单元上的度量与拓扑一致性。
• 证明Hodge分解在连续、离散与有限维层次上均被保持，从而确保数值稳定性与最优收敛性。

提出的方法

• 基于代数拓扑的仿射框架，将微分形式在具有内在方向性（内/外）的胞复形上离散化为余链。
• 采用约化、重构与投影算子，将连续微分形式映射至离散余链，确保离散算子与连续对应算子之间精确交换。
• 应用有界线性投影（πh 与 π̃h）与余投影（π⋆h 与 π̃⋆h），在连续与离散空间之间传递函数，同时保持拓扑与度量结构。
• 通过在曲边单元上使用任意阶次多项式基函数，构建高阶仿射谱元方法，投影基于微分算子（D 与 D̃）。
• 通过度量相容投影引入离散Hodge-⋆算子，确保离散内积与楔积与连续情形一致。
• 证明微分形式的拉回运算与余链映射及离散算子之间交换，从而在曲边网格上实现结构的精确保持。

实验结果

研究问题

• RQ1 如何在曲线四边形单元上构建高阶仿射离散化方法，同时保持微分形式的拓扑与几何结构？
• RQ2 在高阶谱元框架中，离散算子（外微分、Hodge-⋆、共微分）与连续对应算子之间的交换性在多大程度上可以保持？
• RQ3 在该仿射框架中，Hodge分解如何在连续、离散与有限维层次上被保持？
• RQ4 对偶网格与方向性（内/外）在确保曲线网格上仿射离散化的协调性与稳定性方面起到何种作用？
• RQ5 该框架能否确保在曲边映射下，依赖度量的运算（如内积、Hodge星）保持一致性，同时维持拓扑交换关系？

主要发现

• 通过精心构造的投影，该仿射框架确保离散算子（如外微分与Hodge-⋆算子）与连续对应算子之间精确交换。
• Hodge分解在所有层次——连续、离散与有限维——均被保持，为数值稳定性与最优收敛性提供了结构性基础。
• 通过在非线性几何映射下保持拉回与余链映射同微分算子的交换性，该框架在曲边单元上保持了拓扑一致性。
• 对偶网格与方向性（内/外）的使用实现了本构关系与分部积分的物理解释一致处理，解决了交错有限体积方法中常见的问题。
• 通过度量相容投影构造离散内积与楔积，确保离散Hodge-⋆算子与连续情形一致。
• 在数值实验中，该方法实现了最优收敛率，误差按预期随高阶谱元方法缩放，归因于微分结构的精确保持。