Source author record

Liliana Borcea

Liliana Borcea 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

26works
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

26 published item(s)

preprint2022arXiv

Paraxial wave propagation in random media with long-range correlations

We study the paraxial wave equation with a randomly perturbed index of refraction, which can model the propagation of a wave beam in a turbulent medium. The random perturbation is a stationary and isotropic process with a general form of the covariance that may be integrable or not. We focus attention mostly on the non-integrable case, which corresponds to a random perturbation with long-range correlations, that is relevant for propagation through a cloudy turbulent atmosphere. The analysis is carried out in a high-frequency regime where the forward scattering approximation holds. It reveals that the randomization of the wave field is multiscale: The travel time of the wave front is randomized at short distances of propagation and it can be described by a fractional Brownian motion. The wave field observed in the random travel time frame is affected by the random perturbations at long distances, and it is described by a Schroedinger-type equation driven by a standard Brownian field. We use these results to quantify how scattering leads to decorrelation of the spatial and spectral components of the wave field and to a deformation of the pulse emitted by the source. These are important questions for applications like imaging and free space communications with pulsed laser beams through a turbulent atmosphere. We also compare the results with those used in the optics literature, which are based on the Kolmogorov model of turbulence.

preprint2022arXiv

Velocity estimation via model order reduction

A novel approach to full waveform inversion (FWI), based on a data driven reduced order model (ROM) of the wave equation operator is introduced. The unknown medium is probed with pulses and the time domain pressure waveform data is recorded on an active array of sensors. The ROM, a projection of the wave equation operator is constructed from the data via a nonlinear process and is used for efficient velocity estimation. While the conventional FWI via nonlinear least-squares data fitting is challenging without low frequency information, and prone to getting stuck in local minima (cycle skipping), minimization of ROM misfit is behaved much better, even for a poor initial guess. For low-dimensional parametrizations of the unknown velocity the ROM misfit function is close to convex. The proposed approach consistently outperforms conventional FWI in standard synthetic tests.

preprint2021arXiv

Reduced order model approach for imaging with waves

We introduce a novel, computationally inexpensive approach for imaging with an active array of sensors, which probe an unknown medium with a pulse and measure the resulting waves. The imaging function uses a data driven estimate of the "internal wave" originating from the vicinity of the imaging point and propagating to the sensors through the unknown medium. We explain how this estimate can be obtained using a reduced order model (ROM) for the wave propagation. We analyze the imaging function, connect it to the time reversal process and describe how its resolution depends on the aperture of the array, the bandwidth of the probing pulse and the medium through which the waves propagate. We also show how the internal wave can be used for selective focusing of waves at points in the imaging region. This can be implemented experimentally and can be used for pixel scanning imaging. We assess the performance of the imaging methods with numerical simulations and compare them to the conventional reverse-time migration method and the "backprojection" method introduced recently as an application of the same ROM.

preprint2020arXiv

Reduced Order Model Approach to Inverse Scattering

We study an inverse scattering problem for a generic hyperbolic system of equations with an unknown coefficient called the reflectivity. The solution of the system models waves (sound, electromagnetic or elastic), and the reflectivity models unknown scatterers embedded in a smooth and known medium. The inverse problem is to determine the reflectivity from the time resolved scattering matrix (the data) measured by an array of sensors. We introduce a novel inversion method, based on a reduced order model (ROM) of an operator called wave propagator, because it maps the wave from one time instant to the next, at interval corresponding to the discrete time sampling of the data. The wave propagator is unknown in the inverse problem, but the ROM can be computed directly from the data. By construction, the ROM inherits key properties of the wave propagator, which facilitate the estimation of the reflectivity. The ROM was introduced previously and was used for two purposes: (1) to map the scattering matrix to that corresponding to the single scattering (Born) approximation and (2) to image i.e., obtain a qualitative estimate of the support of the reflectivity. Here we study further the ROM and show that it corresponds to a Galerkin projection of the wave propagator. The Galerkin framework is useful for proving properties of the ROM that are used in the new inversion method which seeks a quantitative estimate of the reflectivity.

preprint2019arXiv

