多孔介质二相驱动问题一种修正特征有限元方法的迭代解法及其收敛性分析

0 引言 多孔介质二相驱动问题的数学模型是由压力方程与浓度方程组成的偏微分方程组的初边值问题.关于该问题的数值解问题,已有大量的文献.为了得到最优的L~2-模误差估计,好多方法用混合元方法解压力方程.我们知道,混合元法得到的方程组系数矩阵是非正定的,从而解混合元比解标准元要困难得多,虽然许多人研究了混合元方法的求解问题,但到目前为止,还没有看到令人满意的好的算法.为了避开对混合元的求解,著名学者T.F.Russell考虑了用标准有限元方法解压力方程,用特征有限元方法解浓度方程的求解方法及其迭代解法,对只有分子扩散的二相驱动问题得到了最优的L~2模误差估计,对有机械弥散的一般二相驱动问题得不到最优的L~2模误差估计,同时在收敛性证明中要求压力有限元空间的指数至少是二.

AN ITERATIVE METHOD AND ITS CONVERGENT ANALYSIS OF A NEW MODIFIED CHARACTERISTIC FINITE ELEMENT METHOD FOR MISCIBLE DISPLACEMENT IN POROUS MEDIA

In this paper, an iterative method of a new modified characteristic finite element method is proposed for incompresible miscible two phase displacement in porous media. In this method, standard finite element method is used to solve the presure equation, after postprocessing, a new modified characteristic method is used to solve the concentration equation. The preconditioned conjugate gradient iterative method is used to solve discrete equations, the optimal error estimate of the iterative method, which could not be obtained by Russell1], is obtained and further piecewise linear finite element space can be used to solve pressur.