[論文レビュー] A two-grid method for the $C^0$ interior penalty discretization of the Monge-Ampère equation

Gerard Awanou, Hengguang Li+1 more|arXiv (Cornell University)|Apr 12, 2018

Nonlinear Partial Differential Equations参考文献 2被引用数 1

ひとこと要約

本稿では、Monge-Ampère方程式のC⁰内部罰則有限要素離散化から生じる非線形連立一次方程式を効率的に解くための二重グリッド法を提案する。粗いメッシュ上で非線形連立一次方程式を解き、その解を細かいメッシュ上で1回のニュートン反復の初期推定値として用いることで、標準的な細かいメッシュ上のニュートン法と比較して計算コストを著しく削減するとともに、準最適なW¹,∞誤差推定値を達成する。数値結果により、その効率性と最適収束率が確認されている。

ABSTRACT

No description supplied

研究の動機と目的

• 楕円型Monge-Ampère方程式のC⁰内部罰則離散化から生じる非線形連立一次方程式を解くための効率的な数値法を開発すること。
• 粗いメッシュ上の解と細かいメッシュ上で1ステップのニュートン反復を組み合わせた二重グリッドスキームの収束性と誤差推定値を分析すること。
• 3次以上の有限要素を用いたC⁰内部罰則法に対して、準最適なW¹,∞誤差推定値を確立すること。
• 数値実験を通じて、二重グリッド法と細かいメッシュ上の標準的ニュートン法との比較により、計算効率を示すこと。

提案手法

• 本手法は二重グリッド戦略を用いる：粗いメッシュ（サイズH）上で非線形離散連立一次方程式を解き、その解を細かいメッシュ（サイズh）上で1回のニュートン反復の初期推定値として用いる。
• 粗いメッシュ上の解は、相対残差の許容誤差10⁻⁶を満たすまでニュートン法により得られる。
• 細かいメッシュ上のニュートンステップは1回のみ実行され、完全なニュートン反復と比較して計算コストを著しく削減する。
• 本手法は、離散双一次形式のフレシェ微分構造を活用しており、ヘッセ行列の余因子行列に対応する線形楕円型作用素に対応する。
• 誤差解析は、逆推定、トレース不等式、離散ソボレフ埋め込みを用いて、二重グリッド解と細かいメッシュ解の差を評価する。
• 収束率は、H = h^λ かつ λ > max{(k+2+ε)/(2k), 3/(2k−1)} を満たすものと仮定し、k ≥3 の場合に最適なH¹誤差収束を保証する。

実験結果

リサーチクエスチョン

• RQ1 C⁰内部罰則離散化から生じるMonge-Ampère方程式の非線形連立一次方程式に、二重グリッド法を効果的に適用可能か？
• RQ2 高次有限要素を用いた場合、二重グリッド解のW¹,∞ノルムにおける収束速度と誤差推定値は何か？
• RQ3 二重グリッド法の計算コストは、細かいメッシュ上の標準的ニュートン法と比べてどの程度か？
• RQ4 二重グリッドアプローチは、元の有限要素離散化の最適収束率を保持するか？
• RQ5 粗いメッシュから細かいメッシュへのメッシュ比λの変化が、手法の精度と効率に与える影響は何か？

主な発見

• 二重グリッド法は、3次以上の有限要素を用いたC⁰内部罰則離散化に対して、準最適なW¹,∞誤差推定値を達成する。
• H = h^λ かつ λ > max{(k+2+ε)/(2k), 3/(2k−1)} を満たすとき、k ≥3 の場合に二重グリッド解のH¹誤差は最適なO(h^k)の収束率を示す。
• 二重グリッド法の計算時間は、細かいメッシュ上のニュートン法と比較して著しく低く、すべてのメッシュ細分化において高速化が観察された。
• 数値実験により、P2要素ではH¹ノルムで2次収束、P3要素では最適収束が確認されたが、理論はk ≥3 の場合に焦点を当てている。
• 二重グリッド法は最適収束率を維持しており、細かいグリッド上の完全なニュートン反復と比較して計算的にも優れており、収束率に損なわれることはない。
• 細かいメッシュ上で2回目のニュートンステップを実行すると精度がわずかに向上するが、収束率に影響を与えないことから、1回の細かいメッシュ補正で既に最適な性能が達成されていることが示された。