2差分格式的构.docVIP

§2. 差分格式的构造 差分格式 对于描述流体流动的微分方程定解问题，用差分来近似方程及定解条件中的导数，就可以构造出该定解问题的差分格式。 【例】波动方程初边值定解问题 （注）根据特征线理论，边值函数 和 依赖于初始条件，不能随便给，否则会与初始条件矛盾，导致问题无解，这在数学上叫做定解条件的相容性。这里假定上面所给的定解条件是相容的。 计算区域的离散化：假设求解到 时刻，则问题的求解区域为 ， 。为了进行差分近似，取定自然数 （作为空间网格的数目）和 （作为时间网格的数目），空间 方向的步长取为 ，时间 方向的步长取为 ，记 ，（） ，（） 求解区域里像 这样的点称作网格点。 定解问题的离散化：在网格点 处列出方程 对方程中的二阶时间导数和二阶空间导数，都采用中心差分近似，即 代入原方程，就有 略去误差项，得 用 表示 的近似解（ ，），我们希望这一近似解能够使上述近似表达式还原成精确成立的等式（这样才能确切地使用它），也就是要求 满足 此式是一个差分方程，称为原方程的差分近似。 （注）因为是在网格点 处列出的方程，所以在近似时间导数时，空间自变量 的取值固定在 处，而在近似空间导数时，时间自变量 的取值固定在 时刻。 上述导出差分方程的过程比较繁琐，但我们可以把这一过程总结为： “用近似解的差分表达式代替原方程中的（偏）导数，得到其差分