Développement de stratégies numériques basées sur la DEM et la CFD-DEM améliorant la performance computationnelle de simulations d'écoulements granulaires et polyphasiques

La conception et l’opération des unités industrielles nécessitent une compréhension du com-portement des poudres et des particules, qui peut être contre-intuitif par rapport à nos connaissances basées sur les liquides et les gaz. Les interactions entre les particules, ou avec le fluide dans le cas d’écoulement polyphasique, impliquent plusieurs considérations. La mé-thode des éléments discrets (DEM) est l’une des méthodes numériques les plus populaires et adaptées pour simuler ces écoulements. Le DEM permet un suivi distinct des particules selon la loi du mouvement et modélise les interactions de collision à l’aide de relations force-déplacement, au détriment de la performance. Cette méthode est souvent couplée avec un solveur de mécanique des fluides (CFD-DEM) pour des simulations d’écoulements polypha-siques. Puisque les simulations des systèmes industriels impliquent un nombre considérable de particules, allant jusqu’à plusieurs millions, l’utilisation d’approches améliorant la perfor-mance computationnelle des simulations est justifiée. Des stratégies de parallélisation, des méthodes basées sur les processeurs graphiques (GPU) et des techniques algorithmiques ont été proposées pour améliorer les performances de la DEM et de la CFD-DEM. Toutefois, le matériel informatique de haute performance n’est pas accessible à tous et peut représen-ter un investissement important. De plus, les méthodes algorithmiques actuelles sont parfois limitées à certaines spécificités du cas étudié. L’objectif principal de ce projet est le développement d’une nouvelle méthode pour améliorer les performances de calcul des écoulements granulaires et polyphasiques applicable à la DEM et à la CFD-DEM. L’hypothèse fondamentale repose sur l’identification de régions quasi statiques où les vitesses des particules présentent une faible variance entre les particules avoi-sinantes par une évaluation de la température granulaire. Dans ces zones, les déplacements des particules résultant de collisions n’affecteraient pas le comportement macroscopique de l’écoulement, suggérant l’omission du calcul de leur force de contact et de leur position. Le mécanisme a été développé afin de réduire la charge de calcul du DEM en minimisant les calculs résultant des collisions de particules tout en maintenant le calcul intégral dans les régions critiques. Il est applicable à la DEM et également à la CFD-DEM avec l’ajout d’un terme d’advection moyenné pour prendre en charge le mouvement en vrac. La désactivation adaptative des contacts de particules dans les systèmes dynamiques crée des zones inactives dispersées qui migrent ou disparaissent pendant la simulation. Ce phénomène a conduit au nom du mécanisme de l’algorithme : contacts creux adaptatifs ou Adaptive Sparse Contacts (ASC).