A novel algorithm for DNS and LES simulations of cavitating flows

Un nouvel algoritme pour la simulation DNS et LES des ecoulements cavitants

Anton Znidarcic (1)

1 LML - Laboratoire de Mécanique de Lille - FRE 3723

Cavitation-turbulence interactions are problematic aspect of cavitating flows which imposes limitations in development of better cavitation and turbulence models. DNS simulations with homogeneous mixture approach are proposed to overcome this and offer more insight into the phenomena. As DNS simulations are highly demanding and a variety of cavitation models exists, a tool devoted specifically to them is needed. Such tools usually demand application of highly accurate discretization schemes, direct solvers and multi domain methods enabling good scaling of the codes. As typical cavitating flow geometries impose limits on suitable discretization methods, compact finite differences offer the most appropriate discretization tool. The need for fast solvers and good code scalability leads to request for an algorithm, capable of stable and accurate cavitating flow simulations where solved systems feature multiplication of implicitly treated variables only by constant coefficients. A novel algorithm with such ability was developed in the scope of this work using Concus and Golub method introduced into projection methods, through which the governing equations for homogeneous mixture modeling of cavitating flows can be resolved. Work also proposes an effective and new approach for verification of the new and existing algorithms on the basis of Method of Manufactured Solutions.

Le couplage diphasique-turbulence est une propriété clé des écoulements cavitants, qui est un frein important à l’amélioration des modèles de cavitation et de turbulence. Réaliser des simulations directes (DNS) est le moyen proposé ici pour s’affranchir du modèle de turbulence et obtenir des informations nouvelles sur les phénomènes mis en jeu. Ce type de simulation est exigeant sur le plan numérique, et requiert le développement d’un solveur spécifique intégrant les spécificités des modèles de cavitation. Cela inclue notamment des schémas de discrétisation d’ordre élevé, un solveur direct, et une résolution multi-domaines associée à une parrallélisation efficace. Une discrétisation par différences compactes finies s’avère être le meilleur choix. La contrainte de rapidité et de parrallélisation impose un algorithme où les systèmes résoudre n’impliquent des multiplications des variables implicites que par des coefficients invariants au cours du calcul. Un nouvel algorithme réunissant ces critères a été développé durant cette thèse, à partir de la combinaison de la méthode de Concus & Golub et d’une méthode de projection, qui permet de résoudre les équations associées à la modélisation homogène de la cavitation. Une nouvelle approche de vérification de ce nouvel algorithme est également proposée et mise en œuvre sur la base de la méthode des solutions manufacturées (MMS).

Mots clés: Concus and Golub method; Compact finite differences; Dns; Cavitation; Homogeneous mixture approach; Influence matrix

Mots clés: Dns; Méthode Concus et Golub; Différences finies compactes; Matrice d'influence; Modèle homogène; Cavitation

Domaines: Energie électrique

Anton Znidarcic. Un nouvel algoritme pour la simulation DNS et LES des ecoulements cavitants. Energie électrique. Ecole nationale supérieure d'arts et métiers - ENSAM, 2016. Français. ⟨NNT : 2016ENAM0056⟩. ⟨tel-01535010⟩