Linearized Navier-Stokes for aeroelasticity, shape optimisation and aeroacoustics
Résolution des équations de Navier-Stokes linéarisées pour l'aéroélasticité, l’optimisation de forme et l’aéroacoustique
Résumé
The linearized Navier-Stokes equations are solved at Dassault Aviation within numerical simulations for aerodynamic shape optimisation, flutter calculations and aeroacoustics. In order to improve the robustness and efficiency of the Navier-Stokes solver, this thesis followed two complementary paths. The first is work on the iterative methods used to solve linear systems, and the second is the improvement of the numerical scheme underlying these linear systems. In the first part, the extension to multiple right-hand sides of the GMRES algorithm with spectral deflation was tested on industrial test cases. The use of the ILU(k) preconditioner within an additive Schwarz method led to a reduction of the time needed to solve the systems with GMRES by a factor ten. It also enabled the convergence of some numerically very difficult cases which could not be solved by the software available before this thesis. The second part of the manuscript begins with work on the SUPG method used to stabilise the finite element scheme. A new way of computing the stabilisation matrix gave promising results on non-linear cases, which were however not observed for linear cases. A study on Dirichlet boundary conditions concludes this part. An algebraic method to impose non homogeneous Dirichlet boundary conditions on non-trivial variables is introduced. It enables the use, in an industrial context, of linearized Navier-Stokes for aeroacoustics. Moreover, the transparent behaviour of a homogeneous Dirichlet boundary conditions on all variables is proved to be due to the SUPG stabilisation.
Les équations de Navier-Stokes linéarisées sont utilisées dans l’industrie aéronautique pour l’optimisation de forme aérodynamique, l’aéroélasticité et l’aéroacoustique. Deux axes ont été suivis pour accélérer et rendre plus robuste la résolution de ces équations. Le premier est l’amélioration de la méthode itérative de résolution de systèmes linéaires utilisée, et le deuxième la formulation du schéma numérique conduisant à ce système linéaire. Dans cette première partie, l’extension de l’algorithme GMRES avec déflation spectrale à des systèmes à plusieurs seconds membres a été testée sur des cas tests industriels. L’amélioration du préconditionnement de la méthode GMRES par l’utilisation d'une méthode de Schwarz additive avec préconditionneur ILU(k) a permis une accélération du temps de résolution allant jusqu’à un facteur dix, ainsi que la convergence de cas jusqu’alors impossibles à résoudre. La deuxième partie présente d’abord un travail sur la stabilisation SUPG du schéma élément fini utilisé. La forme proposée de la matrice de stabilisation, dite complète, a donné des résultats encourageants en non-linéaire qui ne se sont pas transposés en linéarisé. Une étude sur les conditions aux limites de Dirichlet clôt cette partie. Une méthode algébrique d’imposition de conditions non homogènes sur des variables non triviales du calcul, qui a permis l’application industrielle à l’aéroacoustique, y est détaillée. De plus, la preuve est apportée que le caractère transparent d’une condition de Dirichlet homogène sur toutes les variables s’explique par le schéma SUPG.
Origine : Version validée par le jury (STAR)