×

zbMATH — the first resource for mathematics

An adaptive scheme for a class of interpolatory model reduction methods for frequency response problems. (English) Zbl 1352.74132
Summary: Frequency sweep problems arise in many structural dynamic, acoustic, and structural acoustic applications. In each case, they incur the evaluation of a frequency response function for a typically large number of frequencies. Because each function evaluation requires the solution of an often large-scale system of equations, frequency sweep problems are computationally intensive. Interpolatory model order reduction is a powerful tool for reducing their cost. However, the performance of this tool depends on the location and number of the interpolation frequency points. It also depends on the number of consecutive frequency derivatives of the response function that are matched at each frequency point. So far, these two choices have been made in the literature in a heuristic manner. In contrast, this paper proposes an automatic adaptive strategy based on monitoring the Euclidean norm of the relative residual associated with the function to be evaluated over the frequency band of interest. More specifically, the number of interpolation points and the number of matched frequency derivatives are adaptively increased until the global Euclidean norm of the relative residual is reduced below a user-specified tolerance. The robustness, accuracy, and computational efficiency of this adaptive strategy are highlighted with the solution of several frequency sweep problems associated with large-scale structural dynamic, acoustic, and structural acoustic finite element models.

MSC:
74H15 Numerical approximation of solutions of dynamical problems in solid mechanics
65D05 Numerical interpolation
74H45 Vibrations in dynamical problems in solid mechanics
76Q05 Hydro- and aero-acoustics
Software:
MUMPS
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Antoulas, Approximation of Large-Scale Dynamical Systems (2005) · Zbl 1112.93002 · doi:10.1137/1.9780898718713
[2] Bai, Numerical Methods in Electromagnetics XIII, in: Handbook of Numerical Analysis pp 825– (2005)
[3] Freund, Model reduction methods based on Krylov subspaces, Acta Numerica 12 pp 267– (2003) · Zbl 1046.65021 · doi:10.1017/S0962492902000120
[4] Hetmaniuk, Review and assessment of interpolatory model order reduction methods for frequency response structural dynamics and acoustics problems, International Journal for Numerical Methods in Engineering 90 pp 1636– (2012) · Zbl 1246.74023 · doi:10.1002/nme.4271
[5] Carlberg, Efficient nonlinear model reduction via a least-squares Petrov-Galerkin projection and compressive tensor approximations, International Journal for Numerical Methods in Engineering 86 pp 155– (2011) · Zbl 1235.74351 · doi:10.1002/nme.3050
[6] Avery, Fast frequency sweep computations using a multi-point Padé-based reconstruction method and an efficient iterative solver, International Journal for Numerical Methods in Engineering 69 (13) pp 2848– (2007) · Zbl 1194.74107 · doi:10.1002/nme.1879
[7] Beattie, Interpolatory projection methods for structure-preserving model reduction, Systems & Control Letters 58 pp 225– (2009) · Zbl 1159.93317 · doi:10.1016/j.sysconle.2008.10.016
[8] Druskin, Solution of large scale evolutionary problems using rational Krylov subspaces with optimized shifts, SIAM Journal on Scientific Computing 31 (5) pp 3760– (2009) · Zbl 1204.65042 · doi:10.1137/080742403
[9] Gugercin, H2 model reduction for large-scale linear dynamical systems, SIAM Journal on Matrix Analysis and Applications 30 (2) pp 609– (2008) · Zbl 1159.93318 · doi:10.1137/060666123
[10] Tuck-Lee, Adaptive frequency windowing for multifrequency solutions in structural acoustics based on the matrix Padé-via-Lanczos algorithm, International Journal for Numerical Methods in Engineering 73 (5) pp 728– (2008) · Zbl 1159.74048 · doi:10.1002/nme.2102
[11] Olsson, Rational Krylov for eigenvalue computation and model order reduction, BIT Numerical Mathematics 46 pp S99– (2006) · Zbl 1105.65036 · doi:10.1007/s10543-006-0085-9
[12] Huynh, High-fidelity real-time simulation on deployed platforms, Computer and Fluids 43 pp 74– (2011) · Zbl 05909599 · doi:10.1016/j.compfluid.2010.07.007
[13] Farhat C Hetmaniuk U Tezaur R A simple adaptive scheme for a class of interpolatory model reduction methods for frequency response problems Proceedings of the 10th International Conference on the Mathematical and Numerical Aspects of Waves 2011 327 330 · Zbl 1246.74023
[14] Djellouli, A fast method for solving acoustic scattering problems in frequency bands, Journal of Computational Physics 168 pp 412– (2001) · Zbl 1153.76394 · doi:10.1006/jcph.2001.6707
[15] Avery, A Padé-based factorization-free algorithm for identifying the eigenvalues missed by a generalized symmetric eigensolver, International Journal for Numerical Methods in Engineering 79 pp 239– (2009) · Zbl 1171.74462 · doi:10.1002/nme.2572
[16] Djellouli, Finite element solution of two-dimensional acoustic scattering problems using arbitrarily shaped convex artificial boundaries, Journal of Computational Acoustics 8 pp 81– (2000) · Zbl 1360.76131 · doi:10.1142/S0218396X00000066
[17] Tezaur, Three-dimensional finite element calculations in acoustic scattering using arbitrarily shaped convex artificial boundaries, International Journal for Numerical Methods in Engineering 53 pp 1461– (2002) · Zbl 0996.76058 · doi:10.1002/nme.346
[18] Berenger, A perfectly matched layer for the absorption of electromagnetic waves, Journal of Computational Physics 114 (2) pp 185– (1994) · Zbl 0814.65129 · doi:10.1006/jcph.1994.1159
[19] Amestoy, A fully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM Journal on Matrix Analysis and Applications 23 pp 15– (2001) · Zbl 0992.65018 · doi:10.1137/S0895479899358194
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.