汪夢雅 曹 蓉 夏杰楨 吳 琪

(西藏大學(xué)理學(xué)院 西藏 拉薩 850000)

靜電場問題是電磁學(xué)理論之后又一個重要的理論知識，是電動力學(xué)里的一個核心問題[1,2].其基本問題就是求解關(guān)于靜電勢的泊松方程的邊值問題，解法主要有分離變量法、格林函數(shù)法、鏡像法和電多極矩法等，其中最常用的兩種解法是分離變量法和鏡像法[3].

由于電動力學(xué)的數(shù)學(xué)要求遠(yuǎn)比普通理工科專業(yè)高，所以為了能將學(xué)生從復(fù)雜的微分方程、偏微分方程的求解過程中解放出來，引入了Mathematica[4].Mathematica是美國WolframResearch公司開發(fā)的數(shù)學(xué)軟件，是現(xiàn)在使用最廣的軟件之一.該軟件除了很好地結(jié)合了數(shù)值和符號計算之外，還具有強(qiáng)大的繪圖功能，操作簡單，界面友好，代碼形式接近于自然語言，廣泛用于數(shù)學(xué)、物理、生物等領(lǐng)域[5～7].

本文利用Mathematica求解靜電場中分離變量法的問題出發(fā)，并討論在這種問題下的一些極限情況來探討Mathematica在靜電場中的應(yīng)用.通過Mathematica的使用消除學(xué)生對復(fù)雜公式的畏懼，使學(xué)生可以擺脫繁瑣的數(shù)學(xué)計算，并且通過Mathematica的繪圖使學(xué)生對物理圖像有個較為直觀準(zhǔn)確的認(rèn)識，激發(fā)學(xué)生的興趣，促進(jìn)課程的學(xué)習(xí).

1 靜電場的標(biāo)勢及其微分方程

1.1 靜電場的標(biāo)勢

眾所周知，相距為dl的兩點(diǎn)的電勢差為[8]

dφ=-E·dl

(1)

然后由于

(2)

因此，電場強(qiáng)度E等于電勢φ的負(fù)梯度

(3)

1.2 靜電場的微分方程和邊值關(guān)系

研究一個電荷對與它鄰近的電場是怎樣作用的，一點(diǎn)上的電場和它鄰近的電場又是怎么聯(lián)系的，即要找出電荷和電場相互作用規(guī)律的微分形式，而在導(dǎo)體表面或其他邊界上場和電荷的相互關(guān)系則由邊界條件反映出來.這種問題稱為邊值問題，即求微分方程滿足邊界條件的解.

在均勻同向的介質(zhì)中，D=εE，再把式(3)代入

(4)

其中ρ是自由電荷密度.式(4)是靜電勢滿足的基本微分方程，也就是泊松方程，給出邊界條件就可以確定電勢φ的解.

在兩介質(zhì)界面上，靜電勢所滿足的邊值關(guān)系分別為

φ1=φ2

(5)

(6)

以上是邊值關(guān)系的一般形式，但是在靜電問題中，常常有一些導(dǎo)體存在，由于導(dǎo)體的特殊性質(zhì)，在導(dǎo)體表面上的邊值關(guān)系有它的特點(diǎn).設(shè)導(dǎo)體表面所帶自由電荷面密度為σ，它外面的介質(zhì)電容率為ε，則由式(5)、(6)和導(dǎo)體靜電條件得到導(dǎo)體表面的邊界條件

φ=常量

(7)

(8)

靜電學(xué)的基本問題就是求出在每個均勻區(qū)域內(nèi)滿足泊松方程，在所有分界面上滿足邊值關(guān)系和在所研究的整個區(qū)域邊界上滿足邊界條件的電勢的解.

2 分離變量法求解的方程和通解

