A model reduction approach to numerical inversion for a parabolic partial differential equation
From MaRDI portal
Publication:2941005
DOI10.1088/0266-5611/30/12/125011zbMATH Open1323.35207arXiv1210.1257OpenAlexW2141409082WikidataQ115293867 ScholiaQ115293867MaRDI QIDQ2941005FDOQ2941005
Authors: L. Borcea, V. Druskin, Alexander V. Mamonov, M. Zaslavsky 
Publication date: 21 January 2015
Published in: Inverse Problems (Search for Journal in Brave)
Abstract: 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.
Full work available at URL: https://arxiv.org/abs/1210.1257
Recommendations
- Solution of the time-domain inverse resistivity problem in the model reduction framework. I. One-dimensional problem with SISO data
- On combining model reduction and Gauss–Newton algorithms for inverse partial differential equation problems
- Multigrid Algorithms for Inverse Problems with Linear Parabolic PDE Constraints
- A fast algorithm for parabolic PDE-based inverse problems based on Laplace transforms and flexible Krylov solvers
- scientific article; zbMATH DE number 665281
rational interpolationmodel order reductionparabolic problemnumerical parameter identificationrational Krylov subspace projection
Cited In (22)
- On extension of the data driven ROM inverse scattering framework to partially nonreciprocal arrays
- Solution of the time-domain inverse resistivity problem in the model reduction framework. I. One-dimensional problem with SISO data
- Untangling the nonlinearity in inverse scattering with data-driven reduced order models
- Transport Reversal for Model Reduction of Hyperbolic Partial Differential Equations
- When data driven reduced order modeling meets full waveform inversion
- Reduced order models for spectral domain inversion: embedding into the continuous problem and generation of internal data
- A fast algorithm for parabolic PDE-based inverse problems based on Laplace transforms and flexible Krylov solvers
- Electromagnetic inverse wave scattering in anisotropic media via reduced order modeling
- A new algorithm to determine the creation or depletion term of parabolic equations from boundary measurements
- A numerical method for an inverse source problem for parabolic equations and its application to a coefficient inverse problem
- A discrete Liouville identity for numerical reconstruction of Schrödinger potentials
- Lippmann-Schwinger-Lanczos algorithm for inverse scattering problems
- Direct, nonlinear inversion algorithm for hyperbolic problems via projection-based model reduction
- A reduced order model approach to inverse scattering in lossy layered media
- Nonlinear parametric inversion using interpolatory model reduction
- Robust nonlinear processing of active array data in inverse scattering via truncated reduced order models
- Distance preserving model order reduction of graph-Laplacians and cluster analysis
- ROM inversion of monostatic data lifted to full MIMO
- A convergent numerical method to recover the initial condition of nonlinear parabolic equations from lateral Cauchy data
- On combining model reduction and Gauss–Newton algorithms for inverse partial differential equation problems
- Computing reduced order models via inner-outer Krylov recycling in diffuse optical tomography
- A nonlinear method for imaging with acoustic waves via reduced order model backprojection
This page was built for publication: A model reduction approach to numerical inversion for a parabolic partial differential equation
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2941005)
