Source author record

Eric Michielssen

Eric Michielssen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

Wigner-Smith Time Delay Matrix for Acoustic Scattering: Computational Aspects

The Wigner-Smith (WS) time delay matrix relates an acoustic system's scattering matrix to its wavenumber derivative. The entries of the WS time delay matrix can be expressed in terms of energy density-like volume integrals, which cannot be efficiently evaluated in a boundary element method framework. This paper presents two schemes for efficiently populating the WS time delay matrix. The direct formulation casts the energy density-like volume integrals into integrals of the incident field and the field and/or its normal derivative over the scatterer surface. The indirect formulation computes the system's scattering matrix and its wavenumber derivative, again via surface integration, and then invokes the WS relationship to compute the WS time delay matrix. Both the direct and the indirect formulations yield equivalent results and can be easily integrated into standard boundary element codes.

preprint2021arXiv

The Design of Dual Band Stacked Metasurfaces Using Integral Equations

An integral equation-based approach for the design of dual band stacked metasurfaces which are invariant in one-dimension is presented. The stacked metasurface will generate collimated beams at desired angles in each band upon reflection. The conductor-backed stacked metasurface consists of two metasurfaces (a patterned metallic cladding supported by a dielectric spacer) stacked one upon the other. The stacked metasurface is designed in three phases. First the patterned metallic cladding of each metasurface is homogenized and modeled as an inhomogeneous impedance sheet. An Electric Field Integral Equation (EFIE) is written to model the mutual coupling between the homogenized elements within each metasurface, and from metasurface to metasurface. The EFIE is transformed into matrix equations by the method of moments. The nonlinear matrix equations are solved at both bands iteratively resulting in dual band complex-valued impedance sheets. In the second phase, optimization is applied to transform these complex-valued impedance sheets into purely reactive sheets suitable for printed circuit board fabrication by introducing surface waves. In the third phase, the metallic claddings of each metasurface are patterned for full-wave simulation of the dual band stacked metasurface. Using this approach, two dual band stacked metasurfaces are designed

preprint2021arXiv

Wigner-Smith Time Delay Matrix for Electromagnetics: Theory and Phenomenology

Wigner-Smith (WS) time delay concepts have been used extensively in quantum mechanics to characterize delays experienced by particles interacting with a potential well. This paper formally extends WS time delay theory to Maxwell's equations and explores its potential applications in electromagnetics. The WS time delay matrix relates a lossless and reciprocal system's scattering matrix to its frequency derivative and allows for the construction of modes that experience well-defined group delays when interacting with the system. The matrix' entries for guiding, scattering, and radiating systems are energy-like overlap integrals of the electric and/or magnetic fields that arise upon excitation of the system via its ports. The WS time delay matrix has numerous applications in electromagnetics, including the characterization of group delays in multiport systems, the description of electromagnetic fields in terms of elementary scattering processes, and the characterization of frequency sensitivities of fields and multiport antenna impedance matrices.

preprint2020arXiv

Butterfly factorization via randomized matrix-vector multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of a $m\times n$ matrix with $m\approx n$ provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of $O(\log n)$ sparse factors, each containing $O(n)$ nonzero entries. The factorization can be attained using $O(n^{3/2}\log n)$ computation and $O(n\log n)$ memory resources. The proposed algorithm applies to matrices with strong and weak admissibility conditions arising from surface integral equation solvers with a rigorous error bound, and is implemented in parallel.

preprint2019arXiv

Magnetic and Combined Field Integral Equations Based on the Quasi-Helmholtz Projectors

Boundary integral equation methods for analyzing electromagnetic scattering phenomena typically suffer from several of the following problems: (i) ill-conditioning when the frequency is low; (ii) ill-conditioning when the discretization density is high; (iii) ill-conditioning when the structure contains global loops (which are computationally expensive to detect); (iv) incorrect solution at low frequencies due to current cancellations; (v) presence of spurious resonances. In this paper, quasi-Helmholtz projectors are leveraged to obtain a magnetic field integral equation (MFIE) formulation that is immune to drawbacks (i)-(iv). Moreover, when this new MFIE is combined with a regularized electric field integral equation, a new quasi-Helmholtz projector combined field integral equation is obtained that also is immune to (v). Numerical results corroborate the theory and show the practical impact of the newly proposed formulations.

preprint2015arXiv

The transmission coefficient distribution of highly scattering sparse random media

We consider the distribution of the transmission coefficients, i.e. the singular values of the modal transmission matrix, for 2D random media with periodic boundary conditions composed of a large number of point-like non-absorbing scatterers. The scatterers are placed at random locations in the medium and have random refractive indices that are drawn from an arbitrary, known distribution. We construct a randomized model for the scattering matrix that retains scatterer dependent properties essential to reproduce the transmission coefficient distribution and analytically characterize the distribution of this matrix as a function of the refractive index distribution, the number of modes, and the number of scatterers. We show that the derived distribution agrees remarkably well with results obtained using a numerically rigorous spectrally accurate simulation. Analysis of the derived distribution provides the strongest principled justification yet of why we should expect perfect transmission in such random media regardless of the refractive index distribution of the constituent scatterers. The analysis suggests a sparsity condition under which random media will exhibit a perfect transmission-supporting universal transmission coefficient distribution in the deep medium limit.

preprint2014arXiv

Backscatter analysis based algorithms for increasing transmission through highly-scattering random media using phase-only modulated wavefronts

Recent theoretical and experimental advances have shed light on the existence of so-called `perfectly transmitting' wavefronts with transmission coefficients close to 1 in strongly backscattering random media. These perfectly transmitting eigen-wavefronts can be synthesized by spatial amplitude and phase modulation. Here, we consider the problem of transmission enhancement using phase-only modulated wavefronts. We develop physically realizable iterative and non-iterative algorithms for increasing the transmission through such random media using backscatter analysis. We theoretically show that, despite the phase-only modulation constraint, the non-iterative algorithms will achieve at least about 25$π$% or about 78.5% transmission assuming there is at least one perfectly transmitting eigen-wavefront and that the singular vectors of the transmission matrix obey a maximum entropy principle so that they are isotropically random. We numerically analyze the limits of phase-only modulated transmission in 2-D with fully spectrally accurate simulators and provide rigorous numerical evidence confirming our theoretical prediction in random media with periodic boundary conditions that is composed of hundreds of thousands of non-absorbing scatterers. We show via numerical simulations that the iterative algorithms we have developed converge rapidly, yielding highly transmitting wavefronts using relatively few measurements of the backscatter field. Specifically, the best performing iterative algorithm yields approx 70% transmission using just 15-20 measurements in the regime where the non-iterative algorithms yield approximately 78.5% transmission but require measuring the entire modal reflection matrix.

preprint2014arXiv

Iterative, backscatter-analysis algorithms for increasing transmission and focusing light through highly-scattering random media

Scattering hinders the passage of light through random media and consequently limits the usefulness of optical techniques for sensing and imaging. Thus, methods for increasing the transmission of light through such random media are of interest. Against this backdrop, recent theoretical and experimental advances have suggested the existence of a few highly transmitting eigen-wavefronts with transmission coefficients close to one in strongly backscattering random media. Here, we numerically analyze this phenomenon in 2-D with fully spectrally accurate simulators and provide rigorous numerical evidence confirming the existence of these highly transmitting eigen-wavefronts in random media with periodic boundary conditions that is composed of hundreds of thousands of non-absorbing scatterers. Motivated by bio-imaging applications where it is not possible to measure the transmitted fields, we develop physically realizable algorithms for increasing the transmission through such random media using backscatter analysis. We show via numerical simulations that the algorithms converge rapidly, yielding a near-optimum wavefront in just a few iterations. We also develop an algorithm that combines the knowledge of these highly transmitting eigen-wavefronts obtained from backscatter analysis, with intensity measurements at a point to produce a near-optimal focus with significantly fewer measurements than a method that does not utilize this information.