流体のエネルギー方程式
理論と物理
概要
先生、流体のエネルギー方程式ってどんなときに必要になるんですか？
温度場を求めたい場合に必須だ。熱交換器の設計、電子機器の冷却、燃焼、自然対流などはすべてエネルギー方程式を解く必要がある。等温流れなら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）を認識しておく必要がある。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
エネルギー方程式の離散化
エネルギー方程式はどうやって数値的に解くんですか？
有限体積法では、セルの面フラックスとして移流項と拡散項を離散化する。
特徴的なのは、エネルギー方程式はスカラー輸送方程式であり、速度場が既知であれば線形問題として解ける点だ。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をずっと大きく取れる。「方程式の形は同じ」でも離散化戦略が変わると計算時間は数十倍変わる——これが数値解法の現実です。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
共役熱伝達（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検証手順とされている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ 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+依存の切り替えが自動化されているツールと手動設定のツールで結果が乖離する。信頼性の高い熱流束予測のためには必ず単純な検証ケース（平板強制対流等）でツールのデフォルト設定を確認することが重要だ。
選定で最も重要な3つの問い
- 「何を解くか」：流体のエネルギー方程式に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピック
エネルギー方程式の先端的な話題にはどんなものがありますか？
温度場モデリングの高度化は産業的にも研究的にも非常に活発な分野だ。
乱流熱フラックスの高度モデリング
標準的な渦拡散率モデル（$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%のエラーが冷却材温度の誤評価につながり、安全評価に直結するため、当時の規制当局は計算コードの検証を義務化しました。流体エネルギー方程式のトラブルは、理論ミスより実装ミスが圧倒的に多い——という教訓が今も生きています。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——流体のエネルギー方程式の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告