Boundary integral equations for bending of thermoelastic plates with transmission boundary conditions

The solution of an initial‐boundary value problem for bending of a piecewise‐homogeneous thermoelastic plate with transverse shear deformation is represented as various combinations of single‐layer and double‐layer time‐dependent potentials. The unique solvability of the boundary integral equations generated by these representations is proved in spaces of distributions. Copyright © 2009 John Wiley & Sons, Ltd.     

Time-dependent boundary integral equations for multiply connected plates

The existence of distributional solutions is discussed for theinitial-boundary value problems associated with the motion ofa thin, elastic, multiply connected plate, and for the boundaryequations arising from integral representations of such solutions.     

A new system of boundary integral equations for plates with free edges

We present a new system of boundary integral equations for the free-edge plate. This system expresses an abstract symmetric variational problem posed on the boundary of the plate. In order to obtain this abstract problem, we must set the exterior boundary value problem corresponding to the free-edge plate in a framework of weighted Sobolev spaces. Finally, we take care of the hypersingular kernels appearing in our system of BIEs, by using an abstract technique of integration by parts.     

Combined boundary integral equations for acoustic scattering‐resonance problems

In this paper, boundary integral formulations for a time‐harmonic acoustic scattering‐resonance problem are analyzed. The eigenvalues of eigenvalue problems resulting from boundary integral formulations for scattering‐resonance problems split in general into two parts. One part consists of scattering‐resonances, and the other one corresponds to eigenvalues of some Laplacian eigenvalue problem for the interior of the scatterer. The proposed combined boundary integral formulations enable a better separation of the unwanted spectrum from the scattering‐resonances, which allows in practical computations a reliable and simple identification of the scattering‐resonances in particular for non‐convex domains. The convergence of conforming Galerkin boundary element approximations for the combined boundary integral formulations of the resonance problem is shown in canonical trace spaces. Numerical experiments confirm the theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

Weak solutions for time‐dependent boundary integral equations associated with the bending of elastic plates under combined boundary data

The existence, uniqueness, stability, and integral representation of distributional solutions are investigated for the equations of motion of a thin elastic plate with a combination of displacement and moment‐stress components prescribed on the boundary. Copyright © 2004 John Wiley & Sons, Ltd.     

伴有边界摄动的Volterra型积分微分方程组的奇摄动

讨论了伴有边界摄动的二阶非线性Volterra型积分微分方程组的奇摄动.在适当的条件下,利用对角化技巧证明了解的存在性,构造出解的渐近展开式并给出余项的一致有效估计.     

Prewavelet approximations for a system of boundary integral equations for plates with free edges on polygons

We consider the plate equation in a polygonal domain with free edges. Its resolution by boundary integral equations is considered with double layer potentials whose variational formulation was given in Reference 25. We approximate its solution (u, (∂u/∂n)) by the Galerkin method with approximated spaces made of piecewise polynomials of order 2 and 1 for, respectively, u and (∂u/∂n). A prewavelet basis of these subspaces is built and equivalences between some Sobolev norms and discrete ones are established in the spirit of References 14, 16, 30 and 31. Further, a compression procedure is presented which reduces the number of nonzero entries of the stiffness matrix from O(N²) to O(N log N), where N is the size of this matrix. We finally show that the compressed stiffness matrices have a condition number uniformly bounded with respect to N and that the compressed Galerkin scheme converges with the same rate than the Galerkin one. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

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.     

Boundary integral equations in time-dependent bending of thermoelastic plates

The initial-boundary value problems with the Dirichlet and Neumann boundary conditions arising in the theory of bending of thermoelastic plates with transverse shear deformation are reduced to time-dependent boundary integral equations by means of layer potentials. The solvability of these equations is then investigated in Sobolev-type spaces.     

A conformal mapping algorithm for the Bernoulli free boundary value problem

We propose a new numerical method for the solution of the Bernoulli free boundary value problem for harmonic functions in a doubly connected domain D in where an unknown free boundary Γ₀ is determined by prescribed Cauchy data on Γ₀ in addition to a Dirichlet condition on the known boundary Γ₁. Our main idea is to involve the conformal mapping method as proposed and analyzed by Akduman, Haddar, and Kress for the solution of a related inverse boundary value problem. For this, we interpret the free boundary Γ₀ as the unknown boundary in the inverse problem to construct Γ₀ from the Dirichlet condition on Γ₀ and Cauchy data on the known boundary Γ₁. Our method for the Bernoulli problem iterates on the missing normal derivative on Γ₁ by alternating between the application of the conformal mapping method for the inverse problem and solving a mixed Dirichlet–Neumann boundary value problem in D. We present the mathematical foundations of our algorithm and prove a convergence result. Some numerical examples will serve as proof of concept of our approach. Copyright © 2015 John Wiley & Sons, Ltd.     

