Incremental Proper Orthogonal Decomposition for PDE Simulation Data
Abstract
We propose an incremental algorithm to compute the proper orthogonal decomposition (POD) of simulation data for a partial differential equation. Specifically, we modify an incremental matrix SVD algorithm of Brand to accommodate data arising from Galerkin-type simulation methods for time dependent PDEs. The algorithm is applicable to data generated by many numerical methods for PDEs, including finite element and discontinuous Galerkin methods. The algorithm initializes and efficiently updates the dominant POD eigenvalues and modes during the time stepping in a PDE solver without storing the simulation data. We prove that the algorithm without truncation updates the POD exactly. We demonstrate the effectiveness of the algorithm using finite element computations for a 1D Burgers' equation and a 2D Navier-Stokes problem.
Recommended Citation
H. Fareed et al., "Incremental Proper Orthogonal Decomposition for PDE Simulation Data," Computers and Mathematics with Applications, vol. 75, no. 6, pp. 1942 - 1960, Elsevier, Mar 2018.
The definitive version is available at https://doi.org/10.1016/j.camwa.2017.09.012
Department(s)
Mathematics and Statistics
Research Center/Lab(s)
Center for High Performance Computing Research
Keywords and Phrases
Eigenvalues and eigenfunctions; Galerkin methods; Navier Stokes equations; Numerical methods; Partial differential equations; Principal component analysis; 2D Navier-Stokes problem; Burgers' equations; Discontinuous Galerkin methods; Finite element computations; Incremental algorithm; Proper orthogonal decompositions; Time dependent PDEs; Weighted norm; Finite element method
International Standard Serial Number (ISSN)
0898-1221
Document Type
Article - Journal
Document Version
Citation
File Type
text
Language(s)
English
Rights
© 2018 Elsevier, All rights reserved.
Publication Date
01 Mar 2018
Comments
J. Singler and Y. Zhang were supported in part by National Science Foundation grant DMS-1217122 . J. Singler and Y. Zhang thank the Institute for Mathematics and its Applications at the University of Minnesota for funding research visits, during which some of this work was completed. The authors thank the referees for their comments, which helped to improve the manuscript.