L'hétérogénéité de la zone non saturée joue un rôle important dans le transfert d'eau et de soluté car elle accentue à la fois le développement des zones de stockage temporelles et les écoulements préférentiels. Par conséquent, la validation des modèles prédictifs nécessite le développement des outils expérimentaux spécifiques afin d'observer et de quantifier les mécanismes de transport impliqués dans un système non saturé hétérogène. Cette thèse vise à étudier l'effet combiné de la vitesse d'infiltration, de la barrière capillaire et l'angle de la pente d'interface entre deux matériaux sur les processus de l'écoulement de l'eau et du transport de soluté dans un modèle physique, le lysimètre de laboratoire 1x1x1.6 m3, nommé LUGH (Lysimeter for Urban Groundwater Hydrology) et un modèle numérique 3D de ce lysimètre. Le lysimètre LUGH est rempli par un sable fin et un mélange bimodal (50 % sable fin et 50 % gravier) en deux configurations: un profil uniforme de matériau bimodal ou un profil avec deux couches avec une pente de 14o. Ces agencements figurent l'hétérogénéité structurale et texturale observée sur un des sites expérimentaux de l'OTHU (Observatoire de Terrain en Hydrologie Urbaine) : le bassin d'infiltration d'eaux pluviales Django Reinhardt géré par la ville de Lyon (France). Le lysimètre est alimenté en eau et avec un traceur inerte (bromure de potassium, KBr) sur une partie de la surface par un système d'arrosage automatique. Les effluents ont été recueillis dans quinze sorties différentes en bas du lysimètre. La forte hétérogénéité des flux des sorties et des courbes de percée souligne la mise en place des écoulements préférentiels résultant à la fois de l'effet de barrière capillaire et de l'effet de fond du lysimètre. A partir des résultats expérimentaux, la modélisation numérique à l'aide de logiciel COMSOL MultiphysicsTM a permis de mieux comprendre les mécanismes responsables de ces transferts hétérogènes. Lorsque le modèle numérique validé, un test de sensibilité a été conduit pour étudier les effets de la vitesse d'infiltration et de la pente de l'interface sur l'écoulement. Les résultats montrent que la diminution de la vitesse d'infiltration ou l'augmentation de la pente de l'interface favorisent le développement des écoulements préférentiels. Notre étude a donné également des renseignements pertinents sur le couplage entre les processus hydrodynamiques et le transfert des solutés dans les sols non saturés hétérogènes en soulignant le rôle de la géométrie des interfaces ainsi que celui des conditions aux limites comme des facteurs clés pour la quantification des écoulements préférentiels.

Les systèmes électrochimiques lithium-air sont des concepts naissants mais exhibent des performances théoriques intéressantes qui laissent espérer une rupture technologique dans le domaine des batteries pour véhicule électrique. La possibilité d'atteindre une densité d'énergie supérieure à 500 Wh kg-1 est effectivement en ligne de mire. A contrario de la technologie lithium-air anhydre, les systèmes lithium-air aqueux n'ont, jusqu'à présent, fait l'objet d'aucune étude approfondie. Ce travail concerne donc le développement d'un système lithium-air aqueux, à trois électrodes, et vise également à améliorer nos connaissances fondamentales dans le domaine. La présente étude se focalise sur le compartiment positif de la cellule, dans lequel les réactions de l'oxygène sont mises en jeu. Dans un premier temps, une électrode spécifiquement dédiée à la réaction de dégagement d'oxygène a été élaborée à partir d'un acier 316L. L'étude de son comportement a révélé une bonne propension à catalyser la réaction de dégagement d'oxygène ainsi qu'une bonne stabilité sur 3 000 heures de fonctionnement. Néanmoins, d'importants problèmes de catalyse ont pu être observés et attribués à la présence des ions lithium dans l'électrolyte de la batterie. Les ions Li+ bloquent les transitions électrochimiques des sites actifs à l'origine des propriétés d'électrocatalyse. Le comportement d'une électrode à air, composée de carbone et d'oxydes de manganèse, a par la suite été caractérisé dans ce milieu. L'étude révèle deux phénomènes importants réduisant la performance de l'électrode et dont l'origine a également été attribuée aux ions lithium : un blocage des transitions (MnIII/MnIV), et une stabilisation des groupements oxygénés à la surface du carbone. Pour finir, il a été proposé d'optimiser le système électrolytique en limitant l'activité des ions Li+ en solution et ainsi d'améliorer le rendement en potentiel de charge/décharge de la batterie.

