Skip to main content

An improved algorithm for balanced POD through an analytic treatment of impulse response tails

Author(s): Tu, Jonathan H; Rowley, Clarence W

To refer to this page use:
Abstract: We present a modification of the balanced proper orthogonal decomposition (balanced POD) algorithm for systems with simple impulse response tails. In this new method, we use dynamic mode decomposition (DMD) to estimate the slowly decaying eigenvectors that dominate the long-time behavior of the direct and adjoint impulse responses. This is done using a new, low-memory variant of the DMD algorithm, appropriate for large datasets. We then formulate analytic expressions for the contribution of these eigenvectors to the controllability and observability Gramians. These contributions can be accounted for in the balanced POD algorithm by simply appending the impulse response snapshot matrices (direct and adjoint, respectively) with particular linear combinations of the slow eigenvectors. Aside from these additions to the snapshot matrices, the algorithm remains unchanged. By treating the tails analytically, we eliminate the need to run long impulse response simulations, lowering storage requirements and speeding up ensuing computations. To demonstrate its effectiveness, we apply this method to two examples: the linearized, complex Ginzburg–Landau equation, and the two-dimensional fluid flow past a cylinder. As expected, reduced-order models computed using an analytic tail match or exceed the accuracy of those computed using the standard balanced POD procedure, at a fraction of the cost.
Publication Date: 20-Jun-2012
Citation: Tu, Jonathan H, Rowley, Clarence W. "An improved algorithm for balanced POD through an analytic treatment of impulse response tails" Journal of Computational Physics,231, 16, 5317 - 5333, doi:10.1016/
DOI: doi:10.1016/
ISSN: 0021-9991
Pages: 5317 - 5333
Type of Material: Journal Article
Journal/Proceeding Title: Journal of Computational Physics
Version: This is the author’s final manuscript. All rights reserved to author(s).

Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.