自然対流のCFD解析 — 収束困難と検証の実際
定常計算が収束しない
自然対流の定常計算で残差が振動して収束しません。
自然対流では物理的に非定常な流れが発生することが多い。残差振動は「定常解が存在しない」ことを意味している場合がある。$Ra > 10^8$（密閉キャビティ）や $Ra > 10^9$（垂直平板）では非定常計算に切り替えるべきだ。
定常計算のまま収束させるテクニックはありますか？
(1) FluentのPseudo Transient設定を有効にする。(2) under-relaxation factorsを大幅に下げる（圧力0.2、運動量0.3程度）。(3) First Order Upwindで初期収束させてからSecond Orderに切り替える。ただしこれらは人工的に定常化しているだけなので、物理的妥当性の確認が必要だよ。
Nu数がde Vahl Davisのベンチマーク値と合わない
差分加熱キャビティで $Ra = 10^6$ のNu数がベンチマークの8.8ではなく6.5になります。
メッシュが粗すぎる可能性が高い。$Ra = 10^6$ では壁面近傍の境界層が非常に薄く、不均一メッシュでの壁面第一層が厚すぎると熱流束を過小評価する。壁面法線方向に最低20層、壁面第一層厚さが境界層の1/10以下になるようにメッシュを調整しよう。
メッシュを細かくしたら8.6になりました。まだ2%ずれています。
Second Order Upwindを使っているか確認。First Orderだと数値拡散でNu数が低く出る。また、Boussinesqモデルの参照温度 $T_0$ が高温壁と低温壁の平均温度に設定されているか確認。参照温度の不適切な設定はNu数に影響する。
輻射モデルを入れると結果が大きく変わる
輻射モデルを追加したらNu数が50%増加しました。妥当ですか？
壁面放射率（emissivity）の設定を確認しよう。金属光沢面なら $\varepsilon \approx 0.1$、塗装面なら $\varepsilon \approx 0.9$。デフォルトの $\varepsilon = 1.0$（完全黒体）のままだと輻射を過大評価する。実際の壁面材料に合わせて設定すべきだ。
実験と比較する場合、輻射は入れるべきですか？
実験データが対流成分のみのNu数を報告しているのか、輻射込みの全熱伝達を報告しているのかを確認する必要がある。多くの実験では輻射を分離できていないので、CFDでも輻射を含めて比較するのが妥当だ。
レイノルズの実験（1883年）——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
トラブル解決の考え方
デバッグのイメージ
CFDのデバッグは「水道管の詰まり修理」に似ている。まず「どこで詰まっているか」（どの残差が下がらないか）を特定し、次に「何が詰まっているか」（メッシュ品質？境界条件？乱流モデル？）を調べ、最後に「どう直すか」（メッシュ修正？緩和係数？）を判断する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——自然対流のCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
自然対流のCFD解析の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →