Modélisation multi-échelle du transport électrocinétique en milieu poreux chargé

As part of the disposal of nuclear waste in deep geological formations (Cigéo project), it is necessary to predict the flow of radionuclides in ionic form in charged porous media such as clays. The clay materials are complex and it is difficult to know their structure, especially at the nanoscale where electrokinetic phenomena become paramount. In this case, the strategy adopted in this work is to represent the porosity with a network of pores interconnected by channels. This requires to know the properties of the coupled transport at the channel scale taking into account the flow of solvent, solutes and charges under the influence of gradients of pressure, salt concentration and electrostatic potential. We used the Poisson-Nernst-Planck model to describe the electrokinetic transport in the channels. It can be shown that the behavior of these properties may vary between the channel scale and the sample scale (random networks constructed from parameters distributions). For example, it is found that the permeability of a network decreases as electrokinetic effects increase, which is not the case at the level of a single channel. These differences arise from the simultaneous presence of coupling between the flows and the heterogeneity of the sample. A numerical sample having the same properties as real clay sample has been created following a process of selection of the channels parameters, it has been used to conduct preliminary studies on clogging and desaturation in clay.Dans le cadre du stockage des déchets nucléaires en couche géologique profonde (projet Cigéo), il est nécessaire de pouvoir prédire l'écoulement des radionucléides sous forme ionique dans des milieux poreux chargés tels que l'argile. Les matériaux argileux sont complexes et il est difficile de connaître leur structure, surtout à l'échelle nanoscopique où les phénomènes électrocinétiques deviennent primordiaux. Dans ce cas, la stratégie adoptée dans ce travail est de représenter la porosité à l'aide d'un réseau de pores connectés entre eux par des canaux. Cela nécessite de connaître les propriétés du transport couplé à l'échelle du canal en prenant en compte les flux de solvant, de solutés et de charges sous l'influence de gradients de pression, de concentration en sel et de potentiel électrostatique. La description du transport électrocinétique utilisée correspond à celle du modèle de Poisson-Nernst-Planck. On peut ainsi montrer que le comportement de ces propriétés peut différer entre l'échelle du canal et celle de l'échantillon (réseaux construits aléatoirement à partir de distributions de paramètres). Par exemple, on constate que la perméabilité d'un réseau diminue si les effets électrocinétiques augmentent, ce qui n'est pas le cas à l'échelle du canal. Ces différences proviennent de la présence simultanée de couplages entre les flux et de l'hétérogénéité du milieu. Un échantillon numérique ayant les mêmes propriétés qu'un échantillon réel d'argile a ainsi été créé au terme d'un processus de sélection des distributions de paramètres des canaux, il a ainsi servi à réaliser des études préliminaires sur le colmatage et la désaturation