Gotowa bibliografia na temat „Écoulements incompressibles diphasiques”

Le modèle de Chavent est utilise pour étudier l'écoulement diphasique de fluides incompressibles, dans un milieu hétérogène, au cours de la récupération secondaire de pétrole. Les équations du milieu homogène sont déduites des équations hétérogènes. Une simulation du déplacement de l'huile par de l'eau dans un gisement hétérogène forme de cellules informement espacées est réalisée et comparée a celle du déplacement de l'huile dans un réservoir homogène. Les simulations numériques sont réalisées avec le code Bidimix (inria et chevron code (Wyoming). Un logiciel de calcul permet d'obtenir les caractéristiques du milieu poreux homogène équivalent. L'étude est étendue aux écoulements diphasiques incompressibles dans un réservoir a deux types de roches

Deux approches au problème de changement d'échelle sont confrontées: les méthodes de calcul de pseudo-fonctions, la prise de moyenne volumique a grande échelle. Cette dernière fondée sur une théorie mathématique plus riche donne des indications précieuses sur les pseudo-fonctions. L’analyse pratique de ces méthodes se fait a l'aide de simulations numériques d'écoulements diphasiques en milieux poreux homogènes et stratifies sous divers régimes d'écoulement. Dans certains cas pratiques, la méthode de prise de moyenne à grande échelle fournit des résultats aussi satisfaisants que ceux obtenus par des méthodes beaucoup plus couteuses et élaborées. Enfin, des règles générales pour le calcul des propriétés effectives en milieux poreux aléatoires sont établies.

On étudie ici, dans un domaine borné, un problème de perturbations singulières pour une classe d'équations paraboliques non linéaires du second ordre lorsque les effets de la diffusion sont négligeables devant ceux de la convection. Ce travail est effectué pour deux types de conditions de bord (Dirichlet, mêlées). On établit, dans un espace approprié, la convergence de la suite de solutions de ces problèmes vers la solution anthropique du problème limite hyperbolique non linéaire du premier ordre. La difficulté essentielle est, compte tenu du caractère non linéaire de l'équation, de trouver des estimations suffisantes pour cette suite de solutions afin d'établir un résultat de convergence. Pour cela, on utilise entre autre des techniques de régularisation, la notion de fonctions à variation bornée et une classe de fonctions, introduite par F. Mignot et J. -P. Puel pour étudier des problèmes hyperboliques linéaires du premier ordre, ayant pour principale propriété celle de s'annuler dans un voisinage de la surface définissant la couche limite. Ce travail est divisé en trois chapitres. Dans les deux premiers, on examine les problèmes non dégénérés qui fournissent une approximation des problèmes dégénérés considérés dans le troisième.

Alors que la simulation aux grandes échelles (L.E.S.) des écoulements monophasiques est largement répandue même dans le monde industriel, ce n'est pas le cas pour la L.E.S. d'écoulements diphasiques avec interface (c'est-à-dire d'écoulements où les deux phases liquide et gazeuse sont séparées par une interface). La difficulté majeure réside dans le développement de modèles de sous-maille adaptés au caractère diphasique de l'écoulement. Le but de ce travail est de générer une base D.N.S. dans le cadre d'écoulements diphasiques turbulents avec interface pour comprendre les interactions entre les petites échelles turbulentes et l'interface. Les différents termes sous-maille proviendront d'une analyse a priori de cette base D.N.S. Pour mener à bien ce travail, différentes techniques numériques sont testées et comparées dans le cadre de configurations turbulentes où de grandes déformations interfaciales apparaissent. Puis, l'interaction interface/turbulence est étudiée dans le cadre où les deux phases, séparées par une interface largement déformée, sont résolues par une approche D.N.S. La configuration retenue est l'interaction entre une nappe initialement plane et une T.H.I. libre. Les rapports de densités et de viscosités sont fixés à 1 pour se concentrer sur l'effet du coefficient de tension de surface. Une étude paramétrique sur le nombre de Weber est menée. Finalement, un filtrage a priori de la base D.N.S. est réalisé et les termes sous-maille qui en découlent sont comparés les uns aux autres.

Ce travail est consacré à l'étude de l'unicité des solutions de systèmes de lois de conservation de masse (équations de continuité) issus de la modélisation des écoulements diphasiques de fluides en milieu poreux. Autant le problème de l'existence semble complètement résolu (chap1), autant celui de l'unicité globale en dimension quelconque posait encore des difficultés principalement a cause du manque de résultats de régularité sur les solutions et du couplage du système. Nous obtenons un premier résultat d'unicité pour un problème pose dans un domaine géométrique particulier (l'espace compris entre deux cylindres de même hauteur et de même axe principal), ce cas pouvant se ramener, dans certaines conditions d'exploitation, au cas des écoulements monodirectionnels pour lesquels le système se découple naturellement par une quadrature. L’unicité s'obtient alors par une technique à la Kruskov, usuelle pour les problèmes hyperboliques du premier ordre et utilisant des approximations lipchitziennes de la fonction de Heaviside (chap2). Une étude générale est ensuite menée sur la régularité des solutions d'équations de type elliptique. Outre une généralisation du théorème de Meyers nous obtenons des résultats permettant particulièrement de pouvoir maitriser et contrôler la régularité du gradient de pression (chap3). Les résultats ainsi obtenus permettent alors de s'affranchir de l'hypothèse que le gradient de pression est borne (hypothèse généralement admise dans la littérature mais pas toujours acquise). On obtient alors des résultats d'unicité pour des systèmes simplifies en utilisant une argumentation usuelle pour des équations de Navier-stokes (chap4) puis d'autres propriétés d'unicité sont obtenus dans le cadre de méthodes de transposition inspirées par S. N. Antontsev (chap5)

