Click on “Download PDF” for the PDF version or on the title for the HTML version.


If you are not an ASABE member or if your employer has not arranged for access to the full-text, Click here for options.

Application of Mixed Finite Element Method to One-Dimensional Unsaturated Flow

Published by the American Society of Agricultural and Biological Engineers, St. Joseph, Michigan www.asabe.org

Citation:  Paper number  032112,  2003 ASAE Annual Meeting . (doi: 10.13031/2013.14005) @2003
Authors:   Debasmita Misra, John L. Nieber
Keywords:   Unsaturated flow, Darcy flux, Richards equation, Mixed finite element method

Flow in the unsaturated porous media is governed by a strongly nonlinear diffusion type equation with a nonlinear, forced convection term due to gravity. These features make it particularly difficult to solve by any means. The flow induced heterogeneity and the specific moisture capacity of the medium make the governing equation of flow highly non-linear due to the dependence of these parameters on the pressure and the moisture content distribution in the medium.

The velocity at which water flows through the unsaturated zone is at all times less than or equal to that in a porous medium saturated with the fluid. As a result, the chemicals transported by the water find more resident time in the medium for reactions, transformations and attenuations. This could have significant impact on the quality of water reaching the ground water table or on the chemicals absorbed by the plants or escaping to the atmosphere in volatile form.

Accurate simulation of Darcy flux is essential to simulate contaminant transport in the unsaturated zone accurately. The mixed finite element method has been used to obtain highly accurate flux distribution in unsaturated zone flow applications. The objective of this paper is to present a simplified conceptual description of the mixed finite element method as applied to unsaturated flow modeling, to compare the solutions obtained from the method to those of the conventional finite element methods and to analyze any special properties of the solutions obtained from the mixed finite element methods. It has been shown in this paper that the solutions obtained from the mixed finite element method are highly accurate in rapidly changing flux distributions and heterogeneous flow distributions even when relatively coarse grids are used to obtain the solution.

(Download PDF)    (Export to EndNotes)