Burman, E;
Hansbo, P;
Larson, MG;
(2018)
A simple finite element method for elliptic bulk problems with embedded surfaces.
Computational Geosciences
10.1007/s10596-018-9792-y.
(In press).
Preview |
Text
Burman 2018_Article_ASimpleFiniteElementMethodForE.pdf - Published Version Download (1MB) | Preview |
Abstract
In this paper we develop a simple finite element method for simulation of embedded layers of high permeability in a matrix of lower permeability using a basic model of Darcy flow in embedded cracks. The cracks are allowed to cut through the mesh in arbitrary fashion and we take the flow in the crack into account by superposition. The fact that we use continuous elements leads to suboptimal convergence due to the loss of regularity across the crack. We therefore refine the mesh in the vicinity of the crack in order to recover optimal order convergence in terms of the global mesh parameter. The proper degree of refinement is determined based on an a priori error estimate and can thus be performed before the actual finite element computation is started. Numerical examples showing this effect and confirming the theoretical results are provided. The approach is easy to implement and beneficial for rapid assessment of the effect of crack orientation and may for example be used in an optimization loop.
Type: | Article |
---|---|
Title: | A simple finite element method for elliptic bulk problems with embedded surfaces |
Open access status: | An open access version is available from UCL Discovery |
DOI: | 10.1007/s10596-018-9792-y |
Publisher version: | https://doi.org/10.1007/s10596-018-9792-y |
Language: | English |
Additional information: | This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. |
Keywords: | Darcy equation, Fracture, Embedded layer, Cut finite element methods |
UCL classification: | UCL UCL > Provost and Vice Provost Offices UCL > Provost and Vice Provost Offices > UCL BEAMS UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences > Dept of Mathematics |
URI: | https://discovery.ucl.ac.uk/id/eprint/1574161 |
Archive Staff Only
View Item |