強制対流のCFD解析
理論と物理
強制対流の基本概念
先生、強制対流って「ファンやポンプで流体を動かして冷やす」ことですよね？ CFDで解くときのポイントは何ですか？
そのとおり。強制対流では外部から与えられた流れ（ファン、ポンプ、配管内圧力差など）が主体で、浮力は無視できるほど小さい。CFDでは速度場・圧力場を解いた上で、壁面の熱流束やNusselt数を正確に予測することが目標になる。
Nusselt数って何を意味するんですか？
Nusselt数 $Nu$ は壁面の対流伝熱と流体の熱伝導の比で、伝熱性能の無次元指標だ。
ここで $h$ は熱伝達係数、$L$ は代表長さ、$k_f$ は流体の熱伝導率。$Nu$ が大きいほど対流伝熱が支配的で、効率的に熱を運んでいるということだ。
代表的なNu数相関式
実験的な相関式にはどんなものがありますか？
管内流れでよく使われるのがDittus-Boelter式とGnielinski式だ。
$n = 0.4$（加熱）、$n = 0.3$（冷却）。適用範囲は $Re_D > 10{,}000$、$0.6 < Pr < 160$、$L/D > 10$。
Gnielinski式はDittus-Boelter式より広い範囲（$2300 < Re_D < 5 \times 10^6$）で精度が良い。CFD結果の検証にはまずこれらの相関式と照合するのが定石だよ。
外部流れの場合は？
平板上の強制対流層流では $Nu_x = 0.332 Re_x^{1/2} Pr^{1/3}$ 、乱流では $Nu_x = 0.0296 Re_x^{0.8} Pr^{1/3}$ が基本だ。円柱まわりの流れではChurchill-Bernstein式が使われる。
Nusselt数相関式は「実験式の化石」——今も設計現場で使い続けられる理由
Dittus-Boelterの式など古典的Nusselt数相関式は1930〜50年代に実験で求められた経験式だ。CFDが普及した今もなお、設計初期段階での熱伝達率の見積もりに第一線で使われている。理由は速さと手軽さ——CFD計算を回す前に「大体この程度」と感覚を掴むのに相関式は最強だ。ただし適用範囲（Re数・Pr数・形状）を外れた使い方をすると50%以上外れることもある。式を使う前に「この条件で使っていいのか」を確認する習慣が、実務エンジニアの基本です。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
壁面処理の選択
強制対流CFDでは壁面メッシュの処理が特に重要だと聞きました。
その通り。壁面の熱伝達を正確に予測するには、温度境界層を十分に解像する必要がある。2つの主要アプローチがある。
(1) Low-Reynolds number approach: 壁面第一セルを $y^+ \approx 1$ に配置し、粘性底層から直接解く。SST k-ωモデルと組合せるのが標準。
(2) 壁面関数（Wall Function）approach: $y^+ \approx 30$〜$300$ に第一セルを配置し、壁面近傍を対数則で近似する。Fluentの Enhanced Wall Treatment は $y^+$ によって自動切り替えするので実用的だ。
熱伝達予測にはどちらがいいんですか？
定量的なNu数予測が必要なら $y^+ \approx 1$ のLow-Re approachを推奨する。壁面関数は流れの傾向把握には使えるけど、Nu数を10〜20%過大/過小評価するリスクがある。特に剥離を伴う流れ（後ろ向きステップ、ブラフボディ）では壁面関数の精度が悪化するよ。
乱流モデルの比較
どの乱流モデルが強制対流に向いていますか？
| 乱流モデル | 管内流れ | 外部流れ | 剥離流れ | 計算コスト |
|---|---|---|---|---|
| 標準k-ε | 可 | 可 | 不可 | 低 |
| Realizable k-ε | 良 | 良 | 可 | 低 |
| SST k-ω | 良 | 良 | 良 | 低〜中 |
| Transition SST | 遷移域で良 | 遷移域で良 | 良 | 中 |
| k-ω BSL RSM | 非等方性重要時 | 良 | 良 | 高 |
管内流れの発達領域なら標準k-εでもいいですか？
完全発達した円管内乱流のNu数予測なら、標準k-εでもDittus-Boelter式と5%以内で一致する。ただし入口領域や急縮小・急拡大のある管路ではSST k-ωのほうが安全だ。遷移領域（$Re \approx 2300$〜$10000$）を扱うならTransition SSTモデルが必要になる。
壁面処理の選択ミスで熱伝達率が2倍ズレる
強制対流CFDで最も影響が大きいのが壁面乱流処理の選択だ。壁関数を使うか、低Re数モデルで壁近傍を解像するかで熱伝達係数が大きく変わる。あるカーエアコン設計案件で、担当者が壁関数（y+≈30を想定）のまま粗いメッシュで解析したところ、実測より熱伝達係数が50%以上高く出た。実際のy+は200を超えており、壁関数の適用範囲を完全に逸脱していた。「とりあえず壁関数」は強制対流で最もコストが高い誤解の一つです。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
ヒートシンクのCFD解析
電子機器のヒートシンク設計にCFDを使うケースが多いと聞きました。
そうだね。アルミ押出しフィン、ピンフィン、液冷コールドプレートなど、いずれも強制対流による放熱が主体だ。設計変数としてはフィンピッチ、フィン高さ、流路幅、流量などがあり、CFDで熱抵抗 $R_{th} = \Delta T / Q$ を評価してパラメトリックスタディを行う。
具体的なワークフローを教えてください。
(1) フィンの対称性を利用して1ピッチ分だけモデル化する（周期境界条件を使用）。(2) 入口に一様流速、出口に圧力出口を設定。(3) フィン底面に一定熱流束を付与。(4) SST k-ωで定常計算。(5) フィン底面の平均温度から熱抵抗を算出。
1ピッチ分で済むのは計算効率がいいですね。流路の発達領域の影響は？
短いヒートシンク（$L/D_h < 20$）では入口効果が無視できないので、全体モデルが必要だ。長いヒートシンクで完全発達流れが主体なら1ピッチ+周期条件で十分。Fluentではperiodic conditionにmass flow rateを指定できるし、STAR-CCM+でもperiodic interfaceで設定可能だ。
液冷コールドプレートの設計
パワー半導体の液冷設計もCFDですか？
液冷コールドプレートは典型的なCHT問題だ。アルミや銅のプレートに微細流路（マイクロチャネル）を加工し、水やクーラントを流す。流路の水力直径 $D_h$ が1mm以下になるとRe数が低くなり層流になることも多い。その場合は層流モデルで十分だ。
圧力損失も重要ですよね？
非常に重要。ポンプの動力制約があるから、熱抵抗と圧力損失のトレードオフを最適化する必要がある。パレートフロント（多目的最適化）の考え方が有効で、STAR-CCM+のDesign ManagerやFluentのDesign ExplorationツールでDOEベースの最適化ができるよ。
Ergun式のような圧損相関もCFD結果の検証に使えますか？
直線流路ならDarcy-Weisbach式 $\Delta p = f \cdot (L/D_h) \cdot (\rho u^2/2)$ で検証できる。層流なら $f = 64/Re$ 、乱流ならMoody線図の値と比較する。CFDの圧損がこれらと10%以上ずれていたらメッシュか設定に問題がある可能性が高い。
ヒートシンクのフィン形状は「最適化し尽くされた」と思ったら……
アルミ押出しのストレートフィンヒートシンクは1980年代からほとんど形が変わらず「枯れた技術」に見える。ところが強制対流CFDを使った形状最適化を行うと、従来のストレートフィンより20〜30%熱抵抗を下げられる波型やオフセットフィン構造が見つかることが多い。加工コストとのトレードオフで採用されていなかっただけで、「最適解」だったわけではない。3Dプリンティングで複雑形状が安く作れるようになった今、ヒートシンク形状は再び進化が始まっています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
Ansys Fluentでの設定
Fluentで強制対流伝熱を解くときの推奨設定を教えてください。
定常・非圧縮性・乱流の標準ケースでの推奨を挙げよう。
| 設定項目 | 推奨値 | 理由 |
|---|---|---|
| Solver | Pressure-Based, Coupled | エネルギー方程式との連成収束性が良い |
| 乱流モデル | SST k-ω | 壁面近傍の精度が高い |
| 壁面処理 | Enhanced Wall Treatment | $y^+$値によらず適応的に切替え |
| 離散化 | Second Order Upwind (全変数) | 精度確保のため |
| 圧力補間 | Second Order | Body Force Weighted (浮力がある場合) |
| 収束基準 | 残差 $10^{-6}$ (Energy), $10^{-4}$ (他) | エネルギーは厳しめに |
CoupledソルバーとSIMPLEソルバーで精度は変わりますか？
最終的な解の精度は同等だけど、Coupledのほうが収束が速い。特にCHT問題や浮力が効く場合はCoupledが安定する。SIMPLEは反復ごとのメモリ使用量が少ないので、メモリ制約がある場合に選択する。
OpenFOAMでの設定
OpenFOAMならどうですか？
ソルバーはbuoyantSimpleFoam（定常、浮力項あり）またはsimpleFoam+エネルギー方程式追加が標準。乱流モデルはconstant/turbulencePropertiesでkOmegaSSTを指定する。壁面境界条件はkqRWallFunction、omegaWallFunction、nutUSpaldingWallFunctionの3点セットだ。
fvSchemesの離散化はどうすればいいですか？
divスキームはbounded Gauss linearUpwindが安定と精度のバランスが良い。Gauss linearは精度が高いけど非定常流れでは振動しやすい。gradスキームはGauss linearで、limiterを入れる場合はcellLimitedを使おう。
STAR-CCM+での設定
STAR-CCM+の場合は？
Physics ModelでSegregated Flowを選択し、エネルギーモデルをSegregated FluidTemperatureで有効化する。乱流はSST k-ω + All y+ Wall Treatment（$y^+$ によらず自動対応）が推奨。Mesh ModelではPolyhedral Mesher + Prism Layer Mesherの組合せが品質と効率の両立に優れているよ。
強制対流CFDツール比較——STAR-CCM+ のコンジュゲートHTモードの実力
強制対流解析では共役熱伝達（CHT: Conjugate Heat Transfer）機能の使いやすさがツール選定の鍵となる。STAR-CCM+はソリッドと流体を単一メッシュ上で連成解析できる「In-situ CHT」が強みで、界面の熱抵抗（接触抵抗）もGUI上で簡単に設定できる。一方ANSYS FluentはFluent-Mechanicsのシステム連成が必要なケースが多く、ライセンスコストが上がる。OpenFOAMはchtMultiRegionFoamで同等の連成が無償で可能だが、初期セットアップに熟練を要する。自動車部品やモーター冷却など量産品の熱設計には、CADとの統合性を含めたトータルワークフローで選ぶべきだ。
選定で最も重要な3つの問い
- 「何を解くか」：強制対流のCFD解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
層流-乱流遷移の扱い
Re数が中間的な領域で遷移が起こる場合、どうモデリングすればいいですか？
$Re \approx 2300$〜$10000$ の遷移領域はRANS乱流モデルの弱点だ。通常のk-εやk-ωでは遷移の位置を予測できない。Ansys FluentのTransition SSTモデル（$\gamma$-$Re_{\theta}$モデル）がこの問題に対応しており、遷移の開始位置と長さを実験的相関式に基づいて予測する。
STAR-CCM+にも同等のモデルはありますか？
STAR-CCM+にもGamma Transition Modelがある。OpenFOAMではkOmegaSSTLM（Langtry-Menter遷移モデル）が利用可能だ。いずれも自由流れ乱流強度 $Tu$ と自由流れ粘性比 $\mu_t/\mu$ の入口条件が遷移位置に大きく影響するので、これらを実験条件に合わせて正確に設定することが重要だよ。
乱流増強構造（タービュレータ）
リブやディンプルで伝熱を促進する設計もCFDで評価するんですか？
もちろん。内燃機関の冷却通路やガスタービンの内部冷却ではリブ付き流路が標準的だ。リブの高さ $e$、ピッチ $P$、角度 $\alpha$ が設計パラメータで、Nu数の増大比（smooth管比）とfriction factorの増大比の関係を評価する。
リブ付き流路のCFDで注意すべきなのは、リブ間の再循環域を正確に捉えること。SST k-ωは再付着点の位置を比較的よく予測するが、標準k-εは再付着長さを過小評価する傾向がある。メッシュはリブ高さ方向に20層以上配置し、リブ間を流れ方向に少なくとも40分割するのが目安だ。
熱流体連成最適化
フィン形状やリブ形状の最適化もできますか？
トポロジー最適化がホットなトピックだ。Ansys Fluent 2023R2以降ではAdjoint Topology Optimizationが使えるようになり、流路形状を目的関数（熱抵抗最小化+圧損制約など）に基づいて自動的に進化させられる。COMSOL MultiphysicsのTopology Optimization Moduleも同様の機能を持っている。得られた形状は有機的で従来設計では思いつかないような構造になることが多いよ。
その形状って製造できるんですか？
3Dプリンティング（SLMやDMLS）なら製造可能だ。実際にAdditive Manufacturingを前提としたトポロジー最適化ヒートシンクの研究事例が増えている。製造制約（最小壁厚、オーバーハング角度など）を最適化に組み込むのが実用上のポイントだよ。
Coffee Break よもやま話
「再層流化」——強制対流で乱流が層流に戻る不思議な現象
強制対流では一度乱流になった流れが、局所的な加速などで層流に戻る「再層流化（relaminarization）」が起きることがある。例えばガスタービン翼の冷却通路内で、コーナー付近の強い加速で乱流が消えてしまい、予測していたより熱伝達が悪化するケースが報告されている。一般的な乱流モデル（k-εなど）はこの現象を捉えにくく、設計マージンとして補正係数を掛けることが多い。「乱流は一方通行で増える」という思い込みが危ない、という教訓を与えてくれます。
Re数が中間的な領域で遷移が起こる場合、どうモデリングすればいいですか？
$Re \approx 2300$〜$10000$ の遷移領域はRANS乱流モデルの弱点だ。通常のk-εやk-ωでは遷移の位置を予測できない。Ansys FluentのTransition SSTモデル（$\gamma$-$Re_{\theta}$モデル）がこの問題に対応しており、遷移の開始位置と長さを実験的相関式に基づいて予測する。
STAR-CCM+にも同等のモデルはありますか？
STAR-CCM+にもGamma Transition Modelがある。OpenFOAMではkOmegaSSTLM（Langtry-Menter遷移モデル）が利用可能だ。いずれも自由流れ乱流強度 $Tu$ と自由流れ粘性比 $\mu_t/\mu$ の入口条件が遷移位置に大きく影響するので、これらを実験条件に合わせて正確に設定することが重要だよ。
リブやディンプルで伝熱を促進する設計もCFDで評価するんですか？
もちろん。内燃機関の冷却通路やガスタービンの内部冷却ではリブ付き流路が標準的だ。リブの高さ $e$、ピッチ $P$、角度 $\alpha$ が設計パラメータで、Nu数の増大比（smooth管比）とfriction factorの増大比の関係を評価する。
リブ付き流路のCFDで注意すべきなのは、リブ間の再循環域を正確に捉えること。SST k-ωは再付着点の位置を比較的よく予測するが、標準k-εは再付着長さを過小評価する傾向がある。メッシュはリブ高さ方向に20層以上配置し、リブ間を流れ方向に少なくとも40分割するのが目安だ。
フィン形状やリブ形状の最適化もできますか？
トポロジー最適化がホットなトピックだ。Ansys Fluent 2023R2以降ではAdjoint Topology Optimizationが使えるようになり、流路形状を目的関数（熱抵抗最小化+圧損制約など）に基づいて自動的に進化させられる。COMSOL MultiphysicsのTopology Optimization Moduleも同様の機能を持っている。得られた形状は有機的で従来設計では思いつかないような構造になることが多いよ。
その形状って製造できるんですか？
3Dプリンティング（SLMやDMLS）なら製造可能だ。実際にAdditive Manufacturingを前提としたトポロジー最適化ヒートシンクの研究事例が増えている。製造制約（最小壁厚、オーバーハング角度など）を最適化に組み込むのが実用上のポイントだよ。
「再層流化」——強制対流で乱流が層流に戻る不思議な現象
強制対流では一度乱流になった流れが、局所的な加速などで層流に戻る「再層流化（relaminarization）」が起きることがある。例えばガスタービン翼の冷却通路内で、コーナー付近の強い加速で乱流が消えてしまい、予測していたより熱伝達が悪化するケースが報告されている。一般的な乱流モデル（k-εなど）はこの現象を捉えにくく、設計マージンとして補正係数を掛けることが多い。「乱流は一方通行で増える」という思い込みが危ない、という教訓を与えてくれます。
トラブルシューティング
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設定がこの用途に適している。
強制対流CFDで「Nu数が実験の半分」になるよくある理由
強制対流CFDの検証で「Nusselt数が実験値の半分しか出ない」というトラブルは珍しくない。原因として多い順番に挙げると：(1) 流入境界条件の乱流強度の設定漏れ、(2) 温度境界条件を断熱にしたまま熱流束を別途入れた二重設定、(3) 流れが十分に発達する前の位置で熱伝達係数を評価している——の3パターンだ。特に(1)は「乱流モデルを使えば乱流が自動的に生まれる」という誤解から来る。乱流は入口で「種を蒔く」必要があります。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——強制対流のCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告