プロペラCFD解析 — 典型的な問題と対処法
$K_T$ が実験値と10%以上ずれる
オープンウォーター計算で $K_T$ が実験値と合いません。
まず以下を確認しよう。(1) 翼面メッシュが前縁曲率を十分に解像しているか(前縁に10セル以上)。(2) $y^+ < 1$ を達成しているか。特にSST k-ωでは $y^+$ が5〜30の中途半端な値だと精度が悪化する。(3) 翼端の後流が十分に解像されているか。メッシュ感度分析で収束を確認する。
MRFとSliding Meshで結果が違うんですが、どちらが正しいですか?
オープンウォーター条件(一様流入)ならMRFとSliding Meshの定常平均値はほぼ一致するはず。差が5%以上あるなら、MRFのインターフェース面でのメッシュ不整合が原因の可能性がある。FluentではMesh InterfaceのInterpolation方法をConservativeに設定し、orphan cellsがないか確認しよう。
Sliding Meshで圧力が振動する
Sliding Mesh計算で圧力が大きく振動して発散しそうです。
時間刻みが大きすぎる可能性がある。プロペラが1ステップで回転する角度を1度以下にするのが目安。$n = 10$ rpsなら $\Delta t < 1/(360 \times 10) \approx 2.8 \times 10^{-4}$s。また、sliding interface上のメッシュサイズが回転域と固定域で大きく異なると補間エラーが出る。interface面のセルサイズを揃えるのが重要だ。
キャビテーション解析が不安定
キャビテーション解析が数回転で発散します。
キャビテーションモデルの追加は数値的に不安定になりやすい。(1) まず非キャビテーション条件で十分に収束した流れ場を初期値にする。(2) 時間刻みをさらに小さくする($\Delta t \sim 10^{-5}$s)。(3) VOFのunder-relaxationを0.3程度に下げる。(4) Courant数を0.5以下に抑える。STAR-CCM+のAutomatic CFL Ramping機能は初期の不安定期を自動的に乗り越えるのに有効だよ。
レイノルズの実験(1883年)——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
トラブル解決の考え方
デバッグのイメージ
CFDのデバッグは「水道管の詰まり修理」に似ている。まず「どこで詰まっているか」(どの残差が下がらないか)を特定し、次に「何が詰まっているか」(メッシュ品質?境界条件?乱流モデル?)を調べ、最後に「どう直すか」(メッシュ修正?緩和係数?)を判断する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——プロペラCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
プロペラCFD解析の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →