Picture of smart phone in human hand

World leading smartphone and mobile technology research at Strathclyde...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by University of Strathclyde researchers, including by Strathclyde researchers from the Department of Computer & Information Sciences involved in researching exciting new applications for mobile and smartphone technology. But the transformative application of mobile technologies is also the focus of research within disciplines as diverse as Electronic & Electrical Engineering, Marketing, Human Resource Management and Biomedical Enginering, among others.

Explore Strathclyde's Open Access research on smartphone technology now...

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.