An Algorithm for the lattice Boltzmann method adapted to sparse geometries and distributed memory systems

Dans cette thèse, une modification de l’algorithme de la méthode de Boltzmann sur réseau est effectuée pour obtenir de meilleures performances et par conséquent moins de temps d’exécution. L’objectif principal du projet est d’améliorer les algorithmes de Boltzmann sur réseau actuels et d’en créer un qui fonctionnera efficacement sur des systèmes à mémoire distribuée ainsi qu’avec un seul processeur, et ce, pour simuler de grands systèmes avec des géométries irrégulières. L’idée suivie dans la thèse est d’éviter d’emmagasiner les données des noeuds solides. L’autre amélioration de l’algorithme de Boltzmann sur réseau consiste à parcourir en une seule passe toutes les fonctions de distribution de probabilité afin d’effectuer les étapes de collision et de propagation en même temps. Ainsi, comme une des boucles principales est sur les fonctions de distribution de probabilité, il est possible d’utiliser un tableau de structures pour simplifier la mise en oeuvre de l’étape de propagation. En considérant les difficultés d’implémentation des conditions aux limites avec la méthode de Boltzmann sur réseau, un nouvel environnement est développé pour simplifier la mise en oeuvre des conditions limites. Une autre idée dans cette thèse est de concevoir un nouvel algorithme pour la communication de l’interface de passage de messages (MPI en anglais) de telle manière qu’au niveau de la programmation, il serait possible de cacher la couche de programmation parallèle des parties du code en série. Cependant, toute l’implémentation et la vérification sont faites pour des cas bidimensionnels. En suivant les idées mentionnées pour atteindre les objectifs du projet, un nouvel algorithme de Boltzmann sur réseau a été introduit qui utilise une disposition de mémoire simple. Comme les étapes de collision et de propagation sont effectuées en même temps, la localité des données est augmentée. Par conséquent, un simple décalage des données est utilisé pour effectuer l’étape de propagation. En évitant d’emmagasiner les noeuds solides, l’algorithme proposé est adapté aux systèmes à faible taux de vide. De cette manière, outre la diminution de l’utilisation de la mémoire, comme aucun calcul n’est effectué pour les noeuds solides, les ressources de calcul sont également utilisées efficacement. Pour vérifier le code implémenté et basé sur l’algorithme mentionné, plusieurs cas tests académiques sont simulés, notamment l’écoulement de Poiseuille, l’écoulement dans une cavité et l’écoulement autour d’un cylindre. L’algorithme proposé utilise la même quantité de mémoire que l’algorithme à deux réseaux, mais plus que l’algorithme à deux étapes. Puisqu’il utilise un schéma d’adressage indirect, il est néanmoins plus efficace en termes d’utilisation de la mémoire, en particulier pour les géométries à faible taux de vide. Comme les données des noeuds solides ne sont pas emmagasinées dans cette méthode, il n’y a pas de calcul pour ces types de noeuds, ce qui entraîne par conséquent de meilleures performances pour les géométries irrégulières avec de faibles taux de vide. En tenant compte des difficultés de mise en oeuvre des conditions limites, un nouvel environnement est introduit pour mettre en oeuvre de nouvelles conditions limites. Dans cet environnement, une condition limite est mise en oeuvre pour un seul côté/coin d’un domaine, puis appliquée aux autres côtés/coins de manière automatique. Ainsi, comme les formulations habituelles d’une condition limite dépendent de l’orientation des noeuds frontière, l’utilisation de cet environnement réduit le nombre d’implémentations, ce qui facilite l’implémentation des conditions limites typiques de la méthode de Boltzmann sur réseau. L’une des principales contributions de ce travail est le concept d’application de l’opération de décalage circulaire et de retournement de données sur un tableau de structure pour projeter les fonctions de distribution de probabilité de n’importe quel bord/coin vers le noeud frontière gauche/inférieur-gauche. Pour vérifier l’environnement, plusieurs des conditions limites communes de la méthode Boltzmann sur réseau ont été mises en oeuvre, puis vérifiées à l’aide de références bien connues dans le domaine de la méthode de Boltzmann sur réseau.