The multigrid solution of coupled porous media and Stokes flow problems is considered. The Darcy equation as the saturated porous medium model is coupled to the Stokes equations by means of appropriate interface conditions. We focus on an efficient multigrid solution technique for the coupled problem, which is discretized by finite volumes on staggered grids, giving rise to a saddle point linear system. Special treatment is required regarding the discretization at the interface. An Uzawa smoother is employed in multigrid, which is a decoupled procedure based on symmetric Gauss–Seidel smoothing for velocity components and a simple Richardson iteration for the pressure field. Since a relaxation parameter is part of a Richardson iteration, local Fourier analysis is applied to determine the optimal parameters. Highly satisfactory multigrid convergence is reported, and, moreover, the algorithm performs very well for small values of the hydraulic conductivity and fluid viscosity, which are relevant for applications.

, , , , , , , ,
SIAM Journal on Scientific Computing
Scientific Computing

Luo, P, Rodrigo, C, Gaspar, F.J, & Oosterlee, C.W. (2017). Uzawa smoother in multigrid for the coupleD porous medium and stokes flow system. SIAM Journal on Scientific Computing, 39(5), S633–S661. doi:10.1137/16M1076514