Дисертації з теми "Two-phase incompressible flows"

L'objectif de cette thèse est de développer une méthode précise d'ordre élevé pour résoudre le problème d'écoulementlaminaire incompressible à deux phases. Trois tâches principales sont à accomplir. Premièrement, la méthode doit être stable en énergie, ce qui signifie que la condition sans divergence de l'équation de Navier-Stokes incompressible est satisfaite partout dans le domaine de calcul. Deuxièmement, les discontinuités locales apparaissant dans le champ d'écoulement diphasique doivent être capturées avec précision. Troisièmement, l'interface matérielle entre les deux fluides doit être représentée avec précision à chaque pas de temps. Dans ce travail, une nouvelle méthode Hybridizable Discontinuous Galerkin (HDG) est utilisée pour la discrétisation spatiale. Cette méthode hybride qui appartient à la famille des méthodes DG-FEM satisfait la condition sans divergence en introduisant des variables de trace de vitesse et de pression du même ordre plus une approximation de vitesse et de pression adaptée à l'intérieur des éléments. Deplus, les concepts de FEM eXtended (X-FEM) sont utilisés pour approximer les discontinuités dans le champ d'écoulement en enrichissant l'approximation FEM standard dans les éléments où deux fluides existent. Enfin, l'interface du matériau en mouvement entre les deux fluides est capturée à l'aide de la méthode Level-Set
The objective of this thesis is to develop a high-order accurate method to solve the two-phase incompressible laminar flowproblem. Three main tasks are to be achieved. First, the method has to be energy-stable meaning that the divergence-free condition of the incompressible Navier-Stokes equation is satisfied everywhere in the computational domain. Second, the local discontinuities arising in the two-phase flow field have to be captured accurately. Third, the material interface betweenthe two fluids has to be represented accurately in each time step. In this work, a novel Hybridizable Discontinuous Galerkin (HDG) method is used for the spatial discretization. This hybrid method that belongs to the family of DG-FEM methods satisfies the divergence-free condition by introducing velocity and pressure trace variables of the same order plus a tailoredvelocity and pressure approximation inside the elements. Furthermore, the concepts of eXtended FEM (X-FEM) are used toapproximate discontinuities in the flow field by enriching the standard FEM approximation in elements where two fluids exist. Finally, the moving material interface between the twofluids is captured using the Level-Set method

Cette thèse est consacrée au développement et à la comparaison des méthodes de suivi d'interface pour les écoulements diphasiques incompressibles. Elle s'intéresse à la sélection de méthodes robustes de suivi d'interface, puis à leur couplage avec le solveur des équations de Navier-Stokes. La méthode level-set est en premier lieu étudiée, en particulier l'influence du schéma d'advection et de l'étape de réinitialisation sur la qualité des résultats du suivi d'interface. Il a été montré que la méthode de réinitialisation avec contrainte de volume est robuste et précise en combinaison avec des schémas conservatifs WENO d'ordre 5 pour l'advection. Il a été constaté que les erreurs du suivi d'interface augmentent de manière abrupte lorsque la condition CFL est trop petite. Comme remède, la réinitialisation du champ level-set effectuée moins souvent réduit la diffusion numérique et le déplacement non-physique de l'interface. La conservation de la masse n'est pas assurée avec les méthodes level-set. Les méthodes VOF (volume-of-fluid) qui conservent naturellement la masse du fluide de référence sont alors étudiées. Une résolution géométrique avec un schéma consistent et conservatif est alors adoptée, ainsi qu'une autre technique alternative plus aisément extensible en 3D. Il a été trouvé que ces deux dernières méthodes donnent des résultats très proches. La méthode MOF (moment-of-fluid), qui reconstruit l'interface en utilisant le centre de masse du fluide de référence, est plus précise que les méthodes VOF. Différentes méthodes couplées entre level-set et VOF sont alors étudiées, notamment: CLSVOF, MCLS, VOSET et CLSMOF. Il a été observé que la méthode level-set tend à épaissir les filaments minces, tandis que VOF et les méthodes couplées les fragmentent en petites particules. Finalement, on a couplé les méthodes level-set et VOF avec le solveur incompressible des équations de Navier-Stokes. On a comparé différentes manières de prise en compte des conditions de saut à l'interface (lisse et raide). Il a été montré que les méthodes VOF sont plus robustes, et donnent d'excellents résultats pour quasiment toutes les simulations. Deux méthodes level-set donnant de très bons résultats, comparables à ceux de VOF, sont aussi identifiées
This thesis is devoted to the development and comparison of interface methods for incompressible two-phase flows. It focuses on the selection of robust interface capturing methods, then on the manner of their coupling with the Navier-stokes solver. The level-set method is first investigated, in particular the influence of the advection scheme and the reinitialization step on the accuracy of the interface capturing. It is shown that the volume constraint method for reinitialization is robust and accurate in combination with the conservative fifth-order WENO schemes for the advection. It is found that interface errors increase drastically when the CFL number is very small. As a remedy, reinitializing the level-set field less often reduces the amount of numerical diffusion and non-physical interface displacement. Mass conservation is, however, not guaranteed with the level-set methods. The volume-of-fluid (VOF) method is then investigated, which naturally conserves the mass of the reference fluid. A geometrical consistent and conservative scheme is adopted, then an alternative technique more easily extended to 3D. It is found that both methods give very similar results. The moment-of-fluid (MOF) method, which reconstructs the interface using the reference fluid centroid, is found to be more accurate than the VOF methods. Different coupled level-set and VOF methods are then investigated, namely: CLSVOF, MCLS, VOSET and CLSMOF. It is observed that the level-set method tends to thicken thin filaments, whereas the VOF and coupled methods break up thin structures in small fluid particles. Finally, we coupled the level-set and volume-of-fluid methods with the incompressible Navier-Stokes solver. We compared different manners (sharp and smoothed) of treating the interface jump conditions. It is shown that the VOF methods are more robust, and provide excellent results for almost all the performed simulations. Two level-set methods are also identified that give very good results, comparable to those obtained with the VOF methods

