Évènements

Dans cet exposé, nous présenterons les schémas de bi-projection récemment développés pour la discrétisation en temps des équations de Navier-Stokes incompressibles isothermes avec une rhéologie visco-plastique (Bingham, Drucker-Prager). La difficulté principale, tant du point de vue de l’analyse mathématique que de l’approximation numérique, est due à  la non-différentiabilité de la partie plastique du tenseur des contraintes dans les régions de l’espace o๠le tenseur des déformations est nul. Une formulation basée sur une réécriture de la définition de la partie plastique du tenseur en terme d’une projection est utilisée. Un nouveau schéma de semi-discrétisation en temps, basé sur un schéma de projection incrémental classique pour les équations de Navier-Stokes dans le cas des fluides newtoniens, est proposé. Le tenseur plastique est traité en implicite dans l’étape de prédiction du schéma de projection et un algorithme de point fixe est utilisé pour son évaluation. Un terme de pseudo-relaxation temporel est ajouté dans la projection de Bingham afin d’assurer une convergence rapide (géométrique) du point fixe. Des analyses de stabilité et d’erreurs du schéma numérique ont été obtenues dans un premier travail dans le cas d’écoulements homogènes puis étendues au cas d’écoulements à  densité, viscosité et seuil de plasticité variables. Des résultats numériques d’écoulements dans une cavité entraînée 2D ainsi que d’instabilités de Rayleigh-Taylor seront présentés afin de montrer l’efficacité de la méthode de bi-projection. Enfin, des applications à  des écoulements géophysiques seront discutées.