Thèses : « Écoulement en milieux poreux – Simulation par ordinateur »

1

Laquerie, Nathalie. « Simulation numérique d'un écoulement diphasique multicomposant en milieu poreux ». Bordeaux 1, 1997. http://www.theses.fr/1997BOR10526.

Texte intégral

Résumé :

Le but de la thèse est de simuler le lessivage par injection d'eau d'un aquifère pollué par un hydrocarbure. L'accent est mis sur la modélisation où l'on prend en compte de façon explicite les échanges entre les phases, sans faire l'hypothèse d'équilibre local. Ce modèle nous a conduit à étudier un système de convection-diffusion-réaction dans lequel les temps caractéristiques régissant les divers mécanismes mis en jeu sont très différents. Dans un premier temps, nous avons fait l'étude monodimensionnelle, afin de cerner les principales difficultés : non linéarité, dégénérescence, caractère hyperbolique des équations. Afin de traiter aussi précisément que possible ces trois phénomènes, un schéma basé sur la méthode des pas fractionnaires est mis en place. Pour résoudre la partie convection, nous avons utilisé les schémas de Roe et de Harten, puis un schéma basé sur la méthode MUSCL introduite par Van Leer. La partie diffusion est traitée de manière semi-implicite par un schéma aux différences finies centrées. La partie réaction se ramène à l'étude d'un système d'équations différentielles ordinaires, qui est traité de façon semi-implicite. Dans ces deux dernières parties, les non-linéarités sont résolues à l'aide de la méthode de Newton. Ensuite, le schéma proposé est étendu à la dimension deux. On présente quelques tests numériques pour prouver la faisabilité de la méthode

Styles APA, Harvard, Vancouver, ISO, etc.

2

Chupin, Olivier. « Écoulement et transport couplés en milieux poreux saturés : application à l'injection des sols ». Nantes, 2005. http://www.theses.fr/2005NANT2096.

Texte intégral

Styles APA, Harvard, Vancouver, ISO, etc.

3

Dang, The Ba. « Contribution à l'étude des écoulements triphasiques en milieux poreux : évaluation des modèles : exemple d'application à un réservoir de pétrole ». Toulouse, INPT, 2003. http://www.theses.fr/2003INPT001H.

Texte intégral

Résumé :

L'étude des écoulements triphasiques en milieu poreux présente un intérêt technico-économique et scientifique évident. Du point de vue technico-économique ils jouent un rôle primordial dans l'exploitation de certains réservoirs de pétrole. Du point de vue scientifique ils offrent l'opportunité de s'interroger sur le rôle conjoint des fractions volumiques et des propriétés interfaciales, des phases en présence, dans la modélisation d'une classe particulière de milieux hétérogènes. Issu des difficultés rencontrées dans l'exploitation des gisements de pétrole vietnamiens, le présent travail s'efforce de répondre principalement aux questions suivantes : où en est-on, aujourd'hui, dans les tentatives de modélisation des écoulements triphasiques en milieux poreux et notamment des lois constitutives : perméabilités relatives et pression capillaires ? Quel crédit peut-on accorder à ces lois ? Et enfin, quelle est l'incidence de leur mise en oeuvre dans la prédiction des écoulements ? Pour répondre à ces questions, un recensement et une analyse aussi exhaustifs que possible des lois de fermeture, proposées dans la littérature, ont tout d'abord été effectué. Puis ces lois ont été évaluées par comparaison avec des données expérimentales existantes. L'élaboration d'un code de calcul, destiné à la résolution du système d'équations triphasiques tridimensionnelles a ensuite permis d'apprécier la sensibilité des solutions à ces différentes lois. Dans un premier temps ce code a été exploité dans une version unidimensionnelle, les résultats des simulations étant alors comparés à des résultats d'expériences effectués en laboratoire, puis dans sa version intégrale pour estimer l'influence des différentes lois sur l'historique de production d'un gisement réel.

Styles APA, Harvard, Vancouver, ISO, etc.

4

Rachmat, Sudjati H. « Modélisation numérique d'écoulements polyphasiques en milieu poreux : application à l'étude des champs de pressions et de saturaitons d'un gisement de pétrole ». Lyon 1, 1987. http://www.theses.fr/1987LYO19002.

Texte intégral

Résumé :

Developpement d'un modele numerique pour representer les ecoulements polyphasiques instationnaires dans un milieu poreux bidimensionnel, en presence de perturbations simulant, un puit de production ou d'injection. Les fluides sont du petrole, de l'eau et du gaz que l'on considere comme des fluides faiblement compressibles

Styles APA, Harvard, Vancouver, ISO, etc.

5

Id, Moulay Mohamed. « Simulation numérique 3D d'Écoulement Multiphasiques Réactifs en Milieux Poreux ». Thesis, Pau, 2019. http://www.theses.fr/2019PAUU3015.

Texte intégral

Résumé :

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
Reactive transport modeling is used in many energy and environmental applications related to subsurface flows. Modeling such problems leads to a highly nonlinear system of PDEs coupled with algebraic or ODEs. Two types of approaches for the numerical solving of reactive transport problems are widely used in the literature. One is the operator-splitting approach which consists in splitting the flow and reactive transport problems. These latter are solved sequentially at each time step. The other strategy is based on the fully coupled approach in which the entire system is solved simultaneously. The goal of the PhD thesis is the development of a fully coupled fully implicit finite volume scheme for numerical modeling of single and two-phase multicomponent flows with reactive transport in porous media. New reactive transport modules will be implemented in DuMuX, a free and open-source simulator for flow and transport processes in porous media. Numerical simulations for 2D and 3D including benchmark tests and high performance computing will be performed to validate the modules

Styles APA, Harvard, Vancouver, ISO, etc.

6

Verdière, Sophie. « Méthodes numériques de double maillage pour la simulation d'écoulements polyphasiques dans les milieux poreux ». Pau, 1997. http://www.theses.fr/1997PAUU3004.

Texte intégral

Résumé :

Les réservoirs pétroliers sont constitués de roches fortement hétérogènes. Les modèles géologiques ainsi générés utilisent un nombre très important d'éléments ou mailles. Pour des raisons de coût de calcul, la simulation numérique des écoulements dans les réservoirs nécessite de travailler sur un nombre de mailles plus réduits. La méthode classique consiste à déterminer le maillage réservoir en mettant à l'échelle les paramètres pétrophysiques. Cette démarche a l'inconvénient de ne pas tenir compte de l'évolution au cours du temps des variables du problème. Pour pallier ce défaut, on propose d'avoir recours à une homogénéisation pendant la résolution du problème. La méthode de double maillage consiste à résoudre, pour un système couple pression-saturation, chacune des équations du système avec une discrétisation en temps et en espace spécifique. L'appliquer à un problème diphasique revient à résoudre l'équation en pression (parabolique) sur un maillage plus grossier que l'équation en saturation (hyperbolique). Par rapport a un schéma Impes classique, il faut : 1) assurer le passage des résultats de la résolution implicite de l'équation en pression pour faire évoluer la saturation sur le maillage fin ; 2) une fois la saturation mise à jour, éventuellement après plusieurs pas de temps locaux, on calcule les paramètres homogénéisés nécessaires pour la prochaine étape du calcul en pression en tenant compte de la distribution au cours du temps des saturations. Le travail de ces trois ans a permis in fine de montrer la validité de la méthode de double maillage non seulement de manière numérique mais aussi théorique. En effet, la méthode a été validée numériquement sur des écoulements diphasiques incompressibles en milieux hétérogènes que les rapports de mobilités soient favorables ou non. De plus, une démonstration de convergence a assuré la validité théorique de la méthode pour un cas simplifié homogène (système elliptique/hyperbolique).

Styles APA, Harvard, Vancouver, ISO, etc.

7

Agnaou, Mehrez. « Une étude numérique des écoulements mono et diphasique inertiels en milieux poreux ». Thesis, Paris, ENSAM, 2015. http://www.theses.fr/2015ENAM0050/document.

Texte intégral

Résumé :

: Ce travail concerne l'écoulement inertiel en milieu poreux rencontré dans diversessituations telles que les écoulements autour des puits pour la récupération pétrolière, lesécoulements dans les réacteurs catalytiques, etc. En régime stationnaire, les différents modèlesmacroscopiques pour décrire ces écoulements inertiels (non-linéaires) demeurent encore sujetsà débat. Ces modèles consistent en une loi de Darcy corrigée de termes dont la dépendancevis à vis de la vitesse de filtration relève du régime d'écoulement. Dans ce travail, une attentionparticulière est portée tout d'abord à l'étude numérique (DNS), sur des structures modèles, de lalimite de stationnarité de l'écoulement monophasique newtonien qui correspond à la premièrebifurcation de Hopf, caractérisée par un nombre de Reynolds critique. La connaissance de cettelimite est cruciale puisqu'elle détermine le domaine de validité des modèles macroscopiquesstationnaires pertinents. Dans un deuxième temps, la dépendance de la déviation (inertielle) àla loi de Darcy par rapport aux propriétés de la structure poreuse (forme des grains, désordre)et à l'orientation de l'écoulement est étudiée dans le cas de structures 2D et 3D. Les propriétéseffectives de la structure à l'échelle macroscopique sont déterminées à partir de la résolutionnumérique des problèmes de fermeture associés au modèle macroscopique obtenu par prisede moyenne des équations de Navier-Stokes. Afin de déceler l'origine de cette déviation et sesdifférentes formes, l'évolution de la structure microscopique de l'écoulement en fonction dunombre de Reynolds est analysée. Plus particulièrement, le rôle des zones de recirculation, etles corrélations avec la courbure des lignes de courant multipliée par l’énergie cinétique localeet la variation de l’énergie cinétique le long de ces lignes sont étudiés. La dernière partie dutravail est consacrée à une étude numérique, toujours dans des situations modèles, de ladéviation à la loi de Darcy généralisée dans le cas de l'écoulement diphasique inertiel
This work focuses on inertial flow in porous media encountered in differentindustrial situations such as flow around wells in oil recovery, flow in filters and in columns ofreactors for chemical engineering, etc. In stationary flow regime, the different macroscopicmodels describing inertial (non-linear) flow are still discussed. These models consist in theDarcy’s law with correction extra terms whose dependence upon the filtration velocity is afunction of the flow regime. In this work, a particular attention is attributed first to the numericalinvestigation (DNS), on model structures, of the limit of one phase Newtonian stationary flowwhich corresponds to the first Hopf bifurcation, characterized by a critical Reynolds number.The knowledge of this limit is crucial since it establishes the ranges of validity of the relevantmacroscopic stationary models. In a second step, the dependence of the deviation (inertial)from Darcy’s law on the properties of the porous structure (grains shape, disorder) and on theorientation of the flow is analyzed in 2D and 3D situations. The effective properties of thestructure and the flow at the macroscopic scale are obtained from the numerical resolution ofthe closure problems associated to the macroscopic model obtained from an up-scalingprocedure (volume averaging) of the Navier-Stokes equations. In order to identify the origin ofthe deviation and its different forms, the variation of the microscopic flow structure with theReynolds number is analyzed. More specifically, the role of the recirculation zones, and thecorrelations with flow streamlines curvature multiplied by the local kinetic energy and thevariation of the kinetic energy along these lines are studied. The last part of the work isdedicated to a numerical investigation of the deviation from the generalized Darcy’s law in thecase of two phase inertial flow

Styles APA, Harvard, Vancouver, ISO, etc.

8

Amir, Laila. « Modèles couplés en milieux poreux : transport réactif et fractures ». Phd thesis, Université Paris Dauphine - Paris IX, 2008. http://tel.archives-ouvertes.fr/tel-00373688.

Texte intégral

Résumé :

Cette thèse porte sur la simulation numérique de modèles couplés pour l'écoulement et le transport dans les milieux poreux. Nous présentons une nouvelle méthode de couplage entre les réactions chimiques et le transport en utilisant une méthode de Newton-Krylov, et nous étudions également un modèle d'écoulement en milieu fracturé qui traite l'intersection des fractures par une méthode de décomposition de domaine.
Ce travail est divisé en trois parties : la première partie contient une analyse de différents schémas numériques pour la discrétisation des problèmes d'advection-diffusion, notamment par une technique de séparation d'opérateurs, ainsi que leur mise en oeuvre informatique, dans un code industriel.
La deuxième partie, qui est la contribution majeure de cette thèse, est consacrée à la modélisation et à l'implémentation d'une méthode de couplage globale pour le transport réactif. Le système couplé transport-chimie est décrit, après discrétisation en temps, par un système d'équations non linéaires. La taille du système sous-jacent, à savoir le nombre de points de grille multiplié par le nombre d'espèces chimiques, interdit la résolution du système linéaire par une méthode directe. Pour remédier à cette difficulté, nous utilisons une méthode de Newton-Krylov qui évite de former et de factoriser la matrice Jacobienne.
Dans la dernière partie, nous présentons un modèle d'écoulement dans un milieu fracturé tridimensionnel, basé sur une méthode de décomposition de domaine, et qui traite l'intersection des fractures. Nous démontrons l'existence et l'unicité de la solution, et nous validons le modèle par des tests numériques.

Styles APA, Harvard, Vancouver, ISO, etc.

9

Dimassi, Mohamed. « Modélisation de la structure poreuse du non-tissé par des techniques de calcul avancé ». Electronic Thesis or Diss., Lille 1, 2009. http://www.theses.fr/2009LIL10017.

Texte intégral

Résumé :

Les non-tissés sont des matériaux dont la structure poreuse est difficilement qualifiable ce qui rend la compréhension de leurs propriétés hydrauliques complexe. Nous avons opté pour l'observation directe du phénomène physique par analyse d'images de la dynamique de l'ascension capillaire dans un non-tissé. Pour cela nous avons mis en place un banc de test à l'aide d'une caméra numérique et d'un éclairage bien adapté pour faire des essais sur les non-tissés. Les algorithmes que j'ai développés et déployés, permettent une extraction des informations nécessaires à la compréhension et à la simulation du phénomène de la montée capillaire dans le réseau poreux du non-tissé: -un algorithme permettant l'extraction des propriétés métriques et géométriques des pores contenus dans un échantillon. -un algorithme permettant l'extraction de la propriété hydraulique de l'échantillon en fonction du temps. Cette propriété dynamique correspond à une propriété fonctionnelle de choix pour ces matériaux destinés à l'hygiène industrielle et domestique. -un algorithme permettant la détermination des paramètres de structure les plus pertinents en se basant sur un modèle d'arbre de décision. -un algorithme permettant la détermination et la simulation du comportement dynamique d'un liquide très mouillant au sein du réseau poreux. Cette approche tire son originalité de l'importance qui a été accordé au modèle expérimental et aux règles tirées de ce dernier afin de construire la simulation. Dans ce travail, nous partons d'observations et de constatations expérimentales pour en déduire une description algorithmique fidèle à la réalité. L'approche est appuyée par des théories existantes
Non-wowen are materials with a complicate porous structure. lndeed, the tangle of fibers and their consolidation, make the comprehension of their hydraulic properties complicated. After having explored three important tracks of research in this field, we chose the direct observation of the physical phenomenon by image analysis and especially the dynamics of the capil!ary rise in a non-woven material. For that, we set up a test bed using a digital camera and a well-adapted lighting to carry out tests on the non-woven materials. I developed and deployed algorithms, which allow an extraction of data in order to understand and simulate the capillary rise in the porous network ofthe non-woven material. These algorithms can be classified in the following way: -an algorithm for the extraction of metric and geometrical properties of pores contained in a non-woven sample.These data represent the structural parameters. -an algorithm which a!lows the extraction of the dynamic hydraulic property. This dynamic property corresponds to a functional property for these materials intended to be used for industrial and domestic hygiene, -an algorithm that determines the most relevant structural parameters based on a decision tree approach, -an algorithm allowing the determination and the simulation of a very wetting Iiquid's dynamic behaviour within the porous network. This approach draws its originality. from the importance attached to the, experimental model and the rules drawn from the Iatter in order.to build the simulation. ln this work, we start from experimental observatIons to deduce an algonthmlc descnptIon faithful to reality.The approach is supported by existing theones