在許多實(shí)際問題中，靜電場是由帶電導(dǎo)體決定的[9].比如，電容器內(nèi)部的電場是由作為電極的兩個導(dǎo)體板上所帶電荷決定的，電子光學(xué)系統(tǒng)的靜電透鏡內(nèi)部，電場是由分布于電極上的自由電荷決定的.這些問題都有一個共同的特點(diǎn)，就是自由電荷只出現(xiàn)在一些導(dǎo)體的表面上，而在空間中沒有其他自由電荷分布.這時，當(dāng)選擇這些導(dǎo)體表面作為區(qū)域V的邊界，則在V內(nèi)部自由電荷密度ρ=0.此時泊松方程可以化為較為簡單的拉普拉斯方程(拉式方程)，即

(9)

因此，這一類問題的解法就是求滿足邊界條件的拉普拉斯方程的解.此拉普拉斯方程的通解可以用分離變量法求出，在求解前我們要先選擇合適的坐標(biāo)系，常用的坐標(biāo)系有球坐標(biāo)系和柱坐標(biāo)系，在這里我們寫出用球坐標(biāo)系得出的通解形式，球坐標(biāo)用(r,θ,φ)表示，其中r為點(diǎn)到原點(diǎn)的距離，θ為極角，φ為方位角，拉式方程在球坐標(biāo)中的通解為

(10)

(11)

其中Pn(cosθ)為勒讓德函數(shù)，an和bn是任意常數(shù)，由邊界條件或者說由邊值關(guān)系確定.

在每一個沒有電荷分布的區(qū)域內(nèi)，φ滿足拉普拉斯方程，其通解已由式(10)和式(11)給出，剩下的問題就是由邊界條件確定這些通解中所含的任意常數(shù)，得到滿足邊界條件的特解.

3 分離變量法的實(shí)際應(yīng)用

【例題】一半徑為a，介電常數(shù)為ε的介質(zhì)球放置在均勻電場E0中，如圖1所示，求介質(zhì)球內(nèi)、外的電位及電場.

圖1 題圖

分析：介質(zhì)球在外電場中極化，在它表面上產(chǎn)生束縛電荷.這些束縛電荷激發(fā)的電場疊加在原外電場E0上，得總電場E.束縛電荷分布和總電場E互相制約，邊界條件正確的反映這種制約關(guān)系.

球的半徑為a，球外為真空，這問題具有軸對稱性，對稱軸為通過球心沿外電場E0方向的軸線，取此軸為極軸.介質(zhì)球的存在使空間分為兩均勻區(qū)域，球外區(qū)域和球內(nèi)區(qū)域，兩區(qū)域內(nèi)部沒有自由電荷，因此電勢φ都滿足拉普拉斯方程，以φ1代表球外區(qū)域的電勢，φ2代表球內(nèi)的電勢，由式(11)得出的兩區(qū)域的通解為

(12)

(13)

式中an,bn,cn,dn是待定常數(shù).

邊界條件有：

(1)無窮遠(yuǎn)處，E→E0，均勻電場中空間任一點(diǎn)P的電勢φ(P)=φ0-E0·x，由于均勻電場可以看作由無窮大平行板電容器產(chǎn)生，其電荷分布不在有限區(qū)域內(nèi)，因此不能選φ(∞)=0.若選φ0=0，則有φ(P)=-E0·x，可得

φ1→-E0rcosθ=-E0rP1(cosθ)

(14)

所以有

a1=-E0an=0(n≠1)

(15)

(2)r=0處，φ2應(yīng)為有限值，因此

dn=0

(16)

(3)在介質(zhì)球面上(r=a)

(17)

把式(12)和(13)代入得

(18)

比較P1的系數(shù)得

(19)

由式(19)解出

(20)

比較式(18)其他Pn項(xiàng)的系數(shù)可解出

bn=cn=0n≠1

(21)

所有常數(shù)已經(jīng)定出，因此本問題的解為

(22)

(23)

為了簡化運(yùn)算，用Mathematica來求解本題，并畫出介質(zhì)球附近的電場線和等位面圖.所需代碼如下所示：

