Name
Abstract | ||
|---|---|---|
The most popular iterative linear solvers in Computational Fluid Dynamics (CFD) calculations are restarted GMRES and BiCGStab. At the beginning of most incompressible flow calculations, the computation time and the number of iterations to converge for the pressure Poisson equation are quite high. In this case, the BiCGStab algorithm, with relatively cheap but non-optimal iterations, may fail to converge for stiff problems. Thus, a more robust algorithm like GMRES, which guarantees monotonic convergence, is preferred. To reduce the large storage requirements of GMRES, a restarted version - GMRES(m) or its variants - is used in CFD applications. However, GMRES(m) can suffer from stagnation or very slow convergence. For this reason, we use the rGCROT method. rGCROT is an algorithm that improves restarted GMRES by recycling a selected subspace of the search space from one restart of GMRES(m) to the next as well as building and recycling this outer vector space from one problem to the next (subsequent time steps in this context). In the current work, we apply both GMRES and Bi-Lanczos based recycling to CFD simulations. For a turbulent channel flow problem with rGCROT we get comparable performance to BiCGStab. The rGCROT performance is still not as cheap as the BiCGStab algorithm in terms of both storage requirements and time to solution. Thus we have used a novel hybrid approach in which we get the benefits of the robustness of the GMRES(m) algorithm and the economical iterations of BiCGStab. For the first few time steps, the algorithm builds an outer vector space using rGCROT and then it recycles that space for the rBiCGStab solver. Time to solution with this approach is comparable to that for the BiCGStab algorithm with savings in the number of iterations and improved convergence properties. Overall, the hybrid approach also performs better than optimized rGCROT and GMRES(m). |
Year | Venue | Field |
|---|---|---|
2015 | arXiv: Numerical Analysis | Convergence (routing),Mathematical optimization,Generalized minimal residual method,Biconjugate gradient stabilized method,Linear subspace,Robustness (computer science),Solver,Computational fluid dynamics,Mathematics,Computation |
DocType | Volume | Citations |
Journal | abs/1501.03358 | 0 |
PageRank | References | Authors |
0.34 | 0 | 5 |
Name | Order | Citations | PageRank |
|---|---|---|---|
| Amit Amritkar | 1 | 0 | 0.34 |
| Eric de Sturler | 2 | 398 | 27.32 |
| Katarzyna Swirydowicz | 3 | 5 | 1.35 |
| Danesh Tafti | 4 | 0 | 0.34 |
| Kapil Ahuja | 5 | 0 | 0.34 |