Premium
A mass conservative numerical solution for two‐phase flow in porous media with application to unsaturated flow
Author(s) -
Celia Michael A.,
Binning Philip
Publication year - 1992
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/92wr01488
Subject(s) - linearization , porous medium , richards equation , flow (mathematics) , mechanics , advection , finite element method , numerical analysis , mathematics , two phase flow , porosity , soil water , geotechnical engineering , mathematical analysis , nonlinear system , physics , geology , thermodynamics , soil science , quantum mechanics
A numerical algorithm for simulation of two‐phase flow in porous media is presented. The algorithm is based on a modified Picard linearization of the governing equations of flow, coupled with a lumped finite element approximation in space and dynamic time step control. Numerical results indicate that the algorithm produces solutions that are essentially mass conservative and oscillation free, even in the presence of steep infiltrating fronts. When the algorithm is applied to the case of air and water flow in unsaturated soils, numerical results confirm the conditions under which Richards's equation is valid. Numerical results also demonstrate the potential importance of air phase advection when considering contaminant transport in unsaturated soils. Comparison to several other numerical algorithms shows that the modified Picard approach offers robust, mass conservative solutions to the general equations that describe two‐phase flow in porous media.