What problem does this paper attempt to address?

The problem that this paper attempts to solve is to achieve a stable MATLAB implementation of the Isogeometric Boundary Element Method (IGABEM) in two - dimensional space. Specifically, the paper focuses on how to reduce the degrees of freedom while maintaining or improving the accuracy of the numerical solution by using smooth splines (i.e., NURBS). In addition, the paper also explores how to optimize the convergence behavior under high - order discretization through a posteriori error estimation and adaptive strategies, especially when the given data and the unknown solution are not smooth. These techniques are particularly useful for dealing with boundary integral equations of partial differential equations (PDEs), especially when dealing with singularities and complex geometries. ### Main contributions of the paper 1. **Implementation of adaptive IGABEM**: - An easy - to - use MATLAB package IGABEM2D is provided, which implements the adaptive Galerkin boundary element method. - This method is carried out in the isogeometric analysis framework, using NURBS as basis functions to improve the accuracy and efficiency of the numerical solution. 2. **A posteriori error estimation**: - Several different nodal - basis error estimators are proposed, including the (h - h/2) estimator, the Faermann estimator, and the weighted residual estimator. - These estimators are used to guide the adaptive mesh refinement process to optimize the convergence rate. 3. **Adaptive algorithm**: - An adaptive IGABEM algorithm is designed. By marking the nodes with large errors and performing mesh refinement, the accuracy of the solution is gradually improved. - The algorithm takes into account the maximum mesh ratio and multiplicity increase to control the mesh quality in the parameter domain. 4. **Numerical experiments**: - Through numerical experiments on different geometries (such as pacman geometry and crack geometry), the effectiveness and superiority of the adaptive IGABEM are verified. - The experimental results show that the adaptive method can significantly improve the convergence rate, especially when dealing with problems with singularities. ### Mathematical formulas - **Single - layer operator \( V \)**: \[ [V\varphi](x) := -\frac{1}{2\pi} \int_{\Gamma} \log|x - y| \varphi(y)\, dy \] - **Hypersingular operator \( W \)**: \[ [Wu](x) := \frac{1}{2\pi} \frac{\partial}{\partial \nu(x)} \int_{\Gamma} \frac{\partial}{\partial \nu(y)} \log|x - y| u(y)\, dy \] - **Double - layer operator \( K \)**: \[ [Ku](x) := -\frac{1}{2\pi} \int_{\Gamma} \frac{\partial}{\partial \nu(y)} \log|x - y| u(y)\, dy \] - **Adjoint double - layer operator \( K' \)**: \[ [K'\varphi](x) := -\frac{1}{2\pi} \int_{\Gamma} \varphi(y) \frac{\partial}{\partial \nu(x)} \log|x - y|\, dy \] - **Weakly singular integral equation**: \[ [V\varphi](x) = f(x)\quad \text{for }x\in \Gamma \] - **Hypersingular integral equation**: \[ [Wu](x) = g(x)\quad \text{for }x\in \Gamma \] ### Conclusion This paper demonstrates the effectiveness and superiority of the adaptive isogeometric boundary element method in dealing with complex geometries and singular problems through providing a stable MATLAB implementation. Through a posteriori error estimation and adaptive strategies, this method can maintain or improve the accuracy of the numerical solution while reducing the degrees of freedom, thus providing a powerful tool for practical applications.

Stable implementation of adaptive IGABEM in 2D in MATLAB

Implementation of Isogeometric Fast Multipole Boundary Element Methods for 2D Half-Space Acoustic Scattering Problems with Absorbing Boundary Condition

Isogeometric Boundary Element Analysis for 2D Transient Heat Conduction Problem with Radial Integration Method

Weight-adjusted discontinuous Galerkin methods: matrix-valued weights and elastic wave propagation in heterogeneous media

IB2d Reloaded: a more powerful Python and MATLAB implementation of the immersed boundary method

Recent Advances of Isogeometric Analysis in Computational Electromagnetics

Adaptive isogeometric analysis on two-dimensional trimmed domains based on a hierarchical approach

New Perspective to Isogeometric Analysis: Solving Isogeometric Analysis Problem by Fitting Load Function

An Efficient Implementation of Brezzi-Douglas-Marini (BDM) Mixed Finite Element Method in MATLAB

Isogeometric Boundary Element Method for Two-Dimensional Steady-State Non-Homogeneous Heat Conduction Problem

Discontinuous Galerkin Methods for Isogeometric Analysis for Elliptic Equations on Surfaces

Multi-resolution Isogeometric Analysis -- Efficient adaptivity utilizing the multi-patch structure

Multi-Adaptive Galerkin Methods for ODEs I

An efficient mass lumping scheme for isogeometric analysis based on approximate dual basis functions

ESPFEM2D: A MATLAB 2D explicit smoothed particle finite element method code for geotechnical large deformation analysis

Structural Shape Optimization Based on Multi-Patch Weakly Singular IGABEM and Particle Swarm Optimization Algorithm in Two-Dimensional Elastostatics

Adaptive Isogeometric analysis for transient dynamics: Space–time refinement based on hierarchical a-posteriori error estimations

FESTUNG: A MATLAB /GNU Octave toolbox for the discontinuous Galerkin method. Part IV: Generic problem framework and model-coupling interface

Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements

Structural shape optimization of three dimensional acoustic problems with isogeometric boundary element methods

An Application of Fast Multipole Method to Isogeometric Boundary Element Method in Two-Dimensional Acoustic Problems