A numerical approach for approximating statistical moments of hydraulic heads of variably saturated flows in multi-dimensional porous media is developed. The approximation relies on a first-order Taylor series expansion of a finite element flow model and an adjoint state numerical method for variably saturated flows to evaluate sensitivities. This approach can be employed to analyze uncertainti...

A computationally simple, numerical algorithm capable of solving a wide variety of twodimensional, variably saturated flow problems is developed. Recent advances in modeling variably saturated flow are incorporated into the algorithm. A physically based form of the general, variably saturated flow equation is solved using finite differences (centered in space, fully implicit in time) employing ...

In this paper, we compare the effectiveness of three preconditioning strategies in simulations of variably saturated flow. Using Richards’ equation as our model, we solve the nonlinear system using a Newton-Krylov method. Since Krylov solvers can stagnate, resulting in slow convergence, we investigate different strategies of preconditioning the Jacobian system. Our work uses a multigrid method ...

We consider the simulation of steady-state variably saturated groundwater flow using Richards’ equation (RE). The difficulties associated with solving RE numerically are well known. Most discretization approaches for RE lead to nonlinear systems that are large and difficult to solve. The solution of nonlinear systems for steadystate problems can be particularly challenging, since a good initial...

Standard Galerkin finite element methods for variably saturated groundwater flow have several deficiencies. For instance, local oscillations can appear around sharp infiltration fronts without the use of mass-lumping, and velocity fields obtained from differentiation of pressure fields are discontinuous at element boundaries. Here, we consider conforming finite element discretizations based on ...

A new computer program, 1DTempPro, is presented for the analysis of vertical one-dimensional (1D) temperature profiles under saturated flow conditions. 1DTempPro is a graphical user interface to the U.S. Geological Survey code Variably Saturated 2-Dimensional Heat Transport (VS2DH), which numerically solves the flow and heat-transport equations. Pre- and postprocessor features allow the user to...

A simple and efficient computational modelling framework is presented for the simulation of variably saturated flows in porous media, which is not based on a conventional numerical solution of the Richards’ equation. The computational domain is discretized with a regular grid and simple rules govern the flow dynamics. Each cell can have distinct hydraulic properties and its state depends only o...