Simulation of non-newtonian flows in macroscopic porous media with the lattice-Boltzmann method

Un fluide non-newtonien est un fluide dont la relation entre son taux de déformation et la contrainte qui lui est imposé n'est pas linéaire. Dans un milieu poreux, la contrainte imposée au fluide dépend de la pression imposée mais aussi de la taille des pores, et donc de la perméabilité à l'échelle macroscopique. Certains fluides peuvent avoir une rhéologie qui présente un changement à partir d'un seuil en contrainte. Si la taille des pores est aléatoire, le fluide va alors changer de régime de façon hétérogène dans le milieu. L'écoulement pourra alors présenter un premier régime où tout le fluide est en dessous du seuil, un régime où tout le fluide est au-dessus du seuil, et un régime intermédiaire pour lequel les deux types de rhéologie coexistent. Nous nous intéressons à ce régime intermédiaire pour des écoulements de fluide non-newtonien en milieu poreux macroscopique, étudiés par des simulations. Plus particulièrement, nous regardons l'écoulement d'un fluide de Bingham et celui d'un fluide de Carreau. Le fluide de Bingham ne s'écoule qu'à partir d'une contrainte seuil. En dessous du seuil, il se comporte comme un solide. Au-delà, sa relation taux de déformation/contrainte suis une loi affine. Les fluides de Carreau ont une relation taux de déformation/contrainte qui change de régime entre une loi linéaire et une loi de puissance. L'étude à l'échelle macroscopique se fait en simulant une loi de Darcy-Brinkman dans un champ de perméabilité hétérogène. Nous utilisons pour nos simulations la méthode de Boltzmann sur réseau, avec une grille de nœuds régulière, plus particulièrement le schéma à deux temps de relaxation d'Irina Ginzburg. Pour chaque fluide, nous regardons la relation débit-pression, ainsi que les propriétés géométriques des différents régimes d'écoulement. Nous caractérisons plus particulièrement les propriétés multi-échelles des régions dans un le même régime d'écoulement (clusters), telles que leur taille ou leur forme. Nous faisons aussi le lien entre ces propriétés géométriques et la théorie de la percolation, qui étudie le comportement de cartes de nœuds s'ouvrant aléatoirement et qui prédit des propriétés fractales.A non-newtonian fluid is a fluid which relation between it's shear rate and the stress under which it is put, is not linear. In a porous medium, the stress imposed to the fluid depends on the imposed pressure, but also on the pores size, and therefore on the macroscopic scale permeability. Some fluids have a rheology such that the fluid show a change of behaviour reaching a yield stress. If the pore size is random, then the fluid will present heterogeneous regime changes in the medium. The flow will then show a first regime where the whole fluid will be under the threshold, a regime where the whole fluid will be far above the threshold, and an intermediate regime for which both rheologies coexists. We are interested in intermediate regime for the flow of non-newtonian fluids in macroscopic porous media, and study it with numerical simulations. More particularly, we look at the flow of a Bingham fluid and that of a Carreau fluid. The Bingham fluid doesn't flow under a yield stress. Under the threshold, it behaves as a solid. Beyond, it's shear-rate/stress relation is an affine law. Carreau fluids have a shear-rate/stress relation that change regime between that of a newtonian fluid, and a power law. The macroscopic scale study is done simulating a Darcy-Brinkman law in a heterogeneous permeability field. We use for our simulations the lattice-Boltzmann method, on a regular node grid, and more specifically Irina Ginzburg two relaxation-time scheme. For each fluid, we study the flow-pressure relationship, as well as the geometric properties and the multi-scale properties in the fluid regions in the same flow regime (clusters), properties such as their size and shape. We also link these geometrical properties to the percolation theory, which studies the behaviour of randomly opening node maps and predicts fractal properties