流量・温度・総括伝熱係数Uを変えながら熱交換量・LMTD・F補正・伝熱面積・NTU・有効度をリアルタイム計算。温度プロファイルをグラフ表示します。
熱交換器の性能を評価するための二つの基本的なアプローチ（LMTD法とNTU法）を支える中心的な式です。
伝熱量Qと必要な伝熱面積Aを求める基本設計式です。
$$ Q = \dot{m}_h C_{ph}(T_{hi}- T_{ho}) = \dot{m}_c C_{pc}(T_{co}- T_{ci}) $$ $$ A = \frac{Q}{U \cdot F \cdot LMTD}$$LMTD（対数平均温度差）: $LMTD = \frac{\Delta T_1 - \Delta T_2}{\ln(\Delta T_1 / \Delta T_2)}$
ここで、対向流の場合 $\Delta T_1 = T_{hi}- T_{co}, \Delta T_2 = T_{ho}- T_{ci}$
Q: 伝熱量 [W], ṁ: 質量流量 [kg/s], Cp: 比熱 [J/kg･K]
U: 総括伝熱係数 [W/m²･K], A: 伝熱面積 [m²], F: 補正係数
熱交換器の形状とサイズが決まっている時、与えられた入口条件から出口温度や伝熱量を予測する方法です。
$$ NTU = \frac{U A}{C_{min}}, \quad \varepsilon = \frac{Q}{Q_{max}}= \frac{C_h (T_{hi}-T_{ho})}{C_{min}(T_{hi}-T_{ci})}= \frac{C_c (T_{co}-T_{ci})}{C_{min}(T_{hi}-T_{ci})}$$NTU（伝熱単位数）: 熱交換器の性能の高さを表す無次元数。値が大きいほど高性能。
ε（有効度）: 実際の伝熱量と、理論上可能な最大伝熱量との比。0から1の値を取る。
C: 水当量 $C = \dot{m} C_p$ [W/K], Cmin: 高温側と低温側の水当量の小さい方
化学プラント・石油精製：プロセス流体の加熱・冷却に不可欠です。例えば、蒸留塔から出てくる高温の留分で、塔へ送入する原料を予熱する「フィード・エフェレント交換器」として多用されます。省エネルギー化の鍵となる設備です。
発電所（火力・原子力）：蒸気タービンで仕事をした後の低温低圧の蒸気を復水器で冷却水（海水や河川水）で凝縮し、再びボイラーに戻す循環システムの中核です。巨大なシェルアンドチューブ熱交換器が使われます。
空調・冷凍システム：ビルの空調システムでは、冷水と冷却水の熱交換に使われます。また、アンモニアやフロンなどの冷媒を凝縮する「凝縮器」や、蒸発する「蒸発器」としても基本的な構造が応用されています。
船舶・車両エンジン：大型船舶のディーゼルエンジンでは、潤滑油を冷却する「オールクーラー」や、過給機で圧縮されて高温になった吸入空気を冷却する「インタークーラー」としてコンパクトなシェルアンドチューブ熱交換器が採用されています。
このシミュレーターを使い始める際、特に初心者がハマりがちな落とし穴がいくつかあります。まず第一に、「総括伝熱係数Uは定数」と思い込んでしまうこと。実はU値は流速や温度、流体の汚れ（スケーリング）によって大きく変動します。例えば、冷却水側の流速を2倍にすると、U値はおよそ2の0.8乗倍（約1.74倍）に増えることが多く、必要な伝熱面積Aが一気に小さくなります。ツールでUを固定して設計しても、実際の運転条件でUが想定より低ければ、性能不足に陥るので要注意です。
次に、「温度差が大きければ何でも良い」という誤解。確かにLMTDが大きいほど小さな面積で済みますが、高温側と低温側の入口温度差（例えばThiとTci）を極端に大きくしすぎると、今度は材料の熱応力やコストの問題が発生します。また、低温側出口温度Tcoを高温側入口温度Thiに極限まで近づけようとすると（「温度接近」）、LMTDが劇的に小さくなり、必要面積が無限大に近づくという現象が起こります。これはツールでTcoを95℃などに設定してみると、面積Aが急激に増大するので体感できるでしょう。
最後に、「補正係数Fは後で考えればいい」という軽視。Fはシェル側・チューブ側のパス数や温度効率P、熱容量比Rから決まりますが、ここで安易な値を選ぶと大変なことになります。例えば、1シェルパス-2チューブパスでPが0.7を超えるような条件では、Fが0.8を切ってしまうことがあります。これは「理想の70%しか温度差が活用できていない」ことを意味し、設計面積が30%以上も過大になる可能性があります。まずはツールでパス数を変えながらFがどう動くか、感覚を掴んでください。
シェルアンドチューブ熱交換器の設計計算は、単なる熱力学の応用ではなく、複数の工学分野が交差する総合技術です。まず深く関わるのが「流体力学」、特に圧力損失の計算です。伝熱面積を大きくするためにチューブ数を増やしたり流速を上げると、必ず圧力損失ΔPが増加します。ポンプやコンプレッサーの動力はΔPに比例するので、伝熱設計と流体設計のトレードオフを考える「熱流体設計」が必須となります。
次に「材料工学」との関連です。総括伝熱係数Uを高めるために伝熱管を薄くするほど、強度や腐食に対する余裕が減ります。また、高温部や低温部で異なる材料を使う場合、熱膨張の差による応力（熱応力）を「構造力学」的に評価する必要があります。特にシェル側で蒸気凝縮のような相変化を伴う場合は、「二相流」というさらに複雑な分野の知見が必要になります。
さらに現代では「システム制御」との連携も重要です。部分負荷運転時や流入温度が変動する条件下で、目標とする出口温度を維持するには、バイパス流路を制御したり、流量を調整したりします。その制御系を設計する際の基礎となるのが、このツールで計算するNTU（伝熱単位数）です。NTUが大きい熱交換器は応答が早く、制御性が良い傾向があります。このように、一つの熱交換器を軸に、伝熱・流動・材料・制御の知識が一本の線で繋がっていくのです。
このシミュレーターでLMTD法とNTU法の関係に慣れたら、次のステップとして「微視的な視点」を加える学習をお勧めします。具体的には、総括伝熱係数Uを構成する個々の熱抵抗を計算してみることです。Uは次の式で表されます：
$$ \frac{1}{U A} = \frac{1}{h_{h} A_{h}} + R_{fouling, h} + \frac{\ln(r_o/r_i)}{2\pi k L} + R_{fouling, c} + \frac{1}{h_{c} A_{c}} $$
ここで$h$は熱伝達率、$R_{fouling}$は汚れ係数、$k$は伝熱管の熱伝導率です。次の課題は、チューブ内の流速から$h$を推定する「相関式」（例えばディタス・ベルトルの式）を自分で計算し、それがUや最終的な面積Aにどう効くかを追ってみることです。これが「詳細設計」の第一歩です。
数学的には、LMTDがなぜ対数平均になるのかを理解するために、熱交換器の微小区間での熱収支を微分方程式でモデル化し、それを積分するプロセスを一度たどってみると、理解が深まります。これは「分布定数系」の考え方の入り口でもあります。
最後に、シェルアンドチューブ以外の熱交換器の形態（プレート式、フィンチチューブ式、スパイラル式）について調べ、なぜその形状が選ばれるのかを比較してみましょう。例えば、高粘度流体にはどういう形式が向くのか、コンパクトさが求められる場合はどうか。これにより、単なる計算ツールの操作から、「最適な熱交換器の選定」という実務的な視点を養うことができます。