第四章格林函数法优秀PPT.ppt

第四章格林函数法第一页，本课件共有39页格林函数又称为格林函数又称为点源函数点源函数或或影响函数影响函数。顾名思。顾名思义，它表示一个点源在一定的边界条件和义，它表示一个点源在一定的边界条件和(或或)初值条初值条件下所产生的场或影响。由于任意分布的源所产生的件下所产生的场或影响。由于任意分布的源所产生的场均可看成许许多多点源产生的场的叠加，因此格林场均可看成许许多多点源产生的场的叠加，因此格林函数一旦求出，就可算出任意源的场。格林函数法以函数一旦求出，就可算出任意源的场。格林函数法以统一的方式处理各类数学物理方程，既可以研究常微统一的方式处理各类数学物理方程，既可以研究常微分方程，又可以研究偏微分方程；既可以研究齐次方分方程，又可以研究偏微分方程；既可以研究齐次方程又可以研究非齐次方程；既可以研究有界问题，又程又可以研究非齐次方程；既可以研究有界问题，又可以研究无界问题。它的内容十分丰富，应用极其广可以研究无界问题。它的内容十分丰富，应用极其广泛。这一章，我们主要介绍用格林函数求解拉普拉斯泛。这一章，我们主要介绍用格林函数求解拉普拉斯方程的边值问题。方程的边值问题。第二页，本课件共有39页4.1 格林公式及其应用4.1.1 基本解基本解对拉普拉斯方程对拉普拉斯方程,其球坐标形式为：其球坐标形式为：(4.1.1)求方程求方程(4.1.1)(4.1.1)的球对称解的球对称解(即与即与和和无关的解无关的解),),则有：则有：其通解为：其通解为：为任意常数为任意常数)。若取若取,则得到特解则得到特解,称此解为三维称此解为三维Laplace 方程的方程的基本解基本解,它在研究三维拉普拉斯方程中起着重要的作用它在研究三维拉普拉斯方程中起着重要的作用.第三页，本课件共有39页对二维拉普拉斯方程对二维拉普拉斯方程,其极坐标形式为其极坐标形式为:(4.1.2)求方程求方程(4.1.2)(4.1.2)的径向对称解的径向对称解(即与即与无关的解无关的解),),则有：则有：其通解为：其通解为：为任意常数为任意常数)。若取若取,则得到特解则得到特解,称此解为二维称此解为二维Laplace方程的方程的基本解基本解.第四页，本课件共有39页4.1.2 格林公式格林公式由高斯公式由高斯公式,则得到则得到格林第一公式格林第一公式：令令 将以上两公式相减，得到将以上两公式相减，得到格林第二公式格林第二公式：调和函数调和函数：具有二阶偏导数并且满足拉普拉斯方程的连续函数。：具有二阶偏导数并且满足拉普拉斯方程的连续函数。第五页，本课件共有39页4.1.3 调和函数的积分表达式由由Green公式可导出调和函数的积分表示。由于函数：公式可导出调和函数的积分表示。由于函数：除在除在 点外处处满足三维点外处处满足三维Laplace方程方程，于是有，于是有 定理：若函数定理：若函数 在在 上有一阶连续偏导数，且在上有一阶连续偏导数，且在 内调和，则内调和，则 调和函数在区域内任一点的值可以通过积分表达式用这个调和函数在区域内任一点的值可以通过积分表达式用这个函数在区域边界上的值和边界上的法向导数来表示。函数在区域边界上的值和边界上的法向导数来表示。第六页，本课件共有39页 若函数若函数 在在 上有一阶连续偏导数，且在上有一阶连续偏导数，且在 内满足内满足Poisson方程方程 ，则同样有，则同样有 4.1.4 调和函数的性质调和函数的性质 性质性质1.设设 是区域是区域 内的调和函数，它在内的调和函数，它在上有一阶连续偏导数，则上有一阶连续偏导数，则其中其中 的外法线方向。的外法线方向。是是证明证明 只要在只要在Green公式中取公式中取 即证。即证。注：此性质表明调和函数的法向导数沿区域边界的积分为零。注：此性质表明调和函数的法向导数沿区域边界的积分为零。对稳定的温度场，流入和流出物体界面的热量相等，否则就对稳定的温度场，流入和流出物体界面的热量相等，否则就不能保持热的动态平衡，而使温度场不稳定。不能保持热的动态平衡，而使温度场不稳定。第七页，本课件共有39页 思考：思考：Laplace方程方程Neumann问题有解的必要条件是什么？问题有解的必要条件是什么？性质性质2(平均值定理平均值定理)设函数设函数在区域在区域 内调和，内调和，是是 内任意一点，若内任意一点，若是以是以 为中心，为中心，a为半径为半径的球面，此球完全落在区域的球面，此球完全落在区域 的内部，则有的内部，则有证明：证明：由调和函数的积分表示：由调和函数的积分表示：及由性质及由性质1，有，有 第八页，本课件共有39页上式称为调和函数的球面平均值公式。上式称为调和函数的球面平均值公式。又因为，在又因为，在 上有上有，所以，所以 性质性质3(极值原理极值原理)设函数设函数在区域在区域 内调和，内调和，它在它在上连续且不为常数，则它的最大值与最小值上连续且不为常数，则它的最大值与最小值只能在边界上达到。只能在边界上达到。推论推论1 设在设在 内有内有在在上连续且在边界上连续且在边界上有上有，则在，则在内有内有推论推论2 Dirichlet问题问题 的解是唯一的。的解是唯一的。第九页，本课件共有39页第十页，本课件共有39页第十一页，本课件共有39页4.2 格林函数格林函数由于调和函数有积分表示由于调和函数有积分表示:又因为又因为Dirichlet边值问题边值问题 的解唯一的解唯一,故希望故希望将此问题的解用积分表示出来。但由于在积分表达示中，将此问题的解用积分表示出来。但由于在积分表达示中，u在边界上的值虽然已知在边界上的值虽然已知,而而 在边界上的值却不知道在边界上的值却不知道.那么那么,能否作为边界条件加上能否作为边界条件加上 的值呢？的值呢？因为因为,此时的解已经是唯一的了此时的解已经是唯一的了.那么只有想办法去掉那么只有想办法去掉 为此，引入格林函数的概念。为此，引入格林函数的概念。显然这是行不通的，显然这是行不通的，(4.2.1)第十二页，本课件共有39页格林函数的物理背景格林函数的物理背景原点处点电荷电量原点处点电荷电量 ，点电荷密度点电荷密度处点电位处点电位即即 处点电荷电量处点电荷电量点电荷密度点电荷密度处点电位处点电位第十三页，本课件共有39页4.2.1 格林函数的定义格林函数的定义设在设在 内有内有在在上有一阶连续上有一阶连续偏导数，则由格林第二公式有偏导数，则由格林第二公式有(4.2.2)将（将（4.2.1）和）和(4.2.2)两式加起来：两式加起来：(4.2.3)选择调和函数选择调和函数v满足满足,于是有：于是有：(4.2.4)第十四页，本课件共有39页记记(4.2.5)则有则有(4.2.6)称称 为为Laplace方程的格林函数。若方程的格林函数。若上有一阶连续偏导数，则当上有一阶连续偏导数，则当Dirichlet问题问题且在且在 上具有一阶连续偏导数的解存在时，解可以表示为上具有一阶连续偏导数的解存在时，解可以表示为在在(4.2.7)存在存在 第十五页，本课件共有39页对对Poisson方程的方程的Dirichlet问题问题 上存在具有一阶连续偏导数的解，则解可以如果在如果在表示为表示为由此可见由此可见,求解求解Dirichlet问题问题,关键是求关键是求Green函数函数(4.2.5),其中其中v满足一个特殊的满足一个特殊的Dirichlet问题：问题：(4.2.8)称由函数称由函数v确定的格林函数为确定的格林函数为第一边值问题的格林函数第一边值问题的格林函数。第十六页，本课件共有39页4.2.2 格林函数的性质格林函数的性质1.格林函数格林函数在除去点在除去点 外处处满足外处处满足 Laplace方程，当方程，当 时，时，其阶数与其阶数与 相同。相同。2.在边界上，格林函数恒等于零：在边界上，格林函数恒等于零：3.在区域在区域 内成立不等式：内成立不等式：（用极值原理证明）（用极值原理证明）4.（由格林第二公式证明）（由格林第二公式证明）5.第十七页，本课件共有39页4.3 格林函数的应用格林函数的应用 用镜象法求特殊区域上的函数。用镜象法求特殊区域上的函数。4.3.1 上半空间内的Green函数及Dirichlet问题 求解上半空间求解上半空间 内的内的Dirichlet问题问题 先求上半空间先求上半空间 内的内的Green函数函数（4.3.1），即求解问题，即求解问题 第十八页，本课件共有39页 在区域外找出区域内一点关于边界的象点，在这两个点放置适当在区域外找出区域内一点关于边界的象点，在这两个点放置适当的电荷，这两个电荷产生的电位在曲面边界上相互抵消。这两个电的电荷，这两个电荷产生的电位在曲面边界上相互抵消。这两个电荷在区域中形成的电位就是所要求的格林函数。荷在区域中形成的电位就是所要求的格林函数。第十九页，本课件共有39页于是，半空间上的格林函数为于是，半空间上的格林函数为(4.3.2)从而，问题从而，问题(4.3.1)的解可表示为的解可表示为:由于平面由于平面z=0上的外法线方向即上的外法线方向即oz轴的负向轴的负向,所以所以 即即 所以，问题所以，问题(4.3.1)的解为的解为:第二十页，本课件共有39页例例2 求解下列定解问题求解下列定解问题解：解：第二十一页，本课件共有39页4.3.2 球域上的球域上的Green函数及函数及Dirichlet问题问题 其中，其中，（4.3.3），即求解问题，即求解问题 求解球域上的求解球域上的Dirichlet问题问题 是以坐标原点是以坐标原点O为球心，为球心，R为半径的球域。为半径的球域。先求球域上的先求球域上的Green函数函数第二十二页，本课件共有39页第二十三页，本课件共有39页球内的格林函数球内的格林函数 M0点处点电荷电量点处点电荷电量 ，M1点处点电荷电量点处点电荷电量 第二十四页，本课件共有39页第二十五页，本课件共有39页从而，问题从而，问题(4.3.3)的解可表示为：的解可表示为：因因其中其中是是与与的夹角，于是：的夹角，于是：（4.3.4）此公式称为球域上的泊松积分公式。如果用球坐标表示，则有此公式称为球域上的泊松积分公式。如果用球坐标表示，则有（4.3.5）其中其中 是点是点 的球坐标，的球坐标，是是 上动点的坐标上动点的坐标，第二十六页，本课件共有39页是是与与的夹角。由于的夹角。由于 所以所以(4.3.6)第二十七页，本课件共有39页例例1.设有一半径为设有一半径为R的均匀球，上半球面的温度保持为的均匀球，上半球面的温度保持为。求球内温度的稳定分布。求球内温度的稳定分布。下半球面的温度保持为下半球面的温度保持为 解：考虑定解问题解：考虑定解问题 由球域上的泊松积分公式由球域上的泊松积分公式(4.3.5)，得，得 第二十八页，本课件共有39页由于此积分的计算很困难，下面我们只考虑一些特殊位置的由于此积分的计算很困难，下面我们只考虑一些特殊位置的温度分布。比如，求温度在球的铅垂直径温度分布。比如，求温度在球的铅垂直径（直径的上（直径的上半部）和半部）和（直径的下半部分）上的分布。（直径的下半部分）上的分布。当当 时，时，(见见(4.3.6)式式)，故有：，故有：第二十九页，本课件共有39页当当 时，时，,故有故有 在以上两个公式中，当在以上两个公式中，当 时，球的温度为时，球的温度为 .第三十页，本课件共有39页4.3.3 4.3.3 四分之一空间的格林函数四分之一空间的格林函数 第三十一页，本课件共有39页4.4 试探法及试探法及Poisson方程的求解方程的求解 4.4.1 试探法试探法 对某些定解问题对某些定解问题,根据问题的物理意义和几何特征，可假设根据问题的物理意义和几何特征，可假设解具有某种特殊形式，将这种形式的解代入方程进行试探直至求解具有某种特殊形式，将这种形式的解代入方程进行试探直至求出特解。这种方法称为试探法。出特解。这种方法称为试探法。第三十二页，本课件共有39页例例1.设有一半径为设有一半径为R的无限均匀圆柱体的无限均匀圆柱体,已知圆柱内无热源已知圆柱内无热源,圆柱圆柱面上的温度分布为面上的温度分布为,试求圆柱内温度的稳定分布试求圆柱内温度的稳定分布.解：因柱面上温度与解：因柱面上温度与z无关无关,则域内温度也应与则域内温度也应与z无关无关,故原问题故原问题可简化为求解圆域上可简化为求解圆域上Laplace方程的第一边值问题方程的第一边值问题,采用极坐标采用极坐标,我们考虑问题我们考虑问题:由由(4.4.2),设设(4.4.1)得得,代入代入,再由再由(4.4.2)得得 由 的任意性得:第三十三页，本课件共有39页例例2 求圆柱域求圆柱域 内的电位内的电位u,使在柱面上有给定的电场强度使在柱面上有给定的电场强度的法向分量的法向分量,即即 解解:由边界条件知，问题可化为平面问题：由边界条件知，问题可化为平面问题：由边界条件由边界条件(4.4.4),设设,显然显然 满足方程满足方程(4.4.3)及条件及条件(4.4.4)，于是问题的解为：，于是问题的解为：第三十四页，本课件共有39页例例3 3 求由两同心球面导体求由两同心球面导体 和和 构成的电容器内构成的电容器内的电位，使内球面的电位，使内球面 保持常电位保持常电位 外球面接地。外球面接地。解解:采用球坐标，考虑定解问题采用球坐标，考虑定解问题 由边界条件知，球内电位的分布仅与由边界条件知，球内电位的分布仅与r有关，即电位有关，即电位函数是球对称的，而电位与函数是球对称的，而电位与r成反比，故可设成反比，故可设 第三十五页，本课件共有39页显然显然 满足满足(4.4.5),这是因为这是因为,是三维是三维Laplace方程方程的基本解。由的基本解。由(4.4.6)于是于是(4.4.5)(4.4.6)的解为：的解为：第三十六页，本课件共有39页如果知道如果知道Poisson方程的一个特解，则通过函数代换，方程的一个特解，则通过函数代换，4.4.2 Poisson方程的求解方程的求解 就可将就可将Poisson方程边值问题化成方程边值问题化成Laplace方程的边值问题。方程的边值问题。例例1 求求 的特解。的特解。解解:设其特解为设其特解为,则则 于是，其解有无穷多个，如于是，其解有无穷多个，如 等等。等等。第三十七页，本课件共有39页例例2 求下列问题的解求下列问题的解 解解:显然方程有一个特解显然方程有一个特解,故令故令,则则 由极值原理，上述问题的解为由极值原理，上述问题的解为,故原问题的解为：第三十八页，本课件共有39页4.4.3 Dirichlet外问题与外问题与Neumann外问题简介外问题简介 Dirichlet内问题内问题Dirichlet外问题外问题Neumann内问题内问题Neumann外问题外问题由于外问题在无穷区域上提出，需附加条件：由于外问题在无穷区域上提出，需附加条件：其中，其中，。从数学角度来讲，此条件。从数学角度来讲，此条件可以保证外问题的解是唯一的。可以保证外问题的解是唯一的。第三十九页，本课件共有39页