Multiscale S-Fraction Reduced-Order Models for Massive Wavefield Simulations

Abstract

We developed a novel reduced-order multiscale method for solving large time-domain wavefield simulation problems. Our algorithm consists of two main stages. During the first “off-line” stage the fine-grid operator (of the graph Laplacian type) is partitioned on coarse cells (subdomains). Then projection-type multiscale reduced-order models (ROMs) are computed for the coarse cell operators. The off-line stage is embarrassingly parallel as ROM computations for the subdomains are independent of each other. It also does not depend on the number of simulated sources (inputs) and it is performed just once before the entire time-domain simulation. At the second “on-line” stage the time-domain simulation is performed within the obtained multiscale ROM framework. The crucial feature of our formulation is the representation of the ROMs in terms of matrix Stieltjes continued fractions (S-fractions). The layered structure of the S-fraction introduces several hidden layers in the ROM representation that results in the block-tridiagonal dynamic system within each coarse cell. This allows us to sparsify the obtained multiscale subdomain operator ROMs and to reduce the communications between the adjacent subdomains, which is highly beneficial for a parallel implementation of the on-line stage. Our approach perfectly suits the high performance computing architectures; however, in this paper we present rather promising numerical results for a serial computing implementation only. These results include three-dimensional acoustic and multiphase anisotropic elastic problems.