Méthodes numériques adaptées à la résolution des équations de Navier-Stokes

Numerical methods suitable for solving the Navier-Stokes equations

• Auteur / Autrice : Yann Guevel
• Direction : Jean-Marc Cadou
• Type : Thèse de doctorat
• Discipline(s) : Sciences pour l'Ingénieur
• Date : Soutenance le 15/01/2016
• Etablissement(s) : Lorient
• Ecole(s) doctorale(s) : École doctorale Santé, information-communication et mathématiques, matière (Brest, Finistère)
• Partenaire(s) de recherche : COMUE : Université Bretagne Loire (2016-2019)
• Partenaire(s) de recherche : Laboratoire : Laboratoire d'Ingénierie des Matériaux de Bretagne / LIMATB
• Jury : Examinateurs / Examinatrices : Grégory Girault, Denis Sipp, Michel Potier-Ferry, Patrice Le Gal
• Jury : Rapporteurs / Rapporteuses : Laurette S. Tuckerman, Jean-Christophe Robinet

Mots clés

Mots clés contrôlés

Méthode asymptotique numérique

Modèles mathématiques

Écoulement instationnaire (dynamique des fluides)

Équations de Navier-Stokes -- Solutions numériques

Mots clés libres

Bifurcation stationnaire

Brisure de symétrie

Écoulements tridimensionnels

Steady bifurcation

Symmetry breaking

Threedimensional flow

Résumé

Le groupe de recherche Instabilités et Méthodes Numériques Spécifiques mène ses activités dans le développement d’outils numé- riques pour la résolution de problèmes non linéaires en utilisant, no- tamment, la Méthode Asymptotique Numérique (MAN). Basée sur le couplage d’une méthode de perturbation et de discrétisation spa- tiale, la MAN est efficace et permet de déterminer précisément les transitions telles que, par exemple, la perte d’unicité de la solution. L’objectif de ce travail de thèse est de proposer des méthodes numé- riques alternatives à la fois robustes, performantes pour la résolution des équations de Navier-Stokes. Nous nous intéressons à l’analyse de bifurcation stationnaire, mais aussi à la simulation d’écoulement dépendant du temps. Dans un premier temps, des techniques d’analyse de bifurcation nu- mérique pour des problèmes stationnaires à très grand nombre de degrés de liberté sont décrites. Nous implémentons ces techniques, basées sur la MAN, dans le logiciel open-source multi-physique ELMER . Nous détaillons l’implémentation des méthodes d’analyse de bifurcation stationnaire telles que la continuation de branches solutions, les techniques de détection des valeurs critiques du pa- ramètre de charge et les changements de branche en un point de bifurcation stationnaire. L’émergence d’une progression géométrique dans les termes de séries MAN à l’approche d’une singularité est dé- crite. Des discussions sont proposées pour le cas de bifurcations par brisure de symétrie. Les méthodes proposées dans ce travail sont validées en utilisant des cas référencés dans la littérature, tels que des écoulements dans des conduites à expansion/contraction sou- daine. Une étude paramétrique permet de présenter de nouveaux ré- sultats pour les écoulements tridimensionnels dans une expansion brusque. L’utilisation de librairies de calculs intensifs rend possible la réalisation d’analyse de bifurcation pour des modèles à très grand nombre de degrés de liberté, en des temps de calcul abordables. Dans un deuxième temps, des solveurs d’ordre élevé sont proposés pour la simulation d’écoulements instationnaires. Une technique d’homotopie à combinaison convexe et une technique de pertur- bation, sont couplées à un schéma d’intégration temporelle pour résoudre les équations instationnaires de Navier-Stokes. Le cas d’un écoulement bidimensionnel autour d’un cylindre fixe est étudié. Ce problème de référence nous permet de valider et discuter des amélio- rations proposées. De cette manière, nous confirmons, au cours des essais numériques, qu’il est possible de réduire les temps de cal- cul en évitant des assemblages d’opérateurs et des résolutions de systèmes linéaires qui n’apportent aucune information supplémen- taire pour la qualité des solutions. De plus, un nouvel éclairage est apporté sur l’utilisation des approximants de Padé par rapport aux travaux antérieurs. L’utilisation de ces solveurs non linéaires nous permet de réduire significativement le nombre de factorisations de matrice en les conservant valides pour un grand nombre de pas de temps, et parfois sur le domaine temporel complet. De nombreuses perspectives sont envisagées, notamment pour l’analyse des séries pour le cas d’un point limite, la bifurcation de Hopf, l’étude d’autre cas d’écoulements tridimensionnels, le couplage fluide-structure. De même, l’association des techniques MAN aux techniques de réductions de modèles et l’analyse de stabilité des orbites périodiques sont envisageables.

The research group "Instabilités et Méthodes Numériques Spéci-fiques" operates in the development of numerical tools for solving nonlinear problems by using, in particluar, the Asymptotic Numer- ical Method (ANM). Based on coupling a perturbation method and a spatial discretization, the ANM is effective and makes it possible to precisely determine the transitions such as, for example, loss of uniqueness of the solution. The objective of this thesis is to offer al- ternative numerical methods both robust and effective, for solving the Navier-Stokes equations. We are interested in steady bifurcation analysis, and in time dependent flow simulation .Initially, numerical bifurcation analysis techniques for steady flow problems in very large number of degrees of freedom are de- scribed. These techniques, based on the ANM, are implemented in the multiphysics ELMER open-source software. We detail the im- plementation of the steady bifurcation analysis methods such as continuation of solution branches, detection of load parameter critical values and branch switching at steady bifurcation point. The emer- gence of a geometric progression in ANM series terms in the vicinity of a singularity is described. Discussions are proposed for the case of symmetry breaking bifurcations. The methods described in this the- sis are validated using reference cases of the literature, such as flow in pipe with sudden expansion/contraction. New results for three- dimensional flow in a sudden expansion, are obtained according to a parametric study. The use of high performance computing libraries makes possible the bifurcation analysis for models with high number of degrees of freedom, in affordable computing times. Secondly, high-order solvers are proposed for the simulation of un- steady flows. Homotopy with convex combination and a perturba- tion technique, are coupled to a time integration scheme in order to solve the unsteady Navier-Stokes equations. The case of two- dimensional flow around a fixed cylinder is studied. This reference problem allows us to validate and discuss proposed improvements. In this way, we confirm, in the numerical tests, that it is possible to reduce the computation time by avoiding operators assembly and resolution of unuseful linear systems in respect to the solution quality. In addition, new lighting is provided on the use of Padé approximants over previous work. The use of these nonlinear solvers allows us to significantly reduce the number of matrix factorization retaining them valid for many time steps, and sometimes on the complete time do- main. Many opportunities are envisaged, in particular the analysis of ANM series for the case of limit point, the Hopf bifurcation, the study of other cases of three-dimensional flow, the fluid-structure interaction. Similarly, the combination of ANM models with reduction techniques f stable periodic orbits are possible.