| 查看: 822 | 回复: 8 | ||||
我要飞铁虫 (正式写手)
|
[求助]
数值解高手 已有2人参与
|
|
如何分析下面三层差分格式的收敛性？ |
回复此楼» 猜你喜欢
求调剂一志愿武汉理工大学材料工程（085601）
已经有7人回复
-
一志愿C9材料与化工专业总分300求调剂
已经有4人回复
-
321求调剂
已经有4人回复
-
307求调剂
已经有4人回复
-
293求调剂
已经有6人回复
-
085602 289分求调剂
已经有3人回复
-
一志愿北京化工大学材料与化工 264分各科过A区国家线
已经有3人回复
-
277分求调剂，跨调材料
已经有3人回复
-
材料292调剂
已经有6人回复
-
材料与化工考研调剂
已经有6人回复
zaq123321
专家顾问 (著名写手)
-
专家经验: +342 - 数学EPI: 6
- 应助: 298 (大学生)
- 贵宾: 0.247
- 金币: 11336.3
- 红花: 29
- 帖子: 1221
- 在线: 538.8小时
- 虫号: 405284
- 注册: 2007-06-17
- 性别: MM
- 专业: 生物大分子结构与功能
- 管辖: 数学
【答案】应助回帖
感谢参与，应助指数 +1
| Suppose this is the Richardson scheme for diffusion equation and it's a consistent scheme. According to Lax equivalence theorem, the stability of this scheme is equivalent to the convergence. It has been proved that the scheme is not stable, thus not convergent. See the example 3.4 in the second edition of Numerical solution of PDE by Lu Jinfu and Guan Zhi. |
2楼2014-11-28 21:36:23
我要飞
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1515.3
- 散金: 144
- 红花: 5
- 帖子: 458
- 在线: 591.2小时
- 虫号: 671362
- 注册: 2008-12-09
- 专业: 计算数学与科学工程计算
3楼2014-11-29 08:38:23
感谢参与，应助指数 +1
feixiaolin: 屏蔽内容, 笔误。修正后见5楼 2014-11-29 13:58:41
feixiaolin: 应助指数-1 2014-11-29 13:59:19
feixiaolin: 屏蔽内容, 笔误。修正后见5楼 2014-11-29 13:58:41
feixiaolin: 应助指数-1 2014-11-29 13:59:19
|
本帖内容被屏蔽 |
4楼2014-11-29 13:40:46
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
5楼2014-11-29 13:54:03
我要飞
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1515.3
- 散金: 144
- 红花: 5
- 帖子: 458
- 在线: 591.2小时
- 虫号: 671362
- 注册: 2008-12-09
- 专业: 计算数学与科学工程计算
6楼2014-11-29 17:07:01
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
7楼2014-11-29 19:40:54
我要飞
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1515.3
- 散金: 144
- 红花: 5
- 帖子: 458
- 在线: 591.2小时
- 虫号: 671362
- 注册: 2008-12-09
- 专业: 计算数学与科学工程计算
8楼2014-11-30 08:22:22
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
|
好吧，你想从收敛的定义来而不是Lax equivalence theorem 判断此格式不收敛。 不去咬文嚼字，收敛性的定义大致就是离散解uh逼近精确解u，当dt和tdx趋于零时。 我们知道原始解析方程的精确解一定是有界的，因为diffusion operator F 的所有的特征值都是纯实数小于零。事实上，对于一类边界条件，u最终会变成一条直线，这里就不详细证明了。可是，leap-frog scheme 的离散特征值lambda显然有大于零的实部，因为我们在楼上实际上是证明了 r = exp (lambda dt) 的绝对值大于零，也就是lambda有大于零的实部。那么就是说，uh(t) 会指数增长; 从而不会收敛到一个有界的精确解u(t)。这就证明了uh(t)是不收敛的。线性数值稳定性分析无非就是看离散谱点如何近似于精确谱点。 事实上，Lax equivalence theorem 有相当的建设性。因为，如果证明了不稳定性，也就是知道不稳定的数值特征值（很多线性稳定性分析都是直接找到增长的离散谱点）。增长率就决定了发散的速率。所以知道不稳定，就知道不收敛的方式，包括空间震荡（特征函数）的具体形式。不知道这是你不是你说的“收敛性的详细过程”。 |
9楼2014-11-30 09:37:12
