TY - THES A1 - Barsukow, Wasilij T1 - Low Mach number finite volume methods for the acoustic and Euler equations T1 - Finite Volumen Methoden für den Grenzwert niedriger Machzahlen der akustischen und der Euler-Gleichungen N2 - Finite volume methods for compressible Euler equations suffer from an excessive diffusion in the limit of low Mach numbers. This PhD thesis explores new approaches to overcome this. The analysis of a simpler set of equations that also possess a low Mach number limit is found to give valuable insights. These equations are the acoustic equations obtained as a linearization of the Euler equations. For both systems the limit is characterized by a divergencefree velocity. This constraint is nontrivial only in multiple spatial dimensions. As the Jacobians of the acoustic system do not commute, acoustics cannot be reduced to some kind of multi-dimensional advection. Therefore first an exact solution in multiple spatial dimensions is obtained. It is shown that the low Mach number limit can be interpreted as a limit of long times. It is found that the origin of the inability of a scheme to resolve the low Mach number limit is the lack a discrete counterpart to the limit of long times. Numerical schemes whose discrete stationary states discretize all the analytic stationary states of the PDE are called stationarity preserving. It is shown that for the acoustic equations, stationarity preserving schemes are vorticity preserving and are those that are able to resolve the low Mach limit (low Mach compliant). This establishes a new link between these three concepts. Stationarity preservation is studied in detail for both dimensionally split and multi-dimensional schemes for linear acoustics. In particular it is explained why the same multi-dimensional stencils appear in literature in very different contexts: These stencils are unique discretizations of the divergence that allow for stabilizing stationarity preserving diffusion. Stationarity preservation can also be generalized to nonlinear systems such as the Euler equations. Several ways how such numerical schemes can be constructed for the Euler equations are presented. In particular a low Mach compliant numerical scheme is derived that uses a novel construction idea. Its diffusion is chosen such that it depends on the velocity divergence rather than just derivatives of the different velocity components. This is demonstrated to overcome the low Mach number problem. The scheme shows satisfactory results in numerical simulations and has been found to be stable under explicit time integration. N2 - Finite Volumen Methoden für die kompressiblen Euler-Gleichungen zeigen übermäßige Diffusion im Grenzwert kleiner Machzahlen. Diese Dissertation beschäftigt sich mit neuen Ansätzen, um dieses Problem zu beheben. Die Analyse eines Systems einfacherer Gleichungen, die ebenso einen Grenzwert niedriger Machzahlen haben, liefert wichtige Einsichten. Diese Gleichungen sind die als Linearisierung der Euler-Gleichungen erhaltenen akustischen Gleichungen. Für beide Gleichungssysteme ist der Grenzwert durch ein divergenzfreies Geschwindigkeitsfeld charakterisiert, was nur in mehreren Raumdimensionen nichttrivial ist. Da die Jacobi-Matrizen des akustischen Systems nicht vertauschen, kann Akustik nicht auf irgendeine Art mehrdimensionaler Advektion zurückgeführt werden. Deswegen wird zunächst eine exakte Lösung in mehreren Raumdimensionen gefunden. Es wird gezeigt, dass sich der Grenzwert kleiner Machzahlen als Grenzwert langer Zeiten interpretieren lässt. Als der Ursprung des Versagens eines Schemas im Grenzwert kleiner Machzahlen wird das Fehlen einer diskreten Entsprechung zum Grenzwert langer Zeiten identifiziert. Numerische Schemata, deren diskrete stationäre Zustände alle analytischen stationären Zustände diskretisieren, werden stationaritätserhaltend genannt. Es zeigt sich, dass für die akustischen Gleichungen stationaritätserhaltende Schemata vortizitätserhaltend sind, und gerade diejenigen sind, die auch den Grenzwert kleiner Machzahlen aufzulösen vermögen. Das zeigt eine neue Verbindung zwischen diesen drei Konzepten auf. Erhaltung der Stationarität wird für lineare Akustik im Detail für Schemata studiert, die nach Raumdimensionen aufgeteilt sind, und auch für multi-dimensionale Schemata. Insbesondere wird ein Grund geliefert, warum die gleichen multi-dimensionalen diskreten Operatoren in der Literatur in sehr unterschiedlichen Kontexten auftauchen: Sie sind Diskretisierungen der Divergenz, für die eine stabilisierende, stationaritätserhaltende Diffusion gefunden werden kann. Auch für nichtlineare Gleichungen, wie die Euler-Gleichungen, kann die Erhaltung der Stationarität verallgemeinert werden. Es werden dazu mehrere Wege der Konstruktion numerischer Schemata gezeigt. Insbesondere im Hinblick auf den Grenzwert kleiner Machzahlen wird ein neuartiges Schema hergeleitet, dessen Diffusion so gewählt ist, dass es von der Divergenz der Geschwindigkeit, und nicht bloß von irgendswelchen Ableitungen der Geschwindigkeitskomponenten abhängt. Es wird gezeigt, dass dieses Schema in der Lage ist, den Grenzwert kleiner Machzahlen aufzulösen. Das Schema zeigt zufriedenstellende Resultate in Simulationen und ist stabil unter Verwendung eines expliziten Zeitintegrators. KW - Finite-Volumen-Methode KW - Machzahl KW - finite volume method KW - Euler equations KW - Acoustic equations KW - low Mach number KW - vorticity preserving Y1 - 2018 UR - https://opus.bibliothek.uni-wuerzburg.de/frontdoor/index/index/docId/15996 UR - https://nbn-resolving.org/urn:nbn:de:bvb:20-opus-159965 ER -