流体のエネルギー方程式
流体のエネルギー方程式の理論基礎
概要
先生、流体のエネルギー方程式ってどんなときに必要になるんですか？
温度場を求めたい場合に必須だ。熱交換器の設計、電子機器の冷却、燃焼、自然対流などはすべてエネルギー方程式を解く必要がある。等温流れならNS方程式と連続の式だけで閉じるが、温度が絡むとエネルギー方程式が加わる。
エネルギー方程式（温度形）
非圧縮性流れの温度場に対するエネルギー方程式は次の通り。
ここで $c_p$ は定圧比熱、$k$ は熱伝導率、$\Phi$ は粘性散逸関数、$\dot{q}$ は内部発熱率だ。
粘性散逸って何ですか？
粘性によって運動エネルギーが熱に変換される項だ。非圧縮の場合、
通常の工学的な流れでは粘性散逸は無視できるほど小さいが、高粘度流体（ポリマー溶融など）や超高速流では無視できなくなる。
エッカート数
粘性散逸の重要性はエッカート数（Eckert number）で評価する。
$Ec \ll 1$ なら粘性散逸は無視可能。例えば空気流（$U = 50$ m/s, $\Delta T = 20$ K）では $Ec \approx 0.12$ で小さい。一方、ポリマー押出（$U = 0.1$ m/s, $\mu = 1000$ Pa·s）では散逸が温度上昇の主因になり得る。
エネルギー方程式（エンタルピー形）
圧縮性流れではどう変わりますか？
全エンタルピー $h_0 = h + \frac{1}{2}|\mathbf{u}|^2$ を用いたエネルギー方程式が使われる。
圧縮性では状態方程式 $p = \rho R T$（理想気体）が加わり、方程式系が閉じる。密度、速度、圧力、温度の4つの場が全て連成する。
無次元パラメータ
| パラメータ | 定義 | 物理的意味 |
|---|---|---|
| Prandtl数 Pr | $\nu/\alpha = \mu c_p/k$ | 運動量拡散と熱拡散の比 |
| Nusselt数 Nu | $hL/k$ | 対流/伝導の熱伝達比 |
| Peclet数 Pe | $Re \cdot Pr$ | 移流/拡散の比 |
| Eckert数 Ec | $U^2/(c_p\Delta T)$ | 運動エネルギー/熱エネルギー |
Pr数が流体の性質で決まるのに対して、Nuは計算結果として得られるんですね。
その通り。空気は $Pr \approx 0.71$、水は $Pr \approx 7$、エンジンオイルは $Pr \approx 100$〜1000。Pr数が大きいほど熱境界層が薄くなり、壁面近傍の温度勾配を正確に解像するために細かいメッシュが必要になる。
エネルギー方程式の歴史——ジュールの熱の仕事当量実験（1843年）から流体力学へ
流体のエネルギー方程式に含まれる熱の仕事当量の概念は、英国のジェームズ・ジュール（James Joule）が1843年に水の攪拌実験で確立した。1cal=4.186Jというジュール定数は今も熱力学の基礎定数だ。これを流体力学に統合したのはフーリエの熱伝導方程式（1822年）とナビエ-ストークス方程式の結合で、エネルギー方程式の標準形は1850年代に確立した。特に「粘性散逸（Viscous Dissipation）項」——速度勾配が熱に変換される効果——は超音速流や高粘性流体で無視できず、現代のCFDでは高精度解析時には必ず有効にする項だ。低速・低粘性の工学CFDでは省略される場合が多いが、その適用限界（Br=ηU²/(kΔT)>0.1）を認識しておく必要がある。
流体のエネルギー方程式の数値計算手法
エネルギー方程式の離散化
エネルギー方程式はどうやって数値的に解くんですか？
有限体積法では、セルの面フラックスとして移流項と拡散項を離散化する。
特徴的なのは、エネルギー方程式はスカラー輸送方程式であり、速度場が既知であれば線形問題として解ける点だ。NS方程式を解いた後に、後処理的にエネルギー方程式を解くアプローチが可能。
移流スキームの選択
温度場の移流スキームは速度場と同じでいいんですか？
Pe数（Peclet数）が大きいと移流支配になり、中心差分では数値振動が発生する。一般に以下の通り。
| Pe数範囲 | 推奨スキーム | 備考 |
|---|---|---|
| Pe < 2 | 中心差分（CD） | 拡散支配で安定 |
| Pe > 2 | 風上差分（Upwind） | 数値拡散あり |
| 高精度 | QUICK, TVD (MUSCL等) | 精度と安定性の両立 |
熱境界条件
壁面での熱境界条件は主に3種類ある。
| 境界条件 | 数学的表現 | 用途 |
|---|---|---|
| 固定温度（第1種） | $T_{wall} = T_0$ | 冷却水壁面、恒温槽 |
| 固定熱流束（第2種） | $-k\frac{\partial T}{\partial n} = q_w$ | ヒーター、発熱面 |
| 対流伝熱（第3種） | $-k\frac{\partial T}{\partial n} = h(T - T_\infty)$ | 外部環境との熱交換 |
乱流での温度場
乱流の場合、温度場にもモデルが必要ですか？
RANSでは乱流熱フラックス $\overline{u_i'T'}$ のモデリングが必要だ。最も一般的なのは渦拡散率モデル。
$Pr_t$（乱流プラントル数）は通常0.85〜0.9に設定する。壁面近傍の温度プロファイル（Jayatilleke の壁関数等）も重要だ。
乱流Pr数ってどれくらい結果に影響しますか？
壁面熱伝達率に10〜20%程度影響する。特に液体金属（$Pr \ll 1$、$Pr_t \approx 1$〜4）では標準値の0.85では不正確。液体金属の場合は専用のモデルが必要だ。
エネルギー方程式の離散化——「陰解法か陽解法か」は温度の急変速度で決まる
流体エネルギー方程式を数値的に解く際、陽解法（explicit）か陰解法（implicit）かの選択は計算効率を大きく左右します。陽解法は実装が簡単ですが、熱拡散の安定性条件 Δt≤ρcΔx²/(2k) に縛られる。断熱材と金属を混在させたモデルでは熱伝導率が100倍違うため、安定なΔtが最も薄い金属セルに制約され、全体の計算時間が爆増します。このとき陰解法に切り替えると、Δtをずっと大きく取れる。「方程式の形は同じ」でも離散化戦略が変わると計算時間は数十倍変わる——これが数値解法の現実です。
流体のエネルギー方程式の実務適用
実践ガイド
共役熱伝達（CHT）解析って最近よく聞くんですけど、どういうものですか？
流体の対流と固体の熱伝導を同時に解く手法だ。電子基板の冷却設計、エンジンブロック、ガスタービン翼の冷却など、流体-固体界面の温度分布が重要な問題で必須になる。
CHT解析のポイント
流体-固体界面では以下の連続条件が満たされる。
| ソルバー | CHT手法 | 設定方法 |
|---|---|---|
| Fluent | Coupled Wall | 固体領域と流体領域をCoupled Wallで接続 |
| STAR-CCM+ | Multi-Region | リージョン間にinterfaceを定義 |
| OpenFOAM | chtMultiRegionFoam | 各リージョンでソルバーを切替 |
Nusselt数の算出
CFDからNusselt数を計算するにはどうすればいいですか？
壁面熱伝達率 $h$ から算出する。
$T_{\text{ref}}$ は問題に応じて入口温度やバルク温度を使う。円管内流れなら Dittus-Boelter式 $Nu = 0.023 Re^{0.8} Pr^{0.4}$ と比較するのが定石だ。
熱解析のメッシュ要件
温度場の精度にはメッシュ品質が直結する。
| 対象 | 追加要件 | 理由 |
|---|---|---|
| 自然対流 | 壁面付近に密なメッシュ | 浮力駆動の温度境界層 |
| 高Pr流体 | $y^+ \approx 1$ | 熱境界層が速度境界層より薄い |
| CHT | 固液界面のメッシュ整合 | 界面での温度補間精度 |
| 放射 | 視野角の解像 | S2S, DOM等のモデルに依存 |
高Pr数の流体だと、速度場は壁関数で十分でも温度場にはLow-Reが必要ってことですか？
まさにそう。水（$Pr \approx 7$）では熱境界層の厚さが速度境界層の約 $Pr^{-1/3} \approx 0.52$ 倍。エンジンオイル（$Pr \approx 500$）ではさらに薄くなるので、壁面第1セルの解像度が極めて重要になる。
自然対流の注意点
Boussinesq近似では密度変化を浮力項にのみ反映する。
Rayleigh数 $Ra = Gr \cdot Pr = g\beta\Delta T L^3/(\nu\alpha)$ が $10^8$ を超えると乱流遷移が始まる。$\beta\Delta T < 0.1$ がBoussinesq近似の目安だ。
ガスタービン燃焼器の総エンタルピー保存——CFDでの入口条件設定の実際
ガスタービン燃焼器のCFD解析では、圧縮機出口の高圧高温ガスを入口境界条件として正確に設定することが精度の基本だ。エネルギー方程式の実装では「静温度（Static Temperature）」と「総温度（Total Temperature）」の混同が最も多い設定ミスで、高速流（Mach>0.3）では両者が数十℃差になる。また燃料の化学エンタルピー（Lower Heating Value: 43.1MJ/kg for CH4）を熱源項として適切に扱わないと、燃焼後温度が500〜800K過小評価される。実際の設計では入口エンタルピー保存の質量・エネルギーバランスを出口条件と照合し、誤差1%以内を収束判定の目安とする手法が国内航空エンジンメーカーの標準的なCFD検証手順とされている。
流体のエネルギー方程式のソフトウェア比較
エネルギー方程式関連の商用ツール比較
熱解析の機能って、各ソルバーで差がありますか？
基本的なエネルギー方程式はどのソルバーでも解けるが、放射モデル、CHT、相変化の扱いに差がある。
放射モデルの比較
| モデル | Fluent | CFX | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| S2S (Surface to Surface) | 対応 | 対応 | 対応 | viewFactorモデル |
| DOM (Discrete Ordinates) | 対応 | 対応 | 対応 | fvDOM |
| P-1 (球面調和近似) | 対応 | 対応 | 対応 | P1 |
| Monte Carlo | Fluent 2024以降 | 非対応 | 対応 | 非標準 |
どの放射モデルを選べばいいですか？
光学的厚さ $\tau = \kappa L$ で判定する。
- $\tau > 3$（光学的に厚い）: P-1で十分
- $\tau < 0.1$（光学的に薄い）: S2Sが効率的
- それ以外: DOMが汎用的
相変化・沸騰モデル
| 機能 | Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 凝固・融解 | Enthalpy-Porosity法 | 同様 | solidificationMelting |
| 蒸発・凝縮 | Lee model, VOF | VOF + 相変化 | 自作UDF |
| 沸騰 | RPI沸騰モデル | RPI沸騰モデル | 研究用コードあり |
ソルバー別の設定のコツ
各ソルバーでのエネルギー方程式有効化の手順を覚えておこう。
- Fluent: Models > Energy をON。CHT は Cell Zone Conditions > Solid に固体材料を設定
- CFX: Domain > Heat Transfer > Thermal Energy or Total Energy
- STAR-CCM+: Physics > Energy をON。Multi-Region CHT は Region Interface で設定
- OpenFOAM: buoyantSimpleFoam（浮力あり定常）、chtMultiRegionFoam（CHT）を選択
FluentのEnergyモデルをONにすると計算時間はどれくらい増えますか？
温度がスカラー方程式として1つ追加されるだけなので、等温計算と比べて10〜20%増程度。放射モデルを追加するとさらに20〜50%増。CHT（固体領域の追加）の場合はセル数に応じて増加するが、固体は対流項がないので流体より軽い。
エネルギー方程式実装のツール差——FluentとOpenFOAMの熱流束計算の違い
CFDソルバーによってエネルギー方程式の離散化と壁面熱流束の計算方法に微妙な差がある。ANSYS Fluentは壁面熱流束をqw = h_wall × (T_wall - T_fluid_ref)で計算する際、T_fluid_refの定義（壁面セル中心値 vs 対数則補正値）をモデル設定で切り替えられる。一方OpenFOAMのbuoyantSimpleFoamでは壁面境界条件にfixedHeatFluxかexternalWallHeatFluxTemperatureを使い分けるが、設定を誤ると局所的な熱流束が2倍以上ずれるケースがある。特に壁関数（Wall Function）と低Re処理の混在設定では、y+依存の切り替えが自動化されているツールと手動設定のツールで結果が乖離する。信頼性の高い熱流束予測のためには必ず単純な検証ケース（平板強制対流等）でツールのデフォルト設定を確認することが重要だ。
流体のエネルギー方程式の先端研究
先端トピック
エネルギー方程式の先端的な話題にはどんなものがありますか？
温度場モデリングの高度化は産業的にも研究的にも非常に活発な分野だ。
乱流熱フラックスの高度モデリング
標準的な渦拡散率モデル（$Pr_t = 0.85$ 固定）は、以下のケースで不十分。
- 液体金属冷却（$Pr \sim 0.01$）: $Pr_t$ が $Pr$ に依存
- 強い逆圧力勾配下: 乱流熱フラックスと温度勾配が非平行
- 噴流衝突冷却: 淀み点近傍での熱伝達予測
高度なモデルとして、代数的熱フラックスモデル（AFM）やフル2方程式熱フラックスモデル（$\overline{T'^2}$-$\varepsilon_\theta$ モデル）がある。
LESでの温度場
LESでは温度のサブグリッドスケール(SGS)フラックスのモデルが必要。
$Pr_{\text{sgs}} \approx 0.4$〜0.6が典型的。Dynamic SGSモデルでは $Pr_{\text{sgs}}$ も動的に計算できる。
超臨界流体の伝熱
超臨界CO2サイクルが注目されてますけど、CFDで扱えるんですか？
超臨界状態では擬臨界点近傍で物性値（$c_p$, $\rho$, $k$, $\mu$）が急激に変化する。NIST REFPROPやCoolPropなどの物性ライブラリと連携して、温度・圧力依存の物性テーブルを使う必要がある。
標準的なRANS乱流モデルは超臨界伝熱の予測精度が低いことが知られている。壁面近傍の密度変化が乱流構造に影響するためだ。
燃焼と化学反応
燃焼ではエネルギー方程式に化学反応の発熱項が加わる。
$h_k^0$ は種 $k$ の生成エンタルピー、$\dot{\omega}_k$ は反応速度。Flamelet model、EDC（Eddy Dissipation Concept）、PDF輸送方程式など、多数の燃焼モデルがある。
エネルギー方程式が関わる問題は、本当に多岐にわたるんですね。
そう。温度が絡まない工学問題はむしろ少数派だ。だからこそエネルギー方程式の物理を深く理解しておくことが重要だよ。
データセンターの「液浸冷却」——流体エネルギー方程式が再び主役に
近年急速に普及する液浸冷却（immersion cooling）では、サーバー基板を絶縁油や冷却液に直接漬ける。このとき流体エネルギー方程式の対流項と熱伝導項のバランスが、CPU温度を決定します。空冷と比べて熱伝達率が20〜50倍になるため、エネルギー方程式の係数がまったく異なるスケールになる。最先端の2相液浸（相変化冷却）では、さらに蒸発・凝縮の潜熱項が加わり、エネルギー方程式の形が大きく変わる。流体エネルギー方程式の先端は今、AI半導体の冷却設計と直結しています。
流体のエネルギー方程式のトラブル対応
トラブルシューティング
温度場の計算でよくあるトラブルを教えてください。
温度が発散する、壁面熱伝達率が合わない、という問題が多い。
1. 温度が非物理的な値に発散
症状: 温度が$10^{10}$ Kなどの異常値になる。
| 原因 | 対策 |
|---|---|
| エネルギー残差が未収束 | 緩和係数を下げる（0.8→0.5） |
| 粘性散逸を有効化したが速度場が不正確 | 速度場を先に収束させてからエネルギーON |
| 放射モデルの吸収係数が不適切 | 光学的性質を見直す |
| 温度依存物性のテーブル範囲外 | 物性テーブルの温度範囲を拡張 |
2. Nusselt数が実験相関と合わない
円管内流れのNuが Dittus-Boelter式と30%もずれるんですが…
| 原因 | 診断 | 対策 |
|---|---|---|
| y+が不適切 | 壁面y+を確認 | 熱伝達には $y^+ < 5$ 推奨 |
| 乱流Pr数 | デフォルト0.85か確認 | Pr > 10の液体では感度解析 |
| 助走区間 | 入口から発達区間を除外しているか | $L/D > 10$ を確保 |
| 参照温度 | $T_{\text{ref}}$ の定義 | バルク温度を使用 |
Dittus-Boelter式自体にも $\pm 25$%の誤差がある。CFDとの比較では $\pm 10$%の一致で良好と判断してよい。
3. 自然対流が発生しない
症状: 加熱壁面があるのに流れが全く発生しない。
原因:
- 重力が未設定（$g = 0$）
- Boussinesq近似をONにしていない（密度が一定のまま）
- Operating Density の設定が不適切
対策: Fluent なら Operating Conditions > Gravity を設定し、材料物性で Density > Boussinesq を選択。Operating Density は参照温度の密度に設定する。
4. CHT界面の温度不連続
固体と流体の境界で温度がジャンプしてるんですけど…
原因: 界面のメッシュが整合していない（conformal でない）、またはinterface condition の設定ミス。
対策:
- 流体と固体のメッシュ界面を一致させる（conformal mesh）
- Non-conformal の場合は interface の Mapped 設定を確認
- 壁面の Thermal Condition が Coupled に設定されているか確認
温度場のトラブルは、流れ場の計算が正しいことが大前提なんですね。
その通り。温度場は速度場に従属するから、まず速度場と圧力場を収束させ、その上でエネルギー方程式のデバッグに移るべきだ。
「エネルギーが保存されない」バグは電力業界で大問題になった
1970年代、原子力発電所の冷却系CFD解析で「計算上のエネルギー収支が現実より5%以上ずれる」問題が頻発しました。原因の多くは粘性散逸項の符号ミスか、対流項と拡散項の離散化不整合。5%のエラーが冷却材温度の誤評価につながり、安全評価に直結するため、当時の規制当局は計算コードの検証を義務化しました。流体エネルギー方程式のトラブルは、理論ミスより実装ミスが圧倒的に多い——という教訓が今も生きています。
関連トピック
なった
詳しく
報告