Clear["Global′*"]

v1=Sum[(An*R^n+Bn/R^(n+1))*Pn,{n,0,m})]

v2=Sum[(Cn*R^n+Dn/R^(n+1))*Pn,{n,0,m})]

v1=Sum[(An*R^n+Bn/R^(n+1))*Pn,{n,0,1})]/.{P0→1,P1→Cos[θ]}

v2=Sum[(Cn*R^n+Dn/R^(n+1))*Pn,{n,0,1})]/.{P0→1,P1→Cos[θ]}

(v1/.{B0→0,B0→0})==-E0*R*Cos[θ]

A0=0;A1=-E0;

D0=0;D1=0;

(v1/.R→a)==(v2/.R→a)

(ε0D[v1,R]/.R→a)==(εD[v2,R]/.R→a)

B0=0;

C0=0;

sol=Solve[{eq1,eq3},{B1,C1}]

v11=v1/.sol//First

v22=v2/.sol//First

E1=Grad[-v11,{R,θ,φ},"Spherical"]//

Simplify

E2=Grad[-v22,{R,θ,φ},"Spherical"]//

Simplify

E2=Grad[-v22,{x,y,z},"Cartesian"]

?VectorFieldPlots′;

GradientFieldPlot[-v3,{z,-2,2},{x,-2,2},PlotPoints→20]

Plot3D[v3,{z,-2,2},{x,-2,2},PlotPoints→30,LabelStyle→(FontSize→20),BoxRatios→{1,1,1}]

Mathematica繪制出的電場線圖和等電位圖如圖2所示.

(a)電場線圖

接下來對這題做一些極限情況下的討論：

(1)第一種極限情況：當(dāng)ε=ε0時，介質(zhì)球與環(huán)境的介電常數(shù)一樣，其對外場的響應(yīng)消失，把其代入式(22)、(23)，此時φ1=φ2=-E0rcosθ空間的電場就是均勻電場.用Mathematica來分析此時的電場，繪制電場線圖的主要代碼如下所示，繪制的電場線圖如圖3所示.

Clear["Global′*"]

v11=-E0RCos[θ]

v22=-E0RCos[θ]

Grad[-v11,{R,θ,φ},"Spherical"]//Simplify

E2=Grad[-v22,{R,θ,φ},"Spherical"]//

Simplify

v22=-E0x;

E2=Grad[-v22,{x,y,z},"Cartesian"]

?VectorFieldPlots′;

GradientFieldPlot[-v3,{z,-2,2},{x,-2,2},PlotPoints→20]

圖3 極限情況1下介質(zhì)球附近的電場線圖

結(jié)合φ1=φ2=-E0rcosθ和圖3進(jìn)一步驗(yàn)證了此種情況下的電場為勻強(qiáng)電場.

(24)

介質(zhì)球的總電偶極矩為

(25)

式(22)φ1中的第二項(xiàng)正是這個電偶極矩所產(chǎn)生的電勢

(26)

由于此時的外電場相當(dāng)于偶極子產(chǎn)生的電場，所以用Mathematica來繪制偶極子的電場來近似這種情況下的電場，Mathematica的代碼如下所示：

Clear["Global′*"];

V=q/(4*Pi*Subscript[ε,0])*(1/r1-1/r2)

cond={r1→R-l/2*Cos[θ],r2→R+l/2*

Cos[θ]}

V=V/.cond//SimplifySeries[V,{l,0,6}]//

Normal

Solve[V1==c,R]

Ee=Grad[-V1,{R,θ,ω},"Spherical"]

D[R[θ],θ]/R[θ]==Ee[[1]]/Ee[[2]]

DSolve[%,R[θ],θ]

