空力学 — CAE用語解説
空力学とは
空力学ってざっくり言うと何ですか? Navier-Stokes方程式を解くのと何が違うんですか?
理論と物理
基本概念と支配方程式
空力解析でよく出てくる「抗力係数」と「揚力係数」って、具体的に何を表しているんですか?単に力の大きさじゃないんですよね?
その通り、力そのものではなく「無次元化された力の大きさ」です。具体的には、機体の形状と迎角が、どれだけ効率的に抗力や揚力を生み出しているかを示す指標です。例えば、同じ速度・空気密度・翼面積でも、抗力係数
支配方程式はナビエ-ストークス方程式と聞きますが、実際のCAEでは全部解いているんですか?それとも簡略化するんですか?
用途によります。最も一般的なのは、時間平均したRANS(Reynolds-Averaged Navier-Stokes)方程式です。これにより、乱流の詳細な時間変動はモデル化し、平均的な流れ場を効率的に解きます。式で書くと、連続の式と運動量式は以下のようになります。
「圧力係数」
物体表面の各点での圧力を、無次元化した指標です。定義は
数値解法と実装
離散化とソルバー設定
CFDで「セル」や「要素」の種類(四面体、六面体など)を選ぶ時、空力解析では何を基準に決めるんですか?
境界層の分解能が最大のポイントです。物体表面に垂直な方向には、流速の勾配が非常に大きい「境界層」が形成されます。これを正確に捉えるには、表面に沿って細長いプリズム(六面体の一種)や多面体セルを積層するのが一般的です。具体的には、第1層セルの高さを無次元壁面距離
ソルバーの「圧力ベース」と「密度ベース」の違いは何ですか?どちらを選べばいいですか?
歴史的には、非圧縮性流れ(マッハ数Ma < 0.3)には圧力ベース、圧縮性流れ(Ma > 0.3)には密度ベースが適していました。しかし現在のソフトウェア(Fluent, STAR-CCM+など)では、その境界は曖昧です。実務的な選択基準はこうです:圧力ベースは、自動車や建築物の低速流れ、自然対流など、広範囲な問題にロバストです。一方、密度ベースは、超音速流れや高度に圧縮性が関与する燃焼解析で、衝撃波の解像度が高い傾向があります。ただし、Ansys Fluentでは密度ベースソルバーもすべての速度領域で使えるように改良されています。迷ったら、まずは圧力ベースの「Coupled」アルゴリズムから始めるのが無難です。
「定常解析」と「非定常解析」は、どういう時に切り替えるんですか?非定常は計算時間がかかるので避けたいですが。
流れ場に「本質的に時間変動する構造」があるかどうかです。例えば、自動車のリアビークル後方の渦(ウィークトラフィック)は、定常解析でも平均的なパターンは得られます。しかし、橋梁の「ギャロッピング」や航空機の「フラッタ」、風力発電ブレードの「ストール現象」は、時間とともに流れ構造が大きく変動するため、非定常解析が必須です。また、カルマン渦列(円柱後流など)も非定常です。判断の一助として、定常解析を実行した後、残差が振動し続けたり、力の係数が一定値に収束しない場合は、非定常現象の可能性が高いです。非定常解析では、時間刻み幅は現象の特徴周波数(例えばストローハル数から推定)の20倍以上になるように設定します。
実践ガイド
ワークフローとチェックリスト
空力解析を始める前の「計算領域(ドメイン)の大きさ」は、どう決めればいいですか?大きすぎると計算コストが、小さすぎると結果がおかしくなりそうで。
良い質問です。これは規格やベストプラクティスが参考になります。自動車業界では、SAE J1252やJ2966といった規格にガイドラインがあります。一般的な目安は、物体の特性長さL(車なら全長)に対して、
・流入面:物体前端から 3L ~ 5L 上流
・流出面:物体後端から 7L ~ 10L 下流
・側面/上面/下面:物体中心から 3L ~ 5L
です。例えば全長4.5mのセダンなら、計算領域は少なくとも幅・高さ約20m、長さ約50mは必要です。建築物の風荷重解析では、AIJ(日本建築学会)の「建築物荷重指針」が参考になります。領域が小さすぎると、境界の影響で流れが歪み、特に後流の圧力分布が不正確になります。
境界条件で「速度入口」と「圧力入口」の使い分けは?「壁面」の設定で注意することは?
「速度入口」は流入速度が既知の場合(風洞実験を模倣する場合など)に使います。一方、「圧力入口」は流入速度が結果として決まる場合や、外部の大気開放条件を模擬する場合に使います。実務では、自動車の走行解析では「速度入口」が、建築物周りの自然風解析では乱流強度と特性長さを指定した「圧力入口」がよく使われます。
「壁面」で最も重要なのは、先述の境界層メッシュと、壁面関数の選択の整合性です。
結果の収束をどう判断すればいいですか?残差が下がればOK?
残差の低下は必要条件ですが、十分条件ではありません。最も重要なのは、求めたい物理量(抗力係数
1. 主要な残差(連続、運動量)が少なくとも3〜4桁低下し、低い値で平坦化している。
2. モニタしている力や係数の値が、最後の500〜1000ステップで変動幅が1%未満になっている。
3. 計算領域の出入口での質量流量や運動量流量のバランスが取れている(誤差1%以内)。
4. 物理的に意味のある値になっているか(例:自動車の
収束しない場合、メッシュ品質や時間刻み(非定常の場合)、ソルバーの緩和係数を疑います。
ソフトウェア比較
Ansys Fluent vs STAR-CCM+ vs OpenFOAM
自動車メーカーではSTAR-CCM+がよく使われていると聞きます。Fluentと比べて、空力解析において決定的に優れている点は何ですか?
最大の特徴は「オールインワン」の統合環境と、強力な自動メッシング機能です。STAR-CCM+は前処理(CAD処理、メッシュ生成)、求解、後処理が単一のGUI内で完結し、ワークフローの自動化(Javaマクロ)が非常に強力です。自動車の外気解析では、「Surface Remesher」でCADの細かい傷を吸収し、「Polyhedral Mesher」と「Prism Layer Mesher」を組み合わせて、数クリックで複雑な車体周りの高品質なメッシュを生成できます。また、車両運動(ピッチング、ヨーイング)との連成解析や、バッテリー熱流体とのマルチフィジックス連成も得意としています。Fluentも同様の機能はありますが、ワークフローの統合性と自動化の容易さでSTAR-CCM+が好まれる傾向があります。
では、Ansys Fluentはどんな場面で使うメリットがあるんですか?
Fluentは、豊富な物理モデルと長い歴史に裏打ちされたロバスト性、そしてAnsysエコシステムとの親和性に強みがあります。特に、
・複雑な化学反応を伴う燃焼解析(航空エンジンなど)
・多相流(スプレー、気泡)
・高度な乱流モデル(DES, LESの各種派生モデル)
・Ansys Mechanicalとの緊密な連成(流体構造連成:FSI)
といった高度な解析において、多くの実績とノウハウが蓄積されています。また、ユーザー定義関数(UDF)によるカスタマイズ性も高いです。航空宇宙分野や、研究開発機関では依然としてFluentのシェアが高いです。
無償のOpenFOAMは、有償ソフトと比べて実務で使えるレベルなんですか?
計算エンジンとしてのコア性能は非常に高く、研究や特定の産業分野(特に風力発電、船舶)で広く使われています。有償ソフトと同等かそれ以上の物理モデルを実装しています。しかし、最大のハードルは前処理と後処理です。GUIが貧弱(または標準では付属しない)で、メッシュ生成や設定はコマンドラインや辞書ファイルで行う必要があり、習得に時間がかかります。そのため、実務で導入するには、社内に専門知識を持つエンジニアがいるか、有償のサポートとGUIを提供するディストリビューション(ESI-OpenCFD, foam-extendなど)を購入する必要があります。予算は限られるが深くカスタマイズしたい研究開発プロジェクトには最適です。
トラブルシューティング
よくあるエラーと対策
計算中に「負の体積」や「負のヤコビアン」というエラーが出て止まります。これは何が原因で、どう直せばいいですか?
メッシュ品質が極端に悪いことが原因です。具体的には、
1. スキュー角が大きすぎる要素:特に六面体メッシュで、要素がひしゃげている。
2. アスペクト比が大きすぎる要素:細長すぎる要素。境界層メッシュは別ですが、流域内部では通常1000を超えると危険。
3. 急激なメッシュサイズ変化:隣接するセルの体積比が急激に変わると、勾配計算が破綻します。
対策は、メッシャー(Fluent Meshing, STAR-CCM+のメッシャーなど)の品質レポートを確認し、「スキュー角」を0.9以下(理想は0.8以下)、「アスペクト比」を流域内部で100以下に収めるようにメッシュ設定を見直します。また、リメッシュやメッシュスムージング機能を使うことも有効です。
定常解析なのに、抗力係数が一定値に収束せず、周期的に振動します。これは非定常現象ですか?
その可能性が高いです。例えば、鈍頭物体(トラックのキャブ後部など)からの渦放出は本質的に非定常です。定常ソルバーは「時間平均化された解」を見つけようとしますが、物理が非定常であるため、解が振動して収束しません。まず、振動の周波数を推定し(FFTなど)、その現象が設計上重要か判断します。重要であれば、非定常解析に切り替えます。重要でない平均値だけ欲しい場合は、定常解析の緩和係数をさらに強くする(デフォルトの0.2から0.1などに下げる)か、RANSモデルをSST k-ωからより拡散性の高いk-ε Realizableモデルに変更することで、擬似的に「平均化」された解が得られることがありますが、根本的な解決にはなりません。
風洞実験のデータとCFDの結果が、特に高迎角で大きくずれます。原因として考えられることは?
高迎角では流れの剥離が大きく、その再現性が鍵です。主な原因は:
1. 乱流モデルの限界:標準的なRANSモデルは剥離点や大規模剥離域の予測が苦手です。DES(Detached Eddy Simulation)やLESへの切り替えを検討します。
2. メッシュ分解能不足:剥離・再付近領域のメッシュが粗すぎます。剥離が予想されるエッジ周辺のメッシュを局所的に細かくします。
3. 流入乱流条件の不一致:風洞には固有の乱流強度とスケールがあります。CFDの流入条件を「乱流強度1%、特性長さ0.01m」などと適当に設定せず、風洞の実測値(可能ならば)に合わせます。
4. モデルの忠実度:CFDモデルは支え棒や床面境界層の影響を無視していませんか?風洞実験にはこれらの干扰があるため、CFDでも同様の条件を再現する必要があります。まずは、剥離が起きない低迎角で実験とCFDを一致させ、そこから検証を積み上げていくのが基本です。
関連トピック
なった
詳しく
報告