A Non-intrusive Approach for the Imposition of Strong Dirichlet Boundary Conditions in Unfitted Boundary Meshes
Abstract
The enforcement of essential boundary conditions is a fundamental challenge in unfitted boundary methods. This paper presents a non-intrusive, black-box strategy for imposing such conditions in unfitted meshes. The approach is intended for situations where the user does not have access to the solver's source code or its mathematical formulation, which is often the case when using commercial software. The proposed algorithm allows solvers originally designed for body-fitted meshes to be used in unfitted cases, provided that four conditions are satisfied: (i) the solver must support user customization by means of scripting, (ii) allow the imposition of Dirichlet boundary conditions at the node level through scripting, (iii) permit the deactivation of elements outside the physical domain, and (iv) provide access to the solution gradient within active elements. The last condition can also be satisfied by externally reconstructing the gradient from nodal values and connectivity information, provided the element formulation is known, making it optional in practice. These requirements are very fair demands and are satisfied by the vast majority of production-ready, possibly commercial, codes. In the current work, we show the application of this non-intrusive algorithm in the context of the Finite Element Method (FEM) and Isogeometric Analysis (IGA) discretizations, demonstrating optimal L2-norm error convergence. This is demonstrated using the Kratos Multiphysics code (release \texttt{v10.1}) \emph{from the user API, simply leveraging the capabilities mentioned above.}
Versions
Related
- A review of weakly enforced Dirichlet boundary conditions in computational flow analysis
- Residual-Based Adaptive Refinement in Hybrid Variational-Collocation PINNs for Singular Elliptic PDEs with Mixed Dirichlet-Neumann Boundary Conditions
- Generalization of Rayleigh's high-frequency theory for the 2D Helmholtz equation in a half-space subject to a radiation condition at infinity and a Dirichlet condition on a 1D periodically-uneven boundary
- Boundary Conditions and Entanglement in Anti-de Sitter Space
- The Arrow of Records: Structured Boundary Conditions, Receptive Capacity, and Why History Accumulates Toward the Future
BibTeX
@misc{camarotti2026nonintrusive,
title = {{A Non-intrusive Approach for the Imposition of Strong Dirichlet Boundary Conditions in Unfitted Boundary Meshes}},
author = {J. I. Camarotti and R. Aristio and R. Rossi and R. Zorrilla and R. Wüchner},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.15398},
eprint = {2609.15398},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.15398},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - A Non-intrusive Approach for the Imposition of Strong Dirichlet Boundary Conditions in Unfitted Boundary Meshes
AU - J. I. Camarotti
AU - R. Aristio
AU - R. Rossi
AU - R. Zorrilla
AU - R. Wüchner
PY - 2026
DA - 2026/09/14
PB - arXiv
DO - 10.48550/arxiv.2609.15398
UR - https://arxiv.org/abs/2609.15398
AB - The enforcement of essential boundary conditions is a fundamental challenge in unfitted boundary methods. This paper presents a non-intrusive, black-box strategy for imposing such conditions in unfitted meshes. The approach is intended for situations where the user does not have access to the solver's source code or its mathematical formulation, which is often the case when using commercial software. The proposed algorithm allows solvers originally designed for body-fitted meshes to be used in unfitted cases, provided that four conditions are satisfied: (i) the solver must support user customization by means of scripting, (ii) allow the imposition of Dirichlet boundary conditions at the node level through scripting, (iii) permit the deactivation of elements outside the physical domain, and (iv) provide access to the solution gradient within active elements. The last condition can also be satisfied by externally reconstructing the gradient from nodal values and connectivity information, provided the element formulation is known, making it optional in practice. These requirements are very fair demands and are satisfied by the vast majority of production-ready, possibly commercial, codes. In the current work, we show the application of this non-intrusive algorithm in the context of the Finite Element Method (FEM) and Isogeometric Analysis (IGA) discretizations, demonstrating optimal L^2-norm error convergence. This is demonstrated using the Kratos Multiphysics code (release \texttt{v10.1}) \emph{from the user API, simply leveraging the capabilities mentioned above.}
ER -