More Like this

1. An algorithm for computing scattering poles based on dual characterization to interior eigenvalues

   https://doi.org/10.1098/rspa.2024.0015  

   Cakoni, Fioralba; Haddar, Houssem; Zilberberg, Dana (June 2024, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences)

   We present an algorithm for the computation of scattering poles for an impenetrable obstacle with Dirichlet or Robin boundary conditions in acoustic scattering. This paper builds upon the previous work of Cakoni et al. (2020) titled ‘A duality between scattering poles and transmission eigenvalues in scattering theory’ (Cakoni et al. 2020 Proc. A. 476, 20200612 (doi:10.1098/rspa.2020.0612)), where the authors developed a conceptually unified approach for characterizing the scattering poles and interior eigenvalues corresponding to a scattering problem. This approach views scattering poles as dual to interior eigenvalues by interchanging the roles of incident and scattered fields. In this framework, both sets are linked to the kernel of the relative scattering operator that maps incident fields to scattered fields. This mapping corresponds to the exterior scattering problem for the interior eigenvalues and the interior scattering problem for scattering poles. Leveraging this dual characterization and inspired by the generalized linear sampling method for computing the interior eigenvalues, we present a novel numerical algorithm for computing scattering poles without relying on an iterative scheme. Preliminary numerical examples showcase the effectiveness of this computational approach. 

   more » « less
2. Multiple-scattering frequency-time hybrid solver for the wave equation in interior domains

   https://doi.org/10.1090/mcom/3872  

   Bruno, Oscar; Yin, Tao (March 2024, Mathematics of Computation)

   Brenner, Susan (Ed.)

   This paper proposes a frequency-time hybrid solver for the time-dependent wave equation in two-dimensionalinterior spatial domains. The approach relies on four main elements, namely, (1) A multiple scattering strategy that decomposes a giveninteriortime-domain problem into a sequence oflimited-durationtime-domain problems of scattering by overlapping open arcs, each one of which is reduced (by means of the Fourier transform) to a sequence ofHelmholtz frequency-domain problems; (2) Boundary integral equations on overlapping boundary patches for the solution of the frequency-domain problems in point (1); (3) A smooth“Time-windowing and recentering”methodology that enables both treatment of incident signals of long duration and long time simulation; and, (4) A Fourier transform algorithm that delivers numerically dispersionless,spectrally-accurate time evolutionfor given incident fields. By recasting the interior time-domain problem in terms of a sequence of open-arc multiple scattering events, the proposed approach regularizes the full interior frequency domain problem—which, if obtained by either Fourier or Laplace transformation of the corresponding interior time-domain problem, must encapsulate infinitely many scattering events, giving rise to non-uniqueness and eigenfunctions in the Fourier case, and ill conditioning in the Laplace case. Numerical examples are included which demonstrate the accuracy and efficiency of the proposed methodology. 

   more » « less
3. Sum rules for energy deposition eigenchannels in scattering systems

   https://doi.org/10.1364/OL.468697  

   Yamilov, Alexey; Bender, Nicholas; Cao, Hui (September 2022, Optics Letters)

   In a random-scattering system, the deposition matrix maps the incident wavefront onto the internal field distribution across a target volume. The corresponding eigenchannels have been used to enhance the wave energy delivered to the target. Here, we find the sum rules for the eigenvalues and eigenchannels of the deposition matrix in any system geometry: including two- and three-dimensional scattering systems, as well as narrow waveguides and wide slabs. We derive a number of constraints on the eigenchannel intensity distributions inside the system as well as the corresponding eigenvalues. Our results are general and applicable to random systems of arbitrary scattering strength as well as different types of waves including electromagnetic waves, acoustic waves, and matter waves. 

   more » « less
4. Convergence of an adaptive finite element DtN method for the elastic wave scattering by periodic structures

   https://doi.org/10.1016/j.cma.2019.112722  

   Li, P.; Yuan, X. (January 2020, Computer methods in applied mechanics and engineering)

   Consider the scattering of a time-harmonic elastic plane wave by a periodic rigid surface. The elastic wave propagation is governed by the two-dimensional Navier equation. Based on a Dirichlet-to-Neumann (DtN) map, a transparent boundary condition (TBC) is introduced to reduce the scattering problem into a boundary value problem in a bounded domain. By using the finite element method, the discrete problem is considered, where the TBC is replaced by the truncated DtN map. A new duality argument is developed to derive the a posteriori error estimate, which contains both the finite element approximation error and the DtN truncation error. An a posteriori error estimate based adaptive finite element algorithm is developed to solve the elastic surface scattering problem. Numerical experiments are presented to demonstrate the effectiveness of the proposed method. 

   more » « less
5. A spectral target signature for thin surfaces with higher order jump conditions

   https://doi.org/10.3934/ipi.2022020  

   Cakoni, Fioralba; Lee, Heejin; Monk, Peter; Zhang, Yangwen (January 2022, Inverse Problems and Imaging)

   In this paper we consider the inverse problem of determining structural properties of a thin anisotropic and dissipative inhomogeneity in \begin{document}$$ {\mathbb R}^m $$\end{document}, \begin{document}$ m = 2, 3 $$\end{document} from scattering data. In the asymptotic limit as the thickness goes to zero, the thin inhomogeneity is modeled by an open \begin{document}$$ m-1 $$\end{document}$ dimensional manifold (here referred to as screen), and the field inside is replaced by jump conditions on the total field involving a second order surface differential operator. We show that all the surface coefficients (possibly matrix valued and complex) are uniquely determined from far field patterns of the scattered fields due to infinitely many incident plane waves at a fixed frequency. Then we introduce a target signature characterized by a novel eigenvalue problem such that the eigenvalues can be determined from measured scattering data, adapting the approach in [20]. Changes in the measured eigenvalues are used to identified changes in the coefficients without making use of the governing equations that model the healthy screen. In our investigation the shape of the screen is known, since it represents the object being evaluated. We present some preliminary numerical results indicating the validity of our inversion approach 

   more » « less