A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems
Journal Article
·
· Computer Methods in Applied Mechanics and Engineering
- NexGen Analytics, Sheridan, WY (United States)
- Sandia National Lab. (SNL-CA), Livermore, CA (United States)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.
- Research Organization:
- NexGen Analytics, Sheridan, WY (United States); Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC04-94AL85000; NA0003525
- OSTI ID:
- 1810771
- Alternate ID(s):
- OSTI ID: 1998972
- Report Number(s):
- SAND--2021-6660J; 696654
- Journal Information:
- Computer Methods in Applied Mechanics and Engineering, Journal Name: Computer Methods in Applied Mechanics and Engineering Vol. 384; ISSN 0045-7825
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
The Adjoint Petrov–Galerkin method for non-linear model reduction
Model reduction of dynamical systems on nonlinear manifolds using deep convolutional autoencoders
Model reduction for hypersonic aerodynamics via conservative LSPG projection and hyper-reduction
Journal Article
·
Wed Mar 18 20:00:00 EDT 2020
· Computer Methods in Applied Mechanics and Engineering
·
OSTI ID:1618094
Model reduction of dynamical systems on nonlinear manifolds using deep convolutional autoencoders
Journal Article
·
Thu Nov 14 19:00:00 EST 2019
· Journal of Computational Physics
·
OSTI ID:1574441
Model reduction for hypersonic aerodynamics via conservative LSPG projection and hyper-reduction
Conference
·
Sat Jan 04 23:00:00 EST 2020
· AIAA Scitech Forum
·
OSTI ID:1619242