A new method for true and spurious eigensolutions of arbitrary cavities using the combined Helmholtz exterior integral equation formulation method

Integral equation methods have been widely used to solve interior eigenproblems and exterior acoustic problems (radiation and scattering). It was recently found that the real-part boundary element method (BEM) for the interior problem results in spurious eigensolutions if the singular (UT) or the hypersingular (LM) equation is used alone. The real-part BEM results in spurious solutions for interior problems in a similar way that the singular integral equation (UT method) results in fictitious solutions for the exterior problem. To solve this problem, a Combined Helmholtz Exterior integral Equation Formulation method (CHEEF) is proposed. Based on the CHEEF method, the spurious solutions can be filtered out if additional constraints from the exterior points are chosen carefully. Finally, two examples for the eigensolutions of circular and rectangular cavities are considered. The optimum numbers and proper positions for selecting the points in the exterior domain are analytically studied. Also, numerical experiments were designed to verify the analytical results. It is worth pointing out that the nodal line of radiation mode of a circle can be rotated due to symmetry, while the nodal line of the rectangular is on a fixed position.     

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

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

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.     

波动方程基于自然边界归化的区域分解算法

提出了无界区域波动方程的区域分解算法,基于自然边界归化,分别研究了重叠型与非重叠型区域分解算法,首先将控制方程对时间进行离散化,得到关于时间步长离散化格式,对每一时间步长给出了Dirichlet-Neumann和Schwartz交替算法,对Schwartz交替算法,给出了算法的收敛性,对圆外区域研究了压缩因子,并给出了数值例子。     

Two transparent boundary conditions for the electromagnetic scattering from two-dimensional overfilled cavities

We consider electromagnetic scattering from two-dimensional (2D) overfilled cavities embedded in an infinite ground plane. The unbounded computational domain is truncated to a bounded one by using a transparent boundary condition (TBC) proposed on a semi-ellipse. For overfilled rectangular cavities with homogeneous media, another TBC is introduced on the cavity apertures, which produces a smaller computational domain. The existence and uniqueness of the solutions of the variational formulations for the transverse magnetic and transverse electric polarizations are established. In the exterior domain, the 2D scattering problem is solved in the elliptic coordinate system using the Mathieu functions. In the interior domain, the problem is solved by a finite element method. Numerical experiments show the efficiency and accuracy of the new boundary conditions.     

Transformation Electrics: Cloaking and Rotating Electric Current

Transformation optics provides great versatility for precisely manipulating electromagnetic waves. It has been extended to other fields including acoustics, thermotics, and electrics. Taking advantage of the transformation electrics method, we demonstrate that the square-shaped cloak can guide electric current around the cloaked region smoothly without perturbing the exterior electric current. And the cylindrical rotator can rotate the electric current.Inside the enclosed domain of the rotator, the electric current from the outside will appear as if it is coming from a different angle. Finally, the related experimental realizations and potential applications are also discussed.     

抛物方程基于自然边界归化的耦合法

将冯康和余德浩提出的自然边界归化方法^[1~4]应用于求解抛物方程初边值外区域问题,提出一种自然边界元与有限元耦合算法。先将控制方程对时间进行离散化,得到关于时间步长的离散化格式,给出圆外域上的自然积分方程,基于此研究抛物方程无界区域问题的自然边界元与有限元耦合法,最后给出相应的数值例子。     

A New Numerical Method for Solving the Acoustic Radiation Problem

A numerical method of solving the problem of acoustic wave radiation in the presence of a rigid scatterer is described. It combines the finite element method and the boundary algebraic equation one. In the proposed method, the exterior domain around the scatterer is discretized, so that there appear an infinite domain with regular discretization and a relatively small layer with irregular mesh. For the infinite regular mesh, the boundary algebraic equation method is used with spurious resonance suppression according to Burton and Miller. In the thin layer with irregular mesh, the finite element method is used. The proposed method is characterized by simple implementation, fair accuracy, and absence of spurious resonances.     

A hybrid finite element approach to modeling sound radiation from circular and rectangular ducts

A numerical model based on a hybrid finite element method is developed that seeks to join sound pressure fields in interior and exterior regions. The hybrid method is applied to the analysis of sound radiation from open pipes, or ducts, and uses mode matching to couple a finite element discretization of the region surrounding the open end of the duct to wave based modal expansions for adjoining interior and exterior regions. The hybrid method facilitates the analysis of ducts of arbitrary but uniform cross section as well the study of conical flanges and here a modal expansion based on spherical harmonics is applied. Predictions are benchmarked against analytic solutions for the limiting cases of flanged and unflanged circular ducts and excellent agreement between the two methods is observed. Predictions are also presented for flanged and unflanged rectangular ducts, and because the hybrid method retains the sparse banded and symmetric matrices of the traditional finite element method, it is shown that predictions can be obtained within an acceptable time frame even for a three dimensional problem.     

Data completion method for the characterization of sound sources

This paper presents a different approach to solve the inverse acoustic problem. This problem is an "ill-posed" problem since the solution is very sensitive to measurement precision. A classical way to solve this problem consists in inversing a propagation operator which relates structure quantities (acoustic pressures or gradients) to near-field quantities (acoustic pressures or gradients). This can be achieved by using near-field acoustical holography (NAH) in separable coordinate systems. In order to overcome this limitation, the inverse boundary element method (IBEM) can be implemented to recover all acoustic quantities in a three-dimensional space and on an arbitrary three-dimensional source surface. In this paper, the data completion method (DCM) is developed: the acoustic gradients and pressures are known on a surface surrounding the source, but are unknown on its structure. The solution is given by the resolution of the Helmholtz formulation applied on the empty domain between the two boundaries made by the measurements quantities and the structure of the source. The conventional method applies directly the integral formulation for the empty domain. Another way of solving this Helmholtz formulation can be achieved by splitting it in two well-posed subproblems in a Steklov-Poincare?'s formulation. The data completion method allows one to solve the problem with acoustic perturbations due to sources on the exterior domain, or due to a confined domain, without altering the results.     

Scattering of Gaussian beam by arbitrarily shaped inhomogeneous particles

A hybrid finite element-boundary integral method is applied to characterize the scattering of an arbitrarily incident-focused Gaussian beam by arbitrarily shaped inhomogeneous particles. Specifically, the Davis–Barton fifth-order approximation in combination with rotation Euler angles is used to represent the arbitrarily incident Gaussian beams. The finite element method is employed to formulate the fields in the interior region of the inhomogeneous particle, while the boundary integral equation is applied to represent the fields in the exterior region. The interior and exterior fields are coupled by means of the field continuity conditions. To reduce the computational burden, the frontal method and the multilevel fast multipole algorithm are adopted to solve the resultant matrix equation. Numerical results for differential scattering cross sections of several selected inhomogeneous particles are presented and can be served as further study on this subject.     

Formal iteration scheme for scalar fields in some spherically symmetric spaces

The characteristic initial-value problem for the scalar wave equation in some spherically symmetrical geometries, including Reissner-Nordström exterior geometry, is considered. A modified Picard iteration making explicit use of the (known) static solutions yields a formula for the fields. The formulas extend in to the outer horizon, and from them are obtained useful expressions for the reflection and transmission amplitudes in terms of the static solutions.     

Incompressible Flow Around a Small Obstacle and the Vanishing Viscosity Limit

In this article we consider viscous flow in the exterior of an obstacle satisfying the standard no-slip boundary condition at the surface of the obstacle. We seek conditions under which solutions of the Navier-Stokes system in the exterior domain converge to solutions of the Euler system in the full space when both viscosity and the size of the obstacle vanish. We prove that this convergence is true assuming two hypotheses: first, that the initial exterior domain velocity converges strongly in L ² to the full-space initial velocity and second, that the diameter of the obstacle is smaller than a suitable constant times viscosity, or, in other words, that the obstacle is sufficiently small. The convergence holds as long as the solution to the limit problem is known to exist and stays sufficiently smooth. This work complements the study of incompressible flow around small obstacles, which has been carried out in [4–6].     

Acoustoelasticity: General theory,acoustic natural modes and forced response to sinusoidal excitation,including comparisons with experiment

A comprehensive theoretical model has been developed for interior sound fields which are created by flexible wall motion resulting from exterior sound fields. Full coupling between the wall and interior acoustic cavity is permitted. An efficient computational method is used to determine acoustic natural frequencies of multiply connected cavities. Simplified formulae are developed for interior sound levels in terms of cavity, wall and external acoustic field parameters. Comparisons of theory and experiment show generally good agreement.     

Transient acoustic radiation from impulsively accelerated bodies by the finite element method

Bodies under impulsive motion, immersed in an infinite acoustic fluid, severely put to test any numerical method for the transient exterior acoustic problem. Such problems, in the context of the finite element method (FEM), are not well studied. FE modeling of such problems requires truncation of the infinite fluid domain at a certain distance from the structure. The volume of computation depends upon the extent of this domain as well as the mesh density. The modeling of the fluid truncation boundary is crucial to the economy and accuracy of solution and various boundary dampers have been proposed in the literature for this purpose. The second order damper leads to unsymmetric boundary matrices and this necessitates the use of an unsymmetric equation solver for large problems. The present paper demonstrates the use of FEM with zeroth, first and second order boundary dampers in conjunction with an unsymmetric, out of core, banded equation solver for impulsive motion problems of rigid bodies in an acoustic fluid. The results compare well with those obtained from analytical methods.     

Acoustic-elastodynamic interaction in isotropic fractal media

This research explores the acoustic-elastodynamic interaction in isotropic fractal media. The analysis discusses the direct coupling of two constitutive models under dynamic loading: a continuous solid and an isotropic fractal medium. We consider two situations where in the first, the fractal medium is enclosed within a thin spherical shell (interior problem), while in the second, the fractal medium extends infinitely outside the shell (exterior problem). The two problems are simulated analytically, and the exact solution for the shell displacement is expressed in closed form in the Laplace domain. The formulation of the radiation condition for infinite fractal media is essential to derive the exterior problem’s solution. This study represents a meaningful idealization of real-application problems involving the interaction of multi-constitutive media, e.g. the human brain, whereby fractal features affect the response of this body under various excitations.     

Relarivistic gravitational fields of rotating bodies

The Einstein equations for the exterior fields of stationarily rotating axisymmetric sources are solved by the inverse (scattering) method, i.e. the metric can be calculated from the “scattering” coefficients via a Zakharov-Shabat equation or a Gelfand-Levitan-Marchenko equation.     

Finite element model for waves guided along solid systems of arbitrary section coupled to infinite solid media

The Semi-Analytical Finite Element (SAFE) method is becoming established as a convenient method to calculate the properties of waves which may propagate in a waveguide which has arbitrary cross-sectional shape but which is invariant in the propagation direction. A number of researchers have reported work relating to lossless elastic waves, and recently the solutions for nonpropagating waves in elastic guides and for complex waves in viscoelastic guides have been presented. This paper presents a further development, addressing the problem of attenuating waves in which the attenuation is caused by leakage from the waveguide into a surrounding material. This has broad relevance to many practical problems in which a waveguide is immersed in a fluid or embedded in a solid. The paper presents the principles of a procedure and then validates and illustrates its use on some examples. The procedure makes use of absorbing regions of material at the exterior bounds of the discretized domain.     

Application of Padé via Lanczos approximations for efficient multifrequency solution of Helmholtz problems

This paper addresses the efficient solution of acoustic problems in which the primary interest is obtaining the solution only on restricted portions of the domain but over a wide range of frequencies. The exterior acoustics boundary value problem is approximated using the finite element method in combination with the Dirichlet-to-Neumann (DtN) map. The restriction domain problem is formally posed in transfer function form based on the finite element solution. In order to obtain the solution over a range of frequencies, a matrix-valued Padé approximation of the transfer function is employed, using a two-sided block Lanczos algorithm. This approach provides a stable and efficient representation of the Padé approximation. In order to apply the algorithm, it is necessary to reformulate the transfer function due to the frequency dependency in the nonreflecting boundary condition. This is illustrated for the case of the DtN boundary condition, but there is no restriction on the approach which can also be applied to other radiation boundary conditions. Numerical tests confirm that the approach offers significant computational speed-up.     

Calculation of electromagnetic fields in cavity resonators with allowance for energy flux through slots

By the example of a cylindrical cavity with an annular slot, resonance conditions in cavity resonators with slots in the perfectly conducting walls are theoretically studied. The stationary electromagnetic fields of eigenmodes are described in terms of the partial domain method with the use of the homogeneous Maxwell equations. The complex wavenumber and other parameters of these fields are determined from the condition that the fields are continuous on the slot boundaries. An energy resonance condition is formulated that requires that the energy flux of the resonance field through the slot be minimal. It is shown that these conditions, taken together, remove infeasible (nonphysical) solutions to the homogeneous problem, which are associated with out-of-cavity excitation sources.