数值传热学3.pdfVIP

数值传热学 第三章 离散方程的数学与物理特性分析 主讲 陶文铨 西安交通大学能源与动力工程学院 热流中心 CFD-NHT-EHT CENTER 2011年9月19 日, 西安 1/61 第3 章教学目录 3.1 离散方程的相容性、收敛性与稳定性 3.2 分析初值问题稳定性的von Neumann方 法 3.3 离散方程的守恒性 3.4 离散方程的迁移特性 2/61 3.1 离散方程的相容性、收敛性与稳定性 3.1.1 截断误差及相容性 3.1.2 离散误差与收敛性 3.1.3 舍入误差与初值问题的稳定性 3.1.4 数值特性分析举例 3/61 3.1 离散方程的相容性、收敛性与稳定性 3.1.1 截断误差及相容性（consistency) 1. 离散方程的精确解 在求解过程中不引入任何舍入误差的解，记为： n  , 假定可以对其做Taylor展开。 i 2. 微分算子与差分算子 （1）微分算子－ L ()i ,n 对函数  在某点（i,n) 做某些微分运算及算术运算的算子， 如： 2         L ( )i,n ( u 2 S )i ,n t x x 4/61 则 L ()i ,n 0 即为一维模型方程。 n 对函数 n 在某点 （i,n)做 (2) 差分算子－ L ( )  x t i  , i 某些差分及算术运算的算子，如 n1 n n n n n n n      2  n L ( )  i i u i1 i1  i1 i i1 S x t i i  , t 2x x2 则 n L , ( ) 0 即为1-D模型方程的显式格式（时 x t i 间向前，空间二阶格式－FTCS）。 3. 离散方程的截断误差 Forward time and central space 同一地点上差分算子