Component-based reduced basis for parametrized symmetric eigenproblems
Name
40323_2015_Article_21.pdf
Size
3.27 MB
Format
Adobe PDF
Checksum (MD5)
601e390039089b02ebeef4658b3dbab4
Author(s)
Vallaghe, Sylvain P.
Huynh, Dinh Bao Phuong
Knezevic, David
Nguyen, Loi
Patera, Anthony T.
Date Issued
May 2015
Journal
Advanced Modeling and Simulation in Engineering Sciences
Publisher
Springer International Publishing
Citation
Vallaghe, Sylvain, Phuong Huynh, David J. Knezevic, Loi Nguyen, and Anthony T. Patera. “Component-Based Reduced Basis for Parametrized Symmetric Eigenproblems.” Advanced Modeling and Simulation in Engineering Sciences 2, no. 1 (May 23, 2015).
Version
Final published version
Abstract
Background:
A component-based approach is introduced for fast and flexible solution of parameter-dependent symmetric eigenproblems.
Methods:
Considering a generalized eigenproblem with symmetric stiffness and mass operators, we start by introducing a “ σ-shifted” eigenproblem where the left hand side operator corresponds to an equilibrium between the stiffness operator and a weighted mass operator, with weight-parameter σ>0. Assuming that σ=λ n >0, the nth real positive eigenvalue of the original eigenproblem, then the shifted eigenproblem reduces to the solution of a homogeneous linear problem. In this context, we can apply the static condensation reduced basis element (SCRBE) method, a domain synthesis approach with reduced basis (RB) approximation at the intradomain level to populate a Schur complement at the interdomain level. In the Offline stage, for a library of archetype subdomains we train RB spaces for a family of linear problems; these linear problems correspond to various equilibriums between the stiffness operator and the weighted mass operator. In the Online stage we assemble instantiated subdomains and perform static condensation to obtain the “ σ-shifted” eigenproblem for the full system. We then perform a direct search to find the values of σ that yield singular systems, corresponding to the eigenvalues of the original eigenproblem.
Results:
We provide eigenvalue a posteriori error estimators and we present various numerical results to demonstrate the accuracy, flexibility and computational efficiency of our approach.
Conclusions:
We are able to obtain large speed and memory improvements compared to a classical Finite Element Method (FEM), making our method very suitable for large models commonly considered in an engineering context.
MIT Department
Massachusetts Institute of Technology. Department of Mechanical Engineering
Terms of Use
Creative Commons Attribution
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1186/s40323-015-0021-0
