Extrapolation of Nyström Solutions of Boundary Integral Equations on Non-Smooth Domains

The interior Dirichlet problem for Laplace's equation on a plane polygonal region $\Omega$ with boundary $\Gamma$ may be reformulated as a second kind integral equation on $\Gamma$. This equation may be solved by the Nyström method using the composite trapezoidal rule. It is known that if the mesh has $O(n)$ points and is graded appropriately, then $O(1/n^2)$ convergence is obtained for the solution of the integral equation and the associated solution to the Dirichlet problem at any $x\in \Omega$. We present a simple extrapolation scheme which increases these rates of convergence to $O(1/n^4)$ .     

An adapted Galerkin method for the resolution of Dirichlet and Neumann problems in a polygonal domain

The Neumann problem for Laplace's equation in a polygonal domain is associated with the exterior Dirichlet problem obtained by requiring the continuity of the potential through the boundary. Then the solution is the simple layer potential of the charge q on the boundary. q is the solution of a Fredholm integral equation of the second kind that we solve by the Galerkin method. The charge q has a singular part due to the corners, so the optimal order of convergence is not reached with a uniform mesh. We restore this optimal order by grading the mesh adequately near the corners. The interior Dirichlet problem is solved analogously, by expressing the solution as a double layer potential.     

Integral equation methods for scattering by infinite rough surfaces

In this paper, we consider the Dirichlet and impedance boundary value problems for the Helmholtz equation in a non‐locally perturbed half‐plane. These boundary value problems arise in a study of time‐harmonic acoustic scattering of an incident field by a sound‐soft, infinite rough surface where the total field vanishes (the Dirichlet problem) or by an infinite, impedance rough surface where the total field satisfies a homogeneous impedance condition (the impedance problem). We propose a new boundary integral equation formulation for the Dirichlet problem, utilizing a combined double‐ and single‐layer potential and a Dirichlet half‐plane Green's function. For the impedance problem we propose two boundary integral equation formulations, both using a half‐plane impedance Green's function, the first derived from Green's representation theorem, and the second arising from seeking the solution as a single‐layer potential. We show that all the integral equations proposed are uniquely solvable in the space of bounded and continuous functions for all wavenumbers. As an important corollary we prove that, for a variety of incident fields including an incident plane wave, the impedance boundary value problem for the scattered field has a unique solution under certain constraints on the boundary impedance. Copyright © 2003 John Wiley & Sons, Ltd.     

A wavelets method for plane elasticity problem

A Neumann boundary value problem of plane elasticity problem in the exterior circular domain is reduced into an equivalent natural boundary integral equation and a Poisson integral formula with the DtN method. Using the trigonometric wavelets and Galerkin method, we obtain a fast numerical method for the natural boundary integral equation which has an unique solution in the quotient space. We decompose the stiffness matrix in our numerical method into four circulant and symmetrical or antisymmetrical submatrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform (FFT) and the inverse fast Fourier transform (IFFT) instead of the inverse matrix. Examples are given for demonstrating our method has good accuracy of our method even though the exact solution is almost singular.     

Boundary Integral Equations for Mixed Boundary Value Problems in R3

Both exterior and interior mixed Dirichlet-Neumann problems in R³ for the scalar Helmholtz equation are solved via boundary integral equations. The integral equations are equivalent to the original problem in the sense that the traces of the weak seolution satisfy the integral equations, and, conversely, the solution of the integral equations inserted into Green's formula yields the solution of the mixed boundary value problem. The calculus of pseudodifferential operators is used to prove existence and regularity of the solution of the integral equations. The regularity results — obtained via Wiener-Hopf technique — show the explicit “edge” behavior of the solution near the submanifold which separates the Dirichlet boundary from the Neumann boundary.     

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.     

平面中Poisson方程的Dirichlet问题

基于向量旋转内积不变的特点,通过对Green积分公式的推广,得到与空间Green三个公式相似的平面Green公式.从而,得到平面中Poisson方程Robin问题的解和平面中Poisson方程Dirichlet问题的解.     

嵌在弹性半空间的刚性变直径圆轴的扭转*

本文用线载荷积分方程法(LLIEM)研究嵌在弹性半空间的刚性变直径圆轴的轴对称扭转.利用将“半空间的点环力偶”(PRCHS)这一虚的基本载荷沿对称轴的圆轴区间中分布就能使本问题归结为一个一维的非奇异的Ferdholm第一种积分方程,且很易用数值求解.文中给出刚性圆锥,圆柱,圆锥柱嵌在弹性半空间的扭转的数值计算的例并和他人的已知结果相比较.文中也给出了刚性半球嵌在弹性半空间的扭转的精确解.     

Alternative derivations for the Poisson integral formula

Poisson integral formula is revisited. The kernel in the Poisson integral formula can be derived in a series form through the direct BEM free of the concept of image point by using the null-field integral equation in conjunction with the degenerate kernels. The degenerate kernels for the closed-form Green's function and the series form of Poisson integral formula are also derived. Two and three-dimensional cases are considered. Also, interior and exterior problems are both solved. Even though the image concept is required, the location of image point can be determined straightforward through the degenerate kernels instead of the method of reciprocal radii.     

The exterior dirichlet problem for the Laplace equation

