第四章 Laplace方程的格林函数法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第四章 Laplace 方程的格林函数法
在第二、三两章,系统介绍了求解数学物理方程的三种常用方法—分离变量法、行波法与积分变换法,本章来介绍Laplace 方程的格林函数法。先讨论此方程解的一些重要性质,在建立格林函数的概念,然后通过格林函数建立Laplace 方程第一边值问题解的积分表达式。 §4.1 Laplace 方程边值问题的提法
在第一章,从无源静电场的电位分布及稳恒温度场的温度分布两个问题推导出了三维Laplace 方程
2
2
2
2
2
2
2
u u u u u x
y
z
∂∂∂∇=∆≡
+
+
=∂∂∂
作为描述稳定和平衡等物理现象的Laplace 方程,它不能提初始条件。至于边界条件,如第一章所述的三种类型,应用得较多的是如下两种边值问题。
(1)第一边值问题 在空间(,,)x y z 中某一个区域Ω的边界Γ上给定了连续函数f ,要求这样一个函数(,,)u x y z ,它在闭域Ω+Γ(或记作Ω)上连续,在Ω内有二阶连续偏导数且满足Laplace 方程,在Γ上与已知函数f 相重合,即
u
f
Γ
= (4.1)
第一边值问题也称为狄利克莱(Dirichlet )问题,或简称狄氏问题,§2.3中所讨论过的问题就是圆域内的狄氏问题。
Laplace 方程的连续解,也就是所,具有二阶连续偏导数并且满足Laplace 方程的连续函数,称为调和函数。所以,狄氏问题也可以换一种说法:在区域Ω内找一个调和函数,它在边界Γ上的值为已知。
(2)第二边值问题 在某光滑的闭曲面Γ上给出连续函数f ,要求寻找这样一个函数(,,)u x y z ,它在Γ内部的区域Ω中是调和函数,在
Ω+Γ
上连续,在Γ上任一点处法向导数
u n
∂∂存在,并且等于已知函数f
在该点的值:
u f
n
Γ
∂=∂ (4.2)
这里n 是Γ的外法向矢量。
第二边值问题也称纽曼(Neumann )问题。
以上两个问题都是在边界Γ上给定某些边界条件,在区域内部要求满足Laplace 方程的解,这样的问题称为内问题。
在应用中我们还会遇到Dirichlet 问题和Neumann 问题的另一种提法。例如,当确定某物体外部的稳恒温度场时,就归结为在区域Ω的外部求调和函数u ,使满足边界条件u
f
Γ
=,这里Γ是Ω的边界,f
表示物体表面的温度分布。像这样的定解问题称为Laplace 方程的外问题。
由于Laplace 方程的外问题是在无穷区域上给出的,定解问题的解是否应加以一定的限制?基于电学上总是假定无穷远处的电位为零,所以在外问题中常常要求附加如下条件:
lim (,,)0(r u x y z r →∞
==
(4.3)
(3)狄氏外问题 在空间(,,)x y z 的某一闭曲面Γ上给定连续函数
f
,要找出这样一个函数(,,)u x y z ,它在Γ的外部区域'Ω内调和,在
'Ω+Γ
上连续,当点(,,)x y z 趋于无穷远时,(,,)u x y z 满足条件(4.3),
并且它在边界Γ上取所给的函数值
u
f
Γ
= (4.4)
(4)纽曼外问题 在光滑的闭曲面Γ上给定连续函数f ,要找出这样一个函数(,,)u x y z ,它闭曲面Γ的外部区域'Ω内调和,在'Ω+Γ上连续,在无穷远处满足条件(4.3),而且它在Γ上任一点的法向导数u n
∂∂存在,满足
u f
n Γ
∂='
∂ (4.5)
这里'n 是边界曲面Γ内内法向矢量。
重点讨论内问题,所用的方法也可以用于外问题。 §4.2 格林公式
为了建立Laplace 方程解的积分表达式,需要先推导出格林公式,而格林公式则是曲面积分中高斯公式的直接推论。
设Ω是以足够光滑的曲面Γ为边界的有界区域,(,,)P x y z ,
(,,)Q x y z ,(,,)R x y z 是在Ω+Γ
上连续,在Ω内具有一阶连续偏导数的任
一函数,则成立如下的高斯公式
(
)[cos(,)cos(,)cos(,)]P Q R dV P x Q y R z dS
x
y
z
Ω
Γ
∂∂∂+
+
=
++∂∂∂⎰⎰⎰
⎰⎰
n n n (4.6)
其中d V 是体积元素,n 是Γ的外法向矢量,d S 是Γ上的面积元素。
下面来推导公式(4.6)的两个推论。
设函数(,,)u x y z 和(,,)v x y z 在Ω+Γ上具有一阶连续偏导数,在Ω内
具有连续的所有二阶偏导数,在(4.6)中令
,,v v v P u
Q u
R u
x y z
∂∂∂===∂∂∂
则有
()(
)u v u v u v v u v dV dV u
dS
x x
y y
z z
n
Ω
Ω
Γ
∂∂∂∂∂∂∂∆+
+
+
=
∂∂∂∂∂∂∂⎰⎰⎰
⎰⎰⎰
⎰⎰
或
()v u v dV u
dS u vdV
n
Ω
Γ
Ω
∂∆=
-
∇∇∂⎰⎰⎰
⎰⎰
⎰⎰⎰
(4.7)
(4.7)式称为第一格林公式。
在公式(4.7)中交换,u v 位置,则得
()u v u dV v
dS v udV
n
Ω
Γ
Ω
∂∆=
-
∇∇∂⎰⎰⎰
⎰⎰
⎰⎰⎰
(4.8)
将(4.7)与(4.8)式相减得到
()()v u u v v u dV u
v
dS
n
n
Ω
Γ
∂∂∆-∆=
-∂∂⎰⎰⎰
⎰⎰
(4.9)
(4.9)式称为第二格林公式。
利用格林公式可以推出调和函数的一些基本性质。 (ⅰ)调和函数的基本表达式
所谓调和函数的积分表达式,就是用调和函数及其区域边界Γ上的法向导数沿Γ的积分来表达调和函数在Ω内某一固定点,现在我们就来求调和函数在这点的值。为此,构造一个函数
1v r
=
=
(4.10)
函数1r
除点0M 外处处满足Laplace 方程,它在研究三维Laplac 方程中起着重要的作用,通常称它为三维Laplace 方程的基本解。由于1
v r
=在
Ω
内有奇异点0M ,我们作一个以0M 为中心,以充分小的正数ε为半