[论文解读] A Hybridized Weak Galerkin Finite Element Method for the Biharmonic Equation

Chunmei Wang, Junping Wang|arXiv (Cornell University)|Feb 5, 2014

Numerical methods in engineering被引用 23

一句话总结

本文通过在单元边界引入拉格朗日乘子以放松连续性约束，提出了一种用于双调和方程的混合型弱伽辽金有限元方法（HWG）。该方法实现了最优误差估计，并可实现舒尔补约化，显著减少了通过消除单元内部自由度而得到的系统规模，从而获得对称正定的降阶系统，实现高效求解算法。

ABSTRACT

This paper presents a hybridized formulation for the weak Galerkin finite element method for the biharmonic equation. The hybridized weak Galerkin scheme is based on the use of a Lagrange multiplier defined on the element boundaries. The Lagrange multiplier is verified to provide a numerical approximation for certain derivatives of the exact solution. An optimal order error estimate is established for the numerical approximations arising from the hybridized weak Galerkin finite element method. The paper also derives a computational algorithm (Schur complement) by eliminating all the unknown variables on each element, yielding a significantly reduced system of linear equations for unknowns on the boundary of each element.

研究动机与目标

• 开发双调和方程的弱伽辽金有限元方法的混合化形式，以提高计算效率。
• 在单元边界引入拉格朗日乘子，以近似精确解的导数。
• 在HWG框架下建立数值近似解的最优阶误差估计。
• 通过消除内部未知量，推导出舒尔补系统，从而减小全局系统规模。
• 通过变量降阶与对称正定的降阶线性系统，实现高效且可并行计算。

提出的方法

• 采用多项式近似 $P_k/P_{k-2}/P_{k-2}$ 的弱伽辽金有限元格式，分别用于解、其迹以及单元边界上的梯度。
• 在单元界面引入拉格朗日乘子，以弱形式强制连续性，从而实现混合化。
• 采用舒尔补技术消除每个单元上的所有内部未知量，将全局系统简化为仅包含边界未知量的形式。
• 推导出仅涉及边界变量的降阶线性系统，该系统具有对称正定性，从而可高效使用迭代求解器。
• 通过两步算法求解全局系统：首先求解降阶的舒尔补系统，然后逐单元重构内部值。
• 利用离散弱海森矩阵与局部投影算子定义变分格式，并确保稳定性。

实验结果

研究问题

• RQ1 混合化方法能否有效应用于双调和方程的弱伽辽金有限元方法，以降低计算成本？
• RQ2 在单元边界引入拉格朗日乘子是否能提供解及其导数的稳定且收敛的数值近似？
• RQ3 舒尔补技术能否成功应用于混合型弱伽辽金格式，从而获得未知量显著减少的降阶系统？
• RQ4 混合型弱伽辽金方法求解双调和方程的最优收敛速率是多少？
• RQ5 在实际应用中，如何高效求解所得降阶系统，特别是采用预条件技术时？

主要发现

• 混合型弱伽辽金方法对解及其导数均实现了最优阶误差估计，验证了理论收敛速率。
• 单元边界上的拉格朗日乘子为精确解的法向导数提供了数值近似，从而能够准确恢复梯度信息。
• 舒尔补系统将全局未知量数量减少至仅内部和外部边界的未知量，显著降低计算成本。
• 降阶系统具有对称正定性，确保了稳定性，并可使用高效迭代求解器配合预条件技术。
• 计算算法允许在每个单元上并行求解局部问题，显著提升了可扩展性与性能。
• 该方法在保持原始弱伽辽金格式一致性和稳定性的同时，通过混合化显著提升了计算效率。