流體的能量守恆方程式推導(Energy Conservation Equation)

在國高中物理裡面，相信大家對能量守恆方程式並不陌生，不過此次要帶大家來看的是用來求解流體問題的能量守恆偏微分方程。

Skiptocontent首頁關於科技解密產品設計CAE法律主廚私藏聯絡Menu首頁關於科技解密產品設計CAE法律主廚私藏聯絡前面推導連續方程式與Navier-Stokes方程式後，下一個要來推導的是能量守恆方程式。

在國高中物理裡面，相信大家對能量守恆方程式並不陌生，不過此次要帶大家來看的是用來求解流體問題的能量守恆偏微分方程。

能量守恆方程式推導圖1控制體積的能量進出示意圖流體內能與動能考慮一個沿著流跡線移動的控制體積(ControlVolume)，如圖1所示。

流體進出控制體積的能量分別包含流體本身內能與動能，基於能量守恆概念，進出控制體積的流體能量總和等於進來能量扣除出去能量。

因此，想要求得流體進出能量總和的話，需要知道如何表示進來能量與出去能量，如式(1)–式(4)。

 E_x=\rhou(e+\frac{1}{2}V\cdotV)dy -(1)E_y=\rhov(e+\frac{1}{2}V\cdotV)dx-(2) E_{x,dx}=\left(\rhou(e+\frac{1}{2}V^2)+\frac{\partial(\rhou(e+\frac{1}{2}V^2))}{\partialx}dx\right)dy -(3)E_{y,dy}=\left(\rhov(e+\frac{1}{2}V^2)+\frac{\partial(\rhov(e+\frac{1}{2}V^2))}{\partialx}dy\right)dx -(4)其中 E_x與 E_y分別為X方向與Y方向進來能量， E_{x,dx}與 E_{y,dy}為X方向與Y方向出去能量，ρ為密度，e為流體單位質量內能，V為速度向量，dx與dy分別是控制體積長度。

藉由式(1)–式(4)，控制體積內部的能量總和可表示成式(5)。

  dE=-\left(\frac{\partial(\rhou(e+\frac{1}{2}V^2))}{\partialx}+\frac{\partial(\rhov(e+\frac{1}{2}V^2))}{\partialy}\right)dxdy-式(5)流體熱傳除了流體本身內能與動能外，流體帶來的熱能也是一個需要考量的因素，而流體熱傳關係式可用傅立葉定律表示，如式(6)。

q=-\lambda\frac{\partialT}{\partialx} -(6)其中q為熱通量，λ為熱傳導係數，T為溫度。

基於傅立葉定律與能量守恆，可以仿造式(5)的形式列出控制體積內部的熱能總和，如式(7)。

dQ=\left(\frac{\partial}{\partialx}\left(\lambda\frac{\partialT}{\partialx}\right)+\frac{\partial}{\partialy}\left(\lambda\frac{\partialT}{\partialy}\right)\right)dxdy-(7)外力做功在Navier-Stokes方程式推導與物理意義一文中曾提到沿著流跡線移動的控制體積會受到壓力、正應力與剪應力作用，這些外力除了會帶來動量變化以外，也會對控制體積做功，X方向做功總和與Y方向做功總和可分別表示成式(8)跟式(9)。

dW_x=\left(-\frac{\partialP\cdotu}{\partialx}+\frac{\partial\sigma_{xx}\cdotu}{\partialx}+\frac{\partial\tau_{xy}\cdotv}{\partialx}\right)dxdy-(8)dW_y=\left(-\frac{\partialP\cdotv}{\partialy}+\frac{\partial\sigma_{yy}\cdotv}{\partialy}+\frac{\partial\tau_{yx}\cdotu}{\partialy}\right)dxdy-(9)能量守恆前面分別探討流體內能與動能、熱能、外力做功這幾種能量來源，而能量守恆即為控制體積內部的能量變化等於上述幾種能量變化總和，如式(10)所示。

\frac{\partial\left[\rho(e+\frac{1}{2}V^2)\right]}{\partialt}=dE+dQ+dW_x+dW_y+K\cdotV+\rho\cdotq_s-(10)式(10)後兩項分別是其他外力做功與熱能，像是上述未討論到的電磁力作功與燃燒熱或輻射熱等等。