Les glissements de terrain sont les risques naturels fréquents dans le monde. D'une part, ils concernent un comportement alterné pour les géomatériaux : le comportement solide décrit classiquement par la géomécanique, et celui fluide après la rupture. D'autre part, ils peuvent être induits par plusieurs facteurs. Par exemple, la géologie, la topographie, la végétation, la variationclimatique, et les activities humaines. Dans les régions tropiques et tempérées, il est bien connu que les glissements de terrain surviennent fréquemment après les pluies intensives, à cause de la saturation des sols non saturés.Comme jusqu'à présent, il n'y a aucun modèle constitutif satisfaisant pour décrire une telle transition pour les géomatériaux non saturés, nos travaux se concentrent d'abord à un modèle unifié, qui associe un modèle hydro-élasto-plastique pour les géomatériaux non saturés, un modèle visqueux de Bingham, et un critère de transition entre les phases solide et fluide. Le modèle unifié permet dedécrire complètement le comportement avant et après la rupture, induite par des conditions hydromécaniques, pour les géomatériaux non saturés dans un cadre unifié.Ensuite, la Méthode des Eléments Finis avec les Points d'Intégration Lagrangiens (MEFPIL), qui bénéficie des avantages des approches Lagrangiennes (suivre les variables internes) et de celles Eulériennes (traiter les grandes transformations), est améliorée et choisie pour effectuer les calculs dans ce cadre. Nous avons introduit une nouvelle formulation hydro-élasto-plastique avec le couplage hydromécanique basée sur la MEFPIL, et l'implanté dans l'Ellipsis (le code basé sur la MEFPIL). En outre, pour inverser la matrice non symétrique de la formulation, un nouveau solver basé sur la méthode stabilisée du gradient bi-conjugué (BiCGSTAB) a été aussi implanté.Enfin, plusieurs benchmarks sont proposés pour valider le comportement mécanique des géomatériaux non saturés décrit par le modèle, et un glissement de terrain heuristique induit par la pluie est simulé, nous pouvons étudier les effets de paramètres divers sur l'infiltration d'eau et sur la rupture. En outre, une simulation du glissement de terrain à Shien (Chine, 2012), avec des paramètres réels est présentée. / Landslides are common natural disasters all over the world. On one hand, they implicate an alternate behaviour for geomaterials: the solid behaviour dealt with classically by geomechanics, and the fluid behaviour after failure. On the other hand, they can be induced by several factors. For example,geology, topography, vegetation, climatic variation, and human activities. In tropic and tempered regions, it is well known that the landslides occur frequently after intense rainfalls, because of saturation effect of unsaturated soils.As there is, up to now, no satisfactory constitutive model to describe such a transition for unsaturated geomaterials, our work is thus focused firstly on a unified model, that associates an hydro-elastoplastic model for unsaturated geomaterials, a Bingham's viscous law, and a transition criterion between solid and fluid states. The model allows to describe both solid and fluid behaviours for unsaturated geomaterials in a unified framework, and is made possible to simulate completely the behaviour of unsaturated geomaterials before and after failure, induced by hydromecanical conditions.Secondly, the Finite Element Method with Lagrangian Integration points (FEMLIP), that benefits from both the Lagrangian approaches (track internal variables) and Eulerian approaches (solve large transformations), is ameliorated and chosen to carry out calculations in this framework. We have introduced a new hydro-elasto-plastic FEMLIP formulation with hydromecanical coupling, and implemented it in Ellipsis (FEMLIP based code). Moreover, in order to inverse nonsymmetric matrixin the formulation, a new solver based on the biconjugate gradient stabilized method (BiCGSTAB) has been also implemented.Finally, several benchmarks are proposed to validate the model for the main features of unsaturated geomaterials, and a heuristic rainfall-induced landslide is simulated, we could study the effects of various parameters on the water infiltration and on the failure. Moreover, a simulation of the Shien landslide in China (2012) with real parameters is presented.

