Multidimensional self-similar analytical solutions of two-phase flow in porous media

Abstract In general, analytical solutions serve a useful purpose to obtain better insights and to verify numerical codes. For flow of two incompressible and immiscible phases in homogeneous porous media without gravity, one such method that neglects capillary pressure in the solution was first developed by Buckley and Leverett (1942). Subsequently, McWhorter and Sunada (1990) derived an exact solution for the one and two dimensional cases that factored in capillary effects. This solution used a similarity transform that allowed to reduce the governing equations into a single ordinary differential equation (ODE) that can be further integrated into an equivalent integral equation. We present a revision to McWhorter and Sunada solution by extending the self-similar solution into a general multidimensional space. Inspired by the derivation proposed by McWhorter and Sunada (1990), we integrate the resulting ODE in the third and higher dimensions into a new integral equation that can be subsequently solved iteratively by means of numerical integration. We developed implementations of the iterative schemes for one- and higher dimensional cases that can be accessed online on the authors’ website.

[1]  Seyyed A. Hosseini,et al.  Role of partial miscibility on pressure buildup due to constant rate injection of CO2 into closed and open brine aquifers , 2011 .

[2]  Effect of NAPL Source Morphology on Mass Transfer in the Vadose Zone , 2015, Ground water.

[3]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[4]  J. Parlange,et al.  n-Dimensional first integral and similarity solutions for two-phase flow , 2003, The ANZIAM Journal.

[5]  Radek Fucík,et al.  Discontinous Galerkin and Mixed-Hybrid Finite Element Approach to Two-Phase Flow in Heterogeneous Porous Media with Different Capillary Pressures , 2011, ICCS.

[6]  T. Illangasekare,et al.  Semianalytical Solution for Two‐Phase Flow in Porous Media with a Discontinuity , 2008 .

[7]  Radek Fučík,et al.  Mixed-hybrid finite element method for modelling two-phase flow in porous media (Czech Japanese Seminar in Applied Mathematics 2010 : Praha/Telc, September 2010) , 2012 .

[8]  Iuliu Sorin Pop,et al.  A New Class of Entropy Solutions of the Buckley-Leverett Equation , 2007, SIAM J. Math. Anal..

[9]  T. Bjørnarå,et al.  A pseudospectral approach to the McWhorter and Sunada equation for two-phase flow in porous media with capillary pressure , 2013, Computational Geosciences.

[10]  van Cj Hans Duijn,et al.  Similarity solution for capillary redistribution of two phases in a porous medium with a single discontinuity , 1998 .

[11]  J. Molenaar,et al.  The effect of capillary forces on immiscible two-phase flow in heterogeneous porous media , 1995 .

[12]  P Scovazzo,et al.  Modeling of two-phase flow in membranes and porous media in microgravity as applied to plant irrigation in space. , 2001, Water resources research.

[13]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[14]  Michal Beneš,et al.  An Improved Semi‐Analytical Solution for Verification of Numerical Models of Two‐Phase Flow in Porous Media , 2007 .

[15]  A. Kalra,et al.  Evaluating changes and estimating seasonal precipitation for the Colorado River Basin using a stochastic nonparametric disaggregation technique , 2011 .

[16]  S. E. Buckley,et al.  Mechanism of Fluid Displacement in Sands , 1942 .

[17]  Z. Chen,et al.  Some invariant solutions to two-phase fluid displacement problems including capillary effect , 1988 .

[18]  J. Parlange,et al.  Exact solutions to radially symmetric two-phase flow for an arbitrary diffusivity , 2005 .

[19]  David B. McWhorter,et al.  Exact integral solutions for two‐phase flow , 1990 .