热传导方程初边值问题紧致格式.pdfVIP

北京建筑工程学院学报 JOURNAL OF BEIJING INSTITUTE OF CIVIL ENGINEERING AND ARCHITECTURE -Vo!．10 No．1 1994 第l0卷第 l期 胡 ／ 热传导方程初边值问题的紧致格式 苎 ! B．Gu“ a-， 噬 出 0 ． 【摘 要】在此文内，我~fl给 出了热传导方程紧致格式的半离散化能量估计。 然后，用统一 的办法分析 了一族奎离散格式的稳定性．在此基础上 ．我ell讨论 了紧 致格式 的收敛速度。最后将我们的结果推广到二维井指 出在高维情形有 同样 的 结 果 。 ㈣言 专黻暴赫式 徊进_拟孳 众所塌知．国民经济和国防建设中的许多问题都归结为偏微分方程初边值问题的数值求 解 ．因此如何构造可靠而又精确的数值方法去求解它们，自然成为人们关注的一十焦点，自 二次大战以来发展起来的一系列二阶方法，成动地解决了许多问题。但是随着A们对客观世 界认识的深入，也 日益暴露出二阶方法的局限．这主要表现在精度不够理想 要提高二阶方 法的精度须增加大量的网格结点．相应地增加计算 的存贮量和计算量．由于要解决的问题 日益复杂化，大型化，尽管现代计算机已经有了飞跃的发展，但在存贮和运算速度上还 能 荫足需要，即使在将来．解决这个矛盾的前景 电 ；容乐观 这就是为什 幺人们开始研究高阶 格式的主要原 因。 一 般的四阶方法，对一维问题而 占，在空间方 向上一般要涉及五点，迭无疑给计算和边 界条静的处理增加很多困难，但在文 (1]内，Kreiss提出了一个四阶精度的紧致差分格式， 在空间方向上只涉及到三点，显然这具有葜大好处。迄夸为止，已有好多文章报导了用紧致 格式计算的大量复杂 问题，结论是 Krciss所 倡 导 的紧致格式是一个非常有效的方法，正有 鉴于此，关于紧致格式的研究方 未艾，源流不息． 我们将针对热传导方程的第一及第二边值 问题： 害 ()∈(0，1)×(0 (0，t)：0，扪 ，)：0 (O·1) ( O)= () 第1期 黄铎等：热侍旱方程初迎值问题的紧盘格式 = ( })∈ (o，1)x(o，T) d … ～ 0， (1，f)：0 (。一) ： - a 1 O l (，O)= () 首先在第一节 内给出紧致格式的半离散化能量估计．在第二节 内，我们使用一种 统 一 的 方 式，对一族垒离散的紧致格式的稳定性进行了分析，这远较逐一分析它们为简单明了．关于 紧致格式收敛速度的估计则在第三节内给出．在第四节 内，我们指出许多重要的差分格式 ， 其稳定性分析 已包括在我们第二节的主要结果之 内．关于将紧致格式应用高维问题的稳定性 和收敛性分析则在第五带内给出．由于在半离散化能量估计、稳定性及收敛性分析中计及了 边界条件且使用了矩阵分析方法，因而所得结果不可再改进 。 1 紧致格式的半离散化能量估计 在本节 内，我们将就热传导方程韧边值 问题 (01)和 (0．2)的半离散化紧致格式逼近I 1