BEM with linear complexity for the classical boundary integral operators

Alternative representations of boundary integral operators corresponding to elliptic boundary value problems are developed as a starting point for numerical approximations as, e.g., Galerkin boundary elements including numerical quadrature and panel-clustering. These representations have the advantage that the integrands of the integral operators have a reduced singular behaviour allowing one to choose the order of the numerical approximations much lower than for the classical formulations. Low-order discretisations for the single layer integral equations as well as for the classical double layer potential and the hypersingular integral equation are considered. We will present fully discrete Galerkin boundary element methods where the storage amount and the CPU time grow only linearly with respect to the number of unknowns.

     

FINITE ELEMENT NONLINEAR GALERKIN COUPLING METHOD FOR THE EXTERIOR STEADY NAVIER-STOKES PROBLEM

1.IntroductionNonlinearGalerkinmethodsaremultilevelschemesforthedissipativeevolutionpartialdifferentialequations.Theycorrespondtothesplittingsoftheunknownu:u=y z)wherethecomponentsareofdifferentorderofmagnitudewithrespecttoaparameterrelatedtothespati...     

Boundary integral approximation of a heat-diffusion problem in time-harmonic regime

In this paper we propose and analyse numerical methods for the approximation of the solution of Helmholtz transmission problems in the half plane. The problems we deal with arise from the study of some models in photothermal science. The solutions to the problem are represented as single layer potentials and an equivalent system of boundary integral equations is derived. We then give abstract necessary and sufficient conditions for convergence of Petrov–Galerkin discretizations of the boundary integral system and show for three different cases that these conditions are satisfied. We extend the results to other situations not related to thermal science and to non-smooth interfaces. Finally, we propose a simple full discretization of a Petrov–Galerkin scheme with periodic spline spaces and show some numerical experiments.     

A dual-dual mixed formulation for nonlinear exterior transmission problems

We combine a dual-mixed finite element method with a Dirichlet-to-Neumann mapping (derived by the boundary integral equation method) to study the solvability and Galerkin approximations of a class of exterior nonlinear transmission problems in the plane. As a model problem, we consider a nonlinear elliptic equation in divergence form coupled with the Laplace equation in an unbounded region of the plane. Our combined approach leads to what we call a dual-dual mixed variational formulation since the main operator involved has itself a dual-type structure. We establish existence and uniqueness of solution for the continuous and discrete formulations, and provide the corresponding error analysis by using Raviart-Thomas elements. The main tool of our analysis is given by a generalization of the usual Babuska-Brezzi theory to a class of nonlinear variational problems with constraints.

     

The coupling of boundary integral and finite element methods for nonmonotone nonlinear problems ∗

In this paper, we apply the coupling of the boundary integral and finite element methods to study the weak solvability of certain nonmonotone nonlinear exterior boundary value problems. In order to convert the original exterior problem into an equivalent nonlocal boundary value problem on a finite region, we employ two different approaches based on the use of one and two integral equations on the coupling boundary. Existence of a solution for the associated weak formulation, and convergence properties of the corresponding Galerkin approximations are deduced from fundamental results in nonlinear functional analysis. Indeed, the main arguments of our proofs are based on a surjectivity theorem for mappings of type (S) and on the Fredholm alternative for nonlinear A-proper mappings.     

On the convergence of the collocation method for nonlinear boundary integral equations

Recently, Galerkin and collocation methods have been analysed for some nonlinear boundary integral equations. For the collocation method it has been assumed that the nonlinearity is asymptotically linear. In this paper we remove this restriction. We shall prove the convergence of the collocation method for nonlinear boundary integral equations, when the nonlinearity has a polynomial growth condition. In addition to this the optimal order error estimates follow in L^q(Γ)-norm.     

Boundary integral equations for screen problems in IR3

Here we present a new solution procedure for Helmholtz and Laplacian Neumann screen or Dirichlet screen problems in IR³ via boundary integral equations of the first kind having as unknown the jump of the field or of its normal derivative, respectively, across the screen S. Under the assumption of local finite energy we show the equivalence of the integral equations and the original boundary value problems. Via the Wiener-Hopf method in the halfspace, localization and the calculus of pseudodifferential operators we derive existence, uniqueness and regularity results for the solution of our boundary integral equations together with its explicit behavior near the edge of the screen. We give Galerkin schemes based on our integral equations on S and obtain high convergence rates by using special singular elements besides regular splines as test and trial functions.     

An augmented galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problems

Here we present a new solution procedure for Helm-holtz and Laplacian Dirichlet screen and crack problems in IR² via boundary integral equations of the first kind having as an unknown the jump of the normal derivative across the screen or a crack curve T. Under the assumption of local finite energy we show the equivalence of the integral equations and the original boundary value problem. Via the method of local Mellin transform in [5]-[lo] and the calculus of pseudodifferential operators we derive existence, uniqueness and regularity results for the solution of our boundary integral equations together with its explicit behaviour near the screen or crack tips.With our integral equations we set up a Galerkin scheme on T and obtain high quasi-optimal convergence rates by using special singular elements besides regular splines as test and trial functions.     

New boundary element formulas for the biharmonic equation

We propose two new boundary integral equation formulas for the biharmonic equation with the Dirichlet boundary data that arises from plate bending problems in ℝ². Two boundary conditions, u and ∂u/∂n, usually yield a 2 × 2 non-symmetric matrix system of integral equations. Our new formulas yield scalar integral equations that can be handled more efficiently for theoretical and numerical purposes. In this paper we supply complete ellipticity and solvability analyses of our new formulas. Numerical experiments for simple Galerkin methods are also provided. This revised version was published online in June 2006 with corrections to the Cover Date.     

L p estimates of boundary integral equations for some nonlinear boundary value problems

Summary Recently, Galerkin and collocation methods have been analyzed for boundary integral equation formulations of some potential problems in the plane with nonlinear boundary conditions, and stability results and error estimates in theH ^1/2-norm have been proved (Ruotsalainen and Wendland, and Ruotsalainen and Saranen). We show that these results extend toL ^p setting without any extra conditions. These extensions are proved by studying the uniform boundedness of the inverses of the linearized integral operators, and then considering the nonlinear equations. The fact that inH ^1/2 setting the nonlinear operator is a homeomorphism with Lipschitz continuous inverse plays a crucial role. Optimal error estimates for the Galerkin and collocation method inL ^p space then follow.This research was performed while the second author was visiting professor at the University of Delaware, spring 1989     

A numerical method for the solution of multi-point problems for ordinary differential equations with integral constraints

We propose a method for the solution of a system of nonlinear ordinary differential equations with integral constraints, by transforming to multi-point boundary value problems. Examples are given.     

A boundary integral equation method for three-dimensional crack problems in elasticity

This paper presents a solution procedure for three-dimensional crack problems via first kind boundary integral equations on the crack surface. The Dirichlet (Neumann) problem is reduced to a system of integral equations for the jump of the traction (of the field) across the crack surface. The calculus of pseudodifferential operators is used to derive existence and regularity of the solutions of the integral equations. With the concept of the principal symbol and the Wiener-Hopf technique we derive the explicit behavior of the densities of the integral equations near the edge of the crack surface. Based on the detailed regularity results we show how to improve the boundary element Galerkin method for our integral equations. Quasi-optimal asymptotic estimates for the Galerkin error are given.     

A meshless local Galerkin method for the numerical solution of Hammerstein integral equations based on the moving least squares technique

In this paper, a computational scheme is proposed to estimate the solution of one- and two-dimensional Fredholm-Hammerstein integral equations of the second kind. The method approximates the solution using the discrete Galerkin method based on the moving least squares (MLS) approach as a locally weighted least squares polynomial fitting. The discrete Galerkin technique for integral equations results from the numerical integration of all integrals in the system corresponding to the Galerkin method. Since the proposed method is constructed on a set of scattered points, it does not require any background meshes and so we can call it as the meshless local discrete Galerkin method. The implication of the scheme for solving two-dimensional integral equations is independent of the geometry of the domain. The new method is simple, efficient and more flexible for most classes of nonlinear integral equations. The error analysis of the method is provided. The convergence accuracy of the new technique is tested over several Hammerstein integral equations and obtained results confirm the theoretical error estimates.     

Boundary element methods for Maxwell's equations on non-smooth domains

Summary. Variational boundary integral equations for Maxwell's equations on Lipschitz surfaces in are derived and their well-posedness in the appropriate trace spaces is established. An equivalent, stable mixed reformulation of the system of integral equations is obtained which admits discretization by Galerkin boundary elements based on standard spaces. On polyhedral surfaces, quasioptimal asymptotic convergence of these Galerkin boundary element methods is proved. A sharp regularity result for the surface multipliers on polyhedral boundaries with plane faces is established. Received January 5, 2001 / Revised version received August 6, 2001 / Published online December 18, 2001 Correspondence to: C. Schwab     

Galerkin methods for nonlinear Sobolev equations

Summary We study Galerkin approximations to the solution of nonlinear Sobolev equations with homogeneous Dirichlet boundary condition in two spatial dimensions and derive optimalL ² error estimates for the continuous Crank — Nicolson and Extrapolated Crank — Nicolson approximations.     

On the numerical analysis of finite element and dirichlet-to-neumann methods for nonlinear exterior transmission problems

We provide the numerical analysis of the combination of finite elements and Dirichlet-to-Neumann mappings (based on boundary integral operators) for a class of nonlinear exterior transmission problems whose weak formulations reduce to Lipschitz-continuous and strongly monotone operator equations. As a model we consider a nonlinear second order elliptic equation in divergence form in a bounded inner region of the plane, coupled with the Laplace equation in the corresponding unbounded exterior part. A discrete Galerkin scheme is presented by using linear finite elements on a triangulation of the domain, and then applying numerical quadrature and analytical formulae to evaluate all the linear, bilinear and semilinear forms involved. We prove the unique solvability of the discrete equations, and show the strong convergence of the approximate solutions. Furthermore, assuming additional regularity on the solution of the continuous operator equation, the asymptotic rate of convergence O(h) is also derived. Finally, numerical experiments are presented, which confirm the convergence results.     

Boundary element preconditioners for a hypersingular integral equation on an interval

We propose an almost optimal preconditioner for the iterative solution of the Galerkin equations arising from a hypersingular integral equation on an interval. This preconditioning technique, which is based on the single layer potential, was already studied for closed curves [11,14]. For a boundary element trial space, we show that the condition number is of order (1 + | log h _min|)², where h _min is the length of the smallest element. The proof requires only a mild assumption on the mesh, easily satisfied by adaptive refinement algorithms.     

Adaptive Galerkin boundary element methods with panel clustering

In this paper, we will propose a boundary element method for solving classical boundary integral equations on complicated surfaces which, possibly, contain a large number of geometric details or even uncertainties in the given data. The (small) size of such details is characterised by a small parameter and the regularity of the solution is expected to be low in such zones on the surface (which we call the wire-basket zones). We will propose the construction of an initial discretisation for such type of problems. Afterwards standard strategies for boundary element discretisations can be applied such as the h, p, and the adaptive hp-version in a straightforward way. For the classical boundary integral equations, we will prove the optimal approximation results of our so-called wire-basket boundary element method and discuss the stability aspects. Then, we construct the panel-clustering and -matrix approximations to the corresponding Galerkin BEM stiffness matrix. The method is shown to have an almost linear complexity with respect to the number of degrees of freedom located on the wire basket.     

Solvability of a nonlinear model Boltzmann equation in the problem of a plane shock wave

We consider a nonlinear system of integral equations describing the structure of a plane shock wave. Based on physical reasoning, we propose an iterative method for constructing an approximate solution of this system. The problem reduces to studying decoupled scalar nonlinear and linear integral equations for the gas temperature, density, and velocity. We formulate a theorem on the existence of a positive bounded solution of a nonlinear equation of the Uryson type. We also prove theorems on the existence and uniqueness of bounded positive solutions for linear integral equations in the space L ₁[?r, r] for all finite r < +∞. For a more general nonlinear integral equation, we prove a theorem on the existence of a positive solution and also find a lower bound and an integral upper bound for the constructed solution.     

On the boundary element method for some nonlinear boundary value problems

Summary Here we analyse the boundary element Galerkin method for two-dimensional nonlinear boundary value problems governed by the Laplacian in an interior (or exterior) domain and by highly nonlinear boundary conditions. The underlying boundary integral operator here can be decomposed into the sum of a monotoneous Hammerstein operator and a compact mapping. We show stability and convergence by using Leray-Schauder fixed-point arguments due to Petryshyn and Neas.Using properties of the linearised equations, we can also prove quasioptimal convergence of the spline Galerkin approximations.This work was carried out while the first author was visiting the University of Stuttgart