有限差分方法FiniteDifferenceMethod-大丰区新丰中学学校.PDF

第 一章 有 限差 分方 法 ( F in it e D iffer en ce M et h o d ) 本章介绍一些有限差分方法的基本知识。先介绍差分格式的 相容性、收敛性与稳定性; 差分格式 的截断误差与精度; 水利水 电 工程 中常用的差分格式; 显式差分格式与隐式差分格式之应用; 判 别差分格式稳定性的方法。然后介绍预测-校正二步法差分格式; 有限体积法之应用。最后介绍两个应用差分方法计算实际工程 问 题的实例 。 1. 1 概 述 对于求解的偏微分方程 的定解 问题, 通 常都可 以用差分方法 得到近似解 。这是因为导数是作为函数的增量与 自变量的增量之 比的极限, 即 u ( x + Δx ) - u ( x ) u ′( x ) = lim Δx →0 Δx 因此, 在求偏微分方程的近似解时, 可 以用差商 u ( x + Δx ) - u ( x ) Δx 代替 ( ) 的导数, 这样就把偏微分方程的求解 问题化为求解代数 u x 方程组的问题, 这就是差分方法。 差分方法的主要步骤, 是将定解区域用网格线进行剖分, 将计 算区域离散, 化为离散点集, 求解网格线结点上微分方程解的近似 值; 用差商代 替偏微分方程的偏导数, 将偏微分方程化为差分方 程, 并将定解条件离散, 这个过程通常叫做差分格式的构造 。建立 ·2 · 水利枢纽河段模拟 与采用差分格式后, 就将偏微分方程的定解 问题化为代数方程组 的求解 。采用不同的差分格式, 得到不同的代数方程组。通过解代 数方程组, 求得定解 问题在网格点上的近似解; 再应用插值方法, 便可从网格点上的近似解, 得到整个定解区域上的近似解 。 有限差分方法是求解各类偏微分方程和方程组的很通用的方 法, 也是求解水利、水 电工程流体 问题常用的基本方法。但在应用 有限差分方法求解偏微分方程定解 问题时, 会遇到许多重要的、有 时还是相当困难的数学问题 。例如, 有限差分的误差与精度 问题、 收敛性与稳定性 问题 以及差分格式的构造 问题等。 1. 2 有 限差分的格式 一、差分格式的构造 差分格式可利用数值微分的方法, 将微分方程各个微商用适 当的差 商 近 似 地 代 替 而 得 到。也 可 以从 泰 勒 公 式 ( T aylor ′s E xp an s ion ) 推得。 ( 1) 数值微分法构造差分格式 通常从微分方程出发, 利用数值微分的方法, 即将各个微商用 适当的差商近似地代替, 从而得出能够任意逼近微分方程的差分 格式, 这是建立差分格式最常用、最直接的一种方法。这种方法 比 较简单, 而且容易直接列出差分格式, 对双 曲型和抛物型方程更为 常用。 现有一维对流方程的初值 问题 u u + c = 0 , x ∈R , t 0 t x ( 1. 1) u ( x , 0) = f ( x ) , x ∈R 式中, c 是常数。 设 u 是式( 1. 1) 的解, 对于任何结 点( i, n) , u 的微商与向前差 第一章 有限差分方法 ·3 · 商之间有如下关系式: u ( x i, t n+ 1 ) - u ( x i, t n)