The evolution of numerical methods and computational facilities allow re- searchers to explore complex physical phenomenons such as multiphase flows. The specific regime of incompressible, turbulent, bubbly two-phase flow (where a car- rier fluid is infused with bubbles or particles) is also receiving increased attention due to it’s appearance in major industrial processes. The main challenges arise from coupling individual aspects of the physics into a unified model and to provide a robust numerical framework. The presented work aimed at to achieve the second part by employing the most frequently used dispersed two-phase flow model and another incompressible, turbulent single phase solver as a base flow provider for coupled Lagrangian or surface tracking tools. Among the numerical techniques, the finite element method is a powerful can- didate when the need arises for multiphysics simulations (for example coupling with an electrochemical module) where the counterpart has a node based ap- proach. Stabilization schemes such as PSPG/SUPG/BULK provide remedies for the pressure decoupling and the inherent instability of the central discretization when applied for convective flow problems. As an alternative to unsteady solvers based upon an explicit or a fully im- plicit nonlinear treatment of the convective terms, a semi-implicit scheme results in a method of second order accurate in both space and time, has absolute linear stability and requires only a single or two linear system solution per time step. The application of the skew symmetric approach to the convective term further stabilizes the solution procedure and in some cases it even prevents divergence. The Eulerian-Eulerian two-phase flow model poses various issues to be over- come. The major difficulty is the density ratio between the phases; for an ordinary engineering problem it is in the order of thousands or more. The seemingly minus- cule differences in the formulation of the stabilizations can cause very different end results and require careful analysis. Volume fraction boundedness is of concern as well, but it is treatable by solving for its logarithm. Since the equations allow jumps (even separation of the phases) in the volume fraction field, discontinuity capturing techniques are also needed. Besides the standard ’spatial’ stabilization temporal smoothing is also necessary, otherwise the limitation in time step size becomes too stringent. Designing a flow solver is one side of the adventure, but verification is equally important. Comparison against analytical solution (such as the single and two- phase Taylor-Green testcase) provides insight and confirmation about the mathe- matical and physical properties. Meanwhile comparing with real life experiments prove the industrialization and usability of a code, dealing with low quality meshes and effective utilization of computer clusters.
Doctorat en Sciences de l'ingénieur et technologie
info:eu-repo/semantics/nonPublished

