Title

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.

Department(s)

Mathematics and Statistics

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.

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.

Share

 
COinS