第二章有限差分基础.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 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 ) 时间向前差分,空间中心差分得

文档评论(0)

汪汪队 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档