為什麼剛倒的熱咖啡，攪一攪就涼得比較快？

為什麼剛倒的熱咖啡，攪一攪就涼得比較快？

你一定有過這個經驗：剛沖好的咖啡燙得無法入口，於是你拿起湯匙不停攪拌，沒幾分鐘它就降到適口的溫度。但如果只是靜靜放著，可能得等上十幾分鐘。同樣一杯咖啡、同樣的室溫，為什麼「攪一攪」會差這麼多？

答案藏在熱傳學（Heat Transfer）的三種機制裡。熱量從杯壁傳到空氣靠的是對流（convection），攪拌讓杯口表面不斷換上溫度較低的液體，等效於提高了對流速率；同時，熱也透過杯子本身的固體傳導（conduction）散到桌面，並以輻射（radiation）的形式向四周發射紅外線。這三條路徑同時運作，決定了任何物體冷卻或加熱的快慢。

熱傳學研究的核心問題只有一個：熱量以多快的速率、沿著什麼路徑、從高溫處流向低溫處？ 熱力學（thermodynamics）告訴我們「能量守恆、熵會增加、熱往低溫流」，但它不談「要花多久」。熱傳學補上了「速率」這塊拼圖，是熱機、散熱片、建築保溫、晶片冷卻乃至太空船熱控制的共同語言。

第一種機制：傳導，熱在固體裡的接力賽

傳導發生在物質內部：高溫區域的分子振動劇烈，透過碰撞把動能一個傳一個地交給鄰居，能量於是「擴散」過去，物質本身並不移動。金屬之所以摸起來特別「冰」，不是因為它真的比較冷，而是因為它把你手上的熱導走得太快。

描述傳導的核心是傅立葉定律（Fourier's law）。一維穩態下：

$$ q_x = -k A \frac{dT}{dx} $$

其中 $q_x$ 是熱流率（W），$k$ 是材料的熱導率（thermal conductivity，W/m·K），$A$ 是截面積，$\frac{dT}{dx}$ 是溫度梯度。負號代表熱往溫度下降的方向流——這其實是熱力學第二定律在微觀層次的具體表現。

對於一面厚度 $L$、兩側溫度分別為 $T_1$、$T_2$ 的平板，穩態熱流是：

$$ q = \frac{k A (T_1 - T_2)}{L} $$

工程上常把這個式子改寫成「電路類比」的形式。我們定義熱阻（thermal resistance）：

$$ R_{\text{cond}} = \frac{L}{kA}, \qquad q = \frac{T_1 - T_2}{R_{\text{cond}}} $$

這與電學的歐姆定律 $I = \Delta V / R$ 完全同構：溫差像電壓、熱流像電流、$R$ 像電阻。多層牆壁（如建築的磚牆＋保溫層＋石膏板）就像串聯電阻，總熱阻直接相加：

$$ R_{\text{total}} = \sum_i \frac{L_i}{k_i A} $$

不同材料的 $k$ 差距驚人：銅約 400 W/m·K、不鏽鋼約 15、玻璃約 1、木材約 0.15、靜止空氣約 0.026。這也是為什麼保溫材料（羽絨、發泡塑膠、雙層玻璃）的設計精髓，其實是「把空氣關起來不讓它流動」——只要空氣不對流，它就是極佳的絕熱體。

更一般的情況，熱傳導由熱傳導方程（heat equation）描述：

$$ \frac{\partial T}{\partial t} = \alpha \nabla^2 T $$

其中 $\alpha = \dfrac{k}{\rho c_p}$ 是熱擴散率（thermal diffusivity，m²/s），$\rho$ 是密度、$c_p$ 是比熱。這條偏微分方程描述溫度場如何隨時間演化，它與物理裡的擴散方程同形，數學上也與量子力學的薛丁格方程有親緣關係。

第二種機制：對流，熱搭著流體的順風車

對流是熱量隨著流體（液體或氣體）的整體流動被「運走」。它不再只是分子間的傳遞，而是一團團帶著能量的流體在移動。對流分兩類：自然對流（natural convection）靠的是溫差造成的密度差（熱空氣上升、冷空氣下沉），強制對流（forced convection）則靠風扇、幫浦等外力。攪拌咖啡、電腦裝風扇，都是把自然對流升級為強制對流。

對流的工程描述是牛頓冷卻定律（Newton's law of cooling）：

$$ q = h A (T_s - T_\infty) $$

$h$ 是對流換熱係數（convective heat transfer coefficient，W/m²·K），$T_s$ 是表面溫度，$T_\infty$ 是遠處流體溫度。$h$ 不是材料常數，而是綜合了流體性質、流速、幾何形狀的「行為參數」。典型量級：自然對流空氣 $h \approx 5\text{–}25$、強制對流空氣 $h \approx 25\text{–}250$、沸騰的水可高達 $h \approx 2500\text{–}100000$。這就解釋了為什麼 100°C 的蒸氣燙傷遠比 100°C 的空氣嚴重——蒸氣凝結時的 $h$ 大上好幾個數量級。

