高等油藏物理第7章差分格式的构造.pptVIP

考虑一维平流问题，设流体以速度a沿x轴正方向流动。流体中某一物理量，如一物质浓度u(x,t)在流动中满足守恒律。 * 第7章 差分格式构造方法 1.泰勒级数展开法 对于一个偏微分方程，可以建立不同的差分方程，它们的解都是原偏微分方程的近似解，可以用不同的方法得到相同的差分方程。 双曲型一维平流方程及初值条件 空间步长。 时间步长。 网格节点 简记为 。函数值记为 对函数u在空间做泰勒级数展开 一阶中心差商 当 趋于零时，截断误差 也趋于零，因此说差商与微商是相容的。 一阶偏导数的中心差商表达。它具有 阶的截断误差，记为 或者说距离有二阶精度。 具有一阶精度，差商与微商是相容的。 一阶向前差商 一阶向后差商 具有一阶精度，差商与微商是相容的。 二阶中心差商 具有二阶精度，差商与微商是相容的。 对时间的一阶向前差分 用差商代替微商，双曲型一维平流方程及初值条件可写为三种形式。 （1）中心差分格式（中心） 写为 格式截断误差 （2）向前差分格式（右偏） 格式截断误差 （3）向后差分格式（左偏） 格式截断误差 只要已知第n层的值就可以计算n+1层上的值，这样从初始条件开始可以逐层计算下去，不必求解方程组。这种格式叫做显式。 构成差分格式的泰勒级数展开法是一种最常用的方法。它简单但