2025年9月26日 星期五

熱傳導方程的三個「關鍵支柱」

對,你已經抓到熱傳導方程的三個「關鍵支柱」。我幫你整理成一個清晰的邏輯鏈條:


🔑 熱傳導 PDE 的三件事

1. 能量守恆 (Conservation of Energy)

  • 控制體積內能量的變化率 = 邊界熱流 + 內部熱源。

  • 數學表達式:

    tVcρu(x,t)dV=SqndS+VQ(x,t)dV\frac{\partial}{\partial t} \iiint_V c \rho \, u(\vec{x},t) \, dV = - \iint_S \vec{q} \cdot \vec{n} \, dS + \iiint_V Q(\vec{x},t) \, dV

👉 這是 宏觀的能量平衡方程,是出發點。


2. Flux 與溫度的關係 (Fourier’s Law)

  • 熱通量(heat flux)與溫度梯度成正比:

    q=ku\vec{q} = -k \nabla u

    其中:

    • kk:材料的熱傳導係數(thermal conductivity)

    • u\nabla u:溫度場的梯度

👉 這一步把「材料性質」帶進來了。


3. 散度定理 (Divergence Theorem)

  • 把邊界熱流(表面積分)換成體積積分:

    SqndS=VqdV\iint_S \vec{q} \cdot \vec{n} \, dS = \iiint_V \nabla \cdot \vec{q} \, dV

👉 這讓我們可以把所有項目都轉換成「體積積分」形式。


🚀 導出 PDE

代入能量守恆:

cρut=q+Qc\rho \frac{\partial u}{\partial t} = - \nabla \cdot \vec{q} + Q

再代入 Fourier’s Law (q=ku\vec{q} = -k \nabla u):

cρut=k2u+Qc\rho \frac{\partial u}{\partial t} = k \nabla^2 u + Q

最後,定義熱擴散係數 α2=kcρ\alpha^2 = \tfrac{k}{c\rho}

ut=α22u+Qcρ\frac{\partial u}{\partial t} = \alpha^2 \nabla^2 u + \frac{Q}{c\rho}

  • 無熱源 (Q=0Q=0):

    ut=α22uu_t = \alpha^2 \nabla^2 u

這就是 熱傳導 PDE(Heat Equation)。

沒有留言: