Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
156 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
45 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

A semi-implicit second-order temporal scheme for solving the pressure head-based form of Richards' and advection-dispersion equations (2404.03603v3)

Published 4 Apr 2024 in math.NA and cs.NA

Abstract: In this study, a novel semi-implicit second-order temporal scheme combined with the finite element method for space discretization is proposed to solve the coupled system of infiltration and solute transport in unsaturated porous media. The Richards equation is used to describe unsaturated flow, while the advection-dispersion equation (ADE) is used for modeling solute transport. The proposed approach is used to linearize the system of equations in time, eliminating the need of iterative processes. A free parameter is introduced to ensure the stability of the scheme. Numerical tests are conducted to analyze the accuracy of the proposed method in comparison with a family of second-order iterative schemes. The proposed numerical technique based on the optimal free parameter is accurate and performs better in terms of efficiency since it offers a considerable gain in computational time compared to the other methods. For reliability and effectiveness evaluation of the developed semi-implicit scheme, four showcase scenarios are used. The first two numerical tests focus on modeling water flow in heterogeneous soil and transient flow in variably saturated zones. The last numerical tests are carried out to simulate the salt and nitrate transport through unsaturated soils. The simulation results are compared with reference solutions and laboratory data, and demonstrate the effectiveness of the proposed scheme in simulating infiltration and solute transport through unsaturated soils.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (68)
  1. A family of new globally convergent linearization schemes for solving Richards’ equation. Applied Numerical Mathematics, 159:281–296.
  2. Quasilinear elliptic-parabolic differential equations. Math. Z, 183(3):311–341.
  3. Discontinuous Galerkin approximation of flows in fractured porous media on polytopic grids. SIAM Journal on Scientific Computing, 41(1):A109–A138.
  4. Implicit-explicit methods for time-dependent partial differential equations. SIAM Journal on Numerical Analysis, 32(3):797–823.
  5. Adaptive multistep time discretization and linearization based on a posteriori error estimates for the Richards equation. Applied Numerical Mathematics, 112:104–125.
  6. Bear, J. (2013). Dynamics of fluids in porous media. Courier Corporation.
  7. An efficient lumped mixed hybrid finite element formulation for variably saturated groundwater flow. Vadose Zone Journal, 8:352–362.
  8. A continuum model of unstable infiltration in porous media endowed with an entropy function. Advances in Water Resources, 144:103684.
  9. Theoretical and numerical analysis of a class of semi-implicit semi-lagrangian schemes potentially applicable to atmospheric models. Monthly Weather Review, 142:4458–4476.
  10. Localized MQ-RBF meshless techniques for modeling unsaturated flow. Engineering Analysis with Boundary Elements, 130:109–123.
  11. Localized RBF methods for modeling infiltration using the Kirchhoff-transformed Richards equation. Engineering Analysis with Boundary Elements, 152:259–276.
  12. Implicit EXP-RBF techniques for modeling unsaturated flow through soils with water uptake by plant roots. Applied Numerical Mathematics. Positive review.
  13. Properties of porous media affecting fluid flow. Journal of the Irrigation and Drainage Division, American Society of Civil Engineers, 92:61–88.
  14. Hydrogeosphere: a fully integrated, physically based hydrological model. Ground water, 50(2):170–176.
  15. A general mass-conservative numerical solution for the unsaturated flow equation. Water Resources Research, 26:1483–1496.
  16. An adaptive strategy for discontinuous Galerkin simulations of Richards’ equation: Application to multi-materials dam wetting. Advances in Water Resources, 151.
  17. Development of a physically based 1d-infiltration model for irrigated soils. Agricultural Water Management, 85(1):165–174.
  18. Davis, T. A. (2004). Algorithm 832: UMFPACK V4.3 - An unsymmetric-pattern multifrontal method. ACM Transactions on Mathematical Software, 30(2):196 – 199. Cited by: 1072.
  19. Durran, D. R. (2010). Numerical methods for fluid dynamics: With applications to geophysics, volume 32. Springer Science & Business Media.
  20. The finite volume method for Richards equation. Computational Geosciences, 3:259–294.
  21. Numerical solution of Richards’ equation: A review of advances and challenges. Soil Science Society of America Journal, 81(6):1257–1269.
  22. Groundwater prentice-hall inc. Eaglewood Cliffs, NJ.
  23. Gardner, W. R. (1958). Some steady state solutions of the unsaturated moisture flow equation with application to evaporation from a water table. Soil Science, 85:228–232.
  24. Finite element methods for Navier-Stokes equations: theory and algorithms, volume 5. Springer Science & Business Media.
  25. Modeling water flow and solute transport in unsaturated soils using physics-informed neural networks trained with geoelectrical data. Water Resources Research, 59(6).
  26. Modeling water flow and solute transport in unsaturated soils using physics-informed neural networks trained with geoelectrical data. https://doi.org/10.5281/zenodo.7558746. [Dataset].
  27. A comparison of numerical simulation models for one-dimensional infiltration. Soil Science Society of America Journal, 41(2):285–294.
  28. Hecht, F. (2012a). New development in freefem++. Journal of Numerical Mathematics, 20(3-4):251 – 265. Cited by: 2275; All Open Access, Green Open Access.
  29. Hecht, F. (2012b). New development in freefem++. J. Numer. Math., 20(3-4):251–265.
  30. A new convergence criterion for the modified Picard iteration method to solve the variably saturated flow equation. Journal of Hydrology, 178(1-4):69–91.
  31. Numerical solution of time-dependent advection-diffusion-reaction equations, volume 33. Springer.
  32. Finite element model of nitrogen species transformation and transport in the unsaturated zone. Journal of Hydrology, 103(3-4):249–274.
  33. Efficient second-order semi-implicit finite element method for fourth-order nonlinear diffusion equations. Computer Physics Communications, 258:107588.
  34. Implicit and semi-implicit second-order time stepping methods for the Richards equation. Advances in Water Resources, 148.
  35. Lambert, J. D. et al. (1991). Numerical methods for ordinary differential systems, volume 146. Wiley New York.
  36. Levrett, M. C. (1941). Capillary behavior of porous solids. Trans. AIME, 142:152–169.
  37. Modeling of water flow and nitrate transport under surface drip fertigation. Transactions of the ASAE, 48(2):627–637.
  38. Water and nitrogen distribution as affected by fertigation of ammonium nitrate from a point source. Irrigation Science, 22:19–30.
  39. A study on iterative methods for solving Richards’ equation. Computational Geosciences, 20:341–353.
  40. An accelerated Picard method for nonlinear systems related to variably saturated flow. Advances in Water Resources, 38:92–101.
  41. Mass-conservative finite volume methods on 2-D unstructured grid for the Richards’ equation. Adv. Water Resour., 27:1199–1215.
  42. Permeability of porous solids. Transactions of the Faraday Society, 57:1200–1207.
  43. A modified L-scheme to solve nonlinear diffusion problems. Computers and Mathematics with Applications, 77:1722–1738.
  44. Mualem, Y. (1976). A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Research, pages 513–522.
  45. Nishikawa, H. (2019). On large start-up error of BDF2. Journal of Computational Physics, 392:456–461.
  46. Philip, J. R. (1969). Theory of infiltration. In Advances in hydroscience, volume 5, pages 215–296. Elsevier.
  47. Finite element simulation in surface and subsurface hydrology. Elsevier.
  48. Finite element analysis in geotechnical engineering: application, volume 2. Thomas Telford London.
  49. Richards, L. A. (1931). Capillary conduction of liquids through porous mediums. Journal of Applied Physics, 1:318–333.
  50. Numerical analysis of one-dimensional unsaturated flow in layered soils. Advances in Water Resources, 21(4):315–324.
  51. Ross, P. J. (2003). Modeling soil water and solute transport–fast, simplified numerical solutions. Agronomy Journal, 95:1352–1361.
  52. Russo, D. (2012). Numerical analysis of solute transport in variably saturated bimodal heterogeneous formations with mobile–immobile-porosity. Advances in water resources, 47:31–42.
  53. Russo, G. (2002). Central schemes and systems of balance laws, in “hyperbolic partial differential equations. Theory, Numerics and Applications”, edited by Andreas Meister and Jens Struckmeier, Vieweg, Göttingen.
  54. Modeling nitrogen transport and transformations in soils: 1. theoretical considerations: 1. Soil Science, 131(4):233–241.
  55. Modeling compensated root water and nutrient uptake. Ecological Modelling, 220(4):505–521.
  56. The HYDRUS-1D software package for simulating the one-dimensional movement of water, heat, and multiple solutes in variably-saturated media. University of California-Riverside Research Reports, 3:1–240.
  57. The HYDRUS software package for simulating two-and three-dimensional movement of water, heat, and multiple solutes in variably-saturated media. Technical manual, version, 1:241.
  58. Recent developments and applications of the HYDRUS computer software packages. Vadose Zone Journal, 15(7).
  59. Minimizing nitrogen leaching from furrow irrigation through novel fertilizer placement and soil surface management strategies. Agric. Water Manag., 115:242–251.
  60. Strategies to reduce nitrate leaching under furrow irrigation. Int. J. Environ. Sci. Dev., 4:431–434.
  61. Modeling of flow and advection dominant solute transport in variably saturated porous media. Journal of Hydrologic Engineering, 14(1):1–14.
  62. A three-dimensional numerical model for water flow and transport of chemically reactive solute through porous media under variably saturated conditions. Advances in water resources, 15(5):275–287.
  63. Tracy, F. T. (2011). Analytical and numerical solutions of Richards’ equation with discussions on relative hydraulic conductivity. BoD—Books on demand.
  64. van Genuchten, M. (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44:892–898.
  65. Studies on soil physics I. The flow of air and water trough soils. J Agri Sci, 4:1–24.
  66. A new mass lumping scheme for the mixed hybrid finite element method. International Journal for Numerical Methods in Engineering, 67:89–107.
  67. A robust upwind mixed hybrid finite element method for transport in variably saturated porous media. Hydrology and Earth System Sciences, 26(20):5227–5239.
  68. Review of numerical solution of Richardson–Richards equation for variably saturated flow in soils. Wiley Interdisciplinary Reviews: Water, 6(5):e1364.

Summary

We haven't generated a summary for this paper yet.