Academic literature on the topic 'Écoulement en milieux poreux – Simulation par ordinateur'

: 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<br>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

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. <br /> 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.<br /> 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. <br /> 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.

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<br>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

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.