An interior–exterior Schwarz algorithm and its convergenceUn algorithme de Schwarz intérieur–extérieur et sa convergence

In this work we study the solution of Laplace's equation in a domain with holes by an iteration consisting of splitting the problem in an exterior one, around the holes, plus an interior problem in the unholed domain. We show the existence of a decomposition of the solution when the exterior problem is represented by means of a single-layer protential. Also, for the three-dimensional case and with some adjustments for the two-dimensional case, we prove convergence of the method by writing the iteration as a Jacobi iteration for an operator equation and studying the spectrum of the iteration operator. To cite this article: R. Celorrio et al., C. R. Acad. Sci. Paris, Ser. I 334 (2002) 923–926.     

A nodal coupling of finite and boundary elements

This article presents and analyzes a simple method for the exterior Laplace equation through the coupling of finite and boundary element methods. The main novelty is the use of a smooth parametric artificial boundary where boundary elements fit without effort together with a straight approximate triangulation in the bounded area, with the coupling done only in nodes. A numerically integrated version of the algorithm is also analyzed. Finally, an isoparametric variant with higher order is proposed. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 555–570, 2003     

Relaxing the hypotheses of Bielak–MacCamy’s BEM–FEM coupling

In this paper we show that the quasi-symmetric coupling of finite and boundary elements of Bielak and MacCamy can be freed of two very restricting hypotheses that appeared in the original paper: the coupling boundary can be taken polygonal/polyhedral and coupling can be done using the normal stress instead of the pseudostress. We will do this by first considering a model problem associated to the Yukawa equation, where we prove how compactness arguments can be avoided to show stability of Galerkin discretizations of a coupled system in the style of Bielak–MacCamy’s. We also show how discretization properties are robust in the continuation parameter that appears in the formulation. This analysis is carried out using a new and very simplified proof of the ellipticity of the Johnson–Nédélec BEM–FEM coupling operator. Finally, we show how to apply the techniques that we have fully developed in the model problem to the linear elasticity system.     

Asymptotic properties of some triangulations of the sphere

In this paper we analyse a method for triangulating the sphere originally proposed by Baumgardner and Frederickson in 1985. The method is essentially a refinement procedure for arbitrary spherical triangles that fit into a hemisphere. Refinement is carried out by dividing each triangle into four by introducing the midpoints of the edges as new vertices and connecting them in the usual ‘red’ way. We show that this process can be described by a sequence of piecewise smooth mappings from a reference triangle onto the spherical triangle. We then prove that the whole sequence of mappings is uniformly bi-Lipschitz and converges uniformly to a non-smooth parameterization of the spherical triangle, recovering the Baumgardner and Frederickson spherical barycentric coordinates. We also prove that the sequence of triangulations is quasi-uniform, that is, areas of triangles and lengths of the edges are roughly the same at each refinement level. Some numerical experiments confirm the theoretical results.     

Convergence analysis of a high-order Nyström integral-equation method for surface scattering problems

In this paper we present a convergence analysis for the Nyström method proposed in [J Comput Phys 169 (1):80–110, 2001] for the solution of the combined boundary integral equation formulations of sound-soft acoustic scattering problems in three-dimensional space. This fast and efficient scheme combines FFT techniques and a polar change of variables that cancels out the kernel singularity. We establish the stability of the algorithms in the $L^2$ norm and we derive convergence estimates in both the $L^2$ and $L^\infty $ norms. In particular, our analysis establishes theoretically the previously observed super-algebraic convergence of the method in cases in which the right-hand side is smooth.     

Full collocation methods for some boundary integral equations

In this paper we propose a fully discretized version of the collocation method applied to integral equations of the first kind with logarithmic kernel. After a stability and convergence analysis is given, we prove the existence of an asymptotic expansion of the error, which justifies the use of Richardson extrapolation. We further show how these expansions can be translated to a new expansion of potentials calculated with the numerical solution of a boundary integral equation such as those treated before. Some numerical experiments, confirming our theoretical results, are given. This revised version was published online in August 2006 with corrections to the Cover Date.     

Symmetric boundary integral formulations for Helmholtz transmission problems

In this work we propose and analyze numerical methods for the approximation of the solution of Helmholtz transmission problems in two or three dimensions. This kind of problems arises in many applications related to scattering of acoustic, thermal and electromagnetic waves. Formulations based on boundary integral methods are powerful tools to deal with transmission problems in unbounded media. Different formulations using boundary integral equations can be found in the literature. We propose here new symmetric formulations based on a paper by Martin Costabel and Ernst P. Stephan (1985), that uses the Calderón projector for the interior and exterior problems to develop closed expressions for the interior and exterior Dirichlet-to-Neumann operators. These operators are then matched to obtain an integral system that is equivalent to the Helmholtz transmission problem and uses Cauchy data on the transmission boundary as unknowns. We show how to simplify the aspect and analysis of the method by employing an additional mortar unknown with respect to the ones used in the original paper, writing it in an appropriate way to devise Krylov type iterations based on the separate Dirichlet-to-Neumann operators.     

Theoretical aspects of the application of convolution quadrature to scattering of acoustic waves

In this paper we address several theoretical questions related to the numerical approximation of the scattering of acoustic waves in two or three dimensions by penetrable non-homogeneous obstacles using convolution quadrature (CQ) techniques for the time variable and coupled boundary element method/finite element method for the space variable. The applicability of CQ to waves requires polynomial type bounds for operators related to the operator Δ − s ² in the right half complex plane. We propose a new systematic way of dealing with this problem, both at the continuous and semidiscrete-in-space cases. We apply the technique to three different situations: scattering by a group of sound-soft and -hard obstacles, by homogeneous and non-homogeneous obstacles.     

Crouzeix–Raviart boundary elements

This paper establishes a foundation of non-conforming boundary elements. We present a discrete weak formulation of hypersingular integral operator equations that uses Crouzeix–Raviart elements for the approximation. The cases of closed and open polyhedral surfaces are dealt with. We prove that, for shape regular elements, this non-conforming boundary element method converges and that the usual convergence rates of conforming elements are achieved. Key ingredient of the analysis is a discrete Poincaré–Friedrichs inequality in fractional order Sobolev spaces. A numerical experiment confirms the predicted convergence of Crouzeix–Raviart boundary elements. Norbert Heuer is supported by Fondecyt-Chile under grant no. 1080044. F.-J. Sayas is partially supported by MEC-FEDER Project MTM2007-63204 and Gobierno de Aragón (Grupo Consolidado PDIE).     

Analysis of a new BEM‐FEM coupling for two‐dimensional fluid‐solid interaction

In this work we present a new numerical method, based on a coupling of finite and boundary elements, to solve a fluid‐solid interaction problem in the plane. The discrete method uses classical Lagrange finite elements adapted to curved boundaries for the field variable and spectral approximation of the unknowns on the artificial boundary. We provide error estimates for this Galerkin scheme and propose a full discretization based on elementary quadrature formulae, showing that the perturbation due to numerical integration preserves the optimal rate of convergence. We also suggest an iterative method to solve the complicated linear systems arising from this type of schemes. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005