BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Vladimir Druskin

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

ResearcherUnclaimed source record
math.NAmath-phmath.MPNumerical Analysis

Catalog footprint

What is connected

14works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

14 published item(s)

preprint2022arXiv

Model order reduction of layered waveguides via rational Krylov fitting

Rational approximation recently emerged as an efficient numerical tool for the solution of exterior wave propagation problems. Currently, this technique is limited to wave media which are invariant along the main propagation direction. We propose a new model order reduction-based approach for compressing unbounded waveguides with layered inclusions. It is based on the solution of a nonlinear rational least squares problem using the RKFIT method. We show that approximants can be converted into an accurate finite difference representation within a rational Krylov framework. Numerical experiments indicate that RKFIT computes more accurate grids than previous analytic approaches and even works in the presence of pronounced scattering resonances. Spectral adaptation effects allow for finite difference grids with dimensions near or even below the Nyquist limit.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On extension of the data driven ROM inverse scattering framework to partially nonreciprocal arrays

Data-driven reduced order models (ROMs) recently emerged as powerful tool for the solution of inverse scattering problems. The main drawback of this approach is that it was limited to the measurement arrays with reciprocally collocated transmitters and receivers, that is, square symmetric matrix (data) transfer functions. To relax this limitation, we use our previous work [14], where the ROMs were combined with the Lippmann-Schwinger integral equation to produce a direct nonlinear inversion method. In this work we extend this approach to more general transfer functions, including those that are non-symmetric, e.g., obtained by adding only receivers or sources. The ROM is constructed based on the symmetric subset of the data and is used to construct all internal solutions. Remaining receivers are then used directly in the Lippmann-Schwinger equation. We demonstrate the new approach on a number of 1D and 2D examples with non-reciprocal arrays, including a single input/multiple outputs (SIMO) inverse problem, where the data is given by just a single-row matrix transfer function.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Direct, nonlinear inversion algorithm for hyperbolic problems via projection-based model reduction

We estimate the wave speed in the acoustic wave equation from boundary measurements by constructing a reduced-order model (ROM) matching discrete time-domain data. The state-variable representation of the ROM can be equivalently viewed as a Galerkin projection onto the Krylov subspace spanned by the snapshots of the time-domain solution. The success of our algorithm hinges on the data-driven Gram--Schmidt orthogonalization of the snapshots that suppresses multiple reflections and can be viewed as a discrete form of the Marchenko--Gel'fand--Levitan--Krein algorithm. In particular, the orthogonalized snapshots are localized functions, the (squared) norms of which are essentially weighted averages of the wave speed. The centers of mass of the squared orthogonalized snapshots provide us with the grid on which we reconstruct the velocity. This grid is weakly dependent on the wave speed in traveltime coordinates, so the grid points may be approximated by the centers of mass of the analogous set of squared orthogonalized snapshots generated by a known reference velocity. We present results of inversion experiments for one- and two-dimensional synthetic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exponential Convergence through Linear Finite Element Discretization of Stratified Subdomains

Motivated by problems where the response is needed at select localized regions in a large computational domain, we devise a novel finite element discretization that results in exponential convergence at pre-selected points. The two key features of the discretization are (a) use of midpoint integration to evaluate the contribution matrices, and (b) an unconventional bending of the mesh into complex space. Named complex-length finite element method (CFEM), the technique is linked to Pade approximants that provide exponential convergence of the Dirichlet-to-Neumann maps and thus the solution at specified points in the domain. Exponential convergence facilitates drastic reduction in the number of elements. This, combined with sparse computation associated with linear finite elements, results in significant reduction in the computational cost. The paper presents the basic ideas of the method as well as illustration of its effectiveness for a variety of problems involving Laplace, Helmholtz and elastodynamic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Near-optimal perfectly matched layers for indefinite Helmholtz problems

A new construction of an absorbing boundary condition for indefinite Helmholtz problems on unbounded domains is presented. This construction is based on a near-best uniform rational interpolant of the inverse square root function on the union of a negative and positive real interval, designed with the help of a classical result by Zolotarev. Using Krein's interpretation of a Stieltjes continued fraction, this interpolant can be converted into a three-term finite difference discretization of a perfectly matched layer (PML) which converges exponentially fast in the number of grid points. The convergence rate is asymptotically optimal for both propagative and evanescent wave modes. Several numerical experiments and illustrations are included.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonlinear seismic imaging via reduced order model backprojection

