Optimal Local Approximation Spaces for Component-Based Static Condensation Procedures
Name
Optimal local.pdf
Size
852.09 KB
Format
Adobe PDF
Checksum (MD5)
0c37ae214e0cf44c8bbb8d028e81cd00
Author(s)
Smetana, Kathrin
Patera, Anthony T
Date Issued
October 2016
Journal
SIAM Journal on Scientific Computing
Publisher
Society for Industrial and Applied Mathematics
Citation
Smetana, Kathrin, and Anthony T. Patera. “Optimal Local Approximation Spaces for Component-Based Static Condensation Procedures.” SIAM Journal on Scientific Computing 38, no. 5 (January 2016): A3318–A3356. © SIAM
Version
Final published version
Abstract
In this paper we introduce local approximation spaces for component-based static condensation (sc) procedures that are optimal in the sense of Kolmogorov. To facilitate simulations for large structures such as aircraft or ships, it is crucial to decrease the number of degrees of freedom on the interfaces, or “ports,” in order to reduce the size of the statically condensed system. To derive optimal port spaces we consider a (compact) transfer operator that acts on the space of harmonic extensions on a two-component system and maps the traces on the ports that lie on the boundary of these components to the trace of the shared port. Solving the eigenproblem for the composition of the transfer operator and its adjoint yields the optimal space. For a related work in the context of the generalized finite element method, we refer the reader to [I. Babuška and R. Lipton, Multiscale Model. Simul., 9 (2011), pp. 373--406]. We further introduce a spectral greedy algorithm to generalize the procedure to the parameter-dependent setting and to construct a quasi-optimal parameter-independent port space. Moreover, it is shown that, given a certain tolerance and an upper bound for the ports in the system, the spectral greedy constructs a port space that yields an sc approximation error on a system of arbitrary configuration which is smaller than this tolerance for all parameters in a rich train set. We present our approach for isotropic linear elasticity, although the idea may be readily applied to any linear coercive problem. Numerical experiments demonstrate the very rapid and exponential convergence both of the eigenvalues and of the sc approximation based on spectral modes for nonseparable and irregular geometries such as an I-beam with an internal crack.
MIT Department
Massachusetts Institute of Technology. Department of Mechanical Engineering
Terms of Use
Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1137/15M1009603