Factorization method versus migration imaging in a waveguide

We present a comparative study of two qualitative imaging methods in an acoustic waveguide with sound hard walls. The waveguide terminates at one end and contains unknown obstacles of compact support, to be determined from data gathered by an array of sensors that probe the obstacles with waves and measure the scattered response. The first imaging method, known as the factorization method, is based on the factorization of the far field operator. It is designed to image at single frequency and estimates the support of the obstacles by a Picard range criterion. The second imaging method, known as migration, works either with one or multiple frequencies. It forms an image by backpropagating the measured scattered wave to the search points, using the Green's function in the empty waveguide. We study the connection between these methods with analysis and numerical simulations.

preprint2019arXiv

Reduced order models for spectral domain inversion: Embedding into the continuous problem and generation of internal data

We generate data-driven reduced order models (ROMs) for inversion of the one and two dimensional Schrödinger equation in the spectral domain given boundary data at a few frequencies. The ROM is the Galerkin projection of the Schrödinger operator onto the space spanned by solutions at these sample frequencies. The ROM matrix is in general full, and not good for extracting the potential. However, using an orthogonal change of basis via Lanczos iteration, we can transform the ROM to a block triadiagonal form from which it is easier to extract $q$. In one dimension, the tridiagonal matrix corresponds to a three-point staggered finite-difference system for the Schrödinger operator discretized on a so-called spectrally matched grid which is almost independent of the medium. In higher dimensions, the orthogonalized basis functions play the role of the grid steps. The orthogonalized basis functions are localized and also depend only very weakly on the medium, and thus by embedding into the continuous problem, the reduced order model yields highly accurate internal solutions. That is to say, we can obtain, just from boundary data, very good approximations of the solution of the Schrödinger equation in the whole domain for a spectral interval that includes the sample frequencies. We present inversion experiments based on the internal solutions in one and two dimensions.

preprint2016arXiv

Derivation of a one-way radiative transfer equation in random media

We derive from first principles a one-way radiative transfer equation for the wave intensity resolved over directions (Wigner transform of the wave field) in random media. It is an initial value problem with excitation from a source which emits waves in a preferred, forward direction. The equation is derived in a regime with small random fluctuations of the wave speed but long distances of propagation with respect to the wavelength, so that cumulative scattering is significant. The correlation length of the medium and the scale of the support of the source are slightly larger than the wavelength, and the waves propagate in a wide cone with opening angle less than $180^o$, so that the backward and evanescent waves are negligible. The scattering regime is a bridge between that of radiative transfer, where the waves propagate in all directions and the paraxial regime, where the waves propagate in a narrow angular cone. We connect the one-way radiative transport equation with the equations satisfied by the Wigner transform of the wave field in these regimes.

preprint2016arXiv

Imaging in random media with convex optimization

We study an inverse problem for the wave equation where localized wave sources in random scattering media are to be determined from time resolved measurements of the waves at an array of receivers. The sources are far from the array, so the measurements are affected by cumulative scattering in the medium, but they are not further than a transport mean free path, which is the length scale characteristic of the onset of wave diffusion that prohibits coherent imaging. The inversion is based on the Coherent Interferometric (CINT) imaging method which mitigates the scattering effects by introducing an appropriate smoothing operation in the image formation. This smoothing stabilizes statistically the images, at the expense of their resolution. We complement the CINT method with a convex ($l_1$) optimization in order to improve the source localization and obtain quantitative estimates of the source intensities. We analyze the method in a regime where scattering can be modeled by large random wavefront distortions, and quantify the accuracy of the inversion in terms of the spatial separation of individual sources or clusters of sources. The theoretical predictions are demonstrated with numerical simulations.

preprint2016arXiv

Polarization effects for electromagnetic wave propagation in random media

