Soutenance mémoire
Notions de base sur les écoulements
Méthodes graphiques
La méthode graphique est appelée aussi méthode des petits carreaux. Elle est le procédée le plus usuel pour résoudre les problèmes d’écoulements en milieu poreux.
Elle consiste à tracer le réseau orthogonal des lignes de courant et des lignes équipotentielles en tenant compte des conditions aux limites de l’écoulement.
La méthode graphique nécessite une expérience de la part de l’utilisateur et des approximations successives du tracé de réseau. Les résultats obtenus par cette méthode donnent un ordre de grandeur des valeurs recherchées, et permettront ainsi de contrôler les résultats d’autres méthodes.
Le procédé est simple et constitué des étapes suivantes :
a) Définir les conditions aux limites concernant l’écoulement, c’est-à-dire les lignes de courants et les lignes équipotentielles ;
b) Tracer quelques lignes de courant en essayant de satisfaire les conditions aux limites. Le tracé est à refaire jusqu’à la satisfaction de ces conditions. Il est commode que les lignes tracées soient équidistantes ;
c) Tracer des lignes équipotentielles de telle façon qu’elles soient perpendiculaires aux lignes de courant et qu’elles vérifient les conditions aux limites. Les lignes sont choisies équidistantes et ayant un pas identique à celui des lignes de courant.
Méthodes numériques
Elles sont largement utilisées dans les études des écoulements à travers et sous les ouvrages en surface libre ou confinée. Elles sont basées sur le calcul numérique. Les méthodes les plus usuelles sont : la méthode des différences finies et la méthode des éléments finis.
Méthode des différences finies
Elle consiste à remplacer les grandeurs différentielles infiniment petites , , , ℎ, , ,., par des grandeurs de dimensions finies que l’on note x,y,z,h,v,t, petites vis à vis du volume à étudier. La précision des résultats du calcul est d’autant meilleure que les i
sont plus petits. Il y a intérêt à choisir des éléments d’aussi grande dimension que possible pour en limiter le nombre. Mais, si les éléments sont trop grands, le calcul diverge. Il n’y a plus de solution. De telles méthodes doivent donc comporter un critère de convergence qui limite les dimensions des éléments.
En différences finies, l’équation de Laplace :
C’est l’équation de la diffusion en régime permanent pour un écoulement dans un plan vertical d’un milieu homogène et isotrope.
La fonction hx,z, la charge hydraulique (figure 1.6), prend les valeurs h0 pour le point x0, h1pour le point x x x 1 0 et h2 pour le point x x x 2 0.Ce qui permet de trouver le champ de charge, donc les équipotentielles (Humbert, 1984). Le même calcul peut être fait pour la fonction de courant et l’on trouve les lignes de courant.
CONCLUSION
L’eau joue un rôle important dans les ouvrages géotechniques son état défini le champ d’application les potentiels à prendre en compte dans les équations régissant le phénomène de l’écoulement dans le sol.
La loi de Darcy avec l’expression de continuité et l’état isotherme de sol donne l’équation qui régit le problème étudié dans ce projet de fin d’études.
L’ingénieur dispose d’un ensemble de méthodes et techniques permettant de résoudre la problématique de l’eau dans le sol. Il est confronté donc à un problème de choix. Son choix peut être bon dans une situation et mauvais dans une autre, selon les moyens disponibles sa mise en œuvre et surtout de l’intuition et l’expérience de géotechnicien.
Ce qui définit par la suite le modèle de calcul par simplification de la structure réelle moyennant certaines hypothèses simplificatrices.Le chapitre suivant donne un aperçu sur les modèles disponibles dans la littérature.
INTRODUCTION
Ce chapitre contient un aperçu bibliographique sur des éléments permettant la compréhension des écoulements en non saturés. Après quelques définitions de bases sur le processus d’infiltration, le chapitre comporte les méthodes de mesure de la succion et de la conductivité dans le cas d’un sol non saturé, le phénomène d’hystérésis, la non linéarité physique de la courbe de rétention font l’objet de ce chapitre. Ces éléments sont nécessaires pour la modélisation numérique pour le chapitre 4 de projet de fin d’études.
APERÇU GENERAL SUR LE PROCESSUS D’INFILTRATION
Définitions
L’infiltration peut se définir comme le processus physique d’entrée d’eau à l’intérieur du sol à travers sa surface, in-situ l’infiltration représente le même phénomène à travers les parois d’un puits peu profond. Nous intéressons dans ce travail au processus d’infiltration unidirectionnel sous charge constante, qui sera étudié à travers une surface de sol compacté (pour simuler l’infiltration des lixiviats issue de l’évacuation dans les décharges de déchets ménagers)
Afin de mieux mener cette étude, nous allons commencer par une description générale du processus d’infiltration ainsi que l’état énergétique de l’eau dans un sol non saturé, ensuite, on présentera les équations de base décrivant l’écoulement de l’eau dans un sol non saturé et en particulier dans le cas qui nous intéresse qui est celui de l’infiltration en considérant l’approche monophasique de l’écoulement. Enfin, on donne quelques approches physico-empiriques d’estimation de l’infiltration souvent exprimées sous forme de loi d’infiltration.
Paramètres descriptifs de l’infiltration
Le processus d’infiltration est souvent engendré par une modification brutale des conditions de pression et de teneur en eau régnant à la surface du sol lorsque celui-ci reçoit une averse ou s’il est soumis à une submersion. Il est conditionné par divers facteurs dont les plus importants relèvent d’une part, du sol par le biais de ses caractéristiques hydrodynamiques ((), D ()), de sa structure, de sa texture et de son uniformité de profil, et d’autre, part des conditions spécifiques dans lesquelles se déroulent le processus à savoir les conditions initiales et le débit d’alimentation (Musy et Soutter, 1991).
En surface, le processus d’infiltration est caractérisé par le flux traversant le sol et s’écoulant dans le profil appelé le régime d’infiltration i(t). Cette grandeur permet de déterminer la lame d’eau infiltrée I(t) donnée par l’intégrale dans le temps de i(t).
On définit également la capacité d’infiltration ou l’infiltrabilité du sol comme le flux maximum que le sol est en mesure d’absorber à travers sa surface lorsque celle-ci est maintenue en contact avec de l’eau à pression atmosphérique (Hillel, 1988 ;Vauclin, 1984; Musy et Soutter, 1991). En général, l’infiltrabilité du sol est élevée au début de l’infiltration (en particulier pour les sols initialement secs) et décroît rapidement durant les premiers instants jusqu’à ce qu’elle atteigne une valeur stable appelée infiltration permanente ( ). Cette décroissance est due avant tout à la diminution du gradient de succion matriciel, conséquence de l’augmentation de la distance sur laquelle se répartie la différence globale de la charge hydraulique. Dans certains cas, cette diminution peut être renforcée par des facteurs secondaires qui provoquent l’augmentation de la résistance à l’écoulement et qui peuvent être résumés de la façon suivante :
Dégradation graduelle de la structure du sol qui entraîne, par migration des particules, un colmatage partiel du système poreux et par la suite, la formation d’une croûte dense en surface.
Foisonnement des argiles.
Piégeage des bulles d’air ou compression de la masse d’air si elle ne peut s’échapper quand elle est déplacée par l’arrivée d’eau.
Si la surface du sol est alimentée en eau avec une vitesse croissante.
l’alimentation dépassera la capacité d’adsorption en eau du sol. Le processus d’infiltration peur être résumé de la manière suivante (Vauclin.1984) : Aussi longtemps que le débit d’alimentation en surface () est inférieur au régime d’infiltration i(t),
l’eau s’infiltre dans le sol aussi vite qu’elle est fournie (la vitesse d’infiltration () = ()). Dans ce cas, c’est le régime d’alimentation qui détermine le régime d’infiltration et le processus se déroule à flux contrôlé. Quand () dépasse l’infiltrabilité du sol, cette dernière détermine la vitesse d’infiltration de l’eau dans le sol (() = ()) et par conséquent, le processus est contrôlé par le profil du sol. La surface du sol devient saturée et l’excès d’eau ( () − ()) stagne ou ruisselle en présence de pente. Quand le débit d’alimentation ()est constante ( () =
) et inférieur à l’infiltrabilité du sol, il s’établit au sein du sol un profil d’humidité q0 uniforme pour lequel la conductivité hydraulique ( ) est égale . La figure 2.1 résume le processus d’infiltration à charge constante et à flux constant représente la conductivité hydraulique à la saturation du sol.
Distribution de l’eau dans le sol durant l’infiltration
La répartition de l’eau dans le sol à un instant donné peut être décrite par le profil hydrique. La forme de celui-ci est liée d’une part aux propriétés physiques du sol et d’autre part aux circonstances météorologiques très variables. L’infiltration de l’eau en surface provoque une perturbation du profil hydrique qui se propage en profondeur à mesure que le processus se poursuit. Cette perturbation lui confère une structure qui est caractérisée selon Musy et Soutter (1991) par quatre zones (figure 2.2). Dans le cas d’un profil homogène soumis à une alimentation constante en surface, ces zones sont arrangées comme suit :
La zone de saturation en surface qui s’étend généralement à une profondeur de 1 cm ; dans le cas des sols argileux, cette zone peut être soumise à un phénomène de gonflement qui conduit à des teneurs en eau très élevées ;
La zone de transmission, presque saturée, caractérisée par une diminution très rapide de la teneur en eau jusqu’à des valeurs de 70 à 80% de la saturation ; cette zone s’allonge au cours de l’infiltration ;
La zone d’humidification dans laquelle le teneur en eau change rapidement avec la profondeur et le temps, elle se termine avec le front d’humidification ;
Le front d’humidification est une région d’un gradient d’humidité très élevé qui représente une limite visible de pénétration d’humidité surtout quand le sol est initialement sec.
Pendant l’infiltration, l’eau entre dans le sol et descend progressivement vers le bas en établissant des profils de teneur en eau successifs (figure 2.3). La distribution de l’eau dans le sol dépend d’une manière générale de la taille des grains. L’examen périodique des profils d’humidité pendant l’infiltration montre que la zone de transmission s’allonge continuellement, et, que la zone d’humidification et le front se déplacent vers le bas. Ce dernier devient moins net en profondeur. Le gradient de succion, initialement élevé notamment pour un sol sec, diminue au cours de l’infiltration jusqu’à des valeurs unitaires (figure 2.3). Ceci explique la diminution de la vitesse d’infiltration avec le temps.
Influence de certains paramètres du sol sur l’infiltration
En général, les paramètres qui influent sur processus d’infiltration sont nombreux. Ils sont liés soit aux propriétés du sol soit aux conditions initiales et aux limites dans lesquelles l’infiltration se déroule. Les paramètres liés aux propriétés du sol sont : la taille des grains, le type et le degré de développement de la structure, la teneur en matière organique, la porosité, la masse volumique et l’uniformité du profil (Hillel, 1988; Vauclin, 1984). La figure 2.4 montre un exemple de l’influence de la masse volumique sèche du sol sur l’avancement du front d’humidité au cours d’une infiltration. Si on considère que le déplacement du front peut être caractérisé par une teneur en eau q de 10%, on constate que ce front avance de 5 fois plus vite pour le sol ayant la densité la plus faible. Les conditions initiales d’humidité et de succion et les conditions aux limites (charge d’eau appliquée, milieu d’écoulement ouvert ou fermé à la base…) peuvent être contrôlées et étudiées au laboratoire. Ainsi l’influence de ces paramètres a été établie par certains chercheurs (Philip, 1957 et Vauclin, 1986).
L’humidité initiale du sol qi a une grande influence sur la vitesse d’infiltration de l’eau dans le sol surtout au début de l’infiltration. La diminution de qi entraîne une augmentation de la vitesse d’infiltration initiale mais son effet diminue pour des temps d’infiltration élevés (figure 2.5). En effet, pour des faibles qi, la succion correspondante dans le sol est élevée et par conséquent, l’écoulement dû à la capillarité est important entraînant ainsi des vitesses d’écoulement considérables au début de l’infiltration. Au fur et à mesure que le sol se sature, l’effet de la succion sur l’écoulement devient faible par rapport à celui dû à la gravité ; l’influence de la teneur en eau initiale ne se fait plus sentir. Par contre, quand le sol est initialement humide, le gradient de succion est faible au niveau du front d’humidification et réduit donc la vitesse d’entrée d’eau. La vitesse
d’avancement du front peut être rapide du fait du déplacement vers le bas de l’eau déjà existante dans le sol.
L’influence de la charge d’eau en surface sur la dynamique d’infiltration a été également étudiée par Philip (1957 f). La figure 2.6 montre cette influence aussi bien sur l’infiltration cumulée que sur la vitesse d’infiltration. L’augmentation de la charge d’eau entraîne une augmentation de la qualité d’eau infiltrée dans le sol et par conséquent de la vitesse d’infiltration. Ce phénomène est surtout important au début mais devient moins marqué pour des temps élevés. L’augmentation de la charge d’eau entraîne une augmentation de la sorptivité du sol défini par Philip (1957) comme étant la capacité du sol à absorber l’eau par capillarité. Ainsi, la sorptivité du sol augmente typiquement de 2% par centimètre de charge d’eau.
L’influence du confinement sur la lame d’eau infiltrée dans le sol et le flux correspondant a été étudiée par Vauclin (1984) et Wang et Feyen (1995) cité par Amraoui (1996). La figure 2.7 montre que ces deux paramètres d’infiltration diminuent quand l’air initialement présent dans la zone non saturée n’est pas libre de circuler. En effet, dans ce type de milieu et au cours de l’infiltration, l’air ne peut pas s’échapper ; sa pression augmente jusqu’à une valeur critique qui lui permet de se libérer vers le haut au profit d’une avancée d’eau vers le bas. Ce phénomène ralentit la vitesse d’avancement du front d’infiltration. L’air est surtout piégé dans des pores isolés et l’infiltration se réduit par diminution de la porosité effective.
|
Table des matières
INTRODUCTION GENERALE
Chapitre I : Notions de base sur les écoulements
I.1 : Introduction
I.2 : L’eau dans le sol
I.2.1 : L’eau libre
I.2.2 : L’eau capillaire
I.2.3 : L’eau de constitution
I.2.4 : L’eau liée ou adsorbée (Eau de rétention)
I.2.5 : Etude du mouvement d’eau dans le sol
I.2.5.1 : Théorie de la consolidation
I.2.5.2 : Théorie hydrique du sol
I.2.6 : Charge hydraulique
I.2.7 : Gradient hydraulique
I.2.8 : Vitesse d’écoulement
I.2.9 : Lignes de courant
I.3 : Loi de Darcy
I.3.1 : Perméabilité
I.3.1.1 : Le coefficient de perméabilité
I.3.1.2 : La perméabilité intrinsèque
I.3.2 : Limite de la loi de Darcy
I.3.3 : Généralisation de la loi de Darcy
I.3.3.1 : Equations des écoulements en régime permanent
I.3.3.2 : Cas particuliers
1.4 : Généralités sur la résolution de ces équations
I.4.1 : Conditions aux limites des écoulements souterrains
I.4.1.1 : Surface imperméable (condition de Newman)
I.4.1.2 : Surfaces équipotentielles (condition de Dirichlet)
I.4.1.3 : Surface de suintement
I.4.1.4 : Surface libre (condition mixte)
I.4.2 : Méthodes de résolution de ces équations
I.4.2.1 : Méthodes analytiques
I.4.2.2 : Méthodes de similitudes
I.4.2.3 : Modèles réduits
I.4.2.4 : Modèles analogiques
I.4.2.5 : Méthodes graphiques
I.4.2.6 : Méthodes numériques
I.5 : Conclusion
Chapitre II: Processus d’infiltration dans les sols non saturés
II.1 : Introduction
II.2 : Aperçu général sur le processus
II.2.1 : Définitions
II.2.2 : Paramètres descriptifs de l’infiltration
II.2.3 : Distribution de l’eau dans le sol durant l’infiltration
II.2.4 : Influence de certains paramètres du sol sur l’infiltration
II.3 : Succion dans le sol et relation avec les paramètres hydrodynamiques
II.3.1 : Succion dans le sol
II.3.2 : Dispositifs de mesure de la succion dans le sol
II.3.3 : Relation teneur en eau – succion ; courbe de rétention d’un sol
II.3.4 : Relation entre la conductivité hydraulique et la succion du sol
II.4 : Conclusion
Chapitre III: Proposition d’une solution analytique
III.1 : Introduction
III.2 : Mouvement horizontal
III.2.1 : Méthode de Boltzmann
III.2.2 : Méthode de séparation de variables
III.2.3 : Méthode exp-function
III.3 : Mouvement vertical
III.4 : Mouvement bidimensionnel
III.4.1 : Revue des travaux antérieurs
III.4.2 : Positionnement de problème à résoudre
III.4.3 : Formulation mathématique et méthodologie
III.4.3.1 : Cas d’une variation exponentielle de paramètres K(u) et C(u)
III.4.3.2 : Cas d’une variation exponentielle en U (exp-function)
III.4.4 : Résultats et interprétation
III.5 : Conclusion
Chapitre IV: Analyses numériques des écoulements en régime transitoire
IV 1 : Introduction
IV.2 : Description globale de la méthodologie
IV.3 : Cas d’une colonne verticale (imbibition 1D) (BENYELLES, 1988)
IV.3.1 : Introduction
IV.3.2 : Conditions aux limites et propriétés des matériaux
IV.3.3 : Paramètres de convergence
IV.3.4 : Résultats et discussion
IV.3.5 : Influence de la perméabilité à saturation sur le phénomène d’imbibition
IV.3.6 : Influence de la teneur en eau à saturation sur l’évolution de phénomène
IV.4 : Réponse hydraulique d’un barrage en terre
IV.4.1 : Géométrie et conditions aux limites
IV.4.1.1 : Phase de réservoir vide
IV.4.1.2 : Phase de remplissage
IV.4.1.3 : Phase d’exploitation de réservoir plein
IV.4.1.4 : Phase de vidange rapide
IV.4.2 : Propriétés de matériau
IV.4.3 : Résultats et discussions
IV.5 : Conclusion
CONCLUSION GENERALE
REFERENCES BIBLIOGRAPHIQUES
Télécharger le rapport complet
