数值计算在热工中应用B第三章要点.ppt

μ = 1 ? 2r + 2r cosθ ? Ic sinθ 圆心 长轴半径 短轴半径 (b) 2r ≤ 1, c 1 c 1 (c) 2r ≤ 1, c ≤ 1 c 1 椭圆右端点处曲率半径 ≤ 1 从数学公式分析：本例椭圆的参数方程为 曲率半径为 右端点为θ = 0 ，要求右端点曲率半径 ≤ 1 ，得 c ≤ 2r c ≤ 2r 一维模型方程的FTCS格式的稳定性条件为 2 0≤ 2r ≤ 1; c ≤ 2r 2 3.2.5 von Neumann 分析方法应用讨论 1.适用于第一类边界条件的线性初值问题，可得出允 许的最大时间步长； 2.对非线性初值问题（瞬态NS方程）采用局部线性 化的近似：按线性方式分析，得出的时间步长打一个 折扣 3. 是一种有效的分析工具：多重网格分析，代数方 程迭代求解的稳定性分析。 3.3 离散方程的守恒性 3.3.1 定义与分析模型 3.3.2 直接求和分析法 3.3.3 保证离散方程具有守恒性的条件 3.3.4 讨论－追求而不唯一 3.3 离散方程的守恒性 3.3.1 定义与分析模型 1. 定义 2. 分析模型 在模型方程的离散中扩散项采用中心差分，其具 有守恒性，分析仅对对流项的离散格式进行－非稳 态项与对流项组成“平流方程”(advection)： 将一个离散方程在任意有限容积内进行求和， 如所得结果仍满足守恒关系，则称其具有守恒性。 （守恒型） （非守恒型） 3.3.2 直接求和分析法 将守恒型的FTCS格式对区间[ ]求和： l 1 , l 2 空间项为中心 差分未加上标 在 Δt 时间间隔内区域 [l1 , l2 ] 中 φ 的增量 是否等于相同时间间隔内流入、流出之差？ 需对上式右端做分析： * 数值计算在热工中应用 B 天津理工大学自动化学院热能系 第三章 离散方程的数学与物理特性分析 3.1 离散方程的相容性、收敛性与稳定性 3.2 分析初值问题稳定性的von Neumann方 法 3.3 离散方程的守恒性 第 3 章教学目录 3.4 离散方程的迁移特性 3.1 离散方程的相容性、收敛性与稳定性 3.1.1 截断误差及相容性 3.1.3 舍入误差与初值问题的稳定性 3.1.4 数值特性分析举例 3.1.2 离散误差与收敛性 3.1 离散方程的相容性、收敛性与稳定性 3.1.1 截断误差及相容性（consistency) 1. 离散方程的精确解 2. 微分算子与差分算子 在求解过程中不引入任何舍入误差的解，记为： 假定可以对其做Taylor展开。 （1）微分算子－ L(φ )i ,n 对函数Φ在某点（i,n) 做某些微分运算及算术运算的算子， 如： 同一地点上差分算子与微分算子之差。 某些差分及算术运算的算子，如 (2) 差分算子－ LΔx ,Δt (φi,n ) 对函数 在某点（i,n)做 则 LΔx , Δt (φi,n ) = 0 即为1-D模型方程的显式格式（时 间向前，空间二阶格式－FTCS）。 3. 离散方程的截断误差 （2）分析方法—对离散方程的精确解 ， 对 点（i,n) 做Taylor展开（含时间与空间，代入离散方程并整理成两个算子相减的形式。 对1-D模型方程的FTCS格式，可有 （1）－定义— 如瞬态项 即 类似地 故得 假设源项不存在截断误差，则： 1-D模型方程的FTCS格式的截断误差即为： O（Δt,Δx2） 4. 离散方程的相容性 当 Δt → 0, Δx → 0 如果离散方程的截断误差趋于 零，则称该离散方程与微分方程相容。 当离散方程的截断误差表示为 O(Δtn , Δxm)(n, m 0)时 离散方程具有相容性；当截断误差含有 项时，只有 当 Δt 比 Δx 更快地趋于零时，才具有相容性。 3.1.2 离散误差与收敛性 1. 数值解的离散误差(discretization error) 微分方程的精确解 离散方程的精确解 2. 影响离散误差的因素 （1）截断误差：截差阶数越高，同样网格下， 越小 （2）网格步长：同样截差下，网格加密，误差减小； 对一般工程计算格式截差选取建议： 扩散项－2阶，对流项－2、3阶。 3. 离散方程的收敛性 当 Δt → 0, Δx → 0 如果 → 0 则离散方程收敛。 收敛性的证明并不容易，对线性问题有Lax原理。 3.1.3 舍入误差(round-off error)与初值问题的稳定性 1. 数值解的舍入误差 从计算机实际得到的解。 2.