Styles APA, Harvard, Vancouver, ISO, etc.

10

Berton, Gäelle. « Modélisations multi-échelles de structures et d'écoulements engendrés à petits nombres de Reynolds en amont et en aval de milieux poreux fibreux ». Poitiers, 1999. http://www.theses.fr/1999POIT2305.

Texte intégral

Résumé :

Des simulations numeriques d'ecoulements de fluides newtoniens incompressibles en regime permanent sont realisees par une methode d'elements finis en presence de milieux poreux de type fibreux. La resolution associe les equations de navier-stokes et la loi de darcy pour plusieurs conditions de frontieres. Suivant les cas, le milieu poreux est homogene, heterogene, ou constitue de bandes poreuses paralleles representant des espacements entre fils ou des endommagements accidentels. L'analyse du champ des vitesses permet de mettre en evidence l'influence du milieu fibreux sur l'ecoulement et son evolution suivant les changements des parametres hydrauliques (permeabilite) et structuraux tels que le nombre de fentes, la porosite d'ouverture. . . L'etude des variations du champ des pressions permet d'evaluer les pertes de pressions occasionnees par la presence du milieu poreux selon ses differentes configurations. A partir du champ des contraintes, la detection de zones de contraintes particulieres comme les zones de cisaillement pur et d'elongation pure, est un moyen pour determiner le comportement de particules deformables chargeant le fluide. Dans le cas de structures fibreuses non tissees, la methode des boites a permis d'identifier pour certaines leur nature fractale. De nouveaux modeles de reseaux de droites sont proposes. Leur construction est soit reguliere, soit fractale, soit de type aleatoire. La determination des parametres hydrauliques de ces modeles permet de les reinjecter dans la simulation numerique decrite ci-dessus.

Styles APA, Harvard, Vancouver, ISO, etc.

11

Ait, Hammou Oulhaj Ahmed. « Conception et analyse de schémas non-linéaires pour la résolution de problèmes paraboliques : application aux écoulements en milieux poreux ». Thesis, Lille 1, 2017. http://www.theses.fr/2017LIL10215/document.

Texte intégral

Résumé :

L'objectif de cette thèse est de concevoir et d'analyser des schémas numériques performants pour la simulation d'écoulements complexes en milieux poreux. Dans un premier temps nous proposons un schéma CVFE (Control Volume Finite Element) non-linéaire pour approcher la solution de l'équation de Richards anisotrope. La mobilité d'arête est gérée à l'aide d'une procédure de décentrement. On montre d'abord que ce schéma est non-linéairement stable, qu'il admet (au moins) une solution discrète et que la saturation est bornée entre 0 et 1. Ce schéma converge sans restriction sur le maillage. Enfin, en vue de mettre en évidence l'efficacité, la stabilité et la robustesse de la méthode, nous réalisons des tests numériques dans des cas isotropes et anisotropes. Dans un second temps on étudie un schéma Volumes finis (avec décentrement des mobilités) pour un modèle d'intrusion saline. Il préserve au niveau discret les principales propriétés du problème continu: l'existence de solutions discrètes positives, la décroissance de l'énergie et le contrôle de l'entropie et sa dissipation. Nous montrons que ce schéma converge. De plus, nous illustrons numériquement le comportement du modèle. Enfin nous étudions le comportement en temps long d'un modèle d'intrusion saline. Il s'agit d'identifier les états stationnaires qui sont les minimiseurs d'une énergie convexe. On montre pour le problème continu l'existence et l'unicité des minimiseurs de l'énergie, que les minimiseurs sont des états stationnaires et que ces états stationnaires sont radiaux et uniques. Nous donnons une illustration numérique des états stationnaires et nous exhibons le taux de convergence
This thesis is focused on the design and the analysis of efficient numerical schemes for the simulation of complex flows in porous media. First, we propose a nonlinear Control Volume Finite Element scheme (CVFE) in order to approximate the solution of Richards equation with anisotropy. This scheme is based on a suitable upwinding of the mobility which allows the negative transmissibility coefficients. We prove the nonlinear stability of the scheme, that there exists (at least) one discrete solution and that the saturation belongs to the interval [0,1]. Moreover, the convergence of the method is proved as the discretization steps tend to 0. We give some numerical experiments on isotropic and anisotropic cases illustrate the efficiency of the method. Second, we propose and analyze a finite volume scheme based on two-point flux approximation with upwind mobilities for a seawater intrusion model. The scheme preserves at the discrete level the main features of the continuous problem, namely the nonnegativity of the solutions, the decay of the energy and the control of the entropy and its dissipation. We show the convergence of this scheme. Numerical results are provided to illustrate the behavior of the model. Finally, the large time behaviour of the seawater intrusion model is studied. The goal is to identify the steady states which are the minimizers of a convex energy. We prove for the continuous problem the existence and uniqueness of the minimizers of the energy, that the minimizers are stationary states and that these stationary states are radial and unique. We give numerical illustrations of the stationary states and we exhibit the convergence rate

Styles APA, Harvard, Vancouver, ISO, etc.

12

Charton, Sophie. « Étude et modélisation en régime transitoire d'un écoulement de gaz à haute pression dans un milieu poreux réactif ». Vandoeuvre-les-Nancy, INPL, 1997. http://www.theses.fr/1997INPL012N.

Texte intégral

Résumé :

Le présent travail consiste en l'étude et la modélisation d'un procédé dans lequel un réacteur hétérogène est traversé, à débit variable, par un gaz sous haute pression (plusieurs centaines de bars) provenant d'un réservoir fermé. Le gaz, qui réagit avec le solide poreux de manière exothermique, est collecté, à sa sortie du réacteur, dans un second réservoir initialement vide, connecté au réacteur par un tube long et de faible section. Un modèle et un outil de simulation du procédé, prenant en compte le couplage du transfert de matière avec la réaction chimique et le transfert thermique, ont été développés. Les paramètres physico-chimiques qui caractérisent le système et interviennent dans le modèle ont été mesurés (perméabilité et conductivité thermique effective du réacteur, cinétiques d'échange de matière, dispersion axiale, construction de la courbe de Van Deemter à vitesse variable). Une maquette instrumentée du procédé, incluant le couplage du réacteur et d'un spectromètre de masse quadripolaire, a été mise au point. Ce dispositif, qui permet de suivre en continu la pression, la température et la composition du système au cours de la détente rapide du gaz, a été utilisé pour valider le modèle. L’outil de simulation a été exploité pour étudier la sensibilité du procédé à certains paramètres opératoires (pression initiale du gaz, porosité du réacteur, etc. . . . ). Les simulations ont de plus permis de prévoir l'influence de la dégradation du solide poreux sur les performances du système.

Styles APA, Harvard, Vancouver, ISO, etc.

13

Ben, Slimane Férid. « Méthodologie d'étude et modélisation du transport de polluants organiques en milieux poreux naturels : application à l'étude de deux cas de pollutions issues de l'industrie du bois ». Vandoeuvre-les-Nancy, INPL, 1997. http://www.theses.fr/1997INPL081N.

Texte intégral

Résumé :

Le présent travail vise à synthétiser l'état des connaissances en matière de modélisation du transfert de polluants organiques en milieu poreux naturel. Il s'articule autour des deux interrogations suivantes : - devant le nombre croissant de modèles de transfert de polluants, de plus en plus spécifiques et complexes, quels sont les critères de sélection à prendre en considération afin de disposer d'outils prédictifs relativement simples (requérant le moins de paramètres) ? - à l'issue d'une étude historique approfondie, d'un diagnostic classique de site pollué et d'une estimation des paramètres d'entrée des modèles de transport et d'écoulement, obtenue à l'issue d'une recherche bibliographique détaillée, peut-on, à l'aide de simulations relativement simples, valider ou approcher les valeurs observées lors du suivi ? Pour ce faire, dans un premier temps, un modèle de transfert a été sélectionné à l'issue d'un inventaire des codes disponibles sur le marché et dans la littérature. Dans un second temps, le code sélectionné a été appliqué à deux cas de pollutions issues de l'utilisation de produits de préservation du bois : la première due au lindane, la seconde due à la créosote. Les résultats des simulations ont été par la suite confrontés aux valeurs de concentrations relevées à l'issue du suivi réalisé sur ces deux sites.

Styles APA, Harvard, Vancouver, ISO, etc.

14

Puaux, Grégory. « Simulation numérique des écoulements aux échelles microscopique et mésoscopique dans le procédé RTM ». Phd thesis, École Nationale Supérieure des Mines de Paris, 2011. http://pastel.archives-ouvertes.fr/pastel-00667411.

Texte intégral

Résumé :

Le procédé " Resin Transfer Molding " consiste à injecter un polymère thermodurcissable à travers des fibres de renfort, qui polymérise ensuite pour former une pièce composite. L'objectif est de modéliser numériquement l'écoulement à différentes échelles, celle des mèches composants le tissu (mésoscopique) et celle des fibres composants les mèches (microscopique). Cette thèse se décompose en deux parties. La première concerne le calcul de perméabilité d'un volume élémentaire représentatif par prise de moyenne des champs de vitesse et de pression, dans le cadre de la méthode d'immersion de domaines. A l'échelle microscopique, le calcul de perméabilité a été validé en utilisant des lois analytiques. A l'échelle mésoscopique, un couplage entre les équations de Stokes et de Darcy, pour les écoulements entre les mèches et à l'intérieur des mèches a été effectué. La seconde partie aborde l'étape d'imprégnation du renfort à l'échelle microscopique, ce qui inclut la modélisation de l'avancée du front de matière entre les fibres avec prise en compte de la tension de surface. Nous avons implémenté des méthodes pour prendre en compte les phénomènes capillaires. Une méthode d'imposition de l'angle de contact statique de l'interface avec une surface solide, et la condition aux limites de Navier permettant d'imposer un glissement à la ligne triple, ont été implémentées et validées. Tous les développements ont été effectués dans une méthode éléments finis mixtes linéaires en vitesse pression, stabilisée par une fonction bulle (MINI-élément), et en utilisant la méthode d'immersion de domaines.

Styles APA, Harvard, Vancouver, ISO, etc.

15

Puaux, Grégory. « Simulation numérique des écoulements aux échelles microscopique et mésoscopique dans le procédé RTM ». Phd thesis, Paris, ENMP, 2011. http://www.theses.fr/2011ENMP0055.

Texte intégral

Résumé :

Le procédé " Resin Transfer Molding " consiste à injecter un polymère thermodurcissable à travers des fibres de renfort, qui polymérise ensuite pour former une pièce composite. L'objectif est de modéliser numériquement l'écoulement à différentes échelles, celle des mèches composants le tissu (mésoscopique) et celle des fibres composants les mèches (microscopique). Cette thèse se décompose en deux parties. La première concerne le calcul de perméabilité d'un volume élémentaire représentatif par prise de moyenne des champs de vitesse et de pression, dans le cadre de la méthode d'immersion de domaines. A l'échelle microscopique, le calcul de perméabilité a été validé en utilisant des lois analytiques. A l'échelle mésoscopique, un couplage entre les équations de Stokes et de Darcy, pour les écoulements entre les mèches et à l'intérieur des mèches a été effectué. La seconde partie aborde l'étape d'imprégnation du renfort à l'échelle microscopique, ce qui inclut la modélisation de l'avancée du front de matière entre les fibres avec prise en compte de la tension de surface. Nous avons implémenté des méthodes pour prendre en compte les phénomènes capillaires. Une méthode d'imposition de l'angle de contact statique de l'interface avec une surface solide, et la condition aux limites de Navier permettant d'imposer un glissement à la ligne triple, ont été implémentées et validées. Tous les développements ont été effectués dans une méthode éléments finis mixtes linéaires en vitesse pression, stabilisée par une fonction bulle (MINI-élément), et en utilisant la méthode d'immersion de domaines
In the Resin Transfer Molding (RTM) process, a thermoset resin is injected into a fibrous woven reinforcement. This resin polymerises to form a composite part. The goal is to model the flow at different scales: the tow scale (mesoscopic) and the fibre scale (microscopic), using numerical simulation. This thesis is divided in two parts. The first one is about permeability computation of a representative elementary volume (REV), applying an averaging method on pressure and velocity fields in the frame of the immersed domains method. At microscopic scale, permeability computation is validated using analytical laws. At mesoscopic scale, Stokes and Darcy equations, for the flows between and into yarns, have been implemented. The second part concerns the impregnation of the reinforcement at microscopic scale. It's includes modelling of the flow front movement between fibres with surface tension forces. We have implemented methods to take into account for the capillary phenomenon. We have implemented and validated methods to impose static contact angle and slip of the contact line. All developments are made using finite elements method with a velocity pressure formulation stabilised by a bubble function (MINI-element), and using the immersed domains method

Styles APA, Harvard, Vancouver, ISO, etc.

16

Ben, Cheikh Mohamed Najeh. « Effet du contraste de la masse volumique sur le déplacement et l'infiltration de substances non réactives en milieu poreux saturé ». Vandoeuvre-les-Nancy, INPL, 1998. http://www.theses.fr/1998INPL106N.

Texte intégral

Résumé :

L’utilisation de méthodes numériques classiques pour la modélisation du transfert de polluant dans un milieu poreux homogène et saturé se heurte à plusieurs difficultés. Afin de contourner ces obstacles, deux méthodes ont été choisies. La première est celle des éléments finis mixtes hybrides appliquée pour la résolution de l'équation de l'écoulement et des termes dispersif et accumulatif de l'équation de transport. La seconde est la méthode des éléments finis discontinus appliquée pour la résolution du terme convectif de l'équation de transport. Ces méthodes sont soumises à des comparaisons, faisant l'objet d'une synthèse bibliographique, avec d'autres méthodes dans le cas ou le contraste de masse volumique n'est pas pris en compte. Le code de calcul a été généralise pour le cas où la variation de la masse volumique en fonction de la concentration n'est plus négligeable. Le nouveau code est testé sur deux exemples devenus classiques. La zone de mélange lors d'un déplacement, dans un milieu poreux saturé et indéformable, de deux fluides miscibles avec contraste de masse volumique a été étudiée à travers deux types de déplacement. Le déplacement horizontal constitue le premier cas où l'inclinaison du front de polluant a été particulièrement étudiée à travers l'influence des paramètres intervenant dans ce phénomène. Le second type de déplacement est celui de l'infiltration ponctuelle d'une solution saline dans un milieu poreux initialement saturé en eau douce au repos. Une solution théorique simplifiée exprimant l'avancement du front de polluant suivant un axe vertical a été comparée aux résultats des simulations numériques.

Styles APA, Harvard, Vancouver, ISO, etc.

17

El, Soueidy Charbel-Pierre. « Éléments finis discontinus multi-domaines en temps pour la modélisation du transport en milieu poreux saturé ». Université Louis Pasteur (Strasbourg) (1971-2008), 2008. https://publication-theses.unistra.fr/public/theses_doctorat/2008/EL_SOUEIDY_Charbel-Pierre_2008.pdf.

Texte intégral

Résumé :

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
This work treats the modeling of mass transport in porous media. The advective part of the transport equation is solved using the discontinuous Galerkin (DG) finite element method. In the first part, the discretization of the spatial operator is considered and two approximation spaces are studied and compared for unstructured meshes in 2-D. The temporal discretization is tackled in the second part of this study. Two alternatives to the traditional explicit scheme are presented : a class of semi-implicit schemes and an explicit local time-stepping procedure which allows spatially varying time steps. Finally, the last part of this work consists of using the developed numerical tools to simulate density coupled flow and transport in porous media. The local time procedure is implemented in a 3-D numerical code and numerical experiments show that the model gives accurate results being highly efficient for this kind of problems

Styles APA, Harvard, Vancouver, ISO, etc.

18

Laurent, Karine. « Étude de nouveaux schémas numériques pour la simulation des écoulements à rapport de mobilités défavorable dans un contexte EOR ». Thesis, Université Paris-Saclay (ComUE), 2019. http://www.theses.fr/2019SACLC081/document.

Texte intégral

Résumé :

En simulation dynamique des réservoirs, un des artéfacts les plus gênants pour la prédiction de production est l’effet de l’orientation du maillage. Bien que celui-ci soit « normal » pour tout schéma numérique, il se trouve amplifié par l’instabilité du modèle physique, ce qui a lieu lorsque le contraste de mobilités entre l’eau (fluide poussant, utilisé dans les procédés de récupération secondaires) et l’huile (fluide poussé, contenant les hydrocarbures) dépasse un certain seuil critique. On parle alors d’écoulements à rapport de mobilités défavorable. Connu depuis longtemps, ce problème a fait l'objet de nombreux travaux dans les années 1980 ayant abouti au schéma dit à neuf points. Actuellement implanté dans PumaFlow, logiciel développé et commercialisé par IFPEN, ce schéma fonctionne relativement bien en maillages carrés et dépend d’un paramètre scalaire dont le réglage varie selon les auteurs sur la base de considérations heuristiques. Dans cette thèse, nous proposons une nouvelle démarche méthodologique afin non seulement d’ajuster ce paramètre libre de manière optimale mais aussi de généraliser le schéma aux maillages rectangulaires. La stratégie que nous préconisons repose sur une analyse d’erreur du problème, à partir de laquelle il est possible de définir une notion d’erreur angulaire et de garantir que le comportement du schéma obtenu soit le « moins anisotrope » possible via une minimisation de son écart par rapport à un comportement idéal. Cette procédure de minimisation est ensuite appliquée à deux autres familles de schémas numériques~ : (1) un schéma multidimensionnel proposé par Kozdon, dans lequel le paramètre libre est une fonction~ ; (2) un autre schéma à neuf points faisant intervenir deux paramètres scalaires. C’est ce dernier qui réduit le mieux l’effet de l’orientation lorsque le rapport des pas de maillage s’éloigne de 1. Enfin, une extension de la méthode à des modèles physiques plus complets est envisagée
In dynamic reservoir simulation, one of the most troublesome artifacts for the prediction of production is the grid orientation effect. Although this normally arises from any numerical scheme, it happens to be amplified by the instability of the physical model, which occurs when the mobility contrast between the water (pushing fluid, used in the processes of secondary recovery) and the oil (pushed fluid, containing the hydrocarbons) exceeds a some critical threshold. We then speak of flows with adverse mobility ratio. This GOE issue has received a lot of attention from the engineers. Numerous works dating back to the 1980s have resulted in the so-called nine-point scheme. Currently implemented in the IFPEN software PumaFlow, this scheme performs relatively well in square meshes and depends on a scalar parameter whose value varies from one author to another, on the grounds of heuristic considerations. In this thesis, we propose a new methodological approach in order not only to optimally adjust this free parameter, but also to extend the scheme to rectangular meshes. The strategy that we advocate is based on an error analysis of the problem, from which it is possible to define a notion of angular error and to guarantee that the behavior of the obtained scheme is the "least anisotropic" possible through a minimization of its deviation from some ideal behavior. This minimization procedure is then applied to two other families of numerical schemes: (1) a multidimensional scheme proposed by Kozdon, in which the free parameter is a function; (2) another nine-point scheme involving two scalar parameters. The latter provides the best results regarding GOE reduction when the ratio of the mesh steps is far away from 1. Finally, an extension of the method to more sophisticated physical models is envisaged

Styles APA, Harvard, Vancouver, ISO, etc.

19

Martin, Vincent. « Simulations multidomaines des écoulements en milieu poreux ». Paris 9, 2004. https://portail.bu.dauphine.fr/fileviewer/index.php?doc=2004PA090006.

Texte intégral

Résumé :

Cette thèse traite des simulations multidomaines des écoulements en milieu poreux. Trois aspects sont abordés. Une étude est d'abord menée concernant une méthode de décomposition de domaines avec des maillages non-raccordés, utilisant des conditions d'interface de type Robin, pour les éléments finis mixtes. D'autre part, la méthode est implémentée en parallèle à l'aide du système parallèle OcamlP3l, écrit par des informaticiens dans le langage Ocaml. En OcamlP3l, l'utilisateur programme et débogue en séquentiel, puis obtient le code parallèle par une simple recompilation. Enfin, nous présentons un nouveau modèle d'écoulement dans un mileu poreux contenant de grandes fractures, avec des perméabilités très grandes et/ou très faibles. Dans ce modèle, les fractures sont traitées comme des interfaces entre sous-domaines. Une analyse théorique prouve l'existence et l'unicité de la solution, fournit une estimation d'erreur, qui est confirmée par des tests numériques
This thesis is mainly concerned with the multidomain simulations of flow in porous media. Three different themes are considered. First, we study a domain decomposition method with non-matching meshes using Robin type interface conditions, for the mixed finite elements. Second, this method is implemented in parallel using the parallel system OcamlP3l, written in Ocaml by computer scientists. In OcamlP3l, the user develops and debugs sequentially, and obtains the parallel code with a mere recompilation. A realistic 3D simulation is given to validate the procedure. Finally, we present a new model for flow in a porous medium containing large fractures that may have very large and/or very small permeabilities. In this model, the fractures are treated as interfaces between subdomains. Existence and uniqueness of the solution is proved and an error estimate is obtained. Some numerical experiments show the quality of the results

Styles APA, Harvard, Vancouver, ISO, etc.

20

Boubaker, Riadh. « Modélisation thermodynamique instationnaire d'une boucle fluide diphasique à pompage capillaire pour la traction ferroviaire : étude du phénomène de changement de phase dans l'évaporateur ». Thesis, Pau, 2014. http://www.theses.fr/2014PAUU3031/document.

Texte intégral

Résumé :

Les boucles diphasiques à pompage capillaire sont des systèmes performants de transfert thermique qui peuvent transporter, de manière passive, de très importantes quantités de chaleur sur de grandes distances : utilisant les phénomènes qui apparaissent lors de la vaporisation d'un liquide dans un corps poreux (la mèche) contenu dans l'évaporateur, ces dispositifs offrent l'avantage de ne pas utiliser d'organe mécanique de puissance pour mettre lefluide frigorigène en mouvement. Vu les performances dont ces systèmes ont fait preuve dans le domaine spatial, leur utilisation dans le domaine gravitaire est aujourd'hui sérieusement étudiée, en particulier dans le domaine ferroviaire. Le travail de cette thèse porte sur la description du comportement thermodynamique instationnaire global d'une boucle diphasique à pompage capillaire, utilisée par Alstom Transport pour refroidir ses composants d'électronique de puissance. La première partie de ce travail consiste en la modélisation du transfert de masse et de chaleur dans la mèche contenu dans l'évaporateur, qui est le composant clé de la boucle. Un modèle mathématique instationnaire 2D a été développé pour décrire l'écoulement diphasique avec changement de phase dans le milieu poreux. Les résultats numériques de ce modèle montrent la formation d'une poche de vapeur au sein de la mèche poreuse. Une étude approfondie est élaborée pour décrire la dynamique de la croissance de cette poche en fonction de plusieurs paramètres (géométrie, sous-refroidissement, température de saturation, porosité). La deuxième partie de ce travail consiste à coupler le modèle de l'évaporateur au reste de la boucle afin d'obtenir un modèle thermodynamique global avec une interface liquide/vapeur mobile dans la mèche poreuse. Ce modèle est validé en régime transitoire par une confrontation de ses résultats aux mesures obtenues lors d'une précédente campagne expérimentale. Le comportement de la boucle au cours de la phase de démarrage est ensuite étudié. Finalement, la réponse transitoire de la boucle globale est analysée, en portant un intérêt particulier au comportement de l'évaporateur soumis à un créneau de puissance appliqué, à la variation de la température de consigne du réservoir de contrôle et à la variation de la température de la source froide en contact avec le condenseur
Capillary pumped loop is a two phase device that uses the phase change phenomena occurred in the porous wick to transport large heat loads over long distances (several meters) without the need of any mechanical pump. CPLs have been developed and successfully employed for the thermal control of satellites. Thanks to their heat transport capacity, their use is now seriously considered in the gravity field especially for cooling electronic devices of a railroad traction chain. The work of this thesis focuses on the description of the unsteady thermodynamic behavior of a capillary pumped loop used by Alstom Transport to cool its power electronics components. The first part of the thesis consists to study the key component of CPL: the capillary evaporator. A 2D unsteady mathematical model has been developed to describe the heat and mass transfer inside the porous wick of the CPL evaporator. The numerical results of this model show the formation of a vapor pocket inside the porous wick. The influences of evaporator geometry, liquid subcooling, saturation temperature and wick porosity on the dynamic growth of the vapor pocket are discussed in detail. The second part of this work consists of coupling the CPL evaporator model with other CPL components in order to obtain a global capillary pumped poop model with a mobile liquid/vapor interface in the porous wick. The proposed model is validated by comparing numerical simulations with experimental results obtained in a previous works realized at Alstom Transport. The CPL startup is then studied. Finally, the influences of the applied power, reservoir temperature and heat sink temperature on the dynamic response of the overall loop are analyzed, with a particular interest to the evaporator behavior

Styles APA, Harvard, Vancouver, ISO, etc.

21

Deville, Lilas. « Modélisation et caractérisation expérimentale des joints à brosse ». Thesis, Poitiers, 2016. http://www.theses.fr/2016POIT2326/document.

Texte intégral

Résumé :

Les joints à brosse sont des joints compliants utilisés dans les turbomachines. Ils sont constitués d'une brosse de fils très fins compactés entre deux plaques, de telle manière que l'ensemble comble l'écart séparant le stator du rotor. Les fils sont le plus souvent constitués d'un alliage à base de cobalt appelé Haynes 25 (d'autres matériaux peuvent être utilisés, tels que le Kevlar). Le comportement de ces joints présente un fort couplage entre l'écoulement du fluide et la déformation de la brosse. Ce type de joints peut être modélisé en considérant la brosse comme un milieu poreux. Cette méthode présente l'avantage de permettre un calcul relativement simple de l'écoulement, mais elle demande l'usage de données expérimentales afin de calibrer ses perméabilités (i.e. sa capacité à laisser un fluide s'écouler à travers lui). L'objectif de cette thèse est de proposer une modélisation par milieu poreux indépendante des données expérimentales. Une simulation itérative entre le calcul de la déformation des fils sous l'action des efforts de pression et la résolution de l'écoulement dans le milieu poreux formé par la brosse déformée a été mise au point. Elle repose sur le calcul numérique des perméabilités obtenues à partir de la perte de charge générée par un réseau de cylindres représentatifs de la brosse.Des joints de faible diamètre en Haynes 25 ont été testés avec de l'air afin de vérifier la validité des résultats donnés par la simulation. L'évolution du débit avec la différence de pression et avec l'interférence des joints avec le rotor a été étudiée. Les déplacements du rotor et la température des joints ont été suivis durant les essais
Brush seals are compliant dynamic seals that are used in the turbomachines. They are made of a brush of fine bristles sandwiched between two plates, in such a way that the assembly fills the gap between the rotor and the stator. The bristles are usually made from a cobalt based alloy called Haynes 25 (other materials can be used, such as Kevlar). The behaviour of these seals shows a strong coupling between the fluid flow and the brush deformation.This kind of seal can be modelled by considering the brush as a porous medium. This method presents the benefit of a rather simple calculation of the flow, but it requires the use of experimental data in order to calibrate its permeabilities (i.e. its ability to let the fluid flow through it). The aim of this thesis is to propose a model using the porous medium approach but without the need of preliminary experimental data. An iterative simulation between the calculation of the bristles' deformation under the influence of the pressure forces and the solving of the flow in the porous medium formed by the deformed brush was developed. It relies on the numerical calculation of the permeabilities obtained from the pressure drop generated by an array of cylinders representative of the brush.Low diameter Haynes 25 seals were tested with air in order to verify the validity of the results given by the simulation. The evolution of the mass flow rate with the pressure difference and the interference between the seal and the rotor was studied. Displacements of the rotor and temperature of the seals were also monitored during the tests

Styles APA, Harvard, Vancouver, ISO, etc.

22

Dugois, Kévin. « Simulation à l’échelle microscopique et analyse macroscopique de l’imprégnation d’un matériau composite par un fluide chargé en particules ». Thesis, Bordeaux, 2017. http://www.theses.fr/2017BORD0011/document.

Texte intégral

Résumé :

Dans le but d’améliorer le comportement thermo-mécanique des aubes de turbine présentes dans les moteurs d’avion développés par SAFRAN, il est nécessaire de mettre au point un nouveau matériau composite. Le procédé de fabrication de ce matériau est complexe et requiert une densification par voie liquide divisée en deux étapes. Cette thèse s’intéresse à la modélisation numérique de la première étape appelée Slurry Cast/APS.Celle-ci consiste en l’injection et le confinement, dans la préforme fibreuse, de particules préalablement mises en suspension. Pour cela, nous avons développé à l’échelle des fibres,un modèle qui utilise les équations de Navier-Stokes incompressibles et monophasiques ,l’équation de Phillips [Phillips et al., 1992] et une loi rhéologique [Krieger, 1972]. Après validation des résultats numériques par comparaison avec des résultats expérimentaux [Hampton et al., 1997] et théoriques [Belfort et al., 1994], le modèle est utilisé pour simuler l’écoulement autour de géométries de tissage proches du matériau étudié
In order to improve thermo-mechanical behavior of tubine blades in SAFRAN engines plane, a new composite material is necessary. The manufacturing process to obtain this composite is intricate and requires a two steps fluid densification process. This thesis focuses on numerical simulation of the first one called Slurry Cats/APS. In this step, suspended particles are introducted and captured in the reinforcement. For that purpose,we carry out a model at fiber scale, using Navier-Stokes equations in incompressible and monophasic formulation, Phillips equations [Phillips et al., 1992] and a rheological law [Krieger, 1972]. After validation step consisting in a comparison of computational results with experiments [Hampton et al., 1997] and theorical law [Belfort et al., 1994], this model has been used to simulate flow around geometries similar to those encountered in our composite material

