A convergent boundary integral method for three-dimensional water waves

We design a boundary integral method for time-dependent, three-dimensional, doubly periodic water waves and prove that it converges with accuracy, without restriction on amplitude. The moving surface is represented by grid points which are transported according to a computed velocity. An integral equation arising from potential theory is solved for the normal velocity. A new method is developed for the integration of singular integrals, in which the Green's function is regularized and an efficient local correction to the trapezoidal rule is computed. The sums replacing the singular integrals are treated as discrete versions of pseudodifferential operators and are shown to have mapping properties like the exact operators. The scheme is designed so that the error is governed by evolution equations which mimic the structure of the original problem, and in this way stability can be assured. The wavelike character of the exact equations of motion depends on the positivity of the operator which assigns to a function on the surface the normal derivative of its harmonic extension; similarly, the stability of the scheme depends on maintaining this property for the discrete operator. With grid points, the scheme can be implemented with essentially operations per time step.

     

Estimating quadrature errors for an efficient method for quasi-singular boundary integral

The numerical resolution of the boundary integral equations applied to the differential equations of Laplace, Helmholtz and Maxwell requires the handling of quasi-singular integrals with different order of singularity. The numerical approximation of the integral equations of different kinds is made by boundary finite elements. In this paper, we present a complete survey for estimating quadrature errors for the numerical techniques proposed by Huang and Cruse [Q. Huang, T.A. Cruse, Some notes on singular integral techniques in boundary element analysis, Int. J. Numer. Methods Eng. 36 (15) (1993) 2643-2659], to calculate the quasi-singular integrals. To validate the accuracy and efficiency of these techniques and approve our study some numerical examples are presented and discussed.     

On the extraction technique in boundary integral equations

In this paper we develop and analyze a bootstrapping algorithm for the extraction of potentials and arbitrary derivatives of the Cauchy data of regular three-dimensional second order elliptic boundary value problems in connection with corresponding boundary integral equations. The method rests on the derivatives of the generalized Green's representation formula, which are expressed in terms of singular boundary integrals as Hadamard's finite parts. Their regularization, together with asymptotic pseudohomogeneous kernel expansions, yields a constructive method for obtaining generalized jump relations. These expansions are obtained via composition of Taylor expansions of the local surface representation, the density functions, differential operators and the fundamental solution of the original problem, together with the use of local polar coordinates in the parameter domain. For boundary integral equations obtained by the direct method, this method allows the recursive numerical extraction of potentials and their derivatives near and up to the boundary surface.

     

A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

Modified boundary integral formulations for the Helmholtz equation

In this paper we describe some modified regularized boundary integral equations to solve the exterior boundary value problem for the Helmholtz equation with either Dirichlet or Neumann boundary conditions. We formulate combined boundary integral equations which are uniquely solvable for all wave numbers even for Lipschitz boundaries Γ=∂Ω. This approach extends and unifies existing regularized combined boundary integral formulations.     

A boundary integral equation method for the transmission eigenvalue problem

We propose a new integral equation formulation to characterize and compute transmission eigenvalues for constant refractive index that play an important role in inverse scattering problems for penetrable media. As opposed to the recently developed approach by Cossonnière and Haddar [1,2] which relies on a two by two system of boundary integral equations our analysis is based on only one integral equation in terms of Dirichlet-to-Neumann or Robin-to-Dirichlet operators which results in a noticeable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further we employ the numerical algorithm for analytic non-linear eigenvalue problems that was recently proposed by Beyn [3] for the numerical computation of transmission eigenvalues via this new integral equation.     

Block boundary value methods for solving Volterra integral and integro-differential equations

Reducible quadrature rules generated by boundary value methods are considered in block version and applied to solve the second kind Volterra integral equations and Volterra integro-differential equations. These extended block boundary value methods are shown to possess both excellent stability properties and high accuracy for Volterra-type equations. Numerical experiments are presented and the efficiency, accuracy and stability of the schemes are confirmed.     

A plane problem of uncoupled thermomagnetoelasticity for an infinite,elliptical cylinder carrying a steady axial current by a boundary integral method

