Compressed absorbing boundary conditions via matrix probing

Absorbing layers are sometimes required to be impractically thick in order to offer an accurate approximation of an absorbing boundary condition for the Helmholtz equation in a heterogeneous medium. It is always possible to reduce an absorbing layer to an operator at the boundary by layer-stripping elimination of the exterior unknowns, but the linear algebra involved is costly. We propose to bypass the elimination procedure, and directly fit the surface-to-surface operator in compressed form from a few exterior Helmholtz solves with random Dirichlet data. The result is a concise description of the absorbing boundary condition, with a complexity that grows slowly (often, logarithmically) in the frequency parameter.Comment: 29 pages with 25 figure

Viscous theory of surface noise interaction phenomena

A viscous linear surface noise interaction problem is formulated that includes noise production by an oscillating surface, turbulent or vortical interaction with a surface, and scattering of sound by a surface. The importance of viscosity in establishing uniqueness of solution and partitioning of energy into acoustic and vortical modes is discussed. The results of inviscid two dimensional airfoil theory are used to examine the interactive noise problem in the limit of high reduced frequency and small Helmholtz number. It is shown that in the case of vortex interaction with a surface, the noise produced with the full Kutta condition is 3 dB less than the no Kutta condition result. The results of a study of an airfoil oscillating in a medium at rest are discussed. It is concluded that viscosity can be a controlling factor in analyses and experiments of surface noise interaction phenomena and that the effect of edge bluntness as well as viscosity must be included in the problem formulation to correctly calculate the interactive noise

Residual based adaptivity and PWDG methods for the Helmholtz equation

We present a study of two residual a posteriori error indicators for the Plane Wave Discontinuous Galerkin (PWDG) method for the Helmholtz equation. In particular we study the h-version of PWDG in which the number of plane wave directions per element is kept fixed. First we use a slight modification of the appropriate a priori analysis to determine a residual indicator. Numerical tests show that this is reliable but pessimistic in that the ratio between the true error and the indicator increases as the mesh is refined. We therefore introduce a new analysis based on the observation that sufficiently many plane waves can approximate piecewise linear functions as the mesh is refined. Numerical results demonstrate an improvement in the efficiency of the indicators

The jump problem for the laplace equation

AbstractThe boundary value problem for the Laplace equation outside several cuts in a plane is studied. The jump of the solution of the Laplace equation and the jump of its normal derivative are specified of the cuts. The problem is studied under different conditions at infinity, which lead to different uniqueness and existence theorems. The solution of this problem is constructed in the explicit form by means of single-layer and angular potentials. The singularities at the ends of the cuts are investigated

Accurate computation of Galerkin double surface integrals in the 3-D boundary element method

Many boundary element integral equation kernels are based on the Green's functions of the Laplace and Helmholtz equations in three dimensions. These include, for example, the Laplace, Helmholtz, elasticity, Stokes, and Maxwell's equations. Integral equation formulations lead to more compact, but dense linear systems. These dense systems are often solved iteratively via Krylov subspace methods, which may be accelerated via the fast multipole method. There are advantages to Galerkin formulations for such integral equations, as they treat problems associated with kernel singularity, and lead to symmetric and better conditioned matrices. However, the Galerkin method requires each entry in the system matrix to be created via the computation of a double surface integral over one or more pairs of triangles. There are a number of semi-analytical methods to treat these integrals, which all have some issues, and are discussed in this paper. We present novel methods to compute all the integrals that arise in Galerkin formulations involving kernels based on the Laplace and Helmholtz Green's functions to any specified accuracy. Integrals involving completely geometrically separated triangles are non-singular and are computed using a technique based on spherical harmonics and multipole expansions and translations, which results in the integration of polynomial functions over the triangles. Integrals involving cases where the triangles have common vertices, edges, or are coincident are treated via scaling and symmetry arguments, combined with automatic recursive geometric decomposition of the integrals. Example results are presented, and the developed software is available as open source