[論文レビュー] Finite element methods respecting the discrete maximum principle for convection-diffusion equations

Gabriel R. Barrenechea, Volker John+1 more|arXiv (Cornell University)|Apr 15, 2022

Advanced Numerical Methods in Computational Mathematics被引用数 6

ひとこと要約

この論文は、対流優勢な状態における対流拡散反応方程式に対して、離散的最大原理（DMP）を満たす有限要素法についてサーベイする。代数的安定化スキームおよびFEM-FCT や AFC といった非線形手法が、現在入手可能な唯一の方法であり、DMP を同時に満たし、高精度な解を提供することを特定している。また、単調性を保ちつつ高次元離散化における最近の進展にも言及している。

ABSTRACT

Convection-diffusion-reaction equations model the conservation of scalar quantities. From the analytic point of view, solution of these equations satisfy under certain conditions maximum principles, which represent physical bounds of the solution. That the same bounds are respected by numerical approximations of the solution is often of utmost importance in practice. The mathematical formulation of this property, which contributes to the physical consistency of a method, is called Discrete Maximum Principle (DMP). In many applications, convection dominates diffusion by several orders of magnitude. It is well known that standard discretizations typically do not satisfy the DMP in this convection-dominated regime. In fact, in this case, it turns out to be a challenging problem to construct discretizations that, on the one hand, respect the DMP and, on the other hand, compute accurate solutions. This paper presents a survey on finite element methods, with a main focus on the convection-dominated regime, that satisfy a local or a global DMP. The concepts of the underlying numerical analysis are discussed. The survey reveals that for the steady-state problem there are only a few discretizations, all of them nonlinear, that at the same time satisfy the DMP and compute reasonably accurate solutions, e.g., algebraically stabilized schemes. Moreover, most of these discretizations have been developed in recent years, showing the enormous progress that has been achieved lately. Methods based on algebraic stabilization, nonlinear and linear ones, are currently as well the only finite element methods that combine the satisfaction of the global DMP and accurate numerical results for the evolutionary equations in the convection-dominated situation.

研究の動機と目的

• 対流優勢な対流拡散問題において、離散的最大原理（DMP）を厳密に満たす有限要素法を特定・分析すること。
• 安定化有限要素スキームにおけるDMP満たしと数値精度のトレードオフを評価すること。
• 線形と非線形安定化技術の両方が、解の物理的境界を保つ上で果たす役割を評価すること。
• 定常状態および時間依存問題の両方における、代数的安定化およびフラックス補正有限要素法の最近の進展を要約すること。
• 正の保存性を保証し、不自然な振動を回避する手法について包括的な概要を提供すること。

提案手法

• 代数的安定化（人工拡散および局所射影法を含む）を通じてDMPを強制する有限要素法に焦点を当てる。
• 定常状態問題における非線形スキーム（Mizukami–Hughes法、Burman–Ern法、代数的フラックス補正（AFC））を分析する。
• 時間依存問題におけるFEM-FCTおよび単調性保証型局所射影安定化（LPS）法を検討し、DMPを保持する能力に注目する。
• 上流差分法やエッジ平均化有限要素法といった線形手法をレビューし、対流優勢下でのDMP満たしの限界を指摘する。
• 制限子と勾配近似の役割を評価し、単調性と線形保存性を保ちながら精度を向上させる点を検討する。
• Q₁、高次元H¹-適合、Crouzeix–Raviart、および非連続ガラーキン要素を含む、さまざまな有限要素タイプの比較を行う。

実験結果

リサーチクエスチョン

• RQ1 対流優勢状態における対流拡散方程式に対して、どの有限要素法が離散的最大原理（DMP）を厳密に満たすことができるか？
• RQ2 代数的安定化スキームは、古典的安定化手法と比較して、DMP満たしと解の精度の両面でどのように異なるか？
• RQ3 フラックス補正や局所射影といった非線形安定化メカニズムが、単調性を保ち、振動を回避する上で果たす役割は何か？
• RQ4 線形有限要素法は、高次元精度を維持しつつ、グローバルDMPを満たすように修正可能か？
• RQ5 3次元におけるDMPを保つ手法の実装における実用的トレードオフ、特に計算コストと非線形ソルバーの効率性についての考察は？

主な発見

• 定常状態の手法においては、非線形スキーム、特に代数的安定化およびフラックス補正手法のみが、DMPを満たし、同時に高精度な解を提供できる。
• Mizukami–Hughes法は、唯一、DMPを満たすことが証明された古典的安定化手法であり、他の多くは対流優勢下で失敗する。
• 最近の代数的安定化スキーム（修正LPSや勾配ベース制限子を含む）は、過度なスメアリングを回避しつつ、テラス効果を低減し、精度が向上している。
• 時間依存問題においては、代数的安定化FEM-FCTおよびAFCスキームが、現在、グローバルDMPを保ちつつ良好な精度を維持できる唯一の有限要素法である。
• 代数的安定化の線形版は、精度と計算効率のバランスに優れており、実用的応用に適している。
• 半離散スキームに単調性保証型局所射影作用素を追加することで、テラス効果が著しく低減され、解の品質が向上するが、DMPは損なわれない。