Imposing nonlocal boundary conditions in Galerkin-type methods based on non-interpolatory functions

Svajūnas Sajavičius,Thomas Takacs
DOI: https://doi.org/10.1016/j.camwa.2020.09.016
2023-08-01
Abstract:The imposition of inhomogeneous Dirichlet (essential) boundary conditions is a fundamental challenge in the application of Galerkin-type methods based on non-interpolatory functions, i.e., functions which do not possess the Kronecker delta property. Such functions typically are used in various meshfree methods, as well as methods based on the isogeometric paradigm. The present paper analyses a model problem consisting of the Poisson equation subject to non-standard boundary conditions. Namely, instead of classical boundary conditions, the model problem involves Dirichlet- and Neumann-type nonlocal boundary conditions. Variational formulations with strongly and weakly imposed inhomogeneous Dirichlet-type nonlocal conditions are derived and compared within an extensive numerical study in the isogeometric framework based on non-uniform rational B-splines (NURBS). The attention in the numerical study is paid mainly to the influence of the nonlocal boundary conditions on the properties of the considered discretisation methods.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the fundamental challenge of imposing non - homogeneous Dirichlet (essential) boundary conditions in Galerkin - type methods using non - interpolating functions (i.e., functions without the Kronecker delta property). Such functions are typically used in various meshless methods and methods based on isogeometric analysis. Specifically, the paper focuses on methods of imposing non - standard boundary conditions, especially non - local Dirichlet and Neumann - type boundary conditions, in the Poisson equation model problem. ### Background of the Paper 1. **Non - interpolating Functions**: The basis functions used in meshless methods and isogeometric analysis usually do not have the interpolation property, which means that these basis functions are not necessarily zero at other nodes. This leads to difficulties in imposing non - homogeneous Dirichlet boundary conditions in these methods. 2. **Existing Techniques**: There are currently many techniques that can be used to impose classical essential boundary conditions in meshless methods, such as the Lagrange multiplier method, modified variational principles, penalty function methods, perturbed Lagrange methods, modified meshless basis functions, and coupling meshless methods with grid - based techniques (such as the finite element method). ### Research Objectives The main objective of the paper is to study how to effectively impose non - local boundary conditions in Galerkin methods based on non - interpolating basis functions and analyze the properties of these methods. Specifically: - **Model Problem**: The paper considers a Poisson equation model problem with non - standard boundary conditions. - **Variational Form**: Derived and compared the variational forms of imposing and weakly imposing non - homogeneous Dirichlet - type non - local boundary conditions. - **Numerical Study**: Conducted extensive numerical studies based on non - uniform rational B - splines (NURBS) within the framework of isogeometric analysis, with a focus on the influence of non - local boundary conditions on the properties of discretization methods. ### Main Contributions 1. **Non - local Boundary Conditions**: Proposed new methods for dealing with non - local boundary conditions, which can express non - local couplings between different parts of the model geometry. 2. **Variational Forms**: Presented several variational forms, including methods for imposing and weakly imposing non - homogeneous Dirichlet - type non - local boundary conditions. 3. **Numerical Experiments**: Verified the effectiveness of these methods through numerical experiments and compared the accuracy and condition numbers of different methods. ### Conclusions Through theoretical analysis and numerical experiments, the paper demonstrates effective methods for imposing non - local boundary conditions in Galerkin methods based on non - interpolating basis functions and points out the potential advantages and limitations of these methods in practical applications. In particular, the Nitsche method has better accuracy and stability than the penalty function method in some cases.