- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PAGE 1
PAGE 13
有限差分基础(finite difference method,FDM)
偏微分方程的一般形式
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 1 )
网格划分
一般有限差分采用结构化网格划分。即节点对应于当地坐标系统的原点。它的轴同网格线一致。即两个同一族的网格线不相交,且没对网格线对应不同的族。每一个节点可用唯一的一个坐标表示,如(?1, ?2)。网格线能用?1=const, ?2=const 表示。
i-1i
i-1
i
i+1
N
1
N
Nj
j+1
j-1
j
1
1
i-1
i
i+1
(i,j)
2D
有限差分法就是要将节点上的偏微分方程用相邻点上的值表示,变成线性代数方程式。
为流体力学的微分方程的数值求解方法之代表。必要条件:
连续领域内的分配有限的网格
领域内的函数分布可用网格点上的值代表
PWEEENNSSSWW
P
W
E
EE
NN
S
SS
WW
N
P
P
E
W
B
T
S
N
P
P
E
W
S
N
5点计算分子 15点计算分子 7点计算分子
这些节点又称为计算分子。方程的个数应与未知数相同,即每个节点有一个方程式。
Taylor展开
例如:一维时间变量?的理论解为?(t,x),它在离散点上的值为投影(projection)的近似值为:,
n: 时间的step数 i:空间的step数
为了求得此近似解,需对微分方程进行差分近似。利用Taylor展开可得到几个差分表示形式,仅考虑空间依存问题:
在?x很小时,位置j?x内的物理量?用?j来表示,则位置(j+1)?x上的值?j+1表示为:
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 2 )
(j-1)?x上的值?j+1表示为:
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 3 )
基本差分格式
一阶导数(first derivative)的近似
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 4 )
向前差分(forward difference,FDS)
利用( 1 ) 式,可得到1阶微分的向前差分形式:
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 5 )
向后差分(backward difference,BDS)
利用( 1 ) 式,可得到1阶微分的向后差分形式
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 6 )
中心差分(central difference,CDS)
(1)-(2) 得1阶微分的2 次精度中心差分法:
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 7 )
上风法、迎风法(upwind difference, UDS)
与速度有关的微分
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 8 )
二阶导数的近似
中心差分(central difference,CDS)
利用(j?1/2)?x的Tayor展开,可得过且1阶微分的2次精度的向前向后差分形式:
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 9 )
将上二式相减,得2阶微分的差分方程式中心差分:(?x相当)
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 10 )
其它还有各种形式。但最常用的是中心差分。对于扩散项常采用
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 11 )
向后差分
向后差分
向前差分
中心 差分
精确值
i+1
i
i-1
?xi
?xi+1
?
x
数值误差
截断误差(truncation error)
例如:热传递方程
( STYLEREF 1 \s 2 SEQ ( \* ARABIC \s 1 12 )
时间向前差分,空间中心差分得
您可能关注的文档
- 企业的战略规划和营销管理超全.ppt
- 土木地质第1章3-.ppt
- 高一化学下学期第一次月考试题4.doc
- 滦县小学2018-2019学年一年级下学期数学3月月考试卷.doc
- tornado调试技巧课件.ppt
- 第一章中国传统社会.ppt
- 社旗县实验小学2018-2019学年一年级下学期数学3月月考试卷.doc
- 2018澳洲489偏远地区技术移民详解.docx
- 数据库新技术介绍2018.pptx
- 安徽省安庆市宿松县九姑中学2019-2020年高一历史下学期期末试题(含解析).doc
- 第九章放射防护法规与标准...ppt
- 高一物理机械能守恒定律2.ppt
- 医疗器械生产质量管理规范无菌医疗器械现场检查要点(37张).ppt
- 2.4.2平面向量数量积的坐标表示、模、夹角(比赛)1.ppt
- 宣化区第三小学2018-2019学年一年级下学期数学3月月考试卷.pdf
- 甘南县第三中学校2018-2019学年高三上学期11月月考数学试卷含答案.pdf
- 名山区小学2018-2019学年二年级下学期数学3月月考试卷.pdf
- 最新语文S版四年级语文下册26《鲸》教案公开课说课稿(教学设计)d.pdf
- 期末复习材料题.ppt
- 人教版高中化学必修一第3章第1节《金属的化学性质》ppt教学课件3.ppt
文档评论(0)