Stable expression and control parameters in the numerical simulation of unsaturated flow

Abstract

The Richards’ equation describes the flow phenomenon in unsaturated porous media and is essential to hydrology and environmental science. This study evaluated the numerical stability of two different forms of the Richards’ equation. Sensitivity analyses were performed to investigate the control parameters of the equation. The results show that the h-form Richards’ equation has better applicability for calculating variable saturation flows than the θ-form Richards’ equation. For the h-form Richards’ equation, the hydraulic conductivity of the soil in the low-suction range and the specific moisture capacity in the high-suction range primarily influenced the solution. In addition, sensitivity analyses indicated that the saturated hydraulic conductivity, initial condition, and air-entry pressure have a higher sensitivity to the simulation results than the saturated water content, rainfall intensity, and decline rate of hydraulic conductivity. Moreover, their correctness needs to be guaranteed first in numerical simulations. The research findings can provide a helpful reference for improving the reliability of numerical simulations of unsaturated flows.

Graphical abstract

Keywords

Richards’ equation / sensitivity analysis / unsaturated soil / hydraulic diffusivity / seepage simulation

References

[1] Park J S R, Cheung S W, Mai T. Multiscale simulations for multi-continuum Richards equations. Journal of Computational and Applied Mathematics, 2021, 397: 113648

[2] Richards L A. Capillary conduction of liquids through porous mediums. Physics, 1931, 1(5): 318–333

[3] Farthing M W, Ogden F L. Numerical solution of Richards’ equation: A review of advances and challenges. Soil Science Society of America Journal, 2017, 81(6): 1257–1269

[4] Li Z, Özgen-Xian I, Maina F Z. A mass-conservative predictor-corrector solution to the 1D Richards equation with adaptive time control. Journal of Hydrology (Amsterdam), 2021, 592: 125809

[5] Paniconi C, Putti M. Physically based modeling in catchment hydrology at 50: Survey and outlook. Water Resources Research, 2015, 51(9): 7090–7129

[6] Klute A. The determination of the hydraulic conductivity and diffusivity of unsaturated soils1. Soil Science, 1972, 113(4): 264–276

[7] Li X, Zhang L M, Fredlund D G. Wetting front advancing column test for measuring unsaturated hydraulic conductivity. Canadian Geotechnical Journal, 2009, 46(12): 1431–1445

[8] Liu Q, Xi P, Miao J, Li X, Wang K. Applicability of wetting front advancing method in the sand to silty clay soils. Soil and Foundation, 2020, 60(5): 1215–1225

[9] Ma D, Wang Q, Shao M. Analytical method for estimating soil hydraulic parameters from horizontal absorption. Soil Science Society of America Journal, 2009, 73(3): 727–736

[10] Su W, Cui Y J, Qin P J, Zhang F, Ye W M, Conil N. Application of instantaneous profile method to determine the hydraulic conductivity of unsaturated natural stiff clay. Engineering Geology, 2018, 243: 111–117

[11] Li X, Zhang Z, Zhang L, Zhang L, Wu L. Combining two methods for the measurement of hydraulic conductivity over a wide suction range. Computers and Geotechnics, 2021, 135: 104178

[12] Bruce R R, Klute A. The measurement of soil moisture diffusivity1. Soil Science Society of America Journal, 1956, 20(4): 458

[13] Morel-Seytoux H J, Khanji J. Prediction of imbibition in a horizontal column. Soil Science Society of America Journal, 1975, 39(4): 613–617

[14] NingLLikos W J. Unsaturated Soil Mechanics. New Jerse: Wiley, 2004

[15] BrooksR HCorey A T. Hydraulic Properties of Porous Media. Fort Collins: Hydrol Pap, 1964

[16] van Genuchten M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 1980, 44(5): 892–898

[17] Fredlund D G, Xing A. Equations for the soil-water characteristic curve. Canadian Geotechnical Journal, 1994, 31(4): 521–532

[18] Zha Y, Yang J, Shi L, Song X. Simulating one-dimensional unsaturated flow in heterogeneous soils with water content-based Richards equation. Vadose Zone Journal, 2013, 12(2): vzj2012.0109

[19] Tai B, Yue Z, Qi S, Wang P. Experimental and numerical investigation on thermal-moisture-mechanical behaviors on a new anti-frost cutting bed of high-speed railway in deep seasonally frozen ground regions under extreme climate. Computers and Geotechnics, 2021, 136: 104251

[20] Li X, Li X, Wu Y, Wu L, Yue Z. Selection criteria of mesh size and time step in fem analysis of highly nonlinear unsaturated seepage process. Computers and Geotechnics, 2022, 146: 104712

[21] Tai B, Wu Q, Zhang Z, Xu X. Study on thermal performance of novel asymmetric crushed-rock-based embankment on the qinghai-tibet railway in permafrost region. International Journal of Thermal Sciences, 2020, 152: 106333

[22] Fidkowski K J, Darmofal D L. Review of output-based error estimation and mesh adaptation in computational fluid dynamics. AIAA Journal, 2011, 49(4): 673–694