將式(10)化簡後可得能量守恆的偏微分方程，如式(11)所示。

\rho\left(\frac{\partiale}{\partialt}+u\frac{\partiale}{\partialx}+v\frac{\partiale}{\partialy}\right)=\left(\frac{\partial}{\partialx}\left(\lambda\frac{\partialT}{\partialx}\right)+\frac{\partial}{\partialy}\left(\lambda\frac{\partialT}{\partialy}\right)\right)-p\cdot\triangledownV+K\cdotV+\rho\cdotq_s+\mu\cdot\phi-(11)式(11)末項為逸散項，該項只有在超音速下產生極大速度梯度時不能忽略，對於一般的不可壓縮流體而言，該項可忽略不計。

此外，由於可透過連續方程式與Navier-Stoke方程式聯合求解，流體速度與壓力為已知數，故式(11)的未知數剩下流體能量與溫度，由於一個方程式只能求解一個未知數，因此將流體能量用溫度來表示，如式(12)。

e=c_p\cdotT-\frac{P}{\rho}-(12)其中c_p為定壓比熱。

接著將式(12)代入式(11)可得未知數為溫度的流體能量守恆方程式，如式(13)所示。

\rhoC_p\left(\frac{\partialT}{\partialt}+u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}\right)=\left(\frac{\partial}{\partialx}\left(\lambda\frac{\partialT}{\partialx}\right)+\frac{\partial}{\partialy}\left(\lambda\frac{\partialT}{\partialy}\right)\right)+\left(\frac{\partialP}{\partialt}+u\frac{\partialP}{\partialx}+v\frac{\partialP}{\partialy}\right)+K\cdotV+\rho\cdotq_s+\mu\cdot\phi-(13)能量守恆方程式經常因應不同物理狀況來求解，而對於不可壓縮流體(incompressibleflow)而言，式(13)後面四項可用Q^T表示，搭配數學符號就能將式(13)表示成更精簡的式(14)，即為此次要跟大家介紹的流體能量守恆方程式。

\frac{\partial}{\partialt}\left(\rhoC_pT\right)+\triangledown\cdot\left(\rhoC_pVT\right)=\triangledown\cdot\left(\lambda\triangledownT\right)+Q^T-(14)主廚結語本次跟大家介紹流體能量守恆方程式如何推導。

雖然花了一段時間，不過總算介紹完計算流體力學的守恆方程式，而在了解三大守恆方程式之後，後續要來看怎麼在電腦上求解這幾個方程式！如果想持續查看最新消息，歡迎追蹤科技雞湯Facebook或Instagram！參考資料JousefMurad,DerivationoftheEnergyEquation,https://www.youtube.com/watch?v=in68iBENe14F.Moukalledetal.,TheFiniteVolumeMethodinComputationalFluidDynamics,2015相關文章:偏微分方程分類對流動特性的探討[CFD]流體的連續方程式(ContinuityEquation)推導文章導覽←Previous文章Next文章→LeaveaCommentCancelReply發佈留言必須填寫的電子郵件地址不會公開。

必填欄位標示為*Typehere..Name*Email*Website在瀏覽器中儲存顯示名稱、電子郵件地址及個人網站網址，以供下次發佈留言時使用。

主廚小叮嚀文章下方留言需經人工審核才會顯示！若尚未顯示的話，請等等喔！若覺得本站雞湯不錯，可點選下方按鈕至臉書按讚或追蹤IG！把好喝的雞湯拿來跟好朋友一起分享吧！FacebookInstagramEnvelope最新雞湯什麼是邊界條件？有限體積法的邊界條件離散什麼是結構色？蝴蝶翅膀的奈米科技！流體拓樸優化概念與原理什麼是渦振現象？塔科馬海峽吊橋崩塌事件！什麼是NFT？帶你看懂NFT爆紅的背後因素！有限體積法離散原理：傳輸方程式雞湯分類科技解密CAE主廚私藏產品設計法律