對於很多數學和工程問題,我們常常需要使用到梯度、散度和旋度方程,而有的時候,在使用這些方程時,我們卻對它們其中的數學、物理意義不甚清楚,結果便是看着很多在此基礎上建立的公式而一頭霧水。這篇文章便從這三大方程的本質入手,推導它們在三大經典座標系下的形式,揭露其”廬山真面目“!
旋度公式
旋度的理解
旋度,單從字面上不難看出,它是個描述旋轉劇烈程度的量,大自然中有很多旋轉的現象,例如水的漩渦、地球的自轉等,因此我們需要有個來描述這些旋轉現象劇烈程度的"工具",於是,旋度就被先輩們提出來了。那麼,接下來的問題是,我們該如何定義某個物理量(必須是向量,理由見下文)的旋度呢?這便是我們接下來需要討論的問題。
我相信,在閱讀本篇文章的讀者應該對環量不會陌生(如果您現在不清楚環量的概念,或者已經忘記環量的概念,筆者強烈建議您先去學習環量的概念後再閱讀此文),所謂環量,直白地講就是對空間中一條封閉的曲線做第二類曲線積分,環量代表着某一物理量在該曲線上的集度,因爲是第二類曲線積分,所以很明顯被積的物理量得是向量。舉個簡單的例子,你晨跑時繞着操場跑了一圈,這裏,操場便能被看成一條封閉曲線,而你繞着跑一圈的過程中,速度在操場上的投影繞操場的積分就是速度環量。這個繞操場跑的事件,我們可以看成是繞操場某一中心點旋轉的過程,想象一下,如果有兩個人在同一操場上的同一條跑道上用時不同地跑完了一圈,誰的速度環量大呢?答案很明顯,那個用時短的人速度環量大,因爲積分的路徑不變,速度提高了,速度積分自然也就增大了。我們也可以想象到跑的快的人繞操場的"旋轉"越快。從這個例子出發,我們很自然地能類推出在其它一些物理現象中,如果積分路徑相同,環量大的自然旋轉越劇烈。那麼,如果我們需要描述空間中的某個點所在區域上的旋轉劇烈程度呢?很簡單,點可以看成是個面積無窮小的曲面,我們只需要把積分曲線無限縮小,最終趨於一個點,我們把這時獲得的物理向量的環量除以該封閉積分曲線所包圍曲面的面積(無窮小),定義爲某一物理向量在該點處的環量密度,也就是旋度,旋度的大小就可以代表出該點所在區域的旋轉劇烈程度。
圖1.繞操場跑圈
笛卡爾座標系下的旋度公式
我們現在已經知道了旋度是咋計算的,接下來就來推導一下它在三大座標系下的公式吧。老規矩,先來推導最經典的笛卡爾座標系下的旋度公式。
圖2.笛卡爾座標系
上圖中,在D點處有個ω的角速度,我們知道,角速度也是個矢量,所以我們可以將它分解到x、y、z軸上,我們先來看繞x軸的旋轉。假設該空間中有個物理向量A,其繞DCGH的環量我們分開來看,DC上,我們發現物理向量A隨着y的改變也是變化的,但因爲我們現在已經把DCGH縮小到無限小,於是我們可以認爲物理向量A沿着DC是單調變化的,我們可以得到以下不等式:
Ay(x,y,z)dy⩽∫DCAydy⩽Ay(x,y+dy,z)dy
根據夾逼準則,因爲dy趨於無窮小,所以最終得
∫DCAydy=Ay(D點處的值)dy=Ay(x,y,z)dy
同理,我們來算CG上(僅z變化)的積分:
∫CGAzdz=Az(C點處的值)dz=Az(x,y+dy,z)dz
GH上(僅y變化)的積分:
∫GHAydy=−Ay(H點處的值)dy=−Ay(x,y,z+dz)dy
HD上(僅z變化)的積分:
∫HDAzdz=−Az(D點處的值)dz=−Az(x,y,z)dz
所以,A在DCGH上的積分爲:
∮DCHGA⋅dl=Ay(x,y,z)dy+Az(x,y+dy,z)dz−Ay(x,y,z+dz)dy−Az(x,y,z)dz
並且易得DCHG的面積爲dydz,我們將式6除以這個面積即可得到繞x軸旋轉引起的旋度:
dydz∮DCHGA⋅dl=dydAz−dzdAy
同理我們可以得到繞y軸旋轉引起的旋度:
dzdx∮DHEAA⋅dl=dzdAx−dxdAz
繞z軸旋轉引起的旋度:
dxdy∮DABCA⋅dl=dxdAy−dydAx
將式7,8,9相加起來,我們即得到D點處的旋度:
curlA=(dydAz−dzdAy)i+(dzdAx−dxdAz)j+(dxdAy−dydAx)k
柱面座標系下的旋度公式
圖3. 柱面座標系
柱面座標系下旋度公式的推導過程和笛卡爾座標系下的旋度公式相同。我們同樣來考察D點處的旋度,將角速度分解到eρ、eϕ、ez方向上,這裏我們以推導繞ρ軸旋轉引起的旋度爲例:
沿DC上的積分:
∫DCAϕρdϕ=Aϕ(D點處的值)ρdϕ=Aϕ(ρ,ϕ,z)ρdϕ
沿CG上的積分:
∫CGAzdz=Az(C點處的值)dz=Az(ρ,ϕ+dϕ,z)dz
沿GH上的積分:
∫GHAϕρdϕ=−Aϕ(H點處的值)ρdϕ=−Aϕ(ρ,ϕ,z+dz)ρdϕ
沿HD上的積分:
∫HDAzdz=−Az(D點處的值)dz=−Az(ρ,ϕ,z)dz
又因爲DCHG圍成的面積等於ρdϕdz,於是得繞ρ軸旋轉引起的旋度爲:
ρdϕdz∮DCHGA⋅dl=ρdϕdAz−dzdAϕ
同理,繞ϕ軸引起的旋度爲:
dzdρ∮DHEAA⋅dl=dzdAρ−dρdAz
繞z軸引起的旋度爲:
ρdρdϕ∮DABCA⋅dl=ρdρd(ρAϕ)−ρdϕdAz
最終我們得到柱面座標系下的D點處的旋度公式:
curlA=(ρdϕdAz−dzdAϕ)eρ+(dzdAρ−dρdAz)eϕ+(ρdρd(ρAϕ)−ρdϕdAz)ez
球面座標系下的旋度公式
有了笛卡爾座標系和柱面座標系的推導經驗,我相信你現在可以推導球面座標系下的旋度公式了,拿出紙筆動動手吧,如果在推導過程中有啥問題,可以在評論區留言哦!這裏我們直接給出球面座標系下的旋度公式,推導完成後你可以校對一下哦😀:
圖4. 球面座標系
球面座標系下的旋度公式爲:
感謝閱讀!
與諸君共前行!
參考文獻
- Del in cylindrical and spherical coordinates - 維基百科(https://en.wikipedia.org/wiki/Del_in_cylindrical_and_spherical_coordinates)
- 高等數學(西工大版本)