垂直平板の膜状凝縮シミュレーター 戻る
熱工学

垂直平板の膜状凝縮シミュレーター

冷たい垂直壁に飽和蒸気が凝縮し、薄い液膜となって重力で流れ落ちる「膜状凝縮」をNusselt理論で解くツールです。流体・飽和温度・壁面温度・平板の高さを変えると、平均熱伝達率・膜厚・熱流束・凝縮量がリアルタイムで分かります。

パラメータ設定
凝縮する流体
液の密度・熱伝導率・粘性・潜熱を自動設定
飽和温度 T_sat
凝縮する蒸気の飽和温度
壁面温度 T_w
冷却された平板の表面温度
平板の高さ L
m
蒸気が凝縮する平板の鉛直方向の長さ
計算結果
平均熱伝達率 h̄ (W/m²K)
底部の膜厚 δ (mm)
熱流束 q (kW/m²)
凝縮量 (kg/m²·h)
凝縮膜レイノルズ数 Re
過冷度 ΔT (K)
凝縮膜のアニメーション

蒸気が冷たい壁に凝縮し、液膜となって重力で流れ落ちます。膜厚 δ(x) は上から下へ x^(1/4) で厚くなり、その液膜が蒸気→壁の熱伝達の支配抵抗になります。

膜厚 δ vs 高さ x
平均熱伝達率 h̄ vs 過冷度 ΔT
理論・主要公式

$$\bar h=0.943\left[\frac{\rho_l(\rho_l-\rho_v)g\,h_{fg}\,k_l^3}{\mu_l\,L\,\Delta T}\right]^{1/4}$$

垂直平板の層流膜状凝縮における平均熱伝達率 h̄(Nusselt理論)。ρl・ρv:液・蒸気の密度、g:重力加速度、hfg:蒸発潜熱、kl:液の熱伝導率、μl:液の粘性係数、L:平板高さ、ΔT=Tsat−Tw:過冷度。h̄ は ΔT と L のそれぞれ −1/4 乗に比例する。

$$\delta(L)=\left[\frac{4\mu_l k_l L\,\Delta T}{\rho_l(\rho_l-\rho_v)g\,h_{fg}}\right]^{1/4}$$

底部(x=L)における凝縮膜の厚さ δ。膜厚は高さ x の 1/4 乗に比例して上から下へ厚くなる。

$$q=\bar h\,\Delta T,\qquad Re=\frac{4\,\dot m\,L}{\mu_l}$$

熱流束 q と凝縮膜レイノルズ数 Re(ṁ:単位面積あたりの凝縮質量流束 q/hfg)。重力で流れ落ちる液膜が支配抵抗であり、Nusselt理論は層流の重力排水膜を前提とする。Re が約1800を超えると膜は乱流となる。

垂直平板の膜状凝縮とは

🙋
「膜状凝縮」って、お風呂の冷たい鏡が湯気で曇って、水の筋がツーッと流れるあれと同じことですか?
🎓
まさにそれだよ。飽和した蒸気が、自分の飽和温度より冷たい壁に触れると凝縮して液になる。その液が壁全体を「膜」のように覆って、重力で下へ流れ落ちる——これが膜状凝縮だ。発電プラントの復水器やエアコンの凝縮器など、蒸気を液に戻すあらゆる熱交換器で起きている現象なんだ。垂直平板はその一番基本的なモデルだよ。
🙋
なるほど。でも、蒸気が壁にくっついて熱を渡すだけなら、すごく熱がよく伝わりそうですよね。なんでわざわざ計算するんですか?
🎓
いい質問だ。ポイントは「壁に直接触れているのは蒸気じゃなくて液膜」ということ。蒸気から壁へ熱が届くには、その液膜を熱伝導で通り抜けなきゃいけない。液膜は薄いけど、熱伝導率はそんなに高くない。だから液膜そのものが一番の熱抵抗になるんだ。ざっくり言うと h̄ ≈ kl/δ。膜が薄いほど熱がよく伝わる。左で壁面温度 T_w を飽和温度に近づけて過冷度を小さくすると、膜が薄くなって h̄ が上がるのが分かるよ。
🙋
え、じゃあ過冷度を大きくして「もっと冷やす」と、かえって熱伝達率は下がっちゃうんですか?
🎓
そう、ここが面白いところ。過冷度ΔTを大きくすると確かに「熱流束 q」は増えるんだけど、たくさん凝縮するぶん液膜が厚くなる。だから単位温度差あたりの伝わりやすさ、つまり h̄ 自体は ΔT の −1/4 乗で下がっていくんだ。右の「h̄ vs 過冷度」のグラフを見ると、ゆるやかに右肩下がりのカーブになっているだろう? 平板の高さ L も同じで、L の −1/4 乗で h̄ が下がる。背の高い壁ほど下のほうに液が積み重なって膜が厚くなるからだよ。
🙋
膜厚がそんなに大事なんですね。膜の厚さは上から下まで同じじゃないんですか?
🎓
違うんだ。一番上では凝縮が始まったばかりで膜はほぼゼロ。下に行くほど、上から流れてきた液が次々と合流して膜が厚くなる。Nusselt理論ではこの膜厚 δ(x) が高さ x の 1/4 乗に比例する。だから上のほうが膜が薄くて熱がよく伝わり、下のほうは膜が厚くて伝わりにくい。アニメーションでも、壁の上端は薄く、下端でぷっくり厚くなっているのが見えるはずだよ。
🙋
「凝縮膜レイノルズ数」っていうのも出てきますけど、これは何を見ているんですか?
🎓
液膜の流れが層流か乱流かを見る指標だよ。流れ落ちる液膜も「流れ」だから、レイノルズ数で評価できる。Re が約1800を超えると液膜が乱流に遷移する。乱流になると膜の中で混合が起きて熱がよく伝わるようになるから、層流前提のNusselt理論は実際より熱伝達率を低く見積もってしまう。背の高い凝縮器や過冷度の大きい運転では乱流膜になりやすいから、そのときは乱流補正を入れた相関式に切り替える必要があるんだ。

