三維情況下,拉普拉斯方程可由下面的形式描述,問題歸結為求解對實自變量x 、y 、z 二階可微 的實函數φ:
使用笛卡爾坐標 ,
Δ f = ∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2 + ∂ 2 f ∂ z 2 = 0 {\displaystyle \Delta f={\frac {\partial ^{2}f}{\partial x^{2}}}+{\frac {\partial ^{2}f}{\partial y^{2}}}+{\frac {\partial ^{2}f}{\partial z^{2}}}=0} 。使用柱坐標 ,
Δ f = 1 r ∂ ∂ r ( r ∂ f ∂ r ) + 1 r 2 ∂ 2 f ∂ ϕ 2 + ∂ 2 f ∂ z 2 = 0 {\displaystyle \Delta f={\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r{\frac {\partial f}{\partial r}}\right)+{\frac {1}{r^{2}}}{\frac {\partial ^{2}f}{\partial \phi ^{2}}}+{\frac {\partial ^{2}f}{\partial z^{2}}}=0} 使用球面坐標 ,
Δ f = 1 ρ 2 ∂ ∂ ρ ( ρ 2 ∂ f ∂ ρ ) + 1 ρ 2 sin θ ∂ ∂ θ ( sin θ ∂ f ∂ θ ) + 1 ρ 2 sin 2 θ ∂ 2 f ∂ φ 2 = 0 {\displaystyle \Delta f={\frac {1}{\rho ^{2}}}{\frac {\partial }{\partial \rho }}\left(\rho ^{2}{\frac {\partial f}{\partial \rho }}\right)+{\frac {1}{\rho ^{2}\sin \theta }}{\frac {\partial }{\partial \theta }}\left(\sin \theta {\frac {\partial f}{\partial \theta }}\right)+{\frac {1}{\rho ^{2}\sin ^{2}\theta }}{\frac {\partial ^{2}f}{\partial \varphi ^{2}}}=0} 。使用曲線坐標 ,
Δ f = ∂ ∂ ξ j ( ∂ f ∂ ξ k g k i ) + ∂ f ∂ ξ j g j m Γ m n n = 0 , {\displaystyle \Delta f={\frac {\partial }{\partial \xi ^{j}}}\left({\frac {\partial f}{\partial \xi ^{k}}}g^{ki}\right)+{\frac {\partial f}{\partial \xi ^{j}}}g^{jm}\Gamma _{mn}^{n}=0,} 或
Δ f = 1 | g | ∂ ∂ ξ i ( | g | g i j ∂ f ∂ ξ j ) = 0 , ( g = d e t { g i j } ) {\displaystyle \Delta f={\frac {1}{\sqrt {|g|}}}{\frac {\partial }{\partial \xi ^{i}}}\!\left({\sqrt {|g|}}g^{ij}{\frac {\partial f}{\partial \xi ^{j}}}\right)=0,\qquad (g=\mathrm {det} \{g_{ij}\})} 。這組方程又經常寫為
∇ 2 φ = 0 {\displaystyle \nabla ^{2}\varphi =0} 或者
div grad φ = 0 {\displaystyle \operatorname {div} \,\operatorname {grad} \,\varphi =0} ;其中,div表示矢量場 的散度 (結果是一個純量場 ),grad表示純量場的梯度 (結果是一個矢量場 )。
這方程又可寫為
Δ φ = 0 {\displaystyle \Delta \varphi =0} ;其中,Δ稱為拉普拉斯算子 。
拉普拉斯方程的解稱為調和函數 。[1] :671-672
如果等號右邊是一個給定的函數f (x , y , z ),即
Δ φ = f {\displaystyle \Delta \varphi =f} ,則該方程稱為泊松方程 。拉普拉斯方程和泊松方程是最簡單的橢圓型微分方程 。偏微分算子∇ 2 {\displaystyle \nabla ^{2}} 或Δ {\displaystyle \Delta } (可以在任意維空間中定義這樣的算子)稱為拉普拉斯算子 。
邊界條件
編輯
二維拉普拉斯方程
編輯
兩個自變量的拉普拉斯方程具有以下形式:
∂ 2 ψ ∂ x 2 + ∂ 2 ψ ∂ y 2 ≡ ψ x x + ψ y y = 0. {\displaystyle {\frac {\partial ^{2}\psi }{\partial x^{2}}}+{\frac {\partial ^{2}\psi }{\partial y^{2}}}\equiv \psi _{xx}+\psi _{yy}=0.} 解析函數
編輯
解析函數 的實部和虛部均滿足拉普拉斯方程。換言之,若 z = x + iy ,並且
f ( z ) = u ( x , y ) + i v ( x , y ) {\displaystyle f(z)=u(x,y)+iv(x,y)} ,那麼 f (z ) 是解析函數的充要條件 是 u(x,y),v(x,y) 可微,且滿足下列柯西-黎曼方程 :[1] :671-672
u x = v y , v x = − u y {\displaystyle u_{x}=v_{y},\quad v_{x}=-u_{y}} 。上述方程繼續求導就得到
u y y = ( − v x ) y = − ( v y ) x = − ( u x ) x {\displaystyle u_{yy}=(-v_{x})_{y}=-(v_{y})_{x}=-(u_{x})_{x}} 。所以 u 滿足拉普拉斯方程。類似的計算可推得 v 同樣滿足拉普拉斯方程。
反之,給定一個由解析函數(或至少在某點及其鄰域內解析的函數)f (z ) 的實部確定的調和函數,若寫成下列形式:
f ( z ) = φ ( x , y ) + i ψ ( x , y ) {\displaystyle f(z)=\varphi (x,y)+i\psi (x,y)} ,則等式
ψ x = − φ y , ψ y = φ x {\displaystyle \psi _{x}=-\varphi _{y},\quad \psi _{y}=\varphi _{x}} 。成立就可使得柯西-黎曼方程得到滿足。上述關係無法確定ψ,只能得到它的微增量表達式:
d ψ = − φ y d x + φ x d y {\displaystyle d\psi =-\varphi _{y}\,dx+\varphi _{x}\,dy} 。φ 滿足拉普拉斯方程意味著 ψ 滿足可積條件:
ψ x y = ψ y x {\displaystyle \psi _{xy}=\psi _{yx}} 。所以可以通過一個線積分來定義 ψ。可積條件和斯托克斯定理 的滿足說明線積分的結果與積分經過的具體路徑無關,僅由起點和終點決定。於是,我們便通過複變函數 方法得到了 φ 和 ψ 這一對拉普拉斯方程的解。這樣的解稱為一對共軛調和函數 。這種構造解的方法只在局部(複變函數 f (z )) 的解析域內)有效,或者說,構造函數的積分路徑不能圍繞有 f (z ) 的奇點 。譬如,在極坐標 平面 (r ,θ ) 上定義函數
φ = log r {\displaystyle \varphi =\log r} ,那麼相應的解析函數為
f ( z ) = log z = log r + i θ {\displaystyle f(z)=\log z=\log r+i\theta } 。在這裡需要注意的是,極角 θ 僅在不包含原點的區域內才是單值的。
拉普拉斯方程與解析函數之間的緊密聯繫說明拉普拉斯方程的任何解都無窮階可導(這是解析函數的一個性質),因此可以展開成冪級數 形式,至少在不包含奇點的圓域內是如此。這與波動方程 的解形成鮮明對照,後者包含任意 函數,其中一些的可微分階數是很小的。
冪級數和傅立葉級數 之間存在著密切的關係。如果我們將函數 f 在複平面上以原點為中心,R 為半徑的圓域內展開成冪級數,即
f ( z ) = ∑ n = 0 ∞ c n z n {\displaystyle f(z)=\sum _{n=0}^{\infty }c_{n}z^{n}} ,將每一項係數適當地分離出實部和虛部
c n = a n + i b n {\displaystyle c_{n}=a_{n}+ib_{n}} 。那麼
f ( z ) = ∑ n = 0 ∞ [ a n r n cos n θ − b n r n sin n θ ] + i ∑ n = 1 ∞ [ a n r n sin n θ + b n r n cos n θ ] {\displaystyle f(z)=\sum _{n=0}^{\infty }\left[a_{n}r^{n}\cos n\theta -b_{n}r^{n}\sin n\theta \right]+i\sum _{n=1}^{\infty }\left[a_{n}r^{n}\sin n\theta +b_{n}r^{n}\cos n\theta \right]} ,這便是 f 的傅立葉級數。這些三角函數自身也可以用倍角公式 展開。
流體動力學
編輯
設 u {\displaystyle u} 、v {\displaystyle v} 分別為滿足定常 、不可壓縮 和無旋 條件的流體速度場的 x {\displaystyle x} 和 y {\displaystyle y} 方向分量(這裡僅考慮二維流場),那麼不可壓縮條件為:[3] :99-101
u x + v y = 0 {\displaystyle u_{x}+v_{y}=0} ,無旋條件為:
v x − u y = 0 {\displaystyle v_{x}-u_{y}=0} 。若定義一個純量函數 ψ {\displaystyle \psi } ,使其微分滿足:
d ψ = − v d x + u d y {\displaystyle d\psi =-v\,dx+u\,dy} ,那麼不可壓縮條件便是上述微分式的可積條件。積分的結果函數 ψ {\displaystyle \psi } 稱為流函數 ,因為它在同一條流線 上各點的值是相同的。ψ {\displaystyle \psi } 的一階偏導為:
ψ x = − v , ψ y = u {\displaystyle \psi _{x}=-v,\quad \psi _{y}=u} ,無旋條件即令 ψ {\displaystyle \psi } 滿足拉普拉斯方程。ψ {\displaystyle \psi } 的共軛調和函數 φ {\displaystyle \varphi } 稱為速度勢 。柯西-黎曼方程要求
φ x = u , φ y = v {\displaystyle \varphi _{x}=u,\quad \varphi _{y}=v} 。所以每一個解析函數都對應著平面內的一個定常不可壓縮無旋流場。解析函數的實部為速度勢函數,虛部為流函數。
根據麥克斯韋方程組 ,二維空間中不隨時間變化的電場 (u ,v ) 滿足:[4] :83
∇ × ( u , v ) = v x − u y = 0 {\displaystyle \nabla \times (u,v)=v_{x}-u_{y}=0} ,和
∇ ⋅ ( u , v ) = ρ {\displaystyle \nabla \cdot (u,v)=\rho } ,其中 ρ 為電荷密度 。第一個麥克斯韋方程便是下列微分式的可積條件:
d φ = − u d x − v d y {\displaystyle d\varphi =-u\,dx-v\,dy} ,所以可以構造電勢函數 φ 使其滿足
φ x = − u , φ y = − v {\displaystyle \varphi _{x}=-u,\quad \varphi _{y}=-v} 。第二個麥克斯韋方程即:
φ x x + φ y y = − ρ {\displaystyle \varphi _{xx}+\varphi _{yy}=-\rho } ,這是一個泊松方程 ,當空間不包含自由電荷時,方程等號右邊變為0,方程變為拉普拉斯方程。
三維拉普拉斯方程
編輯
拉普拉斯方程的基本解滿足
∇ ⋅ ∇ u = u x x + u y y + u z z = − δ ( x − x ′ , y − y ′ , z − z ′ ) {\displaystyle \nabla \cdot \nabla u=u_{xx}+u_{yy}+u_{zz}=-\delta (x-x',y-y',z-z')} ,其中的三維δ函數 代表位於( x ′ , y ′ , z ′ ) {\displaystyle (x',\,y',\,z')} 的一個點源。
由基本解的定義,若對u 作用拉普拉斯算子,再把結果在包含點源的任意體積內積分,那麼
∭ V ∇ ⋅ ∇ u d V = − 1 {\displaystyle \iiint _{V}\nabla \cdot \nabla udV=-1} 。由於坐標軸旋轉不改變拉普拉斯方程的形式,所以基本解必然包含在那些僅與到點源距離r 相關的解中。如果我們選取包含點源、半徑為a 的球形域作為積分域,那麼根據高斯散度定理 [1] :318-322
− 1 = ∭ V ∇ ⋅ ∇ u d V = ∬ S u r d S = 4 π a 2 u r ( a ) {\displaystyle -1=\iiint _{V}\nabla \cdot \nabla u\,dV=\iint _{S}u_{r}dS=4\pi a^{2}u_{r}(a)} 。求得在以點源為中心,半徑為r 的球面上有
u r ( r ) = − 1 4 π r 2 {\displaystyle u_{r}(r)=-{\frac {1}{4\pi r^{2}}}} ,所以
u = 1 4 π r {\displaystyle u={\frac {1}{4\pi r}}} 。經過類似的推導同樣可求得二維形式的解
u = − ln r 2 π {\displaystyle u={\frac {-\ln r}{2\pi }}} 。格林函數
編輯
格林函數 是一種不但滿足前述基本解的定義,而且在體積域V 的邊界S 上還滿足一定的邊界條件的基本解。譬如,G ( x , y , z ; x ′ , y ′ , z ′ ) {\displaystyle G(x,y,z;x',y',z')\,} 可以滿足
∇ ⋅ ∇ G = − δ ( x − x ′ , y − y ′ , z − z ′ ) in V {\displaystyle \nabla \cdot \nabla G=-\delta (x-x',y-y',z-z')\quad {\hbox{in}}\quad V} ,G = 0 if ( x , y , z ) on S {\displaystyle G=0\quad {\hbox{if}}\quad (x,y,z)\quad {\hbox{on}}\quad S} 。現設u 為在V 內滿足泊松方程的任意解:
∇ ⋅ ∇ u = − f {\displaystyle \nabla \cdot \nabla u=-f} ,且u 在邊界S 上取值為g ,那麼我們可以應用格林定理 (是高斯散度定理的一個推論),得到[1] :652-659
∭ V [ G ∇ ⋅ ∇ u − u ∇ ⋅ ∇ G ] d V = ∭ V ∇ ⋅ [ G ∇ u − u ∇ G ] d V = ∬ S [ G u n − u G n ] d S {\displaystyle \iiint _{V}\left[G\,\nabla \cdot \nabla u-u\,\nabla \cdot \nabla G\right]\,dV=\iiint _{V}\nabla \cdot \left[G\nabla u-u\nabla G\right]\,dV=\iint _{S}\left[Gu_{n}-uG_{n}\right]\,dS} 。un 和Gn 分別代表兩個函數在邊界S 上的法向導數。考慮到u 和G 滿足的條件,可將這滿足狄利克雷邊界條件的公式化簡為
u ( x ′ , y ′ , z ′ ) = ∭ V G f d V − ∬ S G n g d S {\displaystyle u(x',y',z')=\iiint _{V}Gf\,dV-\iint _{S}G_{n}g\,dS} 。所以格林函數描述了量f 和g 對( x ′ , y ′ , z ′ ) {\displaystyle (x',y',z')} 點函數值的影響。
圓球殼案例
編輯
格林函數在半徑為a 的球面內的點上得值可以通過鏡像法 求得:距球心ρ的源點P 的通過球面的「反射鏡像」P' 距球心
ρ ′ = a 2 ρ {\displaystyle \rho '={\frac {a^{2}}{\rho }}} 。需要注意的是,如果P 在球內,那麼P' 將在球外。於是可得格林函數為
G = 1 4 π R − a 4 π ρ R ′ {\displaystyle G={\frac {1}{4\pi R}}-{\frac {a}{4\pi \rho R'}}} ;其中,R 表示距源點P 的距離,R' 表示距鏡像點P' 的距離。從格林函數上面的表示式可以推出泊松積分公式 。設ρ、θ和φ為源點P 的三個球坐標 分量。此處θ按照物理學界的通用標準定義為坐標矢徑與豎直軸(z 軸)的夾角(與歐洲習慣相同,與美國習慣不同)。於是球面內拉普拉斯方程的解為:[2] :64-65
u ( P ) = 1 4 π a 3 ( 1 − ρ 2 a 2 ) ∬ g ( θ ′ , φ ′ ) sin θ ′ d θ ′ d φ ′ ( a 2 + ρ 2 − 2 a ρ cos Θ ) 3 / 2 {\displaystyle u(P)={\frac {1}{4\pi }}a^{3}\left(1-{\frac {\rho ^{2}}{a^{2}}}\right)\iint {\frac {g(\theta ',\varphi ')\sin \theta '\,d\theta '\,d\varphi '}{(a^{2}+\rho ^{2}-2a\rho \cos \Theta )^{3/2}}}} ;其中,cos Θ = cos θ cos θ ′ + sin θ sin θ ′ cos ( θ − θ ′ ) {\displaystyle \cos \Theta =\cos \theta \cos \theta '+\sin \theta \sin \theta '\cos(\theta -\theta ')} 。
這個公式的一個顯見的結論是:若u 是調和函數,那麼u 在球心處的取值為其在球面上取值的平均。於是我們可以立即得出以下結論:任意一個調和函數(只要不是常函數)的最大值必然不會在其定義域的內部點取得。
參考文獻
編輯
^ 1.0 1.1 1.2 1.3 1.4 1.5 Boas, Mary. Mathematical Methods in the Physical Sciences 3rd. Wiley. 2005. ISBN 978-0471198260 .
^ 2.0 2.1 2.2 Jackson, John David, Classical Electrodynamic 3rd., USA: John Wiley & Sons, Inc., 1999, ISBN 978-0-471-30932-1
^ Batchelor, George. An Introduction to Fluid Dynamics . Cambridge University Press. 1999. ISBN 978-0521663960 .
^ Griffiths, David J., Introduction to Electrodynamics (3rd ed.), Prentice Hall, 1998, ISBN 0-13-805326-X
嚴鎮軍編,《數學物理方程》,第二版,中國科學技術大學出版社,合肥,2002,ISBN 978-7-312-00799-6 /O•177
L.C. Evans, Partial Differential Equations , American Mathematical Society, Providence, 1998. ISBN 978-0-8218-0772-9
I. G. Petrovsky, Partial Differential Equations , W. B. Saunders Co., Philadelphia, 1967.
A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists , Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 978-1-58488-299-2
A. Sommerfeld, Partial Differential Equations in Physics , Academic Press, New York, 1949.
Pijush K.Kundu, Fluid Mechanics, Academic Press, 2002. 外部連結
編輯