Modélisation numérique de l'écoulement et du transport en milieux poreux

Modélisation numérique de l'écoulement et du transport en milieux poreux

Author: Selim Ahmed

Publisher:

Published: 2009

Total Pages: 142

ISBN-13:

DOWNLOAD EBOOK

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.


Modélisation numérique du transport de masse et de la filtration dans les milieux poreux saturés

Modélisation numérique du transport de masse et de la filtration dans les milieux poreux saturés

Author: Samira Oukfif

Publisher:

Published: 2010

Total Pages: 185

ISBN-13:

DOWNLOAD EBOOK

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.


Modélisations et simulations numériques d'écoulements d'air dans des milieux micro poreux

Modélisations et simulations numériques d'écoulements d'air dans des milieux micro poreux

Author: Thanh Long Vu

Publisher:

Published: 2011

Total Pages: 0

ISBN-13:

DOWNLOAD EBOOK

Ce travail de thèse a pour objectif de simuler numériquement des écoulements de gaz dans des matrices poreuses dont les pores sont de taille micrométrique. On étudie l'influence des phénomènes de glissement hydrodynamique qui apparaissent lorsque la dimension caractéristique de micro- conduites est caractérisée par des nombres de Knudsen compris entre Kn = 0,01 et Kn = 0,1.Le mémoire de thèse est composé de cinq chapitres suivis d'une conclusion dans laquelle nous présentons quelques perspectives pour une suite de ce travail. Le chapitre I constitue le travail préliminaire de thèse qui s'est ensuite orienté vers des approches complémentaires. Le principe des méthodes d'homogénéisation périodique est d'abord exposé. Suit une présentation de deux méthodes de résolution dans l'espace de Fourier : l'approche en déformation et l'approche en contrainte. L'extension de ces méthodes à la résolution d'écoulements régis par l'équation de Stokes est ensuite décrite. Des applications aux cas d'écoulements à travers des réseaux de cylindres, avec condition d'adhérence ou avec condition de glissement, sont ensuite discutées. Deux techniques de modélisation des phénomènes de transport dans des milieux poreux saturés par un fluide monoconstituant sont présentées dans le second chapitre. La première est basée sur la méthode des développements asymptotiques, appelée aussi méthode d'homogénéisation. On explique que le processus consiste en trois étapes : description locale, localisation et description macroscopique. La seconde technique s'appuie sur la méthode de calcul de moyennes à l'échelle d'un VER. Le point de départ de cette méthode est basé sur des théorèmes donnant les expressions des moyennes de tous les opérateurs intervenant dans une équation de transport. Après une brève présentation du logiciel commercialisé que nous avons utilisé, nous exposons les études de convergence spatiale que nous avons effectuées et nous comparons nos solutions avec des résultats de la littérature dans le chapitre III. Diverses géométries sont considérées (allant de géométries planes à des empilements 3D de cubes ou de sphères).L'effet du glissement sur la perméabilité de milieux microporeux est abordé dans le chapitre IV. Le formalisme résultant de l'homogénéisation de structures périodiques est utilisé pour simuler numériquement des écoulements isothermes de gaz dans divers empilements de complexités croissantes. Les perméabilités sont déterminées en calculant les moyennes spatiales des champs de vitesses, solutions des équations de Stokes. Les valeurs obtenues en imposant des conditions d'adhérence sont comparées à celles obtenues avec des conditions de glissement du premier ordre. Dans le chapitre V, nous présentons des solutions pour des écoulements anisothermes et étudions l'effet du glissement sur la conductivité effective de milieux microporeux 2D et 3D. Dans ce chapitre, nous résolvons les équations de Navier-Stokes et de l'énergie en imposant des conditions de symétries dans une ou deux directions. A partir des solutions locales, sont calculées les moyennes intrinsèques des champs de vitesse et de température. Nous considérons des cas pour lesquels la condition d'équilibre thermique local peut être considérée comme satisfaite et d'autres correspondant à un non-équilibre thermique (NTLE). On détermine les conductivités de dispersion en fonction du nombre du Péclet et on montre l'influence du glissement sur les composantes longitudinales et transverses pour différentes porosités et longueur de glissement. Dans les cas NLTE, le coefficient macroscopique de transfert fluide-solide est aussi calculé


