HR: 0830h
AN: H11G-0926    [PDF]
TI: Higher Order Time Integration and Discontinuous Galerkin Methods for Variably Saturated Groundwater Flow
AU: Farthing, M W
EM: matthew_farthing@unc.edu
AF: Center for the Advanced Study of the Environment, Department of Environmental Sciences and Engineering, University of North Carolina, CB #7431, Rosenau Hall, Chapel Hill, NC 27599-7431 United States
AU: * Li, H
EM: huinali@email.unc.edu
AF: Center for the Advanced Study of the Environment, Department of Environmental Sciences and Engineering, University of North Carolina, CB #7431, Rosenau Hall, Chapel Hill, NC 27599-7431 United States
AU: Miller, C T
EM: casey_miller@unc.edu
AF: Center for the Advanced Study of the Environment, Department of Environmental Sciences and Engineering, University of North Carolina, CB #7431, Rosenau Hall, Chapel Hill, NC 27599-7431 United States
AU: Kees, C E
EM: chris_kees@ncsu.edu
AF: Center for Research in Scientific Computation, Department of Mathematics, North Carolina State University, Box 8205, Harrelson Hall, Raleigh, NC 27695-8205 United States
AB: The numerical simulation of groundwater flow in the vadose zone continues to be a challenge for many problems of practical interest. Under commonly used constitutive relations, the governing equations can be highly nonlinear and produce sharp fronts in the solution variables for problems such as wetting phase infiltration into an initially dry medium. For a number of multiphase flow problems, the use of variable order and variable step size temporal discretizations has shown significant advantages. However, the spatial discretizations commonly used for variably saturated flow are dominated by low-order finite difference and finite element methods. Over the last decade discontinuous Galerkin (DG) finite element methods have received significant attention in a number of fields for hyperbolic PDE's and, more recently, for elliptic and parabolic problems. DG approaches are appealing for modeling subsurface flow since they can lead to velocity fields that are locally mass-conserving without the need for auxiliary variables or alternative meshes. Moreover, DG discretizations are inherently local and so well-suited for unstructured meshes and $h$-$p$ adaption strategies. While some work has been done recently for multiphase subsurface flow, there are a range of issues related to the performance of DG methods for highly nonlinear parabolic problems that have not been investigated fully, particularly for air-water systems. In this work, we consider the combination of higher order adaptive time integration with DG spatial discretizations applied to variably saturated groundwater flow. We compare this approach to standard low order methods for a series of test problems and consider a number of issues including the methods' relative accuracy and computational efficiency.
DE: 1829 Groundwater hydrology
DE: 1875 Unsaturated zone
SC: Hydrology [H]
MN: 2003 Fall Meeting