Details
| Original language | English |
|---|---|
| Publication status | E-pub ahead of print - 4 Jun 2018 |
Abstract
Keywords
- math.NA, 65M60
Cite this
- Standard
- Harvard
- Apa
- Vancouver
- BibTeX
- RIS
2018.
Research output: Working paper/Preprint › Preprint
}
TY - UNPB
T1 - On the implementation of a locally modified finite element method for interface problems in deal.II
AU - Frei, Stefan
AU - Richter, Thomas
AU - Wick, Thomas
PY - 2018/6/4
Y1 - 2018/6/4
N2 - In this work, we describe a simple finite element approach that is able to resolve weak discontinuities in interface problems accurately. The approach is based on a fixed patch mesh consisting of quadrilaterals, that will stay unchanged independent of the position of the interface. Inside the patches we refine once more, either in eight triangles or in four quadrilaterals, in such a way that the interface is locally resolved. The resulting finite element approach can be considered a fitted finite element approach. In our practical implementation, we do not construct this fitted mesh, however. Instead, the local degrees of freedom are included in a parametric way in the finite element space, or to be more precise in the local mappings between a reference patch and the physical patches. We describe the implementation in the open source C++ finite element library deal.II in detail and present two numerical examples to illustrate the performance of the approach. Finally, detailed studies of the behavior of iterative linear solvers complement this work.
AB - In this work, we describe a simple finite element approach that is able to resolve weak discontinuities in interface problems accurately. The approach is based on a fixed patch mesh consisting of quadrilaterals, that will stay unchanged independent of the position of the interface. Inside the patches we refine once more, either in eight triangles or in four quadrilaterals, in such a way that the interface is locally resolved. The resulting finite element approach can be considered a fitted finite element approach. In our practical implementation, we do not construct this fitted mesh, however. Instead, the local degrees of freedom are included in a parametric way in the finite element space, or to be more precise in the local mappings between a reference patch and the physical patches. We describe the implementation in the open source C++ finite element library deal.II in detail and present two numerical examples to illustrate the performance of the approach. Finally, detailed studies of the behavior of iterative linear solvers complement this work.
KW - math.NA
KW - 65M60
M3 - Preprint
BT - On the implementation of a locally modified finite element method for interface problems in deal.II
ER -