[論文レビュー] A numerical stability investigation of strong ZND detonations for Majda's model

Blake Barker, Kevin Zumbrun|arXiv (Cornell University)|Nov 6, 2010

Combustion and Detonation Processes参考文献 26被引用数 7

ひとこと要約

本稿は、活性化エネルギーが変化する範囲において、アーレニウス型 igniting 関数を有するMajdaモデルにおける強力なZND炎上波の数値的安定性を調査する。高度な固有値ソルバーを用いて、高活性化エネルギー極限においても振動不安定性は観測されず、Humpherys–Zumbrun法が古典的Lee–Stewart法に比べ100～1,000倍速いことが示された。特にアダプティブメッシュ化により性能向上が顕著であった。

ABSTRACT

We carry out a systematic numerical stability analysis of ZND detonations of Majda's model with Arrhenius-type ignition function, a simplified model for reacting flow, as heat release and activation energy are varied. Our purpose is, first, to answer a question of Majda whether oscillatory instabilities can occur for high activation energies as in the full reacting Euler equations, and, second, to test the efficiency of various versions of a numerical eigenvalue-finding scheme suggested by Humpherys and Zumbrun against the standard method of Lee and Stewart. Our results suggest that instabilities do not occur for Majda's model with Arrhenius-type ignition function, nor with a modified Arrhenius-type ignition function suggested by Lyng--Zumbrun, even in the high-activation energy limit. We find that the algorithm of Humpherys--Zumbrun is in the context of Majda's model $100$-$1,000$ times faster than the one described in the classical work of Lee and Stewart and $1$-$10$ times faster than an optimized version of the Lee--Stewart algorithm using an adaptive-mesh ODE solver

研究の動機と目的

• 高活性化エネルギー極限において、アーレニウス型 igniting 関数を有するMajdaモデルに、ホフプ・分岐を伴う振動不安定性が生じるかどうかを特定すること。
• 炎上波安定性解析の文脈において、Humpherys–Zumbrun固有値探索スキームが古典的Lee–Stewart法に比べてどの程度効率的であるかを検証すること。
• メッシュの適応的生成、座標パrametrization、スキーム構成といった数値的実装選択が、計算性能および収束性に与える影響を評価すること。
• 不安定固有値の厳密な境界を確立し、特にスイッチバック近似を含む特異極限を扱うハイブリッドスキームを開発すること。
• Majdaモデルを簡略化されたテストベッドとして用い、より複雑な反応性Euler系に適用可能な数値的手法のベンチマークを提供すること。

提案手法

• アーレニウス型および修正アーレニウス型 igniting 関数を有するMajdaモデルにおけるZND炎上波プロファイルに対して、体系的な数値的固有値解析を実施した。
• Humpherys–Zumbrunのエヴァンズ関数型スキームを、x=0からx=−∞への後退シューティング法に再定式化し、スペクトルギャップの利点を活用した。
• 空間（x）および周波数（λ）変数の両方において、アダプティブステップのODEソルバーを適用し、収束性の向上と計算コストの低減を図った。
• z座標系（プロファイルがzにおいて明示的に解ける）とx座標系（元の定式化）を用い、その性能への影響を比較した。
• 高周波数漸近展開とハイブリッドスキームを発展させ、特にスイッチバック近似を含む特異極限を扱った。
• 多数のスキーム間で広範な数値的比較を実施した：μ(x)、随伴μ(x)、ハイブリッド、Lee–Stewart（アダプティブおよび固定メッシュ）、極座標随伴、極座標随伴径方向法。

実験結果

リサーチクエスチョン

• RQ1 Majdaモデルにアーレニウス型 igniting 関数を適用した場合、全反応性Euler方程式と同様に高活性化エネルギー極限において振動不安定性が生じるか？
• RQ2 Humpherys–Zumbrun固有値探索スキームは、ZND炎上波安定性解析において、古典的Lee–Stewart法に比べて計算効率が優れているか？
• RQ3 アダプティブメッシュ対固定メッシュ、x座標系対z座標系、スキーム構成といった異なる数値的選択が、計算速度および収束性に与える影響は何か？
• RQ4 Majdaモデルにおける不安定固有値の最大サイズは何か？ その厳密な境界を導出可能か？
• RQ5 ハイブリッドまたは極限スキームは、igniting 関数の特異的・スイッチバック近似を安定性計算で効果的に扱えるか？

主な発見

• 高活性化エネルギー極限においても、アーレニウス型または修正アーレニウス型 igniting 関数を有するMajdaモデルにおいて、振動的または他の種類の不安定性は観測されなかった。
• Humpherys–Zumbrunスキームは、標準的Lee–Stewart法に比べ100～1,000倍速く、最適化されたLee–Stewartバージョンに対しても1～10倍の性能向上を示した。
• x方向の進化におけるアダプティブステップODEソルバーは、固定メッシュソルバーに比べ計算時間を100～1,000倍短縮し、最も顕著な性能向上を示した。
• z座標系法は、x座標系法に比べ一様に6倍の性能向上を示した。主にプロファイル評価コストの低減によるものであった。
• 最適な数値的選択を採用した場合、1パラメータセットあたり数秒で安定性計算が完了し、簡略化モデルにおける良好な数値的条件付けが示された。
• 高活性化エネルギー極限においては性能が低下したが、成長/減衰挙動の因子分解により、劣化したスキームが失敗する状況でも収束が可能となった。