抛物型方程差分格式特性的研究.docVIP

抛物型方程差分格式的特性实验 一、研究对象 一维热传导方程是一种简单的双曲方程，我们研究如下奇次方程： 该定解问题的精确为： 二、研究方法 1 向前差分格式 差分方程： ⑴ 边界处理： 由问题可知： C程序实现： 1 第0层时间步的值 (j=0、1…jmax) 2 当n=0、1…nmax-1时，令 ① ② (j=1、2…jmax-1) 2 向后差分格式 差分方程 (2) 边界处理： 由问题可知： C程序实现： 1 第0层时间步的值 (j=0、1…jmax) 2令， a=-λ，b=2λ+1，c=-λ (j=2、3…jmax-1) 3当n=0、1…nmax-1时， ① ③ (j=2、3…jmax-2) (j=jmax-2、jamx-3…1) 3 Crank-Nicolson 差分方程 (3) 边界处理： 由问题可知： C程序实现： 1 第0层时间步的值 (j=0、1…jmax) 2令， a=-λ，b=2λ+1，c=-λ (j=2、3…jmax-1) 3当n=0、1…nmax-1时， ① ③ (j=2、3…jmax-2) (j=jmax-2、jamx-3…1) 4 Du Fort-Frankel 差分方程 (4) 边界处理： 由问题可知： C程序实现： 1 第0层时间步的值 (j=0、1…jmax) 2 利用Crank-Nicolson计算出第一层时间步的值。 3 当n=1、2…nmax-1时， ① ② (j=jmax-2、jamx-3…1) 三、研究结果 1 精确解 2 向前差分的结果 3 向后差分的结果 4 Crank-Nicolson格式的结果 5 Du Fort-Frankel格式的结果 6当λ=0.1，t=0.1时，准确解、向前、向后、C-N及DFF格式的比较 当λ=0.5，t=0.1时，准确解、向前、向后、C-N及DFF格式的比较 ... ... 图9 Du Fort-Frankel格式的结果λ=0.5，h=0.1 图10 精确解与其它格式在t=0.1的结果比较，其中λ=0.1，h=0.1 图2 向前差分的结果，λ=0.1，h=0.1 图1 精确解的结果，λ=0.5，h=0.1 图3 向前差分的结果，λ=0.5，h=0.1 图4 向后差分的结果，λ=0.1，h=0.1 图5 向后差分的结果，λ=0.5，h=0.1 图11 精确解与其它格式在t=0.1的结果比较，其中λ=0.5，h=0.1 图6 Crank-Nicolson格式的结果λ=0.1，h=0.1 图7 Crank-Nicolson格式的结果λ=0.5，h=0.1 图8 Du Fort-Frankel格式的结果λ=0.1，h=0.1