数学物理方程第四章 格林函数法
格林函数法
(14.2.12)
考虑到格林函数的齐次边界条件,由公式(14.2.9) 可得第一类边值问题的解
u (r0 ) G (r , r0 ) f (r )dV (r )
T
G (r , r0 ) n
dS
(14.2.13)
另一形式的第一类边值问题的解
u (r ) G (r , r0 ) f ( r0 )dV0 ( r0 )
T 中具有连续二阶导数,应用矢量分析的高斯定理
A S d
AdV =
T
divAdV (14.1.1)
T
单位时间内流体流过边界闭曲面S的流量
单位时间内V内各源头产生的流体的总量
将对曲面 的积分化为体积分
uv S uv )dV uvdV u vdV d (
T0
(14.3.1)
选取 u (r ) 和 G(r , r0 ) 分别满足下列方程
u (r ) f (r )
G(r , r0 ) (r - r0 )
(14.3.2) (14.3.3)
14.3.1 三维球对称
对于三维球对称情形,我们选取 对(14.3.3)式两边在球内积分
r0 0
(14.2.4)
(r r0 ) 代表三维空间变量的 函数,在直角坐标系中其形式为
(r r0 ) ( x x0 ) ( y y0 ) ( z z0 )
(14.2.4)式中
函数前取负号是为了以后构建格林函数方便
格林函数的物理意义【2】:在物体内部(T 内) r0 处放置一个单位点电荷,而该物体的界面保持电位为零, 那么 该点电荷在物体内产生的电势分布,就是定解问题(14.2.4)的解 ――格林函数.由此可以进一步理解通常人们为什么称格林函 数为点源函数.
数理方程第四章 格林函数法
则 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章格林函数法
第四章_拉普拉斯方程的格林函数法
边值问题的提法: 二. 边值问题的提法:
1)第一边值问题(Dirichlet问题/狄氏问题) ∇ 2u = 0, in Ω uΓ = f 数学解释: 在Ω内寻求一个调和函数u, 它在边界Γ上与已知
连续函数f 吻合,即u Γ = f .
2)第二边值问题( Neumann 问题/牛氏问题 ) ∇ 2 u = 0, in Ω (其中 n是 Γ 的外法向量, f 是 连续函数 ) ∂u = f ∂n Γ
数学解释: Ω内寻求一个调和函数,它在闭区域Ω上有一阶 在 连续偏导数,即u ∈ C 2 (Ω) I C1 (Ω),且在边界上满足边界条件。 ,且在边界上满足边界条件。
∂u ∫∫ u (M )dS + 4π a ∫∫ ∂n dS Γ Γ 1
∫∫ u (M )dS
Γ
(4)Laplace方程解的唯一性问题 ) 方程解的唯一性问题
定理:狄氏问题在C 2 (Ω) I C 1 (Ω)内解唯一,牛曼问题除相差一个
常数外解也是唯一确定的。
证明:
设u1 , u2为上述两类问题的解,则它们的差v = u1 − u2必是原问题的
第四章 拉普拉斯方程的格林函数法
第一节 第二节 第三节 第四节 拉普拉斯方程边值问题的提法 格林公式 格林函数 两种特殊区域的格林函数及狄氏 问题的解
格 林 函 数 法
格林函数:又称点源影响函数,是数学物理中的一个重要概念 格林函数 又称点源影响函数, 又称点源影响函数
格林函数代表一个点源在一定的边界条件和初始条件下所产生的 格林函数代表一个点源在一定的边界条件和初始条件下所产生的 知道了点源的场就可以用迭加的方法计算出任意源所产生的场。 场,知道了点源的场就可以用迭加的方法计算出任意源所产生的场。 知道了点源的场就可以用迭加的方法计算出任意源所产生的场
格林函数法 数学物理方程
格林函数法
若L 一个带平滑系数的线性微分算子,当求解形如()L u f =的微分方程时,若对于任意的向量y 都存在广义函数()G x,y ,使得
[]()()L G δ=x x,y x-y
(此处下标x 表示L 作用于()G x,y 时将其当做以x 为自变量的广义函数,而y 为参数) 若再令
()()()d u G f =⎰x x,y y y
将上式代入()L u f =则有
[]()()d ()()d ()()d ()L G f L G f f f δ⎡⎤===⎣⎦
⎰⎰⎰x x,y y y x,y y y x -y y y x 故此时()u x 是微分方程()L u f =的解。
采用上述方法求解微分方程的方法称为格林函数法,广义函数()G x,y 也称为格林函数。
数学物理方法知识体系
数学物理方法所要解决的问题:求解(偏)微分方程
本学期学过的求解方法:变量分离法、积分变换法、格林函数法
变量分离法涉及知识点:傅里叶级数、函数的正交系、贝塞尔函数(Chap.2~Chap.5) 积分变换法涉及知识点:傅里叶变换、拉普拉斯变换、广义函数(Chap.7~Chap.9) 格林函数法涉及知识点:格林函数(Chap.10)
例题数量统计。
第四章 Laplace方程的格林函数法
第四章 Laplace方程的格林函数法第四章laplace方程的格林函数法在第二、三两章,系统介绍了求解数学物理方程的三种常用方法―分离变量法、行波法与积分变换法,本章来介绍laplace方程的格林函数法。
先讨论此方程解的一些重要性质,在建立格林函数的概念,然后通过格林函数建立laplace方程第一边值问题解的积分表达式。
§4.1laplace方程边值问题的提法在第一章,从无源静电场的电位原产及稳恒温度场的温度原产两个问题推论出来了三维laplace方程2u2u2uuu2220xyz2做为叙述平衡和均衡等物理现象的laplace方程,它无法加初始条件。
至于边界条件,例如第一章所述的三种类型,应用领域得较多的就是如下两种边值问题。
(1)第一边值问题在空间(x,y,z)中某一个区域?的边界?上给定了连续函数f,要求这样一个函数u(x,y,z),它在闭域(或记作?)上连续,在?内有二阶连续偏导数且满足laplace方程,在?上与已知函数f相重合,即u?(4.1)?f第一边值问题也称为狄利克莱(dirichlet)问题,或简称狄氏问题,§2.3中所讨论过的问题就是圆域内的狄氏问题。
1laplace方程的连续解,也就是所,具有二阶连续偏导数并且满足laplace方程的连续函数,称为调和函数。
所以,狄氏问题也可以换一种说法:在区域?内找一个调和函数,它在边界?上的值为已知。
(2)第二边值问题在某扁平的闭合曲面?上得出连续函数f,建议找寻这样一个函数u(x,y,z),它在?内部的区域?中就是调和函数,在上连续,在?上任一点处法向导数u存有,并且等同于未知函数f?n在该点的值:unf(4.2)这里n就是?的外法向矢量。
第二边值问题也称纽曼(neumann)问题。
以上两个问题都就是在边界?上取值某些边界条件,在区域内部建议满足用户laplace 方程的求解,这样的问题称作内问题。
在应用中我们还会遇到dirichlet问题和neumann问题的另一种提法。
格林函数法
两种边值问题: 两种边值问题:
第一边值问题
u |Γ = f .
这类问题也叫做狄利克雷 问题。 这类问题也叫做狄利克雷(Dirichlet)问题。 狄利克雷 问题
拉普拉斯方程的连续解,也叫调和函数, 拉普拉斯方程的连续解,也叫调和函数,所以 调和函数 狄利克雷问题也可以叙述为: 狄利克雷问题也可以叙述为:在区域 Ω 内找 一个调和函数, 上的值已知。 一个调和函数 它在边界 Γ 上的值已知。 第二边值问题 在光滑的闭曲面 Γ 上给出连续函数 f,寻找函数 , u(x,y,z) :在 Γ 的内部 Ω 是调和函数 在 ( Ω + Γ ) 是调和函数, 上连续, 上连续,在 Γ 上任一点法向导数存在并且等于 已知函数 f ,即: ∂u =f ∂n Γ 这类问题也叫做纽曼 纽曼(Neumann)问题。 问题。 这类问题也叫做纽曼 问题
在球面坐标下, 拉普拉斯方程为: 在球面坐标下, 拉普拉斯方程为:
1 ∂ 2 ∂u 1 ∂ ∂u 1 ∂2u =0 r + 2 sin θ + 2 2 2 2 r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂ϕ
球对称解 u=u(x,y,z)在以原点为中心的同一球面的 在以原点为中心的同一球面的 值为常数。 的函数:u=u(r)。 值为常数。u 仅为半径 r 的函数 。
(
)
两式相减, 两式相减 得
2 2
第二格林公式
∂v ∂u ∫∫∫ (u∇ v − v∇ u)dV = ∫∫ (u ∂n − v ∂n )dS Ω Γ
利用格林公式, 可以得到调和函数的一些性质: 利用格林公式, 可以得到调和函数的一些性质: 1) 纽曼内问题有解的必要条件 内的调和函数, 设 u 是以 Γ 为边界的区域 Ω 内的调和函数, 在 上有一阶连续偏导数, Ω + Γ上有一阶连续偏导数, 则在第二格林公式 为上述调和函数 调和函数, 中取 u 为上述调和函数,取 v ≡ 1, 有
数学物理方程第四章_格林函数
1 ⎧ ⎪∆G (r , r0 ) = − δ (r − r0 ) ε ⎨ ⎪G Γ = 0 ⎩
(4.3.7) (4.3.8)
以 G (r , r0 ) 乘式 (4.3.5), u (r ) 乘式 (4.3.7), 二式相减后在 Ω 上对 r 积分 ,以 dr 表示 r 点处的体积微元,有
∫
Ω
(G∆u − u∆G )dr = −
第 4 章 格林函数
在这一章里,我们介绍数学物理方程中另外一种常用的方法—格林函数法.从物理上看, 一个数学物理方程是表示一种特定的“场”和产生这种场的“源”之间的关系.例如,热传导 方程表示温度场和热源之间的关系,泊松方程表示静电场和电荷分布的关系,等等.这样,当源 被分解成很多点源的叠加时,如果能设法知道点源产生的场,利用叠加原理,我们可以求出同 样边界条件下任意源的场,这种求解数学物理方程的方法就叫格林函数法.而点源产生的场就 叫做格林函数. 4.1
⎧0, T ( x) = ⎨ ⎩∞,
x≠0 x=0
且
∫Байду номын сангаас
所以有
+∞
−∞
cρT ( x)dx = Q
T ( x) =
Q δ ( x) cρ
通过以上两个例题,我们对 δ ( x) 有了进一步的认识.如果将坐标平移 x0 ,即集中量 出现在点 x = x 0 处,则有
δ ( x − x0 ) = ⎨
且
⎧0, ⎩∞,
∫
= ∫ (u∆v)dΩ + ∫ gradu ⋅ gradvdΩ
Ω Ω
=∫u
Γ
∂v dS ∂n
或表示为
∫
Ω
(u∆v)dΩ = ∫ u
Γ
∂v dS − ∫ gradu ⋅ gradvdΩ Ω ∂n
第四章格林函数法课件
特点:除 M0(x0,y0,z0)点外,任一点满足Laplace方程。
同学们自己验证。
PPT学习交流
2
二维Laplace方程的基本解:
1
1
u(x,y)ln ln
rM M 0
(xx0)2(yy0)2
特点:除 M0(x0, y0) 点外,任一点满足Laplace方程。
同学们自己验证。 问题:基本解是否为整个区域内的解?
n
n
从而得证
1
1 1 u (M )
Ò u (M 0) 4
[u (M ) ( )
nrM M 0 rM M 0
]d S n
PPT学习交流
8
4 调和函数的基本性质
性质1:设 u ( x, y , z ) 在有界区域 内为调和函数,且在
上有一阶连续偏导数,则
Ò
u n
dS
0
证:令 v 1 将 u , v 代入第二Green公式即可。
uv
PPT学习交流
11
证明:用反证法
若在 内有 u v ,即 uv0 ,而在边界上 uv0 , 说明 u v 在内部可能取最大值。
推论2:狄利克莱问题 的解唯一。
u0, u f
(x, y,z)
证明:设 u 1 和 u 2 均为该问题的解,则 u u1 u2 满足
由极值原理, u 0
u0, u 0
于是
r rMM0
r2 MM0
2
乙 u n(rM 1 M 0)d S1 2 u d S1 24 2u4 u
乙 rM 1M0 u ndS1 u ndS4 u n
PPT学习交流
7
代入上式,得
Ò [u( 1)1u]dS4 u4 u0
Chapter 4. 格林函数法
20
为此,在第二Green公式
与调和函数的积分表达式相加
China University of Petroleum
21
其中
称为Laplace方程第一边值问题的格林函数/影响函数
China University of Petroleum
22
China University of Petroleum
16
在上一节的基础上,我们直接给出平面域上的一些结论。
China University of Petroleum
17
第一格林公式
第二格林公式
China University of Petroleum
18
4. Green公式的应用
China University of Petroleum
19
China University of Petroleum
9
为了建立三维拉普拉斯方程解的积分表达式,需要用到格林公
式。而格林公式实际上是高斯公式的直接推论。
China University of Petroleum
10
第一格林公式
第二格林公式
China University of Petroleum
11
4. Green公式的应用
利用格林公式可以推出调和函数的一些基本性质。
26
M
M 0
China University of Petroleum
27
China University of Petroleum
28
什么是电象法?
M
M0
China University of Petroleum
29
q M1
第四章格林函数法1
注1:当M 0取在区域之外或边界上,可用同样的方法导出公式
4 u ( M 0 ), 1 1 u [u ( ) ]dS 2 u ( M 0 ), n r r n 0,
M 0在内; M 0在上; M 0在外。
注2:若u不是调和函数,即2u F,只要u C 2 () C1 (), 我们可以得到类似公式
u u ds ds n n r R D
sin Rd 4 R 0 0 4
2
由牛曼内问题有解的必要条件知该问题无解。
3)平均值公式
定理3:设函数u(M )在区域内调和,M 0 ( x0 , y0 , z0 )为其中 任一点,a是以M 0为中心,以a为半径且完全落在内部 的球面,则下面平均值公式成立 1 u(M 0 ) udS 2 4 a a
P Q R ( ) dV Pdydz Qdzdx Rdxdy (1) x y z 其中取外侧位正向.
由两类曲面积分之间的关系得高斯公式的另一 种形式:
P Q R ( ) dv x y z ( P cos Q cos R cos )dS .
2 2
取u为调和函数,并假定且在上有一阶连续偏导数,v 1/ r则有
1 1 u r (u )dS 0 n r n
1 1 1 r r 注意到:在球面 上, 2 n r
1 1 r 因此可得 u dS 2 n 其中u
两式相减可得 v u (u v v u )dV (u v )dS n n
2 2
第二格林公式
二、调和函数的基本性质
1). 调和函数的积分表达式
定义:所谓调和函数的积分表达式,就是用调和函数及其 在区域边界上的法向导数沿的积分来表达调和函数在 内任一点的值。
4第四章格林函数法
则 u ( M 2 ) u ( M 1 ) 。以 M 2 为中心,以小于 d 的数为半径 在 内作球 k 2 ,在 k 2上 u ( M ) u ( M 2 ) u ( M 1 ) 点 N 一定包含在以某点 M n
c1 d 2 dV V (r ) 0 其通解为: (r ) c2 , (r 0, c1 , c2 为任意常数)。 r dr dr 1 1 若取 c1 , c2 0 ,则得到特解 V0 (r ) 4r ,称此解为 4
三维Laplace方程的基本解,它在研究三维拉普拉斯方程中 起着重要的作用. 对二维拉普拉斯方程 u uxx u yy 0,其极坐标形式为:
数学物理方程与特殊函数
第4章格林函数法
4.2.1 格林函数的定义 设在 内有 u 0, v 0; u, v 在 上有一阶连续 1 v u 偏导数,则由格林第二公式有 0 (u n v n )dS (2) 4 将(1)和(2)两式加起来:
u(M 0 ) 1 4 1 1 u u (v ) (v ) dS (3) n rMM 0 rMM 0 n
4.1.4 调和函数的性质
u u 0, | f . n
u n dS f dS 0.
6
下午10时1分
数学物理方程与特殊函数
第4章格林函数法
性质2 (平均值定理) 设函数 u(M ) 在区域 内调和, M 0 是 内任意一点,若 a 是以 M 0 为中心,a为半径 的球面,此球完全落在区域 的内部,则有 1 u(M 0 ) udS(调和函数的球面平均值公式) 2 a 4a 证明: 由调和函数的积分表示:
数学物理方程第四章 格林函数法
为边界的有界连通区域,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
数学物理方程第四章 格林函数法简化
x2 y2 z2 r
u 1 称作拉普拉斯方程的对称解. r
(2)u
2u x2
2u y2
0
(x, y) D
2u 1 u 1 2u
r 2 r r r 2 2 0
x r cos
y
r
sin
可以验证u(x, y) ln
1
ln 1 满足方程(2)
Dirichlet问题u
2u x 2
2u y 2
2u z 2
0
(x, y, z)
R3
(1.1)
u(x, y, z) f (x, y, z),
(1.2)
第二边值问题(牛曼(Neumann)问题)
在某一个光滑的闭曲面 上给出一个连续函数 f (x, y, z),要寻找一个函数u(x, y, z)它在 的内部区域 中是调和的,在 上连续,且在 上的任意一点
(2)
t
2
a2
2w x2
(x ,t )
x
, t
0
(1.1)
w(0, x) (x), w (x) x (1.2)
t t0
t
u(x,t) d f (,)w(x,t,,)d
0
w( x, , t, )称作初值问题(1)格林函数 .
§1 拉普拉斯方程边值问题的提法
• 静态薄膜的横向位移----二维拉普拉斯方程(也称调和方程)
2u x 2
2u y 2
0
(x, y) R2
第四章_拉普拉斯方程的格林函数法
注:对于外问题来说,求解通常都是在无界区域上,
这时需不需要对解加些限制条件呢?看下面一例子。
易知
u 0, r 1,
u 1 r 1
其中r x2 y2 z2
u 1,
u 1/ r
都是上述定解问题的解,即解不唯一.为了保证解的唯一性,
n
的值,所以要想求得狄氏问题的解就要想法消去积分公式中的
u 。故而我们需要引入格林函数。 n
在第二格林公式 (u2v
v2u)dV
(u
v n
v
u )dS, n
中取u, v C1(),并且都是内的调和函数.则
(u
v n
v
u )dS n
P Q R
(
x
y
z
) dV
Pdydz
Qdzdx Rdxdy
其中取外侧位正向.
由两类曲面积分之间的关系得高斯公式的另一种形式:
(
P x
Q y
R z
)dV
(P cos(n, x) Q cos(n, y) R cos(n, z))dS.
Ka表示以M0 (x0, y0, z0 )为中心,以a为半径且完全落在内部的球面,
则成立下面平均值公式
1
u(M0 ) 4 a2 Ka udS
证明: 将调和函数的积分公式应用到Ka可得
u(M 0 )
1
4
(u(M )
n
(1) r
1 r
数学物理方程课件 第四章拉普拉斯方程的格林函数法
S
u
(
v n
1
4
1 (
n rMM0
)
(
1
4
1 rMM0
v)
u n
dS
v为调和函数,且满足v
1
4
1 rMM 0
1 1
u(M0 )
u ( (
v) dS
S n 4 rMM0
11
G(M0 ) 4 rMM0 v
2u(M ) 0, 内 u | f (M )
G
u(M0 ) S u n dS
数学物理方程与特殊函数
第4章格林函数法
3 区域的格林函数和狄氏问题的解
v为调和函数,且满足v 电象法求格林函数
1
4
1 rMM 0
11
G(M0 ) 4 rMM0 v
在区域外找出区域内一点关于边界的象点,在这两个点放
置适当的电荷,这两个电荷产生的电位在曲面边界上相互抵消。
这两个电荷在区域中形成的电位就是所要求的格林函数。
1
4 rMM0 F (M0 )dV0
(t)
线性系统
h(t)
f (t)
线性系统
f (t) h(t)
G(M , M0 ) F(M0)
1
4 rMM0
F(M0 )dM0
1
4 rM rM0 F (M0 )dM0
数学物理方程与特殊函数
第4章格林函数法
2u(M
)
r
r0
,
内
u | 0
G(M , M0)
1
R1
G(M , M0 ) 4 rMM0
4 rOM0
rMM1
M R
o M0
M1
4格林函数法
那么,如何求解某边界条件下的泊松方程呢?
格林函数法求解拉普拉斯方程的边值问题
边值问题与无界空间的问题不同,要受到边界的影响, 边值问题的解 u(r) 与格林函数 G 的关系就更复杂了,需要 用到格林第二公式。 格林第二公式 设函数 u(r), v(r) 在区域 Ω 直到其边界 Γ 上具有连续一阶 微商,而在 Ω 中有连续二阶微商,有格林第二公式
若知道一个点源在一定的边界条件和(或)初值 条件下所产生的场(称为格林函数),就可用叠加 的方法计算出任意源产生的场。这就是格林函数法 的基本思想。故也称点源法。
δ – 函数(狄拉克函数)
用来描述物理中集中分布的量,如点电荷、点热源、质点、单 位脉冲等,这是通常的函数概念不能描述的。
0 x x0 ( x x0 ) x x0 0 x0 ( a , b ) a ( x x0 )dx 1 x0 (a , b)
r r0 o
2
R (x0, y0, z0)
G G x x 0 , x R R
2G 2G x x0 G 1 ( x x0 )2 2 2 3 x R R R R R
2G 2G 同样计算可 2 和 ,于是 2 y z 2G R 2 G 3 R 2 2G 2 G 2G 3 2 2 R R R 2 R R 0 R R R
于是有积分形式的解
拉普拉斯方程的格林函数
根据格林函数法的基本思想,先求解方程
(x, y, z) r r0 o
G ( r , r0 ) ( r r0 )
2
《数学物理方程》第4章 Laplace方程的Green函数法
w 有 w wdv w ds w wdv n 即 0 0 w wdv 向量 w 0
2
w w w 0 u u w C ( 常数) x y z 1 2
狄内问题,由 w 0 得 C 0 , 故 w 0 即 u u . 1 2
u( M 0 )
2 2
z 0
2
1
1
f( x ,y )z 0
2 2 23 2
[( x x ) ( y y ) z ] 0 0 0
dxd
u 0 x , y ,z 0 2 2 例1求 的 u ( 0 , 0 ,a ) 的值 1 x y 1 u ( x ,y , 0 )f 2 2 (常数 a 0) 0 x y 1
二 狄内问题Green函数法的步骤
1 . 任意取给定 M . v ( M , M ) 0 0 0 4 r 2 v 0 M 2 . 由 求出 v ( M ,M ) 0 v v ( M ) 0 则 G ( M , M ) v ( M , M ) v ( M , M ) 0 0 0 0
2
2
v 即 ( u v u v ) dv u dS n v 2 u vdv u dS u v dv n u 2 v udv v dS u v dv 再取 V v u n u v 2 2 两式相减 ( v u u v ) dv ( v u ) dS n n
§3 Green函数法
一 问题的提出
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
x
2
y
2
z
2
0
2 u u 1 u [ r sin ) (sin ) ( )] 0 (2.4) 2 r sin r r sin
(2.4)称作球坐标下的拉普拉斯方程。
如果u( x, y, z )具有球对称,即u与,无关,只与r有关,
拉普拉斯方程的外问题是在无限区域上给出的, 定解问题的解在无穷远处是否应该加以限制?
事实上如果不加以限制,外问题的解不一定是唯一的。
例如:考察以原点为心的单位球面 作为边界曲面的狄利 克雷外问题,并给出边界条件: u ( x, y , z ) 1
上述问题可以表示为
可以证明
2u 2u 2u 2 2 2 u 0 ( x , y , z ) x y z 1 2 2 2 x y z 2 2 2 u( x , y, z ) 1 x y z 1
2.2 拉普拉斯方程的对称解
首先介绍拉普拉斯方程的球对称解,前面我们知道
u( x , y , z ) 1 x2 y2 z2 1 r
x r sin cos 满足拉普拉斯方程。做变换 x r sin sin z r cos 2u 2u 2u
x
x r
u xx ( x, y, z ) 2 ( 3 ) r r x
u xx u yy u zz r 3 3x 2 r r
6
rx
r 3 3xr 2 rx r
6
r 3 3x 2 r r6
r 3 3 y 2r r
6
r 3 3z 2 r r
6
0
在二维平面的圆域中, 2u 2u 2 u 1 u 1 2 u 2 0 2 2 0 (2.7) 2 2 x y r r r r (2.7)为极坐标下的二维拉普拉斯方程。 如果 u ( x, y ) 关于原点对称性, 即 u (r , ) 不依赖 只与 r 有关,(2.7)变成
u f ( x, y , z ) n (1.3)
以上两个边值问题都是在边界上给定某些条件, 在区域 内部求拉普拉斯方程的解,称这样的问题为内问题。
在 应 用 中 我 们 还 会 遇 到 Dirichlet 问 题 、 Neumann 问题的另外一种提法,即在有限区域 外 部求函数 u ( x, y, z ) 使其在 的边界 满足条件 u u ( x, y , z ) f ( x, y , z ) 或 f ( x, y , z ) n 这样的问题相应地称作外问题。
1 1 u 1 u(M 0 ) u ( )]ds [ 4 r n n r (2.9)
证明:利用公式
(uv vu )dV (u
v u v )ds n n
Hale Waihona Puke v v v 令P u ,Q u , R u , x y z P u v 2 v Q u v 2 v R u v 2v u 2 , u 2 , u 2 x x x x y y y y z z z z
P Q R 2 v 2 v 2 v u v u v u v u( 2 2 2 ) x y z x x y y z z x y z
第二边值问题(牛曼(Neumann)问题)
在某一个光滑的闭曲面 上给出一个连续函数 f ( x, y, z ) , 要寻找一个函数 u ( x, y, z ) 它在 的内部区域 中是调和的,在 上连续,且在 上的任意一点 u 沿 的单位外法线方向 n 的方向导数 存在,并且等 n 于已知函数 f ( x, y, z ) 在该点的函数值,即
二维情形要求在无穷远处的极限有界,即
lim u( x, y)
r 0
§2 调和函数
2.1 格林公式
奥-高公式: 设 是以分片光滑的曲面 为边界的有界连通区域, p( x, y, z ) 、 Q( x, y, z ) R( x, y, z ) 是 上连续,而在内具有一阶连 续偏导数的任意函数,则有
P cos(n, x ) Q cos(n, y ) Rcos(n, z )
代入下式,
u( v v v v cos(n, x ) cos(n, y ) cos(n, z )) u x y z n
(
P Q R )dV ( P cos(n, x ) Q cos(n, y ) Rcos(n, z ))ds (2.1) x y z
(
P Q R )dV ( P cos(n, x) Q cos(n, y ) R cos(n, z )) ds x y z
其中 n 是 的单位外法线向量, ds 是 上面积元素。
格林第一公式: 设 u ( x, y, z ) 、 v( x, y, z ) 以及它们的所有一阶偏导数在闭区 域 上是连续的,他们在 内具有连续的 所有二阶偏导数。令
u u v u v u v ds ( )dV n x x y y z z
(2.2' )
(2.2)-(2.2’)可得格林第二公式:
(uv vu)dV (u
v u v )ds (2.3) n n
牛顿 - -莱布尼茨公式
(2.2)
(2.2' )
(vu)dV v
u, v互换
(uv )dV u
v u v u v u v ds ( )dV (2.2) n x x y y z z
(vu)dV v
u1 ( x, y, z ) 1
u2 ( x , y , z )
1 x2 y2 z2
都是解。
1 r
因此在无穷远点不加任何限制时,外问题的解就不一定唯一。 那么无穷远点处的条件应该如何加?一般地在三维情形要求
lim u( x, y, z ) 0
r
r
x2 y2 z2
第四章 格林函数法
主要内容 第一边值问题(狄利克雷(Dirichlet)问 题)第二边值问题(牛曼(Neumann)问题)
格林第一(二)公式
调和函数的基本性质 *格林函数的定义 及特殊区域上格林函 数的求法
§1 拉普拉斯方程边值问题的提法
静态薄膜的横向位移----二维拉普拉斯方程(也称调和方程)
则称问题( 1.1 ) ~ ( 1.2 )为拉普拉斯第一边值问题或狄利克雷 (Dirichlet)问题, u ( x, y, z ) 为此问题的解。
2u 2u 2u 3 u 0 ( x , y , z ) R (1.1) 2 2 2 x y z u( x , y , z ) f ( x , y , z ), (1.2)
u 0.描述的是稳态时(静态 )物理量,与时间无关 ,
故不提初始条件 !只给出边界条件就可以 . 下面看边界条件的提法 .
(1) 第一边值问题(狄利克雷(Dirichlet)问题) 3 设方程(1.1)的空间变量 ( x, y, z ) , 为 R 的开区域。如果 u ( x, y, z ) 满足方程 (1.1) , 且在 边界 上直接给定了 u ( x, y, z ) 的具体函数形式 f ( x, y, z ) ,即 u ( x, y , z ) f ( x, y , z ) (1.2)
则有格林第一公式:
(uv )dV u
v u v u v u v ds ( )dV n x x y y z z
u u v u v u v ds ( )dV n x x y y z z
(1) 第一边值问题(狄利克雷(Dirichlet)
2u 2u 2u u 2 2 2 0 ( x, y, z ) R3 (1.1) x y z
u( x, y, z) f ( x, y, z)
(2)第二边值问题(牛曼(Neumann)问题)
(1.2)
2u 2u 2u u 2 2 2 0 ( x, y, z ) R3 (1.1) x y z
2u 2u 2 0 ( x , y ) R x 2 y 2 (1.1)
u 0.
2u 2u 2u u 2 2 2 0 ( x, y, z ) R 3 (1.1') x y z
2 2 2 2 2 2 称作拉普拉斯算子 . x y z
3
2.3 调和函数的基本性质
在开区域 内,称具有二阶连续偏导数并且满足 拉普拉斯方程的连续函数为调和函数。又称此函数 在 内是调和的。
性质 1 (积分表达式)设 是以分片光滑的曲面 u ( x, y, z ) 在 上有连续 为边界的有界连通区域, 的一阶偏导数, 在 内调和, 定点 M 0 ( x0 , y0 , z0 ) , r 为定点 M 0 到变点 M ( x, y, z ) 距离: 则有
2 u 1 u 0 (2.7') 2 r r r
显然它的解为(如取 c1 1, c2 0 )
1 u c1 ln r c2 u ln r
(2.8)
1 1 综上所述, 除 r 0外, , ln 分别是 (2.5) 、 '7 . 2 ( r r
2 2 2
)的
u( x, y, z ) f ( x, y, z ), n
2u 2u 2u 3 u 2 2 2 0 ( x , y , z ) R (1.1) x y z u( x , y , z ) f ( x , y , z ), (1.3) n