The level set method is a powerful way of tracking surfaces by defining the surface as a zero level set of a continuous function that is usually a signed distance function. The level set method is one of the best methods for simulating multi-phase flow because it can easily handle fast topological changes, as well as splitting and merging of fluids. In this thesis, a standard level set method was implemented in C++, using the finite element method library deal.II, to simulate incompressible two-phase flow on some benchmark problems. The results show a significant change of mass in the simulations, something that should not be allowed to happen when simulating incompressible fluids. The mass changes mainly occur in the reinitialization phase, where the level set function is rebuilt to look more like a signed distance function.

This item is only available in print in the UCF Libraries. If this is your Honors Thesis, you can help us make it available online for use by researchers around the world by following the instructions on the distribution consent form at http://library.ucf.edu/Systems/DigitalInitiatives/DigitalCollections/InternetDistributionConsentAgreementForm.pdf You may also contact the project coordinator, Kerri Bottorff, at kerri.bottorff@ucf.edu for more information.
Bachelors
Engineering and Computer Science
Mechanical Engineering

PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO DE JANEIRO
CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICO
Esta tese propõe um método meshfree adaptativo de advecção para problemas de fluxo bifásico de fluidos incompressíveis e imiscíveis em meios porosos heterogêneos tridimensionais. Este método se baseia principalmente na combinação do método Semi-Lagrangeano adaptativo com interpolação local meshfree usando splines poliharmônicas como funções de base radial. O método proposto é uma melhoria e uma extensão do método adaptativo meshfree AMMoC proposto por Iske e Kaser (2005) para modelagem 2D de reservatórios de petróleo. Inicialmente este trabalho propõe um modelo em duas dimensões, contribuindo com uma melhoria significativa no cálculo do Laplaciano, utilizando os métodos meshfree de Hermite e Kansa. Depois, o método é ampliado para três dimensões (3D) e para um meio poroso heterogêneo. O método proposto é testado com o problema de five spot e os resultados são comparados com os obtidos por sistemas bem conhecidos na indústria de petróleo.
This thesis proposes an adaptive meshfree advection method for two-phase flow problems of incompressible and immiscible fluids through three-dimensional heterogeneous porous media. This method is based mainly on a combination of adaptive semi-Lagrangian method with local meshfree interpolation using polyharmonic splines as radial basis functions. The proposed method is an improvement and extension of the adaptive meshfree advection scheme AMMoC proposed by Iske and Kaser (2005) for 2D oil reservoir modeling. Initially this work proposes a model in two dimensions, contributing to a significant improvement in the calculation of the Laplacian, using the meshfree methods of Hermite and Kansa. Then, the method is extended to three dimensions (3D) and a heterogeneous porous medium. The proposed method is tested with the five spot problem and the results are compared with those obtained by well-known systems in the oil industry.

Dans cette thèse, on s'intéresse à deux problèmes provenant de l'étude mathématique des fluides incompressibles visqueux : la propagation de la régularité tangentielle et le mouvement d'une surface libre.La première question concerne plus particulièrement l'étude qualitative de l'évolution de quantités thermodynamiques telles que la température dans l'équation de Boussinesq sans diffusion et la densité dans le système de Navier-Stokes non homogène. Typiquement, on suppose que ces deux quantités sont, à l'instant initial, discontinues le long d'une interface à régularité h"oldérienne. Comme conséquence de résultats de propagation de régularité tangentielle pour le champ de vitesses, on établit que la régularité des interfaces persiste pour tout temps aussi bien en dimension deux d'espace, qu'en dimension supérieure (avec condition de petitesse). Notre approche suit celle du travail de J.-Y. Chemin dans les années 90 pour le problème des poches de tourbillon dans les fluides incompressiblesparfaits.Dans le cas présent, outre cette hypothèse de régularité tangentielle, nous n'avons besoin que d'une régularité critique sur le champ de vitesses.La démonstration repose sur le calcul para-différentiel et les espaces de multiplicateurs.Dans la dernière partie de la thèse, on considère le problème à frontière libre pour le système de Navier-Stokes incompressible à deux phases. Ce système permet de décrire l'évolution d'un mélange de deux fluides non miscibles tels que l'huile et l'eau par exemple. Différents cas de figure sont étudiés : le cas d'un réservoir borné, d'une goutte ou d'une rivière à profondeur finie.On établit l'existence et l'unicité à temps petit pour ce problème. Notre démonstration repose fortement sur des propriétés de régularité maximale parabolique de type $L_p$-$L_q
This thesis is dedicated to two different problems in the mathematical study of the viscous incompressible fluids: the persistence of tangential regularity and the motion of a free surface.The first problem concerns the study of the qualitative properties of some thermodynamical quantities in incompressible fluid models, such as the temperature for Boussinesq system with no diffusion and the density for the non-homogeneous Navier-Stokes system. Typically, we assume those two quantities to be initially piecewise constant along an interface with H"older regularity.As a consequence of stability of certain directional smoothness of the velocity field, we establish that the regularity of the interfaces persist globally with respect to time both in the two dimensional and higher dimensional cases (under some smallness condition). Our strategy is borrowed from the pioneering works by J.-Y.Chemin in 1990s on the vortex patch problem for ideal fluids.Let us emphasize that, apart from the directional regularity, we only impose rough (critical) regularity on the velocity field. The proof requires tools from para-differential calculus and multiplier space theory.In the last part of this thesis, we are concerned with the free boundary value problem for two-phase density-dependent Navier-Stokes system.This model is used to describe the motion of two immiscible liquids, like the oil and the water. Such mixture may occur in different situations, such as in a fixed bounded container, in a moving bounded droplet or in a river with finite depth. We establish the short time well-posedness for this problem. Our result strongly relies on the $L_p$-$L_q$ maximal regularity theoryfor parabolic equations