[23] Hénon P, Ramet P, Roman J. Pastix: A high-performance parallel direct solver for sparse symmetric positive definite systems. Parallel Computing, 2002, 28(2): 301–321

[24] LiX. An Overview of SuperLU. New York: ACM Transactions on Mathematical Software (TOMS), 2005

[25] List F, Radu F A. A study on iterative methods for solving Richards’ equation. Computational Geosciences, 2015, 20(2): 1–13

[26] Pettway J S, Schmidt J H, Stagg A K. Adaptive meshing in a mixed regime hydrologic simulation model. Computational Geosciences, 2010, 14(4): 665–674

[27] Pop I S, Radu F, Knabner P. Mixed finite elements for the Richards’ equation: Linearization procedure. Journal of Computational and Applied Mathematics, 2004, 168(1−2): 365–373

[28] Solin P, Kuraz M. Solving the nonstationary Richards equation with adaptive HP-FEM. Advances in Water Resources, 2011, 34(9): 1062–1081

[29] Valdez A R, Rocha B M, Chapiro G, Weber dos Santos R. Uncertainty quantification and sensitivity analysis for relative permeability models of two-phase flow in porous media. Journal of Petroleum Science Engineering, 2020, 192: 107297

[30] Sieber A, Uhlenbrook S. Sensitivity analyses of a distributed catchment model to verify the model structure. Journal of Hydrology (Amsterdam), 2005, 310(1−4): 216–235

[31] Berg S, Unsal E, Dijk H. Sensitivity and uncertainty analysis for parameterization of multiphase flow models. Transport in Porous Media, 2021, 140(1): 27–57

[32] Crosetto M, Tarantola S, Saltelli A. Sensitivity and uncertainty analysis in spatial modelling based on GIS. Agriculture, Ecosystems & Environment, 2000, 81(1): 71–79

[33] Morris M D. Factorial sampling plans for preliminary computational experiments. Technometrics, 1991, 33(2): 161–174

[34] Campolongo F, Tarantola S, Saltelli A. Tackling quantitatively large dimensionality problems. Computer Physics Communications, 1999, 117(1−2): 75–85

[35] Herman J D, Kollat J B, Reed P M, Wagener T. Technical note: Method of Morris effectively reduces the computational demands of global sensitivity analysis for distributed watershed models. Hydrology and Earth System Sciences, 2013, 17(7): 2893–2903

[36] Pang M, Xu R, Hu Z, Wang J, Wang Y. Uncertainty and sensitivity analysis of input conditions in a large shallow lake based on the latin hypercube sampling and morris methods. Water (Basel), 2021, 13(13): 1861

[37] King D M, Perera B J C. Morris method of sensitivity analysis applied to assess the importance of input variables on urban water supply yield—A case study. Journal of Hydrology (Amsterdam), 2013, 477: 17–32

[38] Ren J, Zhang W, Yang J. Morris sensitivity analysis for hydrothermal coupling parameters of embankment dam: A case study. Mathematical Problems in Engineering, 2019, 2019: 1–11

[39] Yi X, Zou R, Guo H. Global sensitivity analysis of a three-dimensional nutrients-algae dynamic model for a large shallow lake. Ecological Modelling, 2016, 327: 74–84

[40] Niessner J, Helmig R, Jakobs H, Roberts J E. Interface condition and linearization schemes in the newton iterations for two-phase flow in heterogeneous porous media. Advances in Water Resources, 2005, 28(7): 671–687

[41] Mualem Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Research, 1976, 12(3): 513–522

[42] CrankJ. The Mathematics of Diffusion. New York: Oxford University Press, 1979

[43] Tracy F T. Clean two- and three-dimensional analytical solutions of Richards’ equation for testing numerical solvers. Water Resources Research, 2006, 42(8): W08503

[44] Gibson R E. A note on the constant head test to measure soil permeability in situ. Geotechnique, 1966, 16(3): 256–259

[45] Gibson R E. An extension to the theory of the constant head in situ permeability test. Geotechnique, 1970, 20(2): 193–197

[46] Powrie W. Soil Mechanics: Concepts and Applications. 2nd ed. Oxford: Taylor & Francis, 2013, 22(3): 381–382

[47] Saltelli A, Sobol I M. About the use of rank transformation in sensitivity analysis of model output. Reliability Engineering & System Safety, 1995, 50(3): 225–239

[48] Zádor J, Zsély I G, Turányi T. Local and global uncertainty analysis of complex chemical kinetic systems. Reliability Engineering & System Safety, 2006, 91(10-11): 1232–1240

[49] Francos A, Elorza F J, Bouraoui F, Bidoglio G, Galbiati L. Sensitivity analysis of distributed environmental simulation models: Understanding the model behaviour in hydrological studies at the catchment scale. Reliability Engineering & System Safety, 2003, 79(2): 205–218

[50] Chui T F M, Freyberg D L. Implementing hydrologic boundary conditions in a multiphysics model. Journal of Hydrologic Engineering, 2009, 14(12): 1374–1377

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grant Nos. 51979002 and 52025094).

Open Access

This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/ by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.

RIGHTS & PERMISSIONS

2022 The Author(s) 2022. This article is published with open access at link.springer.com and journal.hep.com.cn