Crank_Nicolson差分格式及其稳定性研究.doc

C ran k - N ico lso n 差分格式 及其稳定性研究 李 华 周维奎 (成都理工学院, 成都 610059) 邓培智 (核工业部西南物理研究院, 成都 610041)  【摘 要】 本文以自己独特的方式, 构造了一维和二维抛物型方程的 C rank - N ico lso n 差分格式。本文不仅详细地给出了离散误差的表达式, 而且论证了它们的稳定性。该差分格式具有精度高, 稳定性好, 计算量和 存储量都比较小的特点, 是一个很理想, 便于应用的差分格式。 【关键词】 抛物型方程 隐格式 离散误差 绝对稳定性 中图法分类号: P 628 在计算机上用有限差分法求解抛物型方程初边值问题, 已有多种数值解法。这些数值解法 的稳定性、离散误差的估计、计算量和存贮量, 都是人们普遍关注的问题。本文研究的 C ran k - N ico lso n 差分格式, 根据论证是绝对稳定的差分格式, 它的离散误差、计算量和存贮量都比较 小, 是稳定性好, 精确度高, 实用性强的差分格式之一。 1 建立差分格式 考虑一维抛物型方程初边值模型问题:  5u 5t =  5 2 u u (x , 0) = f (x ) ,0≤x ≤au (0, t) = g1 ( t) ,u (a , t) = g2 ( t).t 0 5x 2 , 0 x a , t 0  (1)  其中 t 为时间变量, x 为空间变量, f (x ) , g1 ( t) 和 g2 ( t) 为已知函数。 首先, 对求解区域进行网格剖分。设时间步长为 Σ, 空间步长为 h = a , 得网格节点 (x j, tk ) , N , N , k = 0, 1, 2, 。记 tk + 1 = (h + 1 ) Σ, 网比 r= Σ。在点 (x j, tk + 1 ) 用一阶中心差商 j= 0, 1, 2, h 2 2 2 2 逼近 5u?5t, 用加权二阶中心差商逼近 5 2 u ?5x2 , 即得如下加权六点差分格式1 : uk + 1 u j = r [ Η(u j+ 1 - 2u j + u j- 1 ) + k k + 1 k + 1 k + 1 (1- Η) (u j+ 1 - 2u j + u j- 1 ) ] k k k (2) - j 其中 Η为权因子, 0≤Η≤1。 j= 1, 2, , N - 1 k = 0, 1, 2, 由(2) 式不难看出: 1. Η= 0 时, (2) 为古典显格式: uk + 1 = u k + r (uk 1 - 2uk + uk 1 ) j j j+ j j- 2. Η= 1 时, (2) 为古典隐格式: - ruk + 1 + (1+ 2 r) uk + 1 - ruk + 1 = u k j+ 1 j 3. Η= 1?2 时, (2) 为 C ran k - N ico lso n 格式: j- 1 j - r u k + 1  k + 1  r  k + 1  r k k r k 2 j+ 1 + (1+ r) u j - 2 u j- 1 = 2 u j+ 1 + (1- r) u j + 2 u j- 1 (3) C ran k - N ico lso n 格式是一个隐格式。现在, 我们给出用它来求解初边值问题 (1) 的解法。 根据问题 (1) 的初始条件和边值条件, 有 u 0 j = f (x j) , j= 1, 2, , N - 1  u k k (4) 0 = g1 ( tk ) , uN = g2 ( tk ). k = 0, 1, 2, 显然, 若已求得第 k 层的解 u k , 则差分方程 ( 3) 和初边值条件 ( 4) 给出关于求 uk + 1 的线性 j j 代数方程组。 该方程组是三对角线性方程组, 系数矩阵是一个严格对角占优的三对角矩阵, 故 可以用追赶法求解。每计算一层, 就需要用追赶法求解一次线性方程组, 逐层计算, 就可以求得 初边值问题 (1) 的数值解。 这个解题过程, 实际计算表明, 效果是很好的。  k + 1 k + 1  k + 1  + h 2 )  2离散误差的估计 假设问题 (1) 的解 u = u (x , t) 充分光滑, 利用泰劳展式, 在(x j, tk + 1 ) 有 2 + O (Σ2 1 [ u (x , t - u x , t = 5u ]k +