Cette thèse présente la Simulation des Grandes Echelles (LES) de l’injection, de la pulvérisation et de la cavitation dans un injecteur pour les applications liées aux moteurs à combustion interne. Pour la modélisation de l’atomisation, on utilise le modèle ELSA (Eulerian Lagrangian Spray Atomization). Le modèle résout la fraction volumique du combustible liquide ainsi que la densité de surface d’interface liquide-gaz pour décrire le processus complet d’atomisation. Dans cette thèse, l’écoulement à l’intérieur de l’injecteur est également pris en compte pour une étude ultérieure de l’atomisation. L’étude présente l’application du modèle ELSA à un injecteur Diesel typique, à la fois dans le contexte de RANS et de LES.Le modèle est validé à l’aide de données expérimentales disponibles dans Engine Combustion Network (ECN). Le modèle ELSA, qui est normalement conçu pour les interfaces diffuses (non résolues), lorsque l’emplacement exact de l’interface liquide-gaz n’est pas pris en compte, est étendu pour fonctionner avec une formulation de type Volume of Fluid (VOF) de flux à deux phases, où l’interface est explicitement résolu. Le couplage est réalisé à l’aide de critères IRQ (Interface Resolution Quality), qui prennent en compte à la fois la courbure de l’interface et la quantité modélisée de la surface de l’interface. Le modèle ELSA est développé en premier lieu en considérant les deux phases comme incompressibles. L’extension à la phase compressible est également brièvement étudiée dans cette thèse. Il en résulte une formulation ELSA compressible qui prend en compte la densité variable de chaque phase. En collaboration avec l’Imperial College de Londres, la formulation de la fonction de densité de probabilité (PDF) avec les champs stochastiques est également explorée afin d’étudier l’atomisation. Dans les systèmes d’injection de carburant modernes, la pression locale à l’intérieur de l’injecteur tombe souvent en dessous de la pression de saturation en vapeur du carburant, ce qui entraîne une cavitation. La cavitation affecte le flux externe et la formulation du spray. Ainsi, une procédure est nécessaire pour étudier le changement de phase ainsi que la formulation du jet en utilisant une configuration numérique unique et cohérente. Une méthode qui couple le changement de phase à l’intérieur de l’injecteur à la pulvérisation externe du jet est développée dans cette thèse. Ceci est réalisé en utilisant le volume de formulation de fluide où l’interface est considérée entre le liquide et le gaz; le gaz est composé à la fois de vapeur et d’airambiant non condensable
This thesis presents Large Eddy Simulation (LES) of fuel injection, atomization and cavitation inside the fuel injector for applications related to internal combustion engines. For atomization modeling, Eulerian Lagrangian Spray Atomization (ELSA) model is used. The model solves for volume fraction of liquid fuel as well as liquid-gas interface surface density to describe the complete atomization process. In this thesis, flow inside the injector is also considered for subsequent study of atomization. The study presents the application of ELSA model to a typical diesel injector, both in the context of RANS and LES. The model is validated with the help of experimental data available from Engine Combustion Network (ECN). The ELSA model which is normally designed for diffused (unresolved) interfaces, where the exact location of the liquid-gas interface is not considered, is extended to work with Volume of Fluid (VOF) type formulation of two phase flow, where interface is explicitly resolved. The coupling is achieved with the help of Interface Resolution Quality (IRQ) criteria, that takes into account both the interface curvature and modeled amount of interface surface. ELSA model is developed first considering both phases as incompressible, the extension to compressible phase is also briefly studied in this thesis, resulting in compressible ELSA formulation that takes into account varying density in each phase. In collaboration with Imperial College London, the Probability Density Function (PDF) formulation with Stochastic Fields is also explored to study atomization. In modern fuel injection systems, quite oftenthe local pressure inside the injector falls below the vapor saturation pressure of the fuel, resulting in cavitation. Cavitation effects the external flow and spray formulation. Thus, a procedure is required to study the phase change as well as jet formulation using a single and consistent numerical setup. A method is developed in this thesis that couples the phase change inside the injector to the external jet atomization. This is achieved using the volume of fluid formulation where the interface is considered between liquid and gas; gas consists of both the vapor and non condensible ambient air