Modélisation du transport réactif en milieux poreux saturé

Modélisation du transport réactif en milieux poreux saturé

Author: Tony Khalil

Publisher:

Published: 2011

Total Pages: 177

ISBN-13:

DOWNLOAD EBOOK

Les sols et les sous-sols urbains peuvent être contaminés par des polluants divers. En effet, leur mobilité est régie par différents processus physico-chimiques de transport et de rétention/libération (la rétention tend à freiner la propagation des particules tandis que la libération accélère le transport des particules). Ainsi, il est nécessaire de comprendre les mécanismes du transport des polluants dans des sols et de prédire cette migration à long terme. Tout d’abord, nous détaillerons un modèle de transport et d’écoulement en milieu poreux saturé, qui est basé sur les EDPs. De plus, ces équations se réfèrent au transport de masse en milieu poreux saturé qui est régi en général par la classique équation de type convection-dispersion-diffusion. En bref, les équations de l'écoulement et du transport avec des effets de viscosité et la densité sont mises en œuvre dans les codes de différences et éléments finis en utilisant la Librairie d’éléments finis Diffpack. Ensuite, nous proposons de coupler le modèle de transport avec l’érosion interne. Pour atteindre ce but, nous avons proposé un modèle mathématique de l’érosion interne et nous l’avons validé par des tests expérimentaux. De même, nous avons montré que l’érosion influe directement sur le transport des particules. Ceci nous a amené à définir un couplage entre le transport, l’érosion et le transfert de polluants. En outre, le modèle développé est basé sur le transport réactif qui tente de simuler les processus d’écoulement-transport, les réactions physico-chimiques. Enfin, nous avons introduit le comportement des milieux déformables dans le modèle pour ainsi étudier le rôle de la déformation du sol dans le transport.


Modélisation des écoulements dans des milieux poreux fracturés par la méthode des équations aux intégrales singulières

Modélisation des écoulements dans des milieux poreux fracturés par la méthode des équations aux intégrales singulières

Author: Minh Ngoc Vu

Publisher:

Published: 2012

Total Pages: 0

ISBN-13:

DOWNLOAD EBOOK

Cette thèse est consacrée au développement d'une méthode numérique visant à modéliser des écoulements dans des milieux poreux fissurés, ainsi qu'à déterminer leur perméabilité effective à partir des avancements théoriques récents. En parallèle, elle a été aussi l'occasion de continuer sur la voie théorique et d'obtenir de nouveaux résultats sur ce plan.Les équations générales gouvernant l'écoulement dans de tels matériaux sont rappelées, et plus particulièrement, la conservation de la masse à l'intersection entre fissures est établie explicitement [132, 139]. Des solutions générales du potentiel sont proposées sous la forme d'une équation intégrale singulière décrivant l'écoulement dans et autour des fissures entourées par une matrice infinie soumise à un champ lointain [136, 139]. Ces solutions représentent le champ de pression dans le milieu infini en fonction de l'infiltration dans les fissures qui prennent en compte complètement l'interaction et l'intersection entre fissures. En considérant le problème d'une fissure super-conductrice, des solutions analytiques de l'écoulement ont été développées et ces solutions sont utilisées comme cas de référence pour valider la solution numérique. De plus, la solution théorique de ce problème dans le cas 3D permet de comparer le modèle d'écoulement de Poiseuille dans une fissure elliptique d'épaisseur nulle et le modèle d'inclusion ellipsoïdale aplatie soumise à l'écoulement de Darcy [140]. Des outils numériques ont été développés en se basant sur la méthode des équations intégrales singulières afin de résoudre les équations générales du potentiel [132, 180]. Cela permet, d'une part, de modéliser l'écoulement stationnaire dans un domaine poreux contenant un grand nombre de fissures et, d'autre part, de proposer une solution semi-analytique de l'infiltration dans une fissure isolée dépendant de la perméabilité de la matrice, de la conductivité de la fissure et de la variable géométrique de la fissure. Cette dernière est l'élément essentiel pour déterminer de la perméabilité effective d'un milieu poreux fissuré en utilisant des schémas d'homogénéisation. Ce modèle auto-cohérent révèle un seuil de percolation qui est alors applicable pour l'estimation de la perméabilité effective d'un matériau contenant un grand nombre de fissures. L'approche par sous-structuration permet l'extension de la solution générale du potentiel, écrite pour un domaine infini, à celle d'un domaine fini [181]. Une solution analytique de l'écoulement dans et autour d'une fissure partiellement saturée est établie, fondée sur la solution analytique pour la fissure super-conductrice. Celle-ci est alors utilisée pour estimer la perméabilité effective d'un milieu poreux fissuré non-saturé [141]. Le modèle de la perméabilité effective est appliqué dans le contexte du stockage géologique du CO2 en vue d'étudier le comportement d'une zone de faille constituée par un noyau argileux et des zones fissurées. La pression d'injection provoque l'augmentation de la pression interstitielle dans le réservoir. Cette surpression engendrée dans le réservoir peut affecter la perméabilité de zones fissurées ce qui conduit à des phénomènes hydromécaniques couplés. Les résultats de simulations numériques permettent d'évaluer, d'une part, le risque de la remontée de la saumure à l'aquifère supérieur, et d'autre part, le risque de l'initiation d'une rupture sur le plan de la faille.


