[論文レビュー] A comparison of turbulence models in airship steady-state CFD simulations

Vitaly Voloshin, Yong K. Chen+1 more|arXiv (Cornell University)|Oct 10, 2012

Aerospace Engineering and Energy Systems参考文献 17被引用数 10

ひとこと要約

本研究では、さまざまな迎角における blimp の定常状態 CFD シミュレーションにおいて、標準 k-ω、SST k-ω、実現可能 k-ε、Spallart-Allmaras の4つの渦粘性係数ベースの乱流モデルを評価した。風洞データを用いた検証と1500反復のベンチマークを実施した結果、特に低迎角から中迎角において、精度と計算効率のバランスに優れるため、Spallart-Allmaras モデルが全体として最も最適であることが判明した。一方、高迎角では標準 k-ω モデルが精度が最も高く、リソース消費量は高いが、優れた性能を示した。

ABSTRACT

The accuracy and resource consumption of the four different turbulence models based on the eddy viscosity assumption, namely, $k-\varepsilon$, two $k-ω$ and Spallart-Allmaram models, in modeling airships are investigated. The test airship shape is a conventional shape. Three different angles of attack are considered. The results are checked against the wind tunnel experimental data. The resource consumption study is based on the benchmark of 1500 iterations. Based on all data obtained it is evident that Spallart-Allmaras model is the most optimal one in the majority of cases.

研究の動機と目的

• 定常状態 CFD シミュレーションにおける blimp の乱流モデルの精度と計算コストを評価すること。
• 風洞データとの検証に基づき、空飛ぶ飛行機の空力的シミュレーションに最も適した乱流モデルを特定すること。
• 異なる迎角における精度とリソース消費量のバランスを考慮した CFD プロフェッショナル向けのモデル選定ガイドラインを提供すること。
• 近いゼロ、約10°、約35°の3つの迎角におけるモデル性能を評価すること。
• 1500回の反復という固定回数を用いて、モデル間の公平な比較が可能なリソース消費量のベンチマークを設定すること。

提案手法

• 本研究では、商用 CFD ソフトウェア Star-CCM+ 6.04 を用いて、3つの迎角（近いゼロ、約10°、約35°）における一般的な blimp 形状のシミュレーションを実施した。
• 研究では、4つの渦粘性係数ベースの乱流モデル（標準 k-ω、SST k-ω、実現可能 k-ε、Spallart-Allmaras）を用いたレイノルズ平均ナビエ＝ストークス（RANS）方程式を採用した。
• モデルの精度は、特に圧力力が支配的（98％寄与）である法線力係数 $ c_z $ について、実験的風洞データと比較して検証した。
• 計算リソース消費量は、全モデルで1500反復を完了するまでの CPU 時間で測定した。
• Spallart-Allmaras モデルは、1つの変数 $ \tilde{\nu} $ の輸送方程式を用い、$ \nu_t = \tilde{\nu} f_{v1} $ を通じて乱流粘性係数を直接モデル化する。これには、経験的定数とブレンド関数のセットが含まれる。
• 境界層効果の正確な解像度を確保するため、特に粘性力の予測に重要である近壁領域にメッシュの細分化を施した構造的メッシュを用いた。

実験結果

リサーチクエスチョン

• RQ1 低迎角、中迎角、高迎角における blimp の法線力係数 $ c_z $ を予測する際、どの乱流モデルが最も高い精度を示すか？
• RQ2 固定反復数（1500回）における4つの乱流モデルの計算リソース消費量（CPU 時間）は、どのように比較されるか？
• RQ3 Spallart-Allmaras モデルは、全迎角範囲で高い精度と効率を維持しているか、それとも特定の条件下で他のモデルが優れるか？
• RQ4 圧力力が空飛ぶ飛行機の空力的挙動をどれほど支配しているか、各乱流モデルがそれらをどれほど正確に予測できるか？
• RQ5 標準 k-ω モデルのような複雑なモデルを Spallart-Allmaras モデルよりも採用するにあたって、精度と計算コストのトレードオフが正当化されるか？

主な発見

• Spallart-Allmaras モデルは、特に低迎角（近いゼロ）および中迎角（約10°）において、精度と計算効率のバランスが最も優れており、$ c_z $ 予測の誤差が5.8％にとどまった。
• 高迎角（約35°）では、標準 k-ω モデルが最も高い精度を示し、誤差は5.8％であった。Spallart-Allmaras モデルはこの状況で7.7％の誤差を示し、劣っていた。
• 実現可能 k-ε モデルは1500反復に要する CPU 時間が最長で、精度は中程度（誤差8.5％）であり、計算効率が最も低かった。
• 高迎角においても、Spallart-Allmaras モデルは標準 k-ω モデルに比べて約13％の CPU 時間を削減しており、計算効率に顕著な優位性を示した。
• 圧力力は法線力係数 $ c_z $ の98％を占めており、正確な圧力予測がモデル性能の主な決定要因であることが判明した。
• 予備的な空飛ぶ飛行機シミュレーションでは、k-ω モデルによる2％の精度向上と、Spallart-Allmaras モデルによる13％のリソース消費量削減というトレードオフを考えると、コスト上昇を理由に前者を優先する理由は薄く、多くの用途では後者を推奨する。