[論文レビュー] A general approach to enhance slope limiters on non-uniform grids

Xianyi Zeng|arXiv (Cornell University)|Jan 6, 2013

Computational Fluid Dynamics and Aerodynamics参考文献 19被引用数 3

ひとこと要約

本稿では、特にMUSCL法を対象として、非一様メッシュにおける勾配制限関数の性能を向上させる一般化されたフレームワークを提示する。第二階層の精度、TVD安定性、対称性保存のための十分条件を導出することで、古典的な再構成およびハルテンの安定性理論を非一様1次元メッシュへと拡張する。1次元および2次元の数値例を通じて検証され、従来の制限関数では不規則メッシュで失われていた第二階層収束性が回復される。

ABSTRACT

Abstract. A general approach to study and enhance the slope limiter functions on non-uniform grids is presented. Slope limiters are preferred in high-resolutions schemes in general and MUSCL in particular to solve hyperbolic conservation laws. However, most 1D limiters are developed assuming uniform meshes in space, which are shown to be inadequate on non-uniform grids. Especially, second-order convergence is shown to be lost when the conventional limiters are applied on irregular grids in the case of smooth solutions. A methodology based on the classical reconstruct-evolve-project approach and Harten’s stability theory is presented to study the slope limiters on 1D non-uniform computational grids. Sufficient conditions for the limiters to lead to formal second-order spatial accuracy, total-variational-diminishing stability and symmetry-preserving property are derived. The analysis and results extend naturally to cell-centered finite volume methods in multiple dimensions. Several most widely used conventional limiters are enhanced to satisfy these conditions, and their performances are illustrated by various 1D and 2D numerical examples.

研究の動機と目的

• 高解像度スキームにおける非一様メッシュ上に従来の勾配制限関数を適用した場合に生じる第二階層空間精度の損失を是正すること。
• 再構成-発展-投影フレームワークおよびハルテンの安定性理論に基づく、非一様1次元メッシュ上での勾配制限関数の分析のための体系的アプローチを開発すること。
• 非一様離散化において形式的第二階層精度、全変動減少（TVD）安定性、対称性保存特性を保証する十分条件を導出すること。
• 1次元メッシュから2次元セル中心有限体積法への分析を拡張すること。
• 導出された条件を満たすように、広く用いられる従来の制限関数を改良し、その性能を数値的に検証すること。

提案手法

• 非一様メッシュ上での勾配制限関数の分析の基盤として、再構成-発展-投影アプローチを採用する。
• ハルテンの安定性理論を応用し、非一様メッシュの文脈においてTVD安定性および対称性保存のための十分条件を導出する。
• 非一様メッシュ上での形式的第二階層空間精度を保証するため、制限関数に数学的制約を定式化する。
• 従来の制限関数（例：minmod、van Leer、Superbee）を、導出された条件を満たすように制限関数を変更することで変換する。
• 1次元および2次元の有限体積スキームに改良された制限関数を実装し、不規則メッシュ上での精度と安定性を検証する。
• 滑らかな解を用いて収束率を定量的に評価し、非一様メッシュ上での第二階層精度の回復を確認する。

実験結果

リサーチクエスチョン

• RQ1 非一様メッシュ上での1次元双曲型保存則において、勾配制限関数が第二階層空間精度を保持するためには、どのような条件を満たす必要があるか？
• RQ2 ハルテンの安定性理論は、非一様計算メッシュ上での勾配制限関数の分析にどのように適合可能か？
• RQ3 非一様メッシュ上でのTVD安定性および対称性保存を確保するために、従来の制限関数にどのような修正が必要か？
• RQ4 改良された制限関数は、不規則メッシュ上での滑らかな解において、第二階層収束性をどの程度回復できるか？
• RQ5 提案されたフレームワークは、多次元セル中心有限体積法へと拡張可能か？

主な発見

• 従来の勾配制限関数は、非一様メッシュ上に適用された場合でさえ、滑らかな解に対しても第二階層空間精度を失う。
• 導出された十分条件により、非一様1次元メッシュ上での勾配制限関数に対して、形式的第二階層精度、TVD安定性、対称性保存特性が保証される。
• 広く用いられる制限関数（例：minmod、van Leer、Superbee）の改良版は、非一様メッシュ上でも第二階層収束性を効果的に回復する。
• この手法は、多次元セル中心有限体積法へ自然に一般化可能である。
• 1次元および2次元の数値例により、改良された制限関数が不規則メッシュ上でも高解像度特性を維持しながら最適な精度を達成することが確認された。
• 本手法は、実用的な計算流体力学応用分野における非一様メッシュ向けに、強固で高精度な制限関数を設計する体系的フレームワークを提供する。