The coupling of boundary integral and finite element methods for the nonstationary exterior flow

In this article, we represent a new numerical method for solving the nonstationary Navier–Stokes equations in an unbounded domain. The technique consists of coupling the boundary integral and the finite element method. The variational formulation and the well-posedness of the coupling method are obtained. The convergence and optimal error estimates for the approximate solution are provided. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 549–565, 1998     

Dynamics for a nonlocal competition system with a free boundary

In this paper, we investigate a nonlocal reaction–diffusion competition model with a free boundary and discuss the long time behavior of species. The main objective is to understand the effect of the nonlocal term in the form of an integral convolution on the dynamics of competing species. Specially, for the weak competition case, when spreading occurs, we provide some sufficient conditions to prove that two competing species stabilize at a positive constant equilibrium state. Furthermore, for the case of successful spreading, we estimate the asymptotic spreading speed of the free boundary.     

Wavelet-based preconditioners for boundary integral equations

We study simple preconditioners for the conjugate gradient method when used to solve matrix systems arising from some hypersingular and weakly singular integral equations. The preconditioners, which are of the type of hierarchical basis preconditioners, are based on the decomposition of the piecewise-linear (respectively piecewise-constant) functions as the sum of prewavelets (respectively derivatives of prewavelets). We prove that with these preconditioners the preconditioned systems have condition numbers uniformly bounded with respect to the degrees of freedom. Numerical experiments support our analysis. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.

     

Parabolic finite volume element equations in nonconvex polygonal domains

We study spatially semidiscrete and fully discrete finite volume element approximations of the heat equation with homogeneous Dirichlet boundary conditions in a plane polygonal domain with one reentrant corner. We show that, as a result of the singularity in the solution near the reentrant corner, the convergence rate is reduced from optimal second order, similarly to what was shown for the finite element method in the earlier work 2 . Optimal order convergence may be restored by mesh refinement near the corners of the domain. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Shape-explicit constants for some boundary integral operators

Among the well-known constants in the theory of boundary integral equations are the coercivity constants of the single-layer potential and the hypersingular boundary integral operator, and the contraction constant of the double-layer potential. Whereas there have been rigorous studies how these constants depend on the size and aspect ratio of the underlying domain, only little is known on their dependency on the shape of the boundary. In this article, we consider the homogeneous Laplace equation and derive explicit estimates for the above-mentioned constants in three dimensions. Using an alternative trace norm, we make the dependency explicit in two geometric parameters, the so-called Jones parameter and the constant in Poincaré's inequality. The latter one can be tracked back to the constant in an isoperimetric inequality. There are many domains with quite irregular boundaries, where these parameters stay bounded. Our results provide a new tool in the analysis of numerical methods for boundary integral equations and in particular for boundary element based domain decomposition methods.     

Meshless Galerkin algorithms for boundary integral equations with moving least square approximations

In this paper, we first give error estimates for the moving least square (MLS) approximation in the Hk norm in two dimensions when nodes and weight functions satisfy certain conditions. This two-dimensional error results can be applied to the surface of a three-dimensional domain. Then combining boundary integral equations (BIEs) and the MLS approximation, a meshless Galerkin algorithm, the Galerkin boundary node method (GBNM), is presented. The optimal asymptotic error estimates of the GBNM for three-dimensional BIEs are derived. Finally, taking the Dirichlet problem of Laplace equation as an example, we set up a framework for error estimates of the GBNM for boundary value problems in three dimensions.     

Numerical solution of boundary integral equations by means of attenuation factors

We consider first-kind boundary integral equations with logarithmickernel such as those arising from solving Dirichlet problemsfor the Laplace equation by means of single-layer potentials.The first-kind equations are transformed into equivalent equationsof the second kind which contain the conjugation operator andwhich are then solved with a degenerate-kernel method basedon Fourier analysis and attenuation factors. The approximationswe consider, among them spline interpolants, are linear andtranslation invariant. In view of the particularly small kernel,the linear systems resulting from the discretization can besolved directly by fixed-point iteration.     

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.     

Numerical treatment of retarded boundary integral equations by sparse panel clustering

Stefan Sauter Institute for Mathematics, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland ^{We consider the wave equation in a boundary integral formulation.The discretization in time is done by using convolution quadraturetechniques and a Galerkin boundary element method for the spatialdiscretization. In a previous paper, we have introduced a sparseapproximation of the system matrix by cut-off, in order to reducethe storage costs. In this paper, we extend this approach byintroducing a panel clustering method to further reduce thesecosts.}