Simulation numérique 3D d'Écoulement Multiphasiques Réactifs en Milieux Poreux

Simulation numérique 3D d'Écoulement Multiphasiques Réactifs en Milieux Poreux

Author: Mohamed Id Moulay

Publisher:

Published: 2019

Total Pages: 0

ISBN-13:

DOWNLOAD EBOOK

La modélisation du transport réactif est utilisée dans de nombreuses applications énergétiques et environnementales liées aux écoulements souterrains. La modélisation de tels problèmes conduit à un système hautement non linéaire d'EDP couplées à des équations différentielles ordinaires ou algébriques. Deux types d'approches pour la résolution numérique des problèmes de transport réactif sont largement utilisés dans la littérature. L'une est l'approche de séparation des opérateurs qui consiste à découpler les problèmes d'écoulement et de transport réactif. Ces derniers sont résolus séquentiellement à chaque pas de temps. L'autre stratégie est basée sur une approche entièrement couplée dans laquelle le système entier est résolu simultanément. Le but de la thèse de doctorat est le développement d'un schéma implicite en volumes finis pour la modélisation numérique d'écoulements multicomposants monophasiques et diphasiques avec transport réactif en milieu poreux.Deux nouveaux modules de transport réactif seront implémentés dans DuMuX, un simulateur libre pour les problèmes d'écoulements et de transport dans les milieuw poreux. Des simulations numériques bi et tridimensionnels comprenant des benchmarks et du calcul haute performance, seront effectuées pour valider les modules.


Modélisation des écoulements en milieu poreux hétérogènes 2D

Modélisation des écoulements en milieu poreux hétérogènes 2D

Author: Ahmad Al Bitar

Publisher:

Published: 2007

Total Pages: 150

ISBN-13:

DOWNLOAD EBOOK

