熱交換器の伝熱計算の基礎知識

化学工学/伝熱入門伝熱入門の3本目(全3本)公開更新広告

熱交換器の伝熱計算の基礎として、流体の流れ方向と熱交換量の計算、固体壁内と壁付近の伝熱を解説します。

この記事の内容
  1. 2流体の流れ方向
  2. 熱交換量の計算
  3. 固体壁内の伝熱
  4. 固体壁付近の伝熱
  5. 熱交換器の伝熱計算
  6. 関連記事

2流体の流れ方向

熱交換器の高温側流体と低温側流体の流れ方向によって伝熱結果が大きく異なります。

並流

両流体が同じ方向に流れる並流では、低温側の出口温度が高温側の出口温度を超えないのが特徴です。

これは入口付近の温度差が大きく、出口付近の温度差が小さくなるためです。

加熱温度に上限がある場合に向いています。

向流

両流体が対向する方向に流れる向流では、低温側の出口温度が高温側の出口温度を超えられるのが特徴です。

これは低温側の出口が高温側の入口と同じ端にあり、低温側を高温側の入口温度近くまで温められるためです。

熱交換量を増やしたい場合に向いています。

並流と向流の温度分布。並流では両流体の温度が近づき、低温側の出口は高温側の出口を超えない。向流では低温側の出口が高温側の出口より高くなる。入口側と出口側の温度差ΔT1、ΔT2を示す
流れに沿った温度の変化(例:入口 高温100℃・低温20℃)。ΔT₁、ΔT₂ は LMTD の式の温度差

直交流

並流や向流以外にも、両流体が直交する方向に流れる直交流という考え方もあります。

気体の熱交換器でよく用いられます。

多管式熱交換器の流れ

多管式熱交換器は一般に、シェル内をチューブが往復する2パス形式がとられています。

そのためチューブ流体の往路は向流、復路は並流になります。

シェルの中をチューブが往復する2パスの熱交換器。シェル流体は右から左へ流れ、チューブ流体の往路(左から右)は向流、復路(右から左)は並流になる
2パス多管式の流れ。シェル流体は一方向、チューブ流体は往復する

更に内部にはバッフルが設けられています。

実際の流体の流れは非常に複雑です。

熱交換量の計算

まずは熱交換器において、どの程度の熱交換が行われたかの計算方法を解説します。

熱交換量

エネルギー保存則より、高温流体が失う熱量と低温流体が得る熱量は等しくなります。

ただし、これは十分な断熱により外界への放熱が無いこと、機器の不良で液漏れが無いことなど、外的要因を考慮していないことが前提です。

熱交換量

$$ Q=c_{pH}V_{H}\rho_{H}\left(T_{HI}-T_{HO}\right)=c_{pL}V_{L}\rho_{L}\left(T_{LO}-T_{LI}\right)$$

式は横にスクロールできます

Q:熱交換量[W]、cp:比熱[J/(kg・K)]、V:体積流量[m3/s]
ρ:密度[kg/m3]、T:温度[K]
H:高温流体、L:低温流体、I:熱交換器入口、O:熱交換器出口

質量流量の換算

先に紹介した計算式では体積流量を用いていますが、密度をかけて質量流量に換算しています。

これは体積流量は温度によって値が変化してしまうためです。

コリオリ流量計のように質量流量を直接測定できる場合は、以下の換算式を用いて質量流量で計算できます。

質量流量

$$ W=Vρ$$

式は横にスクロールできます

W:質量流量[kg/s]、V:体積流量[m3/s]、ρ:密度[kg/m3]

固体壁内の伝熱

高温側から低温側へ、固体壁を通って熱が伝わる場合を考えます。

その時の固体壁内の伝熱計算を解説します。

フーリエの法則

固体壁内における伝熱量はフーリエの法則から求められます。

フーリエの法則

$$ Q=kA\frac{T_{H}-T_{L}}{t}$$

式は横にスクロールできます

Q:伝熱量[W]、k:熱伝導度[W/(m・K)]、T:温度[K]、t:壁面厚さ[m]
A:伝熱面積[m2]、H:高温側の壁面、L:低温側の壁面

熱伝導度(熱伝導率)kは物質固有の値です。

温度によって熱伝導度は異なりますので、高温側と低温側の壁面温度の平均値から熱伝導度を決定します。

温度勾配

フーリエの法則における、壁面の温度差(TH-TL)と壁面厚さtの比は温度勾配と呼びます。

温度勾配が大きいほど単位時間当たりの伝熱量が増加します。

平板の壁の断面。高温側の面の温度TH、低温側の面の温度TL、壁の厚さt、伝熱量Qを示す。同じ温度差で壁の厚さを半分にすると温度勾配は2倍になり、伝熱量も2倍になる
壁の中の温度は直線で下がる。薄いほど勾配が急で Q は大きい

熱流束

ある断面に対して直交する方向に伝わる、単位時間・単位面積当たりの熱量を熱流束と呼びます。

固体壁面の伝熱量は熱流束を用いて以下のように表せます。

熱流束を用いた伝熱量

$$ Q=qA$$

式は横にスクロールできます

Q:伝熱量[W]、q:熱流束[W/m2]、A:伝熱面積[m2]

固体壁付近の伝熱

高温側から低温側へ、固体壁を通って熱が伝わる場合を考えます。

その時の固体壁付近の伝熱計算を解説します。

境膜

流体には粘性があるため、固体壁面付近では摩擦が発生して流速には勾配が現れます。

この勾配が存在する固体壁面付近の微小な領域が境膜です。

