含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

Existence and uniqueness of solutions for acoustic scattering over infinite obstacles

The paper considers the solution of the boundary value problem (BVP) consisting of the Helmholtz equation in the region D with a rigid boundary condition on ∂D and its reformulation as a boundary integral equation (BIE), over an infinite cylindrical surface of arbitrary smooth cross-section. A boundary integral equation, which models three-dimensional acoustic scattering from an infinite rigid cylinder, illustrates the application of the above results to prove existence of solution of the integral equation and the corresponding boundary value problem.     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.     

A Symmetric Galerkin Boundary Element Method for Linear Poroelasticity

In linear poroelasticity so far only collocation boundary element methods have been available. However, in some applications, e.g., when coupling with finite elements is desired, a symmetric formulation is preferable. Choosing a Galerkin approach which involves the second boundary integral equation, such a formulation is possible. Here, a previously presented integration by part technique for the regularization of the first boundary integral equation is extended to the second boundary integral equation as well. While the weakly singular representation of the double layer operator has been presented before, the emphasis lies here on the so called hyper-singular boundary integral operator. Due to the regularization, this operator can be evaluated numerically and, hence, be used within a numerical scheme for the first time. Different numerical studies will be presented to show the behavior of the established symmetric Galerkin boundary element method, also comparing it with collocation boundary element methods. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of a numerical scheme of the type of the vortex loop method on a closed surface with an approximation to the surface shape

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises when solving the Neumann boundary value problem for the Laplace equation with the use of the representation of the solution in the form of a double layer potential. We study the case in which an exterior or interior boundary value problem is solved in a domain whose boundary is a smooth closed surface and the integral equation is written out on that surface. For the numerical solution of the integral equation, the surface is approximated by spatial polygons whose vertices lie on the surface. We construct a numerical scheme for solving the integral equation on the basis of such an approximation to the surface with the use of quadrature formulas of the type of the method of discrete singularities with regularization. We prove that the numerical solutions converge to the exact solution of the hypersingular integral equation uniformly on the grid.     

Integral Equation Method for the First and Second Problems of the Stokes System

Using the integral equation method we study solutions of boundary value problems for the Stokes system in Sobolev space H ¹(G) in a bounded Lipschitz domain G with connected boundary. A solution of the second problem with the boundary condition $\partial {\bf u}/\partial {\bf n} -p{\bf n}={\bf g}$ is studied both by the indirect and the direct boundary integral equation method. It is shown that we can obtain a solution of the corresponding integral equation using the successive approximation method. Nevertheless, the integral equation is not uniquely solvable. To overcome this problem we modify this integral equation. We obtain a uniquely solvable integral equation on the boundary of the domain. If the second problem for the Stokes system is solvable then the solution of the modified integral equation is a solution of the original integral equation. Moreover, the modified integral equation has a form f?+?S f?=?g, where S is a contractive operator. So, the modified integral equation can be solved by the successive approximation. Then we study the first problem for the Stokes system by the direct integral equation method. We obtain an integral equation with an unknown ${\bf g}=\partial {\bf u}/\partial {\bf n} -p{\bf n}$ . But this integral equation is not uniquely solvable. We construct another uniquely solvable integral equation such that the solution of the new eqution is a solution of the original integral equation provided the first problem has a solution. Moreover, the new integral equation has a form ${\bf g}+\tilde S{\bf g}={\bf f}$ , where $\tilde S$ is a contractive operator, and we can solve it by the successive approximation.     

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

$C^1$ error estimation on the boundary for an exterior Neumann problem in $\mathbb R^3$

Summary.   In this paper we establish a error estimation on the boundary for the solution of an exterior Neumann problem in . To solve this problem we consider an integral representation which depends from the solution of a boundary integral equation. We use a full piecewise linear discretisation which on one hand leads to a simple numerical algorithm but on the other hand the error analysis becomes more difficult due to the singularity of the integral kernel. We construct a particular approximation for the solution of the boundary integral equation, for the solution of the Neumann problem and its gradient on the boundary and estimate their error. Received May 11, 1998 / Revised version received July 7, 1999 / Published online August 24, 2000     

Helmholtz equation outside an open arc in a plane with a mixed boundary condition

A boundary value problem for the Helmholtz equation outside an open arc in a plane is studied with mixed boundary conditions. In doing so, the Dirichlet condition is specified on one side of the open arc and the boundary condition of the third kind is specified on the other side of the open arc. We consider non-propagative Helmholtz equation, real-valued solutions of which satisfy maximum principle. By using the potential theory the boundary value problem is reduced to a system of singular integral equations with additional conditions. By regularization and subsequent transformations, this system is reduced to a vector Fredholm equation of the second kind and index zero. It is proved that the obtained vector Fredholm equation is uniquely solvable. Therefore the integral representation for a solution of the original boundary value problem is obtained. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Behaviour of Elastic Fields and Boundary Integral Mellin Techniques Near Conical Points