Les systèmes électrochimiques lithium-air sont des concepts naissants mais exhibent des performances théoriques intéressantes qui laissent espérer une rupture technologique dans le domaine des batteries pour véhicule électrique. La possibilité d'atteindre une densité d'énergie supérieure à 500 Wh kg-1 est effectivement en ligne de mire. A contrario de la technologie lithium-air anhydre, les systèmes lithium-air aqueux n'ont, jusqu'à présent, fait l'objet d'aucune étude approfondie. Ce travail concerne donc le développement d'un système lithium-air aqueux, à trois électrodes, et vise également à améliorer nos connaissances fondamentales dans le domaine. La présente étude se focalise sur le compartiment positif de la cellule, dans lequel les réactions de l'oxygène sont mises en jeu. Dans un premier temps, une électrode spécifiquement dédiée à la réaction de dégagement d'oxygène a été élaborée à partir d'un acier 316L. L'étude de son comportement a révélé une bonne propension à catalyser la réaction de dégagement d'oxygène ainsi qu'une bonne stabilité sur 3 000 heures de fonctionnement. Néanmoins, d'importants problèmes de catalyse ont pu être observés et attribués à la présence des ions lithium dans l'électrolyte de la batterie. Les ions Li+ bloquent les transitions électrochimiques des sites actifs à l'origine des propriétés d'électrocatalyse. Le comportement d'une électrode à air, composée de carbone et d'oxydes de manganèse, a par la suite été caractérisé dans ce milieu. L'étude révèle deux phénomènes importants réduisant la performance de l'électrode et dont l'origine a également été attribuée aux ions lithium : un blocage des transitions (MnIII/MnIV), et une stabilisation des groupements oxygénés à la surface du carbone. Pour finir, il a été proposé d'optimiser le système électrolytique en limitant l'activité des ions Li+ en solution et ainsi d'améliorer le rendement en potentiel de charge/décharge de la batterie. / The electrochemical lithium-air devices are emerging concepts and their very high theoretical performances have attracted a lot of attention, especially for an application in the electrical vehicle. A target of at least 500 Wh kg-1 is aimed for. The aqueous lithium-air devices have not yet been studied in detail which is not the case for the anhydrous lithium-air technology. This thesis firstly deals with the development of an aqueous lithium-air cell based on a three electrodes setup, and secondly attempts to improve our theoretical knowledge of these systems. This study particularly focuses on the positive compartment of the cell in which oxygen reactions occur. The first section is dedicated to the development of an oxygen evolution electrode made of 316L stainless steel. The study shows its ability to catalyze the oxygen evolution reaction as well its good stability over 3 000 hours of operation. Nevertheless, major problems of catalysis were observed and assigned to the presence of lithium ions in the electrolyte. Li+ ions inhibit the electrochemical oxidation of the active sites, which are needed for the electrocatalytic properties. The behavior of an air cathode was characterized in the same medium. The results show two important phenomena which reduces the electrode performance and which are also attributed to lithium ions : a deactivation of the (MnIII/MnIV) transition and a stabilization of the oxygenated groups at the carbon surface. Finally, an optimization of the electrolytic system was proposed by limiting the Li+ ions activity in solution, which improves the charge/discharge potential efficiency of the battery.

Les acides gras saturés à chaîne moyenne (AGCM) ont été décrits comme potentiellement bénéfiques pour réguler la balance énergétique chez l’homme. Cependant, la récente découverte de l’activation de la ghréline, hormone gastrique orexigène, par l’acide octanoïque (C8:0), ajoute un nouveau niveau de complexité dans la compréhension des effets physiologiques des AGCM. En effet, pour être reconnue par son récepteur hypophysaire, la ghréline doit préalablement être octanoylée par la Ghréline O-Acyltransférase (GOAT). Or, les paramètres régulés par la ghréline ainsi modifiée s’opposent aux effets connus des AGCM, en favorisant la prise alimentaire et la prise de masse. L’élaboration d’inhibiteurs de la GOAT constitue donc une stratégie thérapeutique nouvelle pour lutter contre l’obésité. Les AGCM sont essentiellement apportés par l’alimentation,leur biosynthèse endogène n’ayant été décrite que dans les cellules de la glande mammaire en lactation. Etudier l’impact du C8:0 d’origine alimentaire sur l’octanoylation de la ghréline est donc d’autant plus pertinent, car celui-ci est suspecté de fournir à la GOAT le co-substrat nécessaire à l’octanoylation. Afin de mieux comprendre l’apparente contradiction entre les effets connus des AGCM sur la masse corporelle et la satiété, et l’effet de l’octanoylation de la ghréline, quatre approches ont été mises en œuvre: (1) Etude des effets physiologiques du C8:0 à différentes doses sur l’axe ghréline-GOAT et sur le métabolisme chez le rat adulte, (2) Etude de l’absorption gastrique du C8:0 par inhibition de la lipase préduodénale, chez le / Medium chain fatty acids (MCFA) have been described as potentially beneficial for regulating human energetic balance. However, the recent discovery of the orexigenic ghrelin activation by octanoic acid revealed a new level of complexity in MCFA effects. Ghrelin must be octanoylated by the Ghrelin O-Acyltransferase (GOAT) to bind its pituitary receptor. Interestingly, the effects mediated by the post-translationally modified ghrelin seem opposed to MCFA known effects, increasing food intake and body weight. Development of GOAT inhibitors could then constitute a therapeutic strategy against obesity. MCFA are mainly provided by the diet, as their endogenous synthesis has only been described in lactating mammary glands. Dietary caprylic acid is then suspected to provide the GOAT enzyme with octanoyl- CoA co-substrates necessaryfor the acyl modification of ghrelin. My objective is to understand the discrepancy between the formerly described beneficial effects of dietary MCFAs on body weight loss and the C8:0 newly reported effect on appetite stimulation via ghrelin octanoylation. This thesis is divided into four parts: (1) Study of C8:0 physiological effects on ghrelin-GOAT axis, in rat, (2) Study of C8:0 gastric absorption by inhibition of preduodenal lipase, in young rat, (3) Study of C8:0 metabolism in gastric cells in vitro, (4) Preliminary study of GOAT activity in vitro.

