Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng-Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young-Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young-Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young-Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests are carried out to verify the effectiveness of the proposed multi-scale method.

[1]  Shuyu Sun,et al.  Compositional modeling of three‐phase flow with gravity using higher‐order finite element methods , 2011 .

[2]  W. M. Haynes,et al.  Equations for the viscosity and thermal conductivity coefficients of methane , 1975 .

[3]  M. Wheeler,et al.  Anisotropic and dynamic mesh adaptation for discontinuous Galerkin methods applied to reactive transport , 2006 .

[4]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[5]  Shuyu Sun,et al.  Efficient numerical methods for simulating surface tension of multi-component mixtures with the gradient theory of fluid interfaces☆ , 2015 .

[6]  Tomoaki Kunugi,et al.  Multi-scale modeling of the gas-liquid interface based on mathematical and thermodynamic approaches , 2010 .

[7]  Dong Liang,et al.  An Approximation to Miscible Fluid Flows in Porous Media with Point Sources and Sinks by an Eulerian-Lagrangian Localized Adjoint Method and Mixed Finite Element Methods , 2000, SIAM J. Sci. Comput..

[8]  E. M. Blokhuis,et al.  Thermodynamic expressions for the Tolman length. , 2006, The Journal of chemical physics.

[9]  Yalchin Efendiev,et al.  Multiscale Finite Element Methods: Theory and Applications , 2009 .

[10]  Karen Schou Pedersen,et al.  Phase Behavior of Petroleum Reservoir Fluids , 2006 .

[11]  Shuyu Sun,et al.  An energy stable evolution method for simulating two-phase equilibria of multi-component fluids at constant moles, volume and temperature , 2016, Computational Geosciences.

[12]  Gretar Tryggvason,et al.  Direct Numerical Simulations of Gas–Liquid Multiphase Flows: Preface , 2011 .

[13]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[14]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[15]  D. Peng,et al.  A New Two-Constant Equation of State , 1976 .

[16]  U. Ozen,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[17]  Abbas Firoozabadi,et al.  Higher-order compositional modeling of three-phase flow in 3D fractured porous media based on cross-flow equilibrium , 2013, J. Comput. Phys..

[18]  JISHENG KOU,et al.  Numerical Methods for a Multicomponent Two-Phase Interface Model with Geometric Mean Influence Parameters , 2015, SIAM J. Sci. Comput..

[19]  Gergina Pencheva,et al.  A Global Jacobian Method for Mortar Discretizations of a Fully Implicit Two-Phase Flow Model , 2014, Multiscale Model. Simul..

[20]  Hongxing Rui,et al.  A mixed element method for Darcy–Forchheimer incompressible miscible displacement problem , 2013 .

[21]  J. Glimm,et al.  A multiscale front tracking method for compressible free surface flows , 2007 .

[22]  Dong Liang,et al.  An efficient S-DDM iterative approach for compressible contamination fluid flows in porous media , 2010, J. Comput. Phys..

[23]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[24]  A. Graciaa,et al.  Modeling of the Surface Tension of Multicomponent Mixtures with the Gradient Theory of Fluid Interfaces , 2005 .

[25]  Mary F. Wheeler,et al.  Discontinuous Galerkin methods for simulating bioreactive transport of viruses in porous media , 2007 .

[26]  Mary F. Wheeler,et al.  Coupling Discontinuous Galerkin and Mixed Finite Element Discretizations using Mortar Finite Elements , 2008, SIAM J. Numer. Anal..

[27]  Hussein Hoteit,et al.  Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures , 2008 .

[28]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[29]  Shuyu Sun,et al.  Multiscale Discontinuous Galerkin and Operator-Splitting Methods for Modeling Subsurface Flow and Transport , 2008 .

[30]  Shuyu Sun,et al.  Upwind discontinuous Galerkin methods with mass conservation of both phases for incompressible two‐phase flow in porous media , 2014 .

[31]  Shuyu Sun,et al.  Two-Phase Fluid Simulation Using a Diffuse Interface Model with Peng-Robinson Equation of State , 2014, SIAM J. Sci. Comput..

[32]  Hussein Hoteit,et al.  Modeling diffusion and gas–oil mass transfer in fractured reservoirs , 2013 .

[33]  Jirí Mikyska,et al.  Compositional modeling in porous media using constant volume flash and flux computation without the need for phase identification , 2014, J. Comput. Phys..

[34]  Shuyu Sun,et al.  Unconditionally stable methods for simulating multi-component two-phase interface models with Peng-Robinson equation of state and various boundary conditions , 2016, J. Comput. Appl. Math..

[35]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy and Free Energy of a Nonuniform System. III. Nucleation in a Two‐Component Incompressible Fluid , 2013 .

[36]  Gretar Tryggvason,et al.  Computational Methods for Multiphase Flow: Frontmatter , 2007 .

[37]  Mary F. Wheeler,et al.  Compatible algorithms for coupled flow and transport , 2004 .

[38]  Daniel A. Cogswell.  A phase-field study of ternary multiphase microstructures , 2010 .

[39]  C. Miqueu,et al.  Modelling of the surface tension of binary and ternary mixtures with the gradient theory of fluid interfaces , 2004 .

[40]  Nikolaus A. Adams,et al.  Scale separation for multi-scale modeling of free-surface and two-phase flows with the conservative sharp interface method , 2015, J. Comput. Phys..

[41]  Shuyu Sun,et al.  An adaptive finite element method for simulating surface tension with the gradient theory of fluid interfaces , 2014, J. Comput. Appl. Math..

[42]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[43]  Abbas Firoozabadi,et al.  Thermodynamics of Hydrocarbon Reservoirs , 1999 .

[44]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[45]  Shuyu Sun,et al.  Coupled Generalized Nonlinear Stokes Flow with Flow through a Porous Medium , 2009, SIAM J. Numer. Anal..

[46]  Mary F. Wheeler,et al.  Symmetric and Nonsymmetric Discontinuous Galerkin Methods for Reactive Transport in Porous Media , 2005, SIAM J. Numer. Anal..

[47]  Ng Niels Deen,et al.  Multi-scale modeling of dispersed gas-liquid two-phase flow , 2004 .

[48]  Yuanle Ma,et al.  Computational methods for multiphase flows in porous media , 2007, Math. Comput..

[49]  Multiscale Finite Element Methods for Elliptic Equations , 2010 .

[50]  A. Firoozabadi,et al.  Unified model for nonideal multicomponent molecular diffusion coefficients , 2007 .