Picture of neon light reading 'Open'

Discover open research at Strathprints as part of International Open Access Week!

23-29 October 2017 is International Open Access Week. The Strathprints institutional repository is a digital archive of Open Access research outputs, all produced by University of Strathclyde researchers.

Explore recent world leading Open Access research content this Open Access Week from across Strathclyde's many research active faculties: Engineering, Science, Humanities, Arts & Social Sciences and Strathclyde Business School.

Explore all Strathclyde Open Access research outputs...

A symmetric nodal conservative finite element method for the Darcy equation

Barrenechea, G. R. and Franca, L. P. and Valentin, F. (2009) A symmetric nodal conservative finite element method for the Darcy equation. SIAM Journal on Numerical Analysis, 47 (5). pp. 3652-3677. ISSN 0036-1429

[img]
Preview
Text (strathprints015098)
strathprints015098.pdf - Accepted Author Manuscript

Download (492kB) | Preview

Abstract

This work introduces and analyzes novel stable Petrov-Galerkin EnrichedMethods (PGEM) for the Darcy problem based on the simplest but unstable continuous P1/P0 pair. Stability is recovered inside a Petrov-Galerkin framework where element-wise dependent residual functions, named multi-scale functions, enrich both velocity and pressure trial spaces. Unlike the velocity test space that is augmented with bubble-like functions, multi-scale functions correct edge residuals as well. The multi-scale functions turn out to be the well-known lowest order Raviart-Thomas basis functions for the velocity and discontinuous quadratics polynomial functions for the pressure. The enrichment strategy suggests the way to recover the local mass conservation property for nodal-based interpolation spaces. We prove that the method and its symmetric version are well-posed and achieve optimal error estimates in natural norms. Numerical validations confirm claimed theoretical results.