TY - JOUR
T1 - Robust high-order unfitted finite elements by interpolation-based discrete extension
AU - Badia, Santiago
AU - Neiva, Eric
AU - Verdugo, Francesc
N1 - Funding Information:
This research was partially funded by the Australian Government through the Australian Research Council (project number DP210103092 ), the European Commission under the FET-HPC ExaQUte project (Grant agreement ID: 800898 ) within the Horizon 2020 Framework Programme and the project RTI2018-096898-B-I00 from the “ FEDER/Ministerio de Ciencia e Innovación (MCIN) – Agencia Estatal de Investigación (AEI) ”. F. Verdugo acknowledges support from the “ Severo Ochoa Program for Centers of Excellence in R&D (2019-2023) ” under the grant CEX2018-000797-S funded by MCIN/AEI/10.13039/501100011033. This work was also supported by computational resources provided by the Australian Government through NCI under the National Computational Merit Allocation Scheme.
Publisher Copyright:
© 2022 Elsevier Ltd
PY - 2022/12/1
Y1 - 2022/12/1
N2 - In this work, we propose a novel formulation for the solution of partial differential equations using finite element methods on unfitted meshes. The proposed formulation relies on the discrete extension operator proposed in the aggregated finite element method. This formulation is robust with respect to the location of the boundary/interface within the cell. One can prove enhanced stability results, not only on the physical domain, but on the whole active mesh. However, the stability constants grow exponentially with the polynomial order being used, since the underlying extension operators are defined via extrapolation. To address this issue, we introduce a new variant of aggregated finite elements, in which the extension in the physical domain is an interpolation for polynomials of order higher than two. As a result, the stability constants only grow at a polynomial rate with the order of approximation. We demonstrate that this approach enables robust high-order approximations with the aggregated finite element method. The proposed method is consistent, optimally convergent, and with a condition number that scales optimally for high order approximation.
AB - In this work, we propose a novel formulation for the solution of partial differential equations using finite element methods on unfitted meshes. The proposed formulation relies on the discrete extension operator proposed in the aggregated finite element method. This formulation is robust with respect to the location of the boundary/interface within the cell. One can prove enhanced stability results, not only on the physical domain, but on the whole active mesh. However, the stability constants grow exponentially with the polynomial order being used, since the underlying extension operators are defined via extrapolation. To address this issue, we introduce a new variant of aggregated finite elements, in which the extension in the physical domain is an interpolation for polynomials of order higher than two. As a result, the stability constants only grow at a polynomial rate with the order of approximation. We demonstrate that this approach enables robust high-order approximations with the aggregated finite element method. The proposed method is consistent, optimally convergent, and with a condition number that scales optimally for high order approximation.
KW - Aggregated finite elements
KW - Embedded methods
KW - High-order finite elements
KW - Immersed methods
KW - Unfitted finite elements
UR - http://www.scopus.com/inward/record.url?scp=85139362091&partnerID=8YFLogxK
U2 - 10.1016/j.camwa.2022.09.027
DO - 10.1016/j.camwa.2022.09.027
M3 - Article
AN - SCOPUS:85139362091
VL - 127
SP - 105
EP - 126
JO - Computers and Mathematics with Applications
JF - Computers and Mathematics with Applications
SN - 0898-1221
ER -