一维不可压缩两相渗流驱动问题的特征有限体积元方法.pdfVIP

烟台大学学报{自然科学与工程版) Vo1. 21 No.3 第21 卷第3 期 2∞8 年7 月 Joumal of Yantai Vniversity (Natural Science and Engineering Edition) Jul. 2∞8 文章编号:1ω4-8820(2∞8)03呵。166啪 一维不可压缩两相渗流驱动问题的 特征有限体积元方法 封常荣，陈传军 (烟台大学教学与信息科学学院，山东烟台2ωω5) 摘 要:在初始网格剖分上采取分段线性函数空间作为特征有限体积元方法的试探函数 空间，在相应的对偶网格上采取分段常数函数空间作为其检验函数空间，将特征线方法 与有限体积元方法相结合，构造出特征有限体积元方法，对一维不可压缩两相渗流驱动 1 模误差估计. 问题提出了全离散有特征限体积元方法，并得到最优阶的 H 关键词:特征有限体积元方法;不可压缩;全离散;线性元;误差估计 : 0241.82 文献标识码:A 中图分类号 有限体积元方法作为求解微分方程的-类主要方法，由于其具有局部守恒性质，日益受到普遍关 注[1-4] 文中将特征线方法与有限体积元方法相结合，构造出特征有限体和、元方法，该方法综合了特征 有限差分方法和特征有限元方法的主要优点，与特征有限差分方法相比，它保持了差分方法的计算简 单性，而且具有网格剖分灵活，可以在不规则网格上计算的优点;与特征有限元方法相比，它具有计 算量小，格式直观易于计算的优点，而且可以近似达到有限元的精度[5] 文中把有特征限体和、元方法应用到一维不可压缩两相渗流驱动问题，文献[6J 曾对二维可压缩可 棍溶驱动问题做了一系列分析与研究，本文对。这一区域截断，取横截面[=[α ， b ]， 在这一横截面上 考察一维不可压缩两相海流驱动问题，使数值模拟更加精细，模型可由下面的微分方程组来表示[7] : 生= _ ~ð (k(x ， s) 学) =q , X E [a ,bJ , 0 t 运 T， (1) dX dX dX ðs φ(x) 一 +u 一-一(D 一) = (s - s)q , X E [α ， b ]， Ot 运 T， (2) ðx k(x) 其中:u 是达西速度， p 是压力， s 是饱和度， q 是产量， φ 为孔隙度. k(x , s) = 一一， k(x) 为介质渗透 μ(s) 率， μ(s) 为依赖于水的饱和度s 的蒙古度. s 为注人流体饱和度. D=D(x) 为扩散系数.给定的初始条件 为 s ( 0 , x) =so ( x) , p ( 0 , x) =PO ( x) , x ε [a ， b ]， u( α) =μ(b) = 0 ,