g2=PolarPlot[Evaluate[Table

Identity];

Show[{g1,g2},DisplayFunction→$Display

Function,Epilog→{PointSize[0.03],Hue[2],

Point[{0,-2.5}],PointSize[0.03],Hue[.3],

Point[{0,2.5}]}]

V1[θ_,R_]:=cos[θ]/R2

RevolutionPlot3D[V1[θ,R],{R,0.01,3},{θ,0,2π},PlotRange→{-3,3},PlotPlints→30,

BoxRatios→{1,1,1}]

Mathematica繪制出的第二種極限情況下也就是電偶極子的電場線圖和等位面圖如圖4所示.

從圖4(a)中可以看出電偶極子的整個電場和電勢是一個數(shù)值分布對稱的區(qū)域，整個圖形是以電偶極子的軸線對稱的.再看圖4(b)，發(fā)現(xiàn)在正電荷附近電偶極子的電勢比較高；而在負(fù)電荷附近，電偶極子的電勢比較低，且在電偶極子的中垂線上電勢為零.

(a)電場線圖

(b)等位面圖

(3)當(dāng)ε→-∞時，介質(zhì)球內(nèi)的場為E內(nèi)→0，其效果相當(dāng)于一個導(dǎo)體球.而此時，p=4πε0a3E0，問題就變成了求解導(dǎo)體球附近的電場，所以，我們可以用導(dǎo)體球的情況來近似這種情況下的電場分布，求導(dǎo)體球附近的電場的Mathematica代碼如下所示：

Clear["Global′*"]

v=Sum[(An*R^n+Bn/R^(n+1))*Pn,{n,0,∞}]

v=Sum[(An*R^n+Bn/R^(n+1))*Pn,{n,0,1}]/.{P0→1,P1→Cos[θ]}

(v/.{B0→0,B1→0})==-E0*R*Cos[θ]

cond1={A1→-E0,A0→0};

(v/.cond1/.R→a)==0

cond2={B1→a^3E0,B0→0};

v=v/.cond1/.cond2

GradientFieldPlot[-v1[z,x],{z,-2.,2.},{x,-2.2.},PlotPoints→20,Epilog→{Hue.

[.6],Disk[{0,0},1]}]

Mathematica繪制出的第三種極限情況下也就是導(dǎo)體球附近的電場線如圖5所示.

圖5 導(dǎo)體球附近的電場線圖

從圖中我們發(fā)現(xiàn)，導(dǎo)體球外的電場分布規(guī)律和介質(zhì)球外的電場分布規(guī)律是一樣的，但是在這種情況下，導(dǎo)體球內(nèi)部是沒有電場存在的.因?yàn)楫?dāng)導(dǎo)體置于電場中時，導(dǎo)體內(nèi)的電荷會重新分布，形成電場來抵消外電場，最終就會導(dǎo)致導(dǎo)體內(nèi)的總電場為零.

4 結(jié)論

本文基于Mathematica討論了分離變量法在靜電場的應(yīng)用，通過用分離變量法求解介質(zhì)球附近的電場以及3種極限情況下電場分布的實(shí)例分析展示了Mathematica具有的強(qiáng)大的數(shù)值計算、符號計算和圖形繪制的功能，研究探討了如何通過Mathematica軟件輔助靜電場的學(xué)習(xí).學(xué)生通過學(xué)習(xí)Mathematica軟件，能夠很好地鍛煉動手能力和自主學(xué)習(xí)能力，他們在學(xué)習(xí)過程中可以將專業(yè)課的復(fù)雜數(shù)學(xué)問題交給Mathematica來進(jìn)行，而自己只需掌握基本的數(shù)學(xué)原理，了解相關(guān)知識，配合Mathematica圖形顯示，就能達(dá)到更充分更深層次理解內(nèi)容本質(zhì)的目的，尤其是西藏地區(qū)來自農(nóng)牧區(qū)的學(xué)生，數(shù)理基礎(chǔ)很差，Mathematica的引入解決了這些學(xué)生由于薄弱的數(shù)學(xué)基礎(chǔ)和抽象的專業(yè)概念所引起的在專業(yè)課學(xué)習(xí)上的困難，對這些學(xué)生具有重大意義.