Styles APA, Harvard, Vancouver, ISO, etc.

23

Szydlarski, Mikolaj. « Algebraic Domain Decomposition Methods for Darcy flow in heterogeneous media ». Phd thesis, Université Pierre et Marie Curie - Paris VI, 2010. http://tel.archives-ouvertes.fr/tel-00550728.

Texte intégral

Résumé :

Afin de répondre aux besoins de l'industrie pétrolière d'une description plus fine de la géométrie et des propriétés pétrophysiques des bassins et des réservoirs, la simulation numérique des écoulements en milieux poreux doit évoluer vers des algorithmes plus performants et plus robustes vis à vis de la taille des simulations, de la complexité des maillages et des hétérogénéités du milieu poreux. Les méthodes de décomposition de domaine constituent une alternative aux méthodes multigrilles et pourraient permettre de lever les difficultés précédentes en terme de robustesse et d'efficacité sur architectures parallèles. Elles sont par nature plus adaptées au calcul parallèle et sont plus robustes en particulier lorsque les sous domaines sont résolus par des méthodes directes. Elles permettent aussi de traiter dans un cadre unique les couplages de modèles comme les puits ou les failles conductrices et s'étendent au cas des systèmes couplés. Le travail de thèse traite plus particulièrement de méthodes définies au niveau algébrique. On ne suppose pas avoir une connaissance préalable du problème continu dont la matrice provient. On n'a pas non plus accés aux matrices avant assemblage. Ce manque d'informations a priori rend plus difficile la construction de méthodes efficaces. On propose deux nouvelles méthodes de construction de méthodes de décomposition de domaine au niveau algébrique: la construction de conditions d'interface optimisées et d'une grille grossière. Ce dernier point est particulièrement important pour avoir des méthodes robustes vis à vis du nombre des sous-domaines. Les méthodes sont adaptatives et basées sur l'analyse de l'espace de Krylov généré durant les premières itérations de la méthode de Schwarz classique. A partir des vecteurs de Ritz correspondant aux plus basses valeurs propres, on construit des conditions d'interface et des grilles grossières qui annihilent l'erreur sur ces composantes. Les méthodes ont été testées sur des calculateurs parallèles pour des matrices issues de la simulation de milieux poreux.

Styles APA, Harvard, Vancouver, ISO, etc.

24

Pauthenet, Martin. « Macroscopic model and numerical simulation of elastic canopy flows ». Thesis, Toulouse, INPT, 2018. http://www.theses.fr/2018INPT0072/document.

Texte intégral

Résumé :

On étudie l'écoulement turbulent d'un fluide sur une canopée, que l'on modélise comme un milieu poreux déformable. Ce milieu poreux est en fait composé d'un tapis de fibres susceptibles de se courber sous la charge hydrodynamique du fluide, et ainsi de créer un couplage fluide-structure à l'échelle d'une hauteur de fibre (honami). L'objectif de la thèse est de développer un modèle macroscopique de cette interaction fluide-structure, afin d'en réaliser des simulations numériques. Une approche numérique de simulation aux grandes échelles est donc mise en place pour capturer les grandes structures de l'écoulement et leur couplage avec les déformations du milieu poreux. Pour cela nous dérivons les équations régissant la grande échelle, au point de vue du fluide ainsi que de la phase solide. À cause du caractère non-local de la phase solide, une approche hybride est proposée. La phase fluide est décrite d'un point de vue Eulerien, tandis que la description de la dynamique de la phase solide nécessite une représentation Lagrangienne. L'interface entre le fluide et le milieu poreux est traitée de manière continue. Cette approche de l'interface fluide/poreux est justifiée par un développement théorique sous forme de bilan de masse et de quantité de mouvement à l'interface. Ce modèle hybride est implémenté dans un solveur écrit en C$++$, à partir d'un solveur fluide disponible dans la librairie CFD \openfoam. Un préalable nécessaire à la réalisation d'un tel modèle macroscopique est la connaissance des phénomènes de la petite échelle en vue de les modéliser. Deux axes sont explorés concernant cet aspect. Le premier consiste à étudier les effets de l'inertie sur la perte de charge en milieu poreux. Un paramètre géométrique est proposé pour caractériser la sensibilité d'une microstructure poreuse à l'inertie de l'écoulement du fluide dans ses pores. L'efficacité de ce paramètre géométrique est validée sur une diversité de microstructures et le caractère général du paramètre est démontré. Une loi asymptotique est ensuite proposée pour modéliser les effets de l'inertie sur la perte de charge, et comprendre comment celle-ci évolue en fonction de la nature de la microstructure du milieu poreux. Le deuxième axe d'étude de la petite échelle consiste à étudier l'effet de l’interaction fluide-structure à l'échelle du pore sur la perte de charge au niveau macroscopique. Comme les cas présentent de grands déplacements de la phase solide, une approche par frontières immergées est proposée. Ainsi deux méthodes numériques sont employées pour appliquer la condition de non-glissement à l'interface fluid/solide: l'une par interface diffuse, l'autre par reconstitution de l'interface. Cela permet une validation croisée des résultats et d'atteindre des temps de calcul acceptables tout en maîtrisant la précision des résultats numériques. Cette étude permet de montrer que l'interaction fluide-structure à l'échelle du pore a un effet considérable sur la perte de charge effective au niveau macroscopique. Des questions fondamentales sont ensuite abordées, telles que la taille d'un élément représentatif ou la forme des équations de transport dans un milieu poreux souple
We study the turbulent flow of a fluid over a canopy, that we model as a deformable porous medium. This porous medium is more precisely a carpet of fibres that bend under the hydrodynamic load, hence initiating a fluid-structure coupling at the scale of a fibre's height (honami). The objective of the thesis is to develop a macroscopic model of this fluid-structure interaction in order to perform numerical simulations of this process. The volume averaging method is implemented to describe the large scales of the flow and their interaction with the deformable porous medium. An hybrid approach is followed due to the non-local nature of the solid phase; While the large scales of the flow are described within an Eulerian frame by applying the method of volume averaging, a Lagrangian approach is proposed to describe the ensemble of fibres. The interface between the free-flow and the porous medium is handle with a One-Domain- Approach, which we justify with the theoretical development of a mass- and momentum- balance at the fluid/porous interface. This hybrid model is then implemented in a parallel code written in C$++$, based on a fluid- solver available from the \openfoam CFD toolbox. Some preliminary results show the ability of this approach to simulate a honami within a reasonable computational cost. Prior to implementing a macroscopic model, insight into the small-scale is required. Two specific aspects of the small-scale are therefore studied in details; The first development deals with the inertial deviation from Darcy's law. A geometrical parameter is proposed to describe the effect of inertia on Darcy's law, depending on the shape of the microstructure of the porous medium. This topological parameter is shown to efficiently characterize inertia effects on a diversity of tested microstructures. An asymptotic filtration law is then derived from the closure problem arising from the volume averaging method, proposing a new framework to understand the relationship between the effect of inertia on the macroscopic fluid-solid force and the topology of the microstructure of the porous medium. A second research axis is then investigated. As we deal with a deformable porous medium, we study the effect of the pore-scale fluid-structure interaction on the filtration law as the flow within the pores is unsteady, inducing time-dependent fluidstresses on the solid- phase. For that purpose, we implement pore-scale numerical simulations of unsteady flows within deformable pores, focusing for this preliminary study on a model porous medium. Owing to the large displacements of the solid phase, an immersed boundary approach is implemented. Two different numerical methods are compared to apply the no-slip condition at the fluid-solid interface: a diffuse interface approach and a sharp interface approach. The objective is to find the proper method to afford acceptable computational time and a good reliability of the results. The comparison allows a cross-validation of the numerical results, as the two methods compare well for our cases. This numerical campaign shows that the pore-scale deformation has a significant impact on the pressure drop at the macroscopic scale. Some fundamental issues are then discussed, such as the size of a representative computational domain or the form of macroscopic equations to describe the momentum transport within a soft deformable porous medium

Styles APA, Harvard, Vancouver, ISO, etc.

25

McKee, François. « Etude et mise à l'échelle des écoulements diphasiques en milieux poreux hétérogènes par une approche d'optimisation ». Nantes, 2014. http://archive.bu.univ-nantes.fr/pollux/show.action?id=bd435ce7-5514-401b-9305-92faceb0a7ba.

Texte intégral

Résumé :

Le couplage de l'expertise géologique et de la simulation numérique permet à l'industrie pétrolière de construire des modèles de réservoir toujours plus détaillés. Ces modèles intègrent la totalité des données de production et de sismique disponibles, mais nécessitent l'utilisation de la géostatistique et d'une approche stochastique. Un grand nombre de simulations est requis pour estimer les réserves d'hydrocarbures et pour optimiser la récupération de pétrole. Lors d'une première étape, les ingénieurs conçoivent un modèle géologique respectant la géométrie du réservoir et pouvant contenir plusieurs millions de mailles. Dans un second temps, un modèle de réservoir, plus grossier, est construit. Celui-ci contient un nombre réduit de mailles et assure un temps de simulation raisonnable. L'upscaling est la caractérisation du modèle de réservoir conformément aux données du modèle géologique. L'upscaling d'écoulements multiphasiques est encore un problème d'actualité. En effet, les méthodes disponibles sont souvent appliquées au cas par cas ou par étapes successives. Les courbes de perméabilités relatives, caractéristiques des écoulements multiphasiques, sont au coeur du problème. Nous proposons ici une méthode d'upscaling des courbes de perméabilités relatives par une approche d'optimisation. Le réservoir étudié est découpé en zones distinctes. Les courbes de perméabilités relatives upscalées, obtenues pour chaque zone, permettent d'obtenir un modèle de réservoir approchant le modèle géologique en termes d'écoulement équivalent
Geological and geophysical expertise coupled to numerical simulation allow the petroleum industry to build increasingly detailed reservoir models. These models integrate the whole set of available data (production and seismic) but involve geostatistics and stochastic approach. A large number of simulations is required to estimate hydrocarbon reserves and optimize oil recovery. During a first step, engineers build a geological model respecting the real geometry of the reservoir and containing possibly million of cells. The second step consist in building a coarser model, the reservoir model. It contains less cells in order to ensure that numerical simulation is feasible within a reasonable time. Upscaling is the characterization of the reservoir model accordingly to the geological model properties. Multiphase flows upscaling is still an actual issue. Indeed, methods are applied from case to case or in successive stages. Relative permeability curves, characteristic of multiphase flow equations, are a real issue. The main purpose here is to propose an upscaling method for relative permeability by an optimization approach. The studied reservoir is divided in several zones. The upscaled relative permeability curves obtained for each zone allow the building of a reservoir model. This reservoir model forms a good match for the geological model in terms of equivalent flow

Styles APA, Harvard, Vancouver, ISO, etc.

26

Gueye, Abdoulaye. « Modélisation et simulations numériques des écoulements et instabilités thermiques de fluides non-Newtonien en milieu poreux ». Thesis, Lille 1, 2015. http://www.theses.fr/2015LIL10154/document.

Texte intégral

Résumé :

Ce travail de thèse sur les milieux poreux est axé sur deux parties. La première concerne l'étude numérique de l'écoulement d'un fluide Newtonien ou non-Newtonien au sein d'un système fluide/poreux. L'approche à un seul domaine qui consiste à écrire l'équation de Navier-Stokes incluant le terme de Darcy-Brinkman-Forchheimer est adoptée dans cette étude. La relation entre le gradient de pression et la vitesse débitante linéaire dans le cas de Darcy où le fluide est Newtonien, est obtenue. Cette relation est étendue dans le cas non-Darcy où le fluide est non Newtonien. L'influence des nombres de Darcy et de Forchheimer sur la structure de l'écoulement est montrée. Dans la seconde partie, une étude de stabilité linéaire et numérique de la convection naturelle de fluides viscoélastiques saturant une couche poreuse horizontale chauffée par un flux constant est réalisée. Une étude d'instabilité primaire et secondaire nous a permis de montrer que pour un fluide Newtonien, la convection monocellulaire perd sa stabilité au profit des rouleaux longitudinaux. Dans le cas des fluides viscoélastiques, on trouve que l'élasticité du fluide induit la sélection des rouleaux transversaux propagatifs. Une solution numérique basée sur un schéma aux différences finies est venue conforter ces résultats analytiques
The present thesis on porous media concentrates in two parts. The first concerns the numerical study of the flow of a Newtonian or a non-Newtonian fluid within a fluid/porous system. The approach of a single domain, which consists of/in writing the Navier-Stokes equation including the Darcy-Brinkham-Forchheimer term, is chosen in this study. The linear relation between the pressure gradient and the bulk velocity in the Darcy case, for which the fluid is Newtonian, is obtained. This relation is extended to the non-Darcy case, for which the fluid is non-Newtonian. The influence of Darcy and Forchheimer numbers on the structure of the flow is presented. In the second part, linear stability and numerical analysis of the natural convection of viscoelastic fluids saturating a horizontal porous layer heated by a constant flux is performed. A primary and secondary instability study allowed to show that, for a Newtonian fluid, the unicellular convection loses its stability to the benefit of longitudinal rolls. In the case of viscoelastic fluids, the elasticity of the fluid leads to the selection of propagation transverse rolls. A numerical solution based on a finite difference scheme has reinforced these analytical results

Styles APA, Harvard, Vancouver, ISO, etc.

27

Mathelin, Lionel. « Etude numérique, expérimentale et théorique de l'influence de l'effusion sur la dynamique et la thermique de l'écoulement autour d'un cylindre circulaire poreux ». Lyon, INSA, 2001. http://theses.insa-lyon.fr/publication/2001ISAL0028/these.pdf.

Texte intégral

Résumé :

Ce travail concerne l'étude numérique et expérimentale de l'écoulement et des transferts de chaleur autour d'un cylindre circulaire poreux en écoulement croisé, soumis à de l'effusion non-isotherme (écoulement interne). Les régimes d'écoulement sont stationnaires et instationnaires, laminaires et turbulents. Deux modèles d'effusion sont proposés. Les manipulations expérimentales permettent de valider les simulations numériques et les modèles d'effusion. En particulier, l'anémométrie à fils chauds et fil froid permet de comparer les profils de vitesse et de température avec les simulations numériques. Un bon accord est constaté. Un couplage écoulement interne à la matrice poreuse / écoulement externe au cylindre est réalisé sous la forme d'une fonction externe au code numérique, prenant en compte les transferts radiatifs afin de déterminer la température de paroi à l'équilibre. Les performances de l'injection de liquide (transpiration) sont expérimentalement évaluées avec de l'eau et de l'éthanol et montrent une excellente efficacité du procédé. Enfin, l'influence de l'effusion sur les instabilités est étudiée. Il est constaté que l'effusion fait baisser la fréquence de détachement des tourbillons primaires et secondaires et modifie la structure de la couche de cisaillement. En présence d'effusion, une série de pics supplémentaires apparaît. En utilisant la théorie de stabilité linéaire, il est montré que les taux de croissance sont significativement affectés par l'effusion et que la structure du sillage proche, et notamment la zone absolument instable, est modifiée. Enfin, des simulations numériques 3-D ont permis de conclure que l'effusion réduit les phénomènes tridimensionnels
This work is concerned about the numerical and experimental study of the flow and heat transfer over a porous circular cylinder in cross-flow submitted to non-isothermal blowing (internal flow). Flow regimes are steady and unsteady, laminar and turbulent. Two blowing models are proposed. The experiments allow for numerical simulation and blowing models validation. In particular, using hot- and cold-wire anemometry, velocity and temperature profiles are compared with numerical simulations. A coupling between internal flow within the porous matrix and external flow over the cylinder is achieved, in the form of an user defined function of the code accounting for radiative heat transfer, to determine the equilibrium wall temperature. Performances of liquid injection (transpiration cooling) are experimentally evaluated with water and ethanol and prove excellent thermal protection effectiveness of the process. Moreover, the blowing influence on instabilities influence on instabilities is studied. It is seen that the blowing lowers the primary and secondary vortex shedding frequency and modifies the shear layer structure. In case of blowing, a series of peaks appears. Using the linear stability theory, it is shown that growth rates are significantly affected by blowing and that the near wake structure, and the absolutely convective zone in particular, is modified. Finally, performed 3-D numerical simulations lead to the conclusion that blowing reduces tri-dimensional effects