La construction des biosystèmes d’oxydation passive du méthane (BOPM) est une option économique et durable pour réduire les émissions de méthane des sites d’enfouissement de déchets et des effets subséquents du réchauffement climatique. Les BOPM sont constitués de deux couches principales: la couche d'oxydation du méthane (MOL) et la couche de distribution du gaz (GDL). L'oxydation du méthane se produit dans la MOL par les réactions biochimiques des bactéries méthanotrophes, et la GDL est construite sous la MOL pour intercepter et distribuer les émissions fugitives de biogaz à la base de la MOL. Fondamentalement, l'efficacité d'un BOPM est définie en fonction de l'efficacité d'oxydation du méthane dans la MOL. Par conséquent, il est indispensable de fournir des conditions adéquates pour les activités bactériennes des méthanotrophes. En plus des paramètres environnementaux, l'intensité et la distribution du biogaz influencent l'efficacité des BOPM, et ils peuvent rendre le matériau de la MOL - avec une grande capacité d'accueillir les activités bactériennes - inutilisables en termes d'oxydation du méthane sur place. L'effet de barrière capillaire le long de l'interface entre la GDL et la MOL peut provoquer des émissions localisées de méthane, due à la restriction ou la distribution non uniforme de l’écoulement ascendant du biogaz à la base de la MOL. L'objectif principal de cette étude est d'incorporer le comportement hydraulique non saturé des BOPM dans la conception des BOPM, afin d’assurer la facilité et la distribution adéquates de l'écoulement du biogaz à la base de la MOL. Les fonctions de perméabilité à l'air des matériaux utilisés pour construire la MOL des BOPM expérimentaux au site d’enfouissement des déchets de St Nicéphore (Québec, Canada), ainsi que celles d'autres de la littérature technique, ont été étudiés pour évaluer le comportement d'écoulement non saturé du gaz dans les matériaux et pour identifier le seuil de migration sans restriction du gaz. Ce dernier seuil a été introduit en tant que un paramètre de conception avec lequel le critère de conception recommandé ici, c’est à dire la longueur de la migration sans restriction de gaz (LMSG), a été défini. La LMSG est considérée comme la longueur le long de l'interface entre la GDL et la MOL où le biogaz peut migrer à travers la MOL sans restriction. En réalisant des simulations numériques avec SEEP/W, les effets de la pente de l'interface, des paramètres définissant la courbe de rétention d'eau, de la fonction de la conductivité hydraulique du matériau de la MOL sur la valeur de la LMSG (représentant la facilité d'écoulement du biogaz à l'interface) et de la distribution de l'humidité (et par conséquent celle du biogaz) ont été évalués. Selon les résultats des simulations, la conductivité hydraulique saturée et la distribution des tailles de pores du matériau de la MOL sont les paramètres les plus importants sur la distribution de l'humidité le long de l'interface. Ce dernier paramètre influe également sur la valeur du degré de saturation et donc la facilité du biogaz à la base de la MOL. La densité sèche du matériau de MOL est un autre paramètre qui contrôle la facilité d'écoulement ascendant du biogaz. Les limitations principales de la présente étude sont associées au nombre de matériaux de MOL testés et à l'incapacité de SEEP/W de considérer l'évapotranspiration. Toutefois, compte tenu des hypothèses raisonnables dans les simulations et en utilisant les données de la littérature, on a essayé de réduire ces limitations. En utilisant les résultats des expériences et des simulations numériques, des étapes et des considérations de conception pour la sélection du matériau de MOL et de la pente d'interface ont été proposées. En effet,le comportement hydraulique non saturé des matériaux serait intégré dans les nécessités de conception pour un BOPM efficace, de sorte que la capacité maximale possible d'oxydation du méthane du matériau de la MOL soit exploitée. / Abstract : Implementation of passive methane oxidation biosystems (PMOB) is a cost effective and sustainable solution to eliminate the methane emissions of landfills to the atmosphere and ensuing global warming effects. PMOBs consist of two main layers: methane oxidation layer (MOL) and gas distribution layer (GDL). The oxidation of methane occurs in MOL through the biochemical reactions of methanotrophic bacteria, and GDL is constructed beneath the MOL to intercept and distribute the fugitive biogas emissions at the base of MOL. Basically, the efficiency of a PMOB is defined based on the methane oxidation efficiency in MOL. Therefore, it is indispensable to provide adequate conditions for the bacterial activities of methanotrophs. In addition to the environmental parameters, the intensity and the distribution of the biogas reaching the MOL material influence the efficiency of PMOBs, and they may cause the MOL material possessing great capacity to host the bacterial activities to be unserviceable in terms of in field methane oxidation. The capillary barrier effect along the GDL MOL interface may provoke localized surface methane emissions, resulted from the restricted and/or non uniform distribution of upward flow of biogas at the base of MOL. The main focus of present study is to incorporate the unsaturated hydraulic behavior of PMOBs into the design of PMOBs, providing adequate ease and distribution of upward flow of biogas at the base of MOL. The air permeability functions of the materials used to construct the MOL of experimental PMOBs at the St Nicephore landfill (Quebec, Canada), along with other materials from the technical literature, were studied to evaluate the unsaturated gas flow behavior of the materials and to identify the threshold of unrestricted gas migration. This latter threshold was introduced as a design parameter based on which the recommended design criterion herein, i.e. the length of unrestricted gas migration (LUGM), was defined. LUGM is considered as the length along the GDL MOL interface along which biogas can migrate upwards without restriction. Performing sets of numerical simulations in SEEP/W, the effect of slope of interface and the parameters defining the water retention curve and hydraulic conductivity function of MOL material on value of LUGM (representing the ease of upward flow of biogas at the interface) and distribution of moisture (and therefore biogas) along the GDL MOL interface were assessed. The saturated hydraulic conductivity and the pore size distribution of the MOL material were the most influencing parameters in distribution of moisture along the interface. The latter parameter influences also the value of degree of saturation and therefore, the ease of biogas at the base of MOL. Dry density of MOL material is another parameter that controls the ease of upward flow of biogas. The main limitations of the present study are associated with the number of tested MOL materials and the inability of SEEP/W in considering the evapotranspiration. However, considering reasonable assumptions in simulations and using the data from the literature, it was attempted to reduce the limitations. Based on the results of experiments and numerical simulations, some design steps and considerations for selection of the MOL material and the slope of interface were suggested that incorporate the unsaturated hydraulic behavior into the design necessities for an efficient PMOB so that the maximum possible methane oxidation capacity of MOL material is exploited.

