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

DOI10.1088/1361-6420/AB750BzbMATH Open1478.65106arXiv1909.06460OpenAlexW3104215129WikidataQ118126684 ScholiaQ118126684MaRDI QIDQ5000606FDOQ5000606

Authors: L. Borcea, V. Druskin, Alexander V. Mamonov, Shari Moskow, M. Zaslavsky 

Publication date: 14 July 2021

Published in: Inverse Problems (Search for Journal in Brave)

Abstract: We generate data-driven reduced order models (ROMs) for inversion of the one and two dimensional Schr"odinger equation in the spectral domain given boundary data at a few frequencies. The ROM is the Galerkin projection of the Schr"odinger 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"odinger 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"odinger 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.

Full work available at URL: https://arxiv.org/abs/1909.06460

Recommendations

• Lippmann-Schwinger-Lanczos algorithm for inverse scattering problems
• Reduced order model approach to inverse scattering
• Reduced order models based on POD method for Schrödinger equations
• A reduced order model approach to inverse scattering in lossy layered media
• Direct, nonlinear inversion algorithm for hyperbolic problems via projection-based model reduction

zbMATH Keywords

inversionreduced order modelinternal dataSchrödinger

Mathematics Subject Classification ID

Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Boundary value problems for second-order elliptic equations (35J25) NLS equations (nonlinear Schrödinger equations) (35Q55) Finite difference methods for boundary value problems involving PDEs (65N06) Numerical methods for ill-posed problems for boundary value problems involving PDEs (65N20)

Cites Work

• Optical tomography: forward and inverse problems
• Electrical impedance tomography and Calderón's problem
• Rational Krylov approximation of matrix functions: numerical methods and optimal pole selection
• Solution of the time-domain inverse resistivity problem in the model reduction framework. I. One-dimensional problem with SISO data
• $\mathcal{H}_2$ Model Reduction for Large-Scale Linear Dynamical Systems
• Electrical impedance tomography
• Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain
• Three-point finite-difference schemes, Padé and the spectral Galerkin method. I. One-sided impedance approximation
• Circular resistor networks for electrical impedance tomography with partial boundary measurements
• Resistor network approaches to electrical impedance tomography
• Optimal finite difference grids for direct and inverse Sturm$ndash$Liouville problems
• Robust nonlinear processing of active array data in inverse scattering via truncated reduced order models
• A model reduction approach to numerical inversion for a parabolic partial differential equation
• Direct, Nonlinear Inversion Algorithm for Hyperbolic Problems via Projection-Based Model Reduction
• Untangling the nonlinearity in inverse scattering with data-driven reduced order models
• A discrete Liouville identity for numerical reconstruction of Schrödinger potentials
• A Nonlinear Method for Imaging with Acoustic Waves Via Reduced Order Model Backprojection
• Pyramidal resistor networks for electrical impedance tomography with partial boundary measurements

Cited In (8)

• On extension of the data driven ROM inverse scattering framework to partially nonreciprocal arrays
• Lippmann–Schwinger–Lanczos algorithm for inverse scattering problems
• A distributional Gelfand–Levitan–Marchenko equation for the Helmholtz scattering problem on the line
• Reduced order model approach for imaging with waves
• Waveform Inversion with a Data Driven Estimate of the Internal Wave
• On Learned Operator Correction in Inverse Problems
• Reduced order model based nonlinear waveform inversion for the 1D Helmholtz equation
• ROM inversion of monostatic data lifted to full MIMO

Uses Software

• Loewner

This page was built for publication: Reduced order models for spectral domain inversion: embedding into the continuous problem and generation of internal data

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5000606)