[論文レビュー] Improving Newton's method performance by parametrization: the case of Richards equation

Konstantin Brenner, Clément Cancès|arXiv (Cornell University)|Jul 6, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 42被引用数 54

ひとこと要約

本稿では、不飽和多孔質媒体流れにおける特異な非線形系を解くニュートン法の安定化を目的として、リチャーズ方程式のパrametrizationに基づく再定式化を提案する。新たな主変数（フラックスに基づくパラメータ化）を導入することで、2次収束性と強固な質量保存性を確保し、特にメッシュの細分化や物性値の変動下でも、標準的な圧力-飽和定式化に比べて反復回数と精度の面で顕著に優れている。

ABSTRACT

The nonlinear systems obtained by discretizing degenerate parabolic equations may be hard to solve, especially with Newton's method. In this paper, we apply to Richards equation a strategy that consists in defining a new primary unknown for the continuous equation in order to stabilize Newton's method by parametrizing the graph linking the pressure and the saturation. The resulting form of Richards equation is then discretized thanks to a monotone Finite Volume scheme. We prove the well-posedness of the numerical scheme. Then we show under appropriate non-degeneracy conditions on the parametrization that Newton\\^as method converges locally and quadratically. Finally, we provide numerical evidences of the efficiency of our approach.

研究の動機と目的

• 離散化されたリチャーズ方程式から生じる非線形系を解く際のニュートン法の収束性の悪さを是正すること。
• 移動度関数および飽和関数の特異性が生じる状況下でも、ニュートン法のロバスト性と効率性を向上させること。
• 圧力-飽和関係の新しいパラメータ化を通じて、質量保存性と2次収束性を保証すること。
• 数値実験を通じて、標準的な圧力ベースの（u）定式化に比べ、フラックスベースの（τ）定式化の優位性を示すこと。

提案手法

• ニュートン反復の安定化のため、飽和関数の逆数として定義される新たな主変数 τ を導入してリチャーズ方程式を再定式化する。
• 安定性と一貫性を保証するため、単調な有限体積法をパラメータ化されたリチャーズ方程式の形に適用する。
• パラメータ化に関する適切な非特異性条件の下で、離散スキームの適切な定義（well-posedness）を証明する。
• 滑らかさおよび非特異性の仮定の下で、パラメータ化された系に対するニュートン法の局所的2次収束性を確立する。
• 各反復で質量保存性を維持できるように、線形化されたヤコビ行列に基づくラインサーチを用いた不正確ニュートン法を実装する。
• ニュートン更新サブプロブレムを解くために、パラメータ化された残差のヤコビ行列に基づく線形化システムを用いる。

実験結果

リサーチクエスチョン

• RQ1 圧力-飽和関係のパラメータ化は、リチャーズ方程式に対するニュートン法の収束性を改善できるか？
• RQ2 提案されたパラメータ化は、標準的定式化に比べてより優れた質量保存性を達成できるか？
• RQ3 新しい定式化において、ニュートン法の性能はメッシュの細分化および物性パラメータ β にどのように依存するか？
• RQ4 パラメータ化された系は、適切な定義と安定性を保ちつつ、2次収束性を維持できるか？
• RQ5 反復回数および誤差伝搬の観点から、τ-定式化とu-定式化の相対的効率性はいかほどか？

主な発見

• τ-定式化では、特にメッシュの細分化下で、u-定式化に比べて収束に必要なニュートン反復回数が最大1桁減少する。
• ε = 10⁻⁸ の場合、u-定式化では相対的飽和誤差が質量保存誤差に支配され（最大10⁻⁴）、一方τ-定式化では質量保存誤差が10⁻¹⁵未満に抑えられる。
• τ-定式化はβの値の変化に対してもロバストな性能を維持するが、u-定式化ではβが増加するにつれて反復回数が著しく増加する。
• 理論的分析により、パラメータ化の非特異性条件の下で、パラメータ化された系に対するニュートン法の局所的2次収束性が確認された。
• 数値結果から、選択されたパrameter領域においてs(τ)が線形であるため、τ-定式化はマシン精度に達する質量保存性を達成している。
• ε = 10⁻¹⁶ を用いてτ-定式化で計算された基準解は、質量誤差が無視できるほど小さく、高精度ベンチマークとしての有効性が検証された。