Boundary element tearing and interconnecting methods in unbounded domains
Applied Numerical Mathematics
Multiscale finite element methods for high-contrast problems using local spectral basis functions
Journal of Computational Physics
Quasi-Monte Carlo methods for elliptic PDEs with random coefficients and applications
Journal of Computational Physics
Robust solvers for symmetric positive definite operators and weighted poincaré inequalities
LSSC'11 Proceedings of the 8th international conference on Large-Scale Scientific Computing
Local-global multiscale model reduction for flows in high-contrast heterogeneous media
Journal of Computational Physics
Auxiliary variables for 3D multiscale simulations in heterogeneous porous media
Journal of Computational Physics
Robust multilevel solvers for high-contrast anisotropic multiscale problems
Journal of Computational and Applied Mathematics
Ensemble level multiscale finite element and preconditioner for channelized systems and applications
Journal of Computational and Applied Mathematics
Hi-index | 0.02 |
We consider additive Schwarz domain decomposition preconditioners for piecewise linear finite element approximations of elliptic PDEs with highly variable coefficients. In contrast to standard analyses, we do not assume that the coefficients can be resolved by a coarse mesh. This situation arises often in practice, for example in the computation of flows in heterogeneous porous media, in both the deterministic and (Monte–Carlo simulated) stochastic cases. We consider preconditioners which combine local solves on general overlapping subdomains together with a global solve on a general coarse space of functions on a coarse grid. We perform a new analysis of the preconditioned matrix, which shows rather explicitly how its condition number depends on the variable coefficient in the PDE as well as on the coarse mesh and overlap parameters. The classical estimates for this preconditioner with linear coarsening guarantee good conditioning only when the coefficient varies mildly inside the coarse grid elements. By contrast, our new results show that, with a good choice of subdomains and coarse space basis functions, the preconditioner can still be robust even for large coefficient variation inside domains, when the classical method fails to be robust. In particular our estimates prove very precisely the previously made empirical observation that the use of low-energy coarse spaces can lead to robust preconditioners. We go on to consider coarse spaces constructed from multiscale finite elements and prove that preconditioners using this type of coarsening lead to robust preconditioners for a variety of binary (i.e., two-scale) media model problems. Moreover numerical experiments show that the new preconditioner has greatly improved performance over standard preconditioners even in the random coefficient case. We show also how the analysis extends in a straightforward way to multiplicative versions of the Schwarz method.