Dans ce travail, on considère la modélisation des écoulements dans des hydro-systèmes comprenant des sols et des aquifères géologiquement complexes et hétérogènes. On considèrera par exemple le cas d'un aquifère côtier soumis à l'intrusion saline, avec couplage densitaire (eau douce / eau salée), phénomène auquel peuvent se greffer d'autres couplages (écoulements à saturation variable, couplages surface / souterrain). On choisit une approche ayant les caractéristiques suivantes : - le modèle est spatialement distribué afin de représenter l'hétérogénéité du milieu ; - le modèle est fortement couplé afin d'appréhender les écoulements dans leur complexité physique. On utilise dans ce but un modèle fortement intégré, à une seule équation générique de type EDP, basée sur une loi de Darcy généralisée permettant de décrire différents "régimes" d'écoulements la co-existant dans un même domaine, tout en conservant robustesse et efficacité. Le travail est divisé en trois parties : dans une première partie on élabore un nouveau modèle numérique 3D, pour la modélisation des écoulements en milieux poreux à densité variable dans l’hypothèse d'une interface abrupte. Ce nouveau modèle est basé sur des relations "effectives" non linéaires de saturation et de perméabilité, dans une équation d'écoulement de type Richards modifiée. La seconde partie correspond à l'élaboration et l'implémentation d'un modèle verticalement intégré d'intrusion saline en aquifère côtier, permettant d'étudier l'effet de l'hétérogénéité stochastique de l'aquifère. Le modèle, basé sur l'hypothèse "interface abrupte", est implémenté comme un module "2D" dans le code volumes finis BigFlow2D/3D. Le nouveau module 2D est utilisé pour analyser la variabilité de l"interface eau douce / eau salée par simulations stochastiques de type Monte Carlo à échantillonnage spatial (réalisation unique). Ces résultats sont comparés à nouvelle théorie, où l"interface aléatoire auto-corrélée est analysée par transformation de variable, combinée à une méthode de perturbation et à une représentation spectrale (Fourier / Wiener-Khinchine). Dans la troisième et dernière partie, on présente un modèle de couplage fortement "intégré" pour la modélisation des écoulements de surface et souterrain en hypothèses d'écoulement plan, verticalement hydrostatique. On s'intéresse au cas d'une vallée fluviale avec cours d'eau, plaine d'inondation, et nappe d'accompagnement. L'écoulement en surface est modélisé par l'équation d'onde diffusante et l'écoulement souterrain par l'équation de Dupuit-Boussinesq. Ce modèle couplé est appliqué à la vallée fluviale de la Garonne dans la région de Toulouse - Moissac (France). Cette application a nécessité l'élaboration d'une méthode d'interpolation géostatistique adaptée à l'élaboration d'un Modèle Intégré Numérique de Terrain ("MINT"), de façon à inclure le fond de la rivière au MNT topographique en haute résolution. Enfin, au-delà de cette application particulière, le modèle d'écoulement couplé surface / souterrain est généralisé au cas d'un couplage densitaire eau douce / eau salée, lorsque la nappe est sujette à l'intrusion saline au voisinage d'une embouchure ou d'un estuaire.


Modélisation tridimensionnelle par la méthode des éléments finis de l'écoulement et du transport en milieu poreux saturé

Modélisation tridimensionnelle par la méthode des éléments finis de l'écoulement et du transport en milieu poreux saturé

Author: Franck Jacob

Publisher:

Published: 1993

Total Pages: 184

ISBN-13:

DOWNLOAD EBOOK

LE TRAVAIL PRESENTE, EST RELATIF A LA CREATION D'UN MODELE PAR LA METHODE DES ELEMENTS FINIS TRIDIMENSIONNEL, POUR L'ECOULEMENT ET LE TRANSPORT D'ENERGIE OU DE SOLUTE EN MILIEU POREUX. LE PHENOMENE DE TRANSPORT PREND EN COMPTE LES REACTIONS CHIMIQUES DE PREMIERE ESPECE, LA VISCOSITE ET LA MASSE SPECIFIQUE DU FLUIDE SONT LINEAIREMENT DEPENDANTES DE LA CONCENTRATION OU DE LA TEMPERATURE. APRES UNE PRESENTATION DU PROBLEME PHYSIQUE, ON S'ASSURE QUE LE PROBLEME MATHEMATIQUE EST BIEN POSE, ET QUE DANS LE CADRE LINEAIRE IL EXISTE UNE SOLUTION UNIQUE, ET ON DONNE UNE DISCRETISATION DE CE PROBLEME EN ESPACE ET EN TEMPS. ENFIN, ON PRESENTE LES CAS TESTS MIS EN UVRE POUR S'ASSURER DE LA VALIDITE DU MODELE