We study Maxwell's equations in random media with small fluctuations of the electric permittivity. We consider a setup where the waves propagate toward a preferred direction, called range. We decompose the electromagnetic wave field in transverse electric and transverse magnetic plane waves, called modes, with random amplitudes that model cumulative scattering effects in the medium. Their evolution in range is described by a coupled system of stochastic differential equations driven by the random fluctuations of the electric permittivity. We analyze the solution of this system with the Markov limit theorem and obtain a detailed asymptotic characterization of the electromagnetic wave field in the long range limit. In particular, we quantify the loss of coherence of the waves due to scattering by calculating the range scales (scattering mean free paths) on which the mean amplitudes of the modes decay. We also quantify the energy exchange between the modes, and consequently the loss of polarization induced by scattering, by analyzing the Wigner transform (energy density) of the electromagnetic wave field. This analysis involves the derivation of transport equations with polarization. We study in detail these equations and connect the results with the existing literature in radiative transport and paraxial wave propagation.

preprint2016arXiv

Pulse propagation in time dependent randomly layered media

We study cumulative scattering effects on wave front propagation in time dependent randomly layered media. It is well known that the wave front has a deterministic characterization in time independent media, aside from a small random shift in the travel time. That is, the pulse shape is predictable, but faded and smeared as described mathematically by a convolution kernel determined by the autocorrelation of the random fluctuations of the wave speed. The main result of this paper is the extension of the pulse stabilization results to time dependent randomly layered media. When the media change slowly, on time scales that are longer than the pulse width and the time it takes the waves to traverse a correlation length, the pulse is not affected by the time fluctuations. In rapidly changing media, where these time scales are similar, both the pulse shape and the random component of the arrival time are affected by the statistics of the time fluctuations of the wave speed. We obtain an integral equation for the wave front, that is more complicated than in time independent media, and cannot be solved analytically, in general. We also give examples of media where the equation simplifies, and the wave front can be analyzed explicitly. We illustrate with these examples how the time fluctuations feed energy into the pulse.

preprint2016arXiv

Pulse reflection in a random waveguide with a turning point

We present an analysis of wave propagation and reflection in an acoustic waveguide with random sound soft boundary and a turning point. The waveguide has slowly bending axis and variable cross section. The variation consists of a slow and monotone change of the width of the waveguide and small and rapid fluctuations of the boundary, on the scale of the wavelength. These fluctuations are modeled as random. The turning point is many wavelengths away from the source, which emits a pulse that propagates toward the turning point, where it is reflected. To focus attention on this reflection, we assume that the waveguide supports a single propagating mode from the source to the turning point, beyond which all the waves are evanescent. We consider a scaling regime where scattering at the random boundary has a significant effect on the reflected pulse. In this regime scattering from the random boundary away from the turning point is negligible, while scattering from the random boundary around the turning point results in a strong, deterministic pulse deformation. The reflected pulse shape is not the same as the emitted one. It is damped, due to scattering at the boundary, and is deformed by dispersion in the waveguide. The reflected pulse also carries a random phase.

preprint2016arXiv

Robust imaging with electromagnetic waves in noisy environments

We study imaging with an array of sensors that probes a medium with single frequency electromagnetic waves and records the scattered electric field. The medium is known and homogenous except for some small and penetrable inclusions. The goal of inversion is to locate and characterize these inclusions from the data collected by the array, which are corrupted by additive noise. We use results from random matrix theory to obtain a robust inversion method. We assess its performance with numerical simulations and quantify the benefit of measuring more than one component of the scattered electric field.

preprint2016arXiv

Synthetic Aperture Imaging of Direction and Frequency Dependent Reflectivities

We introduce a synthetic aperture imaging framework that takes into consideration directional dependence of the reflectivity that is to be imaged, as well as its frequency dependence. We use an $\ell_1$ minimization approach that is coordinated with data segmentation so as to fuse information from multiple sub-apertures and frequency sub-bands. We analyze this approach from first principles and assess its performance with numerical simulations in an X-band radar regime.

preprint2015arXiv

Imaging with electromagnetic waves in terminating waveguides

We study an inverse scattering problem for Maxwell's equations in terminating waveguides, where localized reflectors are to be imaged using a remote array of sensors. The array probes the waveguide with waves and measures the scattered returns. The mathematical formulation of the inverse scattering problem is based on the electromagnetic Lippmann-Schwinger integral equation and an explicit calculation of the Green tensor. The image formation is carried with reverse time migration and with $\ell_1$ optimization.

