Semianalytical Solution to Richards' Equation for Layered Porous Media

Traditional finite-difference and finite-element solutions to Richards' equation can exhibit stability problems and mass balance errors when sharp saturation fronts migrate across material interfaces. Several semianalytical solutions to Richards' equation have been developed for layered media, but none allow for arbitrary constitutive relationships between capillary pressure, water saturation, and relative water permeability. This paper develops a more general semianalytical solution to Richards' equation that can simulate unsaturated flow in layered media and utilize arbitrary constitutive relationships. Based on a combination of the Runge-Kutta and shooting methods, initial-boundary-value problems are solved for layered systems without experiencing the types of stability problems commonly associated with numerical models. The proposed numerical scheme is computationally efficient and numerically stable, and it compares favorably with other analytical solutions. Results of two example simulations demonst...

[1]  Forman S. Acton,et al.  Numerical methods that work , 1970 .

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

[3]  R. H. Brooks,et al.  Properties of Porous Media Affecting Fluid Flow , 1966 .

[4]  H. Keller.  Numerical Methods for Two-Point Boundary-Value Problems , 1993 .

[5]  Rajesh Srivastava,et al.  Analytical solutions for one-dimensional, transient infiltration towards the water table in homogeneous and layered soils , 1991 .

[6]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[7]  Per Moldrup,et al.  RAPID AND NUMERICALLY STABLE SIMULATION OF ONE‐DIMENSIONAL, TRANSIENT WATER FLOW IN UNSATURATED, LAYERED SOILS , 1989 .

[8]  F. Dullien.  Porous Media: Fluid Transport and Pore Structure , 1979 .

[9]  W. R. Gardner.  SOME STEADY‐STATE SOLUTIONS OF THE UNSATURATED MOISTURE FLOW EQUATION WITH APPLICATION TO EVAPORATION FROM A WATER TABLE , 1958 .

[10]  John L. Nieber,et al.  Analytic Steady-State Solution to One-Dimensional Unsaturated Water Flow in Layered Soils , 1989 .

[11]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[12]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[13]  A. Corey.  Mechanics of Immiscible Fluids in Porous Media , 1986 .

[14]  Fred Marinelli,et al.  Theory of Soil‐Water Flow , 2015 .

[15]  J. Bear.  Hydraulics of Groundwater , 1979 .

[16]  William H. Press,et al.  Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing , 1992 .

[17]  R. Lenhard,et al.  Measurement and modeling of three-phase saturation-pressure hysteresis , 1992 .

[18]  P. Milly,et al.  A mass-conservative procedure for time-stepping in models of unsaturated flow , 1985 .

[19]  Keith L. Bristow,et al.  Simulating Water Movement in Layered and Gradational Soils Using the Kirchhoff Transform , 1990 .

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

[21]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[22]  P. J. Ross,et al.  Efficient numerical methods for infiltration using Richards' equation , 1990 .