[論文レビュー] Eliminating spurious eigenvalues in the analysis of incompressible fluids and other systems of differential-algebraic equations

M. Lisa Manning, Bassam Bamieh+1 more|arXiv (Cornell University)|May 10, 2007

Numerical methods for differential equations参考文献 30被引用数 2

ひとこと要約

本稿では、非圧縮性流体やその他の微分代数方程式（DAE）の安定性解析における余剰固有値を排除する記述子表記フレームワークを提案する。解析的削減を伴わずに微分および代数的制約を同時に解くことで、従来の数値的手法で一般的に生じる物理的に意味のない不安定性を回避し、流体力学的安定性問題における正確な固有値計算を保証する。

ABSTRACT

We describe a general framework for avoiding spurious eigenvalues -- unphysical unstable eigenvalues that often occur in hydrodynamic stability problems. In two example problems, we show that when system stability is analyzed numerically using {\em descriptor} notation, spurious eigenvalues are eliminated. Descriptor notation is a generalized eigenvalue formulation for differential-algebraic equations that explicitly retains algebraic constraints. We propose that spurious eigenvalues are likely to occur when algebraic constraints are used to analytically reduce the number of independent variables in a differential-algebraic system of equations before the system is approximated numerically. In contrast, the simple and easily generalizable descriptor framework simultaneously solves the differential equations and algebraic constraints and is well-suited to stability analysis in these systems.

研究の動機と目的

• 非圧縮性流体の流れの数値的安定性解析における、長年の余剰固有値の問題に対処すること。
• DAEにおける変数の解析的削減が、物理的に意味のない固有値の出現を引き起こす原因であることを特定すること。
• 数値的解法中に代数的制約を保持する記述子表記フレームワークを提案し、固有値計算の精度を向上させること。
• 記述子アプローチが流体力学的安定性問題における物理的に意味のない不安定性を排除する有効性を実証すること。
• 広範な微分代数系に適用可能な一般的で容易に拡張可能な手法を提供すること。

提案手法

• 本稿では、DAEにおける代数的制約を明示的に保持する一般化固有値問題としての記述子表記を採用する。
• 独立変数の数を解析的に削減することを避け、数値近似中にシステムの完全な構造を保持する。
• 統一された行列フレームワーク内で、微分方程式と代数的制約を同時に解く。
• 安定性は、記述子表記から導かれる一般化固有値問題を介して分析する。
• 本手法の有効性を検証するために、2つの流体力学的安定性問題に適用する。
• 本フレームワークは、流体力学を越えた他のDAEシステムへも容易に拡張可能であるように設計されている。

実験結果

リサーチクエスチョン

• RQ1 なぜ非圧縮性流体系の数値的安定性解析において余剰固有値が出現するのか？
• RQ2 DAEにおける変数の解析的削減が、物理的に意味のない固有値の出現にどのように寄与するのか？
• RQ3 記述子に基づく定式化によって、余剰固有値を排除しつつ安定性解析の正確さを維持できるか？
• RQ4 記述子フレームワークは、流体力学を越えた他の微分代数系にも一般化可能か？
• RQ5 数値近似中に代数的制約を保持することが、固有値計算にどのような影響を与えるか？

主な発見

• 代数的制約が数値近似の前に解析的に削除されると、物理的に意味のない不安定性が生じ、余剰固有値が出現する。
• 記述子フレームワークは、数値的解法中にシステムの完全な構造を保持することで、余剰固有値を効果的に排除する。
• 変数の削減を伴わずに微分方程式と代数的制約を同時に解くことで、安定性解析の精度を維持する。
• 2つのテストケース（非圧縮性流体の流れを含む）において、本手法は頑健性と信頼性を示した。
• 記述子定式化は、DAEシステムにおける安定性解析のための一般的で安定的かつ容易に一般化可能なフレームワークを提供する。
• 制約を明示的かつ一貫的に扱えるため、本フレームワークは数値的安定性解析に適している。