Styles APA, Harvard, Vancouver, ISO, etc.

28

Arcen, Boris. « Etude par simulation numérique directe du comportement et de la dispersion de particules solides en écoulement non homogène isotherme ou anisotherme ». Nancy 1, 2006. http://docnum.univ-lorraine.fr/public/SCD_T_2006_0163_ARCEN.pdf.

Texte intégral

Résumé :

Le travail présenté dans ce mémoire concerne l'étude du mouvement d'inclusions solides en suspension dans un écoulement turbulent de canal isotherme et anisotherme par simulation numérique directe. Grâce à cet outil de simulation, nous avons pu analyser l'influence de l'inertie et du croisement de trajectoires sur les caractéristiques dynamiques et thermiques de la phase dispersée ainsi que sur celles du fluide vu au sein d'une turbulence non homogène. Cela pourra notamment faciliter le développement futur des modélisations euléro-lagrangienne et euléro-eulérienne dans ce type de turbulence. Nous avons essayé d'examiner au mieux la conséquence de ces effets sur les statistiques de la phase dispersée telles que la concentration, la moyenne et l'écart type de la vitesse, les covariances fluide-particules, les corrélations triples de la vitesse des particules. Parallèlement à cela, les caractéristiques du fluide vu par les particules ont été étudiées, nous nous sommes intéressés à la vitesse de dérive, aux tensions de Reynolds du fluide vu, et la décorrélation temporelle des fluctuations de la vitesse du fluide vu. En ce qui concerne la partie thermique, nous présentons les statistiques thermiques de la phase dispersée et du fluide vu par les particules au sein de l'écoulement anisotherme vertical descendant. Tous ces aspects sont développés en gardant à l'esprit le cadre général de cette étude, c'est-à-dire comprendre le comportement thermique de la phase dispersée et fournir des informations concernant des grandeurs importantes intervenant dans la modélisation de tels écoulements
This dissertation is devoted to the investigation of the motion of solid particles in an isothermal and non isothermal inhomogeneous turbulent gas flow using direct numerical simulation. By means of this numerical tool, it was possible to analyse the inertia and crossing trajectory effects on the thermal and dynamical characteristics of the dispersed phase as well as on those of the fluid seen. A better understanding of the influence of these two effects is an important issue in the development of Eulerian-Eulerian and Eulerian-Lagrangian models. The study was conducted by extracting dispersed phase statistics such as concentration, mean and root mean square of the particle velocity, fluid-particle covariances, and triple particle velocity correlations. The characteristics of the fluid seen which were examined are the drift velocity, the Reynolds stresses, and the decorrelation time scales of the fluctuating velocity of the fluid seen. Concerning non isothermal flows, we present the thermal statistics of the dispersed phase as well as those of the fluid seen obtained in a downward channel flow. These aspects are investigated keeping in mind the objective of the present study, i. E. Understand the thermal behaviour of the dispersed phase and provide information about crucial parameters which appear in the modelling of such flows

Styles APA, Harvard, Vancouver, ISO, etc.

29

Cohen, Charles-Édouard. « Modélisation et simulation de la stimulation acide des puits carbonatés ». Toulouse, INPT, 2007. http://ethesis.inp-toulouse.fr/archive/00000601/.

Texte intégral

Résumé :

L'injection d'acide dans un puits au sein de roches carbonatées permet l'amélioration de sa productivité, en particulier si celui-ci est endommagé. Cette opération s'appuie sur un mécanisme instable qui sous certaines conditions produit de longues percées appelées "wormholes". L'objectif pour l'opérateur est d'utiliser cette instabilité afin de reconnecter le puits au réservoir en minimisant le temps du traitement et la quantité d'acide nécessaire. Les travaux présentés dans cette thèse ont pour objectif d'améliorer la compréhension des mécanismes de l'instabilité et optimiser le traitement acide. Dans une première partie nous réalisons une série d'études de stabilité linéaire du front de dissolution d'un milieu poreux. Les résultats expliquent les rôles des nombres adimensionnels dans les transitions entre les différents régimes de dissolution et permettent de prédire ces transitions. Nous utilisons ensuite des simulations d'acidification de carottes sur les figures de dissolution dans des géométries 2D et 3D, cartésiennes et radiales. Dans une dernière partie nous proposons une modélisation à grande échelle de la dissolution basée sur une approche double milieux. Ce modèle est présenté avec un exemple de simulation de traitement acide d'un puits. Enfin, ce simulateur est inclus dans une boucle d'inversion permettant d'optimiser les paramètres du traitement en fonction de la production d'huile du puits
Matrix acidizing of a well in a carbonated reservoir improves its productivity, especially when the well is damaged. This stimulation technique relies on an unstable process which creates, under specific conditions, empty channels called "wormholes". The objective of the operator is to use this instability to reconnect the well to the reservoir while minimizing the duration of the treatment and the volume of acid injected. The aim of the present thesis is to improve the understanding of the mechanisms of instability and to optimize the acid treatment. In the first part, we present a series of linear stability analysis of the dissolution front in a porous media. Results explain the role of dimensionless numbers in transitions between the different dissolution regimes and serve to predict these transitions. We also use simulations of core acidizing to study the wormhole density, the mechanism of wormhole competition and the effect of confinement induced by the core boundaries on dissolution patterns in 2D and 3D, in linear and radial geometries. Then, in the second part, we propose a large scale model to simulate full acid treatments, based on a dual porosity approach. An exemple of simulation of a full acid treatment illustrated the model. Finally, we link this acidizing simulator to a reservoir simulator and an inversion software in order to optimise the treatment parameters, as a function of the oil production

Styles APA, Harvard, Vancouver, ISO, etc.

30

Coelho, Daniel. « Génération, géométrie et propriétés de transport de milieux granulaires ». Poitiers, 1996. http://www.theses.fr/1996POIT2346.

Texte intégral

Résumé :

Ce travail porte sur les proprietes geometriques et de transport d'une classe de milieux poreux constitues par l'empilement ou l'agregation de particules. Une procedure de simulation numerique tridimensionnelle de la deposition sequentielle de grains de forme quelconque a ete developpee. Elle prend en compte la translation et la rotation des particules durant leur chute. Une etude systematique des proprietes geometriques des empilements qui en resultent a ete conduite. Les proprietes de transport de ces lits de grains ont egalement ete caracterisees et reliees a leurs caracteristiques geometriques. On a evalue les tenseurs de permeabilite, de conductivite, de dispersivite, et dans le cas ou la surface des particules est chargee electriquement, les tenseurs decrivant les effets electroosmotiques. Les outils numeriques developpes pour ce dernier point s'appliquent quand les doubles couches electriques ne sont pas d'epaisseur negligeable, c'est-a-dire hors du domaine ou des simplifications permettent un traitement analytique. Il ressort en particulier que la plupart des coefficients de transport de ces materiaux granulaires peuvent etre relies directement a quelques parametres geometriques globaux. Une nouvelle relation directe entre le coefficient electroosmotique et la permeabilite a egalement ete mise en evidence. Finalement, des agregats fractals de particules ont ete etudies suivant la meme approche

Styles APA, Harvard, Vancouver, ISO, etc.

31

Trottier, Nicolas. « Modélisation des écoulement en milieux poreux fracturés : estimation des paramètres par approche inverse multi-échelle ». Phd thesis, Université de Strasbourg, 2014. http://tel.archives-ouvertes.fr/tel-01037933.

Texte intégral

Résumé :

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.

Styles APA, Harvard, Vancouver, ISO, etc.

32

Ben, amara Sami. « Ecoulements et transferts thermiques en convection naturelle dans les milieux macro-poreux alimentaires application aux refrigerateurs menagers ». Phd thesis, INAPG (AgroParisTech), 2005. http://pastel.archives-ouvertes.fr/pastel-00001581.

Texte intégral

Résumé :

L'étude porte sur les écoulements et les transferts en convection naturelle au sein d'empilements de produits alimentaires et plus particulièrement dans des réfrigérateurs ménagers. L'approche est à la fois expérimentale et numérique. Une première étude a été effectuée dans un empilement ordonné de sphères soumis à une convection à faible vitesse d'air (u<0.2 ms-1). L'intensité des échanges convectifs entre les objets et l'air ainsi qu'entre la paroi de la cavité et l'air a été caractérisée en fonction de la position et de la vitesse. Par ailleurs, les résultats ont montré que les autres modes de transfert (conduction et rayonnement entre les objets et avec les parois) sont du même ordre de grandeur que la convection. Des mesures de vitesse d'air (par vélocimétrie par images de particules) et de température ont été effectuées dans une enceinte fermée représentant une maquette d'un réfrigérateur ménager (à l'échelle 1) dans différentes configurations (maquette vide, maquette chargée d'empilement de produits inertes thermiquement et de produits échangeant de la chaleur). Des instationnarités dans le bas de l'enceinte vide ont été mis en évidence (Ra" 3´108). La présence d'empilements d'objets (même inerte thermiquement) modifie de façon importante les écoulements et les transferts. Les simulations numériques effectuées avec un logiciel de mécanique de fluide pour les différentes configurations étudiées sur la maquette et dans le cas d'un réfrigérateur réel ont donné des prédictions satisfaisantes en comparaison avec les mesures expérimentales. Mots clés: milieu macro-poreux, convection naturelle, coefficient de transfert convectif, conduction, rayonnement, réfrigérateur ménager, simulations numérique (CFD).

Styles APA, Harvard, Vancouver, ISO, etc.

33

Barbosa, Mota José Paulo. « Modélisation des transferts couplés en milieux poreux : étude du stockage de gaz naturel par adsorption et de la convection naturelle en espace poreux annulaire ». Vandoeuvre-les-Nancy, INPL, 1995. http://docnum.univ-lorraine.fr/public/INPL_T_1995_BARBOSA_MOTA_J_P.pdf.

Texte intégral

Résumé :

La première partie de la thèse est consacrée à l'étude du stockage de gaz naturel par adsorption. Plusieurs outils théoriques ont été développés, notamment un simulateur bidimensionnel d'un réservoir de stockage permettant d'étudier les phénomènes d'adsorption, l'hydrodynamique et les effets thermiques associés aux remplissages et vidanges des réservoirs avec des géométries et dimensions réalistes. Les simulations ont mis en évidence des variations importantes de la température du stockage lors des remplissages et des décharges rapides. Ces variations de température entrainent des pertes de performance dans le système en réduisant la capacité de stockage. Plusieurs possibilités d'atténuer les changements de température sont discutées. L’influence des résistances diffusionnelles dans le comportement du stockage ainsi que le développement d'une formule pour estimer les temps minimales de remplissage, sont aussi étudiées. Un montage expérimental simple donnant accès à des paramètres physico-chimiques caractéristiques d'un système gaz/adsorbant est proposé. La deuxième partie de la thèse concerne l'étude de la convection naturelle dans un milieu poreux entre deux cylindres isothermes et horizontaux. Pour une couche concentrique et pour des rapports des rayons supérieurs à 1,7, la relation entre le nombre de Nusselt global et le nombre de Rayleigh présente une boucle fermée d'hystérésis associée à la transition entre les régimes stationnaires à 2 et 4 cellules. Pour les rapports des rayons plus petits, la boucle d'hystérésis disparait et des régimes stationnaires avec 2, 4, 6 et 8 cellules sont progressivement obtenus dans la couche à mesure que le nombre de Rayleigh augmente. Pour des nombres de Rayleigh suffisamment élevés, le déplacement judicieux du centre du cylindre interne vers le haut, afin que son centre soit au-dessus du centre du cylindre externe, diminue les pertes thermiques d'environ 10% par rapport au cas concentrique

Styles APA, Harvard, Vancouver, ISO, etc.

34

Arson, Chloé. « Étude théorique et numérique de l'endommagement thermo-hydro-mécanique des milieux poreux non saturés ». Phd thesis, Ecole des Ponts ParisTech, 2009. http://pastel.archives-ouvertes.fr/pastel-00005608.

Texte intégral

Résumé :

Les déchets radioactifs sont stockés dans des milieux poreux multi-phasiques. Un nouveau modèle d'endommagement, formulé en variable indépendantes (contrainte nette, succion et contrainte thermique) est proposé pour ces géomatériaux. La variable d'endommagement est un tenseur d'ordre deux, dont les valeurs principales croissent avec la tension. Les relations contrainte/déformation sont dérivées de l'énergie libre, dont l'expression est postulée. Les rigidités dégradées sont calculées en appliquant le principe de l'Energie Elastique Equivalente pour chaque variable d'état. La fissuration est prise en compte dans les transferts d'humidité, au moyen de longueurs internes. Le modèle d'endommagement a été intégré dans le code d'Eléments Finis Θ-Stock. Le modèle mécanique a été validé en comparant des résultats numériques à des résultats expérimentaux et théoriques. Les tendances données par le modèle dans les études paramétriques sur des configurations réalistes complexes sont bonnes.

Styles APA, Harvard, Vancouver, ISO, etc.

35

Beaudoin, Anthony. « Contribution numérique à l'étude du transport de masse dans les milieux poreux saturés ». Le Havre, 2004. http://www.theses.fr/2004LEHA0002.

Texte intégral

Résumé :

L'objectif de la thèse est la résolution des problèmes de transport dans les milieux poreux à l'aide des méthodes particulaires. Deux problèmes sont résolus avec les méthodes particulaires : la migration d'un soluté dans une nappe phréatique et l'écoulement de l'eau dans un milieu poreux à l'échelle des pores. Pour résoudre l'équation de convection-dispersion dans le cas général d'une dispersion anisotrope, un algorithme basé sur les méthodes particulaires et la méthode de vitesse de dispersion est présenté. Pour simuler l'écoulement dans un milieu poreux à des faibles nombres de Reynolds, un algorithme basé sur une méthode particulaire anisotrope et une méthode intégrale de frontière est présenté. La méthode particulaire anisotrope permet d'améliorer la description de la couche limite se formant autour des grains du sol. La méthode intégrale de frontière permet de satisfaire la condition d'adhérence en écrivant une équation intégrale de frontière sur le contour des grains du sol
The aim of the thesis is to solve transport problems by means of particles methods. Two problems are resolved with particles methods : the migration of a solute in the groundwater and the water flow in a porous medium at the pore scale. In order to solve the advection-dispersion equation in the general case of a non-isotropic dispersion, a numerical model based on the particles methods and the dispersion velocity method is presented. In order to simulate flows in porous media at low Reynolds numbers, a numerical model based on an anisotropic particles method and a boundary integral method is presented. The anisotropic particles method allows to improve the description of the boundary layer appearing around grains of the ground. The boundary integral method allows to satisfy the no-slip condition by writing a boundary integral equation on the wall of grains of the ground