preprint2015arXiv

Resolution analysis of imaging with $\ell_1$ optimization

We study array imaging of a sparse scene of point-like sources or scatterers in a homogeneous medium. For source imaging the sensors in the array are receivers that collect measurements of the wave field. For imaging scatterers the array probes the medium with waves and records the echoes. In either case the image formation is stated as a sparsity promoting $\ell_1$ optimization problem, and the goal of the paper is to quantify the resolution. We consider both narrow-band and broad-band imaging, and a geometric setup with a small array. We take first the case of the unknowns lying on the imaging grid, and derive resolution limits that depend on the sparsity of the scene. Then we consider the general case with the unknowns at arbitrary locations. The analysis is based on estimates of the cumulative mutual coherence and a related concept, which we call interaction coefficient. It complements recent results in compressed sensing by deriving deterministic resolution limits that account for worse case scenarios in terms of locations of the unknowns in the imaging region, and also by interpreting the results in some cases where uniqueness of the solution does not hold. We demonstrate the theoretical predictions with numerical simulations.

preprint2015arXiv

Source estimation with incoherent waves in random waveguides

We study an inverse source problem for the acoustic wave equation in a random waveguide. The goal is to estimate the source of waves from measurements of the acoustic pressure at a remote array of sensors. The waveguide effect is due to boundaries that trap the waves and guide them in a preferred (range) direction, the waveguide axis, along which the medium is unbounded. The random waveguide is a model of perturbed ideal waveguides which have flat boundaries and are filled with known media that do not change with range. The perturbation consists of fluctuations of the boundary and of the wave speed due to numerous small inhomogeneities in the medium. The fluctuations are uncertain in applications, which is why we model them with random processes, and they cause significant cumulative scattering at long ranges from the source. The scattering effect manifests mathematically as an exponential decay of the expectation of the acoustic pressure, the coherent part of the wave. The incoherent wave is modeled by the random fluctuations of the acoustic pressure, which dominate the expectation at long ranges from the source. We use the existing theory of wave propagation in random waveguides to analyze the inverse problem of estimating the source from incoherent wave recordings at remote arrays. We show how to obtain from the incoherent measurements high fidelity estimates of the time resolved energy carried by the waveguide modes, and study the invertibility of the system of transport equations that model energy propagation in order to estimate the source.

preprint2014arXiv

A model reduction approach to numerical inversion for a parabolic partial differential equation

We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss-Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments.

preprint2013arXiv

A quantitative study of source imaging in random waveguides

We present a quantitative study of coherent array imaging of remote sources in randomly perturbed waveguides with bounded cross-section. We study how long range cumulative scattering by perturbations of the boundary and the medium impedes the imaging process. We show that boundary scattering effects can be mitigated with filters that enhance the coherent part of the data. The filters are obtained by optimizing a measure of quality of the image. The point is that there is an optimal trade-off between the robustness and resolution of images in such waveguides, which can be found adaptively, as the data are processed to form the image. Long range scattering by perturbations of the medium is harder to mitigate than scattering by randomly perturbed boundaries. Coherent imaging methods do not work and more complex incoherent methods, based on transport models of energy, should be used instead. Such methods are nor useful, nor needed in waveguides with perturbed boundaries. We explain all these facts using rigorous asymptotic stochastic analysis of the wave field in randomly perturbed waveguides. We also analyze the adaptive coherent imaging method and obtain a quantitative agreement with the results of numerical simulations.

preprint2013arXiv

Electromagnetic wave propagation in random waveguides

We study long range propagation of electromagnetic waves in random waveguides with rectangular cross-section and perfectly conducting boundaries. The waveguide is filled with an isotropic linear dielectric material, with randomly fluctuating electric permittivity. The fluctuations are weak, but they cause significant cumulative scattering over long distances of propagation of the waves. We decompose the wave field in propagating and evanescent transverse electric and magnetic modes with random amplitudes that encode the cumulative scattering effects. They satisfy a coupled system of stochastic differential equations driven by the random fluctuations of the electric permittivity. We analyze the solution of this system with the diffusion approximation theorem, under the assumption that the fluctuations decorrelate rapidly in the range direction. The result is a detailed characterization of the transport of energy in the waveguide, the loss of coherence of the modes and the depolarization of the waves due to cumulative scattering.

