Divergence-free unfitted finite element discretisations for the Darcy problem
Fuente:
arXiv
Guardado en:
| Autores principales: | , , , , , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2026
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866912984964530176 |
|---|---|
| author | Badia, Santiago Boschman, Anne Martín, Alberto F. Nilsson, Erik Ruiz-Baier, Ricardo Zahedi, Sara |
| author_facet | Badia, Santiago Boschman, Anne Martín, Alberto F. Nilsson, Erik Ruiz-Baier, Ricardo Zahedi, Sara |
| contents | We develop an unfitted compatible finite element discretisation for the Darcy problem based on $H(\mathrm{div})$-conforming flux spaces and discontinuous pressure spaces. The method is designed to preserve pointwise discrete mass conservation while remaining robust in the presence of arbitrarily small cut cells arising from unfitted meshes. Robustness is achieved by combining an $L^2$-stabilisation of the flux with an additional mixed-term stabilisation that enhances pressure control without destroying the local conservation structure. We consider both cell-wise (bulk) and face-based ghost-penalty realisations of the stabilisation. Mixed boundary conditions are handled by weak imposition of both flux and pressure traces on unfitted boundaries. We prove stability and optimal-order a priori error estimates with constants independent of the cut configuration, and establish pressure-robust flux error bounds in the case of pure pressure boundary conditions. We also introduce an augmented Lagrangian variant that improves control of the conservation constraint and is amenable to efficient preconditioning strategies. Numerical experiments for a range of cut configurations, boundary-condition regimes and parameter choices confirm the theoretical results, demonstrating optimal convergence, cut-independent conditioning and mass conservation up to solver tolerance. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2603_26212 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Divergence-free unfitted finite element discretisations for the Darcy problem Badia, Santiago Boschman, Anne Martín, Alberto F. Nilsson, Erik Ruiz-Baier, Ricardo Zahedi, Sara Numerical Analysis 65N30, 65N12, 76S05 We develop an unfitted compatible finite element discretisation for the Darcy problem based on $H(\mathrm{div})$-conforming flux spaces and discontinuous pressure spaces. The method is designed to preserve pointwise discrete mass conservation while remaining robust in the presence of arbitrarily small cut cells arising from unfitted meshes. Robustness is achieved by combining an $L^2$-stabilisation of the flux with an additional mixed-term stabilisation that enhances pressure control without destroying the local conservation structure. We consider both cell-wise (bulk) and face-based ghost-penalty realisations of the stabilisation. Mixed boundary conditions are handled by weak imposition of both flux and pressure traces on unfitted boundaries. We prove stability and optimal-order a priori error estimates with constants independent of the cut configuration, and establish pressure-robust flux error bounds in the case of pure pressure boundary conditions. We also introduce an augmented Lagrangian variant that improves control of the conservation constraint and is amenable to efficient preconditioning strategies. Numerical experiments for a range of cut configurations, boundary-condition regimes and parameter choices confirm the theoretical results, demonstrating optimal convergence, cut-independent conditioning and mass conservation up to solver tolerance. |
| title | Divergence-free unfitted finite element discretisations for the Darcy problem |
| topic | Numerical Analysis 65N30, 65N12, 76S05 |
| url | https://arxiv.org/abs/2603.26212 |