数学物理方法12格林函数.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第十二章 格林函数法,格林（Green）函数,，又称为点源影响函数，是数学物理中的一个重要概念格林函数代表一个点源在一定的边界条件下和初始条件下所产生的场知道了点源的场，就可以用叠加的方法计算出任意源所产生的场,格林函数法是解数学物理方程的常用方法之一,一、格林公式,上具有连续一阶导数,在区域,及其边界,和,中具有连续二阶导数，应用矢量分析的高斯定理,(12.1.1),将对曲面,的积分化为体积分,12.1 泊松方程的格林函数法,以上用到公式,称上式为,第一格林公式,同理有,上述两式相减得到,表示沿边界,的外法向偏导数,称为,第二格林公式,进一步改写为,二、泊松方程的格林函数法,1、讨论具有一定边界条件的泊松方程的定解问题,泊松方程,边值条件,是区域边界,上给定的函数.,是第一、第二、第三类边界条件的统一描述,典型的泊松方程（三维稳定分布）边值问题,表示边界面,上沿界面外法线方向的偏导数,泊松方程,第一类边界条件：第一边值问题(狄里希利问题),第二类边界条件：第二边值问题(诺依曼问题),第三类边界条件：第三边值问题,2、格林函数的引入及其物理意义,引入：为了求解泊松方程的定解问题，我们必须定义一个与此定解问题相应的格林函数,它满足如下定解问题，边值条件可以是第一、二、三类条件：,代表三维空间变量的函数，在直角坐标系中其形式为,格林函数的物理意义,：,在区域,T,内部,处放置一个点源，而在该区域,T,的界面上为零的条件下,那么该点点源在区域,T,内r处产生的场，由此可以进一步理解通常人们为什么称格林函数为点源函数,格林函数互易定理,：,因为格林函数,代表,处的点源在,处所产生的影响（或所产生的场）,所以它只能是距离,的函数,故它应该遵守如下的互易定理：,根据格林第二公式,令,得到,根据,函数性质有：,故有,称为,泊松方程的基本积分公式,格林函数满足互易定理,并利用格林函数的对称性则得到,考虑格林函数所满足的边界条件讨论如下：,第一类边值问题：,相应的格林函数,是下列问题的解：,考虑到格林函数的,齐次边界条件,，,第一类边值问题的解,2.第二类边值问题,相应的格林函数,是下列问题的解：,由公式可得,第二类边值问题解,3.第三类边值问题,相应的格林函数,是下列问题的解：,泊松方程的边值条件，两边同乘以格林函,格林函数的边值条件的两边同乘以函数,得,相减得到,代入（14.2.9）得到,第三类边值问题的解,这就是第三边值问题解的积分表示式,右边第一个积分表示区域,中分布的源,在,点产生的场的总和.,第二个积分则代表边界上的状况对,点场的影响的总和两项积分中的格林函数相同这说明泊松方程的格林函数是点源在一定的边界条件下所产生的场,对于,拉普拉斯方程,第一边值问题的解为,第三边值问题的解,为,12.2,无界空间的格林函数 基本解,无界区域这种情形公式中的面积分应为零，故有,选取,和,分别满足下列方程,三维球对称,对于三维球对称情形，我们选取,对上式两边在球内积分,（14.3.4）,（14.3.5）,利用高斯定理（14.1.1）得到,(14.3.6),故有,使上式恒成立，有,因此,，,故得到,对于三维无界球对称情形的格林函数可以选取为,代入（14.3.1）得到三维无界区域问题的解为,上式正是我们所熟知的静电场的电位表达式,二维轴对称情形,用单位长的圆柱体来代替球积分在单位长的圆柱体内进行，即,因为,由于,只是垂直于轴，且向外的分量，所以上式在,圆柱体上、下底的面积分为零，只剩下沿侧面的积分，即,选取的圆柱的高度为单位长，则很容易得到下面的结果,令积分常数为0，得到,因此二维轴对称情形的格林函数为,将（14.3.9）代入式（14.3.1）得到二维无界区域的解为,用电像法确定格林函数,用格林函数法求解的主要困难还在于如何确定格林函数本身,一个具体的定解问题，需要寻找一个合适的格林函数,为了求解的方便，对一些具体问题我们给出构建格林函数的方法,电像法,考虑一个具体的物理模型：设在一接地导体球内的,放置一个单位正电荷，求在体内的电势分布，并满足边界条件为零,点,对于,第一类边值问题,，其格林函数可定义为下列定解问题的解,为了满足边界条件：电势为零，所以还得在,边界外像点（或对称点）放置,一个合适的负电荷，这样才能使这两个电荷在界面上产生的电势之和为零,这方法是基于,静电学的镜像原理来构建格林函数,，所以我们称这种构建方法为电像法（也称为镜像法）,上半平面区域第一边值问题的格林函数构建,拉普拉斯方程的第一边值问题求解,物理模型：,若在,处放置一,正单位点电荷,则虚设的负,单位点电荷,应该在,于是得到这两点,电荷,在,xoy,的上半平面的电位分布也就是本问题的格林函数，即为,（14.4.2）,据上述物理模型可求解下列定解问题,例14.6.1 定解问题：,【解】根据第一边值问题，构建的格林函数满足,处放置于一个正和一个负的点电荷（或点源）,构,建格林函数为,边界外法线方向为负,轴，故有,代入到拉普拉斯第一边值问题解的公式（14.2.13），拉普拉斯方程的自由项,则由,得,(14.4.3),或代入,拉普拉斯方程的第一边值问题的解公式（14.2.22）,得到,(14.4.4),公式（14.4.3）或（14.4.4）称为,上半平面的拉普拉斯积分公式,2.泊松方程的第一边值问题求解,例14.6.2 定解问题：,根据,第一类边值问题的解公式,(14.2.14),得到,（14.4.5）,根据半平面区域第一类边值问题的格林函数(14.4.2)式，得到,(14.4.6),因为边界上的法线为负y轴，故,(14.4.7),将（14.4.6）和(14.4.7)代入（14.4.5）得到泊松方程在半平面区域第一边值问题的解,14.4.2 上半空间内求解拉普拉斯方程的第一边值问题,物理模型:,例14.4.3 在上半空间,内求解拉普拉斯方程的第一边值问题,【,解,】构建格林函数,满足,根据物理模型和无界区域的格林函数可以构建为,(14.4.8),即有,为了把,代入,拉普拉斯,第一边值问题的解的公式（14.2.22），,需要先计算,即为,代入（14.2.22）即得到,这公式叫作,半空间的拉普拉斯积分,（14.4.9）,14.4.3 圆形区域第一边值问题的格林函数构建,物理模型【2】：,在圆内任找一点,放置一个单位,根据图14.2，这两线电荷在圆内任一观察点,所产生的电势为,当观察点,位于圆周上,时，应该有,，即满足第一类齐次边值条件,即为,上式应对任何,值成立，所以上式对,的导数应为零，即,即得到,要求上式对任意的,值要成立，故提供了确定,的方程,联立解得,于是圆形区域,的第一类边值问题的格林函数为,（14.4.10）,即为,(14.4.11),.,其中,例14.4.4 求解如下泊松方程定解问题,根据第一类边值问题解的公式(14.2.14)，并取沿垂直于圆的方向取单位长积分，这样原来的体积分化为面积分，原来的面积分化为线积分故得到,根据构建的圆内第一边值问题的格林函数（14.4.11）,（14.4.12）,代入得到圆内第一边值问题的解为,(14.4.13),例14.4.5,在圆,内求解拉普拉斯方程的第一边值问题,【,解,】根据公式（14.4.14），故有,（14.4.14）,