第4章 格林函数在这一章里,我们介绍数学物理方程中另外一种常用的方法—格林函数法.从物理上看,一个数学物理方程是表示一种特定的“场”和产生这种场的“源”之间的关系.例如,热传导方程表示温度场和热源之间的关系,泊松方程表示静电场和电荷分布的关系,等等.这样,当源被分解成很多点源的叠加时,如果能设法知道点源产生的场,利用叠加原理,我们可以求出同样边界条件下任意源的场,这种求解数学物理方程的方法就叫格林函数法.而点源产生的场就叫做格林函数.4.1δ函数几何学中的点是没有大小的,它仅仅表示空间的一个位置,因此物理学中的质点、点电荷等点源无法用几何中的点来表示.那么,我们用数学语言如何描述这类具有实际背景的点源呢?考虑一根长为l 的直线,其上任一点的坐标⎦⎤⎢⎣⎡−∈2,2l l x .若总电量为Q 的电荷均匀分布在直线上,则直线上的电荷分布的线密度)(x ρ是⎪⎪⎩⎪⎪⎨⎧≤>=2,2,0)(l x lQ lx x ρ (4.1.1) 由定积分的性质可知x x Q d )(∫+∞∞−=ρ (4.1.2)若将上述线段无限缩小,或者说令0→l ,则我们得到了一个物理上常用的点源—点电荷.此时,电荷分布密度用)(0x ρ表示,同时式(4.1.1)变为⎩⎨⎧=∞≠=0,,0)(0x x x ρ (4.1.3) 而此时,电量仍为Q ,则式(4.1.2)仍然成立.为了理解上的方便,我们修改一下问题的叙述:去电量1=Q ,线段长度为ε2,则密度分布函数为⎪⎩⎪⎨⎧≤>=εεεδεx x x ,21,0)(且1d )(d )(===∫∫−+∞∞−εεεεδδx x x x Q由此可见)(x εδ是偶函数,则由积分第一中值定理可得)()(d )()(d )()(εξεξδξδεε<<−==∫∫+∞∞−+∞∞−f x x f x x f x当0→ε时,我们有了新的结果,我们将它定义为δ函数. 我们称符合下述2个条件的函数为δ函数⎩⎨⎧≠=∞=0,00,)(x x x δ (4.1.4)且∫+∞∞−=1d )(x x δ (4.1.5)由极限理论可知,)(x δ是偶函数.∫∫+∞∞−+∞∞−→→===)0(d )()(d )()(lim )(lim 00f x x f x x x f x f δδξεεε (4.1.6))(x δ不是通常意义下的函数,它用来描述集中分布这种常见而又特殊的一类现象的数学工具.δ函数不局限于描述点电荷的分布密度,它可以用来描述任意点量的密度.借助于δ函数,我们可以方便地描述各类点源的分布情况.如电量Q 的点电荷的分布函数为)()(0x Q x δρ=.例1 设有一条张紧静止的无穷长的细弦,其线密度为1=ρ若在0=x 点,在很短的时间内,用大小为F 的力敲一下,使获得的冲量1=∆⋅t F .问弦上的初始速度v 是怎样的?解 若0≠x ,由于时间非常短,扰动尚未传动,所以0=v ;而在0=x 上有∞=v .此外,由于敲打前弦是静止的,所以弦上的动量是1=∆⋅t F ,即∫∫+∞∞−+∞∞−==⋅1d )()(d x x v x v x ρ故初速度)()(x x v δ=.例2 设有一根温度为C 0o度的导热杆,其线密度为ρ,比热为c ,现用火焰在0=x 处以极短的时间烤一下,传给杆的热量为Q ,请分析一下开始一瞬间杆上的温度)(x T 的分布?解 在刚开始一瞬间,我们有⎩⎨⎧=∞≠=0,,0)(x x x T且∫+∞∞−=Q x x T c d )(ρ所以有)()(x c Qx T δρ=通过以上两个例题,我们对)(x δ有了进一步的认识.如果将坐标平移0x ,即集中量出现在点0x x =处,则有⎩⎨⎧=∞≠=−000,,0)(x x x x x x δ且∫+∞∞−=−1d )(0x x x δ这样，我们可以得到δ函数的一个重要性质)(d )(00x f x x x ∫+∞∞−=−δ或者说⎩⎨⎧><<<=−∫bx a x bx a x x x ba0000,0,1d )(或δ⎩⎨⎧><<<=−∫b x a x bx a x f x x x x f b a00000,0),(d )()(或δ4.2 无界域中的格林函数在第1章中,我们推导出了静电场的电势分布u 满足泊松方程ρε1222222−=∂∂+∂∂+∂∂=∆zu y u x u u (4.2.1)式中,ρ是电荷密度,所占区域为Ω,0r 是Ω中任意一个点.如果不考虑其他因素的影响,对于无界空间中的电势u ,可以利用定积分中的微元法的思想求出来.有库仑定律知,位于0r 点的一个正的单位电荷,在无界空间中点r 处产生的电势是041),(r r r r G −=π (4.2.2)则以0r 为中心的小体积Ωd 在r 处产生的电势为Ω=d )(),(d 00r r r G u ρ因此,在r 处产生的电势为∫∫ΩΩΩ−==d 4)(d )(00r r r u r u πρ为了表述上的方便, 0r 处的体积微元Ωd 以后用0d r 表示,则有∫Ω−=000d 4)()(r r r r r u πρ这样,我们没有直接求解方程,而是通过寻找微元,利用积分的方式求出了方程的解.而点源产生的电势),(0r r G 称为泊松方程式(4.2.1)在无界空间中的格林函数,利用它,我们求出了泊松方程在无界空间的解.无界空间中的格林函数又叫做方程的基本解,因此式(4.2.2)又称为泊松方程的基本解.有时也称它为相应的齐次方程(即拉普拉斯方程)的基本解,记为).,(00r r G基本解式(4.2.2)是密度为0ρ的点源在空间产生的电势,因此它在空间除了0r r =点以外,满足方程001ρε−=∆G而在0r r =点有奇异性.由于格林函数是点源函数,因此在空间某一点有奇异性. 在一般的数学物理方程中,我们需要考虑的是满足一定边界条件和初始条件的解,因此相应的格林函数就比刚才所提到的要复杂.在这种情况下,一个点源所产生的场,同时要受到边界条件及初始条件的影响,而这些影响的本身也是待定的. 例如,在一个接地的导体空腔内的点0P 处放置一个正的单位点电荷(如图4-1),则在点P 处的电势不仅是点电荷本身所产生的场41r r −π,并且还要加上这个点电荷在导体内壁上感应电荷所产生的场.而感应电荷在导体内壁上的分布是未知的,我们只知道在边界上电势为零(接地).因此,在一般情况下,格林函数是一个点源在一定的边界条件和(或)初始条件下所产生的场.通过格林函数,我们可以求得任意分布的源所产生的场.4.3 格林公式 有界域上的格林函数为了进一步探讨利用格林公式函数求解数学物理方程,我们先来推出一个重要工具—格林公式,它是曲面积分中高斯公式的直接推论.设Ω是以足够光滑的曲面Γ为边界的有界域,),,(),,,(),,,(z y x R z y x Q z y x P 在Γ+Ω上是连续的,在Ω内具有一阶连续偏导数,则有如下的高斯公式∫∫∑++=Ω⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂ΩS z n R y n Q x n P z R y Q x P d )],cos(),cos(),cos([d (4.3.1) 式中,Ωd 是体积元素;n 是曲面Γ的外法向量;S d 是Γ上的面积元素.设函数),,(),,,(z y x v z y x u 在Γ+Ω上一阶偏导数连续,在Ω内二阶偏导数连续,则在式(4.3.1)中,令z vR yv u Q x v uP ∂∂=∂∂=∂∂=,,则有∫∫∫∫∫∫ΓΩΩΩΩΩ∂∂=Ω⋅+Ω∆=Ω⎟⎟⎠⎞⎜⎜⎝⎛∂∂⋅∂∂+∂∂⋅∂∂+∂∂⋅∂∂+Ω∆=Ω⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂S nvuv u v u z v z u y v y u x v x u v u z R y Q x P d d grad grad d )(d d )(d 或表示为Ω⋅−∂∂=Ω∆∫∫∫ΩΓΩd grad grad d d )(v u S nvuv u (4.3.2)式(4.3.2)称为格林第一公式.在式(4.3.2)中,交换v u ,的位置，则有Ω⋅−∂∂=Ω∆∫∫∫ΩΓΩd grad grad d d )(v u S nuvu v (4.3.3)式(4.3.2)减式(4.3.3)得∫∫ΓΩ⎟⎠⎞⎜⎝⎛∂∂−∂∂=Ω∆−∆S n u v n vu u v v u d d ][ (4.3.4) 式(4.3.4)称为格林第二公式.下面,我们以泊松方程第一类边值问题为例,进一步阐明格林函数的概念.⎪⎩⎪⎨⎧=−=∆Γ)6.3.4()5.3.4(1f u u ε式中, f 是在区域Ω上的边界Γ上给定的函数.在介绍格林函数之前,我们要引进空间的δ函数来表示点源的密度分布,有)()()()(0000z z y y x x r r −−−=−δδδδ⎩⎨⎧=∞≠=−000,,0)(r r r r r r δ )),,((1d )(00000Ω∈=−∫Ωz y x r r r r δ∫Ω=−)(d )()(00r f r r f r r δ用),(0r r G 表示位于0r 点的单位强度的正点源在第一类边界条件下产生的场,则),(0r r G 作为r 的函数满足⎪⎩⎪⎨⎧=−−=∆Γ)8.3.4(0)7.3.4()(1),(00G r r r r G δε以),(0r r G 乘式(4.3.5),)(r u 乘式(4.3.7),二式相减后在Ω上对r 积分,以r d 表示r 点处的体积微元,有∫∫∫ΩΩΩ−+−=∆−∆r r r r u r G r G u u G d )()(1d 1d )(0δερε利用格林第二公式及δ函数的性质,有)9.3.4(d ),()(d )(),(d ),()(d )(),(d ),()()(),(d )(),()(00000000∫∫∫∫∫∫ΓΩΓΩΓΩ∂∂−=∂∂−=⎦⎤⎢⎣⎡∂∂−∂∂+=S nr r G r f r r r r G S nr r G r u r r r r G S n r r G r u n r u r r G r r r r G r u ερερερ但这个表达式中所表示的意义与我们的初衷相矛盾.),(0r r G 表示的是位于0r 点的点源在r 点产生的场.但我们能证明),(),(00r r G r r G =,这样,式(4.3.9)可以改写成)10.3.4(d ),()(d )(),(d ),()(d )(),()(0000000000∫∫∫∫ΓΩΓΩ∂∂−=∂∂−=S nr r G r f r r r r G Snr r G r f r r r r G r u ερερ这样,式(4.3.1)的物理诠释就很清楚了:右方第一个体积分代表在区域Ω中体分布源)(0r ρ在r 点产生的场的总和,第二个面积分则表示了在边界上的源所产生的场. 下面我们来证明),(),(00r r G r r G =,由式(4.3.7)及式(4.3.8),我们有⎪⎩⎪⎨⎧=−−=∆Γ)12.3.4(0),()11.3.4()(1),(111r r G r r r r G δε⎪⎩⎪⎨⎧=−−=∆Γ)14.3.4(0),()13.3.4()(1),(222r r G r r r r G δε×),(2r r G 是式(4.3.11)—×),(1r r G 式(4.3.13),有)(),()(),()],(),(),(),([21122112r r r r G r r r r G r r G r r G r r G r r G −−−==∆−∆δδε两侧同时对r 积分,有∫∫ΩΩ−−−=∆−∆rr r r r G r r r r G r r r G r r G r r G r r G d )(),()(),(d )],(),(),(),([21122112δδε根据格林公式第二公式及δ函数的性质,有),(),(d ),(),(),(),(12212112r r G r r G S n r r G r r G n r r G r r G −=⎥⎦⎤⎢⎣⎡∂∂−∂∂∫Γε 则根据式(4.3.12)及式(4.3.14),有0),(),(),(),(2112=∂∂−∂∂Γnr r G r r G nr r G r r G 所以),(),(1221r r G r r G =这种性质在物理学中称为倒易性,如图4-2所示,即位于1r 点的点源,在一定的边界情况下,在2r 点产生的场等于位于2r 点的同样强度的点源,在相同的边界情况下在1r 点产生的场.我们称这种现象为格林函数的对称性.应当说明，在得式(4.3.9)时,我们利用格林公式把重积分化为曲面积分时,这要求G ∆(及u ∆)在积分区域Ω内连续为前提,由式(4.3.7)可明显看到G ∆不连续,这样的推导请参阅谷超豪等著《数学物理方程》(第二版).4.4 格林函数的应用在第1章里,我们从无源静电场的电位分布及稳恒温度场的温度分布推出了三维拉普拉斯方程0222222=∂∂+∂∂+∂∂=∆zu y u x u u作为描述稳定或平衡等状态的方程,它与初始状态无关,因而不能提初始条件.对于边界条件,常见的是如下两种现象.第一边值问题 在空间),,(z y x 中某一区域Ω的边界Γ上给定了连续函数f ,要找这样的函数),,(z y x u ,它在闭区域Γ+Ω上连续,且满足⎩⎨⎧==∆Γf u u 0第一边值问题也称为狄利克莱(Dirichlet)问题,或简称为狄氏问题.拉普拉斯方程的连续解,即具有二阶连续偏导数,并且满足拉普拉斯方程的连续函数,称为调和函数.因此, 狄氏问题也可以这样叙述:在区域Ω内找一个调和函数,它在边界Γ上的值是已知的.第二边值问题 在空间在空间),,(z y x 中某一区域Ω的边界Γ上给定了连续函数f ,要找这样的函数),,(z y x u ,它在闭区域Γ+Ω上连续,且满足⎪⎩⎪⎨⎧=∂∂=∆Γf nuu 0 式中,n 是曲面Γ的外法向矢量.第二边值问题也称为诺依曼(Neumann)问题.以上两个边值问题都是在边界Γ上给定某些条件,在区域内部求解拉普拉斯方程,这样的问题称为内问题.在应用中,我们还会碰上另一类现象,如确定某物体外部的稳恒温度场时,就归结为在区域Ω的外部求调和函数u 使之满足边界条件f u =Γ,这里Γ是区域Ω的边界,f 表示物体表面的温度分布.这样的问题称为拉普拉斯方程的外问题. 限于篇幅,本书仅讨论如何利用格林函数求解狄利克莱问题⎩⎨⎧==∆Γ)2.4.4()1.4.4(0fu u至于其他的问题,求解的思考方法是想像的,可查阅相关的书籍.由式(4.4.1)知源的分布密度函数0=ρ,所以上节给出的求解公式就变为∫Τ∂∂−=S nr r G r f r u d ),()()(00 (r 在曲面Γ上) 或∫Τ∂∂−=S nr r G r f r u d ),()()(00 (0r 在曲面Γ上) (4.4.3) 此处介电常数1=ε. 这样,对一个由曲面Γ围成的区域Ω来说,只要求出了格林函数),(0r r G ，则这个区域内狄氏问题的解就可以由式(4.4.3)求出.实际上,求解边值问题式(4.3.7)—式(4.3.8)是很困难的,因此有必要对格林函数),(0r r G 作进一步的剖析.在本章中,我们定义了方程的基本解),(00r r G ,它满足方程式(4.3.7))(),(000r r r r G −−=∆δ但不满足边界条件式(4.3.8).于是我们设)(),(),(000r V r r G r r G +=代入式(4.3.7)及边界条件式(4.3.8),则有⎩⎨⎧−==∆ΓΓ00G V V这样,只要找到满足边界条件ΓΓ−=0G V的调和函数V ,那么就可以由基本解得到格林函数),(0r r G .事实上,当区域的边界具有特殊的对称性时,格林函数是用镜像法(静电源像法)求得的.所谓静像法,就是在区域Ω外找出点0M 关于边界Γ的像点(对称点)1M ,然后在1M 上放置适当的负电荷,由它所产生的负电位与点0M 处单位电荷产生的电位在曲面Γ上相互抵消.此时,放置在0M ,1M 两点处的电荷所形成的电场在Ω内的电位就是所要求的格林函数.下面,我们以寻求半空间、球域的格林函数为例来说明镜像法的具体应用.例3 求解上半空间0>z 内的狄利克莱问题⎪⎩⎪⎨⎧+∞<<−∞=>=∂∂+∂∂+∂∂=)5.4.4(),(0)4.4.4()0(00222222y x u z z uy u xu z解 先求出格林函数),(0r r G .为此在上半空间0>z 中任意一点),,(0000z y x r 处置一单位正电荷,在点0x 关于平面0=z 的对称点),,(0001z y x r −处置一单位负电荷,如图4-3所示.由它们所形成的静电场的电势在平面0=z 上恰好为零.因此上半空间的格林函数为⎟⎟⎠⎞⎜⎜⎝⎛−−−=1001141),(r r r r r r G π(4.4.6)为了利用式(4.4.3)求解问题式(4.4.4),式(4.4.5)需要计算边界曲面上的nG∂∂值.由于在平面0=z 上的外法线方向是Oz 轴的负向,所以)7.4.4(])()[(210])()()[(])()()[(4123220200232020200232020200000z y y x x z z z z y y x x z z z z y y x x z z z G nG z z +−+−−==⎪⎭⎪⎬⎫++++++−⎪⎩⎪⎨⎧−+−+−−=∂∂−=∂∂=ππ则定解问题式(4.4.4),式(4.4.5)的解为∫∫+∞∞−+∞∞−+−+−=ηξηξηξπd d ])()[(),(21),,(23222z y x zf z y x u (4.4.8)用同样的方法，我们可以求出球域上的格林函数,并给出球域内的狄利克莱问题的解.设有一球心在原点,半径为R 的球面Γ.在球内任取一点),,(0000z y x r ，在0Or 的延长线上截取线段1Or ,令00ρ=Or ,11ρ=Or ,使210R =⋅ρρ,这样的点1r 称为点0r 关于球面Γ的反演点(或对称点),如图4-4所示.我们在点0r 处放置一单位正电荷,在点1r 处放置一q 单位的负电荷,通过选择恰当的q 值,使得这两个点电荷所产生的电势在球面Γ为零.即P r qP r 10441ππ=或 Pr P r q 01=式中,P 为球面Γ上任意一点.由于三角形△P Or 1与△P Or 0在点O 处有公共角,且夹这个角的两条边成比例1ρρRR=,因此这两个三角形相似.于是得到01ρRP r P r =因此ρRq =即只要在点1r 处放ρR单位的负电荷,则由0r 及1r 处点源产生的电势在球面上为零,这样,球域内的格林函数为⎟⎟⎠⎞⎜⎜⎝⎛−−−=10001141),(r r R r r r r G ρπ(4.4.9) 式中,r 为球域内任意一点,记0ρ=Or .下面,我们利用格林函数来求解球域内的狄利克莱问题⎩⎨⎧==∆Γf u u 0Ω∈),,(z y x 由式(4.3.9)得(介电常数)1=εS nr r G r f r u d ),()()(00∫Γ∂∂−=因此,我们要计算Γ∂∂n G,由 γρρρρcos 21102200−+=−r rγρρρρcos 21112211−+=−r r012ρρ⋅=R式中,γ是向量0Or 与Or 的夹角.所以⎥⎥⎦⎤⎢⎢⎣⎡+−−−+=40222002200cos 21cos 2141),(R M M G γρρρρργρρρρπ在球面Γ上 2302022022340222002202302020)cos 2(41)cos 2()cos ()cos 2(cos 41γρρρπργρρρργρργρρρργρρπρρR R R R RR R R R RG G −+−−==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−−+−−==∂∂=ΓΓ∂∂ 所以狄氏问题的解为S f R R R R r u d )cos 2(41)(23022220∫∫Γ−+−=γρρρπ (4.4.10)为了方便解释物理现象,我们也可以利用格林函数的倒易性,求出球内任一点r 处的电势)(r u .在球面上应用球坐标系,上式变为∫∫−+−=ππθϕϕγρρρθϕπθϕρ202302222000d d sin )cos 2(),,(4),(R R R R f Ru (4.4.11)式中, ),(000θϕρ是点0r 的坐标;),,(θϕR 是球面Γ上点P 的坐标;γcos 是向量0Or 与OP 的余弦.因为向量0Or 与Or 的方向余弦分别是)cos ,sin sin ,sin (cos )cos ,sin sin ,sin (cos 00000ϕϕθϕθϕϕθϕθi所以可得)cos(sin sin cos cos )cos cos sin (sin sin sin cos cos cos 0000000θθϕϕϕϕθϕθθϕϕϕϕγ−+=++=式(4.4.10)及式(4.4.11)称为球的泊松公式.例4 设有一半径为R 的均匀球,球心在坐标原点,上半球面的温度保持为C o0,下半球面的温度保持为C o2,求:(1) 球内温度的稳定分布; (2) 球内z 轴上温度的分布; (3) 球心的温度.解 这个问题的数学描述为⎪⎪⎩⎪⎪⎨⎧⎪⎪⎩⎪⎪⎨⎧<<<<=<=∆=πϕππϕρρ2,220,0)(0R u R u由泊松公式,球内任一点),(0θϕρ处的温度为∫∫∫∫−+−=−+−=ππππθϕϕγρρρπθϕϕγρρρθϕπθϕρ2023020220220023020222000d d sin )cos 2(2d d sin )cos 2(),,(4),(R R R R R R R R f Ru若只考虑z 轴上的温度,即00=ϕ(上半轴)或πϕ=1(下半轴), 可知:当00=ϕ时,ϕγcos cos =,则⎟⎟⎠⎞⎜⎜⎝⎛+−+−===⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−+−−=−+−=∫∫2020202210202022202230222200112)cos 2(d d sin )cos 2(2),0,(ρρρρπϕπϕϕρρρρθϕϕϕρρρρπθρπππR R R R R R R R R R R Ru当πϕ=0时ϕγcos cos −=,故⎟⎟⎠⎞⎜⎜⎝⎛+−−−=202002020011),,(ρρρρθπρR R R u 当00→ρ时,应用洛必达法则有1),,(lim )0,0,0(00000==→θϕρρu u即球心温度为C o1。
