HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Fucík, R. Articles by Illangasekare, T. H. Search for Related Content PubMed Articles by Fucík, R. Articles by Illangasekare, T. H. GeoRef GeoRef Citation Agricola Articles by Fucík, R. Articles by Illangasekare, T. H. Related Collections Analytical Solutions Numerical Solutions Multiphase Fluid Flow

ORIGINAL RESEARCH

An Improved Semi-Analytical Solution for Verification of Numerical Models of Two-Phase Flow in Porous Media

Radek Fuík^a, Jií Mikyka^a, Michal Bene^a,* and Tissa H. Illangasekare^b

^a Department of Mathematics, Faculty of Nuclear Science and Physical Engineering, Czech Technical University in Prague, Trojanova 13, 120 00 Prague, Czech Republic
^b Center for Experimental Study of Subsurface Environmental Processes (CESEP), Division of Environmental Science and Engineering, Colorado School of Mines, 1400 Illinois Street, Golden CO 80401, CO
^* Corresponding author (benes{at}km1.fjfi.cvut.cz)

Received 13 February 2006.

A closed-form solution for one-dimensional two-phase flow through a homogeneous porous medium is presented that is applicable to water flow in the vadose zone and flow of nonaqueous phase fluids. The solution is a significant improvement to the one originally presented by McWhorter and Sunada, allowing the analysis of wetting phase entry saturations ranging from residual to full. Our aims are to provide a detailed analysis of how the solution to the governing partial differential equation of two-phase flow can be obtained from a functional integral equation arising from the analytical treatment of the problems and to present an improved algorithm for the implementation of this solution. The integral functional equation is obtained by imposing a set of assumptions for the boundary conditions. The proposed method can be used to obtain solutions that incorporate a wide range of saturation values at the entry point. The semi-analytical solution will be useful in the verification of vadose zone flow and multi-phase flow codes designed to simulate more complex two-phase flow problems in porous media where capillary effects must be included.