Home

About

Advanced Search

Browse by Discipline

Scientific Societies

E-print Alerts

Add E-prints

E-print Network
FAQHELPSITE MAPCONTACT US


  Advanced Search  

 
A COMPUTATONAL METHOD FOR APPROXIMATING A DARCY-STOKES SYSTEM
 

Summary: A COMPUTATONAL METHOD FOR APPROXIMATING
A DARCY-STOKES SYSTEM
GOVERNING A VUGGY POROUS MEDIUM
TODD ARBOGAST AND DANA S. BRUNSON
Abstract. We develop and analyze a mixed finite element method for the solution of an elliptic
system modeling a porous medium with large cavities, called vugs. It consists of a second order
elliptic (i.e., Darcy) equation on part of the domain coupled to a Stokes equation on the rest of the
domain, and a slip boundary condition (due to Beavers-Joseph-Saffman) on the interface between
them. The tangential velocity is not continuous on the interface. We consider a 2-D vuggy porous
medium with many small cavities throughout its extent, so the interface is not isolated. We use a
certain conforming Stokes element on rectangles, slightly modified near the interface to account for
the tangential discontinuity. This gives a mixed finite element method for the entire Darcy-Stokes
system with a regular sparsity pattern that is easy to implement, independent of the vug geometry,
as long as it aligns with the grid. We prove optimal global first order L2 convergence of the velocity
and pressure, as well as the velocity gradient in the Stokes domain. Numerical results verify these
rates of convergence, and even suggest somewhat better convergence in certain situations. Finally,
we present a lower dimensional space that uses Raviart-Thomas elements in the Darcy domain and
uses our new modified elements near the interface in transition to the Stokes elements.
Key words. mixed finite elements, Darcy-Stokes system, vuggy porous media, Beavers-Joseph
boundary condition, error estimates

  

Source: Arbogast, Todd - Center for Subsurface Modeling & Department of Mathematics, University of Texas at Austin

 

Collections: Mathematics; Geosciences