Styles APA, Harvard, Vancouver, ISO, etc.

36

Ben, Amara Sami. « Ecoulements et transferts thermiques en convection naturelledans les milieux macro-poreux limentaires : application aux réfrigérateurs ménagers ». Paris, Institut national d'agronomie de Paris Grignon, 2005. http://www.theses.fr/2005INAP0006.

Texte intégral

Styles APA, Harvard, Vancouver, ISO, etc.

37

Kostenko, Romaric. « Simulation d'écoulements non-newtoniens en milieu poreux macroscopique par la méthode de lattice-Boltzmann ». Thesis, Université Paris-Saclay (ComUE), 2019. http://www.theses.fr/2019SACLS561.

Texte intégral

Résumé :

Un fluide non-newtonien est un fluide dont la relation entre son taux de déformation et la contrainte qui lui est imposé n'est pas linéaire. Dans un milieu poreux, la contrainte imposée au fluide dépend de la pression imposée mais aussi de la taille des pores, et donc de la perméabilité à l'échelle macroscopique. Certains fluides peuvent avoir une rhéologie qui présente un changement à partir d'un seuil en contrainte. Si la taille des pores est aléatoire, le fluide va alors changer de régime de façon hétérogène dans le milieu. L'écoulement pourra alors présenter un premier régime où tout le fluide est en dessous du seuil, un régime où tout le fluide est au-dessus du seuil, et un régime intermédiaire pour lequel les deux types de rhéologie coexistent. Nous nous intéressons à ce régime intermédiaire pour des écoulements de fluide non-newtonien en milieu poreux macroscopique, étudiés par des simulations. Plus particulièrement, nous regardons l'écoulement d'un fluide de Bingham et celui d'un fluide de Carreau. Le fluide de Bingham ne s'écoule qu'à partir d'une contrainte seuil. En dessous du seuil, il se comporte comme un solide. Au-delà, sa relation taux de déformation/contrainte suis une loi affine. Les fluides de Carreau ont une relation taux de déformation/contrainte qui change de régime entre une loi linéaire et une loi de puissance. L'étude à l'échelle macroscopique se fait en simulant une loi de Darcy-Brinkman dans un champ de perméabilité hétérogène. Nous utilisons pour nos simulations la méthode de Boltzmann sur réseau, avec une grille de nœuds régulière, plus particulièrement le schéma à deux temps de relaxation d'Irina Ginzburg. Pour chaque fluide, nous regardons la relation débit-pression, ainsi que les propriétés géométriques des différents régimes d'écoulement. Nous caractérisons plus particulièrement les propriétés multi-échelles des régions dans un le même régime d'écoulement (clusters), telles que leur taille ou leur forme. Nous faisons aussi le lien entre ces propriétés géométriques et la théorie de la percolation, qui étudie le comportement de cartes de nœuds s'ouvrant aléatoirement et qui prédit des propriétés fractales
A non-newtonian fluid is a fluid which relation between it's shear rate and the stress under which it is put, is not linear. In a porous medium, the stress imposed to the fluid depends on the imposed pressure, but also on the pores size, and therefore on the macroscopic scale permeability. Some fluids have a rheology such that the fluid show a change of behaviour reaching a yield stress. If the pore size is random, then the fluid will present heterogeneous regime changes in the medium. The flow will then show a first regime where the whole fluid will be under the threshold, a regime where the whole fluid will be far above the threshold, and an intermediate regime for which both rheologies coexists. We are interested in intermediate regime for the flow of non-newtonian fluids in macroscopic porous media, and study it with numerical simulations. More particularly, we look at the flow of a Bingham fluid and that of a Carreau fluid. The Bingham fluid doesn't flow under a yield stress. Under the threshold, it behaves as a solid. Beyond, it's shear-rate/stress relation is an affine law. Carreau fluids have a shear-rate/stress relation that change regime between that of a newtonian fluid, and a power law. The macroscopic scale study is done simulating a Darcy-Brinkman law in a heterogeneous permeability field. We use for our simulations the lattice-Boltzmann method, on a regular node grid, and more specifically Irina Ginzburg two relaxation-time scheme. For each fluid, we study the flow-pressure relationship, as well as the geometric properties and the multi-scale properties in the fluid regions in the same flow regime (clusters), properties such as their size and shape. We also link these geometrical properties to the percolation theory, which studies the behaviour of randomly opening node maps and predicts fractal properties

Styles APA, Harvard, Vancouver, ISO, etc.

38

Schmidmayer, Kevin. « Simulation de l'atomisation d'une goutte par un écoulement à grande vitesse ». Thesis, Aix-Marseille, 2017. http://www.theses.fr/2017AIXM0339/document.

Texte intégral

Résumé :

Depuis le début du millénaire, la simulation numérique directe est apparue comme un outil précieux capable d'étudier l’atomisation d’une goutte isolée par un écoulement à grande vitesse. L’atomisation peut être divisée en deux phases distinctes : l'éclatement se produit d'abord sous la forme d'aplatissement de la goutte, formant également des filaments, puis il se poursuit via l'obtention d'une multitude de gouttes de tailles réduites ce qui complète le processus d’atomisation. Les principaux objectifs pour le présent travail étaient donc d’établir un modèle et une méthode numérique capables d’étudier au mieux ces phénomènes. L'atomisation d’une goutte isolée est présentée et est accompagnée d’une comparaison avec l’expérience qui confirme les capacités du modèle et de la méthode à simuler numériquement les différents processus physiques mis en jeu. Des informations essentielles quant aux mécanismes d’atomisation, non exploitables avec l’expérience, sont décrites et l’objectif d’obtenir des gouttes de tailles réduites est atteint
Only at the beginning of the millennium, direct numerical simulation has emerged as a valuable tool capable of studying the atomization of an isolated droplet by a high-speed flow. The atomization can be divided into two distinct phases: the aerobreakup occurs first in the form of flattening of the droplet, also forming filaments, and then it continues via the obtaining of a multitude of reduced sizes droplets what completes the process of atomization. The main objectives of this work were therefore to establish a model and a numerical method able to study these phenomena as well as possible. The atomization of an isolated droplet is presented and is accompanied by a comparison with the experiment which confirms the capacities of the model and the method to numerically simulate the different physical processes involved. Essential information on atomization mechanisms, which cannot be exploited with experiments, is described and the objective of obtaining droplets of reduced sizes is achieved

Styles APA, Harvard, Vancouver, ISO, etc.

39

Nguyen, Khac Tien. « Dégradation thermique de matériaux solides poreux exposés au feu - Simulation numérique avec prise en compte des processus chimiques et mécanismes de transfert ». Thesis, Poitiers, 2014. http://www.theses.fr/2014POIT2344.

Texte intégral

Résumé :

La description de la dégradation des solides est cruciale dans la simulation numérique d'un feu, les gaz de pyrolyse constituant la source combustible amenée au feu. Il importe de décrire avec précision ces processus en fonction des conditions ambiantes. Toutefois, les couplages entre réactions chimiques et mécanismes de transport de masse et de chaleur modifient fortement les comportements. L'objet de ce travail est la prédiction du comportement macroscopique par simulation numérique de la décomposition thermique des matériaux solides. Les propriétés intrinsèques du matériau, parmi lesquelles un modèle chimique incluant un schéma réactionnel et les paramètres thermocinétiques associés, issus d'expérimentation à petite échelle sont considérés ainsi que les mécanismes de transport et les effets limitants qu'ils peuvent introduire. La réponse du matériau est ainsi prédite pour des configurations géométriques et des scénarios d'exposition arbitraires. L’outil spécifique développé au cours de ce doctorat est basé sur une description à l'échelle de Darcy. Il a pour objectif à terme d’être couplé à un code de simulation du feu dans son ensemble, de sorte que les conditions auxquelles chaque élément solide serait exposé seraient décrites tandis que la prédiction de l'évolution du feu dépendrait de la réponse des matériaux fournie par notre module. Un jeu d'applications est présenté pour deux types de bois ayant des schémas réactionnels différents, avec réactions séquentielles ou concurrentielles. Les configurations et scénarios examinés correspondent à des essais normalisés, sous cône calorimètre, et une comparaison entre les résultats numériques et expérimentaux est effectuée
The description of the degradation of solids is crucial in the numerical simulation of a fire, the pyrolysis gases constituting the combustible source. It is important to describe with accuracy these processes, according to the ambient conditions. However, the couplings between chemical reactions and mechanisms of mass and heat transports strongly modify the behavior of the fire.The objective of this work is to predict the macroscopic behavior of the thermal decomposition of solid fuel by numerical simulation. The intrinsic properties of the material, among which a chemical model including a reactive scheme and the associated thermochemical parameters stemming from small-scale experiment are considered as well as the transport mechanisms and the limiting effects that they can introduce. The comportment of the material is then predicted for geometrical configurations and arbitrary scenarios of exposure. The specific tool developed during this Ph-D is based on a description on the Darcy scale. It has in the future to be coupled with a global code of simulation of the fire, so that the conditions undergone by every exposed solid element would be described whereas the prediction of the fire evolution fire would depend on the behavior of the materials supplied by the present module.A set of applications is presented for two wooden materials, featuring different reactive schemes, with sequential or competitive reactions. The configurations and the examined scenarios correspond to normalized tests under cone calorimeter, and a comparison between the numerical and experimental results is made

Styles APA, Harvard, Vancouver, ISO, etc.

40

Lin, Jian. « Modélisation numérique de l’endommagement des matériaux poreux ductiles et application aux géomatériaux ». Thesis, Lille 1, 2009. http://www.theses.fr/2009LIL10143/document.

Texte intégral

Résumé :

Cette étude est consacrée à la modélisation numérique de l'endommagement pour des matériaux poreux ductiles et à des applications aux géomatériaux. On rappelle et analyse d'abord au chapitre 1 différentes approches d'endommagement pour la classe de matériaux étudiés. L'implantation numérique du modèle basé sur un critère macroscopique de type Gurson, récemment proposé par Monchiet et al.(2007) est effectué au chapitre 2. La bonne performance de ce nouveau modèle est vérifiée en comparant ses prédictions numériques à celles d'autres modèles de type Gurson. Dans le chapitre 3, les modélisations numériques sont étendues aux géomatériaux en implantant dans un code éléments finis un nouveau modèle formulé sur la base du critère de Guo et al.(2008) qui prend en compte de la compressibilité de la matrice. L'adéquation du modèle développé aux géomatériaux est démontrée à travers l'étude du comportement d'une craie. L'implantation numérique et l'analyse du modèle basé sur une version anisotrope du critère proposé par Monchiet et al.(2007) sont présentées au chapitre 4. Les résultats numériques issus de cette implantation sont discutés en les comparant aux ceux des autres modèles existants. Ils apparaissent tout à fait satisfaisants
This study is devoted to the numerical modeling of damage for ductile porous materials and its applications to geomaterials. We first recall and analyze in chapter 1 different approaches of ductile damage. We then propose in chapter 2 a model based on the isotropic Gurson-type macroscopic criterion established by Monchietet al.(2007). This model is implemented numerically in ABAQUS Finite Element software. The good performance of this new model is checked by comparing its numerical predictions with the ones provided by other Gurson-type models. In chapter 3, numerical modelings are performed for geomaterials by implementing an appropriate model (based on Guo et al.(2008) criterion) in ABAQUS. The adequacy of the considered model to geomaterials is demonstrated through the study of a chalk. The numerical implementation and analysis of an anisotropic model based on the anisotropic version of the criterion proposed by Monchiet et al.(2007) are presented in chapter 4. The numerical results derived from this implementation are discussed with the results of others existing models. They appear to be quite satisfactory

Styles APA, Harvard, Vancouver, ISO, etc.

41

Chauchat, Julien. « Contribution à la modélisation diphasique du transport sédimentaire en milieux côtiers et estuariens ». Phd thesis, Université de Caen, 2007. http://tel.archives-ouvertes.fr/tel-00173387.

Texte intégral

Résumé :

Cette thèse constitue une contribution à la modélisation numérique diphasique du transport sédimentaire en milieux côtiers et estuariens. L'approche diphasique diffère de l'approche dite classique en résolvant des équations de conservation de la masse, de la quantité de mouvement et de transport pour l'énergie cinétique turbulente séparément pour chaque phase. Les interactions entre le fluide et les particules sont alors représentées par des termes de transfert entre les deux phases. Le travail présenté dans cette thèse porte essentiellement sur la modélisation de la turbulence dans les écoulements diphasiques fluide particule.

En terme de résultats, nous montrons que l'approche développée est capable de reproduire quantitativement les principaux processus mis en jeu dans le transport sédimentaire de particules non cohésives : la sédimentation et la dispersion turbulente des particules en milieu dilué. Le modèle développé confirme l'existence d'une différence de vitesse horizontale entre les particules et l'eau. Il simule les effets de dispersion des particules par le mouvement turbulent du fluide et l'atténuation de l'énergie cinétique turbulente du fluide due à la présence des particules. Une autre originalité de ce travail est de proposer un modèle diphasique à surface libre, bidimensionnel vertical, pour la simulation du transport sédimentaire. Nous avons identifié des lacunes entre les théories et les expériences notamment pour la simulation de la turbulence en écoulement dense. Nous proposons des solutions pour améliorer la simulation du comportement de matériaux cohésifs. Une tentative de simulation hydrosédimentaire sur l'estuaire de la Seine est présentée. Le phénomène de bouchon vaseux est qualitativement reproduit par le modèle sans qu'aucune loi d'érosion ou de dépôt ne soit imposée.

Styles APA, Harvard, Vancouver, ISO, etc.

42

Sciumè, Giuseppe. « Thermo-hygro-chemo-mechanical model of concrete at early ages and its extension to tumor growth numerical analysis ». Phd thesis, École normale supérieure de Cachan - ENS Cachan, 2013. http://tel.archives-ouvertes.fr/tel-00907395.

Texte intégral

Résumé :

The aim of the PhD thesis has been the development of two multi-physics models based on common theoretical basis, but applied to two very different areas: i) the study of the behavior of concrete at early age, essentially for the prevention of early cracking and related issues- ii) the analysis of physical, chemical and biological processes that govern growth and development of cancer. The development of a numerical tool to model concrete at early age is of great importance for the design of durable and sustainable structures. The model has been implemented on the finite element code CAST3M (developed by CEA), also it was validated and nowadays allows multiple applications: study of stresses and cracking phenomena in young concrete, thermal and hygral gradients, autogenous and drying shrinkage, inhibition of hydration caused by drying, creep, stress redistribution, study repairs, etc.. In the fight against cancer, it is clear that the advance of medical strategies based on numerical analysis have a critical scientific interest and can have a great social impact. The equations which govern the thermo-hydro-chemo-mechanical behavior of concrete at early age have may formal analogies with those used to model tumor growth. Hence, these equations have been readapted and a novel mathematical model for tumor growth has been developed. The model was implemented in Cast3M and the first numerical results have been encouraging since very close to the experimental data present in the literature.

