数学物理方法--格林函数法

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第十二章 格林函数法
第一节 泊松方程的格林函数法 第二节 用电像法求格林函数
11
12.1 泊松方程的格林函数法
有源问题
1. 源问题 例 静电场
a.无界空间
定解=通解+边界条件 求通解=积分

定解=积分+边界条件 (格林函数法)
r r'
q 4
r
r'
r 处静电场
(r)
r
1 r '
u0 (r , r ') G(r , r ')
G(r , r0 ) 4 (r r0)
G(r , r0 ) 0
第一边值问 题格林函数
r r0
0
导体球内有一个点电荷 ,导体接 地。求球内电势。 电荷的存在,在导体上感应了电荷。 球内的电势为自由电荷和感应电荷电势之和。
将感应电荷的电势由一 “电像电荷”的电势表示
10
泊松方程与第一类边界条件,构成第一边值问题(狄里希利问题) 泊松方程与第二类边界条件,构成第二边值问题(诺依曼问题) 泊松方程与第三类边界条件,构成第三边值问题
6
4. 泊松方程的基本积分公式
点源泊松方程 G(r , r ') 4 (r r ') 单位负电荷在
(Gu uG)dV
1
4
G(r , r0)

f n
dS.
8
第三边值问题
u(r
,
r0
)



[
u n
u]


f
[
G n

G]


0
第三边值问 题格林函数
u(r0 )

1
4
G(r , r0)(r )dV T

1
4


fG(r , r0 )dS.
9
12.2 电像法求格林函数
R r2 r22 r2 2rr2 cos
11
例 u 0 (z 0)

u ra f (x, y)

M0 (x0, y0, z0 )
半空间第一边值问题
计算格林函数:
M (x, y, z)
按电磁学思维模式, 应当引入镜像电荷 表示平面(z=0)上 的感应电荷。
镜像电荷的作用为使 平面(z=0)上的电势 为零。显然,这个电荷 位于相对于平面(z=0) 对称的几何点,且有 相反的电量。
3
感应电荷 是边界问题
2. 格林公式
第一格林公式:
区域 T,边界
定解=通解+边界条件 求通解=积分

定解=积分+边界条件 (格林函数法)
T



设 u(r ) 和 v(r ) 在 T 中具有连续二阶导数,
在 上有连续一阶导数。由高斯定理

uv dS (uv)dV
0 z0
1
0 z0
4 [(x x0 )2 ( y y0 )2 (z z0 )2 ]3/2 4 [(x x0 )2 ( y y0 )2 (z z0 )2 ]3/2

1 2
[( x

x0 )2

(y

z0 y0 )2

(z

z0 )2 ]3/ 2
u(x, y)
p
M (r)
o
M0 (r0 )
如右图,当导体外 M1 处有电荷 40q 时,镜像电荷
将在球内M0 处。
M1(r1)
像电荷的大小以及位置:

4 0 q
a r1
a2 r2 r1
G(r , r1)
1 R r1
a r1
1 R r2
其中
R r1 r12 r2 2rr1 cos



(u
v n

v
u )dS n


T
(uv

vu)dV
法向导数
5
3. 边值问题 边界条件
泊松方程
u
[
u n

u]

()
() 定义在
0, 0 0, 0
第一类边界条件 第二类边界条件
0, 0 第三类边界条件
格林函数代表一个点源所产生的场,对点源产生的场 利用叠加定理则可知道任意源分布产生的场
2
r r'
b.有界空间 边界上可能出现感应电荷
(r ')
r
r'
r 处静电场是源电荷与感应电荷
的电势之和。

感应电荷是源电荷的结果。

在一般的定解问题中,是要求解满足一定的边界条 件和初始条件的解。相应的GREEN函数的表达式会 复杂一些。因为此时的点源产生的场还是受边界条 件或初始条件的限制。


G(r

r
')
u(r n
)

u(r
)
G(r n
r
')
ds
7
5. 边值问题的格林函数
还需知道点源泊松方程度解的边界条件。
第一边值问题(狄里希利问题)
u(r , r ')
G(r , r ') 0
u f
第一边值问 题格林函数
u(r0 )

1
4
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

M1(x1, y1, z1)
M1(x1, y1, z1) M1(x0, y0,z0 )
G(r , r0 )


1 4
[( x

x0 )2
(y

1 y0 )2
(z

z0 )2 ]1/ 2

1 4
[( x

x0 )2

(y

1 y0 )2

(z

z0 Leabharlann Baidu2 ]1/ 2
12
1
T
G(r0, r ) f
(r0 )dV0



(r0
)
G(r0
n
,
r)
T


1


4

T
u
(r

r
')dV
由Green公式:
Gu
V
uG(r , r
')dv




G
u n
u
G n
ds
以及 (r r '函) 数的性质: u(r ) (r r ')dv u(r )
V
u(r)
1
4
G(r
V
r ')dv 1 4

T
u vdV uvdV
T
T
4
第二格林公式:
交换
u(r )

v(r )


vu dS v udV vudV

T
T
与上式相减

(uv vu) dS (uv vu)dV

T
即 n
相关文档
最新文档