よくある質問

垂直平板上の層流膜状凝縮では、Nusselt理論から平均熱伝達率を h̄ = 0.943·[ρl(ρl−ρv)g·hfg·kl³ /(μl·L·ΔT)]^(1/4) で求めます。ρl・ρvは液・蒸気の密度、gは重力加速度、hfgは蒸発潜熱、klは液の熱伝導率、μlは液の粘性係数、Lは平板の高さ、ΔT=Tsat−Twは過冷度です。本ツールはこの h̄ と膜厚 δ、熱流束 q、凝縮量を同時に計算します。
壁面で凝縮した液は重力で流れ落ちながら平板上に薄い液膜をつくります。蒸気から壁へ熱が伝わるには、必ずこの液膜を熱伝導で通過しなければなりません。液膜は薄いとはいえ熱伝導率が小さいため、液膜そのものが主たる熱抵抗になります。h̄ ≈ kl/δ のイメージで、膜が薄いほど熱伝達率は高くなります。だから過冷度を小さく、平板を短くして膜を薄く保つ設計が有利です。
Nusselt理論では h̄ は ΔT の −1/4 乗、平板高さ L の −1/4 乗に比例します。過冷度ΔTが大きいほど凝縮量が増えて液膜が厚くなり、L が長いほど上から流れ落ちた液が積み重なって下部の膜が厚くなります。膜が厚くなれば熱抵抗が増えるため、h̄ は下がります。膜厚 δ(x) は高さ x の 1/4 乗に比例して上から下へ厚くなります。
凝縮膜レイノルズ数 Re = 4·ṁ·L/μl は液膜の流れの状態を表します。Re が約1800を超えると液膜が乱流に遷移します。乱流膜では混合により熱伝達が促進されるため、層流を前提とするNusselt理論は実際より低い熱伝達率を予測します(過小評価)。背の高い凝縮器や過冷度が大きい場合は乱流膜となりやすく、乱流補正を加えた相関式が必要です。

実世界での応用

発電プラントの復水器:火力・原子力発電のタービンを出た蒸気は、復水器で冷却水に熱を渡して凝縮し、水に戻ってボイラへ循環します。復水器の伝熱管表面でまさに膜状凝縮が起きており、管表面に張る液膜の厚さが熱交換性能を左右します。プラント全体の熱効率に直結するため、膜状凝縮の熱伝達率を正しく見積もることは発電設計の要です。

空調・冷凍機の凝縮器:エアコンやヒートポンプの凝縮器では、圧縮された冷媒蒸気(R134aなど)が空冷・水冷の管表面で凝縮します。冷媒は水蒸気より潜熱が小さく液の物性も異なるため、同じ過冷度でも熱伝達率や凝縮量は大きく変わります。本ツールで流体を切り替えると、その違いを直感的に確認できます。

化学プラントの蒸留・蒸発設備:蒸留塔のリボイラやコンデンサ、多重効用蒸発缶などでは、溶剤蒸気やアンモニア蒸気が伝熱面で凝縮します。アンモニアは潜熱が大きく熱伝達率が高い反面、過冷度設定や材料選定に注意が必要です。プロセス設計では膜状凝縮の理論値を出発点に、汚れ係数や蒸気流速の影響を加えて設計します。