Styles APA, Harvard, Vancouver, ISO, etc.

43

Lavie, Guillaume. « Simulation et analyse numérique de procédés de récupération de pétrole caractérisés par des fronts raides ». Thesis, Pau, 2011. http://www.theses.fr/2011PAUU3028/document.

Texte intégral

Résumé :

Du fait de la raréfaction des gisements de pétrole conventionnel et de l'accroissement de la demande mondiale, les compagnies pétrolières sont amenées à considérer de nouvelles réserves inexploitées car trop coûteuses à l'exploitation il y a quelques années. Ce travail est consacré à la simulation numérique d'un procédé d'extraction de pétrole lourd, le VAPEX, caractérisé par l'injection de solvant facilitant le drainage par gravité de l'huile. L'enjeu principal de cette thèse est de simuler numériquement et avec précision, le mécanisme de pénétration du solvant dans l'huile, cette zone de pénétration étant relativement mince à l'échelle du réservoir. Notre travail a porté sur l'utilisation d’algorithmes de raffinement adaptatif de maillages pour la simulation de ce procédé. Ces considérations nous ont amené à considérer la problématique d'application d'algorithmes de raffinement de maillages en milieu poreux, et plus particulièrement en milieux poreux hétérogènes pour lesquels les indicateurs (ou estimateurs) existants permettant de déterminer les zones à raffiner ne sont pas clairement établis. Dans un premier temps, nous mettons en place l'ensemble des équations que nous sommes amené à résoudre numériquement au cours de l'élaboration des codes de simulation de procédés d'extraction par injection d'eau et de solvant. Nous décrivons ensuite le procédé VAPEX, et nous établissons l'étude semi-analytique de Butler & Mokrys du procédé à laquelle nous ajoutons la prise en compte des termes capillaires. Nous décrivons ensuite deux méthodes mathématiques permettant d'établir des estimateurs a posteriori pour de tels problèmes. Au cours de ce travail, nous mettons en place deux codes de simulation de réservoirs, le premier est un simulateur basé sur une discrétisation des équations par une méthode de discrétisation par éléments finis mixtes hybrides et permettant de simuler des problèmes d'injection d'eau et le second est basé sur une méthode de discrétisation par volumes finis permettant de simuler des procédés d'injection de solvant tels que le VAPEX. Nous utilisons ensuite ces estimateurs par le biais d’un algorithme de raffinement de maillages pour ces deux codes de simulation. L'étude du modèle semi-analytique que nous avons faite du procédé VAPEX montre que sous certaines hypothèses le terme capillaire peut être assimilé à un terme diffusif pour prédire l'avancement du front. Les estimateurs mis en place basés sur les variations des flux permettent de déterminer les zones de front pour des réservoirs hétérogènes. Nous illustrons ce résultat sur une simulation d'un problème d'injection d'eau à l'aide du simulateur mis en place. Les estimateurs mis en place permettent également de capter les fronts des saturations des phases et de concentration des constituants. Nous illustrons ce résultat sur une simulation du procédé VAPEX.Ce travail a amené à la mise en place d’estimateurs permettant de raffiner des maillages autour des fronts de saturation et de concentration pour des milieux poreux hétérogènes au cours de la simulation d'écoulements en milieux poreux. Pour caractériser ces estimateurs, nous avons distingué un front de saturation (resp. concentration), qui correspond est une forte variation de saturation (resp. concentration) engendrant une forte variation du flux, d'un fort gradient de saturation (resp. concentration) qui apparait automatiquement lorsque la perméabilité du milieu poreux varie brutalement. Les extensions de ce travail pourraient permettre de mettre en œuvre plus de simulations faisant varier les termes capillaires et les termes dispersifs affin de mieux apprécier le comportement des estimateurs introduits. Il serait également intéressant d'étudier l'impact des hypothèses faites pour établir les estimateurs empiriques utilisés. Enfin, de tels estimateurs pourraient s'appliquer à des codes de simulation tenant compte de la dispersion en milieu poreux
With rarefaction of conventional petroleum fields and the world supply increase, petroleum companies have to consider new unexploited reserves because of their exploitation cost few years ago. These reserves are called unconventional reserves, mainly heavy oil or extra heavy oil. This work is devoted to the numerical simulation of an extraction process of heavy oil, the VAPEX. VAPEX is characterized by solvent injection in a horizontal well in order to facilitate gravity drainage of oil in a producer well situated in the same way under the injector well. The main stake to simulate this process is to simulate the mechanism of solvent penetration in oil with precision, this penetration zone being very thin drawn to the reservoir scale. We focus our work on application of algorithm of adaptive mesh refinement to simulate this process. This consideration lead us to consider the issue of application of mesh refinement in porous media, especially in heterogeneous porous media for which existing indicators (or estimators) do not allow to determine properly the zone to be refined. The objective of this work is to determine a posteriori estimators for implementation of simulation codes of heterogeneous petroleum reservoirs and their application to the VAPEX process. Firstly, we set the whole equations we have to solve numerically to build simulators of petroleum extraction processes by water injection and solvent injection. Then, we describe the VAPEX process and we establish the semi-analytical study of Butler & Mokrys of this process for which we take in account influence of capillarity. Afterwards, we describe two mathematical methods to establish a posteriori estimators for such problems. One of these methods is based on empirical extrapolation of existing estimators for hyperbolic problems given by Ohlberger. In the course of this work, we set two simulation codes of petroleum reservoirs, the first is a simulator based on equations discretisation by mixed finite element method allowing simulation of water injection problems and the second is based on discretisation by finite volumes method allowing simulation of solvent injection processes like the VAPEX. Then, we apply these estimators to an algorithm of mesh refinement for these two simulation codes. The study of the semi-analytical model set of the VAPEX process show that under hypothesis, capillarity can be likened to a diffusive term to determine the front tracking. Estimators set are based on flux variation and allow determining sharp zone to be refined in heterogeneous porous media. We illustrate this result for a simulation of water injection with the simulator built. The estimators set can also track saturations fronts and concentrations fronts. We illustrate this result on a simulation of the VAPEX process. This work shows it is possible to set estimators allowing mesh refinement to track saturations and concentrations fronts during simulation of flow in heterogeneous porous media. To set such estimators, we make a distinction between saturation (resp. concentration) front and high saturation (resp. concentration) gradient. To our point of view, a front is a high variation of saturation or concentration that has a high impact on the flux variation. It is different from a high gradient that appear automatically where the permeability of the porous media has a high variation. This work could be complemented by making more simulations with variations of capillarity and dispersive terms in order to have a better appreciation of estimators introduced. It would be advisable to study the impact of suppositions done to establish these empirical estimators set. After all, such estimators could be applied to a simulation code taking in account dispersion in porous media

Styles APA, Harvard, Vancouver, ISO, etc.

44

Wanko, Ngnien Adrien. « Etude des mécanismes de transfert et évaluation des capacités d'oxygénation et de traitement des dispositifs de traitement par cultures fixées sur supports granulaires fins alimentés en discontinu-développement d'un modèle de biodégradation ». Université Louis Pasteur (Strasbourg) (1971-2008), 2005. http://www.theses.fr/2005STR13187.

Texte intégral

Résumé :

L'étude des mécanismes de transfert et évaluation des capacités d'oxygénation et de traitement des dispositifs de traitement par cultures fixées sur matériaux granulaires fins alimentés en discontinu. – Développement d'un modèle de biodégradation" est l'objet du travail mené visant une meilleure compréhension des mécanismes majeurs s'établissant au sein des lits d'infiltration percolation et permettant au moyen de certains choix opérationnels, une pérennisation de la filière. Tout d'abord, dans une première partie, nous exposons l'état d'avancement de la recherche dans ce domaine. Après des généralités sur les filières de traitement des eaux usées, il est fait cas de l'infiltration percolation sur sable; de ses spécificités, des contraintes qu'elle impose et des attentes au niveau de la recherche. La problématique du choix de type de sable qui constitue le matériau de garnissage ou support de la biomasse en relation avec les capacités d'oxygénation du milieu poreux et les risques de colmatage demeure un souci permanent. L'intérêt de la seconde partie essentiellement expérimentale porte sur la caractérisation physique et hydrodynamique des sables roulés et concassés constituant les massifs filtrants en milieux non colonisés. L'hydrodynamique des massifs filtrants a permis une analyse comparative des distributions de temps de séjour du fluide percolateur en milieu vierge et en régime continu puis discontinu. Indépendamment des supports étudiés, les temps de résidence du traceur sont d'autant plus élevés que l'alimentation se fait à fréquences espacées et à volume par bâchée considérable. En revanche, pour chacune des fréquences d'alimentation testées, la récupération du traceur s'est effectuée plus rapidement au sein de la colonne constitué des sables alluvionnaires. Le second volet de l'étude hydrodynamique est effectué en milieu colonisé. L'efficacité du traitement biologique du carbone est supérieure pour un sable d'origine alluvionnaire comparativement à un sable concassé ;la tendance s'inversant significativement lorsqu'il s'agit de l'abattement de l'azote. Cet efficacité est directement corrélée à l'oxygénation du massif filtrant autrement dit au renouvellement de la phase gazeuse du milieu poreux. Exploitant une technologie de fluorescence nous permettant de détecter aussi bien l'oxygène gazeux que dissous, nous parvenons à déterminer les différents ordres de grandeur des phénomènes de convection et de diffusion ainsi que leur importance relative dans la ré-oxygénation d'un massif filtrant. Le développement d'un modèle numérique simulant l'infiltration et le traitement des eaux usées à travers un milieu poreux colonisé constitue l'épilogue de ce travail. Deux approches de résolution sont utilisées : une approche globale qui résout simultanément à l'aide d'un unique schéma les parties transport et chimie de chaque équation et une approche par séparation d'opérateurs qui résout séparément et successivement les parties convection, dispersion et cinétique
Aim of this work is to enhance understanding of the main phenomena occurring during infiltration of wastewater within vertical flow sand filter. We also introduce some operating possibilities which may improve efficiency and keep seepage beds in good working order. In the first part we document state of art in this research topic. Following some general remarks about different wastewater treatment techniques, more precise description of seepage beds is made with particular interest on its constraints and needs of further research. A central issue is the choice of porous media. It should be adapted to allow good oxygen feeding for bacteria growing on it and avoid biological clogging. Second part, mainly experimental, is focused on physical and hydraulic characterization of un-colonized crushed and stream sand. We carried out a residence time distribution (RTD) comparison of these two sands under steady and unsteady flow. Apart from characteristics of porous media and with a constant daily loading rate, decreased frequencies involve longer residence times. On the other hand, output recovery of the mass of injected tracer always occurs faster with stream sand. Continuation of hydrodynamic study is made with colonised porous media. Efficiency of biological removing of carbonaceous pollution is higher within stream sand. Opposite trend is observed with nitrogenous pollution. Removing of both nitrogenous and carbonaceous pollutions is linked to oxygenation. In other words it is linked to renewal of gaseous phase in porous medium. We highlight convection and diffusion of oxygen in porous media by using a fluorescence technology, which allows us to detect both gaseous and dissolved oxygen. We compare participation of convection and diffusion phenomena to renewal of oxygen. Achievement of this work consists on the development of a numerical code which simulates flow and pollution degradation through a colonized porous medium. Two different numerical methods have been used. In the first model, also called global approach, chemical and transport parts of each equation are solved by a single numerical scheme. A second model introduces a splitting operator method which solves separately and successively convective, dispersive and kinetic parts of each equation

Styles APA, Harvard, Vancouver, ISO, etc.

45

Dartois, Arthur. « Étude de la macro-dispersion de particules inertes dans des milieux poreux 3D fortement hétérogènes ». Thesis, Poitiers, 2016. http://www.theses.fr/2016POIT2317/document.

Texte intégral

Résumé :

Les milieux poreux font partie des formations géologiques assez répandue dans la nature et son sujet d'études intensives. L'engouement de ce sujet vient des multiples secteurs d'applications de ces recherches et leur importance dans notre société. Que ce soit de la part des sociétés pétrolières qui souhaitent optimiser leurs moyens de productions, les agences de contrôles environnementaux qui veulent prévenir la contamination de nappe phréatique et la fuite de déchets nucléaires ou encore des industriels avec des problèmes de drainages et de réhabilitation de mines, tous ces acteurs dépendent des recherches faites dans ce domaine. Cependant, un des principaux problèmes de ce sujet est l'inaccessibilité des milieux que nous voulons étudier. Pour palier à cela de nombreuses équipes se sont tournées vers la simulation numérique. Cette thèse s'inscrit dans ce cadre et utilise le module PARADIS du logiciel d'hydrogéologie H2olab pour modéliser le transport de particules dans des milieux poreux fortement hétérogènes. Grâce aux données obtenues et à des comparaisons avec la littérature nous montrerons l'effet du passage au 3D sur la topologie de l'écoulement et les répercussions sur le transport de particules ainsi que l'effet de la diffusion moléculaire sur les coefficients de macro-dispersion. Enfin nous proposerons deux lois de transport reliant macrodispersion, variance du champ de perméabilité et diffusion moléculaire
Heterogeneous porous media have been intensively studied these last fifty years. The popularity of this subject come from the multiple areas where these researches can be applied and their importance to our society. Whether from the oil companies that want to optimize their methods of production, environmental control agencies who want to prevent contamination of ground water and leakage of nuclear waste or industrial with drainage issues and mine rehabilitation, all these actors depend on research done in this area. However, one of the main problems of this subject is the accessibility of these porous media which are often several hundred meters underground. To overcome this, many teams have turned to computer simulation. This thesis is among them and uses the PARADIS module from the hydrogeology software H2olab to model particle transport in highly heterogeneous porous media. Thanks to the data obtained and comparisons with the literature, we show the impact of switching from a 2D to a 3D porous media on the ow topology and the repercussions on the particle transport. Furthermore, we also investigated the effect of molecular diffusion coefficients on macro-dispersion. Finally, we will propose two empirical functions linking macro-dispersion variance of the permeability field and molecular diffusion

Styles APA, Harvard, Vancouver, ISO, etc.

46

Bary, Benoît. « Étude du couplage hydraulique-mécanique dans le béton endommagé ». Cachan, Ecole normale supérieure, 1996. http://www.theses.fr/1996DENS0022.

Texte intégral

Résumé :

Les travaux effectues dans le cadre de la thèse portent sur l'étude du couplage hydromécanique dans le béton, avec pour objectif une application au calcul prédictif du comportement des barrages. L'approche proposée est basée sur la mécanique des milieux poreux saturés couplée à la mécanique de l'endommagement anisotrope. Sont ainsi pris en compte les transferts de fluide sous pression, les phénomènes de dégradation orientée de la matrice, et les répercussions de cet endommagement sur les paramètres de couplage hydromécaniques. Une procédure expérimentale autorisant le pilotage simultané de la pression interstitielle et des contraintes va permettre l'identification des paramètres introduits dans le modèle. Les évolutions de ces paramètres sont alors simulées et comparées à d'autres approches proposées pour ce type de problème. L'implantation numérique du modèle ainsi identifié est réalisée dans le code de calcul par éléments finis Castem 2000. Une méthode de régularisation du problème de localisation des déformations est mise en œuvre avec succès pour le modèle d'endommagement. Une simulation d'un essai de fissuration hydraulique sur une structure pré-entaillée est alors effectuée, et prouve que le caractère dégradable des paramètres hydromécaniques influence fortement la réponse post-pic. On procède ensuite a la simulation du comportement d'un barrage poids soumis a une crue. Les résultats produits permettent de retrouver les principaux aspects mis en évidence sur ouvrages réels

