1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。
该方法将 求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限
差分法以 Taylor 级 数展开等方法,把控制方程中的导数用网格节点上的函数值
的差商代替进行离散,从而 建立以网格节点上的值为未知数的代数方程组。该
方法是一种直接将微分问题变为代数 问题的近似数值解法,数学概念直观,表
达简单,是发展较早且比较成熟的数值方法。对于有限差分格式,从格式的精度
来划分,有一阶格式、二阶格式和高阶格式。从差分 的空间形式来考虑,可分
为中心格式和逆风格式。考虑时间因子的影响,差分格式还可 以分为显格式、
隐格式、显隐交替格式等。目前常见的差分格式,主要是上述几种形式 的组合,
不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步 长
一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本
的差分表达 式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分
和二阶中心差分等, 其中前两种格式为一阶计算精度,后两种格式为二阶计算
精度。通过对时间和空间这几 种不同差分格式的组合,可以组合成不同的差分
计算格式。
2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分
为有限个 互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数
的插值点,将微分 方程中的变量改写成由各变量或其导数的节点值与所选用的
插值函数组成的线性表达式 ,借助于变分原理或加权余量法,将微分方程离散
求解。采用不同的权函数和插值函数 形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每
个单 元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计
算域上总体的 基函数可以看为由每个单元基函数组成的,则整个计算域内的解
可以看作是由所有单元 上的近似解构成。
根据所采用的权函数和插值函数的不同 ,有限元方法也分为多种计算格式。从