数理方程第四章 格林函数法

1、关于数理方程第四章 格林函数法第一张，PPT共三十九页，创作于2022年6月25.08.20221格林函数又称为点源函数或影响函数。顾名思义，它表示一个点源在一定的边界条件和(或)初值条件下所产生的场或影响。由于任意分布的源所产生的场均可看成许许多多点源产生的场的叠加，因此格林函数一旦求出，就可算出任意源的场。格林函数法以统一的方式处理各类数学物理方程，既可以研究常微分方程，又可以研究偏微分方程；既可以研究齐次方程又可以研究非齐次方程；既可以研究有界问题，又可以研究无界问题。它的内容十分丰富，应用极其广泛。这一章，我们主要介绍用格林函数求解拉普拉斯方程的边值问题。第二张，PPT共三十九页，创作

2、于2022年6月25.08.202224.1 格林公式及其应用4.1.1 基本解对拉普拉斯方程, 其球坐标形式为：(4.1.1)求方程(4.1.1)的球对称解(即与和无关的解) ,则有： 其通解为： 为任意常数)。 若取 ,则得到特解 ,称此解为三维Laplace 方程的基本解,它在研究三维拉普拉斯方程中起着重要的作用.第三张，PPT共三十九页，创作于2022年6月25.08.20223对二维拉普拉斯方程 ,其极坐标形式为:(4.1.2)求方程(4.1.2)的径向对称解(即与无关的解) ,则有： 其通解为： 为任意常数)。 若取 , 则得到特解 , 称此解为二维Laplace方程的基本解.第四

3、张，PPT共三十九页，创作于2022年6月25.08.202244.1.2 格林公式由高斯公式 ,则得到格林第一公式： 令 将以上两公式相减，得到格林第二公式： 调和函数：具有二阶偏导数并且满足拉普拉斯方程的连续函数。第五张，PPT共三十九页，创作于2022年6月25.08.202254.1.3 调和函数的积分表达式由Green公式可导出调和函数的积分表示。由于函数： 除在 点外处处满足三维Laplace方程 ，于是有 定理：若函数 在 上有一阶连续偏导数，且在 内调和，则 调和函数在区域内任一点的值可以通过积分表达式用这个函数在区域边界上的值和边界上的法向导数来表示。第六张，PPT共三十九页

4、，创作于2022年6月25.08.20226 若函数 在 上有一阶连续偏导数，且在 内满足Poisson方程 ，则同样有 4.1.4 调和函数的性质 性质1. 设 是区域 内的调和函数，它在上有一阶连续偏导数，则其中 的外法线方向。 是证明 只要在Green公式中取 即证。 注：此性质表明调和函数的法向导数沿区域边界的积分为零。对稳定的温度场，流入和流出物体界面的热量相等，否则就不能保持热的动态平衡，而使温度场不稳定。 第七张，PPT共三十九页，创作于2022年6月25.08.20227 思考：Laplace方程Neumann问题有解的必要条件是什么？性质2 (平均值定理) 设函数在区域 内调

5、和，是 内任意一点，若是以 为中心，a为半径的球面，此球完全落在区域 的内部，则有证明： 由调和函数的积分表示： 及由性质1，有 第八张，PPT共三十九页，创作于2022年6月25.08.20228上式称为调和函数的球面平均值公式。又因为，在 上有，所以 性质3 (极值原理) 设函数在区域 内调和，它在上连续且不为常数，则它的最大值与最小值只能在边界上达到。 推论1 设在 内有在上连续且在边界上有，则在内有推论2 Dirichlet问题 的解是唯一的。第九张，PPT共三十九页，创作于2022年6月25.08.20229第十张，PPT共三十九页，创作于2022年6月25.08.202210第十一

6、张，PPT共三十九页，创作于2022年6月25.08.2022114.2 格林函数由于调和函数有积分表示: 又因为Dirichlet边值问题 的解唯一,故希望将此问题的解用积分表示出来。但由于在积分表达示中，u在边界上的值虽然已知,而 在边界上的值却不知道.那么,能否作为边界条件加上 的值呢？ 因为,此时的解已经是唯一的了.那么只有想办法去掉 为此，引入格林函数的概念。 显然这是行不通的，(4.2.1)第十二张，PPT共三十九页，创作于2022年6月25.08.202212格林函数的物理背景原点处点电荷电量 ，点电荷密度处点电位即 处点电荷电量点电荷密度处点电位第十三张，PPT共三十九页，创作

7、于2022年6月25.08.2022134.2.1 格林函数的定义设在 内有在上有一阶连续偏导数，则由格林第二公式有 (4.2.2)将（4.2.1）和 (4.2.2)两式加起来： (4.2.3)选择调和函数v满足 ,于是有： (4.2.4)第十四张，PPT共三十九页，创作于2022年6月25.08.202214记 (4.2.5)则有 (4.2.6)称 为Laplace方程的格林函数。若上有一阶连续偏导数，则当Dirichlet问题且在 上具有一阶连续偏导数的解存在时，解可以表示为在(4.2.7)存在 第十五张，PPT共三十九页，创作于2022年6月25.08.202215对Poisson方程的

