Dans cette thèse, on étudie les écoulements oscillatoires en milieux poreux (non saturés ou partiellement saturés) dus à des oscillations tidales des niveaux d'eau dans des milieux ouverts adjacents aux milieux poreux. L'étude est centrée sur le cas des plages de sable en hydrodynamique côtière, mais les applications concernent, potentiellement et plus généralement, les problèmes d'oscillation et de variation temporelle des niveaux d'eau dans des systèmes couplés, lorsque ceux-ci mettent en jeu des interactions entre les écoulements de sub-surface (milieux poreux) et les eaux de surface (milieux ouverts) : plages naturelles et artificielles; digues portuaires; barrages en terre; berges de fleuves; estuaires. Le forçage tidal des écoulements souterrains est représenté et modélisé ici, tant expérimentalement que numériquement, par une oscillation quasi-statique du niveau d'eau dans un réservoir externe ouvert, connecté au domaine poreux. On s'intéresse plus particulièrement aux écoulements verticaux forcés par une pression oscillatoire imposée au bas d'une colonne de sol. Sur le plan expérimental, ce type de forçage est obtenu par une machine à marée équipée d'un arbre rotatif. Au total, on utilise dans ce travail trois types d'approches (expérimentale, numérique, analytique), l'objectif étant d'étudier le mouvement vertical de la surface "libre" et l'écoulement non saturé sus-jacent, de façon à prendre en compte aussi bien les pertes de charge dans la zone saturée que les gradients de pression capillaire dans la zone non saturée. […] / In this thesis, we study hydrodynamic oscillations in porous bodies (unsaturated or partially saturated), due to tidal oscillations of water levels in adjacent open water bodies. The focus is on beach hydrodynamics, but potential applications concern, more generally, time varying and oscillating water levels in coupled systems involving subsurface / open water interactions (natural and artificial beaches, harbor dykes, earth dams, river banks, estuaries). The tidal forcing of groundwater is represented and modeled (both experimentally and numerically) by quasi-static oscillations of water levels in an open water reservoir connected to the porous medium. Specifically, we focus on vertical water movements forced by an oscillating pressure imposed at the bottom of a soil column. Experimentally, a rotating tide machine is used to achieve this forcing. Overall, we use three types of methods (experimental, numerical, analytical) to study the vertical motion of the groundwater table and the unsaturated flow above it, taking into account the vertical head drop in the saturated zone as well as capillary pressure gradients in the unsaturated zone. Laboratory experiments are conducted on vertical sand columns, with a tide machine to force water table oscillations, and with porous cup tensiometers to measure both positive pressures and suctions along the column (among other measurement methods). Numerical simulations of oscillatory water flow are implemented with the BIGFLOW 3D code (implicit finite volumes, with conjugate gradients for the matrix solver and modified Picard iterations for the nonlinear problem). In addition, an automatic calibration based on a genetic optimization algorithm is implemented for a given tidal frequency, to obtain the hydrodynamic parameters of the experimental soil. Calibrated simulations are then compared to experimental results for other non calibrated frequencies. Finally, a family of quasi-analytical multi-front solutions is developed for the tidal oscillation problem, as an extension of the Green-Ampt piston flow approximation, leading to nonlinear, non-autonomous systems of Ordinary Differential Equations with initial conditions (dynamical systems). The multi-front solutions are tested by comparing them with a refined finite volume solution of the Richards equation. Multi-front solutions are at least 100 times faster, and the match is quite good even for a loamy soil with strong capillary effects (the number of fronts required is small, no more than N≈ to 20 at most). A large set of multi-front simulations is then produced in order to analyze water table and flux fluctuations for a broad range of forcing frequencies. The results, analyzed in terms of means and amplitudes of hydrodynamic variables, indicate the existence, for each soil, of a characteristic frequency separating low frequency / high frequency flow regimes in the porous system.