preprint2012arXiv

Motion Estimation and Imaging of Complex Scenes with Synthetic Aperture Radar

We study synthetic aperture radar (SAR) imaging and motion estimation of complex scenes consisting of stationary and moving targets. We use the classic SAR setup with a single antenna emitting signals and receiving the echoes from the scene. The known motion estimation methods for such setups work only in simple cases, with one or a few targets in the same motion. We propose to extend the applicability of these methods to complex scenes, by complementing them with a data pre-processing step intended to separate the echoes from the stationary targets and the moving ones. We present two approaches. The first is an iteration designed to subtract the echoes from the stationary targets one by one. It estimates the location of each stationary target from a preliminary image, and then uses it to define a filter that removes its echo from the data. The second approach is based on the robust principle component analysis (PCA) method. The key observation is that with appropriate pre-processing and windowing, the discrete samples of the stationary target echoes form a low rank matrix, whereas the samples of a few moving target echoes form a high rank sparse matrix. The robust PCA method is designed to separate the low rank from the sparse part, and thus can be used for the SAR data separation. We present a brief analysis of the two methods and explain how they can be combined to improve the data separation for extended and complex imaging scenes. We also assess the performance of the methods with extensive numerical simulations.

preprint2012arXiv

Paraxial coupling of propagating modes in three-dimensional waveguides with random boundaries

We analyze long range wave propagation in three-dimensional random waveguides. The waves are trapped by top and bottom boundaries, but the medium is unbounded in the two remaining directions. We consider scalar waves, and motivated by applications in underwater acoustics, we take a pressure release boundary condition at the top surface and a rigid bottom boundary. The wave speed in the waveguide is known and smooth, but the top boundary has small random fluctuations that cause significant cumulative scattering of the waves over long distances of propagation. To quantify the scattering effects, we study the evolution of the random amplitudes of the waveguide modes. We obtain that in the long range limit they satisfy a system of paraxial equations driven by a Brownian field. We use this system to estimate three important mode-dependent scales: the scattering mean free path, the cross-range decoherence length and the decoherence frequency. Understanding these scales is important in imaging and communication problems, because they encode the cumulative scattering effects in the wave field measured by remote sensors. As an application of the theory, we analyze time reversal and coherent interferometric imaging in strong cumulative scattering regimes.

preprint2012arXiv

Synthetic Aperture Radar Imaging and Motion Estimation via Robust Principle Component Analysis

We consider the problem of synthetic aperture radar (SAR) imaging and motion estimation of complex scenes. By complex we mean scenes with multiple targets, stationary and in motion. We use the usual setup with one moving antenna emitting and receiving signals. We address two challenges: (1) the detection of moving targets in the complex scene and (2) the separation of the echoes from the stationary targets and those from the moving targets. Such separation allows high resolution imaging of the stationary scene and motion estimation with the echoes from the moving targets alone. We show that the robust principal component analysis (PCA) method which decomposes a matrix in two parts, one low rank and one sparse, can be used for motion detection and data separation. The matrix that is decomposed is the pulse and range compressed SAR data indexed by two discrete time variables: the slow time, which parametrizes the location of the antenna, and the fast time, which parametrizes the echoes received between successive emissions from the antenna. We present an analysis of the rank of the data matrix to motivate the use of the robust PCA method. We also show with numerical simulations that successful data separation with robust PCA requires proper data windowing. Results of motion estimation and imaging with the separated data are presented, as well.

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

Wave propagation in waveguides with random boundaries

We give a detailed analysis of long range cumulative scattering effects from rough boundaries in waveguides. We assume small random fluctuations of the boundaries and obtain a quantitative statistical description of the wave field. The method of solution is based on coordinate changes that straighten the boundaries. The resulting problem is similar from the mathematical point of view to that of wave propagation in random waveguides with interior inhomogeneities. We quantify the net effect of scattering at the random boundaries and show how it differs from that of scattering by internal inhomogeneities.