CAE・伝熱解析の検証基準:凝縮を含む二相流のCFD解析や熱交換器シミュレーションでは、結果が物理的に妥当かを確かめる基準が必要です。Nusselt理論は解析解が存在する数少ない凝縮問題であり、垂直平板の単純条件ならツール値と解析値を直接比較できます。CFD結果がこの理論値と桁違いなら、相変化モデルや境界条件の設定ミスを疑うサニティチェックになります。

よくある誤解と注意点

まず最大の落とし穴が、「膜状凝縮と滴状凝縮を混同する」ことです。本ツールが扱うのは、液が壁を一様な膜で覆う膜状凝縮です。これに対し、壁が液をはじく(撥水性の)表面では、凝縮液が膜にならず無数の小さな水滴となって転がり落ちる「滴状凝縮」が起こります。滴状凝縮は壁面の大部分が液膜で覆われないため、熱伝達率が膜状凝縮の5〜10倍にもなります。Nusselt理論は膜状凝縮専用であり、滴状凝縮には適用できません。実機の凝縮面は時間とともに濡れていき、たいてい膜状凝縮に落ち着くため、保守的には膜状凝縮で設計します。

次に、「過冷度を大きくすれば凝縮性能が上がる」と思い込むこと。過冷度ΔTを大きくすると、確かに熱流束 q = h̄·ΔT は増えます。しかし熱伝達率 h̄ 自体は ΔT の −1/4 乗で低下します。つまり「冷やせば冷やすほど効率よく熱が伝わる」わけではありません。さらに過冷度や平板高さを大きくしすぎると凝縮膜レイノルズ数が1800を超えて乱流膜に遷移し、層流前提のNusselt理論はもはや成立しません。本ツールの判定が「乱流」を示したら、層流の式の値は過小評価であると理解し、乱流補正式に切り替える必要があります。

最後に、「Nusselt理論の値をそのまま実機性能だと考える」こと。Nusselt理論は、層流・重力排水・液膜の慣性無視・蒸気のせん断応力なし・一定物性、といった理想化の上に成り立っています。実際の凝縮器では、蒸気流速が速いと液膜が引きずられて薄くなり熱伝達が向上したり、逆に膜が波立って挙動が変わったりします。また伝熱面に不凝縮ガス(空気など)が混入すると、それが壁近くに溜まって凝縮を著しく阻害します。本ツールの値はあくまで理想化された出発点であり、実機設計では蒸気流速・不凝縮ガス・汚れ・面の傾きの補正を必ず加えてください。

使い方ガイド

  1. 飽和温度Tsat(℃)を50~150℃の範囲で設定し、凝縮対象の蒸気種(水蒸気)の物性値を決定します
  2. 平板壁面温度Tw(℃)をTsatより5~30K低く設定して過冷度ΔTを確保し、凝縮駆動力を与えます
  3. 垂直平板の高さL(mm)を100~2000mmで入力し、Nusselt理論に基づく膜状凝縮の支配方程式を実行します
  4. 出力の平均熱伝達率h̄、底部膜厚δ、熱流束q、凝縮量(kg/m²·h)、膜レイノルズ数Reを確認します

具体的な計算例

飽和温度100℃、壁温70℃、平板高さ1.5mの水蒸気凝縮を想定します。物性値:密度ρL=960kg/m³、潜熱hfg=2257kJ/kg、動粘度νL=3.5×10⁻⁴m²/s、熱伝導率λL=0.68W/m·Kを用いると、Nusselt式より平均熱伝達率h̄≈11,200W/m²K、底部膜厚δ≈0.18mm、熱流束q≈252kW/m²、凝縮量≈404kg/m²·hが得られます。膜レイノルズ数Re≈142となり層流領域を確認できます。

実務での注意点

  1. Nusselt理論は層流膜状凝縮(Re<1800)を前提とするため、高流量条件では波状流への遷移を検討し、h̄の相関式切り替えが必要です
  2. 蒸気の過冷却度ΔTが大きいほどh̄は増加しますが、凝縮器設計ではエネルギー収支と伝熱面積のトレードオフを評価してください
  3. 垂直平板の高さ増加に伴い膜厚δが増大し、h̄が低下するため、複数段または複数ステーションの凝縮器配置を検討します
  4. 実機の熱交換器では非凝縮性ガス(空気など)の混在、液滴エントレインメント、表面粗さが現象を複雑にするため、安全係数0.7~0.85の適用を推奨します