L'objectif de cette recherche est l'étude du transfert d'eau en phase liquide dans les parois du bâtiment et de définir une nouvelle métrologie pour déterminer le coefficient de perméabilité liquide des matériaux du bâtiment. Une partie théorique fait l'état des connaissances en matière de modélisation du transfert d'eau en phase liquide. Elle est suivie par une partie qui relate les méthodes utilisées pour mesurer ou imposer des pressions capillaires (ou succions) avec un intérêt particulier à la méthode osmotique qui est utilisée plus tard pour imposer des conditions aux limites en terme de succion. Cette dernière partie s'ouvre sur une revue des méthodes de mesure de la perméabilité liquide et sera conclue par la définition et la mise au point d'une nouvelle méthode de mesure de la perméabilité liquide, qui a fait l'objet d'un dépôt de brevet d'invention, s'appuyant sur la méthode osmotique et qui fonctionne en régime permanent, mesurant les débits entrant et sortant de l'échantillon. Une campagne d'essais est menée sur du béton cellulaire autoclave avec différents couples de conditions aux limites. Il s'ensuit une recherche d'identification d'une relation entre la perméabilité liquide et la succion, démarche qui mérite plus d'approfondissement. Finalement cette nouvelle méthode de mesure mérite d'être diffusée, et son caractère opérationnel est manifeste.

Cette thèse concerne la caractérisation expérimentale d’un sable limoneux provenant du barrage de Livet – Gavet (38) dans le cadre du projet ANR TerreDurable avec plusieurs objectifs : 1- Caractériser au travers d’essais de laboratoire le comportement hydromécanique d’un sable fin limoneux (sol A1 dans la classification GTR) en fonction de son état de saturation. Lors de cette étude, un accent particulier est porté sur la caractérisation de ce comportement dans le domaine proche de la saturation. 2- Interpréter le comportement hydrique du matériau sur chemin de drainage – imbibition en relation avec l’analyse de sa microstructure. 3- Fournir d’un point de vue général une base de données et d’analyses exhaustive permettant le développement et la calibration de modèles de comportement des sols fins proches de la saturation, en particulier en considérant des chemins de chargement hydromécanique complexes. Pour l’ensemble de cette étude, le matériau est considéré sous deux états : soit à l’état de pâte (matériau normalement consolidé) préparée à une teneur en eau proche de la limite de liquidité, soit sous forme compactée (matériau sur-consolidé) à différentes énergies de compactage et différentes teneurs en eau initiales. / This thesis concerns the experimental characterization of a silty sand from the Livet - Gavet dam (38) as part of the ANR TerreDurable project, for following objectives: 1- Through laboratory tests, characterize the hydro-mechanical behaviour of a fine silty sand (Type A1 in the GTR classification) according to its saturation state. In this study, particular emphasis is placed on the characterization of this behaviour in the near-saturated domain. 2- Interpret the water behaviour of material on the drainage - imbibition cycles, in relation to the analysis of its microstructure. 3- From a general point of view, provide a comprehensive database and analysis allowing the development and calibration of models of near-saturated fine soil's behaviour, in particular, by considering complex hydro-mechanical loading paths. For all tests in this study, the material is considered in two states: either in the state of paste (normally consolidated material) prepared at water content close to the limit of liquidity, or in compacted state (over consolidated material) at different compaction energies and different initial water contents.

