Source author record

Fernando Guevara Vasquez

Fernando Guevara Vasquez 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

17works
10topics
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

17 published item(s)

preprint2022arXiv

Active exterior cloaking for the 2D Helmholtz equation with complex wavenumbers and application to thermal cloaking

We design sources for the two-dimensional Helmholtz equation that can cloak an object by cancelling out the incident field in a region, without the sources completely surrounding the object to hide. As in previous work for real positive wavenumbers, the sources are also determined by the Green identities. The novelty is that we prove that the same approach works for complex wavenumbers which makes it applicable to a variety of media, including media with dispersion, loss and gain. Furthermore, by deriving bounds on Graf's addition formulas with complex arguments, we obtain new estimates that allow to quantify the quality of the cloaking effect. We illustrate our results by applying them to achieve active exterior cloaking for the heat equation.

preprint2019arXiv

Broadband lightweight flat lenses for longwave-infrared imaging

We experimentally demonstrate imaging in the longwave-infrared (LWIR) spectral band (8um to 12um) using a single polymer flat lens based upon multi-level diffractive optics. The device thickness is only 10μm, and chromatic aberrations are corrected over the entire LWIR band with one surface. Due to the drastic reduction in device thickness, we are able to utilize polymers with absorption in the LWIR, allowing for inexpensive manufacturing via imprint lithography. The weight of our lens is less than 100 times those of comparable refractive lenses. We fabricated and characterized two different flat lenses. Even with about 25% absorption losses, experiments show that our flat polymer lenses obtain good imaging with field of view of ~35degrees and angular resolution less than 0.013 degrees. The flat lenses were characterized with two different commercial LWIR image sensors. Finally, we show that by using lossless, higher-refractive-index materials like silicon, focusing efficiencies in excess of 70% can be achieved over the entire LWIR band. Our results firmly establish the potential for lightweight, ultra-thin, broadband lenses for high-quality imaging in the LWIR band.

preprint2019arXiv

Periodic particle arrangements using standing acoustic waves

We determine crystal-like materials that can be fabricated by using a standing acoustic wave to arrange small particles in a non-viscous liquid resin, which is cured afterwards to keep the particles in the desired locations. For identical spherical particles with the same physical properties and small compared to the wavelength, the locations where the particles are trapped correspond to the minima of an acoustic radiation potential which describes the net forces that a particle is subject to. We show that the global minima of spatially periodic acoustic radiation potentials can be predicted by the eigenspace of a small real symmetric matrix corresponding to its smallest eigenvalue. We relate symmetries of this eigenspace to particle arrangements composed of points, lines or planes. Since waves are used to generate the particle arrangements, the arrangement's periodicity is limited to certain Bravais lattice classes that we enumerate in two and three dimensions.

preprint2016arXiv

Kirchhoff migration without phases

We present a simple, frequency domain, preprocessing step to Kirchhoff migration that allows the method to image scatterers when the wave field phase information is lost at the receivers, and only intensities are measured. The resulting imaging method does not require knowing the phases of the probing field or manipulating the phase of the wave field at the receivers. In a regime where the scattered field is small compared to the probing field, the problem of recovering the full-waveform scattered field from intensity data can be formulated as an embarrassingly simple least-squares problem. Although this only recovers the projection (on a known subspace) of the full-waveform scattered field, we show that, for high frequencies, this projection gives Kirchhoff images asymptotically identical to the images obtained with full waveform data. Our method can also be used when the source is modulated by a Gaussian process and autocorrelations are measured at an array of receivers.

preprint2015arXiv

Imaging with power controlled source pairs

Scatterers in a homogeneous medium are imaged by probing the medium with two point sources of waves modulated by correlated signals and by measuring only intensities at one single receiver. For appropriately chosen source pairs, we show that full waveform array measurements can be recovered from such intensity measurements by solving a linear least squares problem. The least squares solution can be used to image with Kirchhoff migration, even if the solution is determined only up to a known one-dimensional nullspace. The same imaging strategy can be used when the medium is probed with point sources driven by correlated Gaussian processes and autocorrelations are measured at a single location. Since autocorrelations are robust to noise, this can be used for imaging when the probing wave is drowned in background noise.

preprint2015arXiv

On the solvability of the discrete conductivity and Schrödinger inverse problems

