A Mass Conservative Numerical Solution for Two-Phase Flow in Porous Media With Application to Unsaturated Flow

Publication: Research - peer-reviewJournal article – Annual report year: 1992

View graph of relations

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.
Original languageEnglish
JournalWater Resources Research
Issue number10
Pages (from-to)2819-2828
StatePublished - 1992
Download as:
Download as PDF
Select render style:
Download as HTML
Select render style:
Download as Word
Select render style:

ID: 6254232