[论文解读] Galerkin Least-Squares Stabilization in Ice Sheet Modeling - Accuracy, Robustness, and Comparison to other Techniques

Christian Helanow, Josefin Ahlkrona|arXiv (Cornell University)|Feb 27, 2017

Cryospheric studies and observations参考文献 44被引用 6

一句话总结

本研究评估了在冰盖流动有限元模拟中使用Galerkin最小二乘法（GLS）稳定化的性能，与内部惩罚法和压力惩罚法等其他方法进行比较。研究发现，GLS对稳定化参数的选择极为敏感，特别是在垂直速度方面，导致冰盖边缘附近出现振荡；相比之下，内部惩罚法在复杂几何结构中表现出更优的鲁棒性和准确性。

ABSTRACT

We investigate the accuracy and robustness of one of the most common methods used in glaciology for the discretization of the $\mathfrak{p}$-Stokes equations: equal order finite elements with Galerkin Least-Squares (GLS) stabilization. Furthermore we compare the results to other stabilized methods. We find that the vertical velocity component is more sensitive to the choice of GLS stabilization parameter than horizontal velocity. Additionally, the accuracy of the vertical velocity component is especially important since errors in this component can cause ice surface instabilities and propagate into future ice volume predictions. If the element cell size is set to the minimum edge length and the stabilization parameter is allowed to vary non-linearly with viscosity, the GLS stabilization parameter found in literature is a good choice on simple domains. However, near ice margins the standard parameter choice may result in significant oscillations in the vertical component of the surface velocity. For these cases, other stabilization techniques, such as the interior penalty method, result in better accuracy and are less sensitive to the choice of the stabilization parameter. During this work we also discovered that the manufactured solutions often used to evaluate errors in glaciology are not reliable due to high artificial surface forces at singularities. We perform our numerical experiments in both FEniCS and Elmer/Ice.

研究动机与目标

• 评估等阶有限元方法结合Galerkin最小二乘法（GLS）稳定化在冰盖模拟中的精度与鲁棒性。
• 研究速度场中的数值误差（尤其是垂直速度）如何传播至冰面位置与体积预测中。
• 在基准问题与真实冰盖问题上，将GLS稳定化与替代稳定化方法（压力惩罚法、内部惩罚法、压力全局投影法和局部投影法）进行比较。
• 评估制造解在冰川学有限元模型误差评估中的可靠性。
• 确定最优稳定化参数，并识别标准GLS公式在真实冰川学场景中的局限性。

提出的方法

• 使用两个有限元代码（FEniCS和Elmer/Ice）开展数值实验，实现结果的交叉验证。
• 在等阶线性有限元中对p-Stokes方程应用GLS稳定化，稳定化参数源自Franca和Frey（1992），取决于单元尺寸与黏度。
• 采用两种误差测量策略：使用制造解进行收敛性测试，使用高分辨率参考解进行二维验证。
• 利用ISMIP-HOM基准和真实出口冰川设置，测试在不同网格细化程度与黏度条件下的性能表现。
• 在相同条件下实现并比较替代稳定化技术（内部惩罚法、压力惩罚法、压力全局投影法和局部投影法）。
• 网格细化研究聚焦于垂直与水平方向的收敛性，特别分析了在奇点与冰盖边缘附近垂直速度分量的行为。

实验结果

研究问题

• RQ1 GLS稳定化参数的选择如何影响冰盖模型中垂直速度与表面位置的精度？
• RQ2 制造解在冰川学有限元模拟中能在多大程度上提供可靠的误差估计？
• RQ3 不同稳定化有限元方法（GLS、内部惩罚法、压力惩罚法等）在精度、鲁棒性及对稳定化参数敏感性方面如何比较？
• RQ4 垂直速度分量中的数值误差对长期冰面演化与体积预测有何影响？
• RQ5 在何种情况下GLS稳定化会导致非物理解振荡？替代方法如内部惩罚法能否缓解此类问题？

主要发现

• 垂直速度分量对GLS稳定化参数的选择比水平分量更敏感，参数选择不当会导致冰盖边缘附近出现振荡。
• 标准GLS稳定化参数（源自Franca和Frey，1992）在简单区域表现良好，当单元尺寸定义为最小边长且黏度呈非线性变化时；但在真实出口冰川模拟中表现失败。
• 制造解在奇点处产生人工表面力，导致误差估计偏高，并在垂直方向出现超收敛现象，从而削弱其在模型验证中的可靠性。
• 内部惩罚法相比GLS表现出更优的稳定性，且对稳定化参数选择的敏感性更低，尤其在具有高黏度对比度的复杂几何结构中。
• 压力全局投影法精度高，但对大规模冰盖模拟计算成本过高。
• 局部投影稳定化性能与GLS相当，但实现更复杂；其非线性变体可能表现更优，但需进一步研究。