凝縮熱伝達 — CAE用語解説
凝縮熱伝達
先生、凝縮熱伝達って沸騰熱伝達と逆のプロセスですよね？ 何が違うんですか？
凝縮熱伝達の理論基礎
凝縮熱伝達の基本概念
凝縮熱伝達って、単純に「蒸気が冷やされて水になる時の熱移動」で合ってますか？
本質的にはそうですが、CAEで扱う際は「相変化に伴う潜熱の放出」と「凝縮液膜による熱抵抗」の2点が鍵になります。例えば、飽和温度100℃の水蒸気が1kg凝縮すると、約2257kJの潜熱を放出します。この熱量が壁面を通して奪われるプロセスをモデル化する必要があります。
「凝縮液膜」が熱抵抗になるのはなぜですか？ 液体の方が気体より熱伝導率が高いのでは？
良い着眼点です。確かに水の熱伝導率（約0.6 W/m·K）は水蒸気（約0.02 W/m·K）より高い。しかし、問題は「膜の厚さ」です。壁面に形成される液膜は、蒸気と壁面の間に「静止した液体の層」として立ちはだかり、熱がこの層を「伝導」で通過しなければなりません。熱伝達率は
支配方程式は、普通の熱伝導や対流とどう違うんですか？
エネルギー保存則に「相変化に伴う潜熱の項」が追加され、また液膜と蒸気の界面で「質量と運動量のジャンプ条件」が必要になります。単相流のエネルギー式に、蒸発・凝縮による質量流量ṁと潜熱h_fgを考慮した項が加わります：
凝縮熱伝達の数値計算手法
VOF法とEulerian多相流
凝縮をCFDでシミュレーションする時、一番よく使われるアプローチは何ですか？
大きく分けて2つ。界面を追跡する「VOF (Volume of Fluid)法」と、各相を連続体として扱う「Eulerian多相流モデル」です。熱交換器のチューブ内凝縮のように、液膜の形状やドレインが重要な場合はVOF法が選ばれます。Ansys Fluentでは「Evaporation-Condensation」モデルをVOFと組み合わせて使います。
VOF法で凝縮をモデル化する時の、具体的なソース項の与え方は？
界面近傍のセルで、温度に基づく質量移動速度を定義します。よく使われるのはLeeモデルで、蒸発・凝縮の質量流量ṁは：
調整係数rを決める基準はありますか？ でたらめに選ぶとまずいですよね。
その通りです。rは「界面での相変化が物理的に妥当な時間スケールで起こる」ように設定します。目安として、rの逆数（時定数）が流体の滞留時間より十分に小さくなる値にします。例えば滞留時間が0.1秒なら、時定数を0.01秒（r=100）程度に設定します。しかし、これは経験則であり、Nusseltの膜凝縮理論など既知の実験結果や理論解と照合して検証することが必須です。
凝縮熱伝達の実務適用
モデル構築と検証のワークフロー
凝縮熱伝達のシミュレーションを始める時、最初に何を確認すべきですか？
まず「凝縮モード」を特定します。「膜凝縮」か「滴状凝縮」か。工業機器の多くは膜凝縮です。次に、物性値、特に飽和温度と圧力の関係（蒸気表）、潜熱、表面張力を正確に設定します。物性が温度依存する場合は多項式で定義するか、NIST REFPROPなどの外部データベースと連携できるか確認します（COMSOLでは可能）。
メッシュはどのように切るべきですか？ 特に気をつける点は？
液膜が存在する・または形成される壁面付近のメッシュを極めて細かくする必要があります。膜厚は多くの場合サブミリメートルオーダーです。例えば、水平管外の凝縮では、Nusselt理論から予想される平均液膜厚さ（例えば0.1mm）を少なくとも3〜5層のメッシュで解像できるようにします。また、界面のカーブを捉えるため、界面が通過する領域全体のメッシュも細かくします。Ansys Meshingで「Inflation Layer」を適用するのが一般的です。
収束判定はどうすれば？ 単相流より難しい気がします。
非常に難しいです。残差のモニターに加え、「質量保存のチェック」が必須です。計算領域全体の蒸気相と液相の質量を積分し、流入する蒸気質量流量と流出する凝縮水の質量流量、そして領域内の蒸気質量減少率がバランスしているかを確認します。また、代表的な壁面の熱流束や熱伝達率が時間平均で定常値に収束するかも見ます。不安定な場合は、まず疑似時間ステップを小さく（例えば1e-4秒）して、ゆっくり始動させます。
凝縮熱伝達のソフトウェア比較
各ソルバーの特徴と選択指針
Ansys Fluent、STAR-CCM+、COMSOLで凝縮モデルを扱う場合、大きな違いは何ですか？
コアとなる物理モデルは似ていますが、実装の容易さとカスタマイズ性に差があります。FluentとSTAR-CCM+は専用の「相変化モデル」がプリセットされており、GUIから比較的簡単に設定できます。特にSTAR-CCM+の「Eulerian Multiphase」に組み込まれた「Thermal Phase Change」モデルは、液滴の形成・合体も考慮でき、強力です。一方、COMSOLは「Two-Phase Flow, Phase Field」や「Level Set」モジュールに「相変化」機能があり、支配方程式を直接いじるカスタマイズ性が高い代わりに、ユーザーが物理を深く理解している必要があります。
Abaqusで凝縮を扱うことはできますか？
Abaqus/Standardの「熱-電気-構造」連成解析や、湿気拡散解析はできますが、流体の流れを伴う凝縮（CFD）そのものをネイティブで解くことはできません。ただし、Abaqus/CFD（旧称：FEATool）や、Dassault Systèmesの別製品であるXFlow（粒子法ベースのCFD）を使う選択肢はあります。あるいは、凝縮による熱流束を経験式（例えばNusseltの理論式）で計算し、それをAbaqusの熱伝導解析の境界条件として与える「簡易的な連成」は可能です。
オープンソースではどうでしょう？ OpenFOAMで凝縮解析は可能ですか？
可能です。`interCondensatingFoam`や`chtMultiRegionFoam`をベースにしたソルバーがコミュニティから提供されています。また、`reactingTwoPhaseEulerFoam`をカスタマイズして凝縮モデルを組み込むことも一般的です。最大の利点は、ソースコードを直接編集して、独自の相変化速度モデル（例えば、界面での詳細な熱・物質平衡に基づくモデル）を実装できる点です。ただし、プリプロ・ポストプロや収束性のチューニングには商用ソフト以上の知識と労力が必要です。
凝縮熱伝達のトラブル対応
よくあるエラーと対策
計算中に界面が激しく振動したり、液膜がちぎれて飛び散る「数値的なスパイク」が発生します。なぜですか？
これは典型的な問題で、主に2つの原因が考えられます。1つは「メッシュが粗すぎて界面を解像できない」。もう1つは「相変化のソース項（質量移動速度）が大きすぎる」です。後者の場合、先述のLeeモデルの係数rが大きすぎると、一つのタイムステップで過剰な質量が移動し、圧力や体積率が暴走します。rを1/10にして計算を再開し、液膜の成長が物理的に妥当な速度（例えば、数ミリ秒で形成される）か確認してください。
壁面での熱伝達率が、教科書にあるNusselt理論値より異常に低く（または高く）出ます。考えられる原因は？
まず、境界条件を疑います。壁温度は一様で正しく設定されていますか？ 蒸気は飽和状態としていますか？ 過冷却はありませんか？ 次に、液膜の熱抵抗を正しく評価できているか確認。メッシュが粗く、液膜厚さを過小評価していると熱伝達率が高く出ます。逆に、壁面の濡れ性（接触角）の設定を誤り、液膜が思ったより厚く滞留していると、熱伝達率は低くなります。水平管の場合、重力方向や管周方向のメッシュも重要です。
非凝縮性ガス（空気など）が少量混入している場合、シミュレーション上でどう扱えばいいですか？
これは実務で非常に重要な問題です。非凝縮性ガスは蒸気の分圧を下げ、実質的な飽和温度を低下させます。さらに、界面付近にガスが蓄積（濃縮）し、蒸気が液面に拡散してくるのを妨げる「拡散抵抗」が生じ、凝縮速度が大幅に低下します。CFDでは、蒸気-非凝縮性ガスを「Species Transport」モデルで扱い、蒸気の質量分率に基づく局所的な飽和温度を計算し、それを凝縮モデルのT_satとして与える必要があります。この効果を無視すると、熱伝達率を数倍も過大評価する危険があります。
なった
詳しく
報告