[論文レビュー] High resolution methods for time dependent problems with piecewise smooth solutions

Eitan Tadmor|ArXiv.org|May 1, 2003

Computational Fluid Dynamics and Aerodynamics参考文献 20被引用数 7

ひとこと要約

本稿では、区分的連続な解を有する時間依存の対流優勢問題に対する高解像度数値法を提示する。局所的手法には中央差分法を、大域的手法には適応的スペクトル粘性（SV）を用いる。エッジ検出、適応的再構築、スペクトル粘性を用いた指数的収束を示し、振動を伴わず衝撃波や不連続性を高精度に解明する。

ABSTRACT

A trademark of nonlinear, time-dependent, convection-dominated problems is the spontaneous formation of non-smooth macro-scale features, like shock discontinuities and non-differentiable kinks, which pose a challenge for high-resolution computations. We overview recent developments of modern computational methods for the approximate solution of such problems. In these computations, one seeks piecewise smooth solutions which are realized by finite dimensional projections. Computational methods in this context can be classified into two main categories, of local and global methods. Local methods are expressed in terms of point-values (-- Hamilton-Jacobi equations), cell averages (-- nonlinear conservation laws), or higher localized moments. Global methods are expressed in terms of global basis functions. High resolution central schemes will be discussed as a prototype example for local methods. The family of central schemes offers high-resolution ``black-box-solvers'' to an impressive range of such nonlinear problems. The main ingredients here are detection of spurious extreme values, non-oscillatory reconstruction in the directions of smoothness, numerical dissipation and quadrature rules. Adaptive spectral viscosity will be discussed as an example for high-resolution global methods. The main ingredients here are detection of edges in spectral data, separation of scales, adaptive reconstruction, and spectral viscosity.

研究の動機と目的

• 非滑らかな特徴（衝撃波や折り返し点）を有する対流優勢問題を計算する課題に対処すること。
• 不連続性の存在下でも精度と安定性を保持する高解像度数値法を開発すること。
• 区分的連続性を活用することで、理論的収束率と実際の計算性能のギャップを埋めること。
• スペクトル精度を維持しながらエントロピー散逸を強制する安定化手法としての適応的スペクトル粘性を導入すること。
• 適応的モリファイアと物理空間およびフーリエ空間におけるフィルタリングを用いた再構築における指数的収束を示すこと。

提案手法

• 点値、セル平均、非振動的再構築に基づく高解像度中央差分法を用い、振動を伴わず衝撃波を捉える。
• フィルタ $\sigma_p(\xi) = e^{\beta \xi^p / (\xi^2 - 1)}$ を用いた適応的スペクトル粘性（SV）を導入し、スペクトル法の安定化とギブス効果の防止を図る。
• 距離関数 $\theta = d(x)$ および $p \sim d(x)N$ を用いた適応的モリファイア $\Psi_{\theta,p}$ を用い、滑らかな領域で指数的収束を達成する。
• ディリクレ核 $D_p(y)$ とゲヴレークラスのモリファイア $\rho_\beta$ を用い、局所化とキャンセレーション効果を確保し、高精度を実現する。
• エッジ検出とスペクトルフィルタリングを組み合わせ、滑らかさの方向に解を再構築し、不連続性付近の誤差を最小化する。
• SV正則化を施した修正ガラーキン法 $\partial_t v_N + \mathcal{P}_N \nabla_x \cdot f(v_N) = -N \sum_{|k|\leq N} \sigma(|k|/N) \hat{v}_k \phi_k(x)$ を導入し、エントロピー散逸を保証する。

実験結果

リサーチクエスチョン

• RQ1 対流優勢問題における区分的連続な解を高精度に解明するための高解像度数値法は、どのように設計可能か？
• RQ2 適応的スペクトル粘性は、スペクトル精度を損なわずにスペクトル法を安定化させる役割を果たすか？
• RQ3 区分的連続関数のフーリエ変換またはスペクトル射影からの再構築において、指数的収束を達成できるか？
• RQ4 モリファイアパラメータ $\theta$ と $p$ の適応的選択は、不連続性付近での分解能をどのように向上させるか？
• RQ5 不連続性からの距離 $d(x)$、モード数 $N$、およびスペクトル法における収束率との関係は何か？

主な発見

• 区分的連続関数の適応的再構築において指数的収束が達成される：$|\Psi_{\{d(x),d(x)N\}}*\mathcal{P}_N v(x) - v(x)| \leq C (d(x)N)^r e^{-\gamma \sqrt{d(x)N}}$。
• $L^1$-誤差は高解像度スキームで $\|v^h - u\|_{L^1} \leq \alpha_t h |\log h|$ の速度で収束し、一般のBV解の理論的 $O(h^{1/2})$ よりも著しく速い。
• 適応的スペクトル粘性により、エントロピー散逸を強制することでスペクトル近似の強収束が保証され、不適切な振動が防止される。
• わずかに局所的に高波数成分に粘性を追加するだけで、形式的スペクトル精度を維持し、$s \in (0, \infty)$ が散逸スケールを制御する。
• 数値実験により、物理空間およびフーリエ空間の両方で指数的収束率が確認され、図5～8で示される。
• 理論的解析により、非線形保存則では $p > 1$ に対して $L^p$-安定性が成立しないことが示され、BVおよび区分的連続正則性フレームワークの必要性が裏付けられる。