有限差分方法.pdfVIP

有限差分方法 许多物理现象或其运动、演化过程可以用一个微分方程的定解问题来描述，例如无限长 细弦的自由振动向题可归结成二阶双曲型方程的初值问题，而弦对平衡位置的偏移就是方程 的解。但是，绝大多数偏微分方程定解问题的解通常并不能用显式的公式来表达，有时即使 可用公式表示，也往往因为过于复杂，而需要采用各种近似方法来计算它的解。差分方法就 是求解（偏）微分方程定解问题的常用近似方法之一。Courant, Friedrichs, Lewy (1928) 首次 对偏微分方程的差分方法作了完整的论述。第二次世界大战以来，快速电子计算机的诞生为 差分方法提供了强有力的工具，从而促使这一数值分析方法迅速地发展起来。 有限差分法的基本思想是用离散的、只含有限个未知数的差分方程去代替连续变量的微 分方程和定解条件。对于求解的偏微分方程定解问题，有限差分方法的主要步骤如下：利用 网格线将定解区域化为离散点集；在此基础上，通过适当的途径将微分方程离散化为差分方 程，并将定解条件离散化，这一过程叫做构造差分格式，不同的离散化途径一般会得到不同 的差分格式；建立差分格式后，就把原来的偏微分方程定解问题化为代数方程组，通过解代 数方程组，得到出定解问题的解在离散点上的近似值组成的离散解，应用插值方法可从离散 解得到定解问题在整个定解区域上的近似解。由此可见，有限差分方法有大体固定的模式， 它有较强的通用性。但是，不能误认为不去了解这种逼近方法的基本知识，只是单纯模仿， 便能轻易获得满意的结果。因为在应用这种逼近方法时会发生许多重要的但有时还是相当困 难的数学问题，包括精度、稳定性与收敛性等。 本章的目的是介绍有限差分方法的一些基本概念和构造差分格式的基本方法，以及给出 一些力学问题有限差分求解的全过程。 §1 有限差分的基本概念 在本节中，我们将首先针对单变量系统，基于 Taylor 展式导出若干有限差分表达式。接 着，将所得结果推广到两个自变量以及多自变量函数系统。 §1.1 函数的表示 首先考虑单变量函数u(x) ，x 为自变量。我们将区域x 离散化为一系列点（或节点，有 时也称作结点）（如图 1.1），使得 u(x ) u(rh) u ，（r 0, 1, 2, ） （1.1） r r u u u(x ) u(rh) r r u(x) h 2h 3h 4h 5h rh x 图1.1 等格距h 时函数u(x) 的离散化 在用rh 代替x 后，节点坐标就仅仅由整数r 和格距h 的乘积给定(这里不妨假设h 为常 r 数并规范化为小于 1)，其中h 称为沿x 方向的步长。整数r 表示节点沿x 坐标的位置，通常， 当x 0 时r 0 。当h 是常数时，u(rh) 就可表示为ur 。 (a) (b) 图1.2 （a ）二维有限差分离散网格；（b ）三维有限差分离散网格 在二维情形下( 图 1.2a)，函数u(x