For the computation of the local singular behaviour of an homogeneous anisotropic clastic field near the three-dimensional vertex subjected to displacement boundary conditions, one can use a boundary integral equation of the first kind whose unkown is the boundary stress. Mellin transformation yields a one - dimensional integral equation on the intersection curve 7 of the cone with the unit sphere. The Mellin transformed operator defines the singular exponents and Jordan chains, which provide via inverse Mellin transformation a local expansion of the solution near the vertex. Based on Kondratiev's technique which yields a holomorphic operator pencil of elliptic boundary value problems on the cross - sectional interior and exterior intersection of the unit sphere with the conical interior and exterior original cones, respectively, and using results by Maz'ya and Kozlov, it can be shown how the Jordan chains of the one-dimensional boundary integral equation are related to the corresponding Jordan chains of the operator pencil and their jumps across γ. This allows a new and detailed analysis of the asymptotic behaviour of the boundary integral equation solutions near the vertex of the cone. In particular, the integral equation method developed by Schmitz, Volk and Wendland for the special case of the elastic Dirichlet problem in isotropic homogeneous materials could be completed and generalized to the anisotropic case.     

无界区域抛物方程自然边界元方法

本文应用自然边界元方法求解无界区域抛物型初边值问题。首先将控制方程对时间进行离散化,得到关于时间步长离散化的椭圆型问题。通过Fourier展开,导出相应问题的自然积分方程和Poisson积分公式。研究了自然积分算子的性质,并讨论了自然积分方程的数值解法,最后给出数值例子。从而解决了抛物型问题的自然边界归化和自然边界元方法。     

分数阶微分方程边值问题解的存在性

应用Gteen函数将分数阶微分方程边值问题可转化为等价的积分方程．近来此方法被应用于讨论非线性分数阶微分方程边值问题解的存在性．讨论非线性分数阶微分方程边值问题,应用Green函数,将其转化为等价的积分方程,并设非线性项满足Caratheodory条件,利用非紧性测度的性质和M6nch’s不动点定理证明解的存在性．     

A new integral equation approach to the Neumann problem in acoustic scattering

We suggest a new approach of reduction of the Neumann problem in acoustic scattering to a uniquely solvable Fredholm integral equation of the second kind with weakly singular kernel. To derive this equation we placed an additional boundary with an appropriate boundary condition inside the scatterer. The solution of the problem is obtained in the form of a single layer potential on the whole boundary. The density in the potential satisfies a uniquely solvable Fredholm integral equation of the second kind and can be computed by standard codes. Copyright © 2001 John Wiley & Sons, Ltd.     

Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

Some integral relations of Hankel transform type and applications to elasticity theory

The dominant part of an integral equation arising in connection with boundary value problems for the circular disc is evaluated in terms of orthogonal polynomials. This relation leads to an efficient method for numerical solution of the complete integral equation even in the presence of a complicated bounded kernel. The static problem of a circular crack in an infinite elastic body under general loads is used to illustrate vector boundary conditions leading to two coupled integral equations, while the problem of a vibrating flexible circular plate in frictionless contact with an elastic half space is solved by use of the associated numerical method.     

Time-space boundary elements for transient heat conduction problems

This paper is concerned with a generalized time-space boundary element formulation for transient heat conduction problems in anisotropic media. A weighted residual form of the governing equation is used to obtain the boundary integral equation in terms of the fundamental solution. The resulting boundary integral equation is discretized by means of a wide variety of boundary elements from constant-elements to higher-order isoparametric elements located both in time and space.     

The Impedance Boundary Value Problem for the Helmholtz Equation in a Half-Plane

We prove unique existence of solution for the impedance (or third) boundary value problem for the Helmholtz equation in a half-plane with arbitrary L_∞ boundary data. This problem is of interest as a model of outdoor sound propagation over inhomogeneous flat terrain and as a model of rough surface scattering. To formulate the problem and prove uniqueness of solution we introduce a novel radiation condition, a generalization of that used in plane wave scattering by one-dimensional diffraction gratings. To prove existence of solution and a limiting absorption principle we first reformulate the problem as an equivalent second kind boundary integral equation to which we apply a form of Fredholm alternative, utilizing recent results on the solvability of integral equations on the real line in [5]. © 1997 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

The coupling of boundary integral and finite element methods for the bidimensional exterior steady stokes problem

In this paper, we represent a new numerical method for solving the steady-state Stokes equations in an unbounded plane domain. The technique consists in coupling the boundary integral and the finite element methods. An artificial smooth boundary is introduced separating an interior inhomogeneous region from an exterior one. The solution in the exterior domain is represented by an integral equation over the artificial boundary. This integral equation is incorporated into a velocitypressure formulation for the interior region, and a finite element method is used to approximate the resulting variational problem. This is studied by means of an abstract framework, well adapted to the model problem, in which convergence results and optimal error estimates are derived. Computer results will be discussed in a forthcoming paper.     

The first initial boundary value problem for a compound type equation in a three-dimensional multiply connected region

The method of boundary integral equations is used for solving the first initial boundary value problem for a compound type equation in a three-dimensional multiply connected region. The problem is reduced to a uniquely solvable integral equation. The solution of the problem is obtained in the form of dynamic potentials whose density satisfies this integral equation. Thus the existence theorem is proved. Moreover, the uniqueness of the solution is also studied. All the results are valid for both interior and exterior regions provided that the corresponding conditions at infinity are taken into account. Translated fromMatematicheskie Zametki, Vol. 68, No. 2, pp. 249–265, August, 2000.     

New indirect scalar boundary integral equation formulas for the biharmonic equation

We derive scalar boundary integral equation formulas for both interior and exterior biharmonic equations with the Dirichlet boundary data. They are based on indirect boundary integral equation formulas, so-called the Chakrabarty and Almansi formulas. The scalar formulas are derived through an unconventional variational approach. The unique solvability results of the formulas are also obtained.