Les routes sont normalement pourvues de systèmes de drainage dimensionnés et exécutés conformément aux règles de l’art en vue d’évacuer, le plus rapidement possible de l’emprise de la route, les eaux des précipitations à considérer sur la base des périodes de retour prises en compte. Toutefois, il subsiste souvent des écoulements d’eau indésirables au niveau des talus et parfois dans les accotements et/ou les surfaces de roulement non imperméabilisées. Une succession de pluies entraîne des infiltrations d’eau qui varient notamment en fonction des conditions climatiques et suivant la texture et la structure du sol. De telles infiltrations ont pour conséquence la réduction des coefficients de sécurité des talus.<p><p>Bien qu’il existe plusieurs publications scientifiques traitant de pluies qui ont conduit à des glissements de terrain (Lim et al. 1996 ;Cho et al. 2001 ;Kim et al. 2004 ;Xue et al. 2007 ;Gavin et al. 2008), les incidences des infiltrations résultant de pluies successives sur le comportement des couches superficielles des sols non saturés ne sont généralement pas prises en compte. Les modèles permettant le calcul de la stabilité de talus des massifs de sols non saturés exigent beaucoup de paramètres parfois difficiles à évaluer et se rapportent ordinairement aux cas d’instabilité provoqués par une remontée du niveau piézométrique des nappes phréatiques.<p><p>Sur la base des essais réalisés en laboratoire, une méthodologie adaptée permettant d’évaluer la variation spatio-temporelle de la teneur en eau du sol en fonction d’une suite de pluies a été développée. Cette méthodologie facilite la prise en compte des effets cumulés des taux d’infiltration associés aux évènements pluvieux et permet d’en déduire le profil de succion ainsi que celui de la cohésion apparente à utiliser en vue de calculer, pour une inclinaison β du talus, l’intervalle de variation du coefficient de sécurité Fs. La méthodologie développée présente un intérêt particulier dans le cas de budgets et infrastructures limités.<p><p>ABSTRACT<p><p>Roads are normally equipped with drainage systems sized and implemented in accordance with the rules of art to evacuate as quickly as possible to the right of way, water precipitation to be considered on the basis of return periods taken into account. However, there are often water flows at the slope side and sometimes in the shoulders and / or running surfaces that are not waterproof. A succession of rain causes a certain amount of water infiltration, which varies according to climatic conditions and depending on the soil texture and structure. Such infiltrations have resulted in reduced safety factor of slopes.<p><p>Although there are several scientific publications on rainfall leading to landslides (Lim et al. 1996; Cho et al. 2001, Kim et al. 2004; Xue et al. 2007; Gavin and al. 2008), impacts resulting from infiltration of successive rains on the behavior of surface layers of unsaturated soils are usually not taken into account. Models for calculating the slope stability of unsaturated soils require many parameters that can be, in certain circumstances, difficult to assess and refer generally to cases of instability caused by a rise in groundwater level.<p><p>Based on laboratory tests, a suitable methodology for assessing the spatial and temporal variation of soil water content induced by a set of rains has been developed. This methodology facilitates the inclusion of the cumulative effects of the infiltration rates associated with rain events and infers from them the profile of suction and that of the apparent cohesion to be used to calculate, for a slope angle β, the range of variation of the safety factor Fs. This methodology is particularly relevant in the case of limited budgets and infrastructures. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Cette recherche doctorale porte sur la modélisation 3D de la dynamique des teneurs en eau dans le sol et le soussol lorsque les écoulements sont à saturation variable. La modélisation est basée sur une version généralisée de la loi de Darcy-Buckingham et de l’équation de Richards multidimensionnelle. Les recherches présentées dans cette thèse concernent différents volets, présentés ci-dessous, dont le fil conducteur est l’analyse de phénomènes d’écoulements en milieux poreux, contenant de fortes hétérogénéités et/ou perturbés par des fluctuations temporelles. Dans le cadre d’un partenariat de recherche entre l’IMFT et l'IRSN sur la problématique du stockage souterrain de déchets radioactifs, nous avons modélisé en 3D la dynamique du front de désaturation d’une couche argileuse autour d’une galerie souterraine ventilée, à l’aide du code volumes finis BIGFLOW 3D. Ce travail a permis de mettre au point une approche d’immersion pour la modélisation des écoulements en milieux composites. Un autre volet de cette thèse concerne une étude analytique et numérique des profils verticaux de succion et de teneur en eau lors d’une infiltration verticale, non-saturante, dans un sol hétérogène finement stratifié. Des solutions analytiques exactes et approchées sont obtenues en régime permanent, à l’aide de transformations de variables, et sont comparées avec des solutions numériques pour différents degrés d’hétérogénéité. De même, les phénomènes de « barrière capillaire » sont étudiés d’une part analytiquement et d’autre part, par expérimentations numériques transitoires d’infiltration sur des systèmes bicouches, en présence d’une nappe plus ou moins profonde. D’autre part, les écoulements nonsaturés sont étudiés, cette fois, en présence d’un forçage transitoire fortement oscillatoire, dans une colonne de sable fin homogène, sous l’effet de fluctuations périodiques du niveau de la surface libre. Une méthode multi-front mise au point lors d’une thèse précédente à l’IMFT est validée numériquement en montrant que cette méthode simule efficacement la dynamique oscillatoire des flux et des profils de pression avec un nombre limité de « fronts ». De plus, nous avons analysé la phénoménologie capillaire/gravitaire des écoulements oscillatoires dans la colonne grâce à un suivi dynamique du plan de flux nul. Enfin, nous présentons, comme extension aux travaux précédents, une étude préliminaire des phénomènes 3D d’infiltration et de redistribution d’eau, notamment lorsque l’hétérogénéité du sol est de type aléatoire. Les premiers essais d’infiltration permettent de tester, en réplique unique, des méthodes de prises de moyennes spatiales des champs de succion et teneur en eau et évaluer la taille minimum du domaine de calcul 3D permettant d’obtenir des profils verticaux moyens représentatifs de l’infiltration dans une réplique unique du sol aléatoire. / This doctoral research, defended at the Institut de Mécanique des Fluides de Toulouse, is devoted to modeling water pressure and water content in soils and in subsurface geologic formations, in the case of variably saturated flow. One of the main scientific objectives of this work is to analyze the response of unsaturated flow systems, and particularly capillary and gravitational effects, in the presence of material heterogeneities, discontinuities, and/or space-time fluctuations. Modeling is based on a generalized version of Darcy- Buckingham’s law, and of Richards’ flow equation. Overall, the topics developed in this PhD thesis focus on several related aspects of variably saturated water flow in the subsurface. These aspects all occur at once in most applications (drying/wetting,heterogeneity, temporal forcing), but they are "decoupled" here for convenience. A preliminary research (collaborative project IMFT / IRSN) was developed to study the 2D/3D drying process at the porous wall of a deep cylindrical excavation in response to hydrometeorological signals. This project has motivated the design and testing of a novel approach to include cavities in the modeling domain. A detailed study of steady state infiltration was developed for the case of finely stratified soils, with parameters that vary continuously and cyclically with depth. Exact and approximate analytical solutions are calculated based on variable transformation methods and on perturbation type approximations, and they are tested numerically using a finite volume code (BIGFLOW 3D). The sensitivity of suction fluctuations vs. stratification wavelength is investigated, as well as the effect of the degree of heterogeneity, and of water table depth. Capillary barrier effects are studied for the case of unsaturated infiltration in multilayer soil systems characterized by a discontinuity of soil properties at interfaces. Numerical experiments are developed for transient infiltration towards a water table through a two-layer system, the goal being to analyze possible capillary barrier effects under various scenarios. On the other hand, we study numerically the case of a partially saturated / unsaturated soil column submitted to highly variable oscillatory pressure at the bottom of the column: this leads to vertical flow oscillations in the unsaturated zone above the water table influenced by tides (coastal beach sand). We analyze the dynamics of this oscillatory flow, where capillary and gravitational effects compete; for this purpose we use a novel method that tracks the positions of the zero flux plane in the unsaturated column. Finally, we also present, as an extension to the previous studies, a preliminary investigation of multidimensional infiltration/redistribution phenomena, particularly for the case of fully 3D random-type soil heterogeneity. The first numerical experiments of 3D infiltration are undertaken based on the single realization approach to soil heterogeneity, and assuming a uniform distribution of wetting at soil surface.