[論文レビュー] Highly Accurate Solutions of the Blasius and Falkner-Skan Boundary Layer Equations via Convergence Acceleration

B. D. Ganapol|arXiv (Cornell University)|Jun 19, 2010

Iterative Methods for Nonlinear Equations参考文献 15被引用数 14

ひとこと要約

この論文では、再帰的係数を有するマクローリン級数展開、ウィン・エプシロン法による収束加速、連続的解析接続を用いた、ファルクナー＝スカンおよびブレイシウス境界層方程式を解く高精度で安定した数値アルゴリズムを提示する。この手法は、拡張精度演算や記号計算を用いずに、倍精度フォートランで10桁以上の精度を達成し、摩擦係数のベンチマークを新たに設定した。

ABSTRACT

A new highly accurate algorithm for the solution of the Falkner-Skan equation of boundary layer theory is presented. The algorithm, based on a Maclaurin series representation, finds its coefficients from recurrence. In addition, Wynn-epsilon convergence acceleration and continuous analytical continuation enable an accurate evaluation. The most accurate skin friction coefficients (shooting angle) to date are presented along with comparisons to past and present values found in the literature. The algorithm, coded in FORTRAN, uses neither enhanced precision arithmetic beyond quadruple precision nor computer algebra to achieve results in a timely fashion. Key Words: Falkner-Skan flow; Blasius flow; Wynn-epsilon acceleration; Romberg acceleration; Continuous analytical continuation

研究の動機と目的

• ファルクナー＝スカンおよびブレイシウス境界層方程式を解く、安定で高精度な数値アルゴリズムの開発。
• 初期値の微調整が必要で、数値的不安定性を示す既存の数値手法の限界を克服すること。
• 四重精度演算やコンピュータ代数システムに依存せずに、高精度（10桁以上）を達成すること。
• 文献値との比較によるスキン摩擦係数のベンチマーク化と、先行研究を上回る精度の延長。
• 前向きおよび逆向きの流れ領域、特に負のβ値を含む、両方の境界層流れ状態における手法の有効性の検証。

提案手法

• ファルクナー＝スカンODEから導出された再帰関係により生成される係数を有するマクローリン級数として解を表現する。
• ファルクナー＝スカン方程式の収束半径を超えて級数の収束を向上させるために、ウィン・エプシロン法による収束加速を適用する。
• 速度プロファイルが1に近づく漸近領域に解を拡張するために、連続的解析接続を用いる。
• 境界条件が無限遠で満たされるように、初期勾配（シューティング角）を繰り返し調整するためのバイセクション法に基づくシューティング法を用いる。
• ニュートン・ラプソン反復を避け、物理的直観とバイセクション法に依存することで、安定性と収束性を確保する。
• すべての計算は標準的な倍精度フォートランで実行され、精度の確認のため一部の四重精度チェックが実施される。

実験結果

リサーチクエスチョン

• RQ1 収束加速と解析接続を併用したマクローリン級数に基づく手法が、拡張精度演算を用いずにファルクナー＝スカン方程式に対して高精度な解を得られるか。
• RQ2 初期値の微調整を要するルンゲ＝クッタ法に依存する既存の数値手法と比較して、提案手法の精度と安定性はどのように異なるか。
• RQ3 前向き流れを示す負のβ値を含む、ブレイシウスおよびファルクナー＝スカン方程式の最も正確なスキン摩擦係数（シューティング角）は何か。
• RQ4 本手法は、β ≈ -0.198837735という困難なケースを含め、前向きおよび逆向きの境界層流れを信頼性高く取り扱えるか。
• RQ5 収束加速と解析接続を用いることで、マクローリン級数の自然な収束半径を超えて、解の有効な評価領域をどの程度まで拡張できるか。

主な発見

• 本アルゴリズムは、倍精度演算のみを用いても、ファルクナー＝スカン方程式の解を10桁以上（有効数字）の精度で達成した。
• ブレイシウスの場合（β = 0）の最も正確なスキン摩擦係数はα = -0.469599988361と報告され、新たなベンチマーク値として-0.46959998836100000000が提示された。
• β = -0.198837735の場合、スキン摩擦係数はα = -0.19883773500000000000として計算され、流れの逆転を示す臨界値の存在が確認された。
• β < 0の場合、特にβ = -0.16およびβ = -0.198837735に対し、一貫性があり正確なプロファイルを有する逆流れ解を効果的に計算できた。
• 本手法は、極めて悪い初期推定値（例：β = 1000に対してα = 200）からでも収束し、α = 36.5171968という高信頼性の結果を得られた。
• 本研究は、ブレイシウスの元々のマクローリン級数法が、現代の収束技術と組み合わせることで、境界層問題に対して依然として強力で正確な手法であることを確認した。