Styles APA, Harvard, Vancouver, ISO, etc.

47

Wise, John Nathaniel. « Inverse modelling and optimisation in numerical groundwater flow models using proper orthogonal decomposition ». Thesis, Saint-Etienne, EMSE, 2015. http://www.theses.fr/2015EMSE0773/document.

Texte intégral

Résumé :

Des simulateurs numériques sont couramment utilisés pour la prédiction et l'optimisation de l'exploitation d'aquifères et pour la détermination de paramètres physiques (e.g perméabilité) par calcul inverse. L'équation de Richards, décrit l'écoulement d'un fluide dans un milieu poreux non saturé. C'est une équation aux dérivées partielles non linéaires, dont la résolution numérique en grande dimension 3D est très coûteuse et en particulier pour du calcul inverse.Dans ce travail, une méthode de réduction de modèle (ROM) est proposée par une décomposition orthogonale propre (POD) afin de réduire significativement le temps de calcul, tout maîtrisant la précision. Une stratégie de cette méthode est de remplacer localement dans l'algorithme d'optimisation, le modèle fin par un modèle réduit type POD. La méthode de Petroc-Galerkin POD est d'abord appliquée à l'équation de Richards et testée sur différents cas, puis adaptée en linéarisant les termes non linéaires. Cette adaptation ne fait pas appel à une technique d'interpolation et réduit le temps de calcul d'un facteur [10;100]. Bien qu'elle ajoute de la complexité du ROM, cette méthode évite d'avoir à ajuster les paramètres du noyau, comme c'est le cas dans les méthodes du POD par interpolation. Une exploration des propriétés d'interpolation et d'extrapolation inhérentes aux méthodes intrusives est ensuite faite. Des qualités d'extrapolation intéressantes permettent de développer une méthode d'optimisation nécessitant de petits plans d'expériences (DOE). La méthode d'optimisation recrée localement des modèles précis sur l'espace des paramètres, en utilisant une classification à vecteurs de support non linéaire pour délimiter la zone où le modèle est suffisamment précis, la région de confiance. Les méthodes sont appliquées sur un cas d'école en milieu non saturé régit par l'équation de Richards, ainsi que sur un aquifère situé dans le "Table Mountain Group" près de la ville du Cap en Afrique du Sud
The Richards equation describes the movement of an unsaturated fluid through a porous media, and is characterised as a non-linear partial differential equation. The equation is subject to a number of parameters and is typically computationnaly expensive to solve. To determine the parameters in the Richards equation, inverse modelling studies often need to be undertaken. As a solution to overcome the computational expense incurred in inverse modelling, the use of Proper Orthogonal Decomposition (POD) as a Reduced Order Modelling (ROM) method is proposed in this thesis to speed-up individual simulations. The Petrov-Galerkin POD approach is initially applied to the Richards equation and tested on different case studies. However, due to the non-linear nature of the Richards equation the method does not result in significant speed up times. Subsquently, the Petrov-Galerkin method is adapted by linearising the nonlinear terms in the equation, resulting in speed-up times in the range of [10,100]., The adaptation, notably, does not use any interpolation techniques, favouring an intrusive, but physics-based, approach. While the use of intrusive POD approaches add to the complexity of the ROM, it avoids the problem of finding kernel parameters typically present in interpolative POD approaches. Furthermore, the interpolative and possible extrapolation properties inherent to intrusive PODROM's are explored. The good extrapolation propertie, within predetermined bounds, of intrusive POD's allows for the development of an optimisation approach requiring a very small Design of Experiments (DOE). The optimisation method creates locally accurate models within the parameters space usign Support Vector Classification. The limits of the locally accurate model are called the confidence region. The methods are demonstrated on a hypothetical unsaturated case study requiring the Richards equation, and on true case study in the Table Mountain Group near Cape Town, South Africa

Styles APA, Harvard, Vancouver, ISO, etc.

48

Hu, Dawei. « Contribution à l’étude du comportement mécanique et hydromécanique d’une roche endommagée ». Thesis, Lille 1, 2009. http://www.theses.fr/2009LIL10150/document.

Texte intégral

Résumé :

Cette thèse est consacrée aux études expérimentales et numériques du comportement mécanique et poromécanique et de l’évolution de la perméabilité du grès en prenant en compte l’anisotropie du matériau. Après une analyse par la microscopie électronique à balayage (MEB) et les rayons X, les échantillons sont utilisés pour réaliser des essais de compression hydrostatique et de compression triaxiale sous différentes pressions de confinement. Les évolutions des coefficients de Biot et de la perméabilité sont mesurées respectivement et étudiées en prenant en compte le comportement mécanique. Un modèle élastoplastique couplé à l’endommagement en se basant sur l'approche discrète est proposé et appliqué aux grés secs. Les essais en compression sont simulés, l’évolution et la distribution de la variable d’endommagement et de l’écrouissage plastique sont analysées. Le modèle proposé est validé de nouveau par les simulations des essais réalisés (y compris le comportement mécanique, les coefficients de Biot et la perméabilité) sur des échantillons saturés. Une très bonne concordance est obtenue entre les prédictions et les données expérimentales. Enfin, on développe un modèle anisotropique poroplastique couplé avec l’endommagement afin de décrire les propriétés poromécaniques du grès saturé. Le modèle proposé est appliqué pour décrire le comportement du matériau dans les essais triaxiaux en condition non-drainée, et dans les essais triaxiaux en appliquant une pression interstitielle. En général, la aussi nous obtenons une bonne concordance entre modélisation et expérimentation
This thesis presents the laboratory and numerical investigations on anisotropic mechanical, poromechanical and fluid transport behavior of sandstone. After the X-ray and SEM tests, the samples were stressed under triaxial loading to study the initial and stress-induced mechanical behavior. The evolution of anisotropic Biot’s coefficients and axial permeability are also investigated and related to mechanical behavior. A plastic damage model is presented based on discrete approach. The proposed model is applied to simulate the mechanical behavior of dry sandstone under different loading conditions, the evolution and distribution of damage variable and plastic hardening variable are also discussed. Furthermore, the proposed model is applied to simulate the previous experimental results (including mechanical behavior, Biot's coefficients and permeability) of saturated sandstones. The comparisons between the numerical result and test data show a good performance of the present model. Finally, an anisotropic poroplastic damage model is developed to describe the poromechanical behavior of saturated sandstone. The comparisons between numerical simulation and experimental data is given for the triaxial compression tests under drained and undrained condition, the tests of strain response to pore pressure increment is also simulated, a good accordance is obtained

Styles APA, Harvard, Vancouver, ISO, etc.

49

Mnasri, Faiza. « Étude du transfert de chaleur et de masse dans les milieux complexes : application aux milieux fibreux et à l’isolation des bâtiments ». Thesis, Université de Lorraine, 2016. http://www.theses.fr/2016LORR0169/document.

Texte intégral

Résumé :

Le contexte énergétique international impose de nouvelles orientations au secteur du bâtiment neuf ou en rénovation. Toute nouvelle solution doit être techniquement efficace et respectueuse pour l’environnement. Il s'agit dans ce travail de thèse de réaliser une étude numérique et expérimentale de matériaux de construction biosourcés liés au contexte transfrontalier Lorrain (France-Belgique- Luxembourg). En effet, ce travail intègre une partie du projet européen « Ecotransfaire » mené pour le développement d'une filière durable propre aux éco-matériaux. La sélection des matériaux selon une liste de critères à la fois scientifiques, géographiques et environnementaux a permis de répondre à notre problématique en s'orientant vers l'intégration des matériaux biosourcés pour leurs aspects favorables à l'environnement et à l’efficacité énergétique du bâtiment. Intégrés au bâtiment, ces matériaux sont sujets à plusieurs phénomènes de transfert de chaleur et de masse. Dans un premier temps et pour mieux appréhender ces phénomènes, un modèle de transfert couplé de chaleur, d'air, d'humidité (HAM transfers) est utilisé pour simuler le comportement hygrothermique d’un matériau en bois massif à structure supposée homogène. Ce modèle, mis en œuvre et résolu par la méthode des éléments finis, a été validé par des résultats analytiques retenus dans la littérature. L'étude de sensibilité du modèle au couplage, aux dimensions dans l'espace, aux conditions aux limites et aux variabilités des paramètres d'entrée est également présentée. Une des difficultés de l’utilisation de ce modèle réside dans la prise en considération de l'aspect fortement hétérogène de certains matériaux. Ainsi, dans ce travail, nous proposons une approche de caractérisation d'un composite lignocellulosique hétérogène de structure poreuse. En effet ce matériau est composé de deux constituants bien connus dans le domaine de l’industrie de construction: Le bois et le ciment. Le bois est incorporé sous forme de granulats avec des formes et des tailles irrégulières et le ciment est utilisé comme un liant. Le travail réalisé permet de remonter aux propriétés intrinsèques équivalentes de ce matériau (conductivité thermique et perméabilité à la vapeur) à l’aide des techniques de micro-tomographie. La méthodologie suivie consiste à la détermination de la structure d'échantillon par une prise d'images à l'échelle microscopique. Une fois la structure de l’échantillon générée, une reconstruction de la représentation bidimensionnelle précède la génération de la structure tridimensionnelle à l'aide d’un outil numérique qui permet de déterminer les propriétés équivalentes des domaines reconstruits en 3D. La perméabilité et la conductivité thermique équivalentes sont les deux propriétés évaluées dans cette configuration. Ces deux propriétés dépendent fortement de la porosité et de la distribution des pores dans la phase continue (la phase solide). De plus la composition de ce matériau et les fractions volumiques de chacun de ses constituants influent sur la formation de sa microstructure et par conséquent sur ses propriétés de transferts thermiques et hydriques. L'ensemble des connaissances développées dans ce travail permet une piste sérieuse pour l'élaboration d'un éco-matériau à propriétés contrôlées pour des usages spécifiques dans la construction et la rénovation
International energy context requires a new orientation to the building sector as in construction or in renovation. Any new solution must be technically efficient and environmentally acceptable. In this thesis, the object is to achieve a numerical and experimental analysis of a building biobased materials. Some of these materials are included from the study of a transborder project to the Lorraine region (France, Belgium and Luxembourg). Indeed an Ecotransfaire project was included in this work. This project has been oriented to the development of a sustainable eco materials chain. A process of analysis has been established in order to select the materials candidates on the basis of scientific, geographical and environmental criteria. The answers are moving towards the integration of bio-based materials. These materials are subject of several heat and mass transfers phenomena. So understanding these mechanisms within a building material has been achieved firstly. This resulted on a coupled model of heat transfer, air, moisture experienced by the HAM model. This model is applied to a wooden building material whose its structure is assumed homogeneous. Then, this model was implemented and solved by the finite element method. Its numerical solution is validated by analytical results available in the literature. The study of sensitivity of the model coupling, dimensions in space, the boundary conditions and the variability of input parameters is also presented. One of the difficulties of using this model is the case of heterogeneous materials. Thus, in this work, we propose an approach of characterization of a heterogeneous lignocellulosic composite material with a porous structure. In fact, this material is composed of two components: Wood and cement. The wood is presented by a shapes aggregates with irregulars sizes and the cement is considered as the binder in the composition. The object was to predict its equivalent intrinsic properties (thermal conductivity and vapor permeability) by using the micro-tomography techniques.The methodology consists to determine the structure of the sample by taking images at the microscopic scale. Once the structure of the sample is generated, we will conduct from a reconstruction of the two-dimensional representation to a three dimensional structure by using a numerical tool which determines the equivalent properties of the 3D reconstructed domain. The permeability as well as the equivalent thermal conductivity are the two properties evaluated in this configuration. These two properties are strongly depend to the porosity and to pore distribution in the continuous phase (the solid one). Moreover the composition of the material and the volume fractions of each components influence the formation of microstructure and consequently the thermal and hydric transfers

Styles APA, Harvard, Vancouver, ISO, etc.

50

Fanzar, Abdelaziz. « Développement d’un code numérique pour la simulation et l’étude de l’hydrodynamique et de la physico-chimie de milieux diphasiques incompressibles. Cas d’une goutte d’eau dans l’huile de paraffine ». Electronic Thesis or Diss., Aix-Marseille, 2014. http://www.theses.fr/2014AIXM4339.

Texte intégral

Résumé :

Depuis plusieurs décennies, une importante activité scientifique se concentre sur la description numérique, théorique ou expérimentale de l'hydrodynamique des écoulements multiphasiques. Ces écoulements sont caractérisés par l'existence d'interfaces, et d'une force à l'interface, la tension superficielle, séparant généralement deux fluides non miscibles. Un cas d'étude dans ce contexte est le problème du drainage d'une unique goutte dans une phase continue, l'ensemble étant soumis à la gravité. Ce système fait apparaître des écoulements récemment décrits pour une goutte d'eau dans l'huile de paraffine. Ce système constitue également un modèle simple pour l'étude des propriétés aux interfaces, Mais d'un point de vue numérique, se pose alors le problème de la stabilité des algorithmes pouvant être utilisés. Les effets aux interfaces impliquent en effet des domaines spatiaux très limités dans lesquels les grandeurs physiques entre les deux fluides sont discontinues. D'importants artéfacts numériques peuvent alors être générés dans les simulations et faire perdre la richesse de la physico-chimie du système considéré. Le problème de la simulation d'écoulements multiphasiques intéresse aussi bien le monde académique que le monde industriel. L'objectif de ce travail de thèse est donc d'implémenter les techniques numériques les plus récentes et de développer un code pour permettre la simulation de l'hydrodynamique de systèmes dispersés. Pour parvenir à ce but, il reste encore des problèmes algorithmiques importants à résoudre comme la prise en compte des effets thermocapillaires et thermosolutaux. Ces deux derniers points sont l'objet de cette thèse
For several decades, an important scientific activity has focused on the numerical, theoretical and experimental hydrodynamics of drops. This work presents numerical results of a single droplet in the gravity field and in non-isothermal conditions. The simulation such a multiphase system is important in both academic and industrial world. This is particularly the case in the field of emulsions, wetting problems and evaporation. To achieve this goal, there are still important algorithmic problems due to the free moving interfaces and the description of capillary effects. Here, a Volume of Fluid technique has been implemented with high order temporal and spatial schemes to preserve the sharpness of the drop interface. The system under consideration is a simplified model consisting in a single water droplet in a continuous paraffin oil phase. These liquids are immiscible and non-compressible and the overall evolution is unsteady. Capillary contributions such as temperature and surfactant dependent surface tension are fully accounted for. This presentation is aimed to show the capabilities of VOF techniques for the simulations of unsteady multiphase systems in non-isothermal configurations. The role of the droplet initial position and temperature field is described with good numerical stability. There are still important problems remaining in the simulation of free interface systems with such a technique. Spurious currents induced by the description of capillarity can in particular come into play. But these latter can be controlled once the droplet average velocity due to drainage becomes large enough

Styles APA, Harvard, Vancouver, ISO, etc.

Thèses sur le sujet « Écoulement en milieux poreux – Simulation par ordinateur »

Créez une référence correcte selon les styles APA, MLA, Chicago, Harvard et plusieurs autres