Skip to main content
QUICK REVIEW

[论文解读] A discrete Weber inequality on three-dimensional hybrid spaces with application to the HHO approximation of magnetostatics

Florent Chave, Daniele A. Di Pietro|arXiv (Cornell University)|Jul 7, 2020
Advanced Numerical Methods in Computational Mathematics参考文献 33被引用 15
一句话总结

本文在多面体网格上的三维混合多项式空间上建立了离散的韦伯不等式,使能够对磁静力学的两种混合高阶(HHO)方法进行收敛性分析。这两种方法在能量误差范数下均实现了最优的 $k+1$ 阶收敛率,且在立方体和四面体网格上的数值验证结果与理论预测一致。

ABSTRACT

We prove a discrete version of the first Weber inequality on three-dimensional hybrid spaces spanned by vectors of polynomials attached to the elements and faces of a polyhedral mesh. We then introduce two Hybrid High-Order methods for the approximation of the magnetostatics model, in both its (first-order) field and (second-order) vector potential formulations. These methods are applicable on general polyhedral meshes, and allow for arbitrary orders of approximation. Leveraging the previously established discrete Weber inequality, we perform a comprehensive analysis of the two methods. We finally validate them on a set of test-cases.

研究动机与目标

  • 在多面体网格上定义的三维混合多项式空间上,建立第一个离散韦伯不等式的版本。
  • 开发并分析两种用于磁静力学问题的混合高阶(HHO)方法:一种采用一阶场形式,另一种采用二阶矢量势形式。
  • 证明两种方法在能量范数和 $L^2$-类范数下均达到最优的 $k+1$ 阶收敛率,且对任意多项式阶数 $k \geq 0$ 均成立。
  • 在立方体和正四面体网格上通过数值实验验证方法,确认理论收敛率。

提出的方法

  • 通过每个单元上多项式向量函数的分解(引理2),在由网格单元和面上的多项式向量张成的混合空间上推导出离散韦伯不等式。
  • 构建两种 HHO 方法:一种用于一阶场形式,另一种用于磁静力学的广义矢量势形式。
  • 利用局部重构和杂交化方法,通过施密特补(Schur complement)消除单元自由度,从而高效求解降阶后的全局系统。
  • 利用离散 inf-sup 条件和新证明的离散韦伯不等式,建立方法的稳定性和收敛性。
  • 采用类似混合有限元的框架,引入拉格朗日乘子,以在离散设置中保证一致性和相容性。
  • 在单位立方体网格上进行数值实验,使用立方体和四面体单元,采用 SparseLU 求解降阶系统,并在能量范数、$L^2$ 范数和 $L^2$-类范数下测量收敛性。

实验结果

研究问题

  • RQ1能否在一般多面体网格上,于三维混合多项式空间上建立类似于连续情形的离散韦伯不等式?
  • RQ2所提出的用于磁静力学的 HHO 方法是否在能量范数下对场形式和矢量势形式均实现 $k+1$ 阶最优收敛率?
  • RQ3在四面体网格上,移除离散双线性形式 $d_h$ 对 HHO 方法的精度和条件数有何影响?
  • RQ4是否可以仅通过离散韦伯不等式证明稳定性与收敛性,而无需依赖离散精确复形或全局协调性?
  • RQ5在不同网格类型和多项式阶数下,HHO 方法的实际性能如何?其收敛率和求解时间表现如何?

主要发现

  • 在多面体网格上的混合多项式空间上,严格证明了离散韦伯不等式,为收敛性分析奠定了基础。
  • 用于一阶场形式的 HHO 方法对所有 $k \geq 0$ 均实现了最优的 $k+1$ 阶能量误差收敛率。
  • 用于矢量势形式的 HHO 方法同样实现了最优的 $k+1$ 阶能量误差收敛率,且磁感应强度的 $L^2$ 误差收敛率为 $k+2$ 阶。
  • 在匹配的四面体网格上,移除双线性形式 $d_h$ 会使磁感应强度的近似略为精确,同时保持 $k+1$ 阶收敛率。
  • 数值结果在立方体和四面体网格上均确认了预测的收敛率,且求解时间与自由度数成适当比例。
  • 该方法在多项式阶数 $k = 0, 1, 2$ 下均表现稳健,误差度量保持良好条件性,并与理论预期一致。

更好的研究,从现在开始

从阅读论文到最终审阅,大幅缩短您的研究时间。

无需绑定信用卡

本解读由 AI 生成,并经人工编辑审核。