Discontinuous Galerkin methods for the unsteady compressible Navier-Stokes equations

Abstract

In this work a new explicit arbitrary high order accurate discontinuous Galerkin finite element solver for the unsteady compressible Navier-Stokes equations is developed. Although the focus is on the compressible Navier-Stokes equations, the developed framework can directly be applied to other pure hyperbolic, pure parabolic or mixed hyperbolic/parabolic time dependent conservation laws.

Discontinuous Galerkin finite element based methods have several important properties. They are locally conservative schemes, despite their affiliation to the class of finite element methods. They allow to use arbitrary unstructured non-conforming meshes, while remaining their formal (high) order of accuracy, even with skewed and anisotropic grid cells. The resolution can be adapted locally by increasing or decreasing the local polynomial degree, without the difficulties of a conforming finite element approach. The (formal) order of accuracy is essential a parameter, as one only has to choose the polynomial degree of the approximation. The most important property of discontinuous Galerkin schemes is that the solution in a grid cell depends only on data from grid cells sharing a face, independent of the approximation order. This compactness yields an excellent parallelizability of the method, which is essential for large scale computations.

Based on the standard spatial discontinuous Galerkin framework several modifications to increase the computational efficiency are proposed. In a first step a novel construction guideline for modal and nodal basis functions on arbitrary shaped grid cells is introduced. For the discretization of problems with high order derivatives a novel weak formulation is introduced and applied to the second order compressible Navier-Stokes equations. For the approximation of the viscous fluxes a new approximation based on local Riemann solutions is used. This spatial discretization is combined with a new time discretization, which allows time consistent local time stepping.

In order to validate the high accuracy and efficiency of the developed method, several test cases including the linearized Euler equations, the non-linear Euler equations and the compressible Navier-Stokes equations are considered. Finally, this method is applied to solve two and three dimensional compressible unsteady flow problems.

In dieser Arbeit wird ein neues explizites unstetiges Finite Element Verfahren von beliebig hoher Genauigkeit für die Lösung der zeitabhängigen kompressiblen Navier-Stokes Gleichungen entwickelt. Obwohl der Fokus dieser Arbeit auf dem Lösen der kompressiblen Navier-Stokes Gleichung liegt, kann man die entwickelte Methode auch direkt für die Lösung von anderen rein hyperbolischen, rein parabolischen oder gemischt hyberbolisch/parabolischen zeitabhängigen Erhaltungsgleichungen verwenden.

Die unstetige Finite Elemente Methode zeichnet sich durch eine Reihe von wichtigen Eigenschaften aus. So ist sie, obwohl sie zur Klasse der Finite Elemente Verfahren gehört, lokal konservativ. Zudem ist es möglich, allgemeine unstrukturierte Netze mit stark verzerrten anisotropen Gitterzellen zu verwenden, ohne die formale Genauigkeit zu verlieren. Die formale Genauigkeitsordnung des Verfahrens lässt sich beliebig einstellen, da man im Prinzip nur den Polynomgrad ändern muss. Die wichtigste Eigenschaft des unstetigen Finite Element Verfahrens ist, dass die Lösung einer Gitterzelle nur von den Lösungen der Nachbargitterzellen abhängt. Diese Kompaktheit erlaubt eine hohe Parallelisierbarkeit der Methode, was für große Rechnungen essentiell ist.

Basierend auf der Standard unstetigen Finite Elemente Methodik werden verschiedene Modifikationen vorgeschlagen, um die Recheneffizienz zu erhöhen. Zuerst wird die Konstruktion von neuartigen modalen und nodalen Basisfunktionen für beliebige Gitterzellen eingeführt. Für die Navier-Stokes Gleichungen wird eine neuartige schwache Formulierung entwickelt. Die neuartige Approximation der viskosen Flüsse basiert auf einer lokalen Lösung des Riemannproblems. Diese räumliche Diskretisierung wird zudem mit einer neuen Zeitdiskretisierung kombiniert, welche zeitgenaue lokale Zeitschritte erlaubt.

Für die Validierung des entwickelten Verfahrens werden verschiedene Testfälle betrachtet, mit Hilfe derer die hohe Genauigkeit und Effizienz des Verfahrens gezeigt wird. Dieses Verfahren wird dann zur Lösung von zwei- und drei-dimensionalen kompressiblen Strömungsproblemen verwendet.