A mass-conservative control volume-finite element method for solving Richards’ equation in heterogeneous porous media
Tóm tắt
We present a mass-conservative vertex-centred finite volume method for efficiently solving the mixed form of Richards’ equation in heterogeneous porous media. The spatial discretisation is particularly well-suited to heterogeneous media because it produces consistent flux approximations at quadrature points where material properties are continuous. Combined with the method of lines, the spatial discretisation gives a set of differential algebraic equations amenable to solution using higher-order implicit solvers. We investigate the solution of the mixed form using a Jacobian-free inexact Newton solver, which requires the solution of an extra variable for each node in the mesh compared to the pressure-head form. By exploiting the structure of the Jacobian for the mixed form, the size of the preconditioner is reduced to that for the pressure-head form, and there is minimal computational overhead for solving the mixed form. The proposed formulation is tested on two challenging test problems. The solutions from the new formulation offer conservation of mass at least one order of magnitude more accurate than a pressure head formulation, and the higher-order temporal integration significantly improves both the mass balance and computational efficiency of the solution.
Tài liệu tham khảo
Bear, J.: Hydraulics of Groundwater. McGraw-Hill, New York (1979)
Brenan, K., Campbell, S., Petzold, L.: Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations. Classics in Applied Mathematics, vol. 14. SIAM, Philadelphia (1996)
Celia, M., Bouloutas, E., Zarba, R.: A general mass-conservative solution for the unsaturated flow equation. Water Resour. Res. 26, 1483–1496 (1990)
Diersch, H., Kolditz, O.: Coupled groundwater flow and transport: 2. Thermohaline and 3D convection systems. Adv. Water Resour. 21(5), 401–425 (1998)
Diersch, H., Kolditz, O.: Variable-density flow and transport in porous media: approaches and challenges. Adv. Water Resour. 25(8–12), 899–944 (2002)
Diersch, H., Perrochet, P.: On the primary variable switching technique for simulating unsaturated-saturated flows. Adv. Water Resour. 23(3), 271–301 (1999)
Fahs, M., Younes, A., Lehmann, F.: An easy and efficient combination of the mixed finite element method and the method of lines for the resolution of Richards’ equation. Environ. Model. Softw. 24(9), 1122–1126 (2009)
Farthing, M.W., Kees, C.E., Miller, C.T.: Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv. Water Resour. 26(4), 373–394 (2003)
Forsyth, P.: A control volume finite element approach to NAPL groundwater contamination. SIAM J. Sci. Stat. Comput. 12(5), 1029–1057 (1991)
Forsyth, P., Unger, A., Sudicky, E.: Nonlinear iteration methods for nonequilibrium multiphase subsurface flow. Adv. Water Resour. 21(6), 433–449 (1998)
Forsyth, P., Wu, Y., Pruess, K.: Robust numerical methods for saturated–unsaturated flow with dry initial conditions in heterogeneous media. Adv. Water Resour. 18, 25–38 (1995)
Geuzaine, C., Remacle, J.: Gmsh: a three-dimensional finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Methods Eng. 79(11), 1309–1331 (2009)
Helmig, R.: Multiphase Flow and Transport in Processes in the Subsurface: A Contribution to the Modeling of Hydrosystems. Springer, Berlin (1997)
Hindmarsh, A., Brown, P., Grant, K., Lee, S., Serban, R., Shumaker, D., Woodward, C.: SUNDIALS: Suite of nonlinear and differential/algebraic equation solvers. ACM Trans. Math. Sci. 31(3), 363–396 (2005)
Huyakorn, P., Andersen, P., Mercer, J., White, H. Jr: Saltwater intrusion in aquifers: development and testing of a three-dimensional finite element method. Water Resour. Res. 23, 292–312 (1987)
Kees, C., Miller, C.: Higher order time integration methods for two-phase flow. Adv. Water Resour. 25(2), 159–177 (2002)
Kelley, C.: Iterative Methods for Linear and Nonlinear Equations. Frontiers in Applied Mathematics, vol. 16. SIAM, Philadelphia (1995)
Kirkland, M., Hills, R., Wierenga, P.: Algorithms for solving Richards’ equation for variably-saturated soils. Water Resour. Res. 28(8), 2049–2058 (1992)
Knoll, D., Keyes, D.: Jacobian-free Newton-Krylov methods: a survey of approaches and applications. J. Comput. Phys. 193(2), 357–397 (2004)
Li, H., Farthing, M., Miller, C.: Adaptive local discontinuous Galerkin approximation to Richards’ equation. Adv. Water Resour. 30(9), 1883–1901 (2007)
Liu, F., Anh, V., Turner, I., Bajracharya, K., Huxley, W., Su, N.: A finite volume simulation model for saturated-unsaturated flow and application to Gooburrum, Bundaberg, Queensland, Australia. Appl. Math. Model. 30, 352–366 (2006)
Manzini, G., Ferraris, S.: Mass-conservative finite volume methods on 2-D unstructured grids for the Richards’ equation. Adv. Water Resour. 27(12), 1199–1215 (2004)
Miller, C., Williams, G., Kelley, C., Tocci, M.: Robust solution of Richards’ equation for nonuniform porous media. Water Resour. Res. 34(10), 2599–2610 (1998)
Moroney, T., Turner, I.: A finite volume method based on radial basis functions for two-dimensional nonlinear diffusion equations. Appl. Math. Model. 30(10), 1118–1133 (2006)
Tocci, M., Kelley, C., Miller, C.: Accurate and economical solution of the pressure-head form of Richards’ equation by the method of lines. Adv. Water Resour. 20(1), 1–14 (1997)
Tocci, M., Kelley, C., Miller, C., Kees, C.: Inexact Newton methods and the method of lines for solving Richards’ equation in two space dimensions. Comput. Geosci. 2(4), 291–309 (1998)
Turner, I., Ferguson, W.: An unstructured mesh cell-centred control-volume method for simulating heat and mass transfer in porous media: application to softwood drying, Part I: The isotropic model. Appl. Math. Model. 19, 654–667 (1995)
Turner, I., Ferguson, W.: An unstructured mesh cell-centred control-volume method for simulating heat and mass transfer in porous media: application to softwood drying, Part II: The anisotropic model. Appl. Math. Model. 19, 668–674 (1995)
Van Genuchten, M.: Closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892–898 (1980)