TY - JOUR
T1 - Stability of discrete empirical interpolation and gappy proper orthogonal decomposition with randomized and deterministic sampling points
AU - Peherstorfer, Benjamin
AU - Drmac, Zlatko
AU - Gugercin, Serkan
N1 - Funding Information:
\ast Submitted to the journal's Methods and Algorithms for Scientific Computing section December 17, 2019; accepted for publication (in revised form) July 13, 2020; published electronically September 21, 2020. https://doi.org/10.1137/19M1307391 \bfF \bfu \bfn \bfd \bfi \bfn \bfg : The work of the first author was partially supported by the Air Force Center of Excellence on Multi-Fidelity Modeling of Rocket Combustor Dynamics through grant FA9550-17-1-0195. The work of the second author was partially supported by the Croatian Science Foundation through grant IP-2019-04-6268. The work of the third author was partially supported by the NSF through grants DMS-1522616 and DMS-1819110.
Funding Information:
The authors thank Karthik Duraisamy (University of Michigan), Cheng Huang (University of Michigan), and David Xu (University of Michigan) for providing the snapshots corresponding to the single-injector combustion process discussed in subsection 6.2.
Publisher Copyright:
© 2020 Benjamin Peherstorfer, Serkan Gugercin, Zlatko Drmac
PY - 2020
Y1 - 2020
N2 - This work investigates the stability of (discrete) empirical interpolation for nonlinear model reduction and state feld approximation from measurements. Empirical interpolation derives approximations from a few samples (measurements) via interpolation in low-dimensional spaces. It has been observed that empirical interpolation can become unstable if the samples are perturbed due to, e.g., noise, turbulence, and numerical inaccuracies. The main contribution of this work is a probabilistic analysis that shows that stable approximations are obtained if samples are randomized and if more samples than dimensions of the low-dimensional spaces are used. Oversampling, i.e., taking more sampling points than dimensions of the low-dimensional spaces, leads to approximations via regression and is known under the name of gappy proper orthogonal decomposition. Building on the insights of the probabilistic analysis, a deterministic sampling strategy is presented that aims to achieve lower approximation errors with fewer points than randomized sampling by taking information about the low-dimensional spaces into account. Numerical results of reconstructing velocity felds from noisy measurements of combustion processes and model reduction in the presence of noise demonstrate the instability of empirical interpolation and the stability of gappy proper orthogonal decomposition with oversampling.
AB - This work investigates the stability of (discrete) empirical interpolation for nonlinear model reduction and state feld approximation from measurements. Empirical interpolation derives approximations from a few samples (measurements) via interpolation in low-dimensional spaces. It has been observed that empirical interpolation can become unstable if the samples are perturbed due to, e.g., noise, turbulence, and numerical inaccuracies. The main contribution of this work is a probabilistic analysis that shows that stable approximations are obtained if samples are randomized and if more samples than dimensions of the low-dimensional spaces are used. Oversampling, i.e., taking more sampling points than dimensions of the low-dimensional spaces, leads to approximations via regression and is known under the name of gappy proper orthogonal decomposition. Building on the insights of the probabilistic analysis, a deterministic sampling strategy is presented that aims to achieve lower approximation errors with fewer points than randomized sampling by taking information about the low-dimensional spaces into account. Numerical results of reconstructing velocity felds from noisy measurements of combustion processes and model reduction in the presence of noise demonstrate the instability of empirical interpolation and the stability of gappy proper orthogonal decomposition with oversampling.
KW - Empirical interpolation
KW - Gappy proper orthogonal decomposition
KW - Model reduction
KW - Noisy observations
KW - Nonlinear model reduction
KW - Randomized model reduction
UR - http://www.scopus.com/inward/record.url?scp=85093107106&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85093107106&partnerID=8YFLogxK
U2 - 10.1137/19M1307391
DO - 10.1137/19M1307391
M3 - Article
AN - SCOPUS:85093107106
SN - 1064-8275
VL - 42
SP - A2837-A2864
JO - SIAM Journal of Scientific Computing
JF - SIAM Journal of Scientific Computing
IS - 5
M1 - 19M1307391
ER -