叼一Il求解泊松方程的紧致高阶差分方法-西北大学学报编辑部.PDF

维普资讯 1996年 4月 西北大学学报 (自然科学版) ADr．1996 第 26卷 第 2期 JournalofNorthwestUniversity(NaturalScienceEdition) Vo1．26 No 2 ⑤ ＼叼一Il 求解泊松方程的紧致高阶差分方法 ‘ (二)Z ．孑2 田 振 夫 (宁夏大学应用数学研究所，?50021，银川；32岁，男，讲师) 摘 要 基于Hermite插值法的基本思想．提出7求解二维泊松 (Poisson)-b-程的紧致高阶差 分方法 ，得到 了一般 形式 的四阶和 六阶差分 紧致 格式 。通过敷值实验证 明了格式的 良好性 志 。 关键词 塑竺查兰；有限遘坌查塞；高趟生；鉴塾垫查 分类号 0241．3 泊松 (Poisson)方程是最为重要的一类椭 圆型偏微分方程，它是诸如流体动力学、弹性力学、传热传 质学、水动力学、电磁学等基础理论研究的基本控制方程之一，也是数值网格生成技术0所遵循的基本 方程。因此，研究其数值求解方法具有重要意义 ：2 ：： ，． ，， 本文所考虑的二维泊松方程为 L“一 + 一f(x，)， (1) o 0 式中f(x，)可视为已知函数，(，y)是待求函数。为保证文中推求格式的完整性 ，不失一般性，假定函 数 ，)和 f(x，)均为充分光滑函数。 数值求解方程(1)最普通的差分方法乃 5点中心差分格式 (CDS)。当求解区域中差分阿格较粗时， 具有二阶精度的CDS所获得的数值结果往往是不精确的。换言之，若使所求数值解具有较高的可信度， 则必须对网格进行精细划分，然而 以此增加计算量和昂贵的费用为代价，是得不偿失的 高精度差分方 法是解决CDS精度低的有效途径之一。从理论上讲，只要有足够的网格点数，总可以实现任意高精度 差分格式的构造，因此，寻求尽可能少的同格点数，以获得尽可能高精度的差分方法是人们所追求的 目 的 。采用拼凑法[2给 出了正方形网格下LapLance方程9结点六阶精度差分格式，文献3则引用解析 特解法导出了与文献 2相 同的结果。尽管实现了求解泊松方程的高精度差分格式 ]。但是所得格式关 于非齐次项f(x，)的约束条件较苛刻。其中Birkhoff“的四阶与六阶精度格式 ，要求f(x，)分别满足 二阶、四阶可导，而Manohar 的方法虽然在四阶精度的格式上放宽了文献4中对f(x，)的约束条件， 但其六阶差分格式仍未从根本上摆脱 f(x，)必须满足四阶可导的约束条件 。 本文基于Hermite插值法的基本思想n，将待定系数与原方程 (1)有机地结合起来，通过在待定逼 近式中增加导数项的技术．提出了数值求解方程 (1)的一般形式的高阶差分紧致格式 ，所得格式不仅将 文献 2～5中的格式包括于其中，而且无论从理论分析上还是从数值实验结果中，均显示出了低代价与 高精度的特点，同时函数 f(x，)的约束条件也较文献 4．5中的宽松。 1 高阶紧致差分方法 方程 (1)在点 ，)的有限近似值 ．可在该点及相邻点的函数值求得，但为了提高这种近似解的准 · 宁夏 自治区教委科研基金资助课题 收稿 日期 ：1995—03—21 维普资讯 110 一 西北大学学报 (自然科学版)