We study the uniqueness question for two inverse problems on graphs. Both problems consist in finding (possibly complex) edge or nodal based quantities from boundary measurements of solutions to the Dirichlet problem associated with a weighted graph Laplacian plus a diagonal perturbation. The weights can be thought of as a discrete conductivity and the diagonal perturbation as a discrete Schrödinger potential. We use a discrete analogue to the complex geometric optics approach to show that if the linearized problem is solvable about some conductivity (or Schrödinger potential) then the linearized problem is solvable for almost all conductivities (or Schrödinger potentials) in a suitable set. We show that the conductivities (or Schrödinger potentials) in a certain set are determined uniquely by boundary data, except on a zero measure set. This criterion for solvability is used in a statistical study of graphs where the conductivity or Schrödinger inverse problem is solvable.

preprint2013arXiv

Characterization and Synthesis of Rayleigh Damped Elastodynamic Networks

We consider damped elastodynamic networks where the damping matrix is assumed to be a non-negative linear combination of the stiffness and mass matrices (also known as Rayleigh or proportional damping). We give here a characterization of the frequency response of such networks. We also answer the synthesis question for such networks, i.e., how to construct a Rayleigh damped elastodynamic network with a given frequency response. Our analysis shows that not all damped elastodynamic networks can be realized when the proportionality constants between the damping matrix and the mass and stiffness matrices are fixed.

preprint2013arXiv

Restarted inverse Born series for the Schrödinger problem with discrete internal measurements

Convergence and stability results for the inverse Born series [Moskow and Schotland, Inverse Problems, 24:065005, 2008] are generalized to mappings between Banach spaces. We show that by restarting the inverse Born series one obtains a class of iterative methods containing the Gauss-Newton and Chebyshev-Halley methods. We use the generalized inverse Born series results to show convergence of the inverse Born series for the Schrödinger problem with discrete internal measurements. In this problem, the Schrödinger potential is to be recovered from a few measurements of solutions to the Schrödinger equation resulting from a few different source terms. An application of this method to a problem related to transient hydraulic tomography is given, where the source terms model injection and measurement wells.

preprint2012arXiv

Pseudospectra of Isospectrally Reduced Matrices and Systems

The isospectral reduction of matrix, which is closely related to its Schur complement, allows to reduce the size of a matrix while maintaining its eigenvalues up to a known set. Here we generalize this procedure by increasing the number of possible ways a matrix can be isospectrally reduced. The reduced matrix has rational functions as entries. We show that the notion of pseudospectrum can be extended to this class of matrices and that the pseudospectrum of a matrix shrinks as the matrix is reduced. Hence the eigenvalues of a reduced matrix are more robust to entry-wise perturbations than the eigenvalues of the original matrix. We also introduce the notion of inverse pseudospectrum (or pseudoresonances), which indicates how stable the poles of a matrix with rational function entries are to certain matrix perturbations. A mass spring system is used to illustrate and give a physical interpretation to both pseudospectra and inverse pseudospectra.

preprint2011arXiv

Mathematical analysis of the two dimensional active exterior cloaking in the quasistatic regime

We design a device that generates fields canceling out a known probing field inside a region to be cloaked while generating very small fields far away from the device. The fields we consider satisfy the Laplace equation, but the approach remains valid in the quasistatic regime in a homogeneous medium. We start by relating the problem of designing an exterior cloak in the quasistatic regime to the classic problem of approximating a harmonic function with harmonic polynomials. An explicit polynomial solution to the problem was given earlier in [Phys. Rev. Lett. 103 (2009), 073901]. Here we show convergence of the device field to the field needed to perfectly cloak an object. The convergence region limits the size of the cloaked region, and the size and position of the device.

preprint2011arXiv

Resistor network approaches to electrical impedance tomography

We review a resistor network approach to the numerical solution of the inverse problem of electrical impedance tomography (EIT). The networks arise in the context of finite volume discretizations of the elliptic equation for the electric potential, on sparse and adaptively refined grids that we call optimal. The name refers to the fact that the grids give spectrally accurate approximations of the Dirichlet to Neumann map, the data in EIT. The fundamental feature of the optimal grids in inversion is that they connect the discrete inverse problem for resistor networks to the continuum EIT problem.

preprint2011arXiv

Study of noise effects in electrical impedance tomography with resistor networks

