[論文レビュー] 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

ひとこと要約

本研究では、アイスシート流れの有限要素モデルにおけるガレルキン最小自乗（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.

研究の動機と目的

• 等位相の有限要素法にガレルキン最小自乗（GLS）安定化を適用した場合の、アイスシートシミュレーションにおける精度とロバスト性を評価すること。
• 特に垂直速度に生じる数値誤差が、氷表面位置および体積予測にどのように伝播するかを調査すること。
• ベンチマーク問題および現実的なアイスシート問題において、GLS安定化と代替安定化手法（圧力ペナルティ法、内部ペナルティ法、圧力グローバル射影法、局所射影法）を同一条件下で比較すること。
• 流氷学的有限要素モデルにおける誤差評価のためのマニュフェイクドソリューションの信頼性を評価すること。
• 現実的な流氷学的状況下における、標準GLS定式化の最適安定化パラメータの特定と限界の同定すること。

提案手法

• FEniCSおよびElmer/Iceの2つの有限要素コードを用いて数値実験を実施し、結果の相互検証を可能にした。
• 等位相の線形有限要素を用いてp-Stokes方程式にGLS安定化を適用し、安定化パラメータはFrancaとFrey（1992）の式に従い、要素サイズと粘性度に依存する。
• 2つの誤差評価戦略を採用：収束テストにはマニュフェイクドソリューションを、2次元における検証には高解像度の基準解を用いた。
• 性能評価にはISMIP-HOMベンチマークおよび現実的なアウトレットガラスパーの設定を用い、メッシュの細分化および粘性度の変動を想定した条件でテストを行った。
• 同一条件下で、内部ペナルティ法、圧力ペナルティ法、圧力グローバル射影法、局所射影法といった代替安定化手法を実装し、比較した。
• メッシュの細分化に関する研究では、垂直および水平方向の収束を重点的に分析し、特に関節や氷の端縁付近における垂直速度成分の挙動に注目した。

実験結果

リサーチクエスチョン

• RQ1 GLS安定化パラメータの選択が、アイスシートモデルにおける垂直速度および表面位置の精度にどのように影響するか？
• RQ2 マニュフェイクドソリューションは、流氷学的有限要素シミュレーションにおける誤差推定に信頼性があるか？
• RQ3 GLS、内部ペナルティ法、圧力ペナルティ法などの異なる安定化有限要素法は、精度、ロバスト性、安定化パラメータへの感受性の観点からどのように比較されるか？
• RQ4 垂直速度成分における数値誤差が、長期的な氷表面の進化および体積予測に与える影響はどの程度か？
• RQ5 どのような状況でGLS安定化が物理的に非実情な振動を引き起こし、内部ペナルティ法などの代替手法がこれらの問題を緩和できるか？

主な発見

• 垂直速度成分は水平成分よりもGLS安定化パラメータの選択に著しく敏感であり、パラメータが不適切な場合、氷の端縁付近で振動が生じる。
• FrancaとFrey（1992）が提示した標準的なGLS安定化パラメータは、最小辺長を要素サイズとして定義し、粘性度が非線形に変化する単純な領域では良好に機能するが、現実的なアウトレットガラスパーのシミュレーションでは失敗する。
• マニュフェイクドソリューションは特異点に人工的な表面力を取り入れ、垂直方向に誤差推定値が著しく高くなり、超収束を示すことがある。このため、モデル検証の信頼性に欠ける。
• 内部ペナルティ法は、特に高粘性度対比が大きい複雑な幾何形状において、GLS法よりも優れた安定性とパラメータ選択への感受性の低さを示した。
• 圧力グローバル射影法は精度に優れるが、大規模なアイスシートシミュレーションには計算コストが高すぎて現実的ではない。
• 局所射影安定化法はGLS法と同等の性能を示すが、実装がより複雑であり、非線形版は性能向上の可能性を秘めているが、今後の研究が求められる。