碩士
國立交通大學
應用數學系所
97
The mathematical model of the waterflood problem which is applied in this paper can be divided into two sections. One is the pressure equation and the other is the saturation equation. And the saturation equation also can be pa- rtitioned into the transport stage saturation and the diffusive stage saturation. However, we will pay more attention to solve the transport stage saturation in this research. Here we construct a meters reservoir system for simu- lation. An efficient numerical method, locally conservative Eulerian-Lagrangian methods (LCELM), is developed to compute the transport equation to improve the conservation of waterflood. From the results of the numerical simulations, we can realize the relation between temporal variation and the flow condition.

碩士
國立交通大學
應用數學系所
100
The primary objective of this thesis is to introduce a multigrid method to solve elliptic equation with strongly discontinuous coefficients. In the beginning, we explain how to use the multigrid method to solve a 3D elliptic equation with strongly discontinuous coefficients, and then show some numerical testing results. Also, we provide some results compared with other numerical methods to show the efficency of the mutigrid method. Furthermore, we apply the multigrid method to solve two mathematical problems, one is for the waterflooding problem and the other is the incompressible Navier-Stokes equations. A locally conservative Eulerian-Lagrangian method (briefly LCELM) is used to compute the transport part of the two models. Some numerical results for the two problems will be presented as well. ii

abstract: The goal of this paper was to do an analysis of two-dimensional unsplit mass and momentum conserving Finite Volume Methods for Advection for Volume of Fluid Fields with interfaces and validating their rates of convergence. Specifically three unsplit transport methods and one split transport method were amalgamated individually with four Piece-wise Linear Reconstruction Schemes (PLIC) i.e. Unsplit Eulerian Advection (UEA) by Owkes and Desjardins (2014), Unsplit Lagrangian Advection (ULA) by Yang et al. (2010), Split Lagrangian Advection (SLA) by Scardovelli and Zaleski (2003) and Unsplit Averaged Eulerian-Lagrangian Advection (UAELA) with two Finite Difference Methods by Parker and Youngs (1992) and two Error Minimization Methods by Pilliod Jr and Puckett (2004). The observed order of accuracy was first order in all cases except when unsplit methods and error minimization methods were used consecutively in each iteration, which resulted in second-order accuracy on the shape error convergence. The Averaged Unsplit Eulerian-Lagrangian Advection (AUELA) did produce first-order accuracy but that was due to a temporal error in the numerical setup. The main unsplit methods, Unsplit Eulerian Advection (UEA) and Unsplit Lagrangian Advection (ULA), preserve mass and momentum and require geometric clipping to solve two-phase fluid flows. The Unsplit Lagrangian Advection (ULA) can allow for small divergence in the velocity field perhaps saving time on the iterative solver of the variable coefficient Poisson System.
Dissertation/Thesis
Masters Thesis Mechanical Engineering 2019