求解椭圆型方程界面问题修正有限体积方法？.pdfVIP

第32卷 第2期 工 程 数 学 学 报 Vo1．32No．2 2015@z04~ CHINESE JOURNAL OFENGINEERING MATHEMATICS Apr·2015 doi：10．3969／j．issn．1005—3085．2015．02．006 文章编号：1005—3085(2015)02—0213—13 求解椭圆型方程界面问题的修正有限体积方法木 续小磊， 冯秀芳 (宁夏大学数学与计算机学院，银川 750021) 摘 要：本文基于经典的有限体积方法，讨论了带有间断系数二维椭圆型方程的求解问题．文 中通过在求解椭圆型方程时截取通量函数的更多项泰勒展开式，同时结合有限差分方 法对与边界相邻的网格点进行特殊处理，改进了间断系数的求解方法，得到了一种修 正的有限体积方法，该方法在界面处是一阶精度，但方法整体可达到二阶精度．数值 实验表明，文中构造的修正有限体积方法是有效的． 关键词：椭圆型方程；有限体积方法；界面问题；调和平均法 分类号：AMS(2000)65N10；65F30 中图分类号：O241．82 文献标识码：A 1 引言 实际生活中很多现象可以由带有间断系数的椭圆型方程来描述，如在材料科学中，不 同传导率的两种材料的拼接 ．在渗流力学中，复杂地质结构或多相流体导致具有间断渗透 率或扩散系数的不可混溶驱动 问题等 ．由界面问题所导 出的偏微分方程 的解在界面处一般 是不连续的，这使得大多数常规的求解偏微分方程的数值解法往往会失效．为解决这个 问题，流体力学学者做出了很多具有重大意义的探索．1977年Peskin提出了浸入边界方 法 【卜3】用于模拟血液在人体心脏 中的流动，该方法兼具数学建模和数值离散的特点．整个 流场计算都使用笛卡尔网格，而不是按照物体形状生成复杂的贴体网格，无需处理从物理 平面到计算平面的坐标和网格转换问题，因此提高了计算效率而且节省了网格生成所需的 时间．但该方法在界面处一般只能达到一阶精度，想要达到二阶精度，则需要谨慎的选 择Delta函数．1994年，Leveque和Li发展了浸入界面方法 【，】，用于求解带有不连续系 数和奇异源项的椭圆型方程．该方法从实际问题的物理背景或带有约束条件的微分方程本 身出发，预先给出跳跃条件 ，推导出诸多界面关系式．在连续区域通常采用高精度的差分 格式，而在界面附近根据据求得的界面关系对数值方法进行修正，从而得到界面处的差分 格式．80年代，Bell和Samarskii等人提出求解椭圆界面问题的调和平均法[0-8]，调和平 均法在求解一维椭圆型界面问题时有很好的计算结果，基本能达到二阶精度，但在求解 二维 问题时不仅需要较大的工作量求解方法 中所涉及的积分项，而且在实际计算 中并不 能达到理论上的二阶精度．针对这个缺陷，很多学者在提高有限体积法精度上作出了努 收稿 日期：2013—04—28．作者简介：续小磊 (1986年5月生)，男，硕士．研究方向：偏微分方程数值解 基金项目：国家 自然科学基金11361045)；宁夏 自然科学基金 (NZ14233)． 214 工 程 数 学 学 报 第32卷 力。【一】，如Ewing等人[12]通过截取通量函数的更多项泰勒展开式以求获得更高的精度 ． 国内的宋淑红和王双虎[13,14l则发展了自适应的有限体积方法．本文在文献 1【2】的基础上 推导出二维椭圆型方程界面问题的差分格式，并结合有限差分方法对网格边界点进行特殊 处理，对扩散系数的求解方法进行改进得到一种新的差分格式． 2 问题描述 在矩形区域 内考虑带有间断系数的二维椭圆型方程 一 )~Ox2+O 、 y2]= ，)， (z，)∈[0，1]×[0j1] (1) 扩散系数 I1(z，Y)，(，Y)∈Q1，