10.2边值问题的格林函数法
数学物理方法第十章 格林函数法

上式给出了泊松方程解的积分表达(biǎodá),但由于G(M,M0)未知 且不同边值条件也需做进一步的分析。
共二十六页
§10 格林函数(háns
2、泊松方程(fāngchéng)边值问题的积分公式
(A)第一类边界条件 0
由
边界条件变为 u 1 g(M ) f (M )
基本(jīběn)公式变为
这里(zhèlǐ)G就相当于 格林第二公式中的v
(G u u G )d (Gu uG)d
n
n
[u(M ) (M M0 ) G(M , M 0 )h(M )]d
若能由此式化简整理得到u(M),则一定(yīdìng)是方程(1)的解
共二十六页
§10 格林函数(hán
共二十六页
§10 格林函数(hánsh
显然,为了解决这一矛盾,或者修改格林函数所满足的方程
G(M , M0 ) (M M0 )
使之与边界条件
G 相0 容,
n
这就要引入所谓的广义格林函数方程;或者修改边界条件使之
与格林函数所满足的方程相容,这里不再详细讨论。
共二十六页
§10 格林函数(hánsh
(C)第三类边界条件 0, 0
积分变换法:无界区域(qūyù)的定解问题, 解一般为无穷积分
共二十六页
§ 10.1
函数(hánshù)
§10 格林函数(háns
共二十六页
2、定义(dìngyì)
(x)
0
x0 x0
(x)dx 1
更普遍的定义为
§10 格林函数(hánsh
—— 函数
(hánshù)
共二十六页
§10 格林函数(háns
u(M )
数理方程第四章 格林函数法

则 u(M 2 ) u(M1 ) 。以 M 2 为中心,以小于 d 的数为半径
在 内作球 k 2 ,在 k 2上u(M ) u(M 2 ) u(M1 ) ,…, n 次后,
点 N 一定包含在以某点 M n 为中心 ,半径小于 d 的球
kn 内 , 因而 u( N ) u(M n ) u(M1 ) , 由 N 的
性质1. 设 u(x, y, z) 是区域 内的调和函数,它在
上有一阶连续偏导数,则
udS n
0,
其中
,
n
是 的外法线方向。
证明 只要在Green公式中取 v 1即证。
注:此性质表明调和函数的法向导数沿区域边界的积分为零。 对稳定的温度场,流入和流出物体界面的热量相等,否则就 不能保持热的动态平衡,而使温度场不稳定。
3
下午9时12分
HUST 数学物理方程与特殊函数
第4章格林函数法
对二维拉普拉斯方程 u uxx uyy 0 ,其极坐标形式为:
2u r 2
1 r
u r
1 r2
2u
2
0
(4.1.2)
求方程(4.1.2)的径向对称解 u V (r) (即与 无关的解) ,则有:
d 2V dr 2
1 r
dV dr
任意性 ,就得到整个 上有 u( N ) u( M 1 ) ,这与 u 不为
常数矛盾.
10
下午9时12分
HUST 数学物理方程与特殊函数
第4章格林函数法
K1 M2 l M1 K2 M3
S1 S2
Kn N Mn Sn
图4.1
11
下午9时12分
HUST 数学物理方程与特殊函数
第4章格林函数法
格林函数法

通过格林公式,把静电边值问题与相应的格林 函数问题联系起来。 一般的处理方法,在物理学领域有着非常广泛 的应用
3
本节主要内容: 1. 格林函数——对应于给定问题的单位点源
的电势解; 2. 格林函数与泊松方程的解之间的关系; 3. 几种简单边界问题的格林函数形式。
10/20/2014
§5 格林函数法
1
几种方法的比较
1. 镜像法只适用于比较简单(点电荷)问题; 2. 分离变量法是精确求解的方法:除了几个高对
称的边界问题以外,一些实际问题往往难以求 解; 3. 多极展开法只适用于求远处的场(最后一节); 4. 格林函数方法
2
1
10/20/2014
格林函数方法: Green函数本身实际上是对应于给定问题所对
4
2
10/20/2014
几个基本公式:Ñ
1 r
=
-
r r3
,
高斯定理:
ò
E
×
dS
=
1 e0
i
Qi
空间一个单位点电荷的电场: E
=
4
1 e0
r r3
若点电荷处于闭合积分面内:
��������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������
数学物理方法--格林函数法

0, 0
泊松方程与第一类边界条件,构成第一边值问题(狄里希利问题)
泊松方程与第二类边界条件,构成第二边值问题(诺依曼问题) 泊松方程与第三类边界条件,构成第三边值问题
6
4. 泊松方程的基本积分公式
点源泊松方程
G(r , r ') 4 (r r ')
3
感应电荷 是边界问题
2. 格林公式 第一格林公式: 区域 T,边界
定解=通解+边界条件 求通解=积分
定解=积分+边界条件 (格林函数法)
u (r ) 和 设
T
在
上有连续一阶导数。由高斯定理 uv dS (uv)dV
T
v(r ) 在 T 中具有连续二阶导数,
7
5. 边值问题的格林函数
第一边值问题(狄里希利问题)
还需知道点源泊松方程度解的边界条件。
u
f
u (r , r ') G(r , r ') 0
第一边值问 题格林函数
1 u (r0 ) 4
1 G (r , r0 )r (r )dV 4 T
1
1
12.1
泊松方程的格林函数法
有源问题
定解=通解+边界条件 求通解=积分
1. 源问题 例 静电场 a.无界空间
定解=积分+边界条件 (格林函数法)
q 4
r'
r r ' r
r
处静电场
1 (r ) u0 (r , r ') G(r , r ') r r '
格林函数法

为第三边值问题的积分表示式
物理意义:右边第一个积分表示区域T中分布的源在r 点产生的场的总和;第二个积分代表边界上的状况对 r点场的影响的总和;两项积分中的格林函数相同。 说明泊松方程的格林函数是点源在一定的边界条件下 所产生的场。
对于拉普拉斯方程,f(r0)=0,因此可得拉普拉斯 方程第一边值问题的解
因此,我们可设想一个等效的点电荷,它位 于球外M1处,且在球面产生的电势与球内点电荷 在球面产生的电势相反。由物理学知识可知,该 设想的点电荷必位于OM0处的延长线上,如图所 示,并记:
OM r, OM0 r0
在∑ε 上的解,该解表示位于球心r=r0处的电量为ε0的 点电荷在半径为ε的球面上产生的电势,根据电磁学 知识,该电势为:
1
G(r, r0 ) 4
因此我们可得∑ε面上的积分
Ò
u(r)
G n
G
u(r) n
dS
Ò
u(r
)
n
(1
4
)
1
4
u(r) n
dS
Ò
u(
r
)
n
(1
4
)
1
4
u(r n
)
2d
(r r0 ) (x x0) ( y y0) (z z0)
格林函数的物理意义:在物体内部(T内)处放置 一个单位点电荷(或热源),而该物体的界面保持 电位为零(或温度为零), 那么该点电荷(或该点 热源)在物体内产生的电势分布(或稳定温度分 布),就是上述定解问题的解――格林函数。
格林函数互易定理: 格林函数代表r0处的点源在r处 所产生的影响,系统不变,则该影响等同于将移至r 处的该点源在r0处产生影响。故格林函数遵守如下 的互易定理:
高等数学-格林公式及其应用.ppt

l D1
O D2
x
1
2π
d
1 2π
π
20
2
l :4x2 y2 2
法二
l
ydx xdy 4x2 y2
l
ydx
2
xdy
1
2
ydx xd y
l
格林公式
D2是由l 所围区域
4x2 y2 2
所以 I 0 π
π.
1
2
1
2
(1
D2
(2)
π
2
1)dxdy
2
π
25
10.3 格林公式及其应用
Pdx Qdy
L
(L1, L2, L3对D来说为正方向)
8
10.3 格林公式及其应用
(3) 对复连通区域证明:
对若复区连域通不区止域由D一, 格条林闭公曲式线
的右所曲端围线应成积 包.添分 括加,沿且直区边线域界段D的的A方全B向,部CE对边.区界 G D
域则DD来的说边都界是曲正线向由. AB, L2 , BA,
2π 0
格林公式
sin d(
2
(Q P )dxdy D1 x y 0
cos ) cos d(
2
2
0 sin
)
24
10.3 格林公式及其应用
l
ydx xdy 4x2 y2
2π
sin
d(
2
cos
)
2
cos
d(
sin
)
0
2
2 0
π
2
2
sin
2
2
2
2
cos2
d
y L: x2 y2 4
数学物理方程第10讲 格林函数法 叶葱

M(x,y,z)
v u (u v u)dV (u n v n )ds
现在的问题是, V(x,y,z)不包含M0这一点!!!! 所以运用公式时我们要挖去M0点(奇异点)
如何去除M0点??
最简单的,以M0为中心, ɛ 为半径作一个球面, 球面为Ƭɛ,球体积为Kɛ,挖去这样一个球。
1 u(M 0 ) 4 1 rMM 0 n 1 u ( M ) )ds rMM 0 n
(u(M )
我们要求区域内一点M0处的u, 要知道这个函数在区域边界Ƭ上的值 以及在Ƭ上的法向导数的值
1 r 1 u )ds 4u 4 ( u ) 0 根据 (u n r n n
0, lim u u(M 0 )
1 u(M 0 ) 4 1 rMM 0 n 1 u ( M ) )ds rMM 0 n
(u(M )
调和函数的积分表达式
M0(x0,y0,z0)
M(x,y,z)
考虑球面Ƭɛ上,即M点在球面,此时r=ɛ
1 1 r r 1 1 n r r2 2
1 r ds 1 u n 2
uds r 1 u )ds ? (u n r n
2 2
第二格林公式
现在我们求解u(x,y,z)
u0
2
Dirichlet 问题
u
f ( x, y , z )
求出调和函数 的积分表达式
首先构造一个辅助函数
M0(x0,y0,z0) r
M(X,Y,Z)
1 1 v( x, y, z) 2 2 2 r ( x x0 ) ( y y0 ) ( z z0 )
数学物理方法--格林函数法

G(r , r0)r(r )dV T
1
4
f
G(r , r0 ) dS. n
第二边值问题(诺依曼问题)
u(r , r ')
u n
f
第二边值问 题格林函数
G(r , r ')ห้องสมุดไป่ตู้n
0
u(r0 )
1
4
G(r , r0)(r )dV T
(u
v n
v
u )dS n
T
(uv
vu)dV
法向导数
5
3. 边值问题 边界条件
泊松方程
u
[
u n
u]
()
() 定义在
0, 0 0, 0
第一类边界条件 第二类边界条件
0, 0 第三类边界条件
3
感应电荷 是边界问题
2. 格林公式
第一格林公式:
区域 T,边界
定解=通解+边界条件 求通解=积分
定解=积分+边界条件 (格林函数法)
T
设 u(r ) 和 v(r ) 在 T 中具有连续二阶导数,
在 上有连续一阶导数。由高斯定理
uv dS (uv)dV
p
M (r)
o
M0 (r0 )
如右图,当导体外 M1 处有电荷 40q 时,镜像电荷
将在球内M0 处。
M1(r1)
像电荷的大小以及位置:
4 0 q
a r1
格林函数法(课堂PPT)

林函数
25.06.2020
.
12
计算电磁学基础
• 性质:
– 设数学物理方程为 L u(x) = f (x) – 而格林函数方程为 L G(x) =δ(x-x’) – 在相同的齐次定解条件下 – 因为: f(x) =∫f (x’)δ(x-x’) dx’ – 所以: u(x) =∫f (x’) G(x-x’) dx’
VPgQPgQdV Ñ SPQgdS
• 若将第一格林定理相减,即得矢量第二格林定理
VQgPPgQdV Ñ SPQQPgdS
25.06.2020
.
27
计算电磁学基础
• 格林定理说明区域中的场与边界上的场之间的关系。因 此,利用格林定理可以将区域中场的求解问题转变为边 界上场的求解问题。
• 此外,格林定理说明了两种标量场或矢量场之间应该满 足的关系。因此,如果已知其中一种场的分布特性,即 可利用格林定理求解另一种场的分布特性。
(x)Q (xx) Q (xx )0 , (x≠x’点)
Q (xx)dVQ ,(积分区域V包含x=x’点)
V
25.06.2020
.
6
计算电磁学基础
2、 格林函数引入
• Green函数是与理想点源相联系的。
– 具体地说,Green函数是理想点源在给定边界条件下微分方程的 解答。
• 用Green函数求解电磁场是场论中的重要方法之一。
(x)VG (x,x)(x)dV 0S [G (x ,x ) n (x ) n G (x ,x )d S ]
这就是用Green函数求解静电问题的一种形式解。
25.06.2020
.
30
计算电磁学基础
在第一类边值问题中,格林函数满足边界条件
格林函数方法

2 0 2
P ,它的坐标为
4 0 2
(它在 OP 连线上,题中b对应这里的
2 0 2
2 R0 x 2 R R 2
0
R R R R 2 P P r x x R 2R cos R R R R0 Q R0 R02 R02 ∵ Q 1 Q (b R ) R R a R
(3)球外空间的格林函数
P’
P
设点电荷Q = 1 坐标为 P ( x, y , z )
观察点为 P( x, y, z )
R x
R x
x2 y2 z 2
x 2 y 2 z 2
R 相当于题中的 a ) R0 R(
PP r x x R 2 R 2 2RR cos
(x)
S
解法: (1)先求第一类边值问题的格林函数
1 G ( x , x ) ( x x )
2
G s 0
0
(2)
(2). 把(2)的解(格林函数)代入下式即可:
( x ) G ( x , x ) ( x )dV 0 ( x ) G ( x , x )dS V S n
§2.5
内容提要
格林函数方法
一、格林函数
二、用格林函数求解一般的边值问题
机动
目录
上页
下页
返回
结束
1. 处于 x 点上的单位点电荷的电荷分布密度: ( x ) ( x x )
回忆:点电荷密度的
函数表示
V
( x )dx ( x x )dV 1
用格林定理来求解静电边值的方法——格林函数

用格林定理来求解静电边值的方法——格林函数 1.什么是格林函数:在r 'ρ处有一点电荷q ,则电荷密度可写为: ⎩⎨⎧'=∞'≠='-=rr r r r r q ρρρρρρ0)(δρ,(1) 该电荷密度激发的空间的势满足的方程为:)()(2r r q r '--=∇ρρρδεϕ,(2) ∴-∇='-=•∇,),(ϕδεE r r q E ρρρρΘ定义有一个负号。
(3) 同理,处于r 'ρ处的单位电荷的电荷密度为)()(r r r '-=ρρρδρ (4)该单位电荷密度激发的空间的势满足的方程为:)(1)(2r r r '--=∇ρρρδεϕ,(5) 定义一个函数——格林函数,用),(r r G 'ρρ来表示,且满足)(1),(2r r r r G '--='∇ρρρρδε。
(8) 显然格林函数的物理意义为在r 'ρ处的一个单位电荷在空间r ρ处所激发的电势。
显然(8)式对于r r 'ρρ和有对称性,故也可以看作是r ρ处单位电荷在空间r 'ρ处所激发的电势。
由于空间电势分布有两种边界条件,分别为:第一类边界条件:0|),(='s r r G ρρ。
确定了边界上的电势分布(将一大的电势为零的导体与之接触)(9)第二类边界条件:εS n G s 1|-=∂∂。
(10) 确定了边界上的场强分布,也即电荷分布(根据)(1r r E '-=•∇ρρδε,积分形式ετδε1)(1)(=''-=•⎰⎰d r r S d r E ρρρρρ,而E nG ρρ=∂∂-,Eρ的最简单的取法(之后详述)为εS 1(E ρ在边界表面不一定是均匀的) (7)由电势和电荷是共轭量,两个中只能确定一个。
2.格林定理:详细推到见第七讲课件3.2.2式的推导[]⎰∑⎰⎥⎦⎤⎢⎣⎡∂∂-∂∂=∇-∇iidS n n d ψφφψτψφφψ22(11)左边是对所有边界面包括的空间积分,右边是对所有边界面积分(求和),其中对n 的微分代表在该面上求被微分函数的梯度。
《格林函数方法》课件

对于流体动力学中的无界问题,例如流体在大气压下的流动,格林函数方法可以提供一种 有效的求解方法,通过引入适当的边界条件来处理无界问题。
流体动力学中的非线性问题
格林函数方法也可以用于求解流体动力学中的非线性问题,例如流体动力学的非线性波动 和湍流等。
结构力学问题中的应用
弹性力学问题
在结构力学中,弹性力学问题是 常见的,格林函数方法可以用于 求解弹性力学中的各种问题,例 如弹性体的应力、应变和位移等 。
结构稳定性问题
结构稳定性问题也是结构力学中 的重要问题,格林函数方法可以 用于求解结构的稳定性,例如结 构的临界载荷和失稳模态等。
结构动力学问题
在结构动力学中,格林函数方法 可以用于求解结构的动态响应和 振动模态等,例如结构的振动频 率、阻尼和响应等。
探究格林函数方法在污染物扩散、水文模型等 领域的应用。
在材料科学中的应用
研究格林函数方法在材料电磁性能、光学性能等领域的应用。
感谢您的观看
THANKS
粒子动力学问题中的应用
描述粒子运动
在粒子动力学问题中,格林函数可以用来描述粒子的运动轨 迹和行为。通过格林函数,可以将粒子动力学问题转化为求 解微分方程或积分方程的问题,从而得到粒子的运动轨迹和 性质。
04
格林函数在工程问题中的应 用
流体动力学问题中的应用
流体力学中的波动和散射问题
格林函数方法可以用于求解流体力学中的波动和散射问题,例如声波在流体中的传播、波 动在流体界面上的反射和折射等。
格林函数方法的优势与局限性
优势
能够将复杂的微分或积分方程问题转 化为简单的代数方程问题,计算效率 高,适用于大规模问题求解。
局限性
数学物理方程第四章 格林函数法

为边界的有界连通区域,u(x, y, z)在 上有连续
的一阶偏导数,在 内调和,定点 M 0 (x0 , y0 , z0 ) , r 为定点M 0到变点 M (x, y, z) 距离: 则有
u(M0 )
1
4
1 [ r
u n
u
(1)]ds n r
(2.9)
故不提初始条件!只给出边界条件就可以. 下面看边界条件的提法.
(1) 第一边值问题(狄利克雷(Dirichlet)问题)
设方程(1.1)的空间变量(x, y, z) , 为 R3的开区域。如果
u(x, y, z)满足方程(1.1),且在 边界 上直接给定了u(x, y, z)
的具体函数形式 f (x, y, z),即
u(x, y, z) f (x, y, z)
(1.2)
则称问题(1.1)~(1.2)为拉普拉斯第一边值问题或狄利克雷
(Dirichlet)问题,u(x, y, z) 为此问题的解。
2u 2u 2u
u
x 2
y 2
z 2
0
u( x, y,z) f ( x, y,z),
u, v互 换
v
u v u v u v
( uv )dV
u
n
ds
(
x
x
y
y
z
z
)dV
(2.2)
u
u v u v u v
(vu)dV
v
n
ds
(
x
x
y
y
z
数学物理方法第十章_格林函数法讲解

=
1 2π
(x
y0 x0 )2
y02
1 π
(x
y0 x0 )2
y02
1 π
(x
y0 x0 )2
y02
代入到拉普拉斯第一边值问题解的公式,拉普拉斯方程的自由 项 ,则由 f 0
u(r0
)
T
G(r,
r0 )
f
(r)dV
(r)
G(r, n
r0
) ]dS
得
因为
T (r)dV 1
T G(r,0)dV T G(r,0)dV S G(r,0) dS
由于
G
G r
er
,G
只是垂直于轴,且向外的分量,所以上式在
圆柱体上、下底的面积分为零,只剩下沿侧面的积分,即
G r
rddz
T
(r)dV
2π
1 S0 f (r0 ) ln | r r0 | dS0
10.4 用电像法确定格林函数
用格林函数法求解的主要困难还在于如何确定格林函数本身 一个具体的定解问题,需要寻找一个合适的格林函数 为了求解的方便,对一些具体问题我们给出构建格林函数的方法
一、电像法定义 考虑一个具体的物理模型:设在一接地导体球内的 M 0 点
故有
S
G r
r2
sin d d
T
G(r , 0)dV
1
使上式恒成立,有 4πr2 G(r,0) 1 r
G(r,0) 1 c 4πr
r ,G 0 因此 c 0 ,故得到
对于三维无界球对称情形的格林函数可以选取为
格林函数法

r
0
dS
(5.2.8)
式(5.2.8)称为泊松方程的基本积分公式。但是它
的物理解释很困难,因此我们根据前面的格林函数互易定
理,并利用跟林函数的对称性,将上式改为:
u r0
G r0,r
u r
G
T0
r0,r
f
r0
dV0
0
G
r0,r
u r0 n0
n0 dS0(5.2.9)
G r0,r
u r
G r,r0
T0
f
r0
dV0 r0
0
n0 dS0
(2)第二类边值问题
(5.2.13)
14
对应下列格林函数的解:
u r f r
u
n
|
rp
(5.2.14)
G r, r0 r r0 G r, r0
n | 0
(5.2.15)
代入基本积分公式可得第二类边值问题的解的积
G(r, r0) G(r0, r)
上式表明,在位于r0处的脉冲(或点源)在一 定边界条件下在r处产生的影响(或产生的场), 等效于把脉冲(或点源)移至r处在同样边界条件 下在r0处算产生的影响(或场),即物理场的互 易性。
10
根据第二格林公式,得到:
u
r
G n
G
u r
n
dS
T
u rG Gr
u
r
0
T
G
r,
r
0
f
r
dV
1
r
G
r,
n
r0
dS
(5.2.19)
利用格林函数的互易性可得到互易后的解的积
格林函数法

( )
24
令积分常数c=0,得到得到二维球对称情形的格 林函数为
vv 1 1 G r, r0 = ln v v 2π r −r0
( )
(5.3.9)
代入(5.3.1)得到二维无界区域的解为
v 1 v 1 u r = ∫∫ f r0 ln v v dS0 2π S0 r−r0
4π r
(5.3.6)
22
当r→∞时,G→0,因此c=0,故得到
v 1 G r, 0 = 4π r
( )
因此三维无界球对称的格林函数可选取
vv G r, r 0 =
(
)
1 v v 4π r − r 0
(5.3.7)
代入(5.3.1),得到三维无界区域问题的解为
v 1 u r = ∫∫∫ 4π T0
8
那么该点电荷(或热源)在物体内产生电视分布 (或稳定温度分布),这就是定解问题的解—格 林函数。 引入格林函数的目的: 1、解的形式便于理论分析和研究; 2、以同一的形式研究各类定解问题; 3、对于线性问题,格林函数一但求出,就可以 算出任意源的场,关键是求格林函数。
9
因为格林函数G(r,r0)代表r0处的点源在r处所产生的 影响(或产生的场),所以它只能是距离|r-r0|的函 数,它应该遵守如下的互易定理:
v v ∆u (r ) = − f (r ) ∂u αu + β ∂n ∑ v = ϕ r∑
( )
(5.2.3)
为了求上述定解问题,我们需要找一个能表示点 源密度分布的函数,而δ函数则正是描述一个密度 分布的函数。因此定义格林函数G(r,r0)它满足 如下定解问题:
数学物理方法第10章格林函数法

2
格林函数,又称为点源影响函数,是数学物理方程中的
一个重要概念,也是求解各类定解问题的另一种常用方法。
若已知点电荷(点 源)产生的场(边 界无限远,无初始 条件) 积分得到
Uq
任意带电体(任意 源)产生的场(边 界无限远,无初 始条件) 任意源在相同初 始和边界条件下 产生的场
19
1 4 r r0
R 0 r0 R2 4 0 r 2 r0 r0
1 R 4 r r0 r0
1 R2 4 r 2 r0 r0
20
G r0 ; r u r r0 dS0 G r0 ; r f r0 dV0 n T
u r f r u u r n
为了求解上面定解问题,我们必须定义一个与此定解问题相应的 格林函数 G(r , r0 ) 它满足如下定解问题,边值条件可以是第一、二、三类条件:
6
G(r , r0 ) (r r0 ) G [ G ] 0 n
13
对狄利克雷问题的格林函数应满足:
G r ; r0 r r0 G r;r 0 0
令 G G0 G1代入上述定解问题有
G0 G1 r r0 G0 G1 0
显然没有考虑边界的影响 (或者说对应着无界空间)
G r;r0 u r r0 dl0 G r;r0 f r0 dS0 n0 l S
1 G0 4 r r0
1 1 G0 ln c0 2 r r0
G1 0 G1 G0
21
10.2边值问题的格林函数法

dΩ
=
ds r2
=
sinθ
dθ
dϕ
-u(M
,
t)在以M
0为中心,r为半径的球面
sM0 r
上的平均值。
(2)显然
u(M
0
,
t0
)
=
limu
r →0
(r,
t0
)
Wuhan University
§10.2泊松方
二、积分公式-格林函数法 程的狄氏问题
1、(三维)狄氏积分公式 M , M 0 ∈τ
⎪⎧ Δ u = − h ( M ) (1) ⎪⎩⎨u σ = f ( M ) ( 2 )
⎧Δu = −h(M ) ⎩⎨u σ = g(M )
(1) 之解。
(2)
Wuhan University
一、格林公式
§10.2泊松方 程的狄氏问题
3、球面平均值公式
(1)定义
∫∫ u(r,t) = 1
u(M ,t)ds
π4 r 2 srM0
∫∫ = 1
u(M ,t)dΩ
4π srM0
r
•
M0
sM0 r
⎧ΔG = −δ (M − M 0 ) (3)
⎩⎨G σ = 0
(4)
则(3) → G = 1 (10.3.4),
4π r
r=
(x − x0 )2 ( y − y0 )2 (z − z0 )2
[(1) ⋅G − (3) ⋅ u]在[τ −τε ]中积分有 :
∫ ∫ ∫ [GΔu −uΔG]dτ
τ −τε
3、格林第二公式
=
∫σ
v
∂udσ
∂n
(4)
∫τuΔvdτ
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
本节作业
§10.2泊松方 程的狄氏问题
习题 10.2:1(2)
Wuhan University
Good-by!
§10.2泊松方 程的狄氏问题
Wuhan University
附:证明∫τ[uΔv− v源自u]dτ=∫σ(u
∂v ∂n
−v
∂u )dσ
∂n
G(M1, M2) = G(M2, M1)
⎧ΔG(M , M1) = −δ (M − M1)
⎧Δu = −h(M ) ⎩⎨u σ = g(M )
(1) 之解。
(2)
Wuhan University
一、格林公式
§10.2泊松方 程的狄氏问题
3、球面平均值公式
(1)定义
∫∫ u(r,t) = 1
u(M ,t)ds
π4 r 2 srM0
∫∫ = 1
u(M ,t)dΩ
4π srM0
r
•
M0
sM0 r
3、格林第二公式
=
∫σ
v
∂udσ
∂n
(4)
∫τuΔvdτ
−
∫τvΔudτ
=
∫σ
(u
∂v ∂n
−v
∂u )dσ
∂n
(5)
意义:(1) 将u, v, Δu, Δv的值与u, v, ∂v , ∂u 的边值联系起来 .
(2) u, v对称
∂n ∂n
(3) 若已知 v, Δv = 0, 及v σ , 则由格林公式可能求得
dσ 0
QG(M , M0) = G(M0, M )
Wuhan University
-狄氏积分公式
二、积分公式-格林函数法
§10.2泊松方 程的狄氏问题
2、狄氏积分公式的物理意义: 第一项:体内源产生的场的和。
第二项: 边界面上源产生的场的和。
3、(二维)狄氏积分公式 ⎧Δu(M) = −h(M) M ∈σ ⎩⎨u |l = f (M)
(4)
∫ ∫ σ
(G
∂u ∂n
−u
∂G )dσ
∂n
− u(M 0)
=
−
G(M
τ
, M 0)h(M
)dτ
将边界条件(2)(4)代入上式有:
∫ ∫ u(M 0 ) =
G(M
τ
,
M
0 )h(M
)dτ
−
σ
f (M ) ∂G dσ
∂n
∫ ∫ u(M ) =
G
τ
(M
,
M
0
)h(M
0
)dτ
0
−
σ
f
(M 0)
∂G ∂n0
⎧ΔG = −δ (M − M 0 ) (3)
⎩⎨G σ = 0
(4)
则(3) → G = 1 (10.3.4),
4π r
r=
(x − x0 )2 ( y − y0 )2 (z − z0 )2
[(1) ⋅G − (3) ⋅ u]在[τ −τε ]中积分有 :
∫ ∫ ∫ [GΔu −uΔG]dτ
τ −τε
= =
−h(M ) f (M )
,M
∈τ
(1) (2)
∫∫∫ ∫∫ u(M ) =
τ G(M , M 0 )h(M 0 )dτ 0 −
σ
f
(M
0
)
∂G ∂n0
dσ 0
(6)
⎧ΔG = −δ (M − M 0 ) (3)
⎩⎨G σ = 0
(4)
G(M , M0) = G(M0, M )
Wuhan University
Methods in Mathematical Physics
第十章格林函数法
Method of Green’s Function
武汉大学物理科学与技术学院
Wuhan University
第十章格林函数法
§10.2泊松方程的狄氏问题
Dirichlet Problem for the Poisoon’s Equations
G(M
τ
,M
0 )h(M
)dτ
Wuhan University
§10.2泊松方
二、积分公式-格林函数法 程的狄氏问题
1、(三维)狄氏积分公式
狄氏格林函数
⎪⎧ Δ u = − h ( M ) (1) ⎪⎩⎨u σ = f ( M ) ( 2 )
⎧ΔG = −δ (M − M 0 ) (3)
⎩⎨G σ = 0
设u( x, y, z), v( x, y, z)在τ 中具有连续的二阶 导数,在 τ 上具有连续的一阶导数 ,则有
如下格林公式:
Wuhan University
一、格林公式
§10.2泊松方 程的狄氏问题
∫ ∫ 2、格林第一公式
uΔvdτ + ∇u ⋅∇vdτ
τ
τ
=
∫σ
u
∂v ∂n
dσ
(3)
∫τvΔudτ + ∫τ∇u ⋅∇vdτ
⎧Δu = −h(M ) , M ∈τ (1)
⎩⎨u σ = g(M )
(2)
Wuhan University
一、格林公式 1、为何引入格林公式
§10.2泊松方 程的狄氏问题
(1)积分公式的起点是通过直接积分或分部积分将 未知函数从微分号下解脱出来
(2).我们要求解的三类数值 方程中均含有 Δ, 格林公式是将未知函数 ,从微分算符Δ下 解脱出来的工具.
Wuhan University
∫∫ ∫ u(M ) =
σ G (M , M 0 )h(M 0 )dσ 0 −
l
f
(M 0)
∂G ∂n0
dl 0
Wuhan University
三、小结:
§10.2泊松方 程的狄氏问题
1、
∫τuΔvdτ
− ∫τvΔudτ
=
∫σ
(u
∂v ∂n
− v ∂u )dσ
∂n
(5)
2、
⎪⎧ Δ u ⎪⎩⎨u σ
= τ −τε uδ (M − M0 )dτ −
Gh(M )dτ
τ −τε
对左边用格林第二公式有:
∫ ∫ ∫ (G ∂u − u ∂G )dσ
σ +σε ∂n
∂n
=
τ −τε uδ (M − M 0 )dτ −
Gh(M )dτ
τ −τ ε
∫ ∫ σ
(G
∂u ∂n
−u
∂G )dσ
∂n
−
u(M 0)
=
−
⎩⎨G(M , M1) σ = 0
(7) (8)
⎧ ΔG(M , M ⎩⎨G(M , M 2 )
2)
σ
= −δ
=0
(M
−
M2)
(9) (10)
∫ ∫∫τ [(7) ⋅G(M , M 2 ) − (9) ⋅G(M , M1)]dτ : ∫ ∫∫τ [G(M , M 2 ) ⋅ ΔG(M , M1) − G(M , M1) ⋅ ΔG(M , M 2 )]dτ
∫∫ ∫ = τ [−G(M , M 2 )δ (M − M1) + G(M , M1)δ (M − M 2 )]dτ
∫∫σ
[G(M
,
M
2
)
∂ ∂n
G
(
M
,
M
1
)
−G
(
M
,
M
1
)
∂ ∂n
G
(M
,
M
2
)
]
= −G(M1, M 2 ) + G(M 2 , M1)
(8), (10)
0 = −G(M1, M 2 ) + G(M 2 , M1) → G(M1, M 2 ) = G(M 2 , M1)
dΩ
=
ds r2
=
sinθ
dθ
dϕ
-u(M
,
t)在以M
0为中心,r为半径的球面
sM0 r
上的平均值。
(2)显然
u(M
0
,
t0
)
=
limu
r →0
(r,
t0
)
Wuhan University
§10.2泊松方
二、积分公式-格林函数法 程的狄氏问题
1、(三维)狄氏积分公式 M , M 0 ∈τ
⎪⎧ Δ u = − h ( M ) (1) ⎪⎩⎨u σ = f ( M ) ( 2 )