Carlos Pérez-Arancibia

High-order integral equation methods for problems of scattering by bumps and cavities on half-planes

This paper presents high-order integral equation methods for the evaluation of electromagnetic wave scattering by dielectric bumps and dielectric cavities on perfectly conducting or dielectric half-planes. In detail, the algorithms introduced in this paper apply to eight classical scattering problems, namely, scattering by a dielectric bump on a perfectly conducting or a dielectric half-plane, and scattering by a filled, overfilled, or void dielectric cavity on a perfectly conducting or a dielectric half-plane. In all cases field representations based on single-layer potentials for appropriately chosen Green functions are used. The numerical far fields and near fields exhibit excellent convergence as discretizations are refined-even at and around points where singular fields and infinite currents exist.

On the Green's function for the Helmholtz operator in an impedance circular cylindrical waveguide

This paper addresses the problem of finding a series representation for the Greens function of the Helmholtz operator in an infinite circular cylindrical waveguide with impedance boundary condition. Resorting to the Fourier transform, complex analysis techniques and the limiting absorption principle (when the undamped case is analyzed), a detailed deduction of the Greens function is performed, generalizing the results available in the literature for the case of a complex impedance parameter. Procedures to obtain numerical values of the Greens function are also developed in this article.

Electromagnetic power absorption due to bumps and trenches on flat surfaces

This paper presents a study of the absorption of electromagnetic power that results from the interaction of electromagnetic waves and cylindrical bumps or trenches on flat conducting surfaces. Configurations are characterized by means of adequately selected dimensionless variables and parameters so that applicability to mathematically equivalent (but physically diverse) systems can be achieved easily. Electromagnetic fields and absorption increments caused by such surface defects are evaluated by means of a high-order integral equation method which resolves fine details of the field near the surface, and which was validated by fully analytical approaches in a range of computationally challenging cases. The computational method is also applied to problems concerning bumps and trenches on imperfect conducting planes for which analytical solutions are not available. Typically, we find that absorption is enhanced by the presence of the defects considered, although, interestingly, absorption can also be significantly reduced in some cases—such as, e.g., in the case of a trench on a conducting plane where the incident electric field is perpendicular to the plane. Additionally, it is observed that, for some small-skin-depths large-wavelengths, the absorption increment is proportional to the increase in surface area. Significant physical insight is obtained on the heating that results from various types of electromagnetic incident fields.

Windowed Green function method for the Helmholtz equation in the presence of multiply layered media

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 the presence of an arbitrary number of penetrable layers. Relying on the use of certain slow-rise windowing functions, the proposed windowed Green function approach 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 multilayer structures. The proposed methodology, whose theoretical basis was presented in the recent contribution (Bruno et al. 2016 SIAM J. Appl. Math. 76, 1871–1898. (doi:10.1137/15M1033782)), 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 the use of Sommerfeld integrals.

Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D

Abstract We present an effective harmonic density interpolation method for the numerical evaluation of singular and nearly singular Laplace boundary integral operators and layer potentials in two and three spatial dimensions. The method relies on the use of Greens third identity and local Taylor-like interpolations of density functions in terms of harmonic polynomials. The proposed technique effectively regularizes the singularities present in boundary integral operators and layer potentials, and recasts the latter in terms of integrands that are bounded or even more regular, depending on the order of the density interpolation. The resulting boundary integrals can then be easily, accurately, and inexpensively evaluated by means of standard quadrature rules. A variety of numerical examples demonstrate the effectiveness of the technique when used in conjunction with the classical trapezoidal rule (to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type quadrature rule (to integrate over surfaces given as unions of non-overlapping quadrilateral patches) in three-dimensions.

Multitrace/singletrace formulations and Domain Decomposition Methods for the solution of Helmholtz transmission problems for bounded composite scatterers

Abstract We present Nystrom discretizations of multitrace/singletrace formulations and non-overlapping Domain Decomposition Methods (DDM) for the solution of Helmholtz transmission problems for bounded composite scatterers with piecewise constant material properties. We investigate the performance of DDM with both classical Robin and optimized transmission boundary conditions. The optimized transmission boundary conditions incorporate square root Fourier multiplier approximations of Dirichlet to Neumann operators. While the multitrace/singletrace formulations as well as the DDM that use classical Robin transmission conditions are not particularly well suited for Krylov subspace iterative solutions of high-contrast high-frequency Helmholtz transmission problems, we provide ample numerical evidence that DDM with optimized transmission conditions constitute efficient computational alternatives for these type of applications. In the case of large numbers of subdomains with different material properties, we show that the associated DDM linear system can be efficiently solved via hierarchical Schur complements elimination.

Windowed Green Function Method for Nonuniform Open-Waveguide Problems

This contribution presents a novel Windowed Green Function (WGF) method for the solution of problems of wave propagation, scattering, and radiation for structures that include open (dielectric) waveguides, waveguide junctions, as well as launching and/or termination sites and other nonuniformities. Based on the use of a “slow-rise” smooth-windowing technique in conjunction with free-space Green functions and associated integral representations, the proposed approach produces numerical solutions with errors that decrease faster than any negative power of the window size. The proposed methodology bypasses some of the most significant challenges associated with waveguide simulation. In particular, the WGF approach handles spatially infinite dielectric waveguide structures without recourse to absorbing boundary conditions, it facilitates proper treatment of complex geometries, and it seamlessly incorporates the open-waveguide character and associated radiation conditions inherent in the problem under consideration. The overall WGF approach is demonstrated in this paper by means of a variety of numerical results for 2-D open-waveguide termination, launching and junction problems.