We present a study of the numerical solution of the two dimensional electrical impedance tomography problem, with noisy measurements of the Dirichlet to Neumann map. The inversion uses parametrizations of the conductivity on optimal grids. The grids are optimal in the sense that finite volume discretizations on them give spectrally accurate approximations of the Dirichlet to Neumann map. The approximations are Dirichlet to Neumann maps of special resistor networks, that are uniquely recoverable from the measurements. Inversion on optimal grids has been proposed and analyzed recently, but the study of noise effects on the inversion has not been carried out. In this paper we present a numerical study of both the linearized and the nonlinear inverse problem. We take three different parametrizations of the unknown conductivity, with the same number of degrees of freedom. We obtain that the parametrization induced by the inversion on optimal grids is the most efficient of the three, because it gives the smallest standard deviation of the maximum a posteriori estimates of the conductivity, uniformly in the domain. For the nonlinear problem we compute the mean and variance of the maximum a posteriori estimates of the conductivity, on optimal grids. For small noise, we obtain that the estimates are unbiased and their variance is very close to the optimal one, given by the Cramer-Rao bound. For larger noise we use regularization and quantify the trade-off between reducing the variance and introducing bias in the solution. Both the full and partial measurement setups are considered.

preprint2011arXiv

Transformation elastodynamics and active exterior acoustic cloaking

This chapter consists of three parts. In the first part we recall the elastodynamic equations under coordinate transformations. The idea is to use coordinate transformations to manipulate waves propagating in an elastic material. Then we study the effect of transformations on a mass-spring network model. The transformed networks can be realized with "torque springs", which are introduced here and are springs with a force proportional to the displacement in a direction other than the direction of the spring terminals. Possible homogenizations of the transformed networks are presented, with potential applications to cloaking. In the second and third parts we present cloaking methods that are based on cancelling an incident field using active devices which are exterior to the cloaked region and that do not generate significant fields far away from the devices. In the second part, the exterior cloaking problem for the Laplace equation is reformulated as the problem of polynomial approximation of analytic functions. An explicit solution is given that allows to cloak larger objects at a fixed distance from the cloaking device, compared to previous explicit solutions. In the third part we consider the active exterior cloaking problem for the Helmholtz equation in 3D. Our method uses the Green's formula and an addition theorem for spherical outgoing waves to design devices that mimic the effect of the single and double layer potentials in Green's formula.

preprint2010arXiv

Exterior cloaking with active sources in two dimensional acoustics

We cloak a region from a known incident wave by surrounding the region with three or more devices that cancel out the field in the cloaked region without significantly radiating waves. Since very little waves reach scatterers within the cloaked region, the scattered field is small and the scatterers are for all practical purposes undetectable. The devices are multipolar point sources that can be determined from Green's formula and an addition theorem for Hankel functions. The cloaking devices are exterior to the cloaked region.

preprint2009arXiv

Active Exterior Cloaking

A new method of cloaking is presented. For two-dimensional quasistatics it is proven how a single active exterior cloaking device can be used to shield an object from surrounding fields, yet produce very small scattered fields. The problem is reduced to finding a polynomial which is approximately one within one disk and zero within a second disk, and such a polynomial is constructed. For the two-dimensional Helmholtz equation, it is numerically shown that three active exterior devices placed around the object suffice to produce very good cloaking.

preprint2009arXiv

Broadband Exterior Cloaking

It is shown how a recently proposed method of cloaking is effective over a broad range of frequencies. The method is based on three or more active devices. The devices, while not radiating significantly, create a ``quiet zone'' between the devices where the wave amplitude is small. Objects placed within this region are virtually invisible. The cloaking is demonstrated by simulations with a broadband incident pulse.

preprint2009arXiv

Complete characterization and synthesis of the response function of elastodynamic networks

The response function of a network of springs and masses, an elastodynamic network, is the matrix valued function $W(ω)$, depending on the frequency $ω$, mapping the displacements of some accessible or terminal nodes to the net forces at the terminals. We give necessary and sufficient conditions for a given function $W(ω)$ to be the response function of an elastodynamic network, assuming there is no damping. In particular we construct an elastodynamic network that can mimic a suitable response in the frequency or time domain. Our characterization is valid for networks in three dimensions and also for planar networks, which are networks where all the elements, displacements and forces are in a plane. The network we design can fit within an arbitrarily small neighborhood of the convex hull of the terminal nodes, provided the springs and masses occupy an arbitrarily small volume. Additionally, we prove stability of the network response to small changes in the spring constants and/or addition of springs with small spring constants.