On désigne habituellement par retrait-gonflement des sols argileux les variations de volume qu’un massif d’argile peut subir sous l’effet des variations de la pression interstitielle négative (succion) de l’eau, lors d’une humidification ou d’un séchage. L’objectif de ce travail est de compléter les connaissances expérimentales sur le comportement hydromécanique des sols gonflants rencontrés dans les travaux de génie civil, de géotechnique et de géotechnique de l’environnement afin de mieux en tenir compte dans les modèles de calcul. De nombreux travaux ont montré l'intérêt de réaliser des essais œdométriques pour caractériser le potentiel et la cinétique de gonflement des sols. Ils ont révélé également l’influence de l’humidification ou du séchage sur le gonflement ou le retrait, avec des pertes ou des augmentations de rigidité des éprouvettes. En revanche, l’influence de plusieurs cycles de séchage/humidification sur le comportement mécanique n’a été que très rarement analysée. Cette thèse rapporte les résultats d’études expérimentales effectuées à l’œdomètre avec imposition de succion par la méthode osmotique sur des matériaux gonflants compactés lâches et naturels denses. Plusieurs cycles de séchage/humidification ont été appliqués sur ces matériaux sous trois faibles charges mécaniques constantes. Pendant ces cycles, les éprouvettes manifestent un retrait cumulé pour le sol lâche et un gonflement cumulé pour le sol dense. Les résultats montrent que dans les deux cas, les déformations volumiques convergent vers un état d’équilibre où le sol présente un comportement réversible. A la fin des cycles de succion, un cycle de chargement/déchargement a été effectué sous les succions constantes. Les valeurs de la pression de préconsolidation p0(s), de l’indice de compression vierge [lambda(s]] et de l’indice de compression élastique [kappa] dépendent directement des chemins de contrainte suivis. L’ensemble des résultats expérimentaux permet de déterminer les surfaces de charge : la limite de séparation de micro/macro (Lm/M) ; la surface de chargement-effondrement (LC : Loading Collapse du modèle BBM) et la surface de comportement saturé (SCS). La succion limite entre la micro- et la macrostructure (Lm/M) dépend parfaitement de la structure interne et du diamètre qui délimite les deux familles de pores. L’évolution de la pression de préconsolidation en fonction de la succion imposée est présentée par la surface LC. Les courbes de compressibilité sous différentes succions convergent vers la courbe correspondant à l’état saturé sous de fortes contraintes appliquées. La pression à partir de laquelle, le sol continue son chemin sur la courbe du comportement normalement consolidé est appelée la pression de saturation (Psat). Plus la succion imposée est élevée, plus la charge nécessaire pour atteindre cette pression de saturation est importante. La surface SCS présente la variation de la pression de saturation en fonction de la succion imposée. Nous pouvons considérer que les surfaces de charge SCS et LC sont uniques pour les sols denses cependant elles se superposent à la fin des cycles de succion pour les sols lâches. Les cycles hydriques augmentent aussi la limite (Lm/M) entre la micro- et la macrostructure pour les deux sols / We usually define the swelling-shrinkage of the swelling soils by the volume variation of a clayey layer exposed to the negative water pore pressure (suction) variations during the wetting and drying periods. In this research, we try to complete our experimental knowledge on the hydromechanical behaviour of the swelling soils used in civil engineering, geotechnical engineering and geoenvironmental engineering, to better modelize these soils in our numerical calculations. Several authors have used the oedometer tests to characterize the capacity and the process of the soil expansion. They have also studied the influence of the wetting and drying on the soil swelling or shrinkage which can decrease or increase the soil rigidity. However, the influence of several hydraulic cycles on the mechanical behaviour has been rarely studied. This thesis presents an experimental study performed on compacted loose and natural dense expansive soils using osmotic oedometers. Several successive cycles were applied under three different low constant vertical net stresses. The loose soil presents a significant shrinkage accumulation while the dense one produces the swelling accumulation during the suction cycles. The suction cycles induced an equilibrium stage which indicates an elastic behaviour of the samples. At the end of suction cycles, a loading/unloading test was performed at the constant suctions for both materials. The mechanical parameters, i.e. the virgin compression index [lambda(s)], the apparent preconsolidation stress p0(s) and the elastic compression index values [kappa] ?are completely dependent on the followed stress paths. The whole experimental results made it possible to define the yielding surfaces: suction limit between micro and macrostructure (Lm/M), loading collapse (LC) and saturation curve (SCS). The suction limit (Lm/M) depends completely to the soil fabrics and to the diameter separating the micro- and macrostructure. The preconsolidation stress variation with suction is represented by the LC surface. The compression curves at different imposed suctions converge towards the saturated state for the high applied vertical stresses. We consider the saturation pressure (Psat) as the necessary pressure to reach the saturated state for an imposed suction. The higher the suction, the higher the saturation pressure. The yielding surface representing this pressure as a function of suction is called the saturation curve (SCS). Generally we can state that the suction cycles unified the LC and SC surfaces and increased the (Lm/M) up to a higher value