Using a standard application of Green's theorem, the exterior Dirichlet problem for the Laplace equation in three dimensions is reduced to a pair of integral equations. One integral equation is of the second kind and the other is of the first kind. It is known that the integral equation of the second kind is not uniquely solvable, however, it has been demonstrated that the pair of integral equations has a unique solution. The present approach is based on the observation that the known function appearing in the integral equation of the second kind lies in a certain Banach space E which is a proper subspace of the Banach space of continuous complex-valued functions equipped with the maximum norm. Furthermore, it can be shown that the related integral operator when restricted to E has spectral radius less than unity. Consequently, a particular solution to the integral equation of the second kind can be obtained by the method of successive approximations and the unique solution to the problem is then obtained by using the integral equation of the first kind. Comparisons are made between the present algorithm and other known constructive methods. Finally, an example is supplied to illustrate the method of this paper.     

Sharp Growth Estimates for Modified Poisson Integrals in a Half Space

For continuous boundary data, including data of polynomial growth, modified Poisson integrals are used to write solutions to the half space Dirichlet and Neumann problems in Rⁿ. Pointwise growth estimates for these integrals are given and the estimates are proved sharp in a strong sense. For decaying data, a new type of modified Poisson integral is introduced and used to develop asymptotic expansions for solutions of these half space problems.     

A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

Half Dirichlet problem for matrix functions on the unit ball in Hermitian Clifford analysis

The simultaneous null solutions of the two complex Hermitian Dirac operators are focused on in Hermitian Clifford analysis, where the Hermitian Cauchy integral was constructed and will play an important role in the framework of circulant (2×2) matrix functions. Under this setting we will present the half Dirichlet problem for circulant (2×2) matrix functions on the unit ball of even dimensional Euclidean space. We will give the unique solution to it merely by using the Hermitian Cauchy transformation, get the solution to the Dirichlet problem on the unit ball for circulant (2×2) matrix functions and the solution to the classical Dirichlet problem as the special case, derive a decomposition of the Poisson kernel for matrix Laplace operator, and further obtain the decomposition theorems of solution space to the Dirichlet problem for circulant (2×2) matrix functions.     

半空间中一类次调和函数的增长性质

在Rn的半空间{x∈Rn,xn>0}中,得到了具有Dirichlet数据的Poisson积分在自然的积分收敛条件下满足增长性质u(x)=o(|x|),这里|x|→∞,这一性质对于半空间中满足一定条件的次调和函数仍然成立.该结果把复平面C中解析函数的增长性质推广到了n-维Euclidean半空间,并且推广了n-维Euclidean半空间中某些经典的结果.     

双曲空间中Laplace-Beltrami方程的一个带位移的边值问题

讨论了双曲空间中Laplace-Beltrami方程的一个带位移的边值问题.首先将双曲空间中的Laplace-Beltrami方程的解转化为Clifford分析中的超正则函数.然后给出了超正则函数的Plemelj公式并讨论了相关奇异积分算子的性质,最后利用积分方程的方法和压缩不动点原理证明了Laplace-Beltrami方程的一个带位移的边值问题的解的存在性和唯一性.     

NATURAL BOUNDARY ELEMENT METHOD FOR THREE DIMENSIONAL EXTERIOR HARMONIC PROBLEM WITH AN INNER PROLATE SPHEROID BOUNDARY

In this paper, we study natural boundary reduction for Laplace equation with Dirichletor Neumann boundary condition in a three-dimensional unbounded domain, which is theoutside domain of a prolate spheroid. We express the Poisson integral formula and naturalintegral operator in a series form explicitly. Thus the original problem is reduced to aboundary integral equation on a prolate spheroid. The variational formula for the reducedproblem and its well-posedness are discussed. Boundary element approximation for thevariational problem and its error estimates, which have relation to the mesh size andthe terms after the series is truncated, are also presented. Two numerical examples arepresented to demonstrate the effectiveness and error estimates of this method.     

Maslov-Poisson measure and Feynman formulas for the solution of the Dirac equation

As the main step, the method used by V. P. Maslov for representing a solution of the initial-value problem for the classical Schrödinger equation and admitting an application to the Dirac equation includes the construction of a cylindrical countably additive measure (that is an analog of the Poisson distribution) on a certain space of functions (= trajectories in the impulse space) whose Fourier transform coincides with the factor in the formula for representation of the solution of the Schrödinger equation by the integral in the so-called cylindrical Feynman (pseudo) measure (in the trajectory space of the configurational space for the classical system). On the other hand, in the Maslov formula for the solution of the Schrödinger equation, the exponential factor is (with accuracy up to a shift) the Fourier transform of the Feynman pseudomeasure. In the case of the Dirac equation, historically, for the first time, there arose formulas for the impulse representation that use countably additive functional distributions of the Poisson-Maslov measure type but with noncommuting (matrix) values. The paper finds generalized measures whose Fourier transforms coincide with an analog of the exponential factor under the integral sign in the Maslov-type formula for the Dirac equation and integrals with respect to which yield solutions of the Cauchy problem for this equation in the configurational space.     

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

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

Integral representation and asymptotic behavior of harmonic functions in half space

Using Carleman's formula of a harmonic function in the half space and Nevanlinna's representation of a harmonic function in the half sphere, we prove that a harmonic function, whose positive part satisfies a slowly growing condition, can be represented by a certain integral. This improves some classical Poisson integrals for harmonic functions.     

OPTIMIZATION APPROACH FOR THE MONGE-AMPèRE EQUATION

In this paper, we introduce and study a method for the numerical solution of the elliptic Monge-Ampère equation with Dirichlet boundary conditions. We formulate the Monge-Ampère equation as an optimization problem. The latter involves a Poisson Problem which is solved by the finite element Galerkin method and the minimum is computed by the conjugate gradient algorithm. We also present some numerical experiments.