Oscar P. Bruno

Faisal Amlani, Oscar P. Bruno, José Carlos López-Vázquez et al.

We study the scattering of transient, high-frequency, narrow-band quasi-Rayleigh elastic waves by through-thickness holes in aluminum plates, in the framework of ultrasonic nondestructive testing (NDT) based on full-field optical detection. Sequences of the instantaneous two-dimensional (2-D) out-of-plane displacement scattering maps are measured with a self-developed PTVH system. The corresponding simulated sequences are obtained by means of an FC(Gram) elastodynamic solver introduced recently, which implements a full three-dimensional (3D) vector formulation of the direct linear-elasticity scattering problem. A detailed quantitative comparison between these experimental and numerical sequences, which is presented here for the first time, shows very good agreement both in the amplitude and the phase of the acoustic field in the forward, lateral and backscattering areas. It is thus suggested that the combination of the PTVH system and the FC(Gram) elastodynamic solver provides an effective ultrasonic inspection tool for plate-like structures, with a significant potential for ultrasonic NDT applications.

Oscar P. Bruno, Stephane K. Lintner

We present a novel approach for the numerical solution of problems of diffraction by open arcs in two dimensional space. Our methodology relies on composition of {\em weighted versions} of the classical integral operators associated with the Dirichlet and Neumann problems (TE and TM polarizations, respectively) together with a generalization to the open-arc case of the well known closed-surface Calderón formulae. When used in conjunction with spectrally accurate discretization rules and Krylov-subspace linear algebra solvers such as GMRES, the new second-kind TE and TM formulations for open arcs produce results of high accuracy in small numbers of iterations and short computing times---for low and high frequencies alike.

Oscar P. Bruno, Emmanuel Garza

This paper introduces a high-order-accurate strategy for integration of singular kernels and edge-singular integral densities that appear in the context of boundary integral equation formulations of the problem of acoustic scattering. In particular, the proposed method is designed for use in conjunction with geometry descriptions given by a set of arbitrary non-overlapping logically-quadrilateral patches---which makes the algorithm particularly well suited for treatment of CAD-generated geometries. Fejér's first quadrature rule is incorporated in the algorithm, to provide a spectrally accurate method for evaluation of contributions from far integration regions, while highly-accurate precomputations of singular and near-singular integrals over certain "surface patches" together with two-dimensional Chebyshev transforms and suitable surface-varying "rectangular-polar" changes of variables, are used to obtain the contributions for singular and near-singular interactions. The overall integration method is then used in conjunction with the linear-algebra solver GMRES to produce solutions for sound-soft open- and closed-surface scattering obstacles, including an application to an aircraft described by means of a CAD representation. The approach is robust, fast, and highly accurate: use of a few points per wavelength suffices for the algorithm to produce far-field accuracies of a fraction of a percent, and slight increases in the discretization densities give rise to significant accuracy improvements.

Oscar P. Bruno, Manuel A. Santana, Lloyd N. Trefethen

This paper presents a novel algorithm, based on use of rational approximants of a randomly scalarized boundary integral resolvent in conjunction with an adaptive search strategy and an exponentially convergent secant-method termination stage, for the evaluation of acoustic and electromagnetic resonances in open and closed cavities. The desired cavity resonances are obtained as the poles of associated rational approximants; both the approximants and their poles are obtained by means of the recently introduced AAA rational-approximation algorithm. In fact, the proposed resonance-search method applies to any nonlinear eigenvalue problem associated with a given function $F: U \to \mathbb{C}^{d\times d}$, wherein, denoting $F(k) = F_k$, a complex value $k$ is sought for which $F_kw = 0$ for some nonzero $w\in \mathbb{C}^d$. For the scattering problems considered in this paper, $F_k$ is taken to equal a spectrally discretized version of a Green function-based boundary integral operator at spatial frequency $k$. In all cases, the scalarized resolvent is given by an expression of the form $u^* F_k^{-1} v$, where $u,v \in \mathbb{C}^d$ are fixed random vectors. The proposed adaptive search strategy relies on use of a rectangular subdivision of the resonance search domain which is locally refined to ensure that all resonances in the domain are captured. The approach works equally well in the case in which the search domain is an interval of the real line, in which case the rectangles used degenerate into subintervals of the search domain. A variety of numerical results are presented, including comparisons with well-known methods based on complex contour integration, and a discussion of the asymptotics that result as open cavities approach closed cavities -- in all, demonstrating the accuracy provided by the method, for low- and high-frequency states alike.

