Journal
ADVANCES IN WATER RESOURCES
Volume 152, Issue -, Pages -Publisher
ELSEVIER SCI LTD
Keywords
Richards equation; Coupled flow and transport; Linearization; Iterative schemes; Global random walk
Categories
Funding
- Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) [SU 415/4-1 - 405338726]
- Norwegian Academy of Science and Letters [6367]
- Equinor [6367]
Ask authors/readers for more resources
The article presents new random walk methods to solve flow and transport problems in porous media, using global random walk algorithms to approximate the solution and employing linearization techniques to handle nonlinearity and degeneracy. The resulting GRW L-schemes converge with the number of iterations and provide numerical solutions that are first-order accurate in time and second-order in space. These methods are practically free of numerical diffusion and are validated through comparisons with other solvers in benchmark problems.
In this article, we present new random walk methods to solve flow and transport problems in saturated/unsaturated porous media, including coupled flow and transport processes in soils, heterogeneous systems modeled through random hydraulic conductivity and recharge fields, processes at the field and regional scales. The numerical schemes are based on global random walk algorithms (GRW) which approximate the solution by moving large numbers of computational particles on regular lattices according to specific random walk rules. To cope with the nonlinearity and the degeneracy of the Richards equation and of the coupled system, we implemented the GRW algorithms by employing linearization techniques similar to the L-scheme developed in finite element/volume approaches. The resulting GRW L-schemes converge with the number of iterations and provide numerical solutions that are first-order accurate in time and second-order in space. A remarkable property of the flow and transport GRW solutions is that they are practically free of numerical diffusion. The GRW solvers are validated by comparisons with mixed finite element and finite volume solvers in one-and two-dimensional benchmark problems. They include Richards' equation fully coupled with the advection-diffusion-reaction equation and capture the transition from unsaturated to saturated flow regimes.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available