2022年格林函数 2.pdf

1 2.4 格林函数法解的积分公式在第七章至第十一章中主要介绍用分离变数法求解各类定解问题，本章将介绍另一种常用的方法 格林函数方法。格林函数，又称点源影响函数，是数学物理中的一个重要概念。格林函数代表一个点源在一定的边界条件和（或）初始条件下所产生的场。知道了点源的场，就可以用迭加的方法计算出任意源所产生的场。一、泊松方程的格林函数法为了得到以格林函数表示的泊松方程解的积分表示式，需要用到格林公式，为此，我们首先介绍格林公式。设 u（r）和 v（r）在区域 T 及其边界上具有连续一阶导数，而在T 中具有连续二阶导数，应用矢量分析的高斯定理将曲面积分Sdvu化成体积积分.)(TTTvdVuvdVudVvuSdvu（12-1-1）这叫作 第一格林公式 。同理，又有.TTvdVuudVvSduv（12-1-2）（12-1-1）与（ 12-1-2）两式相减，得,)()(TdVuvvuSduvvu亦即.)(TdVuvvudSnuvnvu（12-1-3）n表示沿边界的外法向求导数。（12-1-3）叫作 第二格林公式 。现在讨论带有一定边界条件的泊松方程的求解问题。泊松方程是)(),(Trrfu（12-1-4）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 16 页 - - - - - - - - - 2 第一、第二、第三类边界条件可统一地表为),(Munu（12-1-5）其中 （M）是区域边界上的给定函数。0，0 为第一类边界条件，0，0 是第二类边界条件，、 都不等于零是第三类边界条件。泊松方程与第一类边界条件构成的定解问题叫作第一边值问题或狄里希利问题，与第二类边界条件构成的定解问题叫作 第二边值问题或诺依曼问题，与第三类边界条件构成的定解问题叫作第三边值问题 。为了研究点源所产生的场，需要找一个能表示点源密度分布的函数。 5.3 中介绍的 函数正是描述一个单位正点量的密度分布函数。因此，若以v（r，r0）表示位于 r0点的单位强度的正点源在r 点产生的场，即 v（r，r0）应满足方程).(),(00rrrrv（12-1-6）现在， 我们利用格林公式导出泊松方程解的积分表示式。以 v （r， r0） 乘 （12-1-4） ，u（r）乘（ 12-1-6） ，相减，然后在区域T 中求积分，得.)()(0TTTdVrruvfdVdVvuuv（12-1-7）应用格林公式将上式左边的体积分化成面积分。但是，注意到在rr0点，v 具有 函数的奇异性，格林公式不能用。解决的办法是先从区域T 中挖去包含 r0的小体积，例如半径为的小球 K （图 12-1） ，的边界面为。对于剩下的体积，格林公式成立，.)(dSnvunuvdSnvunuvdVvuuvKT（12-1-8）把（12-1-8）代入挖去 K 的（12-1-7） ，并注意 rr0，故 （rr0）0，于是.KTvfdVdSnvunuvdSnvunuv（12-1-9）当10rr，方程（ 12-1-6）的解 v（r，r0） 位于点 r0而电量为0 的点电荷的静电场中的电势，即140rr。令 0，得O y z x T r0K图 12-1名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 16 页 - - - - - - - - - 3 （12-1-9 ）右边,TvfdV左边的0414102rrnudnudnudSnuv左边的).(141141022rudrrudSrrudSnvu（12-1-10）这样， （12-1-7）成为.),()()(),()(),()(0000dSnrrvrunrurrvdVrfrrvruT（12-1-11）（12-1-11）称为 泊松方程的基本积分公式。（12-1-11）将（ 12-1-4）的解 u 用区域 T 上的体积分及其边界上的面积分表示了出来。那么，能否用（12-1-11）来解决边值问题呢？我们看到， （12-1-11）中需要同时知道 u 及nu在边界上的值，但是，在第一边值问题中， 已知的只是 u 在边界 上的值；在第二边值问题中，已知的只是nu在边界上的值。在第三边值问题中，已知的是 u 和nu的一个线性关系在边界上的值，三类边界条件均未同时分别给出 u 和nu的边界上的值。因此，我们还不能直接利用（12-1-11）解决三类边值问题。其实，这里距离问题的解决已经很近了。原来，对于函数v（r，r0） ，我们还只考虑其满足方程（ 12-1-6） 。如果我们对v（r，r0）提出适当的边界条件，则上述困难就得以解决。对于第一边值问题， u 在边界 上的值是已知的函数（M） 。如果要求 v 满足齐次的第一类边界条件,0v（12-1-12）则（12-1-11）中含nu的一项等于零。从而不需要知道nu在边界上的值。满足方名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 16 页 - - - - - - - - - 4 程（12-1-6）及边界条件（ 12-1-12）的解称为 泊松方程第一边值问题的格林函数，用 G（r，r0）表示。这样，（12-1-11）式成为.),()()(),()(000dSnrrGrdVrfrrGruT（12-1-13）对于第三边值问题，令v 满足齐次的第三类边界条件，.0vnv（12-1-14）满足方程（ 12-1-6）及边界条件（ 12-1-14）的解称为 泊松方程第三类边值问题的格林函数 ，也用 G（r，r0）表示。以 G（r，r0）乘（ 12-1-5）式两边，得.GuGnuG又以 u 乘（12-1-14） ，并以 G 代替其中的 v，得.0uGnGu将这两式相减，得.GnGunuG将此式代入（ 12-1-11） ，得.)(),(1)(),()(000dSrrrGdVrfrrGruT（12-1-15）至于第二边值问题，表面看来，似乎可以按上述同样的办法来解决，即令G 为定解问题),(0rrG（12-1-16）0nG（12-1-17）的解，而由（ 12-1-11）得到.)(),()(),()(000dSrrrGdVrfrrGruT（12-1-18）可是，定解问题（ 12-1-16）（12-1-17）的解不存在。这在物理上是容易理解的：不妨把这个格林函数看作温度分布。泛定方程（12-1-16）右边的 函数表明在所围区域 T 中有一个点热源。边界条件（ 12-1-17）表明边界是绝热的。点热源不停地放也热量。而热量又不能经由边界散发出去，T 里的温度必然要不停地升高，其名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 16 页 - - - - - - - - - 5 分布不可能是稳定的。这就需要引入推广的格林函数 。对于三维空间，,1)()()(000TVzzyyxxG.0nG式中 VT 是 T 的体积。对于二维空间，,1)()(00TAyyxxG.0nG式中 AT 是 T 的面积，方程右边添加的项是均匀分布的热汇密度，这些热汇的总体恰好吸收了点热源所放出的热量，不多也不少。（12-1-13）和（12-1-15）的物理解释有一个困难。公式左边u的宗量 r0表明观测点在 r0，而右边积分中的f（r）表示源在 r，可是，格林函数 G（r，r0）所代表的是 r0的点源在 r 点产生的场。这个困难如何解决呢？原来，这个问题里的格林函数具有对称性 G（r，r0）G（r0，r） ，将（ 12-1-13）和（ 12-1-15）中的 r 和 r0对调，并利用格林函数的对称性， （12-1-13）成为,),()()(),()(0000000dSnrrGrdVrfrrGruT（12-1-19）这就是 第一边值问题解的积分表示式。 （12-1-15）成为,)(),(1)(),()(000000dSrrrGdVrfrrGruT（12-1-20）这就是 第三边值问题解的积分表示式。（12-1-19）和（12-1-20）的物理意义就很清楚了，右边第一个积分表示区域T 中分布的源 f（r0）在 r 点产生的场的总和。第二个积分则代表边界上的状况对r 点场的影响的总和。两项积分中的格林函数相同。这正说明泊松方程的格林函数是点源在一定的边界条件下所产生的场。现在来证明格林函数的对称性。 在 T 中任取两个定点 r1和 r2。 以这两点为中心，各作半径为 的球面1和 2。从 T 挖去 1和 2 所围的球 K1和 K2。在剩下的区域 TK1K2上，G（r，r1）和 G（r，r2）并无奇点。以uG（r，r1） ，vG（r，r2）代入格林公式（ 12-1-3）2121)(KKTdVuvvudSnuvnvu名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 16 页 - - - - - - - - - 6 由于 G（r，r1）和 G（r，r2）是调和函数，上式右边为零。又由于格林函数的边界条件，上式左边0。这样.021dSnuvnvudSnuvnvu令0，上式成为 0v（r1）u（r2）00，即 G（r1，r2）G（r2，r1） 。对于拉普拉斯方程， 即 （12-1-4） 式右边的 f （r） 0， 这时，我们只要令（12-1-19）和（12-1-20）两式右边的体积分值等于零，便可得到拉普拉斯方程第一边值问题的解0000),()()(dSnrrGrru（12-1-21）以及第三边值问题的解000)(),(1)(dSrrrGru（12-1-22）我们看到，借助格林公式，也可利用格林函数方法得到齐次方程定解问题的解。二、用电像法求格林函数（一）无界空间的格林函数基本解从 12.1 讨论可知，确定了 G， 就能利用积分表式求得泊松方程边值问题的解。虽然，求格林函数的问题本身也是边值问题，但这是特殊的边值问题，其求解比一般边值问题简单。特别是对于无界区域的情形，常常还可以得到有限形式的解。无界区域的格林函数称为相应方程的基本解 。我们将一个一般边值问题的格林函数G 分成两部分.10GGG（12-2-1）其中 G0是基本解。对于三维泊松方程，即G0满足).(00rrG（12-2-2）G1则满足相应的齐次方程（拉普拉斯方程）01G（12-2-3）及相应的边界条件。例如在第一边值问题中，0G，从而有名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 16 页 - - - - - - - - - 7 .)(001GGGG（12-2-4）拉普拉斯方程（ 12-2-3）的边值问题的求解是熟知的。至于方程（12-2-2） ，它描述的是点 r0的点源在无界空间产生的稳定场。以静电场为例，它描述在点r0电量为0的点电荷在无界空间中所产生电场的r 点的电势，即004/1rrG。现在再给出（ 12-2-2）的一种解法。先假设点源位于坐标原点，由于区域是无界的，点源产生的场应与方向无关，如果选取球坐标（r， ， ） ，则 G0只是 r 的函数，方程（ 12-2-2）变成一个常微分方程，当r0 时，G0满足拉普拉斯方程,01022drdGrdrdr（12-2-5）其解为.210CrCG（12-2-6）令无穷远处 G00，于是 C20。为了求出 C1，将方程（ 12-2-2）在包含 r00 的区域作体积分，这个区域可取为以r00 为球心，半径为 的小球 K ，其边界面为（参见图 12-1） ，.10KdVG利用（ 12-1-3） （令其中的 u1） ，将上式右边体积分化成面积分。120021004sinCddrrCrdSrGdVGK则411C，从而.141)(0rrG若电荷位于任意点r0，则.141),(000rrrrG（12-2-7）类似地，用平面极坐标可求得二维泊松方程的基本解.1ln21),(000rrrrG（12-2-8）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 16 页 - - - - - - - - - 8 （二）用电像法求格林函数让我们来考虑这样一个物理问题。设在一接地导体球内的M0（r 0）点放置一带电量为 0的点电荷。则球内电势满足泊松方程),(0rrG（12-2-9）边界条件是.0球面G（12-2-10）此处 G 便是泊松方程第一边值问题的格林函数。从电磁学知道，在接地导体球内放置电荷时，导体球面上将产生感应电荷。因此，球内电势应为球内电荷直接产生的电势与感应电荷所产生的电势之和。因此，我们可将G 写成两部分之和,10GGG（12-2-11）其中 G0是不考虑球面边界影响的电势，G1则是感应电荷引起的。由前面的讨论可知，G0满足),(00rrG（12-2-12）从而 G1满足01G（12-2-13）以及边界条件.)(001球面球面球面GGGG（12-2-14）这样， G0就是基本解，0004/1),(rrrrG。至于 G1则可从方程（ 12-2-13）及边界条件（ 12-1-14）用分离变数等方法求得。但这样得到的解往往是无穷级数。现在介绍另一种方法 电像法 ，用电像法可以得到有限形式的解。电像法的基本思想是用另一设想的等效点电荷来代替所有的感应电荷，于是可求得G1的类似于 G0的有限形式的解。显然，这一等效点电荷不能位于球内，因为感应电荷在球内的场满足（12-2-13） ，即球内是无源的。又根据对称性，这个等效电荷必位于OM0 的延长线上的某点 M1， 记等效电荷的电量为P M M0r0O M1图 12-2 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 16 页 - - - - - - - - - 9 q，其在空间任意点M（r）引起的电势是10114/),(rrqrrG。若将场点取在球面上的 P 点，如图 12-2 所示，则OPM0和 OM1P 具有公共角 POM1，如果按比例关系 r0aar1（a 为球的半径）选定 M1（这 M1必在球外）， 则 OPM0 跟OM1P 相似，从而球面上01rr0111rrr球面上.1a因此，若取00/raq，则球面上的总电势是. 01411411410011100rarrrrrrrrrarr正好满足边界条件（ 12-2-10） 。这个设想的位于M1点的等效点电荷称为M0点点电荷的电像。这样，球内任一点的总电势是.141141141141),(0202001000rrarrarrrrrarrrrG（12-2-15） 10.1 例 6 求出球外点电荷的电像（在球内） ，读者不妨把这两种情况中的电像加以对比。若 M0（r0）为圆内的一点，则圆内泊松方程第一边值问题的格林函数满足),(0rrG（12-2-16）.0圆周上G（12-2-17）这个问题也可用电像法求解，结果是,ln211ln211ln21),(0100rarrrrrrG（12-2-18）式中 a 为圆的半径。例 1 在球 ra 内求解拉普拉斯方程的第一边值问题).,(),(03fuaruar解前面已用电像法求得球的第一边值问题的格林函数名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 16 页 - - - - - - - - - 10 .141141),(1000rrrarrrrG把它代入第一边值问题的解的积分公式（12-1-13）就行了。为了把 G（r，r0）代入（ 12-1-19） ，还必须先算出nG。引用球坐标系，极点就取在球心。,cos21120020rrrrrr（12-2-19）其中是矢径 r 跟 r0之间的夹角，).cos(sinsincoscoscos000计算法向导数2/32002020020)cos2(coscos211rrrrrrrrrrrrrn分子里的 cos 可利用（ 12-2-19）消去，.221302202030220200rraarrrrrrrrrrrrn同理，.222130202023033002202202043122121010rrarrrararrrararrrarraarrrrarrran于是.412412413020230202023022020rrararrarrrarraarrrnG代入（ 12-1-13） ，得到球的第一边值问题的解的积分公式名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 16 页 - - - - - - - - - 11 0202/32002202020230202000sin)cos2(),(4sin41),(),(ddrararafaddarrarafru作代换),(),(000rr：0200002/322220000sin)cos2(),(4),(ddrararafaru这叫作球的 泊松积分 。例 2 在半空间 z0 内求解拉普拉斯方程的第一边值问题).,(),0(003yxfuzuz解先求格林函数 G（r，r0）.0),()()(00003zGzzyyxxG这相当于接地导体平面z0上方的电势，在点 M0（x，y，z）放置着电量为 0的点电荷。这电势可用电像法求得。设想在 M0的对称点 M1（x0，y0，z0） 放置电量为 0的点电荷，不难验证，在两个点电荷的电场中，平面z0 上的电势确实是零。在点M1的点电荷就是电像。格林函数.)()()(141)()()(141141141),(202020202020100zzyyxxzzyyxxrrrrrrG为了把 G （r，r0）代入第一边值问题的解的积分公式（12-1-13） ，需要先计算0znGz M0(x0, y0, z0)O x y M(x, y, z)图 12-3 M1(x0, y0, - z0)名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 16 页 - - - - - - - - - 12 即0zzG。.)()(21)()()(141)()()(1412/3202020002020202020200zyyxxzzzyyxxzzzyyxxznGzz代入（ 12-1-13）即得半空间的第一边值问题的解的积分公式dxdyzyyxxzyxfzyxu2/32020200000)()(),(21),(（12-2-21）作代换),(),(000zyxzyx002/32202000)()(1),(2),(dydxzyyxxyxfzzyxu这叫作半空间的 泊松积分 。例 3 在圆 a 内求解拉普拉斯方程的第一边值问题).(),(02fuaua答案.)()cos(212),(200020222dfaaau（12-2-22）例 4 在半平面 y0 内求解拉普拉斯方程的第一边值问题).(),0(002xfuyuy答案.)()(1),(00220dxxfyxxyyxu（12-2-23）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 16 页 - - - - - - - - - 13 三、含时间的格林函数 12.1 12.2 讨论的是稳定场问题的格林函数方法。至于波动与输运这类含时间的问题，同样可以运用格林函数方法求解。本节以波动问题为例介绍含时间的格林函数，并导出波动方程定解问题解的积分表式；对于输运问题，亦给出相应的结果。一般强迫振动的定解问题是),(2trfuautt（12-3-1）),(tMunu（12-3-2）).(),(00ruruttt（12-3-3） 5.3 中曾指出，持续作用的力f（r，t）可年作是前后相继的脉冲力f（r， ） （t ）d 的叠加。现在我们再进一步将一个个连续分布于空间的脉冲力看作是鳞次栉比排列在许许多多点上的力的叠加。总之，把持续作用的连续分布力f（r，t）看作是许许多多脉冲点力的叠加.)()(),(),(000Ttdrdtrrrftrf（12-3-4）把单位脉冲点力所引起的振动记作G（r，t；r0，t0） ，称之为波动问题的格林函数。求得了 G，就可用叠加的方法求出任意力f（r，t）所引起的振动。 G 所满足的定解问题是),()(002ttrrGaGtt（12-3-5）,0Gnu（12-3-6）. 0,000tttGG（12-3-7）我们可以用类似于求解泊松方程的方法求得定解问题（12-3-1）（12-3-3）的解的积分表式。需注意的是含时间的格林函数的对称性不同于泊松方程格林函数的对称性，).,;,(),;,(0000trtrGtrtrG（12-3-8）现在证明对称关系（ 12-3-8） 。在定解问题（ 12-3-5）（12-3-7）中将变量 t，r0，t0分别换为 t，r1，t1，而成为名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 16 页 - - - - - - - - - 14 )()(),;,(),;,(1111211ttrrtrtrGatrtrGtt（12-3-9）0),;,(),;,(1111trtrGntrtrG（12-3-10）.0),;,(,0),;,(011011ttttrtrGtrtrG（12-3-11）以 G（r，t；r1，t1）乘方程（12-3-5） 。同时以 G（r，t；r0，t0）乘方程（12-3-9） ，相减，再对 r 在区域 T 上积分，同时对 t 在区间, t（其中 tt0和 t1）上积分，得).,;,(),;,(),;,(),;,(),;,(),;,(),;,(),;,(),;,(),;,(0011110000121100200111100trtrGtrtrGdVdttrtrGtrtrGatrtrGtrtrGatrtrGtrtrGtrtrGtrtrGttTttt（12-3-12）利用第二格林公式（ 12-1-3） ，上式左端成为.),;,(),;,(),;,(),;,(),;,(),;,(),;,(),;,(00111100200111100dSdttrtrGntrtrGtrtrGntrtrGadVtrtrGtrtrGtrtrGtrtrGttttTtt由定解条件（ 12-3-6）（12-3-7）和（12-3-10）（12-3-11）可以看出，上式为零，从而（ 12-3-12）右端也为零。于是有对称关系（12-3-8） 。现在推导定解问题 （12-3-1） （12-3-3）解的积分表式。 考虑到关系式（12-3-8）中时间变数 t 与 t0不能像空间变数那样简单地对调，我们先将定解问题（12-3-1）（12-3-3）中的 r，t 换为 r0，t0，),(),(),(0000020000trftruatrutt（12-3-13）),(),(),(0000000tMtruntru（12-3-14）).(),(),(),(00000000000rtrurtruttt（12-3-15）将 G 的定解问题中的r 与 r0互换，同时将 t 和 t0分别换为 t0 和t，并利用对称关名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 16 页 - - - - - - - - - 15 系（12-3-8） ，得),()(),;,(),;,(0000020000ttrrtrtrGatrtrGtt（12-3-16）, 0),;,(),;,(00000trtrGntrtrG（12-3-17）.0),;,(,0),;,(00000000ttttrtrGtrtrG（12-3-18）以 G（r，t；r0，t0）乘方程（ 12-3-13） ，以 u（r0，t0）乘方程（ 12-3-16） ，相减，再对 r0在区域 T 上积分，同时对 t0在0，t 上积分，并利用第二格林公式及初始条件（12-3-15）及（ 12-3-18） ，可得.)()(),(),;,()()(0000000000000000020000000TtTtTtTtttttdtdVttrrudtdVtrftrtrGdtdVGuuGadtdVuGGu（12-3-19）其中0，积分后取0，引入是为了使含（t t0）的积分值确定（积分区间包含 t0t 在内） ，于是可得.)()(),(),;,(),(00000200000000000000TtTtttttTtdtdVGuuGadtdVuGGudtdVtrftrtrGtru(12-3-20) 右边第二个积分中0/)(000000dtuGGuduGGutttttt因此，可完成对t0的积分，计及 tt0时 G0，00tG，这样得到.),(),;,(),(000000020000000000TttttTtdVuGGudtdSnGunuGadtdVtrftrtrGtru（12-3-21）对于不同类型的边界条件，可令G 满足相应的齐次边界条件，从而得到适用于不同边界条件的以 G 表示的解的积分表式。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 16 页 - - - - - - - - - 16 对于输运问题，),(2trfuaut（12-3-22）),(tMunu（12-3-23）).(0rut（12-3-24）类似上面的讨论，同样可得到其解的积分表式.),(),;,(),(0000000200000000TttTtdVuGdtdSnGunuGadtdVtrftrtrGtru（12-3-25）作业（ P387） ：1，2名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 16 页 - - - - - - - - -