A Locally Conservative Eulerian-Lagrangian Method for a Model Two-Phase Flow Problem in a One-Dimensional Porous Medium

Motivated by possible generalizations to more complex multiphase multicomponent systems in higher dimensions, we develop an Eulerian--Lagrangian numerical approximation for a system of two conservation laws in one space dimension modeling a simplified two-phase flow problem in a porous medium. The method is based on following tracelines, so it is stable independent of any CFL constraint. The main difficulty is that it is not possible to follow individual tracelines independently. We approximate tracing along the tracelines by using local mass conservation principles and self-consistency. The two-phase flow problem is governed by a system of equations representing mass conservation of each phase, so there are two local mass conservation principles. Our numerical method respects both of these conservation principles over the computational mesh (i.e., locally), and so is a fully conservative traceline method. We present numerical results that demonstrate the ability of the method to handle problems with shoc...

[1]  T. F. Russell,et al.  Convergence analysis of an approximation of miscible displacement in porous media by mixed finite elements and a modified method of characteristics , 1984 .

[2]  Hong Wang,et al.  A numerical modeling of multicomponent compressible flows in porous media with multiple wells by an Eulerian-Lagrangian method , 2005 .

[3]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[4]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[5]  Dong Liang,et al.  An ELLAM Approximation for Highly Compressible Multicomponent Flows in Porous Media , 2002 .

[6]  Todd Arbogast,et al.  A Characteristic-Mixed Method for Contaminant Transport and Miscible Displacement , 1992 .

[7]  Hong Wang,et al.  A Comparison of ELLAM with ENO/WENO Schemes for Linear Transport Equations , 2000 .

[8]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[9]  Todd Arbogast,et al.  A Fully Conservative Eulerian-Lagrangian Stream-Tube Method for Advection-Diffusion Problems , 2012, SIAM J. Sci. Comput..

[10]  A. Tits,et al.  User's Guide for FSQP Version 2.0 A Fortran Code for Solving Optimization Problems, Possibly Minimax, with General Inequality Constraints and Linear Equality Constraints, Generating Feasible Iterates , 1990 .

[11]  Shirley Dex,et al.  JR 旅客販売総合システム（マルス）における運用及び管理について , 1991 .

[12]  Béatrice Rivière,et al.  Computational methods for multiphase flows in porous media , 2007, Math. Comput..

[13]  Mary F. Wheeler,et al.  Some improved error estimates for the modified method of characteristics , 1989 .

[14]  Chi-Wang Shu,et al.  Conservative high order semi-Lagrangian finite difference WENO methods for advection in incompressible flow , 2011, J. Comput. Phys..

[15]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[16]  D. W. Peaceman.  Fundamentals of numerical reservoir simulation , 1977 .

[17]  Todd Arbogast,et al.  A Fully Mass and Volume Conserving Implementation of a Characteristic Method for Transport Problems , 2006, SIAM J. Sci. Comput..

[18]  J. Douglas.  The Convergence of a Locally Conservative Eulerian-Lagrangian Finite Difference Method for a Semilin , 2000 .

[19]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

[20]  Hong Wang,et al.  A locally conservative Eulerian‐Lagrangian control‐volume method for transient advection‐diffusion equations , 2006 .

[21]  Cass T. Miller,et al.  An ELLAM approximation for advective-dispersive transport with nonlinear sorption , 2006 .

[22]  A. Tits,et al.  Nonlinear Equality Constraints in Feasible Sequential Quadratic Programming , 1996 .

[23]  Todd Arbogast,et al.  An Eulerian-Lagrangian WENO finite volume scheme for advection problems , 2012, J. Comput. Phys..

[24]  O. Pironneau.  On the transport-diffusion algorithm and its applications to the Navier-Stokes equations , 1982 .

[25]  Jianxian Qiu,et al.  On the construction, comparison, and local characteristic decomposition for high-Order central WENO schemes , 2002 .

[26]  M. Wheeler,et al.  A characteristics-mixed finite element method for advection-dominated transport problems , 1995 .

[27]  Richard E. Ewing,et al.  Eulerian-Lagrangian Localized Adjoint Methods for a Nonlinear Advection-Diffusion Equation , 1994 .

[28]  C. D. Chambers.  On the Construction of οὐ μή , 1897, The Classical Review.

[29]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[30]  Felipe Pereira,et al.  A locally conservative Eulerian–Lagrangian numerical method and its application to nonlinear transport in porous media , 2000 .