Combining spectral and POD modes to improve error estimation of numerical model reduction for porous media

Abstract

Numerical model reduction (NMR) is used to solve the microscale problem that arises from computational homogenization of a model problem of porous media with displacement and pressure as unknown fields. The reduction technique and an associated error estimator for the NMR error have been presented in prior work, where both spectral decomposition (SD) and proper orthogonal decomposition (POD) were used to construct the reduced basis. It was shown that the POD basis performs better w.r.t. minimizing the residual, but the SD basis has some advantageous properties for the estimator. Since it is the estimated error that will govern the error control, the most efficient procedure is the one that results in the lowest error bound. The main contribution of this paper is further development of the previous work with a proposed combined basis constructed using both SD and POD modes together with an adaptive mode selection strategy. The performance of the combined basis is compared to (i) the pure SD basis and (ii) the pure POD basis via numerical examples. The examples show that it is possible to find a combination of SD/POD modes which is improved, i.e. it yields a smaller estimate, compared to the cases of pure SD or pure POD.