Ce travail de recherche se situe dans le cadre du développement de la méthode des éléments de frontière (BEM) pour les milieux poreux multiphasiques. À l'heure actuelle, l'application de la BEM aux pr oblèmes des milieux poreux non-saturés est encore limitée, car l'expression analytique exacte de la solution fondamentale n'a pas été obtenue, ni dans le domaine transformé ni dans le domaine réel. Ceci provient de la complexité du système d'équations régissant le comportement des milieux poreux non-saturés. Les développements de la BEM pour les sols non-saturés effectués au cours de cette thèse sont basés sur les modèles thermo-hydro-mécanique (THHM) et hydro-mécanique (HHM) présentés dans la première partie de ce mémoire. Ces modèles phénoménologiques basés sur la théorie de la poromécanique et les acquis expérimentaux sont obtenus dans le cadre du modèle mathématique présenté par Gatmiri (1997) et Gatmiri et al. (1998). Après avoir présenté les modèles THHM et HHM, on établit pour la première fois les équations intégrales de frontière et les solutions fondamentales associées pour un milieu poreux non-saturé sous chargement quasi-statique pour les deux cas isotherme (2D dans le domaine de Laplace) et non-isotherme (2D et 3D dans les domaines de Laplace et temporel). Aussi, les équations intégrales de frontière ainsi que les solutions fondamentales 2D et 3D (dans le domaine de Laplace) pour le modèle dynamique couplé des sols non-saturés sont obtenues. Ensuite, les formulations d'éléments de frontière (BEM) basées sur la méthode quadrature de convolution (MQC) concernant les milieux poreux saturé et non-saturé sous chargement quasi-statique isotherme et dynamique sont implémentées dans le code de calcul « HYBRID ». Ayant intégrées les formulations de BEM pour les problèmes de propagation d'ondes ainsi que pour les problèmes de consolidation dans les milieux poreux saturés et non-saturés, il semble que nous ayons fourni à l'heure actuelle le premier code de calcul aux éléments de frontière (BEM) qui modélise les différents problèmes dans les sols secs, saturés et non-saturés. Une fois le code vérifié et validé, des études paramétriques portant sur des effets de site sismiques sont effectuées. Le but recherché est d'aboutir à un critère simple, directement exploitable par les ingénieurs, combinant les caractéristiques géométriques et les caractéristiques du sol, permettant de prédire l'amplification du spectre de réponse en accélération dans des vallées sédimentaires aussi bien que vides / The purpose of this dissertation is to develop a boundary element method (BEM) for multiphase porous media. Nowadays, the application of the BEM for solving problems of unsaturated porous media is still limited, because no fundamental solution exists in the published literature, neither in the frequency nor time domain. This fact rises from the complexity of the coupled partial differential equations governing the behaviour of such media. The developments of the BEM for the unsaturated soils carried out during this thesis are based on the thermo-hydro-mechanical (THHM) and hydro-mechanical (HHM) models presented in the first part of this dissertation. These phenomenological models are presented based on the experimental observations and with respect to the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) and Gatmiri et al. (1998). After having presented the THHM and HHM models, for the first time, one establishes the boundary integral equations (BIE) and the associated fundamental solutions for the unsaturated porous media subjected to quasi-static loading for both isothermal (2D in the Laplace transform domain) and non-isothermal (2D and 3D in Laplace transform and time domains) cases. Also, the boundary integral equations as well as the fundamental solutions (2D and 3D in the Laplace transform domain) are obtained for the fully coupled dynamic model of unsaturated soils.In the next step, the boundary element formulations (BEM) based on the convolution quadrature method (CQM) regarding the saturated and unsaturated porous media subjected to isothermal quasi-static and dynamic loadings are implemented via the computer code HYBRID. Having integrated the BEM formulations for the wave propagation, as well as the consolidation problems in the saturated and unsaturated porous media, it seems that now the first boundary element code is obtained that can model the various problems in dry, saturated and unsaturated soils. Once the code is verified and validated, parametric studies on seismic site effects are carried out. The aim is to achieve a simple criterion directly usable by engineers, combining the topographical and geological characteristics of the soil, to predict the amplification of acceleration response spectra in sedimentary as well as hollow valleys

