微分方程 · 熱方程式

熱方程式

熱方程式uₜ = uₓₓ表示:每一點溫度變化的快慢,由該處溫度曲線的彎曲程度決定。在桿上,它會抹平所有起伏,越細的起伏越早消失。

從圖中看熱方程式

上方圖形顯示從0到π\pi的桿上溫度uu:淡色是初始狀態,黃色是時刻tt。藍色箭頭表示各點升溫或降溫的快慢。下方的桿依溫度著色,最下面一排以對數刻度顯示各正弦模式的大小(兩端絕熱時則是餘弦模式)。

隨時間的變化率,就是沿桿方向的二階導數
ut=uxxu_t = u_{xx}

為什麼曲線彎得越厲害的地方,溫度變化越快?取小步長hh,有uxx(x)≈(u(x−h)+u(x+h)−2u(x))/h2u_{xx}(x) \approx \bigl(u(x-h) + u(x+h) - 2u(x)\bigr)/h^2:也就是兩個鄰點的平均值與u(x)u(x)之差,再乘上2/h22/h^2。當uxx>0u_{xx} > 0,曲線向上彎,該點低於附近兩個鄰點的平均,因此升溫;當uxx<0u_{xx} < 0,則降溫。在直線段上,無論多熱都有ut=0u_t = 0。探針可以用附近兩個鄰點與其平均值展示這個關係。跳躍與尖角會立刻變平滑:只要t>0t > 0,溫度就是平滑曲線。

為什麼細小的起伏先消失?單一模式(兩端維持0時的sin⁡kx\sin kx,或兩端絕熱時的cos⁡kx\cos kx)滿足uxx=−k2uu_{xx} = -k^2 u,所以形狀不變,只按因子e−k2te^{-k^2 t}整體縮小。起伏數量變成3倍,同樣高度下的二階導數大小就變成9倍,衰減也快9倍。任何初始溫度都可寫成這些模式的和,所以最下面的長條從右側開始下降:模式kk在tk=ln⁡10/k2t_k = \ln 10/k^2之後,剩不到初始大小的10%。在對數時鐘下,時間每增加至10倍,被抹平的模式就寬約10≈3.2\sqrt{10} \approx 3.2倍。紅尺2t2\sqrt t表示這個尺度:單一模式的半波長不超過尺長時,振幅最多剩e−π2/4≈8.48%e^{-\pi^2/4} \approx 8.48\%。

兩端會改變什麼?維持0的兩端相當於冰,uu以冰的溫度為基準。隨時間推移,初始狀態中存在的最低頻模式主導形狀。若初始溫度處處非負、且不全為零,熱從兩端流出,∫u dx\int u\,dx下降,曲線變成按e−te^{-t}降低的單一正弦拱。絕熱兩端不讓熱通過,因此∫u dx\int u\,dx不變。此時模式為餘弦,k=0k = 0模式就是平均溫度,永遠不衰減:桿最後各處都趨近平均溫度。兩種邊界條件都從t=0t = 0起作用;若初始溫度在端點不符合條件,例如兩塊冰之間溫度均勻的桿,端點就會立即調整至符合條件。

值得留意

  • 所有箭頭共用一個比例,由當下變化最快的點決定;圖例標出參考箭頭代表的變化率。箭頭長度只能在同一時刻比較,不能跨時刻比較。
  • 比冰冷的溫度為負值。兩端維持0時,若初始溫度有正有負,就顯示淨值∫u dx\int u\,dx,而非比例。

操作方式

  • 按「擴散」,時間每增加至10倍,就留下一條淡色快照。
  • 按一下桿,即可將探針移到該處。
  • 選擇「手繪」,在上方圖形拖曳以設定初始溫度;繪圖會讓tt回到0。

相關內容