早期,在科学技术的发展过程中,理论分析和实验研究一直是两种主要的科学方法,但是随着计算机的出现与发展,情况就大大不同了,研究者可以采用计算机去计算过去那些根本不能求解的问题,模拟一些难以观测的现象,第三种科学方法——数值计算诞生了。
有限差分法作为一种最为朴素的计算方法,发展时间也最为悠久,至今依然活跃。有限差分法的朴素思想就是,在求解连续性问题时,将时空域离散后,利用差商代替微商求解问题。
从应用上看,有限差分法常使用在求解时间依赖问题(双曲型和抛物型方程),而有限元主要求解稳态问题(椭圆型方程)。从建立差分格式的方法上看,有限差分法既可以根据Taylor级数展开建立差分格式,也可以根据积分方法建立差分格式,其中,积分方法由于其本身的特殊性,也被单独称为有限体积法。
从微分方程建立差分方程,方法格式层出不穷,但三个基本性质确是核心所在,即差分格式的相容性、收敛性和稳定性。
- 相容性:当时间和空间离散步长足够小时,差分方程将与微分方程充分的接近。
- 收敛性:当时间和空间离散步长足够小时,差分方程的解将逼近微分方程的解。
- 稳定性:有限差分格式是按照时间层逐层推进,如果初始计算引入的误差,随着逐层计算被控制,那么这个差分格式就是稳定的。
上述三个性质,通常是互相独立的,也就是相容的差分格式并不一定收敛。
通常,判断一个差分格式