Ce travail vise à étudier le couplage entre écoulement et interactions physico-chimiques dans les sols, dans différentes conditions de saturation en eau, afin d’améliorer la prédiction du devenir des polluants. Il s’agit de comprendre en quoi le taux de saturation du milieu affecte la réactivité du sol vis-à-vis des polluants, et d’évaluer le pouvoir prédictif du transport de solutés réactifs étudié en milieu saturé sur la réactivité en conditions non saturées. Différents processus sont considérés : l’échange de cations calcium-zinc sur un milieu poreux modèle (sable-kaolinite), la sorption et désorption d’un composé organique sur une terre non contaminée, le transport de polluants prioritaires tels que les HAP sur une terre de friche industrielle. Dans chaque cas, des expériences en colonne de laboratoire ont été conduites en conditions d’écoulement saturé et non saturé permanent, permettant tout d’abord la caractérisation de l’hydrodynamique, puis l’étude du couplage avec la réactivité. Les courbes de percée obtenues ont été ensuite modélisées avec des codes tels que CXTFIT. On a montré l’influence de la teneur en eau du milieu sur le transport réactif, variable suivant le type de réaction considéré, la structure des milieux jouant également un rôle important. L’échange d’ions sur le milieu modèle n’est globalement pas affecté par la teneur en eau, dans une gamme proche de la saturation. En revanche, une plus forte sorption et une plus faible mobilisation des polluants organiques ont été observées en conditions non saturées. Le transport réactif de ces composés ne peut donc pas être prédit en conditions non saturées à partir de mesures en milieu saturé, qui peuvent surestimer le transport / The aim of this work was to study the link between water flow and physical and chemical interactions in soils under variably water flow conditions, in order to improve the prediction of contaminants fate. It deals with understanding how the porous media water content can modify soil reactivity towards contaminants, and assessing the possibility to predict reactivity under unsaturated conditions with reactive solute transport studied in saturated porous media. Various processes were considered: cations exchange calcium-zinc on a model porous media (sand-kaolinite), sorption and desorption of an organic compound on a non polluted soil, transport of priority contaminants such as PAHs on an industrial contaminated soil. In each case, experiments were carried out with soil columns at the laboratory scale under saturated and unsaturated steady-state flow conditions, in order to characterize at first hydrodynamics and then to study the link with reactivity. Modeling of the breakthrough curves was then performed with codes such as CXTFIT. We showed an influence of porous media water content on reactive transport which was different as a function of the interaction. Porous media structure must also be taken into account. Ions exchange on a model porous media was not globally modified by the water content varying in a range close to saturation. On the contrary, higher sorption and lower migration of organic contaminants were observed under unsaturated conditions. Reactive transport of these compounds cannot therefore be predicted under unsaturated conditions with tests performed on saturated porous media which may overestimate transport

