View in EDS

COMPRESSING LARGE-SCALE WAVE PROPAGATION MODELS VIA PHASE-PRECONDITIONED RATIONAL KRYLOV SUBSPACES.

Saved in:
Bibliographic Details
Title: COMPRESSING LARGE-SCALE WAVE PROPAGATION MODELS VIA PHASE-PRECONDITIONED RATIONAL KRYLOV SUBSPACES.
Authors: DRUSKIN, VLADIMIR1 druskin@slb.com, REMIS, ROB F.2 r.f.remis@tudelft.nl, ZASLAVSKY, MIKHAIL1 mzaslavsky@slb.com, ZIMMERLING, JÖRN T.2 jzimmerl@umich.edu
Source: Multiscale Modeling & Simulation. 2018, Vol. 16 Issue 4, p1486-1518. 33p.
Subjects: Theory of wave motion, High performance computing, Green's functions, Krylov subspace, Acoustic models, Rogue waves
Abstract: Rational Krylov subspace (RKS) techniques are well-established and powerful tools for projection-based model reduction of time-invariant dynamic systems. For hyperbolic wavefield problems, such techniques perform well in configurations where only a few modes contribute to the field. RKS methods, however, are fundamentally limited by the Nyquist--Shannon sampling rate, making them unsuitable for the approximation of wavefields in configuration characterized by large travel times and propagation distances, since wavefield responses in such configurations are highly oscillatory in the frequency-domain. To overcome this limitation, we propose to precondition the RKSs by factoring out the rapidly varying frequency-domain field oscillations. The remaining amplitude-functions are generally slowly varying functions of source position and spatial coordinate and allow for a significant compression of the approximation subspace. Our one-dimensional analysis together with numerical experiments for large-scale two-dimensional acoustic models shows superior approximation properties of preconditioned RKS compared with the standard RKS model-order reduction. The preconditioned RKS results in a reduction of the frequency sampling well below the Nyquist--Shannon rate, a weak dependence of the RKS size on the number of inputs and outputs for multiple-input/multiple-output problems, and, most importantly, in a significant coarsening of the finite-difference grid used to generate the RKS. A prototype implementation indicates that the preconditioned RKS algorithm is competitive in the modern high performance computing environment. [ABSTRACT FROM AUTHOR]
Copyright of Multiscale Modeling & Simulation is the property of Society for Industrial & Applied Mathematics and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Engineering Source
Description
Abstract:Rational Krylov subspace (RKS) techniques are well-established and powerful tools for projection-based model reduction of time-invariant dynamic systems. For hyperbolic wavefield problems, such techniques perform well in configurations where only a few modes contribute to the field. RKS methods, however, are fundamentally limited by the Nyquist--Shannon sampling rate, making them unsuitable for the approximation of wavefields in configuration characterized by large travel times and propagation distances, since wavefield responses in such configurations are highly oscillatory in the frequency-domain. To overcome this limitation, we propose to precondition the RKSs by factoring out the rapidly varying frequency-domain field oscillations. The remaining amplitude-functions are generally slowly varying functions of source position and spatial coordinate and allow for a significant compression of the approximation subspace. Our one-dimensional analysis together with numerical experiments for large-scale two-dimensional acoustic models shows superior approximation properties of preconditioned RKS compared with the standard RKS model-order reduction. The preconditioned RKS results in a reduction of the frequency sampling well below the Nyquist--Shannon rate, a weak dependence of the RKS size on the number of inputs and outputs for multiple-input/multiple-output problems, and, most importantly, in a significant coarsening of the finite-difference grid used to generate the RKS. A prototype implementation indicates that the preconditioned RKS algorithm is competitive in the modern high performance computing environment. [ABSTRACT FROM AUTHOR]
ISSN:15403459
DOI:10.1137/17M1156848