Dans ce travail de thèse, on s'est intéressé à la simulation d'écoulements incompressibles multi-échelles et multiphasiques. L'une des principales difficultés numériques est l'introduction d'une diffusion numérique due aux schémas utilisés. Celle-ci étant indépendante du maillage, une possibilité est de simuler ce type d'écoulement avec un très grand nombre points. Cependant, les besoins en ressources informatiques et en temps deviennent rapidement importants. On a donc développé une méthode de raffinement de maillage adaptatif (AMR) dans le but de suivre, soit des interfaces dans un écoulement diphasique, soit des fronts de concentration dans un écoulement monophasique avec transport d'une espèce inerte, de manière précise tout en optimisant le temps CPU et
la taille mémoire. On montre au travers de cas d'étude 2D et 3D, judicieusement choisis, l'efficacité de cette méthode.

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

Dans cette thèse, nous nous intéressons à la simulation numérique d'écoulements instationnaires de deux fluides visqueux non miscibles, séparés par une interface mobile. Plus particulièrement des écoulements sans choc constitués d'une phase gazeuse et d'une phase liquide sont considérés. Pour modéliser de tels écoulements, une approche dans laquelle le gaz est décrit par les équations de Navier-Stokes compressible et le liquide par les équations de Navier-Stokes incompressible est proposée. C'est le couplage de ces deux modèles qui constitue l'originalité et l'enjeu principal de de cette thèse. Pour traiter cette difficulté majeure, une méthode globale (i.e. la même dans chaque phase) et simple à mettre en \oe uvre est élaborée. L'utilisation des équations de Navier-Stokes formulées de façon unifiée pour les inconnues primitives (pression, vitesse et température) constitue le point de départ pour la construction de notre méthode qui repose sur les composants suivants: - une méthode d'éléments finis stabilisés pour la discrétisation spatiale des équations de Navier-Stokes; - une approche Level Set pour représenter précisément l'interface dont l'équation de transport a été résolue par une méthode de type Galerkin Discontinu; \item et des grandeurs moyennées pour traiter les discontinuités à l'interface. Le bon comportement de notre approche est illustré sur différents tests mono et bi-dimensionnels.

Dans cette thèse, nous nous intéressons à la simulation numérique d'écoulements instationnaires de deux fluides visqueux non miscibles, séparés par une interface mobile. Plus particulièrement des écoulements sans choc constitués d'une phase gazeuse et d'une phase liquide sont considérés. Pour modéliser de tels écoulements, une approche dans laquelle le gaz est décrit par les équations de Navier-Stokes compressible et le liquide par les équations de Navier-Stokes incompressible est proposée. C'est le couplage de ces deux modèles qui constitue l'originalité et l'enjeu principal de de cette thèse. Pour traiter cette difficulté majeure, une méthode globale (i.e. la même dans chaque phase) et simple à mettre en oeuvre est élaborée. L'utilisation des équations de Navier-Stokes formulées de façon unifiée pour les inconnues primitives (pression, vitesse et température) constitue le point de départ pour la construction de notre méthode qui repose sur les composants suivants: une méthode d'éléments finis stabilisés pour la discrétisation spatiale des équations de Navier-Stokes; une approche Level Set pour représenter précisément l'interface dont l'équation de transport a été résolue par une méthode de type Galerkin Discontinu; et des grandeurs moyennes pour traiter les discontinuités à l'interface. Le bon comportement de notre approche est illustré sur différents tests mono et bi-dimensionnels
In this work, we are interested in the numerical simulation of instationnary viscous flows of two immiscible fluids, separated by a mobile interface. In particular, flows without shock composed of a gas phase and a liquid phase are considered. In order to modelize such flows, an approach in which the gaz is described by compressible Navier-Stokes equations and the liquid by incompressible Navier-Stokes équations is proposed. The coupling between these two models is the originality and the stake of this thesis. To treat this important difficulty, a global (i.e. the same for each phase) and simple method is elaborated. In our procedure we propose, using the Navier-Stokes equations formulated in set of primitives unknowns (pressure, velocity and temperature), to elaborate a strategy that relies on the follow components: the stabilized finite element method to discretize spatially the Navier-Stokes equations; the Level Set method for tracking the interface precisely with a discontinuous Galerkin method to solve the associated transport equation; and some averaged quantities to treat the discontinuities at the interface. The good behaviour of this approach is performed on both one and two spatial dimensions