Quantized Tensor FEM for Multiscale Problems: Diffusion Problems in Two and Three Dimensions
Vladimir A. Kazeev, Ivan Valer'evich Oseledets, Maxim Rakhuba, Christoph Schwab · Multiscale Modeling and Simulation · 2022
Abstract. Homogenization in terms of multiscale limits transforms a multiscale problem with [Formula: see text] asymptotically separated microscales posed on a physical domain [Formula: see text] into a one-scale problem posed on a product domain of dimension [Formula: see text] by introducing [Formula: see text] so-called fast variables. This procedure allows one to convert [Formula: see text] scales in [Formula: see text] physical dimensions into a single-scale structure in [Formula: see text] dimensions. We prove here that both the original, physical multiscale problem and the corresponding high-dimensional, one-scale limiting problem can be efficiently treated numerically with the recently developed quantized tensor-train finite-element method (QTT-FEM). The QTT-FE approximation consists in restricting approximation and computation to sequences of nested subspaces, each of which is a tensor product of two factors of low dimension (rank), within a vast, but generic “virtual” (background) discretization space. In practice, these subspaces are determined iteratively and data-adaptively at runtime bypassing any “offline precomputation.” For theoretical analysis, low-dimensional subspaces are constructed analytically to bound the tensor ranks against the error tolerance. We consider a model linear elliptic multiscale problem in several physical dimensions and show, theoretically and experimentally, that both (i) the solution of the associated high-dimensional one-scale problem and (ii) the approximation to the solution of the multiscale problem induced thereby admit efficient QTT-FE approximations. These problems can therefore be numerically solved in a scale-robust fashion by standard (low-order) PDE discretizations combined with state-of-the-art general-purpose solvers for tensor-structured linear systems. Specifically, we prove the existence of QTT-FE approximations with an upper bound on the tensor ranks growing no faster than algebraically with respect to [Formula: see text] and independent of the scale parameters, where [Formula: see text] is the target accuracy for the approximation of the solution and of its gradient. In numerical experiments, we verify the theoretical rank bounds and computationally investigate the dependence of the complexity of the solutions on the number [Formula: see text] of microscales.