Related papers: Corrected Trapezoidal Rules for Boundary Integral …
We consider the numerical algorithm for the two-dimensional time-harmonic elastic wave scattering by unbounded rough surfaces with Dirichlet boundary condition. A Nystr\"{o}m method is proposed for the scattering problem based on the…
We establish improved convergence rates for curved boundary element methods applied to the three-dimensional (3D) Laplace and Helmholtz equations with smooth geometry and data. Our analysis relies on a precise analysis of the consistency…
A method for computing singular or nearly singular integrals on closed surfaces was presented by J. T. Beale, W. Ying, and J. R. Wilson [Comm. Comput. Phys. 20 (2016), 733--753, arXiv:1508.00265] and applied to single and double layer…
In this paper, a highly efficient fast boundary element method (BEM) for solving large-scale engineering acoustic problems in a broad frequency range is developed and implemented. The acoustic problems are modeled by the Burton-Miller…
A Nystrom-based high-order (HO) discretization scheme for surface integral equations (SIEs) for analyzing the electroencephalography (EEG) forward problem is proposed in this work. We use HO surface elements and interpolation functions for…
Peskin's Immersed Boundary (IB) model and method are among one of the most important modeling tools and numerical methods. The IB method has been known to be first order accurate in the velocity. However, almost no rigorous theoretical…
Fully symmetric positive interior (f-SPI) quadrature rules are key building blocks for high-order discretizations of partial differential equations, yet high-degree rules with few nodes remain scarce on reference elements commonly used in…
In this paper, we develop highly accurate Nystr\"{o}m methods for the volume integral equation (VIE) of the Maxwell equation for 3-D scatters. The method is based on a formulation of the VIE equation where the Cauchy principal value of the…
The Boundary Element Method (BEM) is a powerful numerical approach for solving 3D elastostatic problems, particularly useful for crack propagation in fracture mechanics and half-space problems. A key challenge in BEM lies in handling…
We investigate the behavior of integral formulations of variable coefficient elliptic partial differential equations (PDEs) in the presence of steep internal layers. In one dimension, the equations that arise can be solved analytically and…
This paper presents an efficient high-order sharp-interface method for solving the three-dimensional (3D) Poisson equation with Dirichlet boundary conditions on a nonuniform Cartesian grid with irregular domain boundaries. The new approach…
We present a simple and effective method for evaluating double-and single-layer potentials for Laplace's equation in three dimensions close to the boundary. The close evaluation of these layer potentials is challenging because they are…
In this contribution we develop a cut finite element method with boundary value correction of the type originally proposed by Bramble, Dupont, and Thomee. The cut finite element method is a fictitious domain method with Nitsche type…
This paper introduces a novel class of indirect boundary integral equation (BIE) formulations for the solution of electromagnetic scattering problems involving smooth perfectly electric conductors (PECs) in three-dimensions. These…
In this paper, we propose new linearly convergent second-order methods for minimizing convex quartic polynomials. This framework is applied for designing optimization schemes, which can solve general convex problems satisfying a new…
By employing non-equispaced grid points near boundaries, boundary-optimized upwind finite-difference operators of orders up to nine are developed. The boundary closures are constructed within a diagonal-norm summation-by-parts (SBP)…
We propose a boundary value correction approach for cases when curved boundaries are approximated by straight lines (planes) and Lagrange multipliers are used to enforce Dirichlet boundary conditions. The approach allows for optimal order…
Over the past years, inverse problems in partial differential equations have garnered increasing interest among scientists and engineers. However, due to the lack of conventional stability, nonlinearity and non-convexity, these problems are…
Multidimensional diagonal-norm summation-by-parts (SBP) operators with collocated volume and facet nodes, known as diagonal-$ \mathsf{E} $ operators, are attractive for entropy-stable discretizations from an efficiency standpoint. However,…
In this paper, we propose a first second-order scheme based on arbitrary non-Euclidean norms, incorporated by Bregman distances. They are introduced directly in the Newton iterate with regularization parameter proportional to the square…