《计算物理学》课件 第三章 偏微分方程.pptx

1第三章偏微分方程的数值求解方法

偏微分方程?2?

双曲型方程3xt?即对流方程，这里a为常数。顾名思义，对流方程描述流体中某种物质的物理量的变化规律，例如传热过程中温度的演化规律，流体中溶质的浓度等的变化规律等。?

有限差分方法在有限差分方法中，我们放弃了微分方程中独立变量可以取连续值的特征，而关注独立变量离散取值后对应的函数值。原则上，这种方法仍然可以达到任意满意的计算精度，因为方程的连续数值解可以通过减小独立变量离散取值的间隔，或者通过离散点上的函数值插值来近似得到。有限差分法具体操作分两部分：1.用差分（差商）代替微分方程中的微分（微商），将连续变化的变量离散化，从而得到差分方程组的数学形式。2.求解差分方程组。4

在第一步中，通过所谓的网络分割法，将函数定义域分成大量相邻而不重合的子区域，通常采用的是规则的分割方式。网络线划分的交点称为节点。数值求解的关键就是要应用适当的计算方法，求得特定问题在所有这些节点上的近似值。一个函数在x点上的一阶和二阶微商，可以近似地用它所临近的两点上的函数值的差商来表示。具体方法5

差分格式的建立hτtxO差分法是处理偏微分方程的一种常用数值解法6??

7其次，分别利用一阶微商的向前和向后差商公式，以及中心差商公式构造偏微分方程的差分格式。?将这些差商公式分别代入方程，即可构造出不同的差分格式，并建立相应的差分方程。

8??k,jk+1,jk,j+1??

9?k-1,jk,jk,j+1k+1,j?k-1,jk,jk,j+1

10???

11?

12差分格式的相容性和收敛性从上面建立的差分格式和他们相应的截断误差可以看出，当步长h,τ?0时，各差分格式的截断误差也?0。这表明：差分格式的极限形式就是相应的微分方程这称为：差分格式和相应的微分方程是相容的通常是希望用差分格式计算出来的解在h,τ?0时，能逼近于相应微分方程的解，这就是收敛性问题。相容性和收敛性的概念是不同的，相容是必备条件，收敛才是目的。一般来说相容并不能保证收敛，但在一定条件下，两者等价。

差分格式的收敛性–依赖区域t=0t=1t=2t=3BCP若用左图差分格式计算P点的值，依赖于初始线段BC上的网络节点的值；对于右图，则依赖于初始线段AC上节点的值。BC就是格式左的依赖区域；AC就是格式右的依赖区域t=0t=1t=2t=3BCPA13依赖区域跟：差分格式、时间层相关差分格式的收敛性与偏微分方程的依赖区域有关。

Courant条件差分格式收敛的一个必要条件是，差分格式的依赖区域应包含微分方程的依赖区域。——Courant条件14?

可见：可通过控制t的步长τ和x的步长h来保证收敛。15?

差分格式的稳定性数值不稳定性：用差分格式求初值问题数值解，每步计算都会引入舍入误差，而初值问题还是逐层计算的，所以误差还会逐层传播。若误差传播越来越大，就可能会淹没真解。16?差分格式的稳定性不但与差分格式有关，还与网格步长的比值有关。如果在一定步长下，差分格式是稳定的，则称该差分格式为条件稳定。如果对任何步长，差分格式都是稳定的，则称该差分格式为无条件稳定。如果差分格式对任何步长都不稳定，则称该差分格式为完全不稳定。

17一种讨论稳定性的常用方法—Fourier方法常用Fourier方法(也称VonNeumann方法)研究稳定性。??

18??

19????

20??精度稳定性条件逆风差分格式FTCS差分格式不稳定蛙跳差分格式Lax差分格式Lax-Wendroff差分格式

Lax等价定理(相容性、稳定性、收敛性的关系)：对于一个适定的初值问题，一个与它相容的差分格式收敛的充分必要条件是这个格式的稳定性。21对策：

22小结有限差分求解偏微分方程：离散化?建立差分方程讨论差分方程的稳定性和收敛性具体地，离散化、构造网格、用网格函数做近似解、用差商代替微商、同Talyor展开建立差分方程、或用特征线构造差分方程截断误差的引入，相容的概念差分方程的稳定性就是舍入误差是否会逐层传播无限增加，判断其稳定性可以用图表，也可以用Fourier方法Lax定理和Courant条件

例：波动方程23略去小量，并代入方程可得，

24

25????

262、利用二阶泰勒公式构造第二行的值。?这样就得到了一个改进的第二行的近似差分公式。

练习：求解两端固定的弦上的振动形式27

28抛物型方程的差分解法一维热传导方程：?????根据给定的初值和边值，该方程有以下两类定解条件：

292.初边值混合问题：也就是求解函数在一个有限区域内既满足初始条件也满足边界条件的定解。???

30????

31?如果第j层的近似值已经求得，通过以上公式可以计算得到网格中第j+1层的近似值。该差分格式称为显式差分格式。?

32??

33