d_f问题分析(ppt).pptVIP

2003-5-24 DuFork-Frankel格式分析----曲径 提纲 算法的提出 实例验证 理论分析 相容性分析 稳定性分析 结论 算法的提出 扩散方程 的 Richardson格式是一个二阶精度的三层显格式： 并且我们已经知道其截断误差为 阶数高于古典显格式和古典隐格式，但是由于它无条件不稳定，不能应用在计算在中。 算法的提出 课本中为我们提供了一种Richardson格式的改进格式DuFork-Frankel格式，用 代替了 ，这仍然是一个三层格式，用到四个点的值: (i-1,j)，(i+1,j)，(i,j+1)，(i,j-1) ，但是我们的课本中只提到了这个差分格式，却没有进行任何分析，计算或证明，为了验证它对Richardson格式改进后的真实效果，下面就来讨论DuFork-Frankel格式。 算法的提出 实例验证 具体的差分格式，我们只考虑最基本的抛物型方程 应用DuFork-Frankel格式进行如下计算。 实例验证 首先，不固定网比 r，而是选用h=1/16， 分别取1/160，1/16和1/4，也就是使 和 h 的比值小于1，等于1和大于1。这时候计算结果发生了如下现象。 实例验证 实例验证 实例验证 实例验证 由于上面三个问题的网比r都大于1，我们考虑这样的不收敛于原抛物型方程的结果是否是因为差分格式不稳定呢？ 那么，现在固定 r（固定时间步长和空间步长平方的比值），再进行一次计算。 实例验证 固定网比 r，考察一下是不是对于任何r,DuFork-Frankel格式都稳定？ 分别取r=0.4，r=0.8，r=1.2并和精确解比较，可得到如下的结果。 实例验证 实例验证 上述计算可以得到相对误差如下： 可见，相对误差随网比缩小而变小。并且，在时间方向并没有发生如前面两个图的不收敛于抛物方程的情况。 实例验证 既然此时DuFork-Frankel格式对于任何的 r 都是稳定的，为什么会有前面不收敛的现象呢？？ 理论分析 若想得到DuFork-Frankel格式的收敛条件 就要从相容性和稳定性共同分析。 相容性 相容性 得到截断误差为 说明DuFork-Frankel格式要求 ,h 趋向0，并且 趋向于0的速度要比 h 快，才能保证格式与最初的抛物型方程相容。 和我们所遇到的其他差分格式不同，DuFork-Frankel格式的相容性还与 和 h 的比值有关。 相容性 不妨设想 时（c为常数），由上述Taylor级数可以得到这时差分在 和 h 趋向于0的时候，DuFork-Frankel格式趋向于了 这是一个非线性的双曲方程。 这里就解释了为什么在 发生变化的时候，差分格式不收敛于原方程了。 实例验证 稳定性 对于稳定性的讨论，类似于讨论Richardson格式，利用传播因子法，引入向量 其中 满足： 并引入 ，得到： 稳定性 得到传播矩阵 其特征值为 则 ，由Von Neumann条件， DuFork-Frankel格式是无条件稳定的。和我们在前面例子中看到的一致。 结论 在考虑差分格式的收敛性时，一定要兼顾相容性和稳定性。 对于我们所接触到显格式，都在或稳定性，或相容性上存在限制条件。 * * 2003-5-24 DuFork-Frankel格式分析----曲径 曲径 0018138 DuFork-Frankel格式为： 由于DuFork-Frankel格式是一个三层格式，在计算中初始条件后的第一层是不能用DuFork-Frankel格式得到的，所以选择了二阶精度的Crand-Nicholson格式先求出第一层，再进行DuFork-Frankel格式的计算。 就象Lax所证明的那样： 相容性 + 稳定性 = 收敛性 2003-5-24 DuFork-Frankel格式分析----曲径 将上式两项分别展开成Taylor级数 1. 2. *