隐式差分方程课件
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2.3.3 加权六点隐式格式
前面,我们已经推导了热传导方程(2.26)的古典显示格式,古
典显示格式及Crank-Nicolson格式等。实际上,它们都可以作为本 节推导的加权六点隐式格式的特殊情形。
由 umn1 exp(k Dx2 )umn
得 exp(kDx2 )umn1 exp((1 )kDx2 )umn ,1 1
2.3 隐式差分格式
与显式差分格式不同,隐式差分格式中
包括了(n+1)时间层上二个或二个以上结
点处的未知值(例如
U
n 1 m 1
,U
n 1 m
,
U
n 1 m 1
),
使用隐式差分格式和使用显式差分格式求
解完全不同。相对而言,使用隐式差分格
式求解,每时间层包含有较多的计算工作
量。从后面对差分格式的稳定性分析可知,
2
)
am 2
代入式(2.48),则 因此得差分方程
1 2h2
2 x
(umn1
umn )
1
n 1
am 2 k
(umn1
umn
)
1 12 r
1 h2
2 x
[
1
n 1
am 2
(umn1
umn )] (h4
k2)
(2.49.1) 1
n 1
(U
n 1 m
U
n m
)
1 2
r
2 x
)U
n1 m
(1
1 2
r
2 x
)U
n m
(2.43)
这是一个隐式差分格式,称为Crank-Nicolson差分格式,
截断误差阶为 (k 2 h2 ) ,也可写为
(2.44) (1
r
)U
n1 m
1 2
r
(U
n1 m1
U
n1 m1
)
(1
r
)U
n m
1 2
r
(1
)(U
n m1
U
n m1
)
r
(U
n1 m1
U
n1 m1
)
[1
2r(1
)]U
n m
,0
1(2.46)
这是一个六点差分格式(如图2.7所示),称为加权六点差分格式。
m-1,n+1
m,n+1
m+1,n+1
图2.7:
m-1,n
m,n
m+1,n
1
显然,当 0时,加权六点格式为古典显示格式;
0.956 821 703 419 0.000 079
0.915 507 772 134 0.000 151
0.643 146 895 793 0.000 531
0.413 637 929 568 0.000 683
0.171 096 336 778 0.000 564
0.629 273 956 459 0.000 194
1 2
2 x
[(1
1
n 1
)U
n 1 m
]
1 2
2 x
[(1
1
n 1
)U
n m
]
ram 2
6ram 2
6ram 2
格式(2.49.1)具有截断误差阶 (h4 k 2 ) ,可写
成更方便的形式
[1
1 12
a (a n1 2
2
n 1 2
m xm
)1
1 2
n
ram
1 2
k Dx2 )umn1
(1
1 2
k Dx2 )umn
由式(2.19.3),可令 Dx2u
n m
1 h2
2 x
(1
1 12
2 x
)umn
则可得
Dx2
1 h2
2 x
(1
1 12
2 x
)
1
代入上式,则有如下差分格式:
[1
1 2
(r
1 6
)
2 x
]U
n1 m
[1
当
1 2
时,加权六点格式为Crank-Nicolson隐式格式;
当 1 时,加权六点格式为古典隐式格式。
加权六点格式亦可直接由差商代替导数得到
U n1 m
U
n m
k
1 h2
U 2 n1
xm
(1
)
1 h2
U 2 n
xm
2.3.4 系数依赖于x,t的一维热传导方程的一 个隐式格式的推导
考虑方程 的差分逼近。
u t
a(x, t)
2u x 2
(2.47)
由其已Ta知ylor展开x式 2可sin得h(
1 2
hDx
)
据此,可得
1 2h2
2 x
(umn1
umn
)
(
2 x
2u x2 )
h2
n 1 2
m
Dx2
h4 12
Dx2
1 12
h
2
(
4u x 4
r
(U
n m1
U
n m1
)
由于格式(2.44)中包括六个结点,故也可 称为六点格式(如图2.6所示)。
图2.6
m-1,n+1
m,n+1
m+1,n+1
m-1,n
m,n
m+1,n
也可将(
u t
)
n m
1 2
u n1 m
umn
k
(
2u t 2
)
n m
1 2
1
[
u n1 m1
2
SCN
x
,ED
2
Sr SD 分别 ,时间层n
它们的值由表2.2给出。
tn t1 t2 t4 t8 t16 t80 t160 t320 t640 t800
表2.2
Sr
0.994 497 915 630
ECN
0.000011
0.489 026 104 192 0.000 022
0.978 172 634 773 0.000 040
0t T
应用(1) Crank-Nicolson差分格式,(2)
Douglas差分格式解上述问题。对每一种情况,
令 h 20, r 1 20 (r的这个值对Douglas格式
有最小的截断误差),由初值条件和边值条件通
过上述二个格式的每一个逐层求出U
n m
的值。一般
而言,当由第n层去求第(n+1)层的解时,二个
1 2
k 2 Dx4
)u
n 1 m
u
n m
式中左边如果仅保留二阶导数项,且以
格式
(1
k h2
2 x
)U
n1 m
U
n m
1 h2
2 x
替代 Dx2
,则得差分
或者
rU
n1 m1
(1
2r
)U
n1 m
rU
n1 m1
Uຫໍສະໝຸດ Baidu
n m
(2.41)
格式用图2.5表示,其截断误差阶为 (k2 h2) ,与古典差分格式相同。
2umn1 h2
u n1 m1
un m1
2umn h2
un m1
]
代入微分方程(2.26),得到Crank-Nicolson格式。
基于如同Crank-Nicolson格式一样的六个网
格结点可获得另一精度较高的差分格式,如在
前式(2.42)中仅保留直到 Dx2的项,即有
(1
即两边去[1掉高kD于x2 二12阶 2导k 2数Dx4的项,]umn且1 用 [1h12(x21代替)k
Dx2
Dx2
1 2
(1
)2
k
2 Dx4
]umn
,则得差分格式
(1r
2 x
)U
n1 m
[1
(1
)r
2 x
]U
n m
或者 (1
2r
)U
n1 m
隐式格式的优点在于,其稳定性要求对步
长比的限制大为放宽,而这正是我们所期
望的。
2.3.1 古典隐式格式
现在对热传导方程
u t
2u x 2
推导其最简单的隐式差分逼近——古典隐式格式。由umn1 exp( kDx2)umn
故
exp(
k
Dx2
)u
n 1 m
u
n m
(1
k Dx2
格式的每一个都需解一线性代数方程组,其系数
是三对角阵,可用追赶法求解(见2.4)。已知上
述定解问题的理论解,记为 Sr , 有 u et sin x
记SCN , SD 分别为用高速数字计算机解出的Crank-
Nicolson格式的解,而 ECN Sr 表示它们对精确解的误差,在
上, tn nk 。
0.012 108 818 740 0.000 100
ED
0.000 000 000 026 -0.000 000 000 051 -0.000 000 000 101 -0.000 000 000 198 -0.000 000 000 379 -0.000 000 000 331 -0.000 000 001 712 -0.000 000 001 417 -0.000 000 000 485 -0.000 000 000 257
图2.5:
m-1,n+1
m,n+1
m+1,n+1
m,n
为了求得第(n+1)时间层上的
U
的值, n1
m
必须通过解线性代数方程组。这是一个隐
式差分格式,必须联合其初边值条件求解。
格式(2.41)通常称为古典隐式格式。
我们也可以通过直接用差分算子代替 Dx, Dx2
的方法,即
( u t
) n 1 m
n
)m
1
2
(h6 )
(h4
k
2
)
( 1 a
u
)n
1 2
t m
h2 12
2 x 2
(1 a
u t
)
n m
1 2
(h4
k 2)
(2.48)
令
2 x2
(1 a
u t
n
)m
1 2
1 h2
x2[
1
n 1
1 k
(umn1
umn
)
(h2
k
1
2
2 x
]U
n1 m
[1
1 12
a (a n1 2
2
n 1 2
m xm
)1
1 2
n
ram
1
2
2 x
]U
n m
(2.49.2)
这是一个隐式差分格式(如图2.8所示)。
m-1,n+1 m-1,n
m,n+1
m+1,n+1
m,n
m+1,n
图2.8
感谢您的下载让小编的努力能帮助到您, 最后一页是小编对你的谢谢哦,提醒一下, 下载好了几个全部自己看一遍,把用不上 的删除哦!包括最后一页!
u n1 m
umn
k
(
2u x 2
) n 1 m
u n1 m1
2umn1 h2
u n1 m1
代入微分方程,得到格式(2.41)。
2.3.2 Crank-Nicolson隐式格式
Crank-Nicolson隐式差分格式是解热传导方程(2.26)
的常用的差分格式,为了推导它,由式(2.24),有
exp(
1 2
k L)umn1
exp(1 2
k L)umn
由 L Dx2
得 [1
1 2
k Dx2
1 2
(1 2
kDx2 )2
]umn1
[1
1 2
k Dx2
1 2
(1 2
kDx2 )2
]umn
(2.42)
两边仅保留前二项,用
1 h2
代替 2
x
Dx2
,则得差分格式
(1
1 2
(r
1 6
)
2 x
]U
n m
(2.45)
它称为Douglas差分格式,具有截断误差阶
(k 2 h4 )。
例2.1 解初边值问题
{u 2u t x 2 u t 0 s in x u(0, t ) u( , t ) 0
0 x ;0 t T 0 x