Spectral recycling strategies for the solution of nonlinear eigenproblems in thermoacoustics

Pablo Salas (1) , Luc Giraud (2) , Yousef Saad (3) , Stéphane Moreau (4)

1 UdeS - Département de génie mécanique [Sherbrooke]
2 HiePACS - High-End Parallel Algorithms for Challenging Numerical Simulations
3 UMN - University of Minnesota [Twin Cities]
4 UdeS - Université de Sherbrooke

In this work we consider the numerical solution of large nonlinear eigenvalue problems that arise in thermoacoustic simulations involved in the stability analysis of large combustion devices. We briefly introduce the physical modeling that leads to a non-linear eigenvalue problem that is solved using a nonlinear fixed point iteration scheme. Each step of this nonlinear method requires the solution of a complex non-Hermitian linear eigenvalue problem. We review a set of state of the art eigensolvers and discuss strategies to recycle spectral informations from one nonlinear step to the next. More precisely, we consider the Implicitly Restarted Arnoldi method, the Krylov-Schur solver and its block-variant as well as the subspace iteration method with Chebyshev acceleration. On a small test example we study the relevance of the different approaches and illustrate on a large industrial test case the performance of the parallel solvers best suited to recycle spectral information.

Domaines: Analyse numérique [math.NA]; Acoustique [physics.class-ph]; Mécanique des fluides [physics.class-ph]

Pablo Salas, Luc Giraud, Yousef Saad, Stéphane Moreau. Spectral recycling strategies for the solution of nonlinear eigenproblems in thermoacoustics. Numerical Linear Algebra with Applications, 2015, 22 (6), pp.1039-1058. ⟨10.1002/nla.1995⟩. ⟨hal-01238263⟩