Weakly imposed Dirichlet boundary conditions for the Brinkman model of porous media flow

We use low order approximations, piecewise linear, continuous velocities and piecewise constant pressures to compute solutions to Brinkman’s equation of porous media flow, applying an edge stabilization term to avoid locking. In order to handle the limiting case of Darcy flow, when only the velocity component normal to the boundary can be prescribed, we impose the boundary conditions weakly using Nitsche’s method [J. Nitsche, Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind, Abhandlungen aus dem Mathematischen Seminar der Universität Hamburg 36 (1971) 9–15]. We show that this leads to a stable method for all choices of material parameters. Finally we present some numerical examples verifying the theoretical predictions and showing the effect of the weak imposition of boundary conditions. © 2008 IMACS. Published by Elsevier B.V. All rights reserved.