格式分 曲径 0018138
曲径 0018138
■算法的提出 实例验证 理论分析 相容性分析 一稳定性分析 结论 DuFork- Franke格式分析一曲径 2
DuFork-Frankel格式分析----曲径 2 提纲 ◼ 算法的提出 ◼ 实例验证 ◼ 理论分析 – 相容性分析 – 稳定性分析 ◼ 结论
扩散方程 的 Richardson格 式是一个二阶精度的三层显格式 1+t21+1 2L!+ 2T h 并且我们已经知道其截断误差为R=Ox2+h2) 阶数高于古典显格式和古典隐格式,但是由于 它无条件不稳定,不能应用在计算在中。 DuFork- Franke格式分析一曲径 3
DuFork-Frankel格式分析----曲径 3 算法的提出 扩散方程 的 Richardson格 式是一个二阶精度的三层显格式： 并且我们已经知道其截断误差为 阶数高于古典显格式和古典隐格式，但是由于 它无条件不稳定，不能应用在计算在中。 2 1 1 1 1 2 2 h u u u a u u j i j i j i j i j i + − + + − + = + 2 2 x u a t u = R O( h ) j 2 2 i = +
课本中为我们提供了一种 Richardson格式的改 进格式 DuFork- Frankel格式,用u+u 代替了2u},这仍然是一个三层格式,用到 四个点的值:(i-1,j),(i+1,j),(i,j+1), (i,j-1),但是我们的课本中只提到了这个差 分格式,却没有进行任何分析,计算或证明, 为了验证它对 Richardson格式改进后的真实效 果,下面就来讨论 DuFork- Frankel格式。 DuFork- Franke格式分析一曲径
DuFork-Frankel格式分析----曲径 4 算法的提出 课本中为我们提供了一种Richardson格式的改 进格式DuFork-Frankel格式，用 代替了 ，这仍然是一个三层格式，用到 四个点的值: (i-1,j)，(i+1,j)，(i,j+1)， (i,j-1) ，但是我们的课本中只提到了这个差 分格式，却没有进行任何分析，计算或证明， 为了验证它对Richardson格式改进后的真实效 果，下面就来讨论DuFork-Frankel格式。 j 1 i j 1 i u u + − + j 2ui
DuFork- Frankel格式为 (u×1 1i+1 -a-1+1-(u fu +u 0 2τ h 由于 DuFork- Frankel格式是一个三层格式,在 计算中初始条件后的第一层是不能用 DuFork Frankel格式得到的,所以选择了二阶精度的 rand- Nicholson格式先求出第一层,再进行 DuFork- Frankel格式的计算。 DuFork- Franke格式分析一曲径 5
DuFork-Frankel格式分析----曲径 5 算法的提出 0 2 h j i 1 ) u j 1 i u j 1 i (u j i 1 u a 2τ ) j 1 i u j 1 i (u = − + + + − − + − − − + DuFork-Frankel格式为： 由于DuFork-Frankel格式是一个三层格式，在 计算中初始条件后的第一层是不能用DuForkFrankel格式得到的，所以选择了二阶精度的 Crand-Nicholson格式先求出第一层，再进行 DuFork-Frankel格式的计算