Abstract
Nonlinearity continuation method, applied to boundary value problems for steady-state Richards equation, gradually approaches the solution through a series of intermediate problems. Originally, the Newton method with simple line search algorithm was used to solve the intermediate problems. In the present paper, other solvers such as Picard and mixed Picard–Newton methods are considered, combined with slightly modified line search approach. Numerical experiments are performed with advanced finite volume discretizations for model and real-life problems.
-
Funding: The reported study was partially funded by the Russian Foundation for Basic Research (RFBR), project number 20-31-90126.
References
[1] I. Aavatsmark, T. Barkve, O. B⊘e, and T. Mannseth, Discretization on unstructured grids for inhomogeneous, anisotropic media, Part I: Derivation of the methods. SIAM J. Sci. Comp. 19 (1998), No. 5,1700-1716.10.1137/S1064827595293582Search in Google Scholar
[2] E. L. Allgower and K. Georg, Introduction to Numerical Continuation Methods. SIAM, 2003.10.1137/1.9780898719154Search in Google Scholar
[3] D. Anuprienko and I. Kapyrin, Nonlinearity continuation method for steady-state groundwater flow modeling in variably saturated conditions. J. Comp. Appl. Math. (2021), 113502.10.1016/j.cam.2021.113502Search in Google Scholar
[4] D. Anuprienko and I. Kapyrin, Modeling groundwater flow in unconfined conditions: numerical model and solvers’ efficiency. Lobachevskii J. Math. 39 (2018), No. 7, 867-873.10.1134/S1995080218070053Search in Google Scholar
[5] L. Armijo, Minimization of functions having Lipschitz continuous first partial derivatives. Pacific J. Math. 16 (1966), No. 1, 1-3.10.2140/pjm.1966.16.1Search in Google Scholar
[6] J. Bear and A. H. D. Cheng, Modeling Groundwater Flow and Contaminant Transport. Springer Science & Business Media, 2010.10.1007/978-1-4020-6682-5Search in Google Scholar
[7] M. A. Celia, E. T. Bouloutas, and R. L. Zarba, A general mass-conservative numerical solution for the unsaturated flow equation. Water Resources Research 26 (1990), No. 7,1483-1496.10.1029/WR026i007p01483Search in Google Scholar
[8] H.-J. G. Diersch and P. Perrochet, On the primary variable switching technique for simulating unsaturated-saturated flows. Advances in Water Resources 23 (1999), No. 3, 271-301.10.1016/S0309-1708(98)00057-8Search in Google Scholar
[9] J. S. Duff and J. Koster, The design and use of algorithms for permuting large entries to the diagonal of sparse matrices. SIAM J. Matrix Anal. Appl. 20 (1999), No. 4, 889-901.10.1137/S0895479897317661Search in Google Scholar
[10] M. W. Farthing, C. Kees, T. S. Coffey, C. T. Kelley, and C. T. Miller, Efficient steady-state solution techniques for variably saturated groundwater flow. Advances in Water Resources 26 (2003), No. 8, 833-849.10.21236/ADA454444Search in Google Scholar
[11] K. R. Fowler and C. T. Kelley, Pseudo-transient continuation for nonsmooth nonlinear equations. SIAM J. Numer. Anal. 43 (2005), No. 4,1385-1406.10.1137/S0036142903431298Search in Google Scholar
[12] GeRa software website. http://en.ibrae.ac.ru/code/gera7. Accessed: 2021-04-14.Search in Google Scholar
[13] INMOST, a toolkit for distributed mathematical modeling. http://www.inmost.org/. Accessed: 2021-04-14.Search in Google Scholar
[14] J. E. Jones and C. S. Woodward, Newton–Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems. Advances in Water Resources 24 (2001), No. 7, 763-774.10.1016/S0309-1708(00)00075-0Search in Google Scholar
[15] R. Kao, A comparison of Newton-Raphson methods and incremental procedures for geometrically nonlinear analysis. Computers & Structures 4 (1974), No. 5,1091-1097.10.1016/0045-7949(74)90025-XSearch in Google Scholar
[16] I. V. Kapyrin, A family of monotone methods for the numerical solution of three-dimensional diffusion problems on unstructured tetrahedral meshes. In: DokladyMath. 76 (2 0 0 7), No. 2, 734-738.10.1134/S1064562407050249Search in Google Scholar
[17] I. V. Kapyrin, V. A. Ivanov, G. V. Kopytov, and S. S. Utkin, Integral code GeRa for RAW disposal safety validation. Gornyi Zhurnal (2015), No. 10, 44-50 (in Russian).10.17580/gzh.2015.10.08Search in Google Scholar
[18] I. Konshin, I. Kaporin, K. Nikitin, and Yu. Vassilevski, Parallel linear systems solution for multiphase flow problems in the INMOST framework. In: Russian Supercomputing Days (2015), 96-103.Search in Google Scholar
[19] C. T. Kelley, Iterative Methods for Linear and Nonlinear Equations. SIAM, 1995.10.1137/1.9781611970944Search in Google Scholar
[20] N. Li, Y. Saad, and E. Chow, Crout versions of ILU for general sparse matrices. SIAM J. Sci. Comp. 25 (2003), No. 2, 716-728.10.1137/S1064827502405094Search in Google Scholar
[21] K. Lipnikov, D. Moulton, and D. Svyatskiy, New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation. Advances in Water Resources 94 (2016), 11-22.10.1016/j.advwatres.2016.04.016Search in Google Scholar
[22] Y. Mualem, A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Research 12 (1976), No. 3, 513-522.10.1029/WR012i003p00513Search in Google Scholar
[23] L. A. Richards, Capillary conduction of liquids through porous mediums. Physics 1 (1931), No. 5, 318-333.10.1063/1.1745010Search in Google Scholar
[24] C. Paniconi and M. Putti, A comparison of Picard and Newton iteration in the numerical solution of multidimensional variably saturated flow problems. Water Resources Research 30 (1994), No. 12, 3357-3374.10.1029/94WR02046Search in Google Scholar
[25] J.-Y. Parlange, W. L. Hogarth, D. A. Barry, M. B. Parlange, R. Haverkamp, P. J. Ross, T. S. Steenhuis, D. A. DiCarlo, and G. Katul, Analytical approximation to the solutions of Richards’ equation with applications to infiltration, ponding, and time compression approximation. Advances in Water Resources 23 (1999), No. 2,189-194.10.1016/S0309-1708(99)00022-6Search in Google Scholar
[26] A. V. Plenkin, A. Yu. Chernyshenko, V. N. Chugunov, and I. V. Kapyrin, Adaptive unstructured mesh generation methods for hydrogeological problems. Numerical Methods and Programming 16 (2015), 518-533 (in Russian).10.26089/NumMet.v16r449Search in Google Scholar
[27] P. I. Polubarinova-Koch, Theory of Ground Water Movement. Princeton University Press, 2015.Search in Google Scholar
[28] I. S. Pop, F. Radu, and P. Knabner, Mixed finite elements for the Richards’ equation: linearization procedure. J. Comp. Appl. Math. 168 (2004), No. 1-2, 365-373.10.1016/j.cam.2003.04.008Search in Google Scholar
[29] J. N. Reddyand D. K. Gartling, The Finite Element Method in Heat Transfer and Fluid Dynamics. CRC press, 2010.10.1201/9781439882573Search in Google Scholar
[30] W. F. Schmidt, Adaptive step size selection for use with the continuation method. Int.J. Numer. Methods Engrg. 12 (1978), No. 4, 677-694.10.1002/nme.1620120411Search in Google Scholar
[31] K. M. Terekhov, B. T. Mallison, and H. A. Tchelepi, Cell-centered nonlinear finite-volume methods for the heterogeneous anisotropic diffusion problem. J. Comp.Phys. 330 (2017), 245-267.10.1016/j.jcp.2016.11.010Search in Google Scholar
[32] K. Terekhov, Parallel multilevel linear solver within INMOST platform. In: Russian Supercomputing Days (2020), 297-309.10.1007/978-3-030-64616-5_26Search in Google Scholar
[33] M. Th. Van Genuchten, A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1. Soil Sci. Soc. America J. 44 (1980), No. 5, 892-898.10.2136/sssaj1980.03615995004400050002xSearch in Google Scholar
[34] Yu. Vassilevski, K. Terekhov, K. Nikitin, and I. Kapyrin, Parallel Finite Volume Computation on General Meshes. Springer Nature, 2020.10.1007/978-3-030-47232-0Search in Google Scholar
© 2021 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Comparison of nonlinear solvers within continuation method for steady-state variably saturated groundwater flow modelling
- A regularized isothermal phase-field model of two-phase solid–fluid mixture and its spatial dissipative discretization equations
- New correlative randomized algorithms for statistical modelling of radiation transfer in stochastic medium
- The suite of Taylor–Galerkin class schemes for ice transport on sphere implemented by the INMOST package
- Stability analysis of implicit semi-Lagrangian methods for numerical solution of non-hydrostatic atmospheric dynamics equations
Articles in the same Issue
- Frontmatter
- Comparison of nonlinear solvers within continuation method for steady-state variably saturated groundwater flow modelling
- A regularized isothermal phase-field model of two-phase solid–fluid mixture and its spatial dissipative discretization equations
- New correlative randomized algorithms for statistical modelling of radiation transfer in stochastic medium
- The suite of Taylor–Galerkin class schemes for ice transport on sphere implemented by the INMOST package
- Stability analysis of implicit semi-Lagrangian methods for numerical solution of non-hydrostatic atmospheric dynamics equations