标题:Mixed volume element combined with characteristic mixed finite volume element method for oil–water two phase displacement problem
作者:Yuan, Yirang ;Sun, Tongjun ;Li, Changfeng ;Liu, Yunxin ;Yang, Qing
作者机构:[Yuan, Yirang ;Sun, Tongjun ;Liu, Yunxin ] Institute of Mathematics, Shandong University, Jinan; 250100, China;[Yang, Qing ] School of Mathematical Sc 更多
通讯作者:Yuan, Yirang
来源:Journal of Computational and Applied Mathematics
出版年:2018
卷:340
页码:404-419
DOI:10.1016/j.cam.2018.02.038
摘要:As the basic of numerical simulation of energy science, the displacement of three-dimensional oil–water two phase in porous media is discussed in this paper. For incompressible miscible displacement, the pressure is described by a flow equation in an elliptic mathematical formulation and the saturation is defined by a convection–diffusion equation. The pressure exists in the saturation equation by Darcy velocity, and controls the whole flow. We develop and improve substantially the work of Arbogast and Wheeler, then put forward a mixed volume element combined with characteristic mixed finite volume element method for two-phase displacement problem. The flow equation is discretized by the conservative mixed volume element method, which could improve the accuracy of an order for computing Darcy velocity. The saturation equation is solved by a characteristic mixed volume element method, where mixed volume element method is used to compute the diffusion term and the method of characteristics is adopted to discretize the convection term. The method of characteristics can confirm high stability of numerical simulation at the fronts, avoid numerical dispersion and nonphysical oscillation, and can adopt large time steps, obtain smaller time truncation error and improve the computation accuracy. Mixed volume element can solve the convection term, confirm conservation of mass at each element and approximate the saturation and the adjoint vector meanwhile. It is most important in numerical computation of seepage mechanics. Optimal second order estimates in L2norm are derived by theory and special techniques of priori estimates. Finally, numerical experiments are shown to illustrate the efficiency and practicability and to solve the international problem successfully.
© 2018 Elsevier B.V.
收录类别:EI
资源类型:期刊论文
TOP