例えば多管式熱交換器においてはチューブの内側と外側に境膜が存在します。

境膜は流速だけでなく温度の変化も著しく、主に伝熱が起こるのは境膜です。

流速が小さいほど境膜は厚く、大きいほど境膜が薄い特徴があります。

流体と壁の境の温度。流体の本体は一様な温度THで、境膜の中で壁面温度TLまで下がる。流速が小さいと境膜が厚く勾配が緩やか、流速が大きいと境膜が薄く勾配が急になる
温度は境膜の中で落ちる(膜モデルの模式図)。境膜が薄いほど h は大きい

ニュートンの冷却法則

固体壁面と流体との伝熱量はニュートンの冷却法則から計算できます。

ニュートンの冷却法則

$$ Q=hA\left(T_{H}-T_{L}\right)$$

式は横にスクロールできます

Q:伝熱量[W]、h:境膜伝熱係数[W/(m2・K)]、A:伝熱面積[m2]
T:温度[K]、H:高温側、L:低温側

境膜伝熱係数

ニュートンの冷却法則で用いられる境膜伝熱係数hは物性値ではありません。

流体の物性や流速、流れの形式などに依存する値です。

ヌセルト数Nu、レイノルズ数Re、グラスホフ数Gr、プラントル数Prなど無次元数との関係式を用いて算出します。

熱交換器の伝熱計算

ここからは熱交換器における伝熱、つまり熱貫流における伝熱計算を解説します。

伝熱速度式

熱交換器においては高温流体から固体壁面、固体壁面内、固体壁面から低温流体と熱が伝わります。

しかし、これまで紹介したフーリエの法則やニュートンの冷却法則を用いるには、固体壁面の温度が無ければ計算できず現実的ではありません。

そのため実用的には伝熱速度式から算出します。

伝熱速度式

$$ Q=UA\Delta T$$

式は横にスクロールできます

Q:伝熱量[W]、U:総括伝熱係数[W/(m2・K)]、A:伝熱面積[m2]、ΔT:温度差[K]

総括伝熱係数

伝熱速度式では熱貫流のしやすさを表す総括伝熱係数Uが用いられています。

固体壁が平板の場合、この値は以下の式から算出されます(円管では内面と外面の面積の違いも考慮します)。

総括伝熱係数

$$ \frac{1}{U}=\frac{1}{h_{H}}+\frac{t}{k}+\frac{1}{h_{L}}$$

式は横にスクロールできます

U:総括伝熱係数[W/(m2・K)]、h:境膜伝熱係数[W/(m2・K)]
t:壁面厚さ[m]、k:熱伝導度[W/(m・K)]
H:高温側、L:低温側

高温流体、境膜、壁、境膜、低温流体を横切る温度の分布と、1/hH、t/k、1/hLを長さで表した熱抵抗の直列。各抵抗に比例して温度が落ち、和が1/Uになる
温度の落ちは各抵抗の大きさに比例する(例:hH=2000、hL=1000 W/(m²・K)、t=3 mm、k=16 W/(m・K)。U≒590)

ただし熱交換器を長期間使用すると、内部の汚れにより伝熱性能が低下します。

実際の総括伝熱係数は汚れの厚みと熱伝導度の比である汚れ係数rを高温側と低温側のそれぞれに加えて計算されます。

総括伝熱係数(汚れを補正)

$$ \frac{1}{U}=\frac{1}{h_{H}}+\frac{t}{k}+\frac{1}{h_{L}}+r_{H}+r_{L}$$

式は横にスクロールできます

U:総括伝熱係数[W/(m2・K)]、h:境膜伝熱係数[W/(m2・K)]
t:壁面厚さ[m]、k:熱伝導度[W/(m・K)]、r:汚れ係数[m2K/W]
H:高温側、L:低温側

円管の断面と側面。内半径ri、外半径ro、管の長さL。内面積Aiは2πriL、外面積Aoは2πroLで、外面のほうが広い。例では外径25mm、内径20mmで面積比1.25
円管の内面と外面(例:内径 20 mm、外径 25 mm)。面積の取り方で U の値が変わる

温度差に注意(対数平均温度差)

熱交換器内では、高温側と低温側の流体が熱交換しながら流れていきます。

そのため伝熱位置によって両流体の温度は変わり、温度差についても単純に入口温度と出口温度の差だけでは説明できません。

そこで、熱交換器の伝熱計算で用いる温度差には対数平均温度差LMTDを用います。

対数平均温度差LMTD

$$ \begin{align} LMTD&=\frac{\Delta T_{1}-\Delta T_{2}}{\ln \left(\Delta T_{1}/\Delta T_{2}\right)}\end{align}$$

式は横にスクロールできます

LMTD:対数平均温度差[K]
ΔT1:高温流体入口部における低温流体との温度差[K]
ΔT2:高温流体出口部における低温流体との温度差[K]

この式はリービッヒ冷却器のような単純な二重管構造において成り立ちます。

実際には補正係数FTをかけて補正したLMTDを用います。

FTはおおよそ0.8~1.0の値です。

様々な熱交換条件における補正係数FTは伝熱工学資料より参照できます。

伝熱工学資料(別サイト)日本機械学会 著 / 日本機械学会

関連記事

伝熱の種類(対流伝熱・伝導伝熱・放射伝熱)

潜熱と顕熱の違い(温度上昇と状態変化)

多管式熱交換器の構造と分類

伝熱入門の前後の記事

関連する計算ツール

化学工学計算支援ツール(ブラウザで使える無料の計算ツール)

関連する書籍

執筆・監修のご依頼は、ご相談のページで受け付けています。