Reduced order models for spectral domain inversion: embedding into the continuous problem and generation of internal data
From MaRDI portal
Publication:5000606
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)
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 . 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.
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
Cites work
- $\mathcal{H}_2$ Model Reduction for Large-Scale Linear Dynamical Systems
- A discrete Liouville identity for numerical reconstruction of Schrödinger potentials
- A model reduction approach to numerical inversion for a parabolic partial differential equation
- A nonlinear method for imaging with acoustic waves via reduced order model backprojection
- Circular resistor networks for electrical impedance tomography with partial boundary measurements
- Direct, nonlinear inversion algorithm for hyperbolic problems via projection-based model reduction
- Electrical impedance tomography
- Electrical impedance tomography and Calderón's problem
- Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain
- Optical tomography: forward and inverse problems
- Optimal finite difference grids for direct and inverse Sturm$ndash$Liouville problems
- Pyramidal resistor networks for electrical impedance tomography with partial boundary measurements
- Rational Krylov approximation of matrix functions: numerical methods and optimal pole selection
- Resistor network approaches to electrical impedance tomography
- Robust nonlinear processing of active array data in inverse scattering via truncated reduced order models
- Solution of the time-domain inverse resistivity problem in the model reduction framework. I. One-dimensional problem with SISO data
- Three-point finite-difference schemes, Padé and the spectral Galerkin method. I. One-sided impedance approximation
- Untangling the nonlinearity in inverse scattering with data-driven reduced order models
Cited in
(10)- On learned operator correction in inverse problems
- On extension of the data driven ROM inverse scattering framework to partially nonreciprocal arrays
- Waveform Inversion with a Data Driven Estimate of the Internal Wave
- Solving inverse scattering problems via reduced-order model embedding procedures
- Lippmann-Schwinger-Lanczos algorithm for inverse scattering problems
- Reduced order model approach for imaging with waves
- Reduced order model based nonlinear waveform inversion for the 1D Helmholtz equation
- Reduced order model approach to inverse scattering
- A distributional Gelfand–Levitan–Marchenko equation for the Helmholtz scattering problem on the line
- ROM inversion of monostatic data lifted to full MIMO
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)
