Mixed transform finite element method for solving the non-linear equation for flow in variably saturated porous media

R. G. Baca, J. N. Chung, D. J. Mulla

Research output: Contribution to journalArticlepeer-review

24 Scopus citations

Abstract

A new computational method is developed for numerical solution of the Richards equation for flow in variably saturated porous media. The new method, referred to as the mixed transform finite element method, employs the mixed formulation of the Richards equation but expressed in terms of a partitioned transform. An iterative finite element algorithm is derived using a Newton-Galerkin weak statement. Specific advantages of the new method are demonstrated with applications to a set of one-dimensional test problems. Comparisons with the modified Picard method show that the new method produces more robust solutions for a broad range of soil-moisture regimes, including flow in desiccated soils, in heterogeneous media and in layered soils with formation of perched water zones. In addition, the mixed transform finite element method is shown to converge faster than the modified Picard method in a number of cases and to accurately represent pressure head and moisture content profiles with very steep fronts.

Original languageEnglish (US)
Pages (from-to)441-455
Number of pages15
JournalInternational Journal for Numerical Methods in Fluids
Volume24
Issue number5
DOIs
StatePublished - Mar 15 1997

Keywords

  • Finite elements
  • Porous flow
  • Richards equation

Fingerprint

Dive into the research topics of 'Mixed transform finite element method for solving the non-linear equation for flow in variably saturated porous media'. Together they form a unique fingerprint.

Cite this