『如果 σ 代表向量函數, ∇σ 就包含有向量與純量兩部分, 它們可以寫成 S∇σ 和 V ∇σ, 我建議純量部分為 σ 的散度 (divergence), 這是向量函數效應的好名字。 但 是 ∇σ, 一般說來也是有向量的部分, 可以稱這種向量為原向量函數的旋度 (curl ) 或變體 (version)1, 它代表了由向量σ 所承載之內容旋轉的量與方向。』
— James Clerk Maxwell (1831∼1879) —
一、 前言
正如微積分是為了解決牛頓力學而發明的, 向量分析則是應運電磁學而產生的。 學向量分 析心中常有電磁學讓你永遠不孤單也不會迷失。 我願意再強調一次
『數學沒有物理是瞎子, 物理沒有數學是跛子。』
一個數學系特別是應用數學系把 『普通物理』 列為選修 (等價於不用修), 等於告訴世人它是一 個爛系。 歷史告訴我們的還有時代的潮流就是跨領域! 一個只懂數學 (本質上不是真的懂), 而 對其它領域無知的人是絕對沒有競爭力。 就算是數學專業的人也必須同時對於非主修的科目有 所涉獵, 如此才可以培養寬廣的胸襟與開放的心靈, 對於別人的研究也有興趣並懂得欣賞。 否則 會越來越孤立以至於枯竭, 片面內省的數學玄想只會導致貧瘠。 百年前德國哥廷根偉大數學傳 統的建立者 Felix Klein (1849∼1925) 就已經有這樣的呼籲:
『我們強烈感受到, 在現代科學思想的快速發展下, 我們的學科正瀕臨孤立的危險。 自 從現代分析興起以來, 讓數學與理論自然科學雙方皆獲益的親密互動關係, 已經面臨 崩潰的威脅。 · · · 今日數學界的當務之急, 是將其純科學部門和其具有最重要應用價 值的物理科學部門再次緊密結合, 這種關係以往所得的豐碩成果可以從 Lagrange 和 Gauss 的研究裡看到。』
1version 是翻譯、 譯本、 敘述或異見, 但這個字應分解為versio-n-其拉丁文是vertere 而vers- 的意思是轉動 (to turn)。 所以它也 翻為a special form or variant of something (變體)。
30
圖解梯度、 散度與旋度 31
他的吶喊現在看來仍然是先知的警語, 令人深思。 看來數學的發展與文明或民主一樣, 不見得都 是微分大於 0(增函數), 永遠是進步的! 如果你不守護它, 深化其內涵並隨時進步, 一旦獨裁的 敗家子出現, 整個數學社會被迫走入歧途, 則幾代人的努力馬上就煙消雲散。
1865 年左右, 馬克士威 (James Clerk Maxwell, 1831∼1879) 綜合電磁現象的安培定 律及法拉第定律, 建立電磁場的數學理論, 他提出了變化的磁場激發電場的概念之後, 發現還缺 少某些重要的東西。 假如對這些方程式進行數學變換, 那麼就得出一種與電磁學方程組體系的 其它公式矛盾的關係。 假如這個方程式藉著一個附加項, 即所謂的位移電流來補充, 那麼這個矛 盾就可以克服。 所以他提出位移電流來克服電磁學方程組之矛盾, 而且整理成完整的電磁場理 論, 因此馬克士威就從法拉第的解釋者, 成為獨創的學者。 馬克士威方程包含了完整的電學、 磁 學和光學的理論, 這是馬克士威對物理學所作的一項偉大的統一, 將電、 磁、 光三者合而為一。
但是馬克士威的成果在當時 (1865) 並沒有立即引起數學界的注意與興趣, 所以並不像牛頓與 Leibniz (1646∼1716) 那個年代承襲古希臘的信念:
『大自然是根據數學設計的。』
他們將數學看成是通往大自然的真理之路。 這個信仰吸引了 Bernoulli 家族、 Euler (1707
∼1783)、 Lagrange (1736∼1813)、 Laplace (1749∼1827)、 · · · 、 一直到 Fourier (1768
∼1830), 這些同時具有科學家身份的大數學家為之獻身, 進而開創了 17∼18 世紀數學創造最 偉大的年代。 更重要的是除了科學之外他們也將理性的精神注入到人文、 社會領域。 如果 19 世紀後半葉的數學家們能夠像 Euler 等人一樣將牛頓的工作視為研究的核心, 也將 Maxwell 方程當為最根本重要的問題來關心, 則他們極有可能早就發現愛因斯坦的狹義相對論 (special relativity)、 拓撲群及其表現理論還有雙曲偏微分程的大部分理論以及泛函分析。 但是數學界 卻失去了這大好的機會與 Maxwell 方程攜手開創數學與科學史的第二個黃金年代。 事實上那 時候數學界並不缺乏有見識的人, 例如, 牛津的數學家 Henry Smith (1826∼1883) 就已經體 會到需要新的數學以了解 Maxwell 方程, 並在 1873 年左右呼籲數學家研究馬克士威的著作。
但 Smith 已經 46 歲而且他的學養不足以成為新領域的開創者, 有能力的年輕人並不聽從他 的建議。 在 20 世紀完成 Smith 先知呼聲的 Hermann Minkowski (1864∼1909) 與 Jaco- ques Hadamard (1865∼1963) 那時分別是 9 歲與 8 歲, Elie Cartan (1869∼1951) 只有 3 歲, 而 Hermann Weyl (1885∼1955), Jean Leray (1906∼1998) 與 Harish Chandra (1923∼1983) 則還沒有出生。 我想最主要的原因是19世紀後半葉數學家的興趣是複變函數論、
解析數論與不變量理論, 而這些豐富的成果已經形成數學界美學與價值觀的標準, 以至於形式 上複雜的 Maxwell 方程被拒絕了。
雖然我們可以說馬克士威是向量分析的創始人, 但他那個時候還沒有向量的符號, 以至於 方程式必須以各個分量寫出所以非常的繁複, 這使得他偉大的創見並沒有馬上為當時的人所認 可。 如果沒有專門且精準的符號, 則現代數學根本不可能產生。 向量分析最重要的進展主要是歸
Heaviside ab V ab ∇ div curl
梯度、 散度與旋度形式上來看都是一次微分, 但本質上卻有極大的差異。 要瞭解它們就必須從 Gauss 散度定理與 Stokes 定理著手, 雖然這兩個定理的內涵都是微積分基本定理, 但是一個 關心的是法向量 (散度) 另一個則是切向量 (旋度)。 在二維平面給定切向量 (a, b) 則朝外法向 量為 (b, −a), 反過來亦然, 知道其中一個向量立即可得另一個向量。 所以二維的 Green 定理同 時具有 Gauss 散度定理與 Stokes 定理兩個身分。 讀者如果會怕三維空間, 那麼就先把 Green 定理弄清楚, 但不能只學一般微積分課本教你的切向量之情形。 這就是主動與被動學習的差異, 你 (妳) 是要當自己生命的主人還是等著讓別人來決定你 (妳) 的未來! 這兩個積分定理的主要 貢獻者為
(1) G. Green (1793∼1841) : 1828 年提出 Green 定理, Green 等式。
(2) 高斯 (Karl Friedrich Gauss, 1777∼1855) : 提出散度定理。
(3) M. Ostogradsky (1801∼1861) : 1831 年研究熱傳導方程 (heat equation) 時也導出散 度定理。
(4) George Gabriel Stokes (1819∼1903) : 1850 年給 Lord Kelvin (1824∼1907) 的信中 提到 Stokes 定理, Stokes 1854 年在劍橋的數學競試以這定理為題目, 其中得獎的學生有 J. C. Maxwell。 由此可見學生 Maxwell 遠比老師 Stokes 厲害的多。
本文主要的目的是想要以直觀、 圖解的角度來介紹曲線座標系統的梯度、 散度與旋度這三 個向量分析最重要的概念。 大部分數學系只活在一維空間, 把向量分析視為是微積分的推廣。 這 是極大的誤解, 同時也帶來錯誤的數學觀。 通常學生對於非垂直座標的梯度、 散度與旋度都非常 之陌生, 頂多就是把它當成公式, 需要的時候再查表。 對此我們並不滿意, 公式並不僅僅是公式, 它會講話, 因此我們特別介紹量綱 (因次) 分析 (dimensional analysis) 這個簡單但重要的觀 念、 直觀且有感覺地認識所有的方程式與 (不) 等式。 讀者也應嘗試能從公式中解釋各項之意 義, 並從中體會公式 (或方程式) 本身之物理或幾何意義。 想更深入瞭解梯度、 散度與旋度, 我 們需要微分形式 (differential form) 這個微分幾何的基本知識, 我們將在下一篇文章介紹。 關
圖解梯度、 散度與旋度 33
於量綱 (或因次) 分析, 請讀者參閱我在數學傳播的文章[6]。 這個觀念是我最有心得也最引以 為豪的方法, 它幫助我直觀地了解方程式, 讓我所學的數學變得有血有淚, 不再是冰冷無情的怪 咖。
二、 曲線座標系統
回顧一下直角座標的散度 (divergence)、 旋度 (curl) 與梯度 (gradient), 這三個算子正 好可以對應到向量的三個重要運算 : 內積 (inner product)、 外積 (cross product) 與直積 (direct product) (已知 F = (F1, F2, F3) 是向量函數、 f 是純量函數)
內積 : 向量 7→ 純量 div F = ∇ · F = ∂F1
∂x + ∂F2
∂y + ∂F3
∂z
外積 : 向量 7→ 向量 curl F = ∇ × F =
i j k
∂
∂x
∂
∂y
∂
∂z
F1 F2 F3
直積 : 純量 7→ 向量 gradf = ∇f = ∂f
∂x,∂f
∂y,∂f
∂z
(2.1)
直角座標之弧長元素與體積元素
ds2 = dr · dr = dx2+ dy2+ dz2, dV = dxdydz (2.2) 還有 i, j, k 三個方向之面積元素為
dSx = dydz, x = c 之平面 dSy = dzdx, y = c 之平面 dSz = dxdy, z = c 之平面
(2.3)
在彎曲空間梯度、 散度與旋度, 這三個算子與弧長元素、 體積元素與面積元素要如何重新定義 呢? 這是我們所關心的問題。
空間上任一位置向量(position vector) 表示為
r= (x, y, z) = xi + yj + zk (2.4) 其中 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) 就是直角座標的標準基底, 而且滿足右手法則 (i → j → k)
i· i = j · j = k · k = 1
i× j = k, j× k = i, k× i = j
dr = ∂r
∂xdx + ∂r
∂ydy + ∂r
∂zdz = dxi + dyj + dzk (2.6) 換句話說, {i, j, k} 或 {ex, ey, ez} 可由位置向量r之一階偏導數來表示。 全微分並不受座標之 限制, (2.6) 可以換成任意的座標, 也因此提供我們推廣至彎曲空間的思考方式。
考慮一般的曲線座標系統 (curvilinear system)。 已知 x1, x2, x3 是 u1, u2, u3 的函數 xi = xi(u1, u2, u3), ui = ui(x1, x2, x3), i = 1, 2, 3 (2.7) 我們要求 (x1, x2, x3) 7−→ (u1, u2, u3) 這個變換是一對一, 即 Jacobian
∂(x1, x2, x3)
∂(u1, u2, u3) =
∂x1
∂u1
∂x1
∂u2
∂x1
∂u3
∂x2
∂u1
∂x2
∂u2
∂x2
∂u3
∂x3
∂u1
∂x3
∂u2
∂x3
∂u3
= ∂r
∂u1 · ∂r
∂u2 × ∂r
∂u3
6= 0 (2.8)
向量 {∂u∂r
1,∂u∂r
2,∂u∂r
3} 正是 {u1, u2, u3} 三個曲線之切向量。 條件 (2.8) 保證這組向量是線性 獨立 (linear independent), 仿 (2.5) 我們可以選取這組來做為基底, 再化為單位向量 (這是 學數學的好習慣)
e1= 1 h1
∂r
∂u1
, e2 = 1 h2
∂r
∂u2
, e3 = 1 h3
∂r
∂u3
(2.9) h1=
∂r
∂u1
, h2 =
∂r
∂u2
, h3 =
∂r
∂u3
(2.10)
因此 {e1, e2, e3} 是 u1, u2, u3-座標曲線之單位向量, 就是我們所要的彎曲空間之基底, 同時 也是 u1, u2, u3-座標曲線的單位切向量。 若 {e1, e2, e3} 互相垂直則稱曲線座標系為正交曲線 座標系。 因為我們的目的是詮釋梯度、 散度與旋度的物理意義, 為避免節外生枝, 在這篇文章我 們只討論正交曲線座標系。
由全微分得 dr = ∂r
∂u1
du1+ ∂r
∂u2
du2+ ∂r
∂u3
du3 = h1du1e1+ h2du2e2+ h3du3e3 (2.11)
圖解梯度、 散度與旋度 35
圖1: 正交曲線座標系統
值得一提的是dui永遠要配合hi, 因為hidui才是長度, 其量綱為[hidui] = L, 這個事實我們一 直會用到。 由畢氏定理容易看出來弧長元素
ds2 = dr · dr = h21du21+ h22du22+ h23du23 (2.12) 體積元素與各個面的表面積則為 (請參考圖 1)
dV = dx1dx2dx3 = ∂(x1, x2, x3)
∂(u1, u2, u3)du1du2du3 = h1h2h3du1du2du3 (2.13) dS1= h2h3du2du3, dS2 = h3h1du3du1, dS3 = h1h2du1du2 (2.14) 讀者很容易可以驗證
[dV ] = L3, [dSi] = L2, i = 1, 2, 3
例題 2.1 (圓柱座標). 空間上任意一點 (x, y, z) 可以表示為圓柱座標 (r, θ, z)
r= (x, y, z) = (r cos θ, r sin θ, z) (2.15) h1 = hr =
∂r
∂ρ
= 1, h2 = hθ =
∂r
∂θ
= r, h3 = hz =
∂r
∂z
= 1 (2.16) 所以
er= 1 h1
∂r
∂r = (cos θ, sin θ, 0) = cos θi + sin θj
圖2: 圓柱座標
eθ= 1 h2
∂r
∂θ = (− sin θ, cos θ, 0) = − sin θi + cos θj (2.17) ez= 1
h3
∂r
∂z = (0, 0, 1) = k
容易驗證 {er, eθ, ez} 是彼此互相垂直, 而且滿足右手法則
er⊥ eθ, eθ ⊥ ez, ez ⊥ er
er× eθ = ez, eθ× ez = er, ez× er= eθ
(2.18)
為何其右手法則之順序是 r → θ → z 呢? 可以這麼看 : 從原點沿著半徑到 r, 而後逆時針轉 θ, 最後再垂直爬升到 z (參考圖 2 圖形路徑 1 → 2 → 3), 是不是正好就是右手法則 !
位置向量之全微分為
dr = drer+ rdθeθ+ dzez, [dr] = [rdθ] = [dz] = L (2.19) 而圓柱座標之弧長元素、 體積元素與各個面的表面積則為 (請參考圖 3)
ds2 = dr2+ r2dθ2+ dz2, dV = rdrdθdz dSr= rdθdz, dSθ = drdz, dSz = rdrdθ
(2.20)
這裡沒有必要透過全微分與複雜的連鎖律來得到弧長元素 ds2, 因為 {er, eθ, ez} 是一正交 系統, 無窮小方塊 (蛋糕形狀) 的長寬高為 dr, rdθ, dz, 所以由畢氏定理容易看出來弧長元素 (2.20)1, 一旦長寬高決定了則體積元素與各個面的表面積可由圖形直觀地看出來。
例題 2.2 (球座標). 空間上任意一點 (x, y, z) 可以表示為球座標 (ρ, ϕ, θ)
r= (x, y, z) = (ρ sin ϕ cos θ, ρ sin ϕ sin θ, ρ cos ϕ) (2.21)
圖解梯度、 散度與旋度 37
圖3: 圓柱座標各個面
圖4: 球座標之右手法則
h1 = hρ=
∂r
∂ρ
= 1, h2 = hϕ =
∂r
∂ϕ
= ρ, h3 = hθ =
∂r
∂θ
= ρ sin ϕ (2.22) 因此
eρ= 1 h1
∂r
∂ρ = (sin ϕ cos θ, sin ϕ sin θ, cos ϕ) eϕ = 1
h2
∂r
∂ϕ = (cos ϕ cos θ, cos ϕ sin θ, − sin ϕ) eθ = 1
h3
∂r
∂θ = (− sin θ, cos θ, 0)
(2.23)
實際上{eρ, eϕ, eθ} 形成一右手正交座標系統。 為何其右手法則之順序是 ρ → ϕ → θ 呢 ? 可 以這麼看 : 從原點順著半徑到 ρ, 而後垂直往下轉高低角 ϕ (緯度!), 最後再水平逆時針轉水平 角 θ (經度!), (請參考圖 4 圖形路徑 1 → 2 → 3), 是不是正好就是右手法則! 位置向量之全 微分為
dr = dρeρ+ ρdϕeϕ+ ρ sin ϕdθeθ (2.24)
圖5: 球座標
因為角度與三角函數不具有量綱, [dϕ] = [sin ϕ] = [dθ] = 1, 所以 [dρ] = [ρdϕ] = [ρ sin ϕdθ] = L
要特別提醒讀者的是第三個邊長是 ρ sin ϕdθ 並不是 ρdθ, 因為球面上在相同高低角(相同緯 度)兩點之弧長投影至赤道這個平面時其對應的圓的半徑等於 ρ sin ϕ, 所以這兩點的距離為 ρ sin ϕdθ, sin ϕ 這個量顯示的是投影的效應。 對於球座標我們仍然先掌握無窮小方塊(像西瓜 切片但要切除尖的部分之形狀) 的三邊長等於 dρ, ρdϕ, ρ sin ϕdθ, 所以球座標之弧長元素、
體積元素與各個面的表面積為
ds2 = dρ2+ ρ2dϕ2+ ρ2sin2ϕdθ2, dV = ρ2sin ϕdρdϕdθ dSρ= ρ2sin ϕdϕdθ, dSϕ = ρ sin ϕdρdθ, dSθ = ρdρdϕ
(2.25)
三、 方向導數 — 梯度
正如本節之標題: 梯度本質上是方向導數。 給定方向 ν 則純量函數 f 之方向導數定義為
∂f
∂ν = ∇f · ν (3.1)
雖然這是一個等式但左右兩邊卻有本質上的差異。 因為方向導數 ∂f∂ν 與座標之選取無關, 所以 要將梯度推廣到曲線座標系統必須從方向導數著手。
圖解梯度、 散度與旋度 39
圖6: (u1, u2, u3) 到 (u1+ △u1, u2+ △u2, u3+ △u3) 之變化
由 (3.1) 直觀而言, 梯度 ∇f 就是方向導數 ∂f∂ν 除以 ν! 考慮 f 從 (u1, u2, u3) 到 (u1 + △u1, u2+ △u2, u3+ △u3) 之變化並利用連接這兩點之弦的向量 △r = h1△u1e1+ h2△u2e2+ h3△u3e3 得 (請參考圖 6)
△f = f (u1+ △u1, u2+ △u2, u3+ △u3) − f (u1, u2, u3)
= ∂f
∂u1
△u1+ ∂f
∂u2
△u2+ ∂f
∂u3
△u3+ · · ·
= 1 h1
∂f
∂u1
h1△u1+ 1 h2
∂f
∂u2
h2△u2+ 1 h3
∂f
∂u3
h3△u3+ · · ·
= 1 h1
∂f
∂u1
e1+ 1 h2
∂f
∂u2
e2+ 1 h3
∂f
∂u3
e3
· △r + · · · 再除以長度 △s = |△r| = ph21(△u1)2+ h22(△u2)2+ h23(△u3)2
△f
△s = 1 h1
∂f
∂u1
e1+ 1 h2
∂f
∂u2
e2+ 1 h3
∂f
∂u3
e3
·△r
△s + · · · 取極限 (△s → 0) 將弦變為單位切向量 ν = drds 得方向導數
∂f
∂ν = 1 h1
∂f
∂u1
e1+ 1 h2
∂f
∂u2
e2+ 1 h3
∂f
∂u3
e3
· ν (3.2)
比較 (3.1)、 (3.2) 兩式得
∇f = 1 h1
∂f
∂u1
e1+ 1 h2
∂f
∂u2
e2+ 1 h3
∂f
∂u3
e3 (3.3)
藉由類比, 實際上 (3.3) 這個公式是很容易猜得出來。 只需要將 {i, j, k} 換為 {e1, e2, e3} 、 {x, y, z} = {x1, x2, x3} 換為 {u1, u2, u3}, 唯一要注意的是 ∂ui 一定要伴隨著 hi, 所以每一 項都有 h1
i 是自然的。 而背後真正的理由則是: 量綱平衡。
圖7: (r, θ, z) 到 (r + △r, θ + △θ, z + △z) 之變化
例題 3.1. 圓柱與球座標之梯度 (請參考圖 7, 8)
∇f = ∂f
∂xi+ ∂f
∂yj+ ∂f
∂zk
=∂f
∂rer+1 r
∂f
∂θeθ+ ∂f
∂zez
=∂f
∂ρeρ+ 1 ρ
∂f
∂ϕeϕ+ 1 ρ sin ϕ
∂f
∂θeθ
簡單的量綱分析可以看出上式的合理性, 我們留給讀者去 看出來。 讀者當然可以由 (3.3) 直接 代入驗證上式, 但是我們並不鼓勵這種被動的學習方式, 反而是從方向導數著手, 類似前面的推 導方式自己走過一遍, 因為數學是非常的存在主義: 沒有人可以代替你相信!
四、 Gauss 散度定理 — 散度
如果從散度的原始定義 (2.1)1, 它只告訴我們微分而已甚麼東西都看不到。 要瞭解散度就 必須從 Gauss 散度定理著手
Z Z Z
V
divAdV = Z Z
S
A· ndS, S = ∂V (4.1) 兩邊同時除以體積 |V |, 然後再令 |V | → 0 得
divA = lim
|V |→0
1
|V | Z Z
S
A· ndS (4.2)
圖解梯度、 散度與旋度 41
圖8: (ρ, ϕ, θ)到(ρ + △ρ, ϕ + △ϕ, θ + △θ) 之變化
所以散度 divA 就是單位體積向量場 A 之通量 (flux)。 這個定義除了物理意義清楚之外, 也 告訴我們散度與所選取之座標無關, 換言之高斯散度定理提供我們絕佳的散度之定義。
散度 ⇐⇒ 通量(曲面積分) ⇐⇒ Gauss 散度定理 簡單的量綱分析得
[divA] = 1 [|V |]
Z Z
S
[A · ndS] = 1
L3[A]L2 = [A]
L
假設 A = A1e1+ A2e2+ A3e3, 我們先考慮 e1-方向之通量。 此時向量場由 u1 到 u1+ △u1, 利用 △S1 = (h2△u2)(h3△u3) 得
Z Z
S1
A· ndS = (A1△S1)(u1+ △u1) − (A1△S1)(u1)
≈ [(A1h2h3)(u1+ △u1) − (A1h2h3)(u1)] △u2△u3
= ∂
∂u1
(A1h2h3)△u1△u2△u3 (均值定理)
= 1 h1h2h3
∂
∂u1
(A1h2h3)h1△u1h2△u2h3△u3
這裡 A1, h2, h3 是 u1-函數 (因為 u2, u3固定)。 e1-方向有兩個面: 前面與後面, 負號代表後面。
同理可以考慮 e2, e3-方向之通量, 由於 △S2= (h3△u3)(h1△u1)、 △S3 = (h1△u1)(h2△u2) 得
Z Z
S2
A· ndS = (A2△S2)(u2+ △u2) − (A2△S2)(u2)
= ∂
∂u3
(A3h1h2)△u3△u1△u2
= 1 h1h2h3
∂
∂u3(A3h1h2)h3△u3h1△u1h2△u2
這裡 A3, h1, h2 視為 u3-函數 (因為 u1, u2 固定)。 將這三部分通量相加之後再除以體積 h1△u1h2△u2h3△u3 得
divA = 1 h1h2h3
∂
∂u1
(A1h2h3) + ∂
∂u2
(A2h3h1) + ∂
∂u3
(A3h1h2)
(4.3) 只要記得方程式會講話則 (4.3) 實際上是非常直觀且容易記憶的。 以第一項 (e1-方向) 而言
1 h1h2h3
∂
∂u1
(A1h2h3)
把 h1h2h3 視為體積、 h2h3 為底面積則 A1h2h3 是通量, 但有兩個面其朝外法向量正好一正一 負兩個通量必須相減, 也就是差分, 在無窮小的時候 (u1 → 0) 差分變微分, 所以出現 ∂u∂1。 同 理第二、 三項也是如此分析。 表面上 (4.3) 雖然只有三項但本質上是六項, 正是無窮小方塊的六 個面。 所以 (4.3) 告訴我們的正是: 單位體積之通量 (三個方向)。
例題 4.1. 圓柱座標之散度: A= Arer+ Aθeθ+ Azez divA = 1
r
∂
∂r(Arr) + ∂
∂θ(Aθ) + ∂
∂z(Azr)
(4.4)
我們利用圖解(請參考圖 10) 來看出三個方向之通量:
er : Ar· r · 1(r: 大) − Ar· r · 1(r: 小) ≈ ∂
∂r(Arr) eθ : Aθ· 1 · 1(θ: 大) − Aθ· 1 · 1(θ: 小) ≈ ∂
∂θ(Aθ) ez : Az · r · 1(z: 大) − Az· r · 1(z: 小) ≈ ∂
∂z(Azr) 再除以體積 h1h2h3 = hrhθhz = 1 · r · 1 = r 就是 (4.4)。
圖解梯度、 散度與旋度 43
圖9: ei, i= 1, 2, 3,之通量
例題 4.2. 球座標之散度: A = Aρeρ+ Aϕeϕ+ Aθeθ divA = 1
ρ2sin ϕ
∂
∂ρ(Aρρ2sin ϕ) + ∂
∂ϕ(Aϕρ sin ϕ) + ∂
∂θ(Aθρ)
(4.5) 此時體積為 h1h2h3 = hρhϕhθ = ρ2sin ϕ, 而三個方向之通量為 (請參考圖 11, 12)
eρ : Aρ· ρ · ρ sin ϕ(ρ: 大) − Aρ· ρ · ρ sin ϕ(ρ: 小) ≈ ∂
∂ρ(Aρρ2sin ϕ) eϕ : Aϕ· ρ sin ϕ · 1(ϕ: 大) − Aϕ· ρ sin ϕ · 1(ϕ: 小) ≈ ∂
∂ϕ(Aϕρ sin ϕ) eθ : Aθ· ρ · 1(θ: 大) − Aθ· ρ · 1(θ: 小) ≈ ∂
∂θ(Aθρ)
圖10: 圓柱座標er, eθ, ez,之通量
圖11: 球座標 eρ 之通量
五、 Stokes 定理 — 旋度
同樣從旋度的原始定義 (2.1)2, 甚麼東西都看不到。 旋度直接聯想就是環流 (circulation)
圖解梯度、 散度與旋度 45
圖12: 球座標 eϕ, eθ之通量
或者是功 (work) 也就是線積分, 所以要瞭解旋度就必須從 Stokes 定理開始 I
C
A· dr = Z Z
S
curl A · ndS, C = ∂S (5.1) 要提醒讀者的是圍道C的方向與朝外法向量 n 必須滿足右手法則 (如果將 n 當作大拇指的方 向則圍道 C 是以逆時針方向為正)。 (5.1) 兩邊同時除以曲面面積 |S|, 然後再令 |S| → 0 得
curl A · n = lim
|S|→0
1
|S|
I
C
A· dr (5.2)
因為 |curl A · n| ≤ |curl A|, 所以旋度 curl A 可以詮釋為單位面積向量場 A 之最大環流。
由此可知旋度 curl A 之量綱 (或因次; dimension) 等於 A 之量綱除以長度 [curl A] = [A · t]L
L2 = [A]
L
因為線積分與座標無關, 所以旋度與所選取之座標無關, 換言之 Stokes 定理提供我們絕佳的旋 度之定義
旋度 ⇐⇒ 環流(線積分) ⇐⇒ Stokes定理
假設 A = A1e1+ A2e2+ A3e3。 我們先考慮 e1-方向, 此時 S1 是前面 e1 ⊥ S1, C1 = ∂S1, |S1| = (h2△u2)(h3△u3) 按照右手法則 C1 有四個邊, 所以線積分有四項
I
C1
A· dr ≈ (A2h2)(u3)△u2+ (A3h3)(u2+ △u2)△u3
h2h3 ∂u2 ∂u3
同理可以得 e2, e3-方向的分量。 令 S1 是右面, 則
e2 ⊥ S2, C2 = ∂S2, |S2| = (h3△u3)(h1△u1) I
C2
A· dr ≈ (A3h3)(u1)△u3+ (A1h1)(u3+ △u3)△u1
−(A3h3)(u1+ △u1)△u3− (A1h1)(u3)△u1
= ∂(A1h1)
∂u3 − ∂(A3h3)
∂u1
△u3△u1
= 1 h3h1
∂(A1h1)
∂u3 − ∂(A3h3)
∂u1
(h3△u3)(h1△u1)
這裡 A1h1, A3h3 分別視為 u3 與 u1 的函數。 所以 ∇ × A 在 e2-方向 (右面) 的分量為
∇ × A|e2 = curl A|e2 = 1 h3h1
∂(A1h1)
∂u3
− ∂(A3h3)
∂u1
(5.4) 令 S3 是上面, 則
e3 ⊥ S3, C3 = ∂S3, |S3| = (h1△u1)(h2△u2) I
C3
A· dr ≈ (A1h1)(u2)△u1+ (A2h2)(u1+ △u1)△u2
−(A1h1)(u2+ △u2)△u1− (A2h2)(u1)△u2
= ∂(A2h2)
∂u1 − ∂(A1h1)
∂u2
△u1△u2
= 1 h1h2
∂(A2h2)
∂u1
− ∂(A1h1)
∂u2
(h1△u1)(h2△u2)
這裡 A2h2, A1h1 分別視為 u1 與 u2 的函數。 所以 ∇ × A 在 e3-方向 (右面) 的分量為
∇ × A|e3 = curl A|e3 = 1 h1h2
∂(A2h2)
∂u1
− ∂(A1h1)
∂u2
(5.5)
圖解梯度、 散度與旋度 47
將 (5.3)–(5.5) 合併得
curl A = 1 h2h3
∂
∂u2
(A3h3) − ∂
∂u3
(A2h2) e1 + 1
h3h1
∂
∂u3(A1h1) − ∂
∂u1(A3h3) e2 + 1
h1h2
∂
∂u1
(A2h2) − ∂
∂u2
(A1h1)
e3 (5.6)
再由 Laplace 降階法結論
curl A = ∇ × A = 1 h1h2h3
h1e1 h2e2 h3e3
∂
∂u1
∂
∂u2
∂
∂u3
A1h1A2h2A3h3
(5.7)
註解:
(i) 在計算過程中雖然我們只考慮三個面, 但一共有十二個邊, 因此這個立方體所有面之邊界 的線積分都已經考慮到了, 所以最後之結果確實是全部之環流。
(ii) 因為hi必須配合∂ui, 為了保持量綱平衡A的旋度, 例如第一項 (e1方向) 1
h2h3
∂
∂u2
(A3h3) − ∂
∂u3
(A2h2) e1 分子分母的h3、h2分別彼此抵消而1
h2、1
h3則正好配合 ∂
∂u2、 ∂
∂u3, 而且可以驗證其量綱 (因次) 等於[A1]/L 與 [curl A] = [A]/L 相同, 至於正負符號則由右手法則 來決定
∂
∂u2
(A3h3)e1=⇒ (2 → 3 → 1) (+)
∂
∂u3
(A2h2)e1=⇒ (3 → 2 → 1) (−)
(iii) 由 Stokes 定理知道旋度是單位面積之環流, 事實上 (5.6) 告訴我們的正是如此, 看第一 項 (e1-方向, e1為法向量的面)
1 h2h3
∂
∂u2
(A3h3) − ∂
∂u3
(A2h2)
此時 e3-方向之線積分 (環流) 為 A3h3(u2大) 與 −A3h3(u2小), 當u2 → 0時這兩項 合併得 ∂u∂
2(A3h3), (差分變微分!)。 同理 e2-方向之線積分 (環流) 為 −A2h2(u3大) 與 A2h2(u2小), 當 u3 → 0時這兩項合併得 −∂u∂3(A2h2), 因此偏導數是說兩個線積分之差 (正負兩個方向), h2h3視為這個長方形之面積, 整個合併起來就是e1-方向的單位面積之環 流。 同理可以討論e2-方向與e3-方向。
圖13: ei, i= 1, 2, 3 三個方向之環流
我們利用圖解 (請參考圖 13) 來看出三個方向之環流:
I
C1
A· dr ≈ A2h2(u3小) + A3h3(u2大) − A2h2(u3大) − A3h3(u2小)
≈ ∂
∂u2(A3h3) − ∂
∂u3(A2h2) I
C2
A· dr ≈ A3h3(u1小) + A1h1(u3大) − A3h3(u1大) − A1h1(u3小)
≈ ∂
∂u3
(A1h1) − ∂
∂u1
(A3h3) I
C3
A· dr ≈ A1h1(u2小) + A2h2(u1大) − A1h1(u2大) − A2h2(u1小)
≈ ∂
∂u1
(A1h1) − ∂
∂u2
(A1h1) 再分別除以 h2h3, h3h1, h1h2 就是 (5.3)–(5.5)。
圖解梯度、 散度與旋度 49
例題 5.1. 圓柱座標之旋度: A= Arer+ Aθeθ+ Azez
curl A = 1 r
∂Az
∂θ − r∂Aθ
∂z
er+ ∂Ar
∂z −∂Az
∂r
eθ+ 1 r
∂
∂r(rAθ) − ∂Ar
∂θ
ez (5.8) 我們利用圖解 (請參考圖 14) 推得三個方向之環流:
I
Cr
A· dr ≈ Aθr(z小) + Az· 1(θ大) − Aθr(z大) − Az(θ小)
≈ ∂
∂θ(Az) − ∂
∂z(Aθr) I
Cθ
A· dr ≈ Az(r小) + Ar(z大) − Az(r大) − Ar(z小)
≈ ∂
∂z(Ar) − ∂
∂r(Az) I
Cz
A· dr ≈ Ar(θ小) + Aθr(r大) − Ar(θ大) − Aθr(r小)
≈ ∂
∂r(Aθr) − ∂
∂θ(Ar) 再分別除以 r, 1, r (視為 Sr, Sθ, Sz 之面積) 就是 (5.8)。
例題 5.2. 球座標之旋度: A = Aρeρ+ Aϕeϕ+ Aθeθ curl A = 1
ρ2sin ϕ
ρ ∂
∂ϕ(sin ϕAθ) − ρ∂Aϕ
∂θ
eρ + 1
ρ sin ϕ
∂Aρ
∂θ − sin ϕ ∂
∂ρ(ρAθ)
eϕ+ 1 ρ
∂
∂ρ(ρAϕ) −∂Aρ
∂ϕ
eθ (5.9) 我們利用圖解 (請參考圖 15) 推得三個方向之環流:
I
Cρ
A· dr ≈ Aϕρ(θ小) + Aθρ sin ϕ(ϕ大) − Aϕρ(θ大) − Aθρ sin ϕ(ϕ小)
≈ ∂
∂ϕ(Aθρ sin ϕ) − ∂
∂θ(Aϕρ) I
Cϕ
A· dr ≈ Aθρ sin ϕ(ρ小) + Aρ(θ大) − Aθρ sin ϕ(ρ大) − Aρ(θ小)
≈ ∂
∂θ(Aρ) − ∂
∂ρ(Aθρ sin ϕ) I
Cθ
A· dr ≈ Aρ(ϕ小) + Aϕρ(ρ大) − Aρ(ϕ大) − Aϕρ(ρ小)
≈ ∂
∂ρ(Aϕρ) − ∂
∂ϕ(Aρ)
再分別除以 ρ2sin ϕ, ρ sin ϕ, ρ (視為 Sρ, Sϕ, Sθ 之面積) 就是 (5.9)。
圖14: 圓柱座標三個方向之環流
六、Laplace 算子
根據 Laplace 算子的定義
∆ = div∇ = ∇ · ∇ (6.1)
所以由梯度與散度的公式 (3.3) 與 (4.3) 得
∆f = div 1 h1
∂f
∂u1e1+ 1 h2
∂f
∂u2e2+ 1 h3
∂f
∂u3e3
= 1 h1h2h3
∂
∂u1
h2h3
h1
∂f
∂u1
+ ∂
∂u2
h3h1
h2
∂f
∂u2
+ ∂
∂u3
h1h2
h3
∂f
∂u3
(6.2) 雖然這樣可以得到答案, 但對於瞭解 Laplace 算子並沒有幫助, 頂多是散度與梯度的組合。 回 到 Gauss 散度定理
Z Z Z
V
∆f dV = Z Z Z
V
div(∇f )dV = Z Z
S
∇f · ndS = Z Z
S
∂f
∂ndS (6.3)
圖解梯度、 散度與旋度 51
圖15: 球座標三個方向之環流
兩邊同時除以體積 |V |, 然後再令 |V | → 0 得
∆f = lim
|V |→0
1
|V | Z Z
S
∂f
∂ndS (6.4)
所以 ∆f 就是單位體積 f 之法向量的方向導數 ∂f∂n 之通量。 量綱分析得 [∆f ] = [f ]
L2 = 1 [|V |]
∂f
∂n
[dS] = 1 L3
[f ]
L L2 = [f ] L2 我們利用圖解來推得 (6.2)。 ∂n∂f 在三個方向之通量為
e1-方向 (S1曲面): 1 h1
∂f
∂u1
h2h3(u1大) − 1 h1
∂f
∂u1
h2h3(u1小) ≈ ∂
∂u1
h2h3
h1
∂f
∂u1
e2-方向 (S2曲面): 1 h2
∂f
∂u2
h3h1(u2大) − 1 h2
∂f
∂u2
h3h1(u2小) ≈ ∂
∂u2
h3h1
h2
∂f
∂u2
h1 ∂u1,
h2 ∂u2,
h3 ∂u3 =⇒ e1, e2, e3-方向之通量
∂
∂ui
, i = 1, 2, 3 =⇒ 差分變微分 例題 6.1. 極座標之 Laplace 算子
∆f = ∂2f
∂x2 +∂2f
∂y2 = ∂2f
∂r2 +1 r
∂f
∂r + 1 r2
∂2f
∂θ2 (6.5)
我們仍然利用圖解 (請參考圖 16) 推 (6.5)
∂f
∂rr(r大) − ∂f
∂rr(r小) ≈ ∂
∂r
r∂f
∂r
1 r
∂f
∂θ(θ大) − 1 r
∂f
∂θ(θ小) ≈ 1 r
∂2f
∂θ2 這兩式相加之後除以 r 得
∆f = 1 r
∂
∂r
r∂f
∂r
+ 1 r2
∂2f
∂θ2 = ∂2f
∂r2 + 1 r
∂f
∂r + 1 r2
∂2f
∂θ2 例題 6.2. 圓柱座標之 Laplace 算子
∆f = ∂2f
∂x2 +∂2f
∂y2 +∂2f
∂z2 = ∂2f
∂r2 +1 r
∂f
∂r + 1 r2
∂2f
∂θ2 +∂2f
∂z2 (6.6) 圓柱座標幾乎就是極座標我們留給讀者練習。
例題 6.3. 球座標之 Laplace 算子
∆f = 1 ρ2sin ϕ
∂
∂ρ
∂f
∂ρρ2sin ϕ
+ ∂
∂ϕ
∂f
∂ϕsin ϕ
+ ∂
∂θ
1 sin ϕ
∂f
∂θ
= 1 ρ2
∂
∂ρ
ρ2∂f
∂ρ
+ 1 ρ2sin ϕ
∂
∂ϕ
sin ϕ∂f
∂ϕ
+ 1
ρ2sin2ϕ
∂2f
∂θ2 (6.7)
圖解梯度、 散度與旋度 53
圖16: 極座標之 Laplace 算子
讀者可以驗證一下量綱關係 1 ρ2
∂
∂ρ
ρ2∂f
∂ρ
≈ 1 L2
1 L
L2[f ]
L
= [f ] L2 1
ρ2sin ϕ
∂
∂ϕ
sin ϕ∂f
∂ϕ
≈ 1 L2· 1
1 1
1[f ]
1
= [f ] L2 1
ρ2sin2ϕ
∂2f
∂θ2 ≈ 1 L212
[f ] 12 = [f ]
L2 首先由例題 3.1 得
∇f = ∂f
∂ρeρ+ 1 ρ
∂f
∂ϕeϕ+ 1 ρ sin ϕ
∂f
∂θeθ 所以 (6.7) 圖解如下:
∂f
∂ρ · ρ · ρ sin ϕ(ρ大) − ∂f
∂ρ · ρ · ρ sin ϕ(ρ小) ≈ ∂
∂ρ
∂f
∂ρρ2sin ϕ
1 ρ
∂f
∂ϕ · ρ sin ϕ · 1(ϕ大) −1 ρ
∂f
∂ϕ· ρ sin ϕ · 1(ϕ小) ≈ ∂
∂ϕ
∂f
∂ϕsin ϕ
1
ρ sin ϕ
∂f
∂θ · 1 · ρ(θ大) − 1 ρ sin ϕ
∂f
∂θ · 1 · ρ(θ小) ≈ ∂
∂θ
1 sin ϕ
∂f
∂θ
這三式相加再除以 ρ2sin ϕ 就是 (6.7)。
例題 6.4. n 維球座標之 Laplace 算子
∆f = 1 ωnrn−1
∂
∂r
∂f
∂rωnrn−1
+ · · ·
= ∂2f
∂r2 +n − 1 r
∂f
∂r + 1
r2∆∗f (6.8) 這裡 ωn 是 n 維空間的單位球面之表面積
ωn= |Sn−1(1)| = 2πn2
Γ(n2) (6.9)
圖17: n 維球座標之 Laplace 算子
而 ∆∗ 則是單位球面 Sn−1(1) 上的 Laplace 算子。
n 維球座標其實不難想像。 與極座標的扇形方塊類比我們可以想像此時是一個圓錐形的冰 淇淋半徑方向的長度還是 dr, 弧長 rdθ 則取代為半徑等於 r 之球面的一部分其表面積可以設 為 kωnrn−1 (k 與頂角有關), 因此半徑方向之法向方向導數 ∂f∂r 的通量 (請參考圖 17)
∂f
∂rkωnrn−1(r大) −∂f
∂rkωnrn−1(r小) ≈ k ∂
∂r
∂f
∂rωnrn−1
再除以體積 1 · kωnrn−1 就是 (6.8) 第一個等式右邊的第一項。
註解:
(i) 從計算過程中可以看得出來在 (6.8) 出現的 n − 1 是半徑為 r 的球表面積對 r 微分的結 果, 其實就是球面的維數 (dimension)。
(ii) 當 n = 2, 3 時由 (6.5), (6.7) 得
∆∗f = ∂2f
∂θ2, ∆∗f = ∂2f
∂ϕ2 + cos ϕ sin ϕ
∂f
∂ϕ + 1 sin2ϕ
∂2f
∂θ2 (6.10)
誌謝
本文是作者於 2009年暑假在中研院數學所 《數學名題及其故事》 一系列演講之一。 後來更 在中研院數學所星期一的分析討論班中詳細地以圖解的方式講述了一次, 在此特別謝謝李志豪 教授的鼓勵與安排。
圖解梯度、 散度與旋度 55