三维Helmholtz方程外问题的自然积分方程及其数值解

用文[2,3]提出的自然边界归化方法来处理三维Helmholtz方程的外边值问题。在简要介绍如何用球谐展开的方法得到Helmholtz问题在外球域上的自然积分方程后,给出求解该自然积分方程的一种数值方法及相应的数值算例。     

Approximation of the parabolic equation and the wavefield of a point source in a layered random medium

This paper studies the wavefield of a source in a multidimensional randomly layered medium. They obtained asymptotical expressions of the wave statistical characteristics for different boundary conditions both in the framework of the parabolic equation approximation and the exact formulation of the boundary problem for the Helmholtz equation. It is shown that the presence of a small but finite absorption γ is most important for the statistics. The diffraction effects turn out to be like those of absorption, but γ cannot tend to zero in this problem. In an appendix they give the factorization formulae of the wave equation solution in a layered medium.     

Numerical solution of inverse scattering for near-field optics

A novel regularized recursive linearization method is developed for a two-dimensional inverse medium scattering problem that arises in near-field optics, which reconstructs the scatterer of an inhomogeneous medium located on a substrate from data accessible through photon scanning tunneling microscopy experiments. Based on multiple frequency scattering data, the method starts from the Born approximation corresponding to weak scattering at a low frequency, and each update is obtained by continuation on the wavenumber from solutions of one forward problem and one adjoint problem of the Helmholtz equation.     

Approximation of the parabolic equation and the wavefield of a point source in a layered random medium

Abstract

This paper studies the wavefield of a source in a multidimensional randomly layered medium. They obtained asymptotical expressions of the wave statistical characteristics for different boundary conditions both in the framework of the parabolic equation approximation and the exact formulation of the boundary problem for the Helmholtz equation. It is shown that the presence of a small but finite absorption γ is most important for the statistics. The diffraction effects turn out to be like those of absorption, but γ cannot tend to zero in this problem. In an appendix they give the factorization formulae of the wave equation solution in a layered medium.     

Simulation of laser beam propagation with a paraxial model in a tilted frame

We study the Schrödinger equation which comes from the paraxial approximation of the Helmholtz equation in the case where the direction of propagation is tilted with respect to the boundary of the domain. In a first part, a mathematical analysis is made which leads to an analytical formula of the solution in the simple case where the refraction index and the absorption coefficients are constant. Afterwards, we propose a numerical method for solving the initial problem which uses the previous analytical expression. Numerical results are presented. We also sketch an extension to a time dependent model which is relevant for laser–plasma interaction.     

Scattering from a rough layer of a random medium

This paper deals with scattering from a random-medium layer with rough boundaries. The fluctuations of the surface heights and medium permittivity are assumed to be small and smooth. All random quantities are assumed to be stationary and independent of each other. After the introduction of approximate boundary conditions, the system of partial differential equations is transformed into an integral equation where the fluctuations of the problem are represented as a zero-mean random operator. Employing smoothing, integral equations for the coherent fields are obtained. Use of the Helmholtz operator leads to solution for the coherent propagation constant while the boundary operators lead to coherent Fresnel coefficients. The characteristics of the results are illustrated by considering several examples.     

The factorization method in a rough-surface scattering problem: I

In the present paper the wave scattering problem on rough surface is considered for the Helmholtz equation with the Dirichlet boundary condition. An approximate solution is derived with using a factorization approach to the original Helmholtz equation. As a result, the system of two equations of parabolic type appears. The first system equation has an exact analytical solution whereas for the second one, an approximate solution, is considered in terms of perturbation series. It is shown that the obtained approximate solution is the modified classical small perturbation series with respect to small Rayleigh parameter. In Appendix A it is demonstrated that, when the derived perturbation series is converged, it is possible to summarize it and to represent the exact solution of original boundary problem in an analytical symbolical form.     

The factorization method in a rough-surface scattering problem: I

In the present paper the wave scattering problem on rough surface is considered for the Helmholtz equation with the Dirichlet boundary condition. An approximate solution is derived with using a factorization approach to the original Helmholtz equation. As a result, the system of two equations of parabolic type appears. The first system equation has an exact analytical solution whereas for the second one, an approximate solution, is considered in terms of perturbation series. It is shown that the obtained approximate solution is the modified classical small perturbation series with respect to small Rayleigh parameter. In Appendix A it is demonstrated that, when the derived perturbation series is converged, it is possible to summarize it and to represent the exact solution of original boundary problem in an analytical symbolical form.     

Computation of discrete transparent boundary conditions for the 2D Helmholtz equation

We present a family of non-local transparent boundary conditions for the 2D Helmholtz equation. The whole domain, on which the Helmholtz equation is defined, is decomposed into an interior and an exterior domain. The corresponding interior Helmholtz problem is formulated as a variational problem in a standard manner, representing a boundary value problem, whereas the exterior problem is posed as an initial value problem in the radial variable. This problem is then solved approximately by means of the Laplace transformation. The derived boundary conditions are asymptotically correct, model inhomogeneous exterior domains and are simple to implement.     

Multiple monopole scattering of sound waves from spherical particles in liquid and elastic media

A solution to the problem of the mean sound field in liquid and elastic media with spherical particles causing monopole scattering of sound is proposed. The integral equation obtained for the field allows passage to the Helmholtz equation with an effective wave number. The characteristic features of the solution are the absence of radiation loss in the mean field wave and the absence of limitations on the particle concentration. The integral equation is used as the basis for solving the problem of the incidence of a plain sound wave at an arbitrary angle on a plane layer of a medium with particles.     

Equation for wave processes in inhomogeneous moving media and functional solution of the acoustic tomography problem based on it

The paper considers the derivation of the wave equation and Helmholtz equation for solving the tomographic problem of reconstruction combined scalar-vector inhomogeneities describing perturbations of the sound velocity and absorption, the vector field of flows, and perturbations of the density of the medium. Restrictive conditions under which the obtained equations are meaningful are analyzed. Results of numerical simulation of the two-dimensional functional-analytical Novikov–Agaltsov algorithm for reconstructing the flow velocity using the the obtained Helmholtz equation are presented.     

Phase velocities of plane waves in a pipe with a flow whose velocity varies along the radius

The phase velocities of plane waves in a pipe filled with a moving acoustic medium are studied for different laws of flow velocity variation along the pipe radius. The wave equation is solved by the discretization method, which breaks the entire pipe volume into individual cylinders under the assumption that, within each of the cylinders, the flow velocity of the medium is constant. This approach makes it possible to reduce the solution to the wave problem to solving Helmholtz equations for individual cylinders. Based on boundary conditions satisfied at the boundaries between neighboring cylinders, a homogeneous system of linear algebraic equations is obtained. From this system, with the use of the scattering matrices, a simple dispersion equation is derived for determining the phase velocities of plane waves. The stability of the numerical solution to the dispersion equation with respect to the number of cylinders is investigated. The phase velocities of quasi-homogeneous and inhomogeneous waves in a pipe are numerically calculated and analyzed for different velocities of a moving medium and different laws of flow velocity variation along the radius. It is shown that the variation that occurs in the phase velocity of a homogeneous plane wave in a pipe due to the motion of the medium is identical to the mean flow velocity for different laws of flow velocity variation along the radius. For inhomogeneous plane waves, the phase velocity increment exceeds the mean flow velocity several times and depends on both the law of wave amplitude distribution along the radius and the law of the flow velocity variation along the radius.     

用变分法近似计算椭球马克中的倾斜模不稳定性

本文利用变分法的基本概念,提出了一种在任意边界形状的情况下近似求解偏微分方程的方法。把求解偏微分方程的问题近似地化为求解常微分方程组的问题。尤其是当边界形状较为复杂时,其优越性更为显著。作者采用这种方法,对大椭圆度的球形马克(Spheromak)中的倾斜模(Tilting Mode)不稳定性进行了研究。得到的结论是:在扁椭球马克中,不会出现倾斜模不稳定性。 关键词：     

Acoustic cavity modes in lens-shaped structures

A method for solving exactly the Helmholtz equation in parabolic rotational coordinates is presented using separability of the eigenfunctions and the Frobenius power series expansion technique. Two examples of interest in acoustics are considered and analyzed quasianalytically: The acoustic pressure in a cavity defined by two paraboloids (forming a lens-shaped structure) with (I) rigid wall boundary conditions and (II) pressure-release boundaries. The rigid-wall (pressure-release) acoustic enclosure problem is a Neumann (Dirichlet) boundary condition problem. In both cases, eigenfunctions and eigenmodes are calculated and the shape dependence of the eigenvalue for the ground state is examined.     

小尺度阻尼边界封闭空间声传递函数的有限元素法计算

如何求解阻尼边界封闭空间中声源点到接收点的低频声传递函数已成为目前小尺度封闭空间可听化技术研究的关键技术，能处理任意形状及复杂边界条件的有限元素法可作为求解该问题的适合方法，以室内声声有源Helmholtz方程及其相应边界方程为基础，本文推导出了用于小尺度阻尼边界封闭空间声传递函数的有限元素求解方法，并编制了相应的计算机程序，在算例中，首先通过与模态叠加法计算结果进行比较，验证了该方法的正确性。最后计算了某型车体内腔中任意两点间声传递函数。     

Solution of Schrödinger and related equations for irregular and composite regions

The Schrödinger equation with a potential is mathematically equivalent to the Helmholtz equation with a spatially variable propagation constant. A new method is presented for solving certain standard problems associated with these equations. In the Schrödinger language these are the ones in which a potential ν(|r|) acts inside an irregular (aspherical) boundary, where r has its origin inside the boundary. In terms of the Helmholtz equation, these include problems in which a region of constant index of refraction, but arbitrary shape, is embedded in a second uniform region with a different index. It is shown how bound state and scattering problems for such a potential (or region) can be treated in a way that avoids the usually intractable problem of matching solutions across the irregular boundary. The method requires, in general, the truncation of an infinite set of equations for partial wave amplitudes. The special case is discussed of a potential that becomes infinite throughout a region, so the wave amplitude must vanish inside the region (and, hence, on its boundary). For a long wave length this becomes a problem with the Laplace equation, and the general technique is illustrated by a calculation of the free charge on a perfectly conducting spheroid. The theory is extended from a single potential to an ensemble of such potentials, and in particular to an ensemble of potentials with spherical boundaries. In the special case that the potentials are arranged in a periodic lattice the formulas resemble those obtained by the KKR method, but are simpler in some ways. The method is extended to an ensemble of irregular potentials, and these results are shown to be applicable to the special case of an ensemble of finite range, but overlapping, spherical potentials.     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势,     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势, 具有良好的工程应用前景.     

Helmholtz equation and non-singular boundary elements applied to multi-disciplinary physical problems

The famous scientist Hermann von Helmholtz was born 200 years ago. Many complex physical wave phenomena in engineering can effectively be described using one or a set of equations named after him: the Helmholtz equation. Although this has been known for a long time, from a theoretical point of view, the actual numerical implementation has often been hindered by divergence-free and/or curl-free constraints. There is further a need for a numerical method that is accurate, reliable and takes into account radiation conditions at infinity. The classical boundary element method satisfies the last condition, yet one has to deal with singularities in the implementation. We review here how a recently developed singularity-free three-dimensional boundary element framework with superior accuracy can be used to tackle such problems only using one or a few Helmholtz equations with higher order (quadratic) elements which can tackle complex curved shapes. Examples are given for acoustics (a Helmholtz resonator among others) and electromagnetic scattering.     

Bistatic electromagnetic scattering from a three-dimensional perfect electric conducting object above a Gaussian rough surface based on the Kirchhoff–Helmholtz and electric field integral equation

A hybrid integral equation is developed to solve the problem of electromagnetic (EM) scattering from a three-dimensional (3D) perfect electric conducting (PEC) object above a two-dimensional (2D) PEC or dielectric Gaussian rough surface. Firstly, the Kirchhoff–Helmholtz (KH) equation is adopted to describe the wave reflection on the rough surface; only one integral operation on the rough surface is needed, and the scattering from the object can be described by solving the electric field integral equation (EFIE) on the surface of the object. Moreover, according to scattering theory, the KH equation and the EFIE are coupled together (KH-EFIE) to describe wave propagation between the object and the rough surface. Then method of moments (MoM) is adopted to solve the KH-EFIE, and the current is obtained to calculate the scattering field. Finally, compared with other methods, the accuracy of the proposed approach is validated, and its efficiency is proved to be much higher than numerical solutions. Furthermore, by calculating the statistic composite radar cross-section (RCS) of the object/surface and the difference radar cross-section (DRCS) of the object, the influence of the rough surface root mean square (rms) height, the correlation length, the medium permittivity, the shape of the object, and the altitude of the object on the scattering characteristic is investigated.