We introduce a novel nonlinear seismic imaging method based on model order reduction. The reduced order model (ROM) is an orthogonal projection of the wave equation propagator operator on the subspace of the snapshots of the solutions of the wave equation. It can be computed entirely from the knowledge of the measured time domain seismic data. The image is a backprojection of the ROM using the subspace basis for the known smooth kinematic velocity model. The implicit orthogonalization of solution snapshots is a nonlinear procedure that differentiates our approach from the conventional linear methods (Kirchhoff, RTM). It allows for the removal of multiple reflection artifacts. It also enables us to estimate the magnitude of the reflectors similarly to the true amplitude migration algorithms.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Extended Krylov Subspace Model-Order Reduction Technique to Simulate Wave Propagation in Unbounded Domains

In this paper we present a novel extended Krylov subspace reduced-order modeling technique to efficiently simulate time- and frequency-domain wavefields in open complex structures. To simulate the extension to infinity, we use an optimal complex-scaling method which is equivalent to an optimized perfectly matched layer in which the frequency is fixed. Wavefields propagating in strongly inhomogeneous open domains can now be modeled as a non-entire function of the complex-scaled wave operator. Since this function contains a square root singularity, we apply an extended Krylov subspace technique to construct fast converging reduced-order models. Specifically, we use a modified version of the extended Krylov subspace algorithm as proposed by Jagels and Reichel [Linear Algebra Appl., Vol. 434, pp. 1716 - 1732, 2011], since this algorithm allows us to balance the computational costs associated with computing powers of the wave operator and its inverse. Numerical experiments from electromagnetics and acoustics illustrate the performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

S-fraction multiscale finite-volume method for spectrally accurate wave propagation

We develop a method for numerical time-domain wave propagation based on the model order reduction approach. The method is built with high-performance computing (HPC) implementation in mind that implies a high level of parallelism and greatly reduced communication requirements compared to the traditional high-order finite-difference time-domain (FDTD) methods. The approach is inherently multiscale, with a reference fine grid model being split into subdomains. For each subdomain the coarse scale reduced order models (ROMs) are precomputed off-line in a parallel manner. The ROMs approximate the Neumann-to-Dirichlet (NtD) maps with high (spectral) accuracy and are used to couple the adjacent subdomains on the shared boundaries. The on-line part of the method is an explicit time stepping with the coupled ROMs. To lower the on-line computation cost the reduced order spatial operator is sparsified by transforming to a matrix Stieltjes continued fraction (S-fraction) form. The on-line communication costs are also reduced due to the ROM NtD map approximation properties. Another source of performance improvement is the time step length. Properly chosen ROMs substantially improve the Courant-Friedrichs-Lewy (CFL) condition. This allows the CFL time step to approach the Nyquist limit, which is typically unattainable with traditional schemes that have the CFL time step much smaller than the Nyquist sampling rate.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Krylov Stability-Corrected Coordinate-Stretching Method to Simulate Wave Propagation in Unbounded Domains

The Krylov subspace projection approach is a well-established tool for the reduced order modeling of dynamical systems in the time domain. In this paper, we address the main issues obstructing the application of this powerful approach to the time-domain solution of exterior wave problems. We use frequency independent perfectly matched layers to simulate the extension to infinity. Pure imaginary stretching functions based on Zolotarev's optimal rational approximation of the square root are implemented leading to perfectly matched layers with a controlled accuracy over a complete spectral interval of interest. A new Krylov-based solution method via stability-corrected operator exponents is presented which allows us to construct reduced-order models (ROMs) that respect the delicate spectral properties of the original scattering problem. The ROMs are unconditionally stable and are based on a renormalized bi-Lanczos algorithm. We give a theoretical foundation of our method and illustrate its performance through a number of numerical examples in which we simulate 2D electromagnetic wave propagation in unbounded domains, including a photonic waveguide example. The new algorithm outperforms the conventional finite-difference time domain method for problems on large time intervals.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access