Divergence-Free Mixed Finite Elements for the Incompressible Navier-Stokes Equation

Language
en
Document Type
Doctoral Thesis
Issue Date
2008-04-08
Issue Year
2007
Authors
Linke, Alexander
Editor
Abstract

This dissertation is concerned with the numerical approximation of the incompressible Navier-Stokes equation. As our main contribution, we propose and analyze a new stabilized mixed finite element scheme for computing incompressible laminar flows, and investigate several properties of this new scheme. The finite element analysis is presented for the Oseen model problem. The scheme is based on the Scott-Vogelius mixed finite element and on symmetric stabilization operators for dominant convection, proposed in the last recent years. In particular, the scheme delivers divergence-free pointwise velocity approximations, with an approximation quality that is completely independent of the pressure. Therefore, the cumbersome grad-div stabilization drops out from the discrete equations, and standard multi-grid approaches are possible for the solution of the evolving linear systems. We propose such a multi-grid method and deliver a multi-grid analysis for the full symmetric part of the discrete Oseen equation, including the symmetric stabilization operator for dominant convection. We further propose a coupled FEM-FVM scheme for convection-diffusion problems in an incompressible flow field. For the numerical computation of a stationary incompressible Navier-Stokes equation we propose the above stabilized Scott-Vogelius scheme, and for a stationary or nonstationary scalar convection-diffusion equation, we propose a Voronoi-box-based finite volume scheme on boundary conforming Delaunay meshes. Since Scott-Vogelius finite element approximations are divergence-free pointwise, we can prove a discrete maximum principle for the discrete convection-diffusion equation. Finally, we present several numerical examples, in order to illustrate in which situations the proposed stabilized Scott-Vogelius scheme delivers more accurate numerical approximations than other approaches. Here, a 2D example of a colliding flow might deserve further interest, since it seems to have some physical relevance. Further, we present a 3D application from electrochemistry, where the coupled FEM-FVM scheme seems to be promising due to the establishment of a discrete maximum principle.

Abstract

Diese Dissertation befasst sich mit der numerischen Approximation der inkompressiblen Navier-Stokes-Gleichung. Als Hauptbeitrag der Arbeit führen wir eine neue, stabilisierte gemischte Finite-Elemente-Diskretisierung für die numerische Behandlung von inkompressiblen laminaren Strömungen ein, analysieren sie, und untersuchen einige Eigenschaften dieser Diskretisierung. In der Finite-Element-Analysis untersuchen wir das lineare Oseen-Modell-problem. Die Diskretisierung basiert auf dem Scott-Vogelius-Element, und benutzt symmetrische Stabilisierungsoperatoren, um die diskreten Oseen-Gleichungen gegenüber dominanter Konvektion zu stabilisieren. Eine Besonderheit der Diskretisierung besteht darin, dass sie punktweise divergenz-freie Geschwindigkeitsapproximationen liefert, und damit die Approximationsgüte der Geschwindigkeiten von der Druckapproximation völlig entkoppelt. Aufgrund der Divergenzfreiheit fällt die lästige Grad-Div-Stabilisierung aus den diskreten Oseen-Gleichungen heraus, und Standard-Mehrgitter-Methoden sind anwendbar, um die entstehenden linearen Gleichungssysteme zu lösen. Wir stellen eine solche Standard-Mehrgitter-Methode vor, und entwickeln eine Mehrgitter-Analysis für den gesamten symmetrischen Teil der diskreten Oseen-Gleichung, inklusive dem symmetrischen Stabilisierungsoperator. Weiter stellen wir eine gekoppelte Finite-Volumen/Finite-Element-Diskretisierung vor, mit der Konvektions-Diffusions-Probleme in einem inkompressiblen Strömungsfeld behandelt werden können. Für die numerische Approximation der Navier-Stokes-Gleichung wird das vorgeschlagene stabilisierte Scott-Vogelius-Element verwendet, und für die gekoppelte skalare Konvektions-Diffusions-Gleichung setzen wir eine Voronoi-Box-basierte Finite-Volumen-Methode auf randkonformen Delaunaygittern ein. Da das Scott-Vogelius-Element punktweise divergenz-freie Approximationen liefert, können wir ein lokales diskretes Maximumprinzip für die angekoppelte skalare Konvektions-Diffusions-Diskretisierung beweisen. Endlich stellen wir einige numerische Beispiele vor, um veranschaulichen zu können, in welchen Situationen das vorgestellte Scott-Vogelius-Element numerisch genauere Approximation liefert als andere Ansätze. Besonders scheint ein 2D-Beispiel eines kollidierenden Fluids Beachtung zu verdienen, da es vergleichsweise physikalisch bedeutsam ist. Weiter stellen wir Ergebnisse einer 3D-Anwendungsrechnung aus der Elektrochemie vor, bei der das gekoppelte Finite-Volumen/Finite-Element-Schema aufgrund des diskreten Maximumsprinzips eine interessante Diskretisierungsvariante zu sein scheint.

DOI
Document's Licence
Faculties & Collections
Zugehörige ORCIDs