原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
北航计算流体力学第9课
对于迁移方程: (5-1)
1) Euler显式格式
(5-2)
其等价微分方程为: (设)
2) MacCormack格式
(5-6)
其等价微分方程为:
3) Euler隐式格式
(5-7)
其等价微分方程为:
4) Crank-Nicolson格式
(5-8)
其等价微分方程为:
一. 对于迁移方程
(5-1)
设通解为:
(5-10)
式中 , k——波数， ——频率
将通解(5-10)代入方程(5-1),即求得
,
代入(5-1), 整理得:
即:
代入通解式(5-10),得
整理得 (5-11)
式(5-11)表示各种不同波数(k不同)的波都以同样的波速向下游迁移,如图1所示.
所以把方程(5-1)称为迁移方程.
二. 对于
(5-12)
它的通解仍为
(5-10)
求出
将它们代入方程(5-12), 得
整理得
将上式代入通解(5-10), 得到
(5-13)
与迁移方程的通解
相比较,可见方程(5-12)的解表示不同波数的波以不同的波速向下游迁移, 波数越大, 迁移速度就越小, 如图2所示
这种现象称为频散.
三. 对于
(5-14)
它的通解仍为
(5-10)
求出 , 和代入方程(5-14), 整理得
代入通解(5-10), 整理得
(5-15)
上式表明, 不同波数波都以同一波速a向下游迁移, 但在迁移的过程中, 由于, 初始扰动会越来越小. 这种现象称为耗散. 如图3所示.
需要指出的是对于方程 , 其通解为
当时, , 即, 这种现象称为负耗散.
由Taylor公式
（1）
（2）
式（1）（2）得
中心差: (3)
由式(1)得
前差: (4)
由式(2)得:
后差: (5)
由以上式子可见:
1) 中心差分的截断误差项中只包含奇次导数项, 所以由中心差分构成的差分格式具有频散特性;
2) 前差和后差的截断误差项中包含有所有导数项, 且第一项为, 因此由前差和后差构成的格式具有耗散特性,且
a) 当扰动从上游传向下游时, 即, 后差格式为正耗散;
b) 当扰动从下游传向上游时, 即, 前差格式为正耗散.
下面分析几个常见格式频散、耗散特性(以为例)
1. Euler显式格式
等价微分方程为:
显然, 它具有耗散特性(当时)
2. MacCormack格式
等价微分方程为:
显然, 它是频散格式.
3. Crank-Nicolson格式
等价微分方程为:
它是频散格式.
对于Euler中差格式
其等价微分方程为:
显然, 这是负耗散格式，解发散.
作业5
5.1 试写出下列微分方程的通解并分析其频散, 耗散特性
1.
2.
3.
4.
5.2 对于迁移方程() 试推导跳点格式的等价微分方程, 并指出它属于哪种格式.
图1. 的解
图2. 的解
图4. 的解
图3. 的解
文档评论（0）