Les propriétés des argilites du Callovo-Oxfordien les désignent comme une barrière naturelle sûre pour le stockage de déchets radioactifs. Afin d’optimiser la prédiction de leur comportement, leur réponse à diverses sollicitations à court et long termes est étudiée. Notamment, lors de la phase d’exploitation, l’air ventilé dans les galeries n’est pas à l’équilibre hydrique avec la roche et peut perturber ses propriétés hydromécaniques (HM). Il semble alors essentiel de caractériser la réponse HM de la roche à des sollicitations hydriques.La démarche adoptée consiste à proposer un modèle physique adéquat, sur la base d’essais au laboratoire. L’essai de séchage est retenu puisqu’il permet d’étudier la réponse des échantillons en conditions non-saturées. A partir des résultats HM, un modèle élastique 2D isotrope est proposé, puis est élevé à la 3D avec un module de Young isotrope transverse. Ensuite, les données expérimentales servent à estimer certains paramètres poroélastiques et de transfert du modèle. Cette étape est accomplie par procédure inverse (minimisation de l’erreur mesures-calculs). Enfin, une modélisation 2D du comportement in situ est proposée, et compare les prédictions de modèles plastique isotrope et élastique isotrope transverse.Les simulations d’essais au laboratoire reproduisent assez bien les données expérimentales. Pour la modélisation in situ, une bonne corrélation est obtenue entre les prédictions et les mesures, et ce sans ajustement préalable des paramètres. Toutefois, le modèle souligne une influence limitée de la plasticité à l’échelle du laboratoire, alors que les phénomènes dissipatifs sont marqués in situ. Les modélisations 3D au laboratoire ne donnent pas de résultats plus fins qu’en 2D, mais reproduisent plus de données expérimentales (variations de masse, déformations axiales et latérales). De plus, injecter plusieurs types de données dans la formulation inverse permet d’améliorer la précision de l’algorithme. Par ailleurs, une meilleure stabilité de l’algorithme est obtenue en adoptant une convergence en deux étapes (simplex, puis méthode de type gradient). Les estimations numériques des paramètres corroborent les mesures expérimentales obtenues par ailleurs. / The Callovo-Oxfordian claystones’ properties make them reliable as a geological barrier for the confinement of radioactive wastes. In order to optimally predict their behavior, how they respond to various short and long terms loadings has to be studied. Particularly during the exploitation phase, air is continuously ventilated throughout the galleries. The climatic properties of this air are not balanced with those of the rock, and may perturb its hydromechanical (HM) attributes. Thus, assessing the HM response of the rock under hydric loading seems to be a priority.This dissertation begins with laboratory tests to propose an appropriate physical model. Drying tests were studied as they focus on the HM response of samples undergoing hydric loadings. A first 2D isotropic model is proposed, and then enhanced to 3D by considering a transversely isotropic Young modulus. Secondly, experimental results provide relevant data to estimate poroelastic and transport parameters involved in the model. Estimation is achieved according to an inverse procedure, which minimizes the error between measurements and model predictions. Finally, a real-size test is simulated using 2D models: an isotropic plastic one and a transversely isotropic elastic one.Model predictions reproduce well the laboratory tests data. When simulating the in situ behavior, a rather good agreement is obtained between the numerical and experimental results (although using the parameters estimated at the laboratory scale). However, the model highlights a limited influence of plasticity in the laboratory tests, while dissipative phenomena obviously occur in situ. 3D laboratory simulations do not improve the precision of 2D results, but reproduce more experimental data (mass variations, axial and lateral strains). Moreover, the inversion process is more efficient when ran over various kinds of data. Furthermore, stability of the algorithm is improved when adopting a two-phase convergence (simplex, followed by a gradient-like method). Numerical estimates of the parameters are in agreement with the direct experimental measurements obtained through other tests.