Schémas de discrétisation Compatible Discrete Operator pour les équations de Navier–Stokes d’un fluide incompressible en régime instationnaire

Compatible Discrete Operator schemes for the unsteady Navier–Stokes equations for incompressible flows

• Auteur / Autrice : Riccardo Milani
• Direction : Alexandre Ern
• Type : Thèse de doctorat
• Discipline(s) : Mathématiques
• Date : Soutenance le 16/12/2020
• Etablissement(s) : Paris Est
• Ecole(s) doctorale(s) : École doctorale Mathématiques, Sciences et Technologies de l'Information et de la Communication
• Partenaire(s) de recherche : Laboratoire : Centre d'enseignement et de recherche en mathématiques et calcul scientifique (Champs-sur-Marne, Seine-et-Marne)
• Jury : Président / Présidente : Roland Becker
• Jury : Examinateurs / Examinatrices : Alexandre Ern, Raphaèle Herbin, Peter Minev, Pierre Cantin, Stella Krell, Jérôme Bonelle
• Jury : Rapporteurs / Rapporteuses : Raphaèle Herbin, Peter Minev

Mots clés

Mots clés contrôlés

Maillages polyédriques

Navier-Stokes, Équations de

Fluides incompressibles

Mots clés libres

Schémas CDO

Maillages polyédriques

Navier–Stokes

Compressibilité Artificielle

Condition inf-Sup

Convection

CDO schemes

Polyhedral meshes

Navier–Stokes

Artificial Compressibility

Inf-Sup condition

Convection

Résumé

Nous développons des schémas dits face-based Compatible Discrete Operator (CDO-Fb) pour les équations de Stokes et Navier–Stokes incompressibles en régime instationnaire. Des opérateurs pour la reconstruction du gradient, de la divergence et un autre pour le terme de convection sont proposés. On montre que l'opérateur de divergence discret permet de satisfaire une condition inf-sup, tandis que l'opérateur de convection discret est dissipatif, propriété cruciale pour le bilan d'énergie. Le schéma de discrétisation est d'abord testé dans le cas stationnaire sur des maillages généraux mais aussi déformés, afin d'illustrer la flexibilité et la robustesse de la discrétisation CDO-Fb. Dans un deuxième temps, l'attention est placée sur les techniques de marche en temps. En particulier, nous étudions l'approche monolithique traditionnelle qui consiste à résoudre directement le système de point-selle, et la méthode de Compressibilité Artificielle (AC), qui permet de ne plus avoir un système de point-selle à résoudre au prix d'une relaxation du bilan de masse. Trois stratégies classiques pour le traitement du terme non linéaire dû à la convection sont examinées: l'algorithme de Picard, la linéarisation et l'explicitation. Des résultats numériques utilisant des schémas temporels du premier ordre d'abord, puis du deuxième ordre, montrent que la méthode AC constitue une alternative précise et efficace à l'approche monolithique traditionnelle

We develop face-based Compatible Discrete Operator (CDO-Fb) schemes for the unsteady, incompressible Stokes and Navier–Stokes equations. We introduce operators discretizing the gradient, the divergence, and the convection term. It is proved that the discrete divergence operator allows one to recover a discrete inf-sup condition. Moreover, the discrete convection operator is dissipative, a paramount property for the energy balance. The scheme is first tested in the steady case on general and deformed meshes in order to highlight the flexibility and the robustness of the CDO-Fb discretization. The focus is then moved onto the time-stepping techniques. In particular, we analyze the classical monolithic approach, consisting in solving saddle-point problems, and the Artificial Compressibility (AC) method, which allows one to avoid such saddle-point systems at the cost of relaxing the mass balance. Three classic techniques for the treatment of the convection term are investigated: Picard iterations, the linearized convection and the explicit convection. Numerical results stemming from first-order and then from second-order time-schemes show that the AC method is an accurate and efficient alternative to the classical monolithic approach