一维变系数对流扩散方程第三边值问题的紧有限体积方法.pdfVIP

第31卷 第6期 工 程 数 学 学 报 Vo1．31No．6 2014年12月 CHINESEJOURNALOFENGINEERING MATHEMATICS Dec．2014 doi：10．3969／j．issn．1005—3085．2014．06．011 文章编号：1005—3085(2014)06—0889—14 一 维变系数对流扩散方程第三边值问题的 紧有限体积方法木 陈宏霞，王 同科 (天津师范大学数学科学学院，天津 3003871 摘 要：对流扩散方程在工程计算中具有广泛应用．本文研究一维变系数对流扩散方程第三边 值问题的高精度有限体积方法．通过在控制体积上积分导出了方程的积分守恒形式， 然后对积分守恒形式利用泰勒公式和二次埃尔米特插值进行离散得到了紧有限体积格 式．该格式导出的线性代数方程组具有三对角性质，因此可使用追赶法求解．进而， 通过分析截断误差，采用能量方法证明了格式按照几种标准的离散范数四阶收敛．最 后，数值算例验证了格式的正确性和有效性，这与理论分析结果是一致的． 关键词：对流扩散方程第三边值问题；紧有限体积方法；误差估计；四阶精度 分类号：AMS(2000)65N08；65N12 中图分类号：O241．82 文献标识码：A 1 引言 对流扩散方程出现在应用数学的许多分支，其数值算法研究一直非常活跃[1-9]．对流 扩散方程的数值求解方法有多种，其 中包括有限差分方法[卜]，有限体积方法和一些其 他方法8【19】．文献 [2]利用Richardson外推方法和插值技巧将定常型对流扩散方程的紧差 分格式离散精度从4阶提高到了6阶；文献f51对一类常系数对流扩散方程进行转化，给 出了一种可以达到任意阶精度的三点紧致差分格式；文献f61将对流项采用四五阶组合迎 风格式离散，得到了一种求解非定常对流扩散方程的高精度组合紧致有限差分格式；文 献 [8]给出一种新的移动有限点无网格方法，具有2阶收敛率；文献f91提出了一种松弛方 法；文献 [101给出了一类反应扩散方程的二阶有限体积格式．最近，文献f11，121使用有限 体积方法离散具有第三边界条件的扩散方程，得到一类具有三对角性质的四阶有限体积格 式，称为紧有限体积格式．借鉴文献f12]中对二阶导数项的离散思想，本文对于一维变系 数对流扩散方程第三边值问题的扩散项使用同样的方法离散，构造了一类具有四阶精度的 有限体积格式．该格式形成的线性代数方程组具有三对角性质，可用追赶法求解 ．在收敛 性分析部分，文献 1【2]没有达到最优阶，本文使用新的方法证明了格式具有4阶精度，做 到了理论与计算的统一． 收稿 日期：2013—05—15．作者简介：陈宏霞 (1988／f12~生)，女，硕士．研究方向：偏微分方程数值解 基金项目：国家自然科学基金． 890 工 程 数 学 学 报 第31卷 本文第2节针对一维变系数对流扩散方程第三边值问题给出了紧有限体积格式的构 造方法，第 3节给出了该格式的截断误差分析，利用能量方法证明了格式按照离散 。范 数，H 半范数和最大范数达到了4阶精度 ．最后，通过具体算例验证 了格式的高精度 ， 并说明了格式具有较大的实用性． 2 格式构造 考虑区间I=fa，6J上的一维定常型对流扩散方程第三边值问题 一 rf上+p()u+q(X)U=f()，X∈(r上，b)， (1) 一 (口)+O~lU(a)=91， f“(b)+~2u(b)=g2， (2) 其中p()，g()，f(x)适当光滑， ，2为常数，且满足限制条件： (a)