Oscar P. Bruno, Stephen P. Shipman, Catalin Turc et al.

This paper, Part I in a two-part series, presents (i) A simple and highly efficient algorithm for evaluation of quasi-periodic Green functions, as well as (ii) An associated boundary-integral equation method for the numerical solution of problems of scattering of waves by doubly periodic arrays of scatterers in three-dimensional space. Except for certain "Wood frequencies" at which the quasi-periodic Green function ceases to exist, the proposed approach, which is based on use of smooth windowing functions, gives rise to lattice sums which converge superalgebraically fast--that is, faster than any power of the number of terms used--in sharp contrast with the extremely slow convergence exhibited by the corresponding sums in absence of smooth windowing. (The Wood-frequency problem is treated in Part II.) A proof presented in this paper establishes rigorously the superalgebraic convergence of the windowed lattice sums. A variety of numerical results demonstrate the practical efficiency of the proposed approach.

Oscar P. Bruno, Victor Dominguez, Francisco-Javier Sayas

In this paper we present a convergence analysis for the Nystrom method proposed in [Jour. Comput. Phys. 169 pp. 2921-2934, 2001] for the solution of the combined boundary integral equation formulations of sound-soft acoustic scattering problems in three-dimensional space. This fast and efficient scheme combines FFT techniques and a polar change of variables that cancels out the kernel singularity. We establish the stability of the algorithms in the $L^2$ norm and we derive convergence estimates in both the $L^2$ and $L^\infty$ norms. In particular, our analysis establishes theoretically the previously observed super-algebraic convergence of the method in cases in which the right-hand side is smooth.

Oscar P. Bruno, Agustin G. Fernandez-Lado

This paper presents a full-spectrum Green function methodology (which is valid, in particular, at and around Wood-anomaly frequencies) for evaluation of scattering by periodic arrays of cylinders of arbitrary cross section-with application to wire gratings, particle arrays and reflectarrays and, indeed, general arrays of conducting or dielectric bounded obstacles under both TE and TM polarized illumination. The proposed method, which, for definiteness is demonstrated here for arrays of perfectly conducting particles under TE polarization, is based on use of the shifted Green-function method introduced in the recent contribution (Bruno and Delourme, Jour. Computat. Phys. pp. 262--290 (2014)). A certain infinite term arises at Wood anomalies for the cylinder-array problems considered here that is not present in the previous rough-surface case. As shown in this paper, these infinite terms can be treated via an application of ideas related to the Woodbury-Sherman-Morrison formulae. The resulting approach, which is applicable to general arrays of obstacles even at and around Wood-anomaly frequencies, exhibits fast convergence and high accuracies. For example, a few hundreds of milliseconds suffice for the proposed approach to evaluate solutions throughout the resonance region (wavelengths comparable to the period and cylinder sizes) with full single-precision accuracy.

Oscar P. Bruno, Stephen P. Shipman, Catalin Turc et al.

This work presents an efficient method for evaluation of wave scattering by doubly periodic diffraction gratings at or near "Wood anomaly frequencies". At these frequencies, one or more grazing Rayleigh waves exist, and the lattice sum for the quasi-periodic Green function ceases to exist. We present a modification of this sum by adding two types of terms to it. The first type adds linear combinations of "shifted" Green functions, ensuring that the spatial singularities introduced by these terms are located below the grating and therefore outside of the physical domain. With suitable coefficient choices these terms annihilate the growing contributions in the original lattice sum and yield algebraic convergence. Convergence of arbitrarily high order can be obtained by including sufficiently many shifts. The second type of added terms are quasi-periodic plane wave solutions of the Helmholtz equation which reinstate certain necessary grazing modes without leading to blow-up at Wood anomalies. Using the new quasi-periodic Green function, we establish, for the first time, that the Dirichlet problem of scattering by a smooth doubly periodic scattering surface at a Wood frequency is uniquely solvable. We also present an efficient high-order numerical method based on the this new Green function for the problem of scattering by doubly periodic three-dimensional surfaces at and around Wood frequencies. We believe this is the first solver in existence that is applicable to Wood-frequency doubly periodic scattering problems. We demonstrate the proposed approach by means of applications to problems of acoustic scattering by doubly periodic gratings at various frequencies, including frequencies away from, at, and near Wood anomalies.

