Global random walk solvers for fully coupled flow and transport in saturated/unsaturated porous media (extended version)

by   Nicolae Suciu, et al.

In this article, we present new random walk methods to solve flow and transport problems in unsaturated/saturated 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 solutions are validated by comparisons with mixed finite element and finite volume solutions 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. For completeness, we also consider decoupled flow and transport model problems for saturated aquifers.



There are no comments yet.


page 13

page 17

page 18


Global random walk solvers for reactive transport and biodegradation processes in heterogeneous porous media

Flow and multicomponent reactive transport in saturated/unsaturated poro...

Benchmark for numerical solutions of flow in heterogeneous groundwater formations

This article presents numerical investigations on accuracy and convergen...

Random walk approximation for irreversible drift-diffusion process on manifold: ergodicity, unconditional stability and convergence

Irreversible drift-diffusion processes are very common in biochemical re...

Branching random walk solutions to the Wigner equation

The stochastic solutions to the Wigner equation, which explain the nonlo...

Space-time upscaling of reactive transport in porous media

Reactive transport in saturated/unsaturated porous media is numerically ...

A block-coupled Finite Volume methodology for problems of large strain and large displacement

A nonlinear block-coupled Finite Volume methodology is developed for lar...

Numerical Simulations of Electrohydrodynamics flow model based on Burgers' Equation with Transport of Bubbles

In this paper we present numerical models for electrodynamical flows wit...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.