對流同樣可寫成熱阻形式：

$$ R_{\text{conv}} = \frac{1}{hA} $$

於是「固體導熱 + 表面對流」可串聯成一條完整的熱路徑。

要預測 $h$，工程師用無因次數（dimensionless numbers）。最關鍵的是努塞爾數（Nusselt number）：

$$ \mathrm{Nu} = \frac{hL}{k_f} $$

它代表「對流相對於純傳導強了多少倍」。$\mathrm{Nu}$ 通常表達成另外兩個數的函數——雷諾數（Reynolds number） $\mathrm{Re} = \dfrac{\rho u L}{\mu}$（慣性對黏性之比，判斷層流或紊流）與普朗特數（Prandtl number） $\mathrm{Pr} = \dfrac{\mu c_p}{k} = \dfrac{\nu}{\alpha}$（動量擴散對熱擴散之比）。例如管內紊流的經典關係式 Dittus–Boelter：

$$ \mathrm{Nu} = 0.023\, \mathrm{Re}^{0.8} \mathrm{Pr}^{n} $$

（加熱時 $n=0.4$，冷卻時 $n=0.3$）。這種「用無因次數整理實驗數據」的策略，是熱流領域把複雜現象壓縮成可用公式的核心方法論。

第三種機制：輻射，不需要介質的熱

前兩種機制都需要物質當媒介，但輻射不用。所有溫度高於絕對零度的物體，都會以電磁波的形式向外發射能量；太陽的熱穿過真空的太空抵達地球，靠的就是輻射。室溫物體主要輻射紅外線，溫度越高、波長越短，鐵燒到發紅、燈絲發白，都是輻射隨溫度往可見光偏移的結果。

輻射的基石是史蒂芬–波茲曼定律（Stefan–Boltzmann law）。理想黑體（blackbody）的輻射功率密度為：

$$ E_b = \sigma T^4 $$

$\sigma = 5.67 \times 10^{-8}\ \mathrm{W/m^2K^4}$ 是史蒂芬–波茲曼常數。請特別注意這裡是 $T^4$ 而且 $T$ 必須用絕對溫度（K）。真實物體不是完美黑體，需乘上輻射率（emissivity） $\varepsilon \in [0,1]$：

$$ E = \varepsilon \sigma T^4 $$

兩個表面之間的淨輻射交換，常用的簡化式為：

$$ q_{\text{rad}} = \varepsilon \sigma A (T_s^4 - T_{\text{surr}}^4) $$

$T^4$ 這個非線性是輻射最棘手也最迷人之處：溫度翻倍，輻射功率變成 16 倍。因此在高溫場合（火爐、燃氣輪機、太空船），輻射往往主宰一切；而在室溫附近，輻射雖然存在但常被對流掩蓋。值得一提的是，史蒂芬–波茲曼定律可由量子物理的普朗克黑體輻射定律（Planck's law）對全波長積分推導而來——熱傳學的這條工程公式，根源其實在量子力學。

把三者合起來：總熱傳與熱路徑

真實問題幾乎都是三種機制並存。以一面外牆為例：室內空氣對流加熱牆內側 → 熱傳導穿過牆體 → 牆外側對流＋輻射散到戶外。工程師把它組成串聯熱阻，定義總熱傳係數（overall heat transfer coefficient） $U$：

$$ q = U A \,\Delta T, \qquad \frac{1}{UA} = \sum R_i = \frac{1}{h_i A} + \frac{L}{kA} + \frac{1}{h_o A} $$

這個 $U$ 值就是建築能源法規裡規範窗戶、牆體保溫性能的指標。冷氣機的散熱片、汽車的水箱、發電廠的冷凝器，本質上都是在「用最小體積換到最大的 $UA$」。

看一個例子

一片矽晶片運作時產生 $q = 15\ \mathrm{W}$ 的廢熱，必須透過頂部的鋁散熱片排掉。已知散熱片底座與晶片接觸面積 $A = 4\ \mathrm{cm^2} = 4\times10^{-4}\ \mathrm{m^2}$，底座是厚 $L = 5\ \mathrm{mm}$ 的鋁（$k = 200\ \mathrm{W/m\cdot K}$），散熱片表面以強制對流散熱，$h = 80\ \mathrm{W/m^2\cdot K}$，散熱片有效散熱面積 $A_{\text{fin}} = 0.05\ \mathrm{m^2}$，環境空氣 $T_\infty = 30°\mathrm{C}$。問晶片表面溫度約為多少？

步驟一：傳導熱阻。

