偏微分方程离散差分格式差分方法等共57页文档

满足稳定性要求的
amplificat ion factor G
A n1
G An 1
10
3.1.4 差分方法的理论基础（续）
• Fourier (Von Neumann) 稳定性分（续）
G1eikx 1(coksxisinkx) 1(1coskx)isinkx
G2 1(1coskx)2 2sin2 kx14(1)sin2kx,
2
G1 if 1
• 1 称为CFL条件 (Courant, Friedrichs, Levy)
11
3.1.5 守恒型差分格式
• 流体力学方程组描述物理量的守恒性；守恒律组：
• 定义
u d f 0
t i1 xi
对于一维单个守恒律：
u f (u) 0 t x 其差分格式如果具有如
下形式
u n1 j
n 1
xJ 1/ 2 u ( x,0)dx
x J 1 / 2
t 0
n
1
u
(
x
J
1
,
t
)
dt
( 1 ) p 2 p 1 k 2 p 1 p0
格式稳定的充分必要条
件是
( 1 ) p 2 p k 2 p 0 , k p 1
偶次项系数 满足 : ( 1 ) p 2 p 0
对于（ 2）：
1
c，
2
0, 3
1 6
c(c 2 t2
x2)
4
1 c2t 8
(3c 2 t2 x 2 )
t
(e t
1)4
t4
4u t4
可以将
t 表示成
t
(e t
1 ) l 的级数
t
t
t
l1
b l
e
t t
1 l , b 1
1, b 2
1 2
,b3
1 3
,b3
3 8
最后得到
t
t
b
l
e
t
t
1 l
l1
l
bl
1 2
e
x
x
x
e x
1 2
2
e
x
x
2
x
e x
即有
3
3.1.1 模型方程的差分逼近
4
3.1.2 差分格式的构造
5
3.1.3 差分方程的修正方程
• 差分方程所精确逼近的微分方程称为修正方程 • 对于时间发展方程，利用展开的方程逐步消去带时间的高阶导数，只留空间导数。 • Warming-Hyett方法：
u c u 0
(1)
t x
u
n 1 j
u
t
k
k 1
ku xk
p0
2 p1
2 p1u x 2 p1
2p
p 1
2pu x2p
(5)
7
3.1.3 差分方程的修正方程(续)
u
t
k
k 1
ku xk
基本解为
p0
2 p1
2 p1u x 2 p1
2p
p 1
2pu x2p
e e ( i ) t ikx
( 1 ) p 2 p k 2 p p 1
2
x
e x
u
(4)
记算子
t
(e t
1)
t
u t
1 t2 2!
2u t2
1 t3 3!
3u t3
则
t
(e t
1)2
t2
2u t2
1 2
1 2
t 3
3u t3
2
1 6
1 2
1 2
t 4
4u t4
t
(e t
1)3
t3
3u t3
1
1 2
t 4
4u t4
偏微分方程离散差分格式差分方法等
11、用道德的示范来造就一个人，显然比用法律来约束他更有价值。—— 希腊
12、法律是无私的，对谁都一视同仁。在每件事上，她都不徇私情。—— 托马斯
13、公正的法律限制不了好的自由，因为好人不会去做法律不允许的事 情。——弗劳德
14、法律是为了保护无辜而制定的。——爱略特 15、像房子一样，法律和法律都是相互依存的。——伯克
（三）偏微分方程的数值离散方法
• 3.1 有限差分法 • 3.2 有限体积法 • (有限元，谱方法，谱元，无网格，有限
解析，边界元，特征线）
2
3.1 有限差分法
• 3.1.1 模型方程的差分逼近 • 3.1.2 差分格式的构造 • 3.1.3 差分方程的修正方程 • 3.1.4 差分方法的理论基础 • 3.1.5 守恒型差分格式 • 3.1.6 偏微分方程的全离散方法
u
n j
1 2
u j1 u j1
1 2
2
u j1
2u j
u j1
(2)
Taylor 展开
u n1 j
u
n j
t
u t
1 t2 2!
2u t 2
1 t3 3!
3u t 3
t
(e t 1)u
u
j 1
u
j
x
u x
1 2!
x2
2u x 2
1 3!
x3
3u x3
x
(e x 1)u
x
u j1 (e x 1)u ( 2 ) 等价于：
u 1 t t 2
2u t 2
1 t2 6
3u t 3
c
u x
1 x2 6
3u x3
c
2
t
1 2
2u x 2
1 24
4u x 4
(3)
• 差分方程(2)写成算子的形式：
6
3.1.3 差分方程的修正方程 (续)
t
(e t
1)u
1 2
x
(e x
x
e x )u
1 2
2
e
x
x
• 守恒性质：
守恒型差分格式对 j求和 :
j J
u
n j
1
x
jJ
jJ
u
n j
x
jJ
~f n J
1
t
2
~f n J
1
t
2
再对 n求和 :
jJ
u
n j
1
x
jJ
jJ
u
0 j
x
jJ
N k 0
~ f
k
J
1
t
2
N k 0
~ f
k
J
1
t
2
可以看成是积分
u ( x, t )dx xJ1/ 2
x J 1 wenku.baidu.com 2
符合 War min g Hyett 稳定性判别条件
. why CFL 1 for scheme (2) ?
8
3.1.4 差分方法的理论基础
• 相容性，稳定性，收敛性 • 等价性定理 • Fourier稳定性分析
9
3.1.4 差分方法的理论基础（续）
• Fourier (Von Neumann) 稳定性分析
uin 1 tuinc 1 x(uin 1uin)0 , c0(1)
设 c t x
误差的基本解
代入 (1 ) :
u
n i
A n e ikx i
u n1 i
A e n 1 ikx i
u A e n
n ikx i1
i1
u
n 1
u
n i
(
u
n i
u
n i
1
)
A e n 1 ikx i A n e ikx i ( A en ikx i A e ) n ikx i1
u n1 j
t x
~f
n j
1
2
~f
n j
1
2
则为守恒型差分格式。
其中
~f
n j
1
称为数值通量，它是
2 l 个变量的多变量函数：
2
~f
n j
1
~f
(u
n jl
1
,
u
n j
l
2,
,
u
n jl
),
2
~f 满足相容性条件 :
~f (u , u , u ) f (u )
12
3.1.5 守恒型差分格式（续）