Download Free Elements Finis Discontinus Multi Domaines En Temps Pour La Modelisation Du Transport En Milieu Poreux Sature Book in PDF and EPUB Free Download. You can read online Elements Finis Discontinus Multi Domaines En Temps Pour La Modelisation Du Transport En Milieu Poreux Sature and write the review.

Ce travail s’intéresse à la modélisation du transport de polluants en milieu poreux. L’équation hyperbolique caractérisant le transport convectif est résolue à l’aide de la méthode des éléments finis discontinus de Galerkin. La discrétisation spatiale est étudiée dans un premier temps et deux espaces d’approximation sont comparés pour des maillages déstructurés en 2-D. La discrétisation temporelle est abordée dans un second temps et deux alternatives au schéma explicite classique sont présentées afin de s’affranchir de la condition CFL : le schéma semi-implicite et le schéma explicite multi-domaines où différents pas de temps peuvent être utilisés localement. La dernière partie de ce travail porte sur la simulation des problèmes de transfert avec prise en compte du contraste de densité. L’équation de transport est alors couplée à celle de l’hydrodynamique. La discrétisation temporelle multi-domaines est implémentée en 3-D et s’avère redoutablement efficace dans ce genre de situations.
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
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.
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.
LE BUT DE CETTE THESE EST LA MISE AU POINT D'OUTILS ROBUSTES DE SIMULATION DU TRANSFERT DE POLLUANT EN MILIEUX POREUX COMPLEXES. CETTE SIMULATION NECESSITE AU PREALABLE UNE BONNE APPROXIMATION DU CHAMP DE VITESSE. LES MODELES D'ECOULEMENT DEVELOPPES UTILISENT LA METHODE DES ELEMENTS FINIS MIXTES HYBRIDES, DONT LA PARTICULARITE EST D'APPROCHER SIMULTANEMENT LA CHARGE PIEZOMETRIQUE ET LA VITESSE. LES DIFFERENTS TESTS REALISES EN MILIEUX HETEROGENES MONTRENT LA NETTE SUPERIORITE DE CETTE APPROXIMATION. UNE METHODE BASEE SUR LA COMBINAISON D'ELEMENTS FINIS DISCONTINUS (UTILISEE POUR LA DISCRETISATION DU TERME CONVECTIF) ET D'ELEMENTS FINIS MIXTES HYBRIDES (UTILISEE POUR L'APPROXIMATION DU TERME DISPERSIF), A ETE DEVELOPPEE POUR RESOUDRE L'EQUATION DE TRANSPORT. CETTE METHODE EST STABLE QUELLE QUE SOIT LA VALEUR DU NOMBRE DE PECLET. LES MODELES D'ECOULEMENT ET DE TRANSPORT SONT VALIDES A L'AIDE D'EXPERIENCES DE TRANSPORT DE CHLORURE DE SODIUM EN MILIEU HETEROGENE ET EN MILIEU HOMOGENE AVEC CONTRASTE DE MASSE VOLUMIQUE. AFIN DE POUVOIR UTILISER CES OUTILS SUR DES CAS REELS DE POLLUTION, LE MODELE D'ECOULEMENT A ETE COUPLE AVEC UN ALGORITHME D'ESTIMATION DE PARAMETRES. CET ALGORITHME EST BASE SUR UNE PARAMETRISATION MULTI-ECHELLE APPLIQUEE AUX TRANSMISSIVITES ET AUX CONDITIONS AUX LIMITES. CETTE PARAMETRISATION PERMET PAR RAFFINEMENT SUCCESSIF D'ADAPTER LA DISCRETISATION DU MILIEU AUX INFORMATIONS CONTENUES DANS LES MESURES (CHARGES PIEZOMETRIQUES, TRANSMISSIVITES)
Cette étude traite de la modélisation des écoulements d’eau en milieux poreux non saturés. Après avoir caractérisé le domaine d’étude, le modèle déterministe fondé sur l’équation non linéaire de Richards est décrit. Ses limites et d’autres alternatives de modélisation sont évoquées. La méthode des éléments finis mixtes hybrides est présentée et une technique originale de condensation de la masse est proposée pour éviter l’apparition d’oscillations non physiques, notamment lors de la simulation de fronts raides d’infiltration. L’estimation de la conductivité relative, à l’échelle de chaque maille du domaine discrétisé, est un aspect essentiel pour améliorer la précision des résultats. Différentes formulations sont décrites, et les comparaisons effectuées sur de nombreux cas tests conduisent à privilégier la moyenne pondérée. Une partie est également consacrée à la gestion des non linéarités. La forme mixte de l’équation de Richards est conservée, et les investigations aboutissent à un algorithme de sélection de la variable de résolution. Concernant la gestion du temps, les méthodes heuristiques traditionnelles peuvent être avantageusement remplacées par des méthodes a priori, qui relient le pas de temps à un calcul d’erreur temporelle. La méthode d’extrapolation sur l’ordre peut alors être efficacement combinée à un contrôleur intégral. Dans la dernière partie de ce mémoire, toutes les optimisations étudiées dans un contexte monodimensionnel sont généralisées dans un algorithme adapté aux écoulements bidimensionnels. Des cas tests sont proposés pour étudier son efficacité par rapport à des méthodes standards.
L'objet de cette thèse est l'étude et la mise en œuvre d'une méthode d'éléments finis multi-échelles pour la simulation d'écoulements miscibles en milieux poreux. La définition des fonctions de base multi-échelles suit l'idée introduite par F. Ouaki. La nouveauté de ce travail consiste à combiner cette approche multi-échelle avec des éléments finis de type Galerkine Discontinus (DG) de façon à pouvoir utiliser ces nouveaux éléments sur des maillages non-conformes composés de mailles de formes diverses. Nous rappelons, dans un premier temps, le principe des méthodes DG et montrons comment ces méthodes peuvent être utilisées pour discrétiser une équation de convection-diffusion instationnaire identique à celle rencontrée dans le problème d'écoulement considéré dans ce travail. Après avoir vérifié l'existence et l'unicité d'une solution à ce problème, nous redémontrons la convergence des méthodes DG vers cette solution en établissant une estimation d'erreur a priori. Nous introduisons, ensuite, les éléments finis multi-échelles non conformes et détaillons leur mise en œuvre sur ce problème de convection-diffusion. En supposant les conditions aux limites et les paramètres du problème périodiques, nous montrons une nouvelle estimation d'erreur a priori pour cette méthode. Dans une seconde partie, nous considérons le problème d'écoulement complet où l'équation considérée dans la première partie est résolue de manière couplée avec l'équation de Darcy. Nous introduisons différents cas tests inspirés de modèles d'écoulements rencontrés en géosciences et comparons les solutions obtenues avec les deux méthodes DG, à savoir la méthode classique utilisant un seul maillage et la méthode étudiée ici. Nous proposons de nouvelles conditions aux limites pour la résolution des problèmes de cellule qui permettent, par rapport à des conditions aux limites linéaires plus classiquement utilisées, de mieux reproduire les variations des solutions le long des interfaces du maillage grossier. Les résultats de ces tests montrent que la méthode multi-échelle proposée permet de calculer des solutions proches de celles obtenues avec la méthode DG sur un seul maillage et de réduire, de façon significative, la taille du système linéaire à résoudre à chaque pas de temps.
LES TRANSFERTS THERMOCONVECTIFS SONT ETUDIES DANS LES MILIEUX POREUX CONSTITUES PAR DES EMPILEMENTS ALEATOIRES DE FIBRES OU DE PARTICULES DE DIVERS MATERIAUX SATURES PAR DIFFERENTS FLUIDES INCOMPRESSIBLES. LA MODELISATION DU PROBLEME DYNAMIQUE EST BASEE SUR LE MODELE DE DARCY-FORCHHEIMER-BRINKMAN ET LE MODELE D'ERGUN DE LA PERMEABILITE EN TENANT COMPTE D'UNE LOI DE VARIATION DE LA POROSITE. LA MODELISATION DU PROBLEME THERMIQUE POUR DES PHASES SOLIDE ET FLUIDE EN EQUILIBRE THERMIQUE LOCAL EST DERIVEE DU MODELE DE ZEHNER ET SCHLUENDER POUR LA CONDUCTIVITE THERMIQUE MOYENNE ET DU MODELE DE HSU ET CHENG POUR LA CONDUCTIVITE THERMIQUE DISPERSIVE. LA SIMULATION NUMERIQUE EST DEVELOPPEE A L'AIDE DE LA METHODE DES ELEMENTS FINIS EN UTILISANT UN MAILLAGE LAGRANGIEN COMPORTANT DES ELEMENTS QUADRATIQUES COMPLETS. LES EQUATIONS DE CONSERVATION DE LA DYNAMIQUE ET DE LA THERMIQUE SONT RESOLUES A L'AIDE DES VARIABLES PRIMITIVES POUR UNE GEOMETRIE BIDIMENSIONNELLE CARTESIENNE. L'ETUDE DE REFERENCE DES PHENOMENES DE DISPERSION THERMIQUE EST REALISEE DANS LE CAS D'UN MILIEU POREUX CONSTITUE PAR DES PARTICULES SPHERIQUES (VERRE) DELIMITE PAR DEUX SURFACES PLANES PARALLELES ET SATURE PAR DE L'EAU. LES RESULTATS NUMERIQUES PERMETTENT DE VALIDER LA PROGRAMMATION DE LA METHODE D'APRES D'AUTRES RESULTATS NUMERIQUES ET EXPERIMENTAUX. L'ETUDE PRELIMINAIRE DES ECOULEMENTS ET DES TRANSFERTS AUTOUR D'UNE SURFACE CYLINDRIQUE INSEREE DANS UN MILIEU DE POROSITE CONSTANTE EST REALISEE POUR UNE VALIDATION COMPLEMENTAIRE DE LA METHODE NUMERIQUE D'APRES CERTAINS RESULTATS ANALYTIQUES. L'ETUDE GENERALE DES PHENOMENES DE DISPERSION THERMIQUE AUTOUR D'UNE SURFACE CYLINDRIQUE EST EFFECTUEE DANS LE CAS D'UN MILIEU POREUX CONSTITUE PAR UN EMPILEMENT DE PARTICULES SPHERIQUES (VERRE OU D'ALUMINE) SATURE PAR DE L'AIR. LES RESULTATS NUMERIQUES PERMETTENT DE RETROUVER LES PHENOMENES CARACTERISTIQUES OBSERVES AU COURS DES PLUS RECENTS TRAVAUX EXPERIMENTAUX. L'INFLUENCE DU DIAMETRE DES PARTICULES ET DE LA CONDUCTIVITE THERMIQUE DU MATERIAU SUR LES VARIATIONS DU NOMBRE DE NUSSELT MOYEN EST ETUDIEE DANS UN LARGE DOMAINE DU NOMBRE DE PECLET
Une étape primordiale dans la conception et l'optimisation des structures complexes est l'établissement d'un modèle numérique de base, affiné successivement par des essais expérimentaux pour être finalement validé. Cette phase de modélisation, essentielle pour une compréhension future du comportement du système sous différentes sollicitations, suppose le recours à un outil d'analyse numérique performant et maîtrisable, s'appuyant généralement sur la méthode des éléments finis. Cet ouvrage a pour dessein d'exposer les fondements de la méthode des éléments finis et de montrer les qualités - mais aussi les limites - de ce procédé qui constitue à l'heure actuelle la technique la plus répandue de discrétisation spatiale. Son originalité réside dans l'analyse méthodique des problèmes elliptiques du second ordre monodimensionnels, bidimensionnels à variable d'état scalaire et tridimensionnels à variable d'état vectorielle, depuis leur formulation forte classique jusqu'à l'approche locale par la méthode des éléments finis. Mathématiquement rigoureux sans sacrifier les aspects pratiques, l'ouvrage passe systématiquement en revue les formes intégrale, faible et discrète des classes de problèmes couramment rencontrés en mécanique appliquée pour aboutir à une élaboration unifiée d'un modèle d'éléments finis. Comme en témoignent les nombreux exemples et exercices simples qui jalonnent l'exposé, le livre s'adresse en priorité aux étudiants de début de deuxième cycle. Bien que conçu à la base comme support d'enseignement, il est aussi destiné aux chercheurs et ingénieurs praticiens qui désirent s'initier à la méthode des éléments finis.