Numerical modeling of the leak through semipermeable walls for 2D/3D Stokes flow: Experimental scalability of dual algorithms

Abstract

The paper deals with the Stokes flow subject to the threshold leak boundary conditions in two and three space dimensions. The velocity-pressure formulation leads to the inequality type problem that is approximated by the P1-bubble/P1 mixed finite elements. The resulting algebraic system is nonsmooth. It is solved by the path-following variant of the interior point method, and by the active-set implementation of the semi-smooth Newton method. Inner linear systems are solved by the preconditioned conjugate gradient method. Numerical experiments illustrate scalability of the algorithms. The novelty of this work consists in applying dual strategies for solving the problem.