A boundary integral method earlier proposed by two of the authors is used to solve a problem of uncoupled magnetothermoelasticity for an infinite, elliptical cylindrical conductor carrying a steady axial, uniform electric current. The cylinder is placed in a variable ambient temperature and is allowed to exchange heat with the surrounding medium.     

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.     

A note on a singular integral equation arising in water wave scattering

A quick method of solution of a singular integral equationof the first kind involving both logarithmic singularity aswell as Cauchy-type singularity is explained.     

Adaptive methods for boundary integral equations: Complexity and convergence estimates

This paper is concerned with developing numerical techniques for the adaptive application of global operators of potential type in wavelet coordinates. This is a core ingredient for a new type of adaptive solvers that has so far been explored primarily for PDEs. We shall show how to realize asymptotically optimal complexity in the present context of global operators. ``Asymptotically optimal' means here that any target accuracy can be achieved at a computational expense that stays proportional to the number of degrees of freedom (within the setting determined by an underlying wavelet basis) that would ideally be necessary for realizing that target accuracy if full knowledge about the unknown solution were given. The theoretical findings are supported and quantified by the first numerical experiments.

     

Boundary integral methods for singularly perturbed boundary value problems

In this paper we consider boundary integral methods appliedto boundary value problems for the positive definite Helmholtz-typeproblem –U + ²U = 0 in a bounded or unbounded domain,with the parameter real and possibly large. Applications arisein the implementation of space–time boundary integralmethods for the heat equation, where is proportional to 1/(t),and t is the time step. The corresponding layer potentials arisingfrom this problem depend nonlinearly on the parameter and havekernels which become highly peaked as , causing standard discretizationschemes to fail. We propose a new collocation method with arobust convergence rate as . Numerical experiments on a modelproblem verify the theoretical results.     

The boundary integral method for the linearized rotating Navier-Stokes equations in exterior domain

In this paper, we apply the boundary integral method to the linearized rotating Navier-Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and an infinite domain, we obtain a coupled problem by the linearized rotating Navier-Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence and uniqueness of solution. Finally, we study the finite element approximation for the coupled problem and obtain the error estimate between the solution of the coupled problem and its approximation solution.     

Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.     

The method of boundary integral equations for a mixed scattering problem

In this paper we consider the scattering of an electromagnetic time-harmonic plane wave by an infinite cylinder having an open arc and a bounded domain in R² as cross section. To this end, we solve a scattering problem for the Helmholtz equation in R² where the scattering object is a combination of a crack Γ and a bounded obstacle D, and we have Dirichlet-impedance type boundary condition on Γ and Dirichlet boundary condition on ∂D (∂D∈C²). Applying potential theory, the problem can be reformulated as a boundary integral system. We establish the existence and uniqueness of a solution to the system by using the Fredholm theory.     

A new analytical technique to solve Volterra's integral equations

The purpose of the present paper is to show that the well‐known homotopy analysis method for solving ordinary and partial differential equations can be applied to solve linear and nonlinear integral equations of Volterra's type with high accuracy as well. Comparison of the present method with Adomian decomposition method (ADM), a well‐known method to solve integral equations, reveals that the ADM is only especial case of the present method. Furthermore, some illustrating examples such as linear, nonlinear and singular integral equations of Volterra's type are given to show high efficiency with reliable accuracy of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

An extrapolation method for a class of boundary integral equations

Boundary value problems of the third kind are converted into boundary integral equations of the second kind with periodic logarithmic kernels by using Green's formulas. For solving the induced boundary integral equations, a Nyström scheme and its extrapolation method are derived for periodic Fredholm integral equations of the second kind with logarithmic singularity. Asymptotic expansions for the approximate solutions obtained by the Nyström scheme are developed to analyze the extrapolation method. Some computational aspects of the methods are considered, and two numerical examples are given to illustrate the acceleration of convergence.

     

On the numerical solution of a boundary integral equation for the exterior Neumann problem on domains with corners

The authors propose a “modified” Nyström method to approximate the solution of a boundary integral equation connected with the exterior Neumann problem for Laplace's equation on planar domains with corners. They prove the convergence and the stability of the method and show some numerical tests.