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

飞101. 21 No. 1 嬉台大学学报{自然科学与工程摄) 第21 卷第 1 期 2008 年 1 月 Jo泣rnal of Yantai University (Natural Science and Engineering Edition) Jan. 2∞8 文章编号:1∞4-8820(200的01兮∞H)5 一维不可压缩两相渗流驱动问题的有限体积元方法 封常荣，稳传军 (烟台大学数学与信息科学学院，山东施台 264∞5 ) 摘 要:在初始玛格剖分上采取分手支线性函数空间件为有限体和、元方法的试探函数空 间，在相应的对偶网格上采取分技常数函数空间作为其栓验函数空窍，对一维不可压缩两， 2 相渗流驱动问题提出了全离散有限体和、元方法，并得到 L _ 模误差估计. 关键词:有限体积元方法:不可压缩;全离散:线性元;误差估计 中图分类号: 0241.82 文献标识码:A Lipschitz 条件. 将有限体积元方法应黑到一维不可压缩两梧 渗流驱动问题，文献[1 - 3 ]曾对二维两相渗琉驱 (2)q 光滑分布. 动问题做了一系列分析与研究，本文对。这一这 (3)0 :::三 K. ~ k(x , s) 王三 K , 0 :::三 D. 主三 域截靡，取横截面1 = [α ， b ]， 在这一横截面上考 D(x) 运 D* ， 0:::三 φ* 主三 φ(X) 运 φ拿. 察两梧渗混驱动闰题，使数值模拟更加精细，模 其中:儿 ， K* ， D. ， D事， φ. ， φ* 表示正的常数. 型可串下002 个方程藕合而成的方程组来表 在以下的分析中，用C表示一般的正常数，不同之 示[IJ. 处含义不同. 字= - ~ (k怡， s) 学) = q , dX dX dX 1 全离散有黑体积元格式 X E [α ， bJ ， 0 t 运 T， (1) 用11 表示1 = [α ， bJ 上的r I资Sobolev 空间， φ(劝学 +U 孚- ~ (D 孚) = 2 时r dX dX I1X 卜，寸表示L (1) 上的内积. 11. Il r 表示11 的起 (S-S)q ， XE[α ， b]， Ot~T. (2) 1 数.又 H~ (I) = 1u ε H (1) ;u( α) =叭的=