[論文レビュー] Variational structure of the optimal artificial diffusion method for the advection-diffusion equation

K. B. Nakshatrala, Albert J. Valocchi|arXiv (Cornell University)|May 29, 2009

Advanced Numerical Methods in Computational Mathematics参考文献 8被引用数 3

ひとこと要約

本稿は、変分原理から安定な弱形式を導出し、対流拡散問題における最適な人工拡散法の変分的基礎を確立する。この方法は、定数係数および定数項を伴う1次元問題において正確なノード解を生成し、最適な人工拡散法が指数関数的重み関数を用いた重み付き内積に基づく変分形式と等価であることを証明する。

ABSTRACT

In this research note we provide a variational basis for the optimal artificial diffusion method, which has been a cornerstone in developing many stabilized methods. The optimal artificial diffusion method produces exact nodal solutions when applied to one-dimensional problems with constant coefficients and forcing function. We first present a variational principle for a multi-dimensional advective-diffusive system, and then derive a new stable weak formulation. When applied to one-dimensional problems with constant coefficients and forcing function, this resulting weak formulation will be equivalent to the optimal artificial diffusion method. We present representative numerical results to corroborate our theoretical findings.

研究の動機と目的

• 最適な人工拡散法に変分的枠組みを提供すること。これは、これまでにそのような変分的定式化が存在しなかった。
• 対流優勢問題における古典的ガレルキン法の不安定性および振動を是正すること。
• 変分原理を用いて、対流拡散方程式の安定な弱形式を導出すること。
• 導出された弱形式が、1次元の場合に最適な人工拡散法と同一の差分方程式を生成することを示すこと。
• 非自己随伴作用素（例：対流拡散作用素）に対して変分原理が果たす理論的意義を強調すること。

提案手法

• 対流拡散方程式を非自己随伴作用素として定式化し、重み関数 α(x) = exp(−vx/k) を用いた重み付き内積を用いて変分原理を導出する。
• ヴァインバーグの定理を適用し、重み付き残差表現をスカラー汎関数に結びつけることで、変分的定式化を可能にする。
• 重み付き内積を含む汎関数を最小化することで、対流優勢問題に対しても安定な弱形式を導出する。
• 3ノードステンシル上で定義された線形有限要素と形状関数を用いて、離散差分方程式を導出する。
• ペクレ数に依存する重み関数 α(x) を導入し、1次元問題における安定性およびノード正確性を保証する。
• 得られた離散系を古典的な最適な人工拡散法と比較し、1次元において等価であることを確認する。

実験結果

リサーチクエスチョン

• RQ1 最適な人工拡散法は変分原理から導出可能か？
• RQ2 定数係数を伴う1次元問題において、ノード正確解を生成する安定な弱形式が対流拡散方程式に存在するか？
• RQ3 ヴァインバーグの定理を用いて、対流拡散作用素に対してスカラー汎関数を構成可能か？
• RQ4 導出された変分的定式化は、1次元において最適な人工拡散法と同一の離散系を生成するか？
• RQ5 非一様重み関数を用いた変分的定式化が、非自己随伴問題に対して果たす理論的意義は何か？

主な発見

• 変分原理から導出された安定な弱形式は、定数係数および定数項を伴う1次元対流拡散問題において、最適な人工拡散法と同一の差分方程式を生成する。
• 得られた離散系の係数は、最適な人工拡散法のものと完全に一致する：γ−1 = −v/(2h)(1 + coth(Ph/e))、γ0 = v/h coth(Ph/e)、γ+1 = v/(2h)(1 − coth(Ph/e))。
• 数値結果により、本新規定式化が、1次元においてすべてのペクレ数に対してノード正確解を達成することが確認され、理論的等価性が裏付けられる。
• 最適な人工拡散法は、導出された汎関数および重み関数 α(x) = exp(−vx/k) を通じて、明確な変分的基礎を有することが示された。
• 本手法は理論的に堅牢であり、新たな安定化定式化の開発に基盤を提供するが、大規模な2次元および3次元問題では計算コストが高くなる可能性がある。
• 本研究は、非自己随伴系に対する変分原理の理論的枠組み（トンティの理論的業績）が、安定化有限要素法の文脈においていかに重要であるかを強調している。