Gespeichert in:
| Hauptverfasser: | , , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2025
|
| Schlagworte: | |
| Online-Zugang: | https://arxiv.org/abs/2510.22796 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866915579375386624 |
|---|---|
| author | Zhao, Bingying Song, Yin Deng, Quanling Li, Xin |
| author_facet | Zhao, Bingying Song, Yin Deng, Quanling Li, Xin |
| contents | The Generalized Finite Element Method (GFEM) is an effective unfitted numerical method for handling interface problems. By augmenting the standard FEM space with an appropriate enrichment space, GFEM can accurately capture C^0 solutions across the interfaces. While numerous GFEMs for interface problems have been studied, establishing a stable high-order GFEM with optimal convergence rates and robust system conditioning remains a challenge. The highest known order of two was established by Zhang and Babuška (SGFEM2, Comput. Methods Appl. Mech. Engrg. 363 (2020), 112889). In this paper, we propose a unified enrichment space construction and establish arbitrary high-order stable GFEMs (HoSGFEM) for elliptic interface problems. The main idea distinguishes itself from Zhang and Babuška's SGFEM2 substantially and it is twofold: a) we construct dimensionality-reduced auxiliary locally supported piecewise polynomials that satisfy the partition of unity property for elements containing interfaces; b) we construct the enrichment scheme based on d{1,(x-x_c^e),...,(y-y_c^e)^{p-1}} (d is the distance function; (x_c^e, y_c^e) is the center of the element containing interface, thus element-based) for arbitrary p-th order elements instead of d, d{1,x,y} or d{1,x,y,x^2,xy,y^2} (global functions) for p=1,2 in the literature. This idea results in an enrichment space that has a large angle with the standard FEM space, leading to the stability of the method with system condition number growing in order O(h^{-2}). We establish optimal convergence rates for HoSGFEM solutions under the proposed construction. Various numerical experiments with both straight and curved interfaces demonstrate the optimal convergence, FEM-comparable system condition number with O(h^{-2}) growth, and robustness as element boundaries approach interfaces. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2510_22796 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | HoSGFEM: High-order stable generalized finite element method for elliptic interface problem Zhao, Bingying Song, Yin Deng, Quanling Li, Xin Numerical Analysis The Generalized Finite Element Method (GFEM) is an effective unfitted numerical method for handling interface problems. By augmenting the standard FEM space with an appropriate enrichment space, GFEM can accurately capture C^0 solutions across the interfaces. While numerous GFEMs for interface problems have been studied, establishing a stable high-order GFEM with optimal convergence rates and robust system conditioning remains a challenge. The highest known order of two was established by Zhang and Babuška (SGFEM2, Comput. Methods Appl. Mech. Engrg. 363 (2020), 112889). In this paper, we propose a unified enrichment space construction and establish arbitrary high-order stable GFEMs (HoSGFEM) for elliptic interface problems. The main idea distinguishes itself from Zhang and Babuška's SGFEM2 substantially and it is twofold: a) we construct dimensionality-reduced auxiliary locally supported piecewise polynomials that satisfy the partition of unity property for elements containing interfaces; b) we construct the enrichment scheme based on d{1,(x-x_c^e),...,(y-y_c^e)^{p-1}} (d is the distance function; (x_c^e, y_c^e) is the center of the element containing interface, thus element-based) for arbitrary p-th order elements instead of d, d{1,x,y} or d{1,x,y,x^2,xy,y^2} (global functions) for p=1,2 in the literature. This idea results in an enrichment space that has a large angle with the standard FEM space, leading to the stability of the method with system condition number growing in order O(h^{-2}). We establish optimal convergence rates for HoSGFEM solutions under the proposed construction. Various numerical experiments with both straight and curved interfaces demonstrate the optimal convergence, FEM-comparable system condition number with O(h^{-2}) growth, and robustness as element boundaries approach interfaces. |
| title | HoSGFEM: High-order stable generalized finite element method for elliptic interface problem |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2510.22796 |