[2018年最新整理]4-有限差分法-20161116.docxVIP

第4章 有限差分法目录第4章 有限差分法14.1有限差分法的原理24.1.1定解区域离散化34.1.2微分方程离散化34.1.3初边值条件离散化44.1.4解差分方程组44.2椭圆型方程的有限差分法54.2.1 微分方程的离散化61、五点差分网格62、九点差分网格74.2.2 边值条件的离散化81、第一边界条件离散化82、第三边界条件离散化94.3抛物型方程的有限差分法104.3.1一维抛物型方程的差分格式101、向前差分格式112、向后差分格式113、六点对称格式124、Richardson格式124.3.2高维方程的差分格式131、二维问题差分方法132、分裂格式算法144.4 双曲型方程的有限差分法164.4.1 一阶线性双曲型方程的差分格式161、一阶线性双曲型方程初值问题162、一阶线性双曲型方程初边值问题184.4.2 二阶线性双曲型方程的差分格式191、波动方程及其特征192、二阶方程的差分格式203、定解条件的处理21参考文献22第4章 有限差分法刚刚介绍的变分法和加权余量法，都是将实际问题的控制方程转化为积分形式进行求解，而有限差分方法不同，它是用网格覆盖求解区域和时间域，用差分近似替代控制方程中的微分，把求解微分方程的问题改换成求解代数方程问题，进行近似求解。有限差分法是一种经典数值方法，当网格细化或节点较多时，近似解的精度可以得到改进，对于几何形状规则且边界条件简单的问题，差分法有很高的求解精度和收敛性，但用于求解几何形状和边界条件复杂的问题时，它的精度将会降低，甚至产生困难。有限差分法在流体力学和爆炸力学中得到广泛应用，也曾用以求解弹性力学边值问题。目前，弹性静力学一般都不用差分法进行求解了，但是对于弹性动力学等初值问题，仍然要借助于有限差分法予以解决。有限差分法的关键是构造差分格式，保证解的收敛性和解的稳定性。这里简要介绍椭圆型方程、抛物型方程和双曲型方程的差分格式构造方法。4.1有限差分法的原理差分法是微分方程的一种近似数值解法。它不是去寻求函数式的解答，而是去求出函数在一些网格结点上的数值。具体地讲，差分法就是先将求解域离散化，在这些离散点上用差商近似地代替导数，从而把基本方程和边界条件近似地改用差分方程来表示，把求解微分方程的问题改换为求解代数方程的问题。考虑一维热传导方程的第一类初边值问题 （4-1）其中，，a为正常数。假定和在相应区域光滑，并且在满足相容条件，使上述问题有唯一充分光滑的解。用差分法求解上述问题时，通常情况下，首先采用矩形网格对求解区域进行离散化，接着，利用差商或Taylor展式对微分方程和初边值条件等连续性方程，在网格节点处做离散化处理，最后解差分方程组，获得所有网格节点处的解的近似值，用近似研究微分方程初边值问题的解。4.1.1定解区域离散化由于有限差分法求解偏微分方程问题必须把连续问题进行离散化，因此首先要对求解区域进行网格剖分。对于不同的求解问题，求解区域也各不相同。对于问题（4-1）而言，如4-1所示，设N是一给定的正整数，取空间步长=l/N，再取一定的时间步长=T/M，做两族平行线，用平行直线族把区域分成若干个小矩形。位于区域D内部的网格节点（简记为）称为内节点，所有内节点用表示。位于边界上的节点称为界点，所有界点记为 。图4-1 差分法矩形网格剖分4.1.2微分方程离散化由微分学知道，函数的导数是函数的增量与自变量增量之比的极限，即当很小时，可以近似地用差商或代替，可以近似地用二阶差商代替，从而一个微分方程可近似地用一个差分方程来代替。因此，对于问题（4-1）而言，其一维热传导方程可以用方程 （4-2）近似代替，其截断误差为。差分格式（4-2）称为解热传导方程（4-1）的古典显格式，它是一个两层4点格式，所用到的节点图式如图4-2所示。图4-2 古典显式差分格式节点图在有限差分法中，为简单起见，通常将节点处的函数值简记为，依次，节点处的函数值简记为，节点处的函数值简记为。这样的话，式（4-2）就可以写成 （4-3）4.1.3初边值条件离散化对于问题（4-1）中的初边值条件，在网格节点处取值，有4.1.4解差分方程组对于问题（4-1），将微分方程和初边值条件的离散化方程联立，得差分方程组 （4-4）这样，利用k=0层上的数值，即可逐点算出k=1层上的全部离散点处解的近似值，再利用k层上的节点值算出k+1层上的节点值。这样逐层逐点地计算数值解是非常方便的。同样地，若给定的是初值问题，仍采用前面的古典显格式和初值离散方法，也容易写出其差分方程组，但其计算求解范围与初边值问题情形有所不同。此时，利用k=0时间层上的节点近似值只能计算出k=1时间层上关于j=1,2，N-1处的节点值近似值，依次下去，纯初值问题的差分解被局限于图4-1中0A1上的网格节点处。4.2椭圆型方