$$ R_{\text{cond}} = \frac{L}{kA} = \frac{0.005}{200 \times 4\times10^{-4}} = 0.0625\ \mathrm{K/W} $$

步驟二：對流熱阻。（用散熱片的有效面積）

$$ R_{\text{conv}} = \frac{1}{h A_{\text{fin}}} = \frac{1}{80 \times 0.05} = 0.25\ \mathrm{K/W} $$

步驟三：總熱阻與溫升。 兩者串聯：

$$ R_{\text{total}} = 0.0625 + 0.25 = 0.3125\ \mathrm{K/W} $$

$$ \Delta T = q\, R_{\text{total}} = 15 \times 0.3125 = 4.69\ \mathrm{K} $$

$$ T_{\text{chip}} = T_\infty + \Delta T = 30 + 4.69 \approx 34.7°\mathrm{C} $$

晶片表面約 $35°\mathrm{C}$，遠低於危險的 $85°\mathrm{C}$ 上限，散熱設計合格。從計算也看得出：對流熱阻是傳導的 4 倍，是整條路徑的瓶頸——這告訴我們，與其把鋁底座做得更薄，不如加大風扇或散熱片面積，效益更高。這正是「熱阻串聯，最大的那一項決定一切」的工程直覺。

重點回顧

• 三種機制各有公式：傳導用傅立葉定律 $q=-kA\,dT/dx$、對流用牛頓冷卻定律 $q=hA\Delta T$、輻射用史蒂芬–波茲曼定律 $q=\varepsilon\sigma A T^4$。
• 熱阻電路類比是工程主力工具：$R_{\text{cond}}=L/kA$、$R_{\text{conv}}=1/hA$，串聯相加，最大熱阻是散熱瓶頸。
• $k$ 是材料常數，$h$ 不是：$h$ 取決於流速、幾何與流體，需透過 $\mathrm{Nu}=f(\mathrm{Re},\mathrm{Pr})$ 的經驗式推估。
• 輻射對溫度極度敏感（$\propto T^4$，須用絕對溫度），高溫時主宰、室溫時常被對流掩蓋。
• 熱力學談「能不能、往哪走」，熱傳學談「有多快」，兩者互補而不可替代。

深入探討（研究所視角）

研究所層級的熱傳學會在幾個方向深化。第一，對流的邊界層理論。 牛頓冷卻定律的 $h$ 並非憑空而來，而是源自流體在固體表面形成的速度邊界層（velocity boundary layer）與熱邊界層（thermal boundary layer）。事實上，$h$ 由表面處的溫度梯度決定：

$$ h = \frac{-k_f \left.\dfrac{\partial T}{\partial y}\right|_{y=0}}{T_s - T_\infty} $$

對流問題因而本質上是流體力學（Navier–Stokes 方程）與能量方程的耦合求解。普朗特數 $\mathrm{Pr}$ 正是兩個邊界層厚度比的量度：$\mathrm{Pr}\gg1$（如油）熱邊界層薄、$\mathrm{Pr}\ll1$（如液態金屬）熱邊界層厚。

第二，暫態與非穩態傳導。 真實系統常在升溫降溫的過渡中。當物體內部溫度可視為均勻（集總容法 lumped capacitance，判據為畢歐數 $\mathrm{Bi}=hL_c/k < 0.1$），冷卻服從指數律：

$$ \frac{T(t)-T_\infty}{T_i-T_\infty} = \exp\!\left(-\frac{hA}{\rho V c_p}\,t\right) $$

時間常數 $\tau = \rho V c_p / (hA)$ 決定快慢。若 $\mathrm{Bi}$ 不小，就得回到熱傳導偏微分方程，用分離變數、Heisler 圖或數值方法（有限差分、有限元素 FEM、計算流體力學 CFD）求解溫度場。

第三，相變與沸騰傳熱。 沸騰、凝結、熔化涉及潛熱，$h$ 隨熱通量呈現高度非線性的「沸騰曲線」，存在臨界熱通量（CHF）這類安全極限——核反應爐與高功率電子冷卻都繞不開它。

第四，輻射的光譜與方向性。 工程簡化常假設灰體與漫射，但真實輻射率隨波長（光譜選擇性表面，如太陽能集熱器的選擇性塗層）與角度變化，且多表面之間需引入視角因子（view factor）做輻射網路分析。

最後，把熱傳學放回大圖像：它與統計力學（$k$、$c_p$ 的微觀起源）、量子物理（普朗克輻射、聲子導熱）、流體力學（對流的 Navier–Stokes 耦合）共享根基。在奈米尺度，當特徵長度逼近聲子平均自由程，傅立葉定律甚至會失效，需用波茲曼輸運方程（BTE）描述——這正是當代微奈米熱傳（micro/nano-scale heat transfer）的前沿，與半導體散熱、熱電材料息息相關。