数学物理方法第十章格林函数法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
T
S
G(r ,0) dS
G 2 S r r sin d d
故有
G 2 S r r sin d d T G(r ,0)dV 1
使上式恒成立，有
G(r ,0) 4πr 1 r
2
1 G(r ,0) c 4πr
r ， G 0
G (r , r0 ) u (r0 ) G (r , r0 ) f (r )dV (r ) ]dS T n
得
y0 u ( x0 , y0 ) π
( x)
( x x0 ) y
2 2 0
dx
或代入拉普拉斯方程的第一边值问题的解公式
G (r , r0 ) u (r ) (r0 ) ]dS0 n 0 得到
10.4 用电像法确定格林函数
用格林函数法求解的主要困难还在于如何确定格林函数本身 一个具体的定解问题，需要寻找一个合适的格林函数 为了求解的方便，对一些具体问题我们给出构建格林函数的方法 一、电像法定义 考虑一个具体的物理模型：设在一接地导体球内的 M 0 点 放置一个单位正电荷，求在体内的电势分布，并满足边界条件为零
边界外法线方向为负 y 轴，故有
y0 y0 y0 G G 1 1 1 | | y 0 = 2 2 2 2 2 n y 2π ( x x0 ) y0 π ( x x0 ) y0 π ( x x0 )2 y0
代入到拉普拉斯第一边值问题解的公式，拉普拉斯方程的自由 项 ,则由 f 0
因此
c0
,故得到
对于三维无界球对称情形的格林函数可以选取为
G(r , r0 )
代入
1 4π | r r0 |
u (r ) G(r , r0 ) f (r0 )dV0
T0
得到三维无界区域问题的解为
f (r0 ) 1 u (r ) dV0 T 0|r r | 4π 0
上式正是我们所熟知的静电场的电势表达式
二、二维轴对称情形
用单位长的圆柱体来代替球．积分在单位长的圆柱体内进行，即
G(r ,0)dV (r )dV
T T
因为
(r )dV 1
T
G(r ,0)dV G(r ,0)dV
T T
一、三维球对称
对于三维球对称情形，我们选取
r0 0
G(r , r0 ) (r - r0 )
两边在球内积分
G(r ,0)dV (r )dV
T T
(r )dV 1
T
利用高斯定理得到
T
G(r ,0)dV G(r ,0)dV
G ( x, y | x0 , y0 )
( x x0 ) 2 ( y y0 ) 2 1 ln[ ] 2 2 4π ( x x0 ) ( y y0 )
据上述物理模型可求解下列定解问题 例1 定解问题：
u xx u yy 0, ( y 0) u | y 0 ( x)
根据第一类边值问题的解公式得到
u ( x, y )
解： 根据第一边值问题，构建的格林函数满足
2G Gxx Gyy ( x x0 ) ( y y0 )
G | y 0 0
( x0 , y0 ),( x0 , y0 ) 处放置于一个正和一个负的点电荷（或点源）
构建格林函数为
( x x0 )2 ( y y0 )2 1 G ( x, y | x0 , y0 ) ln[ ] 2 2 4π ( x x0 ) ( y y0 )
对于第一类边值问题，其格林函数可定义为下列定解问题的解
G(r , r0 ) (r - r0 ) G(r , r0 ) | 0
为了满足边界条件：电势为零，所以还得在边界外像 点（或对称点）放置一个合适的负电荷，这样才能使这两 个电荷在界面上产生的电势之和为零
这方法是基于静电学的镜像原理来构建格林函数，所 以我们称这种构建方法为电像法（也称为镜像法）．
令积分常数为0，得到
G(r ,0)
1 1 G(r ,0) ln c 2π r
1 1 ln 2π r
因此二维轴对称情形的格林函数为
G(r , r0 ) 1 1 ln 2π | r r0 |
得到二维无界区域的解为
u (r ) 1 1 f ( r )ln dS0 0 S 2π 0 | r r0 |
二、 上半平面区域第一边值问题的格林函数构建
拉普拉斯方程的第一边值问题求解
物理模型：若在
M 0 ( x0 , y0 ) 处放置一正单位点电荷
则虚设的负单位点电荷应该在 M1 ( x0 , y0 ) 于是得到这两点电荷在 xoy 的上半平面的电位分 布．也就是本问题的格林函数，即为
G (r , r0 ) 1 1 1 1 ln ln 2π | r r0 | 2 π | r r1 | 1 1 1 1 ln ln 2π ( x x0 ) 2 ( y y0 ) 2 2π ( x x0 ) 2 ( y y0 ) 2
g ( x0 ) y u ( x, y ) dx0 2 2 π ( x x0 ) y
称为上半平面的拉普拉斯积分公式．
三、 泊松方程的第一边值问题求解
例2
定解问题：
u xx u yy f ( x, y ) u ( x,0) ( x)
( <x <+, y 0) ( <x<+, y 0)
S
G(r ,0) dS
由于
G
源自文库
G er , G r
只是垂直于轴，且向外的分量，所以上式在
圆柱体上、下底的面积分为零，只剩下沿侧面的积分，即
G r rddz T (r )dV 1
选取的圆柱的高度为单位长，则很容易得到下面的结果
G 1 r 2πr
S
G(r ,0) dS
G 2 S r r sin d d
故有
G 2 S r r sin d d T G(r ,0)dV 1
使上式恒成立，有
G(r ,0) 4πr 1 r
2
1 G(r ,0) c 4πr
r ， G 0
G (r , r0 ) u (r0 ) G (r , r0 ) f (r )dV (r ) ]dS T n
得
y0 u ( x0 , y0 ) π
( x)
( x x0 ) y
2 2 0
dx
或代入拉普拉斯方程的第一边值问题的解公式
G (r , r0 ) u (r ) (r0 ) ]dS0 n 0 得到
10.4 用电像法确定格林函数
用格林函数法求解的主要困难还在于如何确定格林函数本身 一个具体的定解问题，需要寻找一个合适的格林函数 为了求解的方便，对一些具体问题我们给出构建格林函数的方法 一、电像法定义 考虑一个具体的物理模型：设在一接地导体球内的 M 0 点 放置一个单位正电荷，求在体内的电势分布，并满足边界条件为零
边界外法线方向为负 y 轴，故有
y0 y0 y0 G G 1 1 1 | | y 0 = 2 2 2 2 2 n y 2π ( x x0 ) y0 π ( x x0 ) y0 π ( x x0 )2 y0
代入到拉普拉斯第一边值问题解的公式，拉普拉斯方程的自由 项 ,则由 f 0
因此
c0
,故得到
对于三维无界球对称情形的格林函数可以选取为
G(r , r0 )
代入
1 4π | r r0 |
u (r ) G(r , r0 ) f (r0 )dV0
T0
得到三维无界区域问题的解为
f (r0 ) 1 u (r ) dV0 T 0|r r | 4π 0
上式正是我们所熟知的静电场的电势表达式
二、二维轴对称情形
用单位长的圆柱体来代替球．积分在单位长的圆柱体内进行，即
G(r ,0)dV (r )dV
T T
因为
(r )dV 1
T
G(r ,0)dV G(r ,0)dV
T T
一、三维球对称
对于三维球对称情形，我们选取
r0 0
G(r , r0 ) (r - r0 )
两边在球内积分
G(r ,0)dV (r )dV
T T
(r )dV 1
T
利用高斯定理得到
T
G(r ,0)dV G(r ,0)dV
G ( x, y | x0 , y0 )
( x x0 ) 2 ( y y0 ) 2 1 ln[ ] 2 2 4π ( x x0 ) ( y y0 )
据上述物理模型可求解下列定解问题 例1 定解问题：
u xx u yy 0, ( y 0) u | y 0 ( x)
根据第一类边值问题的解公式得到
u ( x, y )
解： 根据第一边值问题，构建的格林函数满足
2G Gxx Gyy ( x x0 ) ( y y0 )
G | y 0 0
( x0 , y0 ),( x0 , y0 ) 处放置于一个正和一个负的点电荷（或点源）
构建格林函数为
( x x0 )2 ( y y0 )2 1 G ( x, y | x0 , y0 ) ln[ ] 2 2 4π ( x x0 ) ( y y0 )
对于第一类边值问题，其格林函数可定义为下列定解问题的解
G(r , r0 ) (r - r0 ) G(r , r0 ) | 0
为了满足边界条件：电势为零，所以还得在边界外像 点（或对称点）放置一个合适的负电荷，这样才能使这两 个电荷在界面上产生的电势之和为零
这方法是基于静电学的镜像原理来构建格林函数，所 以我们称这种构建方法为电像法（也称为镜像法）．
令积分常数为0，得到
G(r ,0)
1 1 G(r ,0) ln c 2π r
1 1 ln 2π r
因此二维轴对称情形的格林函数为
G(r , r0 ) 1 1 ln 2π | r r0 |
得到二维无界区域的解为
u (r ) 1 1 f ( r )ln dS0 0 S 2π 0 | r r0 |
二、 上半平面区域第一边值问题的格林函数构建
拉普拉斯方程的第一边值问题求解
物理模型：若在
M 0 ( x0 , y0 ) 处放置一正单位点电荷
则虚设的负单位点电荷应该在 M1 ( x0 , y0 ) 于是得到这两点电荷在 xoy 的上半平面的电位分 布．也就是本问题的格林函数，即为
G (r , r0 ) 1 1 1 1 ln ln 2π | r r0 | 2 π | r r1 | 1 1 1 1 ln ln 2π ( x x0 ) 2 ( y y0 ) 2 2π ( x x0 ) 2 ( y y0 ) 2
g ( x0 ) y u ( x, y ) dx0 2 2 π ( x x0 ) y
称为上半平面的拉普拉斯积分公式．
三、 泊松方程的第一边值问题求解
例2
定解问题：
u xx u yy f ( x, y ) u ( x,0) ( x)
( <x <+, y 0) ( <x<+, y 0)
S
G(r ,0) dS
由于
G
源自文库
G er , G r
只是垂直于轴，且向外的分量，所以上式在
圆柱体上、下底的面积分为零，只剩下沿侧面的积分，即
G r rddz T (r )dV 1
选取的圆柱的高度为单位长，则很容易得到下面的结果
G 1 r 2πr