Deux outils mathématiques pour modéliser l'écoulement et le transport de polluants dans un milieu poreux fracture PDF Download
Are you looking for read ebook online? Search for your book and save it on your Kindle device, PC, phones or tablets. Download Deux outils mathématiques pour modéliser l'écoulement et le transport de polluants dans un milieu poreux fracture PDF full book. Access full book title Deux outils mathématiques pour modéliser l'écoulement et le transport de polluants dans un milieu poreux fracture by Clarisse Alboin. Download full books in PDF and EPUB format.
Author: Clarisse Alboin
Publisher:
ISBN:
Category :
Languages : fr
Pages : 198
Get Book Here
Book Description
Les fractures, qui sont des chenaux de circulation privilégiés, modifient grandement l'écoulement et le transport de polluants dans un milieu. Nous nous intéressons aux fractures colmatées que nous assimilons a des milieux poreux ayant une grande perméabilité. Leur présence génère de fortes hétérogénéités. Notre motivation a été d'étudier des techniques permettant de modéliser les transferts de fluide et de contaminants entre la roche et les fractures. Dans cette thèse, nous proposons deux méthodes numériques adaptées à deux échelles de fracturation. La présence d'un réseau de petites fractures est prise en compte par une modèle double porosité. Ce modèle est obtenu via une homogénéisation du milieu fracture, il est compose des équations standard d'écoulement et de transport dans un milieu poreux auxquelles on a ajoute un terme de couplage qui gère les échanges entre la roche et les fractures. Ce terme de couplage est calcule via la résolution de problèmes auxiliaires qui sont des problèmes de diffusion définis en chaque point de l'espace. Vu le grand nombre de ces problèmes, il est déterminant de les résoudre efficacement. Nous avons montre qu'une résolution analytique est préférable a une discrétisation classique. Pour traiter les fractures plus importantes, nous avons développe une méthode de décomposition de domaine. L’ouverture des fractures étant faible devant les dimensions du site de stockage, nous avons assimile les fractures a des interfaces a travers lesquelles la pression et la concentration sont continues et les flux sont discontinus, contrairement a une décomposition de domaine classique, il s'agit, ici, d'une condition de transmission non locale. La différence des flux provenant des sous-domaines intervient comme un terme source dans les fractures. Les modèles d'écoulement et de transport dans les sous-domaines et dans les fractures ont été obtenus via une analyse asymptotique.
Author: Clarisse Alboin
Publisher:
ISBN:
Category :
Languages : fr
Pages : 198
Get Book Here
Book Description
Les fractures, qui sont des chenaux de circulation privilégiés, modifient grandement l'écoulement et le transport de polluants dans un milieu. Nous nous intéressons aux fractures colmatées que nous assimilons a des milieux poreux ayant une grande perméabilité. Leur présence génère de fortes hétérogénéités. Notre motivation a été d'étudier des techniques permettant de modéliser les transferts de fluide et de contaminants entre la roche et les fractures. Dans cette thèse, nous proposons deux méthodes numériques adaptées à deux échelles de fracturation. La présence d'un réseau de petites fractures est prise en compte par une modèle double porosité. Ce modèle est obtenu via une homogénéisation du milieu fracture, il est compose des équations standard d'écoulement et de transport dans un milieu poreux auxquelles on a ajoute un terme de couplage qui gère les échanges entre la roche et les fractures. Ce terme de couplage est calcule via la résolution de problèmes auxiliaires qui sont des problèmes de diffusion définis en chaque point de l'espace. Vu le grand nombre de ces problèmes, il est déterminant de les résoudre efficacement. Nous avons montre qu'une résolution analytique est préférable a une discrétisation classique. Pour traiter les fractures plus importantes, nous avons développe une méthode de décomposition de domaine. L’ouverture des fractures étant faible devant les dimensions du site de stockage, nous avons assimile les fractures a des interfaces a travers lesquelles la pression et la concentration sont continues et les flux sont discontinus, contrairement a une décomposition de domaine classique, il s'agit, ici, d'une condition de transmission non locale. La différence des flux provenant des sous-domaines intervient comme un terme source dans les fractures. Les modèles d'écoulement et de transport dans les sous-domaines et dans les fractures ont été obtenus via une analyse asymptotique.
Author: Amel Sboui
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
L'objectif de cette thèse est de modéliser et de développer des outils numériques adaptés à l'étude de l'écoulement des eaux souterraines ainsi que la propagation des polluants en milieux poreux. La motivation de ce travail est un benchmark du GDR Momas et de l'Andra pour la simulation de la propagations 3-D des radionucléides autour d'un stockage profond de déchets nucléaires. Premièrement on a construit une nouvelle méthode d'éléments finis mixtes sur un maillage formé d'hexaèdres généraux. La convergence de la méthode est prouvée et confirmée par des tests numériques. Deuxièment, nous présentons une méthode de discrétisation en temps pour une équation d'advection telle que des pas de temps différents sont utilisés dans différents sous-domaines afin de prendre en compte les hétérogèneités. Enfin une méthode numérique pour le calcul de transport de contaminants est proposée. Les techniques précédentes sont implémentées en 3-D et des résultats numériques sont présentés sur le benchmark 3-D champ lointain du GDR Momas et de l'Andra.
Author: Selim Ahmed
Publisher:
ISBN:
Category :
Languages : fr
Pages : 142
Get Book Here
Book Description
L’objectif principal de ce travail est de développer des outils numériques efficaces et fiables pour la modélisation de l’écoulement et du transport de masse en milieux poreux saturé. Dans une première partie, nous avons développé une technique pour améliorer la méthode EFMH sur les quadrangles. Cette technique a permis d’éliminer les oscillations rencontrées en utilisant la méthode EFMH standard pour des solutions à front raide et de gagner en précision et en temps de calcul. Dans la deuxième partie de la thèse un nouveau modèle numérique pour résoudre le systèmes d’équations couplées écoulement–transport a été développé. Ce modèle couple des discrétisations spatiale et temporelle efficaces en utilisant la méthode des lignes. Pour la discrétisation spatiale, nous avons combiné la méthode EFMHc pour l’écoulement, la méthode MPFA pour le terme dispersif et la méthode DG pour le terme convectif de l’équation de transport. Pour la discrétisation temporelle du système EDO résultant, deux solvers basés sur le schéma à ordre et pas de temps adaptatifs, BDF, sont employés. DLSODIS (resp, DASPK), utilise une méthode directe (resp, méthode itérative) pour résoudre les systèmes linéaires surgissant à chaque itération de Newton. Dans un premier temps, deux cas tests numériques standards aussi bien qu’une expérience à l’échelle du laboratoire sont simulés pour mettre en relief l’efficacité et la robustesse du modèle basé sur DLSODIS. Dans un second temps, le solver DASPK a été utilisé pour traiter les systèmes de grande échelle. Le nouveau modèle permet de contrôler l’erreur temporelle et de réduire plus que 30 fois le temps de calcul par rapport au modèle standard.
Author: Samira Oukfif
Publisher:
ISBN:
Category :
Languages : fr
Pages : 185
Get Book Here
Book Description
L’objectif principal de ce travail est le développement d’un modèle numérique qui permet de simuler le transport de masse dans les milieux poreux homogènes et hétérogènes. Cet outil numérique sera à terme un moyen de prédiction de la progression d’une pollution dans un sol et d’évaluation des risques de contamination des nappes phréatiques. Le modèle est basé sur l’équation de convection-dispersion avec une cinétique de dépôt/relargage de premier ordre. La résolution de cette équation en 1D et 2D est basée sur une méthode Lagrangienne, appelée méthode particulaire, qui utilise la technique de vitesse de dispersion. Les conditions aux limites (entrée du domaine) sont traitées avec une technique de particules images. La variation de porosité dans le milieu poreux due à la rétention et au détachement des particules est évaluée à partir de la relation de Kozeny-Carman. L’étude de sensibilité du modèle, effectuée dans différentes configurations lorsque les solutions analytiques sont fournies, montre une précision suffisante pour un choix adéquat de paramètres numériques. La validation du code de calcul est obtenue par l’ajustement d’essais de traçage en colonnes de laboratoire dans des conditions d’écoulement à débit constant ou à charge constante. Dans ce dernier cas, un couplage entre l’équation de transport et celle de l’écoulement (Equation de Darcy) est réalisé. La résolution de l’équation de l’écoulement est effectuée à l’aide d’un schéma numérique en différences finies sur une grille fixe. L’équation de transport est d’abord résolue pour calculer la concentration sur les particules, puis à l’aide d’un schéma séquentiel non itératif, on résout l’équation de l’écoulement pour évaluer la charge hydraulique et la vitesse sur les nœuds de la grille. Les échanges entre la grille et les particules sont assurés au moyen de fonctions d’interpolation. La simulation d’essais en injection instantanée ou continue a montré un bon ajustement entre les courbes de restitution calculées et mesurées, et notamment les essais où le dépôt et le relargage agissent simultanément. Les profils de porosité obtenus le long de la colonne, montrent une diminution de la porosité à l’entrée de la colonne, et qui est plus marquée dans le cas d’un écoulement à charge constante. Le modèle numérique est adapté à la simulation de l’érosion (suffusion) d’un sol en considérant uniquement le relargage, et l’ajustement d’essais d’érosion totale au laboratoire a montré un bon accord. Afin de simuler le transport de soluté dissous dans deux milieux homogène et hétérogène supposé infinis, la cinétique de dépôt/relargage est remplacée par un partage non linéaire de Freundlich. Pour tenir compte des hétérogénéités du champ de perméabilité aléatoire qui engendre des vitesses d’écoulement préférentielles, plusieurs simulations dite monte Carlo sont réalisées et montrent que la vitesse de dispersion offre une alternative intéressante (convergence plus rapide). Les résultats obtenus sont en accord avec les résultats de la littérature. Cette étude numérique a permis la mise en œuvre de la méthode particulaire pour simuler le transport, dépôt et relargage dans un milieu poreux fini. Néanmoins, le modèle de dépôt /relargage de premier ordre adopté peut être amélioré afin de prendre en compte le couplage entre les deux processus et notamment le seuil de détachement des particules.
Author: Loïc Dijoux
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
La représentation des écoulements et du transport de masse en milieu poreux trouve des applications diverses en hydrogéologie dont notamment l'étude de la propagation des polluants dans les sols ou encore l'étude des transferts d'eau salée dans les nappes phréatiques côtières. Ce manuscrit de thèse propose donc d'étudier ces phénomènes physiques à travers la modélisation numérique. Plusieurs ensembles de méthodes d'éléments finis sont présentées et discutées. L'accent est essentiellement mis sur la représentation mathématique des milieux fortement hétérogènes et anisotropes. Deux nouvelles méthodes de résolution numérique sont ainsi introduites, à savoir les méthodes H-RTm et H-RTp. Ces méthodes sont toutes deux issues de la technique d'hybridation appliquée aux éléments finis mixtes et aux éléments finis discontinus de Galerkin. Les bénéfices qu'elles apportent dans la représentation des phénomènes d'écoulement et de transport de masse en milieu poreux sont alors illustrés à travers des exemples numériques retrouvés dans la littérature.
Author: Fanilo Heninkaja Ramasomanana
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
Le devenir des polluants dans les sols constitue un enjeu environnemental majeur. Dans ce travail, nous apportons une contribution à quelques méthodes numériques pour la simulation de l'écoulement et du transfert de polluants en milieu poreux variablement saturés. La propagation d'un contaminant dans les milieux souterrains dépend en premier lieu des caractéristiques de l'écoulement qui le transporte. Dans la première partie de ce travail, nous présentons la méthode des éléments finis mixtes hybrides pour la résolution de l'équation de Richards. Une procédure de condensation de la masse est proposée pour éviter l'apparition d'oscillations non physiques, notamment lors de la simulation de problèmes d'infiltration dans un milieu initialement sec.Dans la deuxième partie de ce travail, la méthode ELLAM (Eulerian-Lagrangian Localized Adjoint Method) est utilisée pour la modélisation du transport réactif en milieux fortement hétérogènes. En effet, les résultats obtenus pour le transport linéaire, décrit par l'équation d'advection-dispersion, avec les ELLAM sont très encourageants. La méthode ELLAM permet (i) de s'affranchir des contraintes de discrétisations spatiale ettemporelle imposées avec les méthodes eulériennes classiques, (ii) de conserver la masse et (iii) de traiter toutes les conditions aux limites. Par ailleurs, nous proposons une nouvelle formulation des ELLAM (C_ELLAM) permettant d'éviter les oscillations numériques et de limiter la diffusion numérique générées parla formulation standard.Dans la dernière partie, le code de calcul élaboré avec la formulation C_ELLAM est utilisé pour la caractérisation de la macrodispersion dans les milieux hétérogènes. Pour ce faire, il est indispensable de disposer d'outils de simulation précis et efficaces car cette étude est basée sur une méthode Monte Carlo nécessitant la réalisation d'un très grand nombre de simulations sur des grilles de calcul de l'ordre du million de mailles. Les résultats obtenus sont comparés avec une étude antérieure basée sur le Random WalkParticle Method.
Author: Van-Linh Nguyen
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
L'effet de serre lié à l'émission de CO2 a conduit à des projets de stockage de ce gaz dans des formations réservoirs. Ces formations peuvent être traversées de failles et l'examen de la sûreté du stockage nécessite alors l'étude du risque de réactivation et de propagation de ces failles. Cette étude passe par des investigations approfondies portant sur des conditions de propagation des fractures sous sollicitations hydromécaniques. Cette thèse a pour objectif l'étude théorique et numérique de ces conditions ainsi que la simulation numérique de la propagation. La modélisation numérique des processus thermo-hydro-mécaniques dans les milieux poreux fracturés par la méthode des éléments finis (MEF) permet de simuler des phénomènes complexes et non linéaires. Les difficultés liées à l'intégration des équations d'échanges de fluide entre la fracture et la matrice environnante avec la MEF ont été résolues dans des travaux récents et nos simulations numériques ont pu être basées sur cette méthode. Dans un premier temps, nous avons modélisé l'écoulement transitoire dans et au voisinage d'une fracture soumise à une injection de fluide et nous avons étudié le facteur d'intensité des contraintes (FIC) à l'extrémité de la fracture dans le cadre de la théorie de la poroélasticité linéaire. Si les conditions d'injection sont maintenues constantes et la fracture n'évolue pas, l'écoulement tend vers un état stationnaire. Le FIC évolue au cours de la phase transitoire pour atteindre une valeur limite dans l'état stationnaire. La modélisation de l'écoulement transitoire est très coûteuse en temps de calcul et il est intéressant de trouver un moyen d'exploiter au mieux les résultats d'un calcul en état stationnaire. L'analyse théorique et les résultats des simulations numériques montrent en effet que le FIC calculé à l'état stationnaire peut fournir certaines bornes pour la propagation des fractures sous l'écoulement transitoire. Dans le cadre de la poroélasticité linéaire et de l'écoulement de Poiseuille dans les fractures, des expressions semi-analytiques pour le FIC à l'état d'écoulement stationnaire ont pu être dérivées. Pour des géométries simples, ces formules approximatives se révèlent efficaces pour discuter des conditions de propagation des fractures pour des cas typiques et simples de géométrie de la fracture et des conditions d'injection de fluide. Dans un deuxième temps, un Modèle de Fracture Cohésive (MFC) a été utilisé pour modéliser la propagation de fracture sur la base de l'endommagent. Ce modèle, basé sur un critère de rupture de Mohr-Coulomb modifié, permet de simuler l'endommagement de l'interface à la fois sous sollicitations en mode I et II. Une relation d'équivalence entre les paramètres de ce modèle et du modèle de Mécanique Linéaire de la Rupture (MLR) a été établie sur la base de la longueur de propagation de fracture sous des charges similaires. Cette relation permet l'extension de l'équivalence théorique entre MLR et MFC établie pour les matériaux fragiles et sur la base de critères énergétiques, à des matériaux quasi-fragiles et ductiles. On a d'ailleurs montré que le MFC permet de simuler certains phénomènes spécifiques tels qu'instabilités de propagation en mode I et II et le branchement de la fracture en mode II. Enfin, la prise en compte de la pression de fluide dans la fracture a permis d'obtenir un modèle de MFC couplé avec l'hydraulique qui a été implémenté dans un code numérique aux éléments finis en vue d'étudier la propagation des fractures sous sollicitations hydromécaniques. Des simulations numériques ont été réalisées afin d'étudier le risque de réactivation et de propagation des failles dans le contexte de stockage du CO2 en particulier dans une configuration de formation réservoir du type Bassin de Paris.
Author: Anis Younes
Publisher:
ISBN:
Category :
Languages : fr
Pages : 181
Get Book Here
Book Description
LA PREMIERE PARTIE DE CE TRAVAIL CONCERNE LA RESOLUTION NUMERIQUE DES EQUATIONS AUX DERIVEES PARTIELLES QUI REGISSENT LE TRANSFERT DE MASSE EN MILIEU POREUX. POUR LA RESOLUTION DE L'EQUATION DE L'ECOULEMENT EN MILIEU HETEROGENE, UNE NOUVELLE FORMULATION POUR LA METHODE DES ELEMENTS FINIS MIXTES AVEC MOINS D'INCONNUES ET SANS AUCUNE APPROXIMATION EST ETABLIE POUR UNE TRIANGULATION QUELCONQUE. DANS LE BUT DE LIMITER LA DISPERSION NUMERIQUE, LA METHODE DES ELEMENTS FINIS DISCONTINUS EST DEVELOPPEE POUR RESOUDRE LA PARTIE CONVECTION DE L'EQUATION DE TRANSPORT. LA SECONDE PARTIE DE CE TRAVAIL EST CONSACREE A L'ETUDE DES PROBLEMES D'ECOULEMENT ET DE TRANSPORT AVEC PRISE EN COMPTE DU CONTRASTE DE MASSE VOLUMIQUE ET/OU DE VISCOSITE. DANS CE BUT, UN MODELE NUMERIQUE FONDE SUR LES METHODES DES ELEMENTS FINIS MIXTES ET DISCONTINUS A ETE DEVELOPPE. UNE PREMIERE VERIFICATION DU MODELE EST REALISEE EN SIMULANT LE PROBLEME D'HENRY POUR LEQUEL UNE SOLUTION SEMI-ANALYTIQUE EXISTE. D'AUTRES SIMULATIONS DE CAS THEORIQUES LARGEMENT EXPLOITES DANS LA LITTERATURE (ELDER, DOME DE SEL) SONT EFFECTUEES ET PERMETTENT UNE BONNE ANALYSE DES PERFORMANCES DES SCHEMAS NUMERIQUES UTILISES. LA DERNIERE PARTIE DE CE TRAVAIL EST CONSACREE A LA SIMULATION D'EXPERIENCES DE TRANSFERT DE SOLUTES AVEC CONTRASTE DE MASSE VOLUMIQUE ET/OU DE VISCOSITE A L'ECHELLE DU LABORATOIRE EN DEUX ET TROIS DIMENSIONS.
Author: Ali Zidane
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
La circulation de l'eau souterraine dans les aquifères contenants des roches évaporitiques peut provoquer des problèmes géo-mécaniques tels que l'affaissement du sol ou l'effondrement. Dans ce travail, nous nous intéressons au développement de modèles numériques permettant de simuler les écoulements dans les milieux poreux et non poreux ainsi que les phénomènes de dissolution et d'évolution des fractures dans le temps. La première partie de cette thèse est consacrée au développement de nouvelles solutions pour la validation des modèles numériques simulant les écoulements densitaires en milieux poreux. La nouvelle procédure consiste à résoudre simultanément les deux systèmes d'écoulement et de transport en utilisant l'algorithme de Levenberg-Marquardt. L'utilisation de cette technique nous a permis de développer, pour la première fois, des solutions semi-analytiques d'intrusion d'eau salée dans le cas de faible diffusion ainsi que dans le cas d'un grand contraste de densité. Dans la deuxième partie de ce travail, nous nous sommes intéressés aux écoulements dans les fractures des roches évaporitiques. Un code de calcul a été développé pour résoudre ce système non linéaire en utilisant des méthodes numériques adaptées. Pour valider ce nouveau modèle, nous avons développé une solution semi-analytique pour les écoulements densitaires de Stokes. La troisième partie de ce travail est consacrée au transport avec dissolution de la roche salée. Dans un premier temps, nous avons étudié l'influence de différents paramètres sur la dissolution du sel dans l'Adler tunnel situé à une profondeur de 160 m dans la région de Bâle en Suisse. Dans un second temps, nous nous sommes intéressés à la simulation de l'évolution dune fracture sous l'effet de la dissolution. Le modèle numérique développé prend en compte les écoulements de Stokes ainsi que le transport de masse avec effets densitaires et la dissolution des parois de la fracture.
Author: Nicolas Trottier
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
Ce travail a pour objectif de développer et d'appliquer une méthode originale permettant de simuler l'écoulement dans un milieu poreux fracturé. Cette méthode repose sur une approche multicouches double continuum permettant de séparer le comportement des différents aquifères présents sur un site. La résolution des écoulements, basée sur la méthode des Eléments Finis de Crouzeix-Raviart, est associée à une méthode inverse (minimisation de type Quasi-Newton combinée à la méthode de l'état adjoint) et à une paramétrisation multi-échelle.La méthode est appliquée dans un premier temps sur l'aquifère fracturé du site expérimental de Poitiers. Les résultats montrent une bonne restitution du comportement de l'aquifère et aboutissent à des champs de transmissivité plus réguliers par rapport à ceux de l'approche simple continuum. L'application finale est réalisée sur le site de Cadarache (taille plus importante et données d'entrée moins denses). Le calage des deux aquifères présents sur le site est satisfaisant et montre que ceux-ci se comportent globalement de façon indépendante. Ce calage pourra être amélioré localement grâce à données de recharge plus fines.