8、Dirichlet问题 上存在具有一阶连续偏导数的解，则解可以如果在表示为由此可见,求解Dirichlet问题,关键是求Green函数(4.2.5),其中v满足一个特殊的Dirichlet问题： (4.2.8)称由函数v确定的格林函数为第一边值问题的格林函数。第十六张，PPT共三十九页，创作于2022年6月25.08.2022164.2.2 格林函数的性质1. 格林函数在除去点 外处处满足 Laplace方程，当 时，其阶数与 相同。 2. 在边界上，格林函数恒等于零：3. 在区域 内成立不等式： （用极值原理证明） 4. （由格林第二公式证明） 5. 第十七张，PPT共三十九页，创作于202

9、2年6月25.08.2022174.3 格林函数的应用 用镜象法求特殊区域上的函数。 4.3.1 上半空间内的Green函数及Dirichlet问题 求解上半空间 内的Dirichlet问题 先求上半空间 内的Green函数 （4.3.1） ，即求解问题 第十八张，PPT共三十九页，创作于2022年6月25.08.202218 在区域外找出区域内一点关于边界的象点，在这两个点放置适当的电荷，这两个电荷产生的电位在曲面边界上相互抵消。这两个电荷在区域中形成的电位就是所要求的格林函数。第十九张，PPT共三十九页，创作于2022年6月25.08.202219于是，半空间上的格林函数为(4.3.2)

10、从而，问题(4.3.1)的解可表示为: 由于平面z=0上的外法线方向即oz轴的负向,所以 即 所以，问题(4.3.1)的解为: 第二十张，PPT共三十九页，创作于2022年6月25.08.202220例2 求解下列定解问题解：第二十一张，PPT共三十九页，创作于2022年6月25.08.2022214.3.2 球域上的Green函数及Dirichlet问题 其中， （4.3.3） ，即求解问题 求解球域上的Dirichlet问题 是以坐标原点O为球心，R为半径的球域。 先求球域上的Green函数第二十二张，PPT共三十九页，创作于2022年6月25.08.202222第二十三张，PPT共三十九

11、页，创作于2022年6月25.08.202223球内的格林函数 M0点处点电荷电量 ，M1点处点电荷电量 第二十四张，PPT共三十九页，创作于2022年6月25.08.202224第二十五张，PPT共三十九页，创作于2022年6月25.08.202225从而，问题(4.3.3)的解可表示为： 因其中是与的夹角，于是： （4.3.4） 此公式称为球域上的泊松积分公式。如果用球坐标表示，则有 （4.3.5） 其中 是点 的球坐标， 是 上动点的坐标， 第二十六张，PPT共三十九页，创作于2022年6月25.08.202226是与的夹角。由于 所以 (4.3.6) 第二十七张，PPT共三十九页，创作

12、于2022年6月25.08.202227例1. 设有一半径为R的均匀球，上半球面的温度保持为 。求球内温度的稳定分布。 下半球面的温度保持为 解：考虑定解问题 由球域上的泊松积分公式(4.3.5)，得 第二十八张，PPT共三十九页，创作于2022年6月25.08.202228由于此积分的计算很困难，下面我们只考虑一些特殊位置的温度分布。比如，求温度在球的铅垂直径 （直径的上半部）和（直径的下半部分）上的分布。 当 时， (见(4.3.6)式)，故有： 第二十九张，PPT共三十九页，创作于2022年6月25.08.202229当 时， ,故有 在以上两个公式中，当 时，球的温度为 . 第三十张，

13、PPT共三十九页，创作于2022年6月25.08.2022304.3.3 四分之一空间的格林函数 第三十一张，PPT共三十九页，创作于2022年6月25.08.2022314.4 试探法及Poisson方程的求解 4.4.1 试探法 对某些定解问题,根据问题的物理意义和几何特征，可假设解具有某种特殊形式，将这种形式的解代入方程进行试探直至求出特解。这种方法称为试探法。 第三十二张，PPT共三十九页，创作于2022年6月25.08.202232例1. 设有一半径为R的无限均匀圆柱体,已知圆柱内无热源,圆柱面上的温度分布为 ,试求圆柱内温度的稳定分布. 解：因柱面上温度与z无关,则域内温度也应与z

14、无关,故原问题可简化为求解圆域上Laplace方程的第一边值问题,采用极坐标,我们考虑问题: 由(4.4.2),设 (4.4.1)得 ,代入, 再由(4.4.2)得 由 的任意性得: 第三十三张，PPT共三十九页，创作于2022年6月25.08.202233例2 求圆柱域 内的电位u,使在柱面上有给定的电场强度的法向分量,即 解: 由边界条件知，问题可化为平面问题：由边界条件(4.4.4),设 , 显然 满足方程(4.4.3)及条件(4.4.4)，于是问题的解为： 第三十四张，PPT共三十九页，创作于2022年6月25.08.202234例3 求由两同心球面导体 和 构成的电容器内的电位，使内球面 保持常电位 外球面接地。 解: 采用球坐标，考虑定解问题 由边界条件知，球内电位的分布仅与r有关，即电位函数是球对称的，而电位与r成反比，故可设 第三十五张，PPT共三十九页，创作于2022年6月25.08.202235显然 满足(4.4.5), 这是因为, 是三维Laplace方程的基本解。由(4.4.6) 于是(4.4.5) (4.4.6)的解为：第三十六张，PPT共三十九页，创作于2022年6月25.08.202236如果知道Poisson方程的一个特解，则通过函数代换，4.4.2 Poisson方程的求解 就可将Poisson