A strongly conservative finite element method for the coupling of darcy and stokes flow JCP 2016.pdfVIP

Journal of Computational Physics 229 2010 5933–5943 Contents lists available at ScienceDirect Journal of Computational Physics journal home A strongly conservative ?nite element method for the coupling of Stokes and Darcy ?ow G. Kanschat a,*,1, B. Rivière b,2 a Department of Mathematics, Texas AM University, 3368 TAMU, College Station, TX 77843-3368, United States b Department of Computational and Applied Mathematics, Rice University, 6100 Main Street, MS-134, Houston, TX 77005-1892, United States a r t i c l e i n f o a b s t r a c t Article history: We consider a model of coupled free and porous media ?ow governed by Stokes and Darcy Received 12 August 2009 equations with the Beavers–Joseph–Saffman interface condition. This model is discretized Received in revised form 6 April 2010 using divergence-conforming ?nite elements for the velocities in the whole domain. Dis- Accepted 13 April 2010 continuous Galerkin techniques and mixed methods are used in the Stokes and Darcy sub- Available online 27 April 2010 div domains, respectively. This discretization is strongly conservative in H X and we show convergence. Numerical results validate our ?ndings and indicate optimal convergence Keywords: orders. Stokes ?ow 2010 Elsevier Inc. All rights reserved. Darcy ?ow Beavers–Joseph–Saffman condition Finite element method Mixed ?nite elements Discontinuous Galerkin Multiphysics 1. Introduction The coupling of Stokes and Darcy equations arises from the modeling of groundwater contamination through streams and ?ltration problems [25,19]. In this work, a new numerical method is proposed, that employs divergence-conforming velocity spaces, i.e. spaces included in Hdiv X . The Darcy ?ow is discretized by a mixed ?nite element method and the Stokes ?ow by a mixed velocity–pressure discontinuous Galerkin DG method. The two types of ?ow are coupled by appropriate interface conditions, namely mass conservation, balance of f