Deflation of conjugate gradients with applications to boundary value problems
SIAM Journal on Numerical Analysis
A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
SIAM Journal on Scientific Computing
A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method)
Journal of Computational Physics
Multigrid
On the Construction of Deflation-Based Preconditioners
SIAM Journal on Scientific Computing
A Comparison of Deflation and Coarse Grid Correction Applied to Porous Media Flow
SIAM Journal on Numerical Analysis
A Comparison of Deflation and the Balancing Preconditioner
SIAM Journal on Scientific Computing
Recycling Subspace Information for Diffuse Optical Tomography
SIAM Journal on Scientific Computing
Recycling Krylov Subspaces for Sequences of Linear Systems
SIAM Journal on Scientific Computing
An accurate conservative level set/ghost fluid method for simulating turbulent atomization
Journal of Computational Physics
Fast and robust solvers for pressure-correction in bubbly flow problems
Journal of Computational Physics
Journal of Scientific Computing
Journal of Computational Physics
Hi-index | 31.45 |
The discretization of Partial Differential Equations often leads to the need of solving large symmetric linear systems. In the case of the Navier-Stokes equations for incompressible flows, solving the elliptic pressure Poisson equation can represent the most important part of the computational time required for the massively parallel simulation of the flow. The need for efficiency that this issue induces is completed with a need for stability, in particular when dealing with unstructured meshes. Here, a stable and efficient variant of the Deflated Preconditioned Conjugate Gradient (DPCG) solver is first presented. This two-level method uses an arbitrary coarse grid to reduce the computational cost of the solving. However, in the massively parallel implementation of this technique for very large linear systems, the coarse grids generated can count up to millions of cells, which makes direct solvings on the coarse level impossible. The solving on the coarse grid, performed with a Preconditioned Conjugate Gradient (PCG) solver for this reason, may involve a large number of communications, which reduces dramatically the performances on massively parallel machines. To this effect, two methods developed in order to reduce the number of iterations on the coarse level are introduced, that is the creation of improved initial guesses and the adaptation of the convergence criterion. The design of these methods make them easy to implement in any already existing DPCG solver. The structural requirements for an efficient massively parallel unstructured solver and the implementation of this solver are described. The novel DPCG method is assessed for applications involving turbulence, heat transfers and two-phase flows, with grids up to 17.8billion elements. Numerical results show a two- to 12-fold reduction of the number of iterations on the coarse level, which implies a reduction of the computational time of the Poisson solver up to 71% and a global reduction of the proportion of communication times up to 53%. As a result, the weak scaling of the LES solver is shown to be clearly improved for massively parallel uses.