Full text loading...
-
Scalable Algebraic Multiscale Linear Solver for Large-scale Reservoir Simulation
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, ECMOR XIV - 14th European Conference on the Mathematics of Oil Recovery, Sep 2014, Volume 2014, p.1 - 20
Abstract
The scalability of the Algebraic Multiscale Solver (AMS) ( Wang et al., 2014 ) for the heterogeneous pressure system that arises from incompressible flow in porous media is analyzed and experimentally demonstrated in parallel computing environments. The major steps of AMS are highly parallel, but the solver overall scalability is strongly tied to the choice of parameters and algorithms involved in each step. These choices additionally impact the convergence properties of the solver. The balance between computational scalability and convergence rate is carefully considered, to ensure high overall performance while maintaining robustness.
The basis-function kernel, which dominates the setup phase, and the local smoother, which dominates the solution phase, are studied in detail. In addition, the balance between convergence rate and scalability as a function of the coarsening ratio, Cr, is investigated.
Based on this analysis, the performance of a scalable AMS implementation is tested using highly heterogeneous problems derived from the SPE10 benchmark ( Christie et al., 2001 ) and geostatistically generated. The problems range in size from a few million to more than a 100 million cells. The solver robustness and scalability is demonstrated on modern multi-core systems and compared with the widely used parallel SAMG solver ( Stüben, 2012 ).