格林函数法求解稳定场问题1 格林函数法求解稳定场问题(Green ’s Function) Green ’s Function, 又名源函数，或影响函数，是数学物理中的一个重要概念。
从物理上看，一个数学物理方程表示一种特定的场和产生这种场的源之间关系：Heat Eq.：()2222 ,ua u f r t t∂-∇=∂ 表示温度场u 与热源(),f r t 之间关系 Poission ’s Eq.： ()20u f r ρε∇=-=-表示静电场u 与电荷分布()f r 之间的关系场可以由一个连续的体分布源、面分布源或线分布源产生，也可以由一个点源产生。
但是，最重要的是连续分布源所产生的场，可以由无限多个电源在同样空间所产生的场线性叠加得到。
例如，在有限体内连续分布电荷在无界区域中产生的电势：()''04r dV r rρφπεΩ=-⎰这就是把连续分布电荷体产生的电势用点电荷产生的电势叠加表示。
或者说，知道了一个点源的场，就可以通过叠加的方法算出任意源的场。
所以，研究点源及其所产生场之间的关系十分重要。
这里就引入Green ’s Functions 的概念。
Green ’s Functions ：代表一个点源所产生的场。
普遍而准确地说，格林函数是一个点源在一定的边界条件和初始条件下所产生的场。
所以，我们需要在特定的边值问题中来讨论 Green ’s Functions.下面，我们先给出Green ’s Functions 的意义，再介绍如何在几个典型区域求出格林函数，并证明格林函数的对称性，最后用格林函数法求解泊松方程的边值问题。
实际上，只限于讨论泊松方程的第一类边值问题所对应的 Green ’s Functions 。
2 泊松方程的格林函数静电场中常遇到的泊松方程的边值问题：()()()()()201 f s u r r u r u r r nρεαβϕ⎧∇=-⎪⎪⎨∂⎡⎤⎪+=⎢⎥⎪∂⎣⎦⎩ 这里讨论的是静电场()u r ， ()f r ρ代表自由电荷密度。
格林函数(),G r r '：位于'r 的单位正电荷在r 处所激发的满足齐次边界条件的电势。
三维 Green ’s Functions 定解问题为：① ()()()()2301,,0S G r r r r G r r G r n δεαβ⎧''∇=--⎪⎪⎨'∂⎡⎤⎪+=⎢⎥⎪∂⎣⎦⎩这里()3r r δ'-表述了单位正电荷的体密度。
注意：对于第二类齐次边界条件且对于有限的研究区域，这个定解问题无解。
这是因为，虽然方程说明V 内有单位正电荷存在，而边界条件(),0SG r r n'∂=∂ 说明点源产生的场在边界S 上电场的法向分量(),n G r r E n'∂=-∂处处为零，说明边界条件与方程不相容。
另外，可以对方程作积分()()21,VG r r dv r r dv δε''∇=--⎰⎰这时要包含r '点，用高斯定理得1SG ds n ε∂=-∂⎰这就矛盾了！！！ 注： 高斯定理 VSAdv n Ads ∇⋅=⋅⎰⎰这时引入广义格林函数② ()()()2301,,0S G r r r r C G r r nδε⎧''∇=--+⎪⎪⎨'∂⎪=⎪∂⎩其中C 为常数，还要增加一个条件，以保证解的唯一性。
求解上面方程组①或②，可得在给定区域V 的泊松方程的各类边值问题的格林函数。
3 镜像法求G . F.用Green ’s Functions 去求解数理方程的定解问题，首先要求出相同边界、同类边值问题的Green ’s Function.3.1 镜像法的基本概念很多物理问题没有一个普遍奏效的解法，人们会发展许多方法，但往往每一种方法只能解决一部分问题。
人们熟知的一种办法是所谓“猜解”，即“尝试解”，这要有所谓的“唯一性定理”做保证。
唯一性定理：某些物理问题（如静电场边值问题）存在唯一解。
可以通过并不唯一的方法找到这个唯一解，这样就意味着解决物理方法上的多样性和灵活性。
静电镜像法是一种特殊的猜解方法，其基本思想：利用点电荷模拟边界面上的感应电荷或极化电荷。
可用于镜像法解决的问题包括：在点或线电荷与导体（或介质）共同存在的系统中，空间任一点的场是由点（或线）电荷与界面上感应（或极化）电荷共同产生的，而感应（或极化）电荷事先并不知道。
通过分析边界条件可以找到一个（或多个）像电荷来等效地代替导体面（或介质面）上的感应（或极化）电荷，从而把点（或线）电荷与界面上感应（或极化）电荷在待求区域产生场的求解问题转化为真实点电荷和虚像电荷在待求区域所产生场的简单叠加。
镜像法求边值问题的一般步骤为（以静电场为例）：1） 列出定解问题：电势在待求区域所满足的微分方程和边界条件； 2） 根据边界条件分析镜像电荷的个数和位置； 3） 写出电势分布的形式表达式（尝试解）；4） 把边界条件带入形式表达式以确定像电荷的量值和位置； 5） 把已求出的像电荷带入形式解以得到真实的电势分布； 6） 根据题意要求可由电势求场强、电荷分布及受力等问题。
静电镜像法分为：反射镜像法：平面镜法 球面镜法半透镜法：平面镜法球面镜法3．2 无界空间定解问题()()()2301,,0r G r r r r G r r δε→∞⎧''∇=--⎪⎨⎪'=⎩ 对应物理问题：单位正电荷1q =置于r '，求空间任一点(),,r x y z =处的电势(),?G r r '=库仑定律给出的解――无界区域的Green ’s Function ：()01,4G r r r r πε'='-=又叫基本解。
3．3 上半空间定解问题()()()()23001,, 0;,0, ,0z r G r r r r z G r r G r r δε=→∞⎧''∇=-->⎪⎨⎪''==⎩这里实际上可以给出满足第一类边界条件的 G . F. of the first kind. 物理问题：在0z = 处，有一无限大接地金属板，在r '处有一单位正电荷q ，求金属板上方任一点r 处的电势分布(),?G r r '=q '=1q =σrr 'r +r -1rr '(),G r r ''镜像法的基本思想用在这里：当电荷q 置于导体板的上方时，由于静电感应，板上出现异号电荷，空间电场是由电荷q 及感应电荷共同激发的，即01G G G =+。
而且满足静电平衡条件时，电力线垂直于导体面。
格林等效层定理：带电导体面上的电荷分布在导体外产生的电势，可以用导体面内的一定的等效电荷分布来代替。
我们通过电场分布分析，引进像电荷――假想电荷――来代替感应电荷作用。
在这里，我们在电荷q 相对于0z =平面的镜像位置引进1q '=-，那么q '和q 激发电场与q 和真实感应电荷激发的电场相同。
这里q '要满足q '和q 共同在导体面上产生的电势为零。
这样引进的像电荷和原电荷一起产生的场，就是要求的由原电荷和感应电荷产生的场，这样我们只需要求出像电荷的位置和大小。
像电荷的正确引进要符合：① 像电荷用在求解区域之外引入，因为感生电荷在上半空间的场1G 处处满足Laplace ’s Eq. 210G ∇=, 即在上半平面内是无源的。
② 像电荷的电量q '和位置要满足边界条件：()0,0r G r r ='= 和(),0r G r r →∞'=。
Then, 1q =和1q '=-激发的电势是待求的格林函数。
()00011,4414G r r r r πεπεπε+-'=-⎡⎤=金属板上的面电荷密度00?z G z σε=∂⎛⎫=-= ⎪∂⎝⎭ 应能证明：金属板上总电荷 1dS σ=-⎰这说明金属板上总感应电荷等于像电荷。
这是因为接地的导体平面相当于一面镜子，而1q '=-则是1q =的像，1q '=-称像电荷。
3．4 球外空间这里还是考虑第一类G.F.函数的求解问题。
定解问题()()()()023001,, ;,0, ,0r r r G r r r r r r G r r G r r δε=→∞⎧''∇=-->⎪⎨⎪''==⎩对应物理问题：接地金属球外r '处，有一单位正电荷1q =，求球外空间任一点r 处的电势(),?G r r '=首先引进像电荷q '，要不违反泊松方程，也就是让q '产生的电势1G 满足 Laplace ’s Eq., 210G ∇=, q '必须在求解区域之外一球内，考虑到对称性，q '还必须在r '上，放在r ''处。
为了保证球面电势为零，即()0,0r r G r r ='= 成立，q '为负电荷。
q '＝？，r ''＝？ 应由边界条件定。
考虑球面上一点0P ，由边界条件得：()001,04r r r r q q G r r r r r r πε==⎛⎫''=+=⎪ ⎪'''--⎝⎭ 也就是00 q q r r r r '=-'''--要求1q =orr 'r ''0rαpq '0P0000 const r r r q r q r r r r '-'''⇒-===='''-注意，这里考虑了若有两个相似三角形'0OQ P OP Q ∆∆，必有'const PQ PQ=。
由此确定了像电荷的位置和电量200r r r r q qr ''=''=-'这样，q 和q '激发的电势就是 Green ’s Function()01,4q q G r r r r r r πε⎛⎫''=+ ⎪ ⎪'''--⎝⎭用球坐标表示：场点： (),,r r θϕ=，q 电荷所在位置：(),,r r θϕ''''=，像电荷q '所在位置：()20,,,,r r r r θϕθϕ⎛⎫''''''''''== ⎪'⎝⎭，（这里,θθϕϕ''''''==） 2r r r'-=+（余弦定理）2r r r ''-=+= （余弦定理） where()()()()cos sin cos sin cos sin sin sin sin cos cos sin sin cos cos cos cos sin sin cos cos cos r r rr xx yy zz rr rr αθϕθϕθϕθϕθθθθϕϕθθαθθϕϕθθ'''''⋅==++''''''=++''''=-+'''⇒=-+（()cosϕϕ'-――加法公式）在考虑1q =，0r q q r'=-'， 我们得 ()01,4G r r πε⎛⎫⎪ ⎪'= ⎝场强： (),E G r r '=-∇ 球面上电荷分布： 00r r r φσε=∂⎛⎫=- ⎪∂⎝⎭ 球面上总电荷：1ds σ=-⎰由于球面上感应电荷在球外的场与像电荷－q '的场等效，所以电荷q 受感应电荷的力为 ()2014qq i F r r πε'=-'''-4 Green ’s Function ’s 对称性 ()(),,G r r G r r ''=重要物理意义：p '点的点源，在一定边界条件下，在p 产生的场等于：在p 置同样强度点源，在相同边界条件下在p '产生的场。