《抛物型方程的差分方法ppt课件.ppt》由会员分享，可在线阅读，更多相关《抛物型方程的差分方法ppt课件.ppt（39页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、非线性动力学数值方法第二章第二章 抛物型方程的差分方法抛物型方程的差分方法2.1 定解问题的离散定解问题的离散变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分2.1 定解问题的离散定解问题的离散 一维热传导方程为： 或 对这样一个问题的求解，分为以下三个步骤来离散。0(, 0 )()tx xL uuuuxx01,0 xtT 220uuLutx（1）变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分在x-t平面上，取 和 分
2、别为函数 的自变量x和t的改变量，由 （ j=0,1,N, h= , n=0,1,M, ）两组平行线构成的矩形网格覆盖x-t平面。h为空间步长， 为时间步长。txnj（j-1,n）（j ,n）（j+1,n）（j,n+1）（j,n-1）hx t ( , )u x t1NTM2.1.1定解区域的离散定解区域的离散X方向节点编号T方向节点编号变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分 定义：定义： 网格节点上的值： 半网格节点上的值： 网格节点上的函数值 简记为u (j, n) 。,0,1,jxjhjN1/
3、21,0,1,12ntnnM1/21,0,1,12jxjhjN,0,1,ntnnM(,)jnu xt变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分 在有限差分离散化时应该注意以下几点：在有限差分离散化时应该注意以下几点：根据问题求解的需要，在x，t方向上离散网格时 和 可以是等分，也可以是不等分，既可以按一定规律来离散，也可以对网格进行局部加密。！对于双曲型和抛物型等发展方程在有限差分离散化时，网格的 和 不能随意选取， 需要满足一定的条件，如稳定性的CFL条件等。！为了保证边界上的计算精度，在网格边界外
4、可设置若干虚拟网络，以保证差分格式在边界处的计算精度和内点精度保持一致。！有限差分离散网格一般选取四边形网格，对于复杂物体外形问题，也可以选择三角形网格或其他形状网格。近年来，发展了一种无结构网格和无网格的有限差分算法，它们的计算网格就更为复杂。txtxxt变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分对于复杂外形飞行器流场的计算，一般需要通过坐标变换，可以把物理平面上的复杂的、非正交的网格转换成在计算平面上的简单、而正交的网格，这就是网格生成技术。特别要指出的是，网格生成技术在网格设计和编程中往往占有很
5、大的工作量，网格生成技术好坏直接影响到数值计算结果的精度，网格生成技术已成为计算流体力学中的一个重要分支。变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分节点(j, n+1)的函数值在(j, n)点作泰勒展开：2.1.2 控制方程的离散控制方程的离散2233232,( ,)( ,)( ,)( ,1)( ,)2!3!12tttj nu j nu j nu j nu j nu j ntttuuu（2）2,1( ,1)2tttj nu j nuuu（3）同理，对于节点(j, n-1)有：变电站电气主接线是指变电站
6、的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分 由（2）得： 其中 表示 一次和一次以 上的小量项22( , )( ,1)( , )( , )2u j nu j nu j nu j ntt( )O22( , )( , )( ,1)( , )2u j nu j nu j nu j ntt224224( ,1) 2 ( , )( ,1)2( , )( , )4!u j nu j nu j nu j nu j ntt （4）（5）（6）（4）（5）得：由（3）得：变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成
7、输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分224224(1, )2 ( , )(1, )2( , )( , )4!u jnu j nu jnhu j nu j nxhx（7）同理：变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分下面引入几个概念：下面引入几个概念：（1）向前差分）向前差分 （forward difference） （2）向后差分）向后差分 (backward difference) 1nnjjnhjuuuh1nnjjnjuuu1nnjjnhjuuuh1nnjjnjuuu(
8、 forward space difference)(forward time difference)(backward space difference)(backward time difference)变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分（3）一阶中心差分（）一阶中心差分（central difference）1122nnjjnhjuuuh1122nnjjnjuuu1122nnjjnhjuuuh1122nnjjnjuuu11222nnnjjjnhjuuuuh2222224nnnjjjnhj
9、uuuuh11222nnnjjjnjuuuu（）二阶中心差分（）二阶中心差分（central difference）变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分1）格式）格式 I显示格式显示格式 (FTCS格式格式)由（4）、（7）代入（1），有 式中： 称为截断误差（Truncation error），它不仅反映了差分算子对微分算子的逼近，也反映了差分解和方程解的误差。截断误差的阶数：就是截断误差中最低阶导数项中 或h的幂次数。（此截断误差，时间1阶，空间2阶）！22224224( , )( , )(
10、 ,1)( , )(1, )2 ( , )(1, )2( , )( , )24!0u j nLuu j ntxu j nu j nu jnu j nu jnuhj nu j nhtx ,nnnhjjjLuLuR2422422( , )( , )()4!2njhRu j nu j nOhxt（8）变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分用 表示u(j ,n）的近似值；用差商近似代替式（1）中的微商后，可得相应的差分方程（通过泰勒展开法最后推出的） (9)记号： 表示微分方程的解 在结点 处的准确值；
11、表示差分方程的解，它是 的近似值； 表示微分方程左端项在结点 处的准确值； 表示用准确值 构造差商； 表示用近似值 构造差商； 表示差商近似微商所产生的截断误差。nju nju, j n111I,220nnnnnjjjjjnhjuuuuuLuh nju( , )u x t, j nnjLunjunju ,nhjLu(,)jnu x t,nhjL unjR变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分注意：注意：由（10）可知，当第n层u已知时，可以直接求出第n+1层上的值，故称之为显式格式。2()njRO
12、hn+1nj-1j+1j111(1 2 )()nnnnjjjjur ur uu（10）2rh令，则(9)式可化为： 差分格式变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分2）格式）格式（BTCS）隐式格式）隐式格式对时间向后差分，对空间用中心差分，得： 2,0nnnhjjhjLuuu 111220nnnnnjjjjjuuuuuh111(12 )nnnnjjjjr ururuu注意：由（12）式不能直接计算出解，而要联立求解代数方程，故称之为隐式格式。2()njROhnn-1j-1j+1j（）变电站电气主接
13、线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分3）格式）格式 CrankNicoLson 格式（格式（CTCS）对时间和空间都用中心差分时间和空间都用中心差分，在 点对u作泰勒展开，得： 1( ,)2j n231,211( ,1)()22! 43! 2ttttttj nuj nuuuu231,211( ,)()22! 43! 2ttttttj nuj nuuuu (13)(14)变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分下面来
14、求 。在 对 点作泰勒展开： 上两式相加， 21( ,1)( , )1,( ,)2242ttttu j nu j nuj nuj n1( ,)2xxuj n1( ,)2j nxxu2,1,2( ,1)()28xxxxxx txx ttj nuj nuuu2,1,2( , )()28xxxxxx txx ttj nuj nuuu211( ,),1), )()22xxxxxxuj nuj nuj nO（(15)(16)(17)(18)变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分而： （19）（20）由（15
15、）式和（18）式得（21）或 ： （22） 22(1,1 -2 ( ,1)(1,1)( ,1)()xxu jnu j nu jnuj nO hh）22(1,-2 ( , )(1, )( , )()xxu jnu j nu jnuj nO hh）111111-11-1III2,2222102nnnnnnnnnjjjjjjjjhjuuuuuuuuLuhh111111112(2)2nnnnnnnnjjjjjjjjruuuuuuuun+1nj-1j+1j变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分注意：泰勒展开
16、点在格边上，不是在结点上，但在格式中未出现格边量。 全二阶精度。在 点展开时，用到了周围6个结点上的量，该格式又称为六点格式。隐式格式。 idea：是将微分方程中的 项以 在第n层和第n+1层上关于x的二阶中心差商的算术平均值来逼近，这一思想已被广泛地应用于一般微分方程，以建立其差分格式。 1222()njROh1( ,)2j n22ux( , )u xt变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分注意： 全二阶精度格式。三层显示格式。11112 (-2)nnnnnjjjjjur uuuu22()njR
17、Ohn+1nj-1j+1jn-1 4）格式（）格式（IV）（）（CTCS）（）（Richardson 格式）格式）对时间中心差分步长放大一倍，空间也中心差分。111-1IV,2202nnnnnjjjjjnhjuuuuuLuh(23)(24)变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分1）推广）推广CrankNicolson （格式（格式III）格式III将差分格式建立在 和 的中点 基础上的。现进一步推广，将差分格式建立在 和 中间任意一点上，即 ，其中 是一参数。按照格式III同样的方法进行差分。(
18、,1)j n( , )j n1( ,)2j n( ,1)j n( , )j n( ,)j n01( ,1)u j n( ,(1)u j n2233,11(1)(1)(1)2!3!ttttttj nuuuu ( , )u j n ( ,)u j n2233,11=( )2!3!ttttttj nuuuu 变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分上两式相减得同样，也可对 作类似格式III的处理。最后 （25） （26） ( ,1)( , )(12 )( ,)( ,)2!tttu j nu j nuj n
19、uj n 3232(1)( ,)( ,)3!3!ttttttuj nuj n 22ux11111-11-12222(1)0nnnnnnnnjjjjjjjjuuuuuuuuhh111111112(1)(2)nnnnnnnnjjjjjjjjuuruuuuuu 变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分注意：此格式利用了在第n层和第n+1层关于 的二阶中心差商的加权平均值，故 称之为加权六点格式。当 =0时，该格式变为格式I; =1时，该格式变为格式II; 时，该格式变为格式III;截断误差 ，当 时为 ；
20、当 时，为 ，故精度与 有关。2)修正修正Richardson格式格式(Dufort-Frankel格式格式)Richardson格式：x121212R22()O h2()O h11112202nnnnnjjjjjuuuuuh变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分而Dufort, Frankel给出的格式为： （27）实际上是将 换为 (由(2),(3)式相加可得) 经过这样的修改可将完全不稳定的Richardson格式变成无条件稳定格式。注意：截断误差： 要求 时速度比 的速度快，才能保证该格式
21、步长趋向于零时，逼近热传导方程，否则 为常数时，该格式逼近双曲型方程。 11111122nnnnnnjjjjjjuuuuuuh1111112 ()nnnnnnjjjjjjuur uuuunju111()2nnnjjjuuu242222( , )()()()ttRuj nO hOOhh00h /h格式名称格式表达式截断误差显式格式（格式I）隐式格式（格式II）CrankNicolson格式（隐式格式）220uutx小结：小结：模型方程 的若干典型差分格式2()Oh2()Oh22()Oh111(1 2 )()nnnnjjjjur ur uu2/rh111111(2)nnnnnjjjjjuur uu
22、u11111111(2) (2)2nnnnnnnnjjjjjjjjruuuuuuuuRichardson格式（显式格式）组合格式六点加权格式（显、隐与 有关） 当 时为 当 时为 Dufort-Frankel格式 复杂1111112 ()nnnnnnjjjjjjuur uuuu11111111() (1)(2nnnnnnnnjjjjjjjjuuruuuuuu 11112 (2)nnnnnjjjjjur uuuu22()Oh121222()O h2()O h变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分(2
23、) 直接差分逼近法直接差分逼近法tuutunjnjt10lim10limnnjjxuuuxx 利用中值定理，可得到：tnjnjnjnjtuttutuu22121nxjnjnjnjxuxxuxuu221212.1.3 差分格式构造方法差分格式构造方法众所周知，微商定义是： (1) Taylor展开法展开法变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分 其中 为 之间的常数。把这些表达式带入微分方程后，在一阶近似条件下略去小量 ，得到相应的差分方程。例如，采用上述方法可得到： nxjnjnjnjnjxuxxu
24、xuuu44222211121210njRnxjtnjnjnjnjnjnjxuxtutRxuutuu22221122 变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分FTBS)(11njnjnjnjuuruuFTCS)(2111njnjnjnjuuruu略去 后，得到差分格式： FTFS)(11njnjnjnjuuruu采用相同的直接差分逼近法，我们也可得到其他差分格式，例如： njR变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个
25、重要组成部分控制体积法不是从连续的微分方程出发，而是从物理量守恒规从物理量守恒规律出发建立离散的差分方程律出发建立离散的差分方程。考虑一维平面流动问题，设流体以速度 沿 轴正方向流动。流体中某一物理量，例如某物质浓度 ，应满足守恒定律。在空间位置 附近画出一个控制体积单元 ，如果控制体积内无源项，且不考虑扩散作用，在 内所含物质总量M应满足守恒定律，即：在控制体积 中物质总增量=由流进 内物质M的净通量。把物质浓度 近似看作是控制体积 内点j上的平均值，在时间 内控制体积 内总增量为：ux,Cx tjxdVdVdVdVdV,Cx tdVtdV,jjC x ttC x tx y z (3) 控制
26、体积法控制体积法变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分从左界面流入控制体积 内物质的平均通量速率为：在 时间内通过左边界面流入控制体积 内总通量为：同理，在 时间内通过右边界面流出控制体积 总通量是：流进控制体积 物质M净通量为流入通量与流出通量差： 由守恒定律得到：dV,2jxuC xty z tdV,2jxuC xty z t tdV,2jxuC xty z t dV,22jjxxu C xtC xty z t ,22jjjjxxC x ttC x txu C xtC xtt 变电站电气主接线
27、是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分在界面上函数值可以取相邻两个节点值的算术平均值，即：把它代入守恒方程，整理后得到 若用 标记，则可以得到差分方程： 1,22jjjxC xtC x tC xx t1,22jjjxC xtC x tC xx t,2jjjjutC x ttC x tC xx tC xx tx, j n1112nnnnjjjjutCCCCx变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分在上述推导中如果我们在界
28、面上函数值C采用不同的近似方式，则可得到不同的差分方程。一般而言，可取： 其中 。当 时为中心差分方程，当 时为前差方程，当 时为后差方程。用控制体积法构造差分方程总是守恒型差分方程。,1,2jjjxC xtC xx tC x t或 ,1,2jjjxC xtC x tC xx t011201变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分考虑对流方程，在矩形网格中控制体积 内，对时间和空间都取前差，即对流方程从 到 ，从 到 进行积分： dVntnttjxjxx 0txuau dxdt通过积分运算： 0jn
29、njjnnjxxttttxxtxxttxu dt dxau dx dt则得到：110jnjnxxttnnjjxtuu dxauu dt 采用积分方法构造差分方程基本思想是把微分方程在一定的控制体积内进行积分得到相应的差分方程。(4) 积分方法积分方法变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分 整理后得到： 这就是对流方程FTFS差分方程。如果采用不同的积分方式，就形成不同的差分方程。但是，如果积分区域不是矩形，则采用积分方法构造差分方程就不会这么简单，但基本原理是一样的。11()0nnnnjjjjuu
30、xa uut 11()nnnnjjjjuur uu把上述方程用数值积分近似表达式来表示，即得到：变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分除上述介绍的构造差分方程的方法外，还有其他构造差分方程的方法，这里不再一一介绍。综上所述，可以看出，对于同一个微分方程可以构造出多种多样的差分方程，而且，即使采用同一种方法也可以得到不同差分方程，不同的方法也可以导出同样的差分方程，这反映了构造差分格式的灵活性和灵活性和多样性。多样性。 变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电
31、任务。变电站的主接线是电力系统接线组成中一个重要组成部分初始条件： 第一边界条件： 第二边界条件： 第三边界条件 ：0( )tux(01)x12(0, )( )(1, )( )xxuttutt0tT 101212( )( )( )( )xxxxut utut ut0tT 12(0, )( )(1, )( )uttutt0tT 2.1.4 初边界条件的离散初边界条件的离散变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分下面仅考虑第三类边界条件的离散。 对初条件，可取：对边条件，有三种三种离散方法：方法一：方法
32、一： 在 用前差， 用后差，则： （28） 为 ，是一阶的。 由于控制方程差分格式的截断误差常达到 ，因此需要构造二阶精度的边界条件离散格式。0()jujh0,1,jN(0, ) t(1, ) t10101122()()()()nnnnnnNNNuununhuununh2()RO h( )O h变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分111011122()()2()()2nnnnnnNNNuununhuununh （29） 区域外的值 ， 需求出来 设 在边界上也成立，由格式I知：（30） 将（29
33、）式代入（30）式，消去 和 ，得 1nu1nNu0Lu 111(1 2 )()nnnnjjjjur ur uu10011111(1 2 )()(1 2 )()nnnnnnnnNNNNur ur uuur ur uu1nu1nNu1010111212122()22()122()22()nnnnnnNNNurhrnururhnurhrnururhn方法二：方法二：用两倍步长的中心差分离散，则 ：变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分方法三：方法三：采用中心差分，则有： 用 处的 代替 处的 ； 用 处的 代替 处的 ； 用 处的 代替 处的 ； 用 处的 代替 处的 。 进而有： 都为 ，二阶的。 2hx xu0 x xu(1)2hxNhxu1x xu2hx u0 x u(1)2hxNhu1x u1010111122()()2()()2nnnnnnnnNNNNuuuunnhuuuunnh2()RO h