离散方程的误差与物理特性的分析.pptxVIP

离散方程的误差与物理特性的分析

2、离散方程得相容性当空间、时间步长趋于零时,如果离散方程得截断误差趋于零,则称离散方程与微分方程相容。截差,方程具有相容性。如截差中出现Dt/Dx时,相容性仅条件满足。

二、离散误差及收敛性1、离散误差离散方程得精确解偏离该点上相应微分方程精确解得值。2、影响因素(1)截断误差一般地说,截断误差得阶数越高,离散误差越小。(2)网格疏密对同一离散格式,网格加密离散误差也会减小。但:网格太密也会引起很大得舍入误差。网格独立得解就是国际杂志接受数值计算论文得基本要求。层次:解

3、离散方程得收敛性当时间与空间步长均趋于零时,如果各节点上得离散误差都趋于零,则称离散方程收敛。

三、舍入误差与初值问题得稳定性1、舍入误差由于计算机存储数据得方式而引入得误差,取决于所采用得计算方法及所用计算机得字长。2、数值解误差得组成总误差=离散误差+舍入误差3、初值问题离散格式得稳定性如果在任一时层计算中所引入得误差都不会在以后各时层得计算中被不断地放大,则称离散格式就是稳定得。稳定与不稳定就是一个离散格式得固有属性！

4、联系收敛性与稳定性得Lax原理对于适定得线性初值问题所建立起得相容格式,稳定性就是收敛性得充分必要条件。工程传热问题中,仅常物性无源项得非稳态导热问题,才能应用Lax原理分析。对于其她问题,相容性与稳定性仅就是获得收敛解得必要条件。

格式、网格疏密、稳定性得影响

5、说明对于工程对流换热问题得数值计算,一般认为:扩散项采用二阶精度截差,对流项采用二阶到三阶得离散格式就是比较合适得。数学特性:相容性、收敛性、稳定性物理特性:守恒性、迁移性、人工粘性(假扩散)在数学上满足稳定性条件得解,不能保证一定得出具有物理意义得结果。

§3、2离散方程得守恒性一、离散方程守恒性定义如果离散方程在定义域内得任一有限空间内作求与得运算,所得表达式满足该区域上物理量守恒得关系,则称该离散格式具有守恒特性。二、对流项中心差分得守恒特性扩散项中心差分具有良好得物理特性与计算精度,就是对流项离散格式导致了各种不良特性。将一维对流扩散方程简化成纯对流方程(平流方程)

二、对流项中心差分得守恒特性(续)离散化:在[l1,l2]内作积分或:

二、对流项中心差分得守恒特性(续)互相抵消

大家学习辛苦了，还是要坚持继续保持安静

二、对流项中心差分得守恒特性(续)满足守恒特性！

三、保证离散方程具有守恒性得条件1、导出离散方程得控制方程就是守恒型得非守恒型守恒型

三、保证离散方程具有守恒性得条件(续)不能抵消,表明非守恒型方程不具有守恒特性。

三、保证离散方程具有守恒性得条件(续)2、在同一界面上各物理量(f及有关物性)及f得一阶导数就是连续得。连续:从界面两侧得两个控制容积来写出得该界面上得值就是相等得(局部守恒性)。表示为:表示对P点写出得e界面上得f值。

三、保证离散方程具有守恒性得条件(续)界面连续性示意图

三、保证离散方程具有守恒性得条件(续)中心差分满足连续性条件,具有守恒特性。四、讨论1、椭圆型问题最希望具有守恒特性。2、抛物型及双曲型问题未必具有守恒特性。

§3、3离散方程得迁移性一、对流与扩散现象在物理本质上得区别扩散项应使扰动向四周传递。对流项应使扰动沿流动方向传递。迁移性！

二、扩散项中心差分可以将扰动均匀地向四周传递一维非稳态扩散方程采用显式格式离散采用离散扰动分析法:设n时层i节点上出现一个扰动e,其余各点均没有。对i点:整理得:,按稳定性要求

扩散项中心差分得迁移性对i+1点:整理得:对i-1点:整理得:

扩散项中心差分得迁移性时扩散作用下扰动得传递扩散项中心差分既具有守恒特性,又能使扰动均匀地向四周传递,符合物理问题得本质,就是一个理想得离散格式。