Proper generalized decomposition solution of the parameterized Helmholtz problem: application to inverse geophysical problems

Abstract

This is the peer reviewed version of the following article: [Signorini, M., Zlotnik, S., and Díez, P. (2017) Proper generalized decomposition solution of the parameterized Helmholtz problem: application to inverse geophysical problems. Int. J. Numer. Meth. Engng, 109: 1085–1102. doi: 10.1002/nme.5313], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5313/full. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The identification of the geological structure from seismic data is formulated as an inverse problem. The properties and the shape of the rock formations in the subsoil are described by material and geometric parameters, which are taken as input data for a predictive model. Here, the model is based on the Helmholtz equation, describing the acoustic response of the system for a given wave length. Thus, the inverse problem consists in identifying the values of these parameters such that the output of the model agrees the best with observations. This optimization algorithm requires multiple queries to the model with different values of the parameters. Reduced order models are especially well suited to significantly reduce the computational overhead of the multiple evaluations of the model.\ud \ud In particular, the proper generalized decomposition produces a solution explicitly stating the parametric dependence, where the parameters play the same role as the physical coordinates. A proper generalized decomposition solver is devised to inexpensively explore the parametric space along the iterative process. This exploration of the parametric space is in fact seen as a post-process of the generalized solution. The approach adopted demonstrates its viability when tested in two illustrative examples.Peer ReviewedPostprint (author's final draft

Generalized parametric solutions in Stokes flow

“…The strategy to deal with geometric parameters in the PGD solver was devised in [12] for Poisson problems and then combined with material parameters in [13] and [14] for heat and wave propagation problems. The fundamental idea is using a reference domain T and a parametric mapping to the physical domain Ω(µ).…”

Certified real‐time shape optimization using isogeometric analysis, PGD model reduction, and a posteriori error estimation

“…An approximate solution of the model can then be obtained for any value of the parameters and be used in an online phase, with cheap and fast computations on light computing platforms, in order to perform real‐time parametric or sensitivity analysis for optimization, inverse identification, or optimal control purposes. The PGD technique was successfully implemented for many applications with parametric analyses, including parameterized geometry . Here, we go one step further by coupling PGD with IGA.…”

“…The PGD technique was successfully implemented for many applications with parametric analyses, including parameterized geometry. 8,35,[50][51][52] Here, we go one step further by coupling PGD with IGA. PGD extra coordinates are then related to the coordinates and weights of control points (at the coarsest-mesh level), and separation of variables is introduced between space and geometry variables.…”

Transport Map sampling with PGD model reduction for fast dynamical Bayesian data assimilation

“…23,24 The construction of the PGD solution is made in an offline phase and leads to a cost-effective evaluation of the numerical model depending on parameters in the online inversion phase. 25,26 The issue of the model evaluation cost is emphasized in Bayesian inference where the posterior density has to be explored to derive useful information on parameter pdfs such as mean, standard deviation, maximum, and marginals. Those quantities require the computation of multidimensional integrals over the parametric space.…”