<p>These days, infiltration techniques are frequently employed to control storm water in cities. Due to their advantages, which include lessening the adverse effects of urbanization, lowering storm water flow in sewage systems, and replenishing groundwater, these methods are employed and acknowledged globally. Richards’ equation (RE) is the governing equation for flow of water through unsaturated porous media. Our goal in this study is to create an effective numerical model for using the Crank-Nicolson finite difference method to solve the 3D mixed form of the RE in a heterogeneous layer soil medium. The system is linearized using a modified version of Picard’s iteration method. The suggested numerical model is validated and simulated through a series of 1D, 2D, and 3D numerical examples. Comparing the simulated results with analytic and Hydrus-1D software solution shows that the current approach is reliable, mass-conserving, and useful with quick convergence. The obtained results unequivocally show how well the suggested numerical model predicts the dynamics of soil moisture in a heterogeneous layer soil medium. The model’s encouraging outcomes enable its use in other contexts, including controlled aquifer recharge, plant water uptake, and the prediction of solute movement when combined with the advection-dispersion transport equation.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Efficient numerical modeling for the 3D mixed form of Richards’ equation in a heterogeneous layered soil with a modified Picard iteration approach

  • Sanjay L. Gosiya,
  • Vikas H. Pradhan

摘要

These days, infiltration techniques are frequently employed to control storm water in cities. Due to their advantages, which include lessening the adverse effects of urbanization, lowering storm water flow in sewage systems, and replenishing groundwater, these methods are employed and acknowledged globally. Richards’ equation (RE) is the governing equation for flow of water through unsaturated porous media. Our goal in this study is to create an effective numerical model for using the Crank-Nicolson finite difference method to solve the 3D mixed form of the RE in a heterogeneous layer soil medium. The system is linearized using a modified version of Picard’s iteration method. The suggested numerical model is validated and simulated through a series of 1D, 2D, and 3D numerical examples. Comparing the simulated results with analytic and Hydrus-1D software solution shows that the current approach is reliable, mass-conserving, and useful with quick convergence. The obtained results unequivocally show how well the suggested numerical model predicts the dynamics of soil moisture in a heterogeneous layer soil medium. The model’s encouraging outcomes enable its use in other contexts, including controlled aquifer recharge, plant water uptake, and the prediction of solute movement when combined with the advection-dispersion transport equation.