Oscar P. Bruno, Liwei Xu, Tao Yin

We present a novel approach for the numerical solution of problems of elastic scattering by open arcs in two dimensions. Our methodology relies on the composition of weighted versions of the classical operators associated with Dirichlet and Neumann boundary conditions in conjunction with a certain "open-arc elastic Calderón relation" whose validity is demonstrated in this paper on the basis of numerical experiments, but whose rigorous mathematical proof is left for future work. Using this Calderón relation in conjunction with spectrally accurate quadrature rules and the Krylov-subspace linear algebra solver GMRES, the proposed overall open-arc elastic solver produces results of high accuracy in small number of iterations---for low and high frequencies alike. A variety of numerical examples in this paper demonstrate the accuracy and efficiency of the proposed methodology.

Oscar P. Bruno, Carlos Pérez-Arancibia

This paper presents a new methodology for the solution of problems of two- and three-dimensional acoustic scattering (and, in particular, two-dimensional electromagnetic scattering) by obstacles and defects in presence an arbitrary number of penetrable layers. Relying on use of certain slow-rise windowing functions, the proposed Windowed Green Function approach (WGF) efficiently evaluates oscillatory integrals over unbounded domains, with high accuracy, without recourse to the highly expensive Sommerfeld integrals that have typically been used to account for the effect of underlying planar multi-layer structures. The proposed methodology, whose theoretical basis was presented in the recent contribution (SIAM J. Appl. Math. 76(5), p. 1871, 2016), is fast, accurate, flexible, and easy to implement. Our numerical experiments demonstrate that the numerical errors resulting from the proposed approach decrease faster than any negative power of the window size. In a number of examples considered in this paper the proposed method is up to thousands of times faster, for a given accuracy, than corresponding methods based on use of Sommerfeld integrals.

Oscar P. Bruno, Manuel A. Santana

Fourier transform-based methods enable accurate, dispersion-free simulations of time-domain scattering problems by evaluating solutions to the Helmholtz equation at a discrete set of frequencies sufficient to approximate the inverse Fourier transform. However, in the case of scattering by trapping obstacles, the Helmholtz solution exhibits nearly-real complex resonances -- which significantly slows the convergence of numerical inverse transform. To address this difficulty this paper introduces a frequency-domain singularity subtraction technique that regularizes the integrand of the inverse transform and efficiently computes the singularity contribution via a combination of a straightforward and inexpensive numerical technique together with a large-time asymptotic expansion. Crucially, all relevant complex resonances and their residues are determined via rational approximation of integral equation solutions at real frequencies. An adaptive algorithm is employed to ensure that all relevant complex resonances are properly identified.

Oscar P. Bruno, Allen Yang

This paper presents a fast "two-dimensional Fourier Continuation" (2D-FC) method for the construction of biperiodic extensions of smooth, non-periodic functions defined over general two-dimensional (2D) domains, including domains with corners. The algorithm operates with an O(N log N) computational cost, for an N-point discretization grid, and it achieves a user-prescribed d-th order of accuracy. The methodology can be generalized to non-smooth domains of arbitrary dimensionality, but such extensions are not considered in the present work. The usefulness and performance of the 2D-FC method are demonstrated through applications to the Poisson problem posed on bounded 2D domains with corners. One illustrative application concerns a Poisson problem on a non-smooth drop-shaped domain with a highly oscillatory forcing term; employing a discretization containing approximately four million degrees of freedom, the method produces the Poisson solution with accuracies of the order of machine precision in computing times on the order of one second on a single core of a present-day computer.