Nonconforming Galerkin methods for the Helmholtz equation

Nonconforming Galerkin methods for a Helmholtz‐like problem arising in seismology are discussed both for standard simplicial linear elements and for several new rectangular elements related to bilinear or trilinear elements. Optimal order error estimates in a broken energy norm are derived for all elements and in L² for some of the elements when proper quadrature rules are applied to the absorbing boundary condition. Domain decomposition iterative procedures are introduced for the nonconforming methods, and their convergence at a predictable rate is established. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 475–494, 2001     

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R~2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.     

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

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

Legendre Gauss Spectral Collocation for the Helmholtz Equation on a Rectangle

A spectral collocation method with collocation at the Legendre Gauss points is discussed for solving the Helmholtz equation –u+(x,y)u=f(x,y) on a rectangle with the solution u subject to inhomogeneous Robin boundary conditions. The convergence analysis of the method is given in the case of u satisfying Dirichlet boundary conditions. A matrix decomposition algorithm is developed for the solution of the collocation problem in the case the coefficient (x,y) is a constant. This algorithm is then used in conjunction with the preconditioned conjugate gradient method for the solution of the spectral collocation problem with the variable coefficient (x,y).     

环形域上半线性椭圆方程正解的存在性

本文利用连续性方法,得到了一类半线性椭圆方程第一边值问题在环形域上径向对称正解的存在性．     

A Pseudo-Kinetic Approach for Helmholtz Equation

A lattice Boltzmann type pseudo-kinetic model for a non-homogeneous Helmholtz equation is derived in this paper. Numerical results for some model problems show the robustness and efficiency of this lattice Boltzmann type pseudo-kinetic scheme. The computation at each site is determined only by local parameters, and can be easily adapted to solve multiple scattering problems with many scatterers or wave propagation in non-homogeneous medium without increasing the computational cost.     

A Radiation Condition for the 2-D Helmholtz Equation in Stratified Media

We study the 2-D Helmholtz equation in perturbed stratified media, allowing the existence of guided waves. Our assumptions on the perturbing and source terms are not too restrictive. We prove two results. Firstly, we introduce a Sommerfeld–Rellich radiation condition and prove the uniqueness of the solution for the studied equation. Then, by careful asymptotic estimates, we prove the existence of a bounded solution satisfying our radiation condition.     

Fourier Moment Method with Regularization for the Cauchy Problem of Helmholtz Equation

In this paper, we consider the reconstruction of the wave field in a bounded domain. By choosing a special family of functions, the Cauchy problem can be transformed into a Fourier moment problem. This problem is ill-posed. We propose a regularization method for obtaining an approximate solution to the wave field on the unspecified boundary. We also give the convergence analysis and error estimate of the numerical algorithm. Finally, we present some numerical examples to show the effectiveness of this method.     

Limiting Absorption Principle for the Dissipative Helmholtz Equation

Adapting Mourre's commutator method to the dissipative setting, we prove a limiting absorption principle for a class of abstract dissipative operators. A consequence is the uniform resolvent estimate for the high-frequency Helmholtz equation when trapped classical trajectories meet the region where the absorption coefficient is non-zero. We also give the resolvent estimate in Besov spaces.     

A Fourth-Order Modified Method for the Cauchy Problem of the Modified Helmholtz Equation

This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0 < x≤ 1, y ∈ R. The Cauchy data at x = 0 is given and the solution is then sought for the interval 0 < x ≤1. This problem is highly ill-posed and the solution (if it exists) does not depend continuously on the given data. In this paper, we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution. Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.     

A nonconforming domain decomposition approximation for the Helmholtz screen problem with hypersingular operator

We present and analyze a nonconforming domain decomposition approximation for a hypersingular operator governed by the Helmholtz equation in three dimensions. This operator appears when considering the corresponding Neumann problem in unbounded domains exterior to open surfaces. We consider small wave numbers and low‐order approximations with Nitsche coupling across interfaces. Under appropriate assumptions on mapping properties of the weakly singular and hypersingular operators with Helmholtz kernel, we prove that this method converges almost quasioptimally, that is, with optimal orders reduced by an arbitrarily small positive number. Numerical experiments confirm our error estimate. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 125–141, 2017     

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.     

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

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

求解修正的Helmholtz方程关于非光滑区域问题

研究了修正的(纯虚波数)Helmholtz方程在阻尼边界条件下,求解含单个角点的闭区域问题.通过采用单双层混合位势来表示其解,进而对其角型区域进行求解.最后,通过数值例子来说明此方法的可行性与可靠性.     

Comparing numerical methods for Helmholtz equation model problem

In this article, we implement a relatively new numerical technique, Adomian’s decomposition method for solving the linear Helmholtz partial differential equations. The method in applied mathematics can be an effective procedure to obtain for the analytic and approximate solutions. A new approach to a linear or nonlinear problems is particularly valuable as a tool for Scientists and Applied Mathematicians, because it provides immediate and visible symbolic terms of analytic solution as well as its numerical approximate solution to both linear and nonlinear problems without linearization [Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic Publishers, Boston, 1994; J. Math. Anal. Appl. 35 (1988) 501]. It does also not require discretization and consequently massive computation. In this scheme the solution is performed in the form of a convergent power series with easily computable components. This paper will present a numerical comparison with the Adomian decomposition and a conventional finite-difference method. The numerical results demonstrate that the new method is quite accurate and readily implemented.     

六次对角型方程的最小解

本文综合运用三种迭代方法证明了方程当ｓ取最小可能的值３７时有非零解，并给出了解的一个上界估计。这在变量个数方面改进了Ｊ．Ｐｉｔｍａｎ的结果。     

Stable Determination of Polyhedral Interfaces from Boundary Data for the Helmholtz Equation

We study an inverse boundary value problem for the Helmholtz equation using the Dirichlet-to-Neumann map. We consider piecewise constant wave speeds on an unknown tetrahedral partition and prove a Lipschitz stability estimate in terms of the Hausdorff distance between partitions.     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Noniterative Reconstruction Method for an Inverse Potential Problem Modeled by a Modified Helmholtz Equation

This paper deals with an inverse potential problem posed in two dimensional space whose forward problem is governed by a modified Helmholtz equation. The inverse problem consists in the reconstruction of a set of anomalies embedded into a geometrical domain from partial measurements of the associated potential. Since the inverse problem, we are dealing with, is written in the form of an ill-posed boundary value problem, the idea is to rewrite it as a topology optimization problem. In particular, a shape functional is defined to measure the misfit of the solution obtained from the model and the data taken from the partial measurements. This shape functional is minimized with respect to a set of ball-shaped anomalies using the concept of topological derivatives. It means that the shape functional is expanded asymptotically and then truncated up to the desired order term. The resulting expression is trivially minimized with respect to the parameters under consideration which leads to a noniterative second-order reconstruction algorithm. As a result, the reconstruction process becomes very robust with respect to noisy data and independent of any initial guess. Finally, some numerical experiments are presented to show the effectiveness of the proposed reconstruction algorithm.     

Application of Adapted-Bubbles to the Helmholtz Equation with Large Wavenumbers in 2D

An adapted-bubbles approach which is a modification of the residual-free bubbles (RFB) method, is proposed for the Helmholtz problem in 2D. A new two-level finite element method is introduced for the approximations of the bubble functions. Unlike the other equations such as the advection-diffusion equation, RFB method when applied to the Helmholtz equation, does not depend on another stabilized method to obtain approximations to the solutions of the sub-problems. Adapted-bubbles (AB) are obtained by a simple modification of the sub-problems. This modification increases the accuracy of the numerical solution impressively. We provide numerical experiments with the AB method up to $ch = 5$ where $c$ is the wavenumber and $h$ is the mesh size. Numerical tests show that the AB method is better by far than higher order methods available in the literature.