強制対流のCFD解析 — 熱伝達予測のトラブルシューティング
Nu数がDittus-Boelterと合わない
管内流れのCFDでNu数を計算したら、Dittus-Boelter式より30%も高いです。
まずメッシュと流れの状態を確認しよう。完全発達しているか？ 入口から十分な助走距離（乱流なら $L > 10D$ 程度）をとっているか？ 助走区間ではNu数が発達値より高くなるのは正常だ。
発達区間で比較してます。
次に壁面処理を確認しよう。壁面関数を使っていて $y^+$ が5〜30の「バッファ層」に入っていないか？ バッファ層は壁面関数の適用範囲外で、熱伝達の予測精度が悪化する。Fluentなら$y^+$のコンター図を確認して、30以上（壁面関数用）か1以下（Low-Re用）のどちらかに統一すべきだ。
温度分布に非対称性が出る
軸対称のはずの管内流れで温度分布が非対称になります。
メッシュの非対称性が原因であることが多い。テトラメッシュでは面の向きがランダムなので、粗いメッシュだと非対称な数値拡散が生じる。ヘキサメッシュに変更するか、メッシュを十分に細かくすれば改善する。また、Second Order Upwindを使っているか確認しよう。First Orderだと数値拡散が大きく非対称になりやすい。
物性値の温度依存性
冷却水の物性値って温度依存性を入れるべきですか？
温度差が大きい場合（$\Delta T > 20$度C程度）は入れるべきだ。特に粘度 $\mu$ は温度依存性が大きく、壁面近傍の速度勾配ひいてはNu数に直結する。Fluentではmaterial propertiesでpolynomialやpiecewise-linearで温度依存性を入力できる。OpenFOAMではtransportProperties内のtransportModelをSutherlandやpolynomialに設定する。
理想気体近似はどこまで使えますか？
空気や窒素を冷媒とする場合、Ma数が0.3以下（圧縮性効果が5%以下）なら非圧縮性近似で十分。ただし大きな温度差がある場合はBoussinesq近似ではなく理想気体として密度の温度依存性を扱ったほうがよい。FluentのIncompressible Ideal Gas設定がこの用途に適している。
レイノルズの実験（1883年）——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
トラブル解決の考え方
デバッグのイメージ
CFDのデバッグは「水道管の詰まり修理」に似ている。まず「どこで詰まっているか」（どの残差が下がらないか）を特定し、次に「何が詰まっているか」（メッシュ品質？境界条件？乱流モデル？）を調べ、最後に「どう直すか」（メッシュ修正？緩和係数？）を判断する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——強制対流のCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
CAEの未来を、実務者と共に考える
Project NovaSolverは、強制対流のCFD解析における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →