[論文レビュー] A WENO-type slope-limiter for a family of piecewise polynomial methods

Darren Engwirda, Maxwell Kelley|arXiv (Cornell University)|Jun 27, 2016

Numerical methods for differential equations参考文献 21被引用数 12

ひとこと要約

本稿では、特徴的な滑らかさと単調性制限付き再構成の非線形結合を用いて、衝撃近傍での振動を抑える高次精度を実現するWENO型勾配制限子を、区分的放物線法（PPM）および区分的四次多項式法（PQM）に導入する。この手法により、滑らかな問題においてPPMとPQMはそれぞれ3次および5次収束を達成でき、従来の単調制限子と比較して数値的拡산を顕著に低減する。

ABSTRACT

A new, high-order slope-limiting procedure for the Piecewise Parabolic Method (PPM) and the Piecewise Quartic Method (PQM) is described. Following a Weighted Essentially Non-Oscillatory (WENO)-type paradigm, the proposed slope-limiter seeks to reconstruct smooth, non-oscillatory piecewise polynomial profiles as a non-linear combination of the natural and monotone-limited PPM and PQM interpolants. Compared to existing monotone slope-limiting techniques, this new strategy is designed to improve accuracy at smooth extrema, while controlling spurious oscillations in the neighbourhood of sharp features. Using the new slope-limited PPM and PQM interpolants, a high-order accurate Arbitrary-Lagrangian-Eulerian framework for advection-dominated flows is constructed, and its effectiveness is examined using a series of one- and two-dimensional benchmark cases. It is shown that the new WENO-type slope-limiting techniques offer a significant improvement in accuracy compared to existing strategies, allowing the PPM- and PQM- based schemes to achieve fully third- and fifth-order accurate convergence, respectively, for sufficiently smooth problems.

研究の動機と目的

• 大気および海洋モデルに用いられる高次有限体積スキームにおける数値的拡散を低減すること。
• 滑らかな極値付近で単調制限子が局所的極値を不自然に平坦化する精度低下問題を解決すること。
• 局所的極値を保存しつつ単調性を維持する高次非振動的再構成戦略を開発すること。
• ALEおよび半ラグランジュフレームワークにおいて、WENOに類似した非線形重み付けを区分的多項式法（PPMおよびPQM）に拡張すること。
• スムーズな問題においてPPMおよびPQMが完全に3次および5次収束を達成できるように、不適切な斜温混合（diapycnal mixing）を最小限に抑えること。

提案手法

• 自然再構成と単調制限付きPPM/PQM補間の非線形重み付けを用いたWENO型非線形再構成を提案する。
• 再構成多項式の形状を制御するためのターゲット曲率点位置を用いることで、単調性を保ちながら振動を回避する。
• 多項式係数に対する2次制約を導出し、一貫性があり単調なプロファイルを強制する、修正された端点勾配および端点値推定値を求める。
• 2段階の制限処理を実施する：まず、係数制約を用いて曲率点をセル境界に配置する。次に、隣接する値および勾配を修正することで、一貫性のない端点勾配を是正する。
• 滑らかな領域では自然PPM/PQM再構成を優遇し、不連続付近では単調制限付きバージョンを優先する非線形重み付け戦略を採用する。
• 移流優勢流体における高次精度を実現するため、任意ラグランジュ・オイラー（ALE）フレームワークに新しい勾配制限子を統合する。

実験結果

リサーチクエスチョン

• RQ1 WENO型非線形再構成戦略は、PPMおよびPQMにおいて滑らかな極値付近での精度を向上させるとともに、衝撃近傍での非振動的挙動を維持できるか？
• RQ2 提案された勾配制限子により、PPMおよびPQMが滑らかな問題において理論的な3次および5次収束率に達するか？
• RQ3 長時間の移流シミュレーションにおいて、従来の単調制限子と比較して、新しい制限子は数値的拡散および斜温混合の観点で優れているか？
• RQ4 点ごとの再構成に限らず、WENOに類似した非線形重み付けは、区分的多項式法へも効果的に適応可能か？
• RQ5 再構成多項式プロファイルの単調性と一貫性を保証するために、端点勾配および端点値にどのような修正が必要か？

主な発見

• 提案されたWENO型勾配制限子により、PPMに基づくスキームは十分に滑らかな問題において3次収束を達成でき、従来の単調制限子で見られる精度低下が解消された。
• PQMに基づくスキームも同様の条件下で5次収束を達成し、新しい制限子が完全な形式的精度を保持していることが示された。
• 滑らかな極値付近では、局所的極値の不自然な平坦化を回避することで、従来の単調制限子と比較して数値的拡散が顕著に低減された。
• 再構成プロファイルは単調かつ非振動的であり、曲率点がセル境界に制限されることで安定性と一貫性が保証された。
• 修正された端点勾配および端点値補正により、勾配推定値の不一致が解消され、最終的な補間関数の正確な単調性が保証された。
• 1次元および2次元のベンチマークテストにより、移流優勢流体における新しい制限子の向上した精度と頑健性が確認された。