A Fast Algorithm for Two-Dimensional Elliptic Problems

In this paper, we extend the work of Daripa et al. [14–16,7] to a larger class of elliptic problems in a variety of domains. In particular, analysis-based fast algorithms to solve inhomogeneous elliptic equations of three different types in three different two-dimensional domains are derived. Dirichlet, Neumann and mixed boundary value problems are treated in all these cases. Three different domains considered are: (i) interior of a circle, (ii) exterior of a circle, and (iii) circular annulus. Three different types of elliptic problems considered are: (i) Poisson equation, (ii) Helmholtz equation (oscillatory case), and (iii) Helmholtz equation (monotone case). These algorithms are derived from an exact formula for the solution of a large class of elliptic equations (where the coefficients of the equation do not depend on the polar angle when written in polar coordinates) based on Fourier series expansion and a one-dimensional ordinary differential equation. The performance of these algorithms is illustrated for several of these problems. Numerical results are presented.     

竖直薄板浸没在表面覆盖冰层海水中时的水波散射

研究在一片均匀薄冰所覆盖的深水中,浸没其间的竖直平板引起的水波散射,冰层看作弹性薄板.通过对障碍物前方的势函数微分,问题被归结为一个超奇异的积分方程,应用适当的Green积分定理,应用一个包含Chebyshev多项式的有限级数配置法,求解该积分方程.得到反射系数和透射系数的数值结果,并在不同的波数和覆盖冰层参数下,用图形表示出来.     

A Second-Order Method for the Electromagnetic Scattering from a Large Cavity

In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.     

二维位势问题的快速多极展开式

快速多极边界元法已经成功地应用于大规模二维三维弹性静力学问题中,有效地减少了计算时间和存储需求.将基于Taylor展式地快速多极边界元法应用到二维位势问题中,提出了二维位势问题地快速多极边界元格式,建立了二维位势问题的快速多极展开式.     

一类奇性积分方程的多尺度配置快速算法

This paper develops fast multiscale collocation methods for a class of Fredholm integral equations of the second kind with singular kernels. A truncation strategy for the coefficient matrix of the corresponding discrete system is proposed, which forms a basis for fast algorithms. The convergence, stability and computational complexity of these algorithms are analyzed.     

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.     

Diagonal forms of the translation operators in the fast multipole algorithm for scattering problems

The integral equations of acoustic and electromagnetic scattering generate large dense systems of linear equations. These systems are efficiently solved with iterative methods where the matrix-vector multiplication is computed using a special fast method, such as the fast Fourier transform or the fast multipole method (FMM). In this paper, the so called diagonal forms of the translation operators for the fast multipole method are derived starting from integral representations of certain special functions. Error analysis of the FMM is given, considering both the truncation error of potential expansions and the errors from the use of numerical integration in the diagonal translation theorem. The implications of the error bounds on the FMM algorithm are discussed.This work has been financially supported by the Jenny and Antti Wihuri Foundation and by the Cultural Foundation of Finland.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.     

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 Comparison of Splittings and Integral Equation Solvers for a Nonseparable Elliptic Equation

Iterative numerical schemes for solving the electrostatic partial differential equation with variable conductivity, using fast and high-order accurate direct methods for preconditioning, are compared. Two integral method schemes of this type, based on previously suggested splittings of the equation, are proposed, analyzed, and implemented. The schemes are tested for large problems on a square. Particular emphasis is paid to convergence as a function of geometric complexity in the conductivity. Differences in performance of the schemes are predicted and observed in a striking manner. AMS subject classification (2000) 31A10, 35C15, 65R20.Received May 2004. Accepted September 2004. Communicated by Anders Szepessy.Johan Helsing: This work was supported by the Swedish Science Research Council under contract 621-2001-2799.     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.     

Simple error estimators for the Galerkin BEM for some hypersingular integral equation in 2D

A posteriori error estimation is an important tool for reliable and efficient Galerkin boundary element computations. For hypersingular integral equations in 2D with a positive-order Sobolev space, we analyse the mathematical relation between the (h???h/2)-error estimator from [S. Ferraz-Leite and D. Praetorius, Simple a posteriori error estimators for the h-version of the boundary element method, Computing 83 (2008), pp. 135–162], the two-level error estimator from [M. Maischak, P. Mund, and E. Stephan, Adaptive multilevel BEM for acoustic scattering, 585 Comput. Methods Appl. Mech. Eng. 150 (1997), pp. 351–367], and the averaging error estimator from [C. Carstensen and D. Praetorius, Averaging techniques for the a posteriori bem error control for a hypersingular integral equation in two dimensions, SIAM J. Sci. Comput. 29 (2007), pp. 782–810]. All of these a posteriori error estimators are simple in the following sense: first, the numerical analysis can be done within the same mathematical framework, namely localization techniques for the energy norm. Second, there is almost no implementational overhead for the realization.     

局部扰动半平面上的散射问题

该文讨论半平面上有局部扰动情况下的散射问题．通过位势理论,应用边界积分方程的方法研究了该问题解的存在与唯一性．主要方法是运用对称反射,使该无界区域上的散射问题变成一个有界区域上的散射问题,只是这一有界区域的边界不光滑．通过仔细分析相应的边界积分算子,作者得到了其解的存在与唯一性．     

含圆形孔口和水平直裂纹的平面问题的应力分析

利用复变函数方法和积分方程理论研究了既含有圆形孔口又含有水平裂纹的无限大平面的平面弹性问题,将复杂的解析函数的边值问题化成了求解只在裂纹上的奇异积分方程的问题.此外,还给出了裂纹尖端附近的应力场和应力强度因子的公式.     

Approximation of potential integral by radial bases for solutions of Helmholtz equation

For a Helmholtz equation Δu(x) + κ ² u(x) = f(x) in a region of R ^s , s ≥ 2, where Δ is the Laplace operator and κ = a + ib is a complex number with b ≥ 0, a particular solution is given by a potential integral. In this paper the potential integral is approximated by using radial bases with the order of approximation derived.      

关于Helmholtz外问题的边界积分方程解的唯一性问题

本文用能量分析的观点探讨了用边界积分方程描述Helmholtz外问题时,解的唯一性不能保持的原因.文中证明了,当利用积分方程来描述问题时,实际上将无穷远处的Sommerfeld条件改成了既适合于外向波(辐射波),又适合于内向波(吸收波),即整个系统的能量保持守恒.并根据此观点解释了保持唯一性的算法.     

A Method for Solving the Inverse Scattering Problem for Shape and Impedance of Crack

The inverse problem considered in this paper is to determine the shape and the impedance of crack from a knowledge of the time-harmonic incident field and the corresponding far field pattern of the scattered waves in two-dimension. The combined single- and double-layer potential is used to approach the scattered waves. As an important feature, this method does not require the solution of $u$ and $\partial u / \partial \nu$ at each iteration. An approximate method is presented and the convergence of this method is proven. Numerical examples are given to show that this method is both accurate and simple to use.     

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.     

A hybrid method for inverse cavity scattering problem for shape

In this paper, we give a hybrid method to numerically solve the inverse open cavity scattering problem for cavity shape, given the scattered solution on the opening of the cavity. This method is a hybrid between an iterative method and an integral equations method for solving the Cauchy problem. The idea of this hybrid method is simple, the operation is easy, and the computation cost is small. Numerical experiments show the feasibility of this method, even for cases with noise.     

求解变系数非齐次亥姆霍茨方程的边界单元法

本文利用拉普拉斯方程的基本解作为权函数，给出求解交系数非齐次亥姆霍茨方程的迭代格式，进而得到求解这类方程的边界元迭代法．文中给出了算例．最后，把本文给出的边界元迭代法与作者早些时候提出的边界元耦合法进行了比较．