Time-harmonic Acoustic Waves Research Articles

Simulating time-harmonic acoustic wave effects induced by periodic holes/inclusions on surfaces

This paper introduces a localized meshless method to analyze time-harmonic acoustic wave propagation on curved surfaces with periodic holes/inclusions. In particular, the generalized finite difference method is used as a localized meshless technique to discretize the surface gradient and Laplace-Beltrami operators defined extrinsically in the governing equations. An absorbing boundary condition is introduced to reduce reflections from boundaries and accurately simulate wave propagation on unclosed surfaces with periodic inclusions. Several benchmark examples demonstrate the efficiency and accuracy of the proposed method in simulating acoustic wave propagation on surfaces with diverse geometries, including complex shapes and periodic holes or inclusions.

Learning scattering waves via coupling physics-informed neural networks and their convergence analysis

Time harmonic acoustic wave scattering problems modeled by the unbounded Helmholtz equations are widely used in medical and military fields. A fast and flexible solver for the Helmholtz equation is expected in engineering applications. Fortunately, data-driven neural network algorithms show extraordinary promise in scientific computing and artificial intelligence. We solve the Helmholtz equation by decomposing the equation into a coupling equation system that explicitly includes the real and imaginary parts of the solution. The real and imaginary parts are approximated by the different neural networks, respectively. The unbounded domain is reduced to a problem on a bounded domain via the Dirichlet-to-Neumann operator. Thus, a coupling physics-informed neural network (cPINN) is constructed. On the other hand, in exploring the convergence analysis for PINNs, we obtain the convergence results for cPINNs through the universal approximation and Sobolev space theory instead of using the neural tangent kernel theory. We prove that the cPINNs converge to the solution of the Helmholtz equation in H01 space as the number of neurons approaches infinity. Finally, we propose some numerical results.

Fast rapidly convergent penetrable scattering computations

We present a fast high-order scheme for the numerical solution of a volume-surface integro-differential equation. Such equations arise in problems of scattering of time-harmonic acoustic and electromagnetic waves by inhomogeneous media with variable density wherein the material properties jump across the medium interface. The method uses a partition of unity to segregate the interior and the boundary regions of the scattering obstacle, enabling us to make use of specially designed quadratures to deal with the material discontinuities in a high-order manner. In particular, the method uses suitable changes of variables to resolve the singularities present in the integrals in conjunction with a decomposition of Green’s function via the addition theorem. To achieve a reduced computational cost, the method employs a Fast Fourier Transform (FFT) based acceleration strategy to compute the integrals over the boundary region. Moreover, the necessary offgrid evaluation of the density and the inter-grid transfer of data is achieved by applying an FFT-based refined-grid interpolation strategy. We validate the performance of the method through multiple scattering simulations. In particular, the numerical experiments demonstrate that the proposed method can handle high-contrast material properties without any adverse effect on the number of GMRES iterations.

The factorization method for a penetrable cavity scattering with interior near-field measurements

Abstract This paper is concerned with the inverse scattering problem of time-harmonic acoustic waves from a penetrable cavity bounded by a layered structure and seeks to determine the shape and location of the cavity from interior near-field measurements. Of particular interest is that the near-field operator does not satisfy the main theorem of the factorization method, so we introduce a modified near-field operator and prove that it can be used to reconstruct the cavity. Numerical examples demonstrate the feasibility and effectiveness of our algorithm.

On the Far Field Pattern for Acoustic Scattering by a Piecewise Homogeneous Obstacle in Two Dimensions

A time-harmonic plane acoustic wave is scattered by a piecewise homogeneous obstacle with a penetrable or impenetrable core. We construct in the close form an integral representation for the far field pattern in which we have incorporated the physical and geometrical characteristics of the scatterer. Through this representation, we obtain the far field pattern for this scatterer. We prove scattering relations between the far field patterns of two scattering problems due to two distinct incident waves on the same scatterer. In particular, we prove reciprocity and general scattering theorems. The optical theorem, connecting the total power that the scatterer extracts from the incident plane wave either by radiation or by absorption with the corresponding far field pattern of an incident plane wave, is recovered as a corollary of the general scattering theorem. Moreover, if we consider incident waves to be both a plane and a spherical, we derive a mixed reciprocity theorem. We define the corresponding far field operators and using these relations, we prove some properties that can be used for solving inverse scattering problems.

Topological sensitivity analysis revisited for time-harmonic wave scattering problems. Part II: recursive computations by the boundary integral equation method

PurposeThe purpose of this paper is to revisit the recursive computation of closed-form expressions for the topological derivative of shape functionals in the context of time-harmonic acoustic waves scattering by sound-soft (Dirichlet condition), sound-hard (Neumann condition) and isotropic inclusions (transmission conditions).Design/methodology/approachThe elliptic boundary value problems in the singularly perturbed domains are equivalently reduced to couples of boundary integral equations with unknown densities given by boundary traces. In the case of circular or spherical holes, the spectral Fourier and Mie series expansions of the potential operators are used to derive the first-order term in the asymptotic expansion of the boundary traces for the solution to the two- and three-dimensional perturbed problems.FindingsAs the shape gradients of shape functionals are expressed in terms of boundary integrals involving the boundary traces of the state and the associated adjoint field, then the topological gradient formulae follow readily.Originality/valueThe authors exhibit singular perturbation asymptotics that can be reused in the derivation of the topological gradient function in the iterated numerical solution of any shape optimization or imaging problem relying on time-harmonic acoustic waves propagation. When coupled with converging Gauss−Newton iterations for the search of optimal boundary parametrizations, it generates fully automatic algorithms.

Topological sensitivity analysis revisited for time-harmonic wave scattering problems. Part I: the free space case

PurposeIn this paper, the authors revisit the computation of closed-form expressions of the topological indicator function for a one step imaging algorithm of two- and three-dimensional sound-soft (Dirichlet condition), sound-hard (Neumann condition) and isotropic inclusions (transmission conditions) in the free space.Design/methodology/approachFrom the addition theorem for translated harmonics, explicit expressions of the scattered waves by infinitesimal circular (and spherical) holes subject to an incident plane wave or a compactly supported distribution of point sources are available. Then the authors derive the first-order term in the asymptotic expansion of the Dirichlet and Neumann traces and their surface derivatives on the boundary of the singular medium perturbation.FindingsAs the shape gradient of shape functionals are expressed in terms of boundary integrals involving the boundary traces of the state and the associated adjoint field, then the topological gradient formulae follow readily.Originality/valueThe authors exhibit singular perturbation asymptotics that can be reused in the derivation of the topological gradient function that generates initial guesses in the iterated numerical solution of any shape optimization problem or imaging problems relying on time-harmonic acoustic wave propagation.

A Source Problem for the Helmholtz Equation via a Dirichlet-to-Neumann Map

In this paper, we consider a source problem for a time harmonic acoustic wave in two-dimensional space. Based on the boundary integral equation method, a Dirichlet-to-Neumann map in terms of boundary integral operators on the boundary of the source is constructed to transform this problem into two boundary value problems for the Helmholtz equation.

A Bayesian method for an inverse transmission scattering problem in acoustics

In this paper, we study an inverse transmission scattering problem of a time-harmonic acoustic wave from the viewpoint of Bayesian statistics. In Bayesian inversion, the solution of the inverse problem is the posterior distribution of the unknown parameters conditioned on the observational data. The shape of the scatterer will be reconstructed from full-aperture and limited-aperture far-field measurement data. We first prove a well-posedness result for the posterior distribution in the sense of the Hellinger metric. Then, we employ the Markov chain Monte Carlo method based on the preconditioned Crank-Nicolson algorithm to extract the posterior distribution information. Numerical results are given to demonstrate the effectiveness of the proposed method.

Design of an acoustic energy distributor using thin resonant slits

We consider the propagation of time-harmonic acoustic waves in a device made of three unbounded channels connected by thin slits. The wavenumber is chosen such that only one mode can propagate. The main goal of this work is to present a device which can serve as an energy distributor. More precisely, the geometry is first designed so that for an incident wave coming from one channel, the energy is almost completely transmitted in the two other channels. Additionally, adjusting slightly two geometrical parameters, we can control the ratio of energy transmitted in the two channels. The approach is based on asymptotic analysis for thin slits around resonance lengths. We also provide numerical results to illustrate the theory.

On corners scattering stably and stable shape determination by a single far-field pattern

In this paper, we establish two sharp quantitative results for the direct and inverse time-harmonic acoustic wave scattering. The first one is concerned with the recovery of the support of an inhomogeneous medium, independent of its contents, by a single far-field measurement. For this challenging inverse scattering problem, we establish a sharp stability estimate of logarithmic type when the medium support is a polyhedral domain in $\mathbb{R}^n$, $n=2,3$. The second one is concerned with the stability for corner scattering. More precisely if an inhomogeneous scatterer, whose support has a corner, is probed by an incident plane-wave, we show that the energy of the scattered far-field possesses a positive lower bound depending only on the geometry of the corner and bounds on the refractive index of the medium there. This implies the impossibility of approximate invisibility cloaking by a device containing a corner and made of isotropic material. Our results sharply quantify the qualitative corner scattering results in the literature, and the corresponding proofs involve much more subtle analysis and technical arguments. As a significant byproduct of this study, we establish a quantitative Rellich's theorem that continues smallness of the wave field from the far-field up to the interior of the inhomogeneity. The result is of significant mathematical interest for its own sake and is surprisingly not yet known in the literature.

The method of finite spheres in acoustic wave propagation through nonhomogeneous media: Inf-sup stability conditions

When the method of finite spheres is used for the solution of time-harmonic acoustic wave propagation problems in nonhomogeneous media, a mixed (or saddle-point) formulation is obtained in which the unknowns are the pressure fields and the Lagrange multiplier fields defined at the interfaces between the regions with distinct material properties. Then certain inf-sup conditions must be satisfied by the discretized spaces in order for the finite-dimensional problems to be well-posed. We discuss in this paper the analysis and use of these conditions. Since the conditions involve norms of functionals in fractional Sobolev spaces, we derive ‘stronger’ conditions that are simpler in form. These new conditions pave the way for the inf-sup testing, a tool for assessing the stability of the discretized problems.

The fluid-solid interaction scattering problem with unknown buried objects

Abstract Consider the direct and inverse scattering problem of time harmonic acoustic waves by a two-dimensional elastic obstacle which contains an unknown impenetrable object inside. We apply the boundary integral equation method to solve the direct scattering problem. Since the obtained boundary integral system is a mixed form of scalar and vector equations, we consider the existence of the solution in the sense of a weak formulation. The factorization method is employed to solve the corresponding inverse scattering problem, that is, we try to recover the shape of the interaction surface from the knowledge of acoustic far-field patterns. With the help of a special data to pattern operator G, we can reconstruct the shape of the interaction surface by using the classical factorization method. Furthermore, a new mixed reciprocity relation is proposed to show the relation between the interior elastic field and unitary matrix 𝒮 {\mathcal{S}} .

Nonlinear Interaction of Acoustic Waves with a Spheroidal Particle: Radiation Force and Torque Effects

The nonlinear interaction of a time-harmonic acoustic wave with an anisotropic particle gives rise to the radiation force and torque effects. These phenomena are at the heart of the acoustofluidics technology, where microparticles such as cells and microorganisms are acoustically manipulated. We present a theoretical model considering a generic acoustic beam interacting with a subwavelength spheroidal particle in a nonviscous fluid. Concise analytical expressions of the radiation force and torque are obtained in the scattering dipole approximation. The radiation force is given in terms of a gradient and scattering force; while the radiation torque has two fundamental contributions, namely, the momentum arm and acoustic spin (spin-torque effect). As a practical example, we use the theory to describe the interaction of two crossed plane waves and a prolate spheroidal particle. The results reveal the particle is transversely trapped in a pressure node and is axially pushed by the radiation force. Also, the momentum arm aligns the particle in the axial direction. At certain specific positions, only the spin-torque occurs. Our findings are remarkably consistent with finite-element simulations. The success of our model enables its use as an investigation tool for the manipulation of anisotropic microparticles in acoustofluidics.

Inverse random source problems for time-harmonic acoustic and elastic waves

This paper concerns the random source problems for the time-harmonic acoustic and elastic wave equations in two and three dimensions. The goal is to determine the compactly supported external force from the radiated wave field measured in a domain away from the source region. The source is assumed to be a microlocally isotropic generalized Gaussian random function such that its covariance operator is a classical pseudo-differential operator. Given such a distributional source, the direct problem is shown to have a unique solution by using an integral equation approach and the Sobolev embedding theorem. For the inverse problem, we demonstrate that the amplitude of the scattering field averaged over the frequency band, obtained from a single realization of the random source, determines uniquely the principle symbol of the covariance operator. The analysis employs asymptotic expansions of the Green functions and microlocal analysis of the Fourier integral operators associated with the Helmholtz and Navier equations.

Unique identification of a multi-layered fluid\u2013solid medium

This paper is concerned with the inverse scattering of time-harmonic acoustic plane waves by a multi-layered fluid–solid medium in the three dimensional space. We establish the global uniqueness in identifying the embedded penetrable solid obstacle, the surrounding fluid medium and its wave number from the acoustic far-field pattern for all incident plane waves at a fixed frequency. The proof depends on constructing different kinds of interior transmission problems in appropriate small domains and the a priori estimates derived for both the elastic wave fields in the embedded solid obstacle and the acoustic wave fields in the surrounding fluid medium.

Least-Squares Pad\xe9 approximation of parametric and stochastic Helmholtz maps

The present work deals with rational model order reduction methods based on the single-point Least-Square (LS) Padé approximation techniques introduced in Bonizzoni et al. (ESAIM Math. Model. Numer. Anal., 52(4), 1261–1284 2018, Math. Comput. 89, 1229–1257 2020). Algorithmical aspects concerning the construction of rational LS-Padé approximants are described. In particular, we show that the computation of the Padé denominator can be carried out efficiently by solving an eigenvalue-eigenvector problem involving a Gramian matrix. The LS-Padé techniques are employed to approximate the frequency response map associated with two parametric time-harmonic acoustic wave problems, namely a transmission-reflection problem and a scattering problem. In both cases, we establish the meromorphy of the frequency response map. The Helmholtz equation with stochastic wavenumber is also considered. In particular, for Lipschitz functionals of the solution and their corresponding probability measures, we establish weak convergence of the measure derived from the LS-Padé approximant to the true one. 2D numerical tests are performed, which confirm the effectiveness of the approximation methods.

The factorization method to reconstruct a combination of penetrable obstacle and arc

ABSTRACT This paper considers the scattering of time-harmonic acoustic plane waves by a mixed scatterer. Such a scatterer is given as the union of a penetrable homogeneous medium with compact support and an open arc. Firstly, the well-posedness of the direct problem is established by the boundary integral method. Then, we study the factorization method for recovering the shape of the mixed scatterer. Furthermore, some numerical experiments are presented to demonstrate the feasibility and effectiveness of our approach.

An inverse acoustic-elastic interaction problem with phased or phaseless far-field data

Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper is concerned with an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to deduce the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integral equations is proposed for the phaseless inverse problem. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed methods.

Determining a Multi-layered Fluid-solid Medium from the Acoustic Measurements

Consider the inverse problem of recovering a multi-layered fluid-solid medium from many acoustic measurements corresponding to time-harmonic acoustic plane waves. We prove that both the supports of the embedded solid obstacle and the surrounding layered fluid medium can be uniquely identified by means of acoustic far-field pattern for all incident wave fields at a fixed frequency. Our proof is based on the constructions of some well-posed partial differential equation systems in sufficiently small domains combined with the a priori estimates for the solutions of the forward scattering problem.