TY - THES A1 - Klotzky, Jens T1 - Well-posedness of a fluid-particle interaction model T1 - Existenz und Eindeutigkeit von Entropielösungen eines Partikel-Fluid-Modells N2 - This thesis considers a model of a scalar partial differential equation in the presence of a singular source term, modeling the interaction between an inviscid fluid represented by the Burgers equation and an arbitrary, finite amount of particles moving inside the fluid, each one acting as a point-wise drag force with a particle related friction constant. \begin{align*} \partial_t u + \partial_x (u^2/2) &= \sum_{i \in N(t)} \lambda_i \Big(h_i'(t)-u(t,h_i(t)\Big)\delta(x-h_i(t)) \end{align*} The model was introduced for the case of a single particle by Lagoutière, Seguin and Takahashi, is a first step towards a better understanding of interaction between fluids and solids on the level of partial differential equations and has the unique property of considering entropy admissible solutions and the interaction with shockwaves. The model is extended to an arbitrary, finite number of particles and interactions like merging, splitting and crossing of particle paths are considered. The theory of entropy admissibility is revisited for the cases of interfaces and discontinuous flux conservation laws, existing results are summarized and compared, and adapted for regions of particle interactions. To this goal, the theory of germs introduced by Andreianov, Karlsen and Risebro is extended to this case of non-conservative interface coupling. Exact solutions for the Riemann Problem of particles drifting apart are computed and analysis on the behavior of entropy solutions across the particle related interfaces is used to determine physically relevant and consistent behavior for merging and splitting of particles. Well-posedness of entropy solutions to the Cauchy problem is proven, using an explicit construction method, L-infinity bounds, an approximation of the particle paths and compactness arguments to obtain existence of entropy solutions. Uniqueness is shown in the class of weak entropy solutions using almost classical Kruzkov-type analysis and the notion of L1-dissipative germs. Necessary fundamentals of hyperbolic conservation laws, including weak solutions, shocks and rarefaction waves and the Rankine-Hugoniot condition are briefly recapitulated. N2 - Diese Arbeit befasst sich mit dem Modell einer skalaren partiellen Differentialgleichung mit singulärem Quellterm, das die Interaktion zwischen einem reibungsfreiem Fluid, dargestellt durch die Burgers Gleichung, und einer gegebenen, endlichen Menge von sich in dem Fluid bewegenden Partikeln beschreibt, die eine punktweise Zugkraft auf das Fluid auswirken und durch eine entsprechende Reibungskonstante charakterisiert sind. \begin{align*} \partial_t u + \partial_x (u^2/2) &= \sum_{i \in N(t)} \lambda_i \Big(h_i'(t)-u(t,h_i(t)\Big)\delta(x-h_i(t)) \end{align*} Das Modell wurde für den Fall der Interaktion mit einem einzelnen Partikel durch Lagoutière, Seguin and Takahashi eingeführt, stellt einen ersten Schritt zu einem besseren Verständnis der Interaktion zwischen einem Fluid und Festkörpern auf dem Level der partiellen Differentialgleichungen dar und hat die einzigartige Eigenschaft, dass Entropielösungen und die Interaktion mit Schockwellen berücksichtigt werden. Das Modell wird zu einer beliebigen, endlichen Anzahl von Partikeln erweitert und Interaktionen wie das Verschmelzen und Spaltung von Partikeln werden behandelt. Existierende Theorie der Entropie-Zulässigkeit im Hinblick auf Interfaces und Erhaltungsgleichungen mit unstetiger Flussfunktion wird zusammengefasst, die Resultate werden verglichen und für die Regionen mit Partikelinteraktionen angepasst. Zu diesem Zweck wird die Theorie der Germs, eingeführt von Andreianov, Karlsen und Risebro, auf den vorliegenden Fall eines nicht-erhaltenden Interfaces erweitert. Für das Riemann Problem von auseinanderdriftenden Partikeln werden die exakten Lösungen berechnet und eine Analyse des Verhaltens von Entropielösungen über die von den Partikeln erzeugten Interface wird genutzt, um ein physikalisch sinnvolles und mit der Theorie eines einzelnen Partikels konsistentes Verhalten beim Verschmelzen und Spalten von Partikeln herzuleiten. Mit Hilfe einer expliziten Konstruktionsmethode, hergeleiteten L-infinity Beschränkungen, einer Approximation der Partikelpfade und Kompaktheitsargumenten wird gezeigt, dass das entsprechende Cauchy Problem wohlgestellt ist. Eindeutigkeit im Raum der schwachen Entropielösungen wird mit beinahe klassischen Argumenten der Theorie von Kruzkov sowie der Theorie von L1-dissipativen Germs gezeigt. Notwendige Grundlagen zu hyperbolischen Erhaltungsgleichungen, unter anderem die Theorie schwacher Lösungen, Schock- und Verdünnungswellen sowie die Rankine-Hugoniot Bedingung, werden in Grundzügen am Anfang der Arbeit wiederholt. KW - Hyperbolische Differentialgleichung KW - Entropielösung KW - Fluid-Partikel-Strömung KW - Burgers-Gleichung KW - Korrekt gestelltes Problem KW - Existenz und Eindeutigkeit KW - Entropiebedingung KW - Well-posedness KW - Entropy admissibility condition Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-169009 ER - TY - THES A1 - Gallego Valencia, Juan Pablo T1 - On Runge-Kutta discontinuous Galerkin methods for compressible Euler equations and the ideal magneto-hydrodynamical model T1 - Runge-Kutta Discontinuous-Galerkin Verfahren für die kompressiblen Euler Gleichungen und das ideale magnetohydrodynamische Modell N2 - An explicit Runge-Kutta discontinuous Galerkin (RKDG) method is used to device numerical schemes for both the compressible Euler equations of gas dynamics and the ideal magneto- hydrodynamical (MHD) model. These systems of conservation laws are known to have discontinuous solutions. Discontinuities are the source of spurious oscillations in the solution profile of the numerical approximation, when a high order accurate numerical method is used. Different techniques are reviewed in order to control spurious oscillations. A shock detection technique is shown to be useful in order to determine the regions where the spurious oscillations appear such that a Limiter can be used to eliminate these numeric artifacts. To guarantee the positivity of specific variables like the density and the pressure, a positivity preserving limiter is used. Furthermore, a numerical flux, proven to preserve the entropy stability of the semi-discrete DG scheme for the MHD system is used. Finally, the numerical schemes are implemented using the deal.II C++ libraries in the dflo code. The solution of common test cases show the capability of the method. N2 - Ein explizite Runge-Kutta discontinous Galerkin (RKDG) Verfahren wird angewendet, um numerische Diskretisierungen, sowohl für die kompressiblen Eulergleichungen der Gasdynamik, als auch für die idealen Magnetohydrodynamik (MHD) Gleichungen zu entwickeln. Es ist bekannt, dass diese System von Erhaltungsgleichungen unstetige Lösungen besitzen. Unstetigkeiten sind die Quelle von störenden Oszillationen im Lösungsprofil der numerischen Näherung, wenn ein numerisches Verfahren von hoher Ordnung verwendet wird. Verschiedene Techniken werden miteinander verglichen um störende Oszillationen zu kontrollieren, die bei der Approximation von Unstetigkeiten in der Lösung auftreten. Ein Verfahren zur Lokalisierung von Schockwellen wird vorgestellt und es wird gezeigt, dass dieses Verfahren nützlich ist um Regionen, in denen störende Oszillationen auftreten, zu bestimmen, so dass ein Limiter verwendet werden kann um diese numerischen Artefakte zu eliminieren. Um die Positivität spezieller Variablen, wie die Dichte und den Druck, zu bewahren, wird ein spezieller „positivitätserhaltender“ Limiter verwendet. Des Weiteren wird ein numerischer Fluss, für den bewiesenermaßen das semi-diskrete DG Verfahren für das MHD System Entropie-Stabil ist, verwendet. Abschließend werden die numerischen Verfahren unter Verwendung der deal.II C++ Bibliotheken im dflo code implementiert. Simulationen bekannter Testbeispiele zeigen das Potential dieses numerischen Verfahrens. KW - explicit discontinuous Galerkin KW - conservation laws KW - numerical methods KW - Euler equations KW - MHD KW - Eulersche Differentialgleichung KW - Galerkin-Methode KW - Numerisches Verfahren Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-148874 ER - TY - THES A1 - Schnücke, Gero T1 - Arbitrary Lagrangian-Eulerian Discontinous Galerkin methods for nonlinear time-dependent first order partial differential equations T1 - Arbitrary Lagrangian-Eulerian Discontinous Galerkin-Methode für nichtlineare zeitabhängige partielle Differentialgleichungen erster Ordnung N2 - The present thesis considers the development and analysis of arbitrary Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) methods with time-dependent approximation spaces for conservation laws and the Hamilton-Jacobi equations. Fundamentals about conservation laws, Hamilton-Jacobi equations and discontinuous Galerkin methods are presented. In particular, issues in the development of discontinuous Galerkin (DG) methods for the Hamilton-Jacobi equations are discussed. The development of the ALE-DG methods based on the assumption that the distribution of the grid points is explicitly given for an upcoming time level. This assumption allows to construct a time-dependent local affine linear mapping to a reference cell and a time-dependent finite element test function space. In addition, a version of Reynolds’ transport theorem can be proven. For the fully-discrete ALE-DG method for nonlinear scalar conservation laws the geometric conservation law and a local maximum principle are proven. Furthermore, conditions for slope limiters are stated. These conditions ensure the total variation stability of the method. In addition, entropy stability is discussed. For the corresponding semi-discrete ALE-DG method, error estimates are proven. If a piecewise $\mathcal{P}^{k}$ polynomial approximation space is used on the reference cell, the sub-optimal $\left(k+\frac{1}{2}\right)$ convergence for monotone fuxes and the optimal $(k+1)$ convergence for an upwind flux are proven in the $\mathrm{L}^{2}$-norm. The capability of the method is shown by numerical examples for nonlinear conservation laws. Likewise, for the semi-discrete ALE-DG method for nonlinear Hamilton-Jacobi equations, error estimates are proven. In the one dimensional case the optimal $\left(k+1\right)$ convergence and in the two dimensional case the sub-optimal $\left(k+\frac{1}{2}\right)$ convergence are proven in the $\mathrm{L}^{2}$-norm, if a piecewise $\mathcal{P}^{k}$ polynomial approximation space is used on the reference cell. For the fullydiscrete method, the geometric conservation is proven and for the piecewise constant forward Euler step the convergence of the method to the unique physical relevant solution is discussed. N2 - Die vorliegende Arbeit beschäftigt sich mit der Entwicklung und Analyse von arbitrar Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) Methoden mit zeitabhängigen Testfunktionen Räumen für Erhaltungs- und Hamilton-Jacobi Gleichungen. Grundlagen über Erhaltungsgleichungen, Hamilton-Jacobi Gleichungen und discontinuous Galerkin Methoden werden präsentiert. Insbesondere werden Probleme bei der Entwicklung von discontinuous Galerkin Methoden für die Hamilton-Jacobi Gleichungen untersucht. Die Entwicklung der ALE-DG Methode basiert auf der Annahme, dass die Verteilung der Gitterpunkte zu einem kommenden Zeitpunkt explizit gegeben ist. Diese Annahme ermöglicht die Konstruktion einer zeitabhängigen lokal affin-linearen Abbildung auf eine Referenzzelle und eines zeitabhängigen Testfunktionen Raums. Zusätzlich kann eine Version des Reynolds’schen Transportsatzes gezeigt werden. Für die vollständig diskretisierte ALE-DG Methode für nichtlineare Erhaltungsgleichungen werden der geometrischen Erhaltungssatz und ein lokales Maximumprinzip bewiesen. Des Weiteren werden Bedingungen für Limiter angegeben. Diese Bedingungen sichern die Stabilität der Methode im Sinne der totalen Variation. Zusätzlich wird die Entropie-Stabilität der Methode diskutiert. Für die zugehörige semi-diskretisierte ALE-DG Methode werden Fehlerabschätzungen gezeigt. Wenn auf der Referenzzelle ein Testfunktionen Raum, der stückweise Polynome vom Grad $k$ enthält verwendet wird, kann für einen monotonen Fluss die suboptimale Konvergenzordnung $\left(k+\frac{1}{2}\right)$ und für einen upwind Fluss die optimale Konvergenzordnung $\left(k+1\right)$ in der $\mathrm{L}^{2}$-Norm gezeigt werden. Die Leistungsfähigkeit der Methode wird anhand numerischer Beispiele für nichtlineare Erhaltungsgleichungen untersucht. Ebenso werden für die semi-diskretisierte ALE-DG Methode für nichtlineare Hamilton-Jacobi Gleichungen Fehlerabschätzungen gezeigt. Wenn auf der Referenzzelle ein Testfunktionen Raum, der stückweise Polynome vom Grad k enthält verwendet wird, kann im eindimensionalen Fall die optimale Konvergenzordnung $\left(k+1\right)$ und im zweidimensionalen Fall die suboptimale Konvergenzordnung $\left(k+\frac{1}{2}\right)$ in der $\mathrm{L}^{2}$-Norm gezeigt werden. Für die vollständig diskretisierte ALE-DG Methode werden der geometrischen Erhaltungssatz bewiesen und für die stückweise konstante explizite Euler Diskretisierung wird die Konvergenz gegen die eindeutige physikalisch relevante Lösung diskutiert. KW - Galerkin-Methode KW - Numerische Strömungssimulation KW - Kontinuitätsgleichung KW - Hamilton-Jacobi-Differentialgleichung KW - Arbitrary Lagrangian-Eulerian KW - Discontinuous Galerkin method KW - Moving mesh method KW - conservation law Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-139579 N1 - zu dieser Arbeit gibt es einen Artikel, der in "Mathematics of Computation" veröffentlicht wurde unter folgendem Link: http://dx.doi.org/10.1090/mcom/3126 ER - TY - THES A1 - Warnecke, Sandra T1 - Numerical schemes for multi-species BGK equations based on a variational procedure applied to multi-species BGK equations with velocity-dependent collision frequency and to quantum multi-species BGK equations T1 - Numerische Verfahren für multispezies BGK Gleichungen mittels Variationsansatz angewandt auf multispezies BGK Gleichungen mit geschwindigkeitsabhängiger Stoßfrequenz sowie auf quantenmechanische multispezies BGK Gleichungen N2 - We consider a multi-species gas mixture described by a kinetic model. More precisely, we are interested in models with BGK interaction operators. Several extensions to the standard BGK model are studied. Firstly, we allow the collision frequency to vary not only in time and space but also with the microscopic velocity. In the standard BGK model, the dependence on the microscopic velocity is neglected for reasons of simplicity. We allow for a more physical description by reintroducing this dependence. But even though the structure of the equations remains the same, the so-called target functions in the relaxation term become more sophisticated being defined by a variational procedure. Secondly, we include quantum effects (for constant collision frequencies). This approach influences again the resulting target functions in the relaxation term depending on the respective type of quantum particles. In this thesis, we present a numerical method for simulating such models. We use implicit-explicit time discretizations in order to take care of the stiff relaxation part due to possibly large collision frequencies. The key new ingredient is an implicit solver which minimizes a certain potential function. This procedure mimics the theoretical derivation in the models. We prove that theoretical properties of the model are preserved at the discrete level such as conservation of mass, total momentum and total energy, positivity of distribution functions and a proper entropy behavior. We provide an array of numerical tests illustrating the numerical scheme as well as its usefulness and effectiveness. N2 - Wir betrachten ein Gasgemisch, das aus mehreren Spezies zusammengesetzt ist und durch kinetische Modelle beschrieben werden kann. Dabei interessieren wir uns vor allem für Modelle mit BGK-Wechselwirkungsoperatoren. Verschiedene Erweiterungen des Standard-BGK-Modells werden untersucht. Im ersten Modell nehmen wir eine Abhängigkeit der Stoßfrequenzen von der mikroskopischen Geschwindigkeit hinzu. Im Standard-BGK-Modell wird diese Abhängigkeit aus Gründen der Komplexität vernachlässigt. Wir nähern uns der physikalischen Realität weiter an, indem wir die Abhängigkeit von der mikroskopischen Geschwindigkeit beachten. Die Struktur der Gleichungen bleibt erhalten, allerdings hat dies Auswirkungen auf die sogenannten Zielfunktionen im Relaxationsterm, welche sodann durch einen Variationsansatz definiert werden. Das zweite Modell berücksichtigt Quanteneffekte (für konstante Stoßfrequenzen), was wiederum die Zielfunktionen im Relaxationsterm beeinflusst. Diese unterscheiden sich abhängig von den jeweils betrachteten, quantenmechanischen Teilchentypen. In dieser Doktorarbeit stellen wir numerische Verfahren vor, die auf oben beschriebene Modelle angewandt werden können. Wir legen eine implizite-explizite Zeitdiskretisierung zu Grunde, da die Relaxationsterme für große Stoßfrequenzen steif werden können. Das Kernstück ist ein impliziter Löser, der eine gewisse Potenzialfunktion minimiert. Dieses Vorgehen imitiert die theoretische Herleitung in den Modellen. Wir zeigen, dass die Eigenschaften des Modells auch auf der diskreten Ebene vorliegen. Dies beinhaltet die Massen-, Gesamtimpuls- und Gesamtenergieerhaltung, die Positivität von Verteilungsfunktionen sowie das gewünschte Verhalten der Entropie. Wir führen mehrere numerische Tests durch, die die Eigenschaften, die Nützlichkeit und die Zweckmäßigkeit des numerischen Verfahrens aufzeigen. N2 - Many applications require reliable numerical simulations of realistic set-ups e.g. plasma physics. This book gives a short introduction into kinetic models of gas mixtures describing the time evolution of rarefied gases and plasmas. Recently developed models are presented which extend existing literature by including more physical phenomena. We develop a numerical scheme for these more elaborated equations. The scheme is proven to maintain the physical properties of the models at the discrete level. We show several numerical test cases inspired by physical experiments. KW - Kinetische Gastheorie KW - Simulation KW - Numerisches Verfahren KW - Gasgemisch KW - Plasma KW - multi-fluid mixture KW - kinetic model KW - entropy minimization KW - IMEX KW - velocity-dependent collision frequency Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-282378 SN - 978-3-95826-192-1 SN - 978-3-95826-193-8 N1 - Parallel erschienen als Druckausgabe in Würzburg University Press, ISBN 978-3-95826-192-1, 32,90 EUR. PB - Würzburg University Press CY - Würzburg ET - 1. Auflage ER - 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 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-159965 ER - TY - THES A1 - Berberich, Jonas Philipp T1 - Fluids in Gravitational Fields – Well-Balanced Modifications for Astrophysical Finite-Volume Codes T1 - Fluide in Gravitationsfeldern - Wohl-Balancierte Modifikationen für Astrophysikalische Finite-Volumen-Codes N2 - Stellar structure can -- in good approximation -- be described as a hydrostatic state, which which arises due to a balance between gravitational force and pressure gradient. Hydrostatic states are static solutions of the full compressible Euler system with gravitational source term, which can be used to model the stellar interior. In order to carry out simulations of dynamical processes occurring in stars, it is vital for the numerical method to accurately maintain the hydrostatic state over a long time period. In this thesis we present different methods to modify astrophysical finite volume codes in order to make them \emph{well-balanced}, preventing them from introducing significant discretization errors close to hydrostatic states. Our well-balanced modifications are constructed so that they can meet the requirements for methods applied in the astrophysical context: They can well-balance arbitrary hydrostatic states with any equation of state that is applied to model thermodynamical relations and they are simple to implement in existing astrophysical finite volume codes. One of our well-balanced modifications follows given solutions exactly and can be applied on any grid geometry. The other methods we introduce, which do no require any a priori knowledge, balance local high order approximations of arbitrary hydrostatic states on a Cartesian grid. All of our modifications allow for high order accuracy of the method. The improved accuracy close to hydrostatic states is verified in various numerical experiments. N2 - Die Struktur von Sternen kann in guter Näherung als hydrostatischer Zustandbeschrieben werden, der durch ein Gleichgewicht zwischen Gravitationskraft undDruckgradient gegeben ist. Hydrostatische Zustände sind statische Lösungen dervollständigen komprimierbaren Euler-Gleichungen mit Gravitationsquellenterm, diezur Modellierung des Sterninneren verwendet werden können. Um Simulationendynamischer Prozesse in Sternen durchführen zu können, ist es wichtig, dass dieverwendete numerische Methode den hydrostatischen Zustand über einen langenZeitraum genau aufrechterhalten kann. In dieser Arbeit stellen wir verschiedene Me-thoden vor, um astrophysikalische Finite-Volumen-Codes so zu modifizieren, dasssie diewell-balancing-Eigenschaft erhalten, d.h., dass sie keine signifikanten Diskre-tisierungsfehler nahe hydrostatischer Zustände verursachen. Unsere well-balancing-Modifikationen sind so konstruiert, dass sie die Anforderungen für Methoden er-füllen, die im astrophysikalischen Kontext angewendet werden: Sie können beliebi-ge hydrostatische Zustände mit jeder Zustandsgleichung, die zur Modellierung derthermodynamischen Beziehungen angewendet wird, balancieren und sind einfach invorhandene astrophysikalische Finite-Volumen-Codes zu implementieren. Eine un-serer well-balancing Modifikationen erhält bekannte Lösungen exakt und kann aufjede Gittergeometrie angewendet werden. Die anderen Methoden, für die keine A-priori-Kenntnisse erforderlich sind, balancieren lokale Approximationen beliebigerhydrostatischer Zustände mit hoher Fehlerordnung auf einem kartesischen Gitter.Alle unsere Modifikationen erlauben eine hohe Fehlerordnung der Methode. Dieverbesserte Genauigkeit nahe an hydrostatischen Zuständen wird in verschiedenennumerischen Experimenten verifiziert. KW - well-balancing KW - Euler equations KW - finite volume methods KW - Fluid KW - Gravitationsfeld KW - Finite-Volumen-Methode Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-219679 ER - TY - THES A1 - Kanbar, Farah T1 - Asymptotic and Stationary Preserving Schemes for Kinetic and Hyperbolic Partial Differential Equations T1 - Asymptotische und Stationäre Erhaltungsverfahren für Kinetische und Hyperbolische Partielle Differentialgleichungen N2 - In this thesis, we are interested in numerically preserving stationary solutions of balance laws. We start by developing finite volume well-balanced schemes for the system of Euler equations and the system of MHD equations with gravitational source term. Since fluid models and kinetic models are related, this leads us to investigate AP schemes for kinetic equations and their ability to preserve stationary solutions. Kinetic models typically have a stiff term, thus AP schemes are needed to capture good solutions of the model. For such kinetic models, equilibrium solutions are reached after large time. Thus we need a new technique to numerically preserve stationary solutions for AP schemes. We find a criterion for SP schemes for kinetic equations which states, that AP schemes under a particular discretization are also SP. In an attempt to mimic our result for kinetic equations in the context of fluid models, for the isentropic Euler equations we developed an AP scheme in the limit of the Mach number going to zero. Our AP scheme is proven to have a SP property under the condition that the pressure is a function of the density and the latter is obtained as a solution of an elliptic equation. The properties of the schemes we developed and its criteria are validated numerically by various test cases from the literature. N2 - In dieser Arbeit interessieren wir uns für numerisch erhaltende stationäre Lösungen von Erhaltungsgleichungen. Wir beginnen mit der Entwicklung von well-balanced Finite-Volumen Verfahren für das System der Euler-Gleichungen und das System der MHD-Gleichungen mit Gravitationsquell term. Da Strömungsmodelle und kinetische Modelle miteinander verwandt sind, untersuchen wir asymptotisch erhaltende (AP) Verfahren für kinetische Gleichungen und ihre Fähigkeit, stationäre Lösungen zu erhalten. Kinetische Modelle haben typischerweise einen steifen Term, so dass AP Verfahren erforderlich sind, um gute Lösungen des Modells zu erhalten. Bei solchen kinetischen Modellen werden Gleichgewichtslösungen erst nach langer Zeit erreicht. Daher benötigen wir eine neue Technik, um stationäre Lösungen für AP Verfahren numerisch zu erhalten. Wir finden ein Kriterium für stationär-erhaltende (SP) Verfahren für kinetische Gleichungen, das besagt, dass AP Verfahren unter einer bestimmten Diskretisierung auch SP sind. In dem Versuch unser Ergebnis für kinetische Gleichungen im Kontext von Strömungsmodellen nachzuahmen, haben wir für die isentropen Euler-Gleichungen ein AP Verfahren für den Grenzwert der Mach-Zahl gegen Null, entwickelt. Unser AP Verfahren hat nachweislich eine SP Eigenschaft unter der Bedingung, dass der Druck eine Funktion der Dichte ist und letztere als Lösung einer elliptischen Gleichung erhalten wird. Die Eigenschaften des von uns entwickelten und seine Kriterien werden anhand verschiedener Testfälle aus der Literatur numerisch validiert. N2 - In this thesis, we are interested in numerically preserving stationary solutions of balance laws. We start by developing finite volume well-balanced schemes for the system of Euler equations and the system of Magnetohydrodynamics (MHD) equations with gravitational source term. Since fluid models and kinetic models are related, this leads us to investigate Asymptotic Preserving (AP) schemes for kinetic equations and their ability to preserve stationary solutions. In an attempt to mimic our result for kinetic equations in the context of fluid models, for the isentropic Euler equations we developed an AP scheme in the limit of the Mach number going to zero. The properties of the schemes we developed and its criteria are validated numerically by various test cases from the literature. KW - Angewandte Mathematik KW - Hyperbolische Differentialgleichung KW - Kinetische Gleichung KW - Euler-Lagrange-Gleichung KW - Magnetohydrodynamische Gleichung KW - Euler equations KW - isentropic Euler equations KW - MHD equations KW - kinetic equations KW - well-balanced scheme KW - asymptotic preserving KW - stationary preserving KW - hyperbolic partial differential equations Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-301903 SN - 978-3-95826-210-2 SN - 978-3-95826-211-9 N1 - Parallel erschienen als Druckausgabe in Würzburg University Press, ISBN 978-3-95826-210-2, 29,80 EUR. PB - Würzburg University Press CY - Würzburg ET - 1. Auflage ER - TY - THES A1 - Zenk, Markus T1 - On Numerical Methods for Astrophysical Applications T1 - Über numerische Methoden für astrophysikalische Anwendungen N2 - Diese Arbeit befasst sich mit der Approximation der Lösungen von Modellen zur Beschreibung des Strömungsverhaltens in Atmosphären. Im Speziellen umfassen die hier behandelten Modelle die kompressiblen Euler Gleichungen der Gasdynamik mit einem Quellterm bezüglich der Gravitation und die Flachwassergleichungen mit einem nicht konstanten Bodenprofil. Verschiedene Methoden wurden bereits entwickelt um die Lösungen dieser Gleichungen zu approximieren. Im Speziellen geht diese Arbeit auf die Approximation von Lösungen nahe des Gleichgewichts und, im Falle der Euler Gleichungen, bei kleinen Mach Zahlen ein. Die meisten numerischen Methoden haben die Eigenschaft, dass die Qualität der Approximation sich mit der Anzahl der Freiheitsgrade verbessert. In der Praxis werden deswegen diese numerischen Methoden auf großen Computern implementiert um eine möglichst hohe Approximationsgüte zu erreichen. Jedoch sind auch manchmal diese großen Maschinen nicht ausreichend, um die gewünschte Qualität zu erreichen. Das Hauptaugenmerk dieser Arbeit ist darauf gerichtet, die Qualität der Approximation bei gleicher Anzahl von Freiheitsgrade zu verbessern. Diese Arbeit ist im Zusammenhang einer Kollaboration zwischen Prof. Klingenberg des Mathemaitschen Instituts in Würzburg und Prof. Röpke des Astrophysikalischen Instituts in Würzburg entstanden. Das Ziel dieser Kollaboration ist es, Methoden zur Berechnung von stellarer Atmosphären zu entwickeln. In dieser Arbeit werden vor allem zwei Problemstellungen behandelt. Die erste Problemstellung bezieht sich auf die akkurate Approximation des Quellterms, was zu den so genannten well-balanced Schemata führt. Diese erlauben genaue Approximationen von Lösungen nahe des Gleichgewichts. Die zweite Problemstellung bezieht sich auf die Approximation von Strömungen bei kleinen Mach Zahlen. Es ist bekannt, dass Lösungen der kompressiblen Euler Gleichungen zu Lösungen der inkompressiblen Euler Gleichungen konvergieren, wenn die Mach Zahl gegen null geht. Klassische numerische Schemata zeigen ein stark diffusives Verhalten bei kleinen Mach Zahlen. Das hier entwickelte Schema fällt in die Kategorie der asymptotic preserving Schematas, d.h. das numerische Schema ist auf einem diskrete Level kompatibel mit dem auf dem Kontinuum gezeigten verhalten. Zusätzlich wird gezeigt, dass die Diffusion des hier entwickelten Schemas unabhängig von der Mach Zahl ist. In Kapitel 3 wird ein HLL approximativer Riemann Löser für die Approximation der Lösungen der Flachwassergleichungen mit einem nicht konstanten Bodenprofil angewendet und ein well-balanced Schema entwickelt. Die meisten well-balanced Schemata für die Flachwassergleichungen behandeln nur den Fall eines Fluids im Ruhezustand, die so genannten Lake at Rest Lösungen. Hier wird ein Schema entwickelt, welches sich mit allen Gleichgewichten befasst. Zudem wird eine zweiter Ordnung Methode entwickelt, welche im Gegensatz zu anderen in der Literatur nicht auf einem iterativen Verfahren basiert. Numerische Experimente werden durchgeführt um die Vorteile des neuen Verfahrens zu zeigen. In Kapitel 4 wird ein Suliciu Relaxations Löser angepasst um die hydrostatischen Gleichgewichte der Euler Gleichungen mit einem Gravitationspotential aufzulösen. Die Gleichungen der hydrostatischen Gleichgewichte sind unterbestimmt und lassen deshalb keine Eindeutigen Lösungen zu. Es wird jedoch gezeigt, dass das neue Schema für eine große Klasse dieser Lösungen die well-balanced Eigenschaft besitzt. Für bestimmte Klassen werden Quadraturformeln zur Approximation des Quellterms entwickelt. Es wird auch gezeigt, dass das Schema robust, d.h. es erhält die Positivität der Masse und Energie, und stabil bezüglich der Entropieungleichung ist. Die numerischen Experimente konzentrieren sich vor allem auf den Einfluss der Quadraturformeln auf die well-balanced Eigenschaften. In Kapitel 5 wird ein Suliciu Relaxations Schema angepasst für Simulationen im Bereich kleiner Mach Zahlen. Es wird gezeigt, dass das neue Schema asymptotic preserving und die Diffusion kontrolliert ist. Zudem wird gezeigt, dass das Schema für bestimmte Parameter robust ist. Eine Stabilität wird aus einer Chapman-Enskog Analyse abgeleitet. Resultate numerische Experimente werden gezeigt um die Vorteile des neuen Verfahrens zu zeigen. In Kapitel 6 werden die Schemata aus den Kapiteln 4 und 5 kombiniert um das Verhalten des numerischen Schemas bei Flüssen mit kleiner Mach Zahl in durch die Gravitation geschichteten Atmosphären zu untersuchen. Es wird gezeigt, dass das Schema well-balanced ist. Die Robustheit und die Stabilität werden analog zu Kapitel 5 behandelt. Auch hier werden numerische Tests durchgeführt. Es zeigt sich, dass das neu entwickelte Schema in der Lage ist, die Dynamiken besser Aufzulösen als vor der Anpassung. Das Kapitel 7 beschäftigt sich mit der Entwicklung eines multidimensionalen Schemas basierend auf der Suliciu Relaxation. Jedoch ist die Arbeit an diesem Ansatz noch nicht beendet und numerische Resultate können nicht präsentiert werden. Es wird aufgezeigt, wo sich die Schwächen dieses Ansatzes befinden und weiterer Entwicklungsbedarf besteht. N2 - This work is concerned with the numerical approximation of solutions to models that are used to describe atmospheric or oceanographic flows. In particular, this work concen- trates on the approximation of the Shallow Water equations with bottom topography and the compressible Euler equations with a gravitational potential. Numerous methods have been developed to approximate solutions of these models. Of specific interest here are the approximations of near equilibrium solutions and, in the case of the Euler equations, the low Mach number flow regime. It is inherent in most of the numerical methods that the quality of the approximation increases with the number of degrees of freedom that are used. Therefore, these schemes are often run in parallel on big computers to achieve the best pos- sible approximation. However, even on those big machines, the desired accuracy can not be achieved by the given maximal number of degrees of freedom that these machines allow. The main focus in this work therefore lies in the development of numerical schemes that give better resolution of the resulting dynamics on the same number of degrees of freedom, compared to classical schemes. This work is the result of a cooperation of Prof. Klingenberg of the Institute of Mathe- matics in Wu¨rzburg and Prof. R¨opke of the Astrophysical Institute in Wu¨rzburg. The aim of this collaboration is the development of methods to compute stellar atmospheres. Two main challenges are tackled in this work. First, the accurate treatment of source terms in the numerical scheme. This leads to the so called well-balanced schemes. They allow for an accurate approximation of near equilibrium dynamics. The second challenge is the approx- imation of flows in the low Mach number regime. It is known that the compressible Euler equations tend towards the incompressible Euler equations when the Mach number tends to zero. Classical schemes often show excessive diffusion in that flow regime. The here devel- oped scheme falls into the category of an asymptotic preserving scheme, i.e. the numerical scheme reflects the behavior that is computed on the continuous equations. Moreover, it is shown that the diffusion of the numerical scheme is independent of the Mach number. In chapter 3, an HLL-type approximate Riemann solver is adapted for simulations of the Shallow Water equations with bottom topography to develop a well-balanced scheme. In the literature, most schemes only tackle the equilibria when the fluid is at rest, the so called Lake at rest solutions. Here a scheme is developed to accurately capture all the equilibria of the Shallow Water equations. Moreover, in contrast to other works, a second order extension is proposed, that does not rely on an iterative scheme inside the reconstruction procedure, leading to a more efficient scheme. In chapter 4, a Suliciu relaxation scheme is adapted for the resolution of hydrostatic equilibria of the Euler equations with a gravitational potential. The hydrostatic relations are underdetermined and therefore the solutions to that equations are not unique. However, the scheme is shown to be well-balanced for a wide class of hydrostatic equilibria. For specific classes, some quadrature rules are computed to ensure the exact well-balanced property. Moreover, the scheme is shown to be robust, i.e. it preserves the positivity of mass and energy, and stable with respect to the entropy. Numerical results are presented in order to investigate the impact of the different quadrature rules on the well-balanced property. In chapter 5, a Suliciu relaxation scheme is adapted for the simulations of low Mach number flows. The scheme is shown to be asymptotic preserving and not suffering from excessive diffusion in the low Mach number regime. Moreover, it is shown to be robust under certain parameter combinations and to be stable from an Chapman-Enskog analysis. Numerical results are presented in order to show the advantages of the new approach. In chapter 6, the schemes developed in the chapters 4 and 5 are combined in order to investigate the performance of the numerical scheme in the low Mach number regime in a gravitational stratified atmosphere. The scheme is shown the be well-balanced, robust and stable with respect to a Chapman-Enskog analysis. Numerical tests are presented to show the advantage of the newly proposed method over the classical scheme. In chapter 7, some remarks on an alternative way to tackle multidimensional simulations are presented. However no numerical simulations are performed and it is shown why further research on the suggested approach is necessary. KW - Strömung KW - Numerical Methods KW - Hyperbolic Partial Differential Equations KW - Well-Balanced KW - Asymptotic Preserving KW - Atmosphäre KW - Mathematisches Modell KW - PDE Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-162669 ER - TY - THES A1 - Pirner, Marlies T1 - Kinetic modelling of gas mixtures T1 - Kinetische Modellierung von Gasgemischen N2 - This book deals with the kinetic modelling of gas mixtures. It extends the existing literature in mathematics for one species of gas to the case of gasmixtures. This is more realistic in applications. Thepresentedmodel for gas mixtures is proven to be consistentmeaning it satisfies theconservation laws, it admitsanentropy and an equilibriumstate. Furthermore, we can guarantee the existence, uniqueness and positivity of solutions. Moreover, the model is used for different applications, for example inplasma physics, for fluids with a small deviation from equilibrium and in the case of polyatomic gases. N2 - Die vorliegende Arbeit beschäftigt sich mit der Modellierung von Gasgemischen mittels einer kinetischen Beschreibung. Es werden Grundlagen über die Boltzmanngleichung für Gasgemische und die BGK-Aproximation präsentiert. Insbesondere wird auf deren Erweiterung auf Gasgemische eingegangen. Es wird ein Gasgemisch bestehend aus zwei Sorten von Gasen ohne chemische Reaktionen betrachtet. Das Gemisch wird mittels eines Systems kinetischer BGK-Gleichungen modelliert, welches je zwei Wechselwirkungsterme enthält, die den Impuls- und Energieaustausch berücksichtigen. Das hier vorgestellte Modell enthält einige von Physikern und Ingenieuren vorgeschlagene Modelle als Spezialfälle. Es wird gezeigt, dass das hier vorgeschlagene Modell die wesentlichen physikalischen Eigenschaften, wie Erhaltungseigenschaften, Positivität aller Temperaturen, das H-Theorem und Maxwellverteilungen im Gleichgewicht, erfüllt. Des Weiteren können die üblichen makroskopischen Gleichungen daraus hergeleitet werden. In der Literatur gibt es ein weiteres vorgeschlagenes Modell für Gasgemische mit nur einem Wechselwirkungsterm von Andries, Aoki und Perthame. In dieser Arbeit werden die Vorteile dieses Modells aus der Literatur und des hier vorgeschlagenen Modells diskutiert. Es wird die Nützlichkeit des hier vorgeschlagenen Modells illustriert, indem es dazu benutzt wird eine unbekannte Funktion in dem makroskopischen Modell für Gasgemische von Dellacherie herzuleiten. Des Weiteren wird für jedes dieser beiden Modelle Existenz, Eindeutigkeit und Positivität der Lösungen gezeigt. Dann wird das hier vorgeschlagene Modell auf bestimmte physikalische Situationen angewandt: auf Elektronen und Ionen in einem Plasma, auf ein Gasgemisch, welches sich nicht im Gleichgewicht befindet und ein Gasgemisch bestehend aus Molekülen mit zusätzlichen inneren Freiheitsgraden. Als erste Anwendung wird das Modell für geladene Teilchen erweitert und auf ein Gemisch aus Elektronen und Ionen angewandt, welches sich teilweise im Gleichgewicht befindet, teilweise nicht. Man findet solch eine Konstellation zum Beispiel bei der Fusion in einem Tokamak. Das Modell, welches hier vorgestellt wird, wird hier benutzt, da es die Wechselwirkungen zwischen Teilchen von der gleichen Sorte und Wechselwirkungen zwischen Teilchen verschiedener Sorten separiert. Dann wird ein neues Modell mithilfe der Mikro-Makro-Zerlegung hergeleitet, welches numerisch in einem Regime angewandt wird, in dem Gase teilweise im Gleichgewicht sind, teilweise nicht. Es werden theoretische Ergebnisse vorgestellt, zum einen Konvergenzraten gegen das Gleichgewicht im räumlich homogenen Fall, zum anderen die Landau-Dämpfung für Gasgemische, um sie mit Ergebnissen aus numerischen Simulationen vergleichen zu können. Als zweite Anwendung wird ein Gasgemisch betrachtet, welches eine Abweichung vom Gleichgewichtszustand hat und makroskopisch mithilfe der Navier-Stokes-Gleichungen beschrieben wird. In dieser makroskopischen Beschreibung erwartet man vier physikalische Größen, die das physikalische Verhalten eines Gases beschreiben, den Diffusionskoeffizienten, den Viskositätskoeffizienten, die Wärmeleitfähigkeit und den thermischen Diffusionsparameter. Es wird eine Chapman-Enskog-Entwicklung des hier vorgestellten Modells durchgeführt, um drei dieser vier physikalischen Größen zu bestimmen. Zusatzlich werden mehrere mögliche Erweiterungen zu einem ES-BGK-Modell für Gasgemische vorgeschlagen um die vierte physikalische Größe zu bestimmen. Es wird eine Erweiterung präsentiert, die möglichst einfach gehalten ist, eine intuitive Erweiterung, die den Fall einer Gassorte ähnelt und eine Erweiterung, die die physikalische Motivation des Physikers Holway, der das ES-BGK-Modell erfunden hat, berücksichtigt. Es wird gezeigt, dass die Erweiterungen die Erhaltungseigenschaften erfüllen, alle Temperaturen positiv sind und das H-Theorem erfüllt ist. Als dritte Anwendung wird das hier vorgestellte Modell zu einem Modell für Moleküle mit zusätzlichen inneren Freiheitsgraden erweitert. Die zwei Gassorten dürfen dabei eine unterschiedliche Anzahl an inneren Freiheitsgraden haben und werden beschrieben durch ein System von kinetischen ES-BGK-Gleichungen. Es wird gezeigt, dass das Modell die Erhaltungseigenschaften erfülllt, dass alle Temperaturen positiv sind und dass das H-Theorem erfüllt ist. Für numerische Zwecke wird die Chu-Reduktion angewandt um die Komplexität des Modells zu reduzieren und eine Anwendung gezeigt, bei dem eine Gassorte keine inneren Freiheitsgrade hat und die andere Sorte zwei Rotationsfreiheitsgrade besitzt. Als letztes wird der Grenzwert des hier vorgestellten Modells zu den dissipativen Eulergleichungen bewiesen. N2 - The present thesis considers the modelling of gas mixtures via a kinetic description. Fundamentals about the Boltzmann equation for gas mixtures and the BGK approximation are presented. Especially, issues in extending these models to gas mixtures are discussed. A non-reactive two component gas mixture is considered. The two species mixture is modelled by a system of kinetic BGK equations featuring two interaction terms to account for momentum and energy transfer between the two species. The model presented here contains several models from physicists and engineers as special cases. Consistency of this model is proven: conservation properties, positivity of all temperatures and the H-theorem. The form in global equilibrium as Maxwell distributions is specified. Moreover, the usual macroscopic conservation laws can be derived. In the literature, there is another type of BGK model for gas mixtures developed by Andries, Aoki and Perthame, which contains only one interaction term. In this thesis, the advantages of these two types of models are discussed and the usefulness of the model presented here is shown by using this model to determine an unknown function in the energy exchange of the macroscopic equations for gas mixtures described in the literature by Dellacherie. In addition, for each of the two models existence and uniqueness of mild solutions is shown. Moreover, positivity of classical solutions is proven. Then, the model presented here is applied to three physical applications: a plasma consisting of ions and electrons, a gas mixture which deviates from equilibrium and a gas mixture consisting of polyatomic molecules. First, the model is extended to a model for charged particles. Then, the equations of magnetohydrodynamics are derived from this model. Next, we want to apply this extended model to a mixture of ions and electrons in a special physical constellation which can be found for example in a Tokamak. The mixture is partly in equilibrium in some regions, in some regions it deviates from equilibrium. The model presented in this thesis is taken for this purpose, since it has the advantage to separate the intra and interspecies interactions. Then, a new model based on a micro-macro decomposition is proposed in order to capture the physical regime of being partly in equilibrium, partly not. Theoretical results are presented, convergence rates to equilibrium in the space-homogeneous case and the Landau damping for mixtures, in order to compare it with numerical results. Second, the model presented here is applied to a gas mixture which deviates from equilibrium such that it is described by Navier-Stokes equations on the macroscopic level. In this macroscopic description it is expected that four physical coefficients will show up, characterizing the physical behaviour of the gases, namely the diffusion coefficient, the viscosity coefficient, the heat conductivity and the thermal diffusion parameter. A Chapman-Enskog expansion of the model presented here is performed in order to capture three of these four physical coefficients. In addition, several possible extensions to an ellipsoidal statistical model for gas mixtures are proposed in order to capture the fourth coefficient. Three extensions are proposed: An extension which is as simple as possible, an intuitive extension copying the one species case and an extension which takes into account the physical motivation of the physicist Holway who invented the ellipsoidal statistical model for one species. Consistency of the extended models like conservation properties, positivity of all temperatures and the H-theorem are proven. The shape of global Maxwell distributions in equilibrium are specified. Third, the model presented here is applied to polyatomic molecules. A multi component gas mixture with translational and internal energy degrees of freedom is considered. The two species are allowed to have different degrees of freedom in internal energy and are modelled by a system of kinetic ellipsoidal statistical equations. Consistency of this model is shown: conservation properties, positivity of the temperature, H-theorem and the form of Maxwell distributions in equilibrium. For numerical purposes the Chu reduction is applied to the developed model for polyatomic gases to reduce the complexity of the model and an application for a gas consisting of a mono-atomic and a diatomic gas is given. Last, the limit from the model presented here to the dissipative Euler equations for gas mixtures is proven. KW - Polyatomare Verbindungen KW - Gasgemisch KW - Transportkoeffizient KW - Plasma KW - multi-fluid mixture KW - kinetic description of gases KW - entropy inequality KW - transport coefficients KW - Modellierung KW - well posedness KW - plasma modelling KW - polyatomic molecules KW - hydrodynamic limits Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-161077 SN - 978-3-95826-080-1 (Print) SN - 978-3-95826-081-8 (Online) N1 - Parallel erschienen als Druckausgabe in Würzburg University Press, ISBN 978-3-95826-080-1, 27,80 EUR. PB - Würzburg University Press CY - Würzburg ET - 1. Auflage ER - TY - THES A1 - Herrmann, Marc T1 - The Total Variation on Surfaces and of Surfaces T1 - Die totale Variation auf Oberflächen und von Oberflächen N2 - This thesis is concerned with applying the total variation (TV) regularizer to surfaces and different types of shape optimization problems. The resulting problems are challenging since they suffer from the non-differentiability of the TV-seminorm, but unlike most other priors it favors piecewise constant solutions, which results in piecewise flat geometries for shape optimization problems.The first part of this thesis deals with an analogue of the TV image reconstruction approach [Rudin, Osher, Fatemi (Physica D, 1992)] for images on smooth surfaces. A rigorous analytical framework is developed for this model and its Fenchel predual, which is a quadratic optimization problem with pointwise inequality constraints on the surface. A function space interior point method is proposed to solve it. Afterwards, a discrete variant (DTV) based on a nodal quadrature formula is defined for piecewise polynomial, globally discontinuous and continuous finite element functions on triangulated surface meshes. DTV has favorable properties, which include a convenient dual representation. Next, an analogue of the total variation prior for the normal vector field along the boundary of smooth shapes in 3D is introduced. Its analysis is based on a differential geometric setting in which the unit normal vector is viewed as an element of the two-dimensional sphere manifold. Shape calculus is used to characterize the relevant derivatives and an variant of the split Bregman method for manifold valued functions is proposed. This is followed by an extension of the total variation prior for the normal vector field for piecewise flat surfaces and the previous variant of split Bregman method is adapted. Numerical experiments confirm that the new prior favours polyhedral shapes. N2 - Die vorliegende Arbeit beschäftigt sich mit der Anwendung der totalen Variation (TV) als Regularisierung auf Oberflächen und in verschiedenen Problemen der Formoptimierung. Die daraus entstehenden Optimierungsprobleme sind aufgrund der TV-Seminorm nicht differenzierbar und daher eine Herausforderung. Allerdings werden dadurch, im Gegensatz zu anderen Regularisierungen, stückweise konstante Lösungen favorisiert. Dies führt bei Problemen der Formoptimierung zu stückweise flachen Geometrien. Der erste Teil dieser Arbeit widmet sich der Erweiterung des Ansatzes zur mathematischen Bildverarbeitung [Rudin, Osher, Fatemi (Physica D, 1992)] von flachen Bildern auf glatte Oberflächen und deren Texturen. Für das damit verbundene Optimierungsproblem wird das Fenchel präduale Problem hergeleitet. Dies ist ein quadratisches Optimierungsproblem mit Ungleichungsrestriktionen für dessen Lösung ein Innere-Punkte-Verfahren in Funktionenräumen vorgestellt wird. Basierend auf einer Quadraturformel, wird im Anschluss eine diskrete Variante (DTV) der TV-Seminorm für global unstetige und stetige Finite- Elemente-Funktionen auf triangulierten Oberflächen entwickelt. (DTV) besitzt positive Eigenschaften, wie eine praktische duale Darstellung. Im letzten Teil wird zuerst ein TV-Analogon für die Oberflächennormale von glatten Formen in 3D gezeigt und mit Hilfe von Differentialgeometrie analysiert. Danach wird eine mögliche Erweiterungen für stückweise glatte Oberflächen vorgestellt. Zur Lösung von beiden Regularisierungen wird eine Variante des Split-Bregman-Verfahrens für Funktionen mitWerten auf Mannigfaltigkeiten benutzt. KW - Gestaltoptimierung KW - optimization KW - total variation KW - Formoptimierung KW - Shape Optimization KW - Optimierung KW - Totale Variation KW - Finite-Elemente-Methode Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-240736 ER - TY - THES A1 - Birke, Claudius B. T1 - Low Mach and Well-Balanced Numerical Methods for Compressible Euler and Ideal MHD Equations with Gravity T1 - Low Mach und Well-Balanced Numerische Verfahren für die kompressiblen Euler und idealen MHD Gleichungen mit Gravitation N2 - Physical regimes characterized by low Mach numbers and steep stratifications pose severe challenges to standard finite volume methods. We present three new methods specifically designed to navigate these challenges by being both low Mach compliant and well-balanced. These properties are crucial for numerical methods to efficiently and accurately compute solutions in the regimes considered. First, we concentrate on the construction of an approximate Riemann solver within Godunov-type finite volume methods. A new relaxation system gives rise to a two-speed relaxation solver for the Euler equations with gravity. Derived from fundamental mathematical principles, this solver reduces the artificial dissipation in the subsonic regime and preserves hydrostatic equilibria. The solver is particularly stable as it satisfies a discrete entropy inequality, preserves positivity of density and internal energy, and suppresses checkerboard modes. The second scheme is designed to solve the equations of ideal MHD and combines different approaches. In order to deal with low Mach numbers, it makes use of a low-dissipation version of the HLLD solver and a partially implicit time discretization to relax the CFL time step constraint. A Deviation Well-Balancing method is employed to preserve a priori known magnetohydrostatic equilibria and thereby reduces the magnitude of spatial discretization errors in strongly stratified setups. The third scheme relies on an IMEX approach based on a splitting of the MHD equations. The slow scale part of the system is discretized by a time-explicit Godunov-type method, whereas the fast scale part is discretized implicitly by central finite differences. Numerical dissipation terms and CFL time step restriction of the method depend solely on the slow waves of the explicit part, making the method particularly suited for subsonic regimes. Deviation Well-Balancing ensures the preservation of a priori known magnetohydrostatic equilibria. The three schemes are applied to various numerical experiments for the compressible Euler and ideal MHD equations, demonstrating their ability to accurately simulate flows in regimes with low Mach numbers and strong stratification even on coarse grids. N2 - Physikalische Regime mit sehr niedrigen Machzahlen und starken Abschichtungen stellen konventionelle Finite Volumen Verfahren vor erhebliche Herausforderungen. In dieser Arbeit präsentieren wir drei neue Verfahren, die in der Lage sind, die Herausforderungen zu bewältigen. Die neuen Verfahren sind speziell an kleine Machzahlen angepasst und können (magneto-)hydrostatische Gleichgewichte exakt erhalten. Diese Eigenschaften sind essentiell für eine effiziente Berechnung präziser Lösungen in den betrachteten Regimen. Zunächst konzentrieren wir uns auf die Konstruktion eines approximativen Riemannlösers innerhalb von Godunov-artigen Finite Volumen Verfahren. Ein neues Relaxationssystem führt zu einem Relaxationslöser für die Euler Gleichungen mit Gravitation, der zwei Relaxationsgeschwindigkeiten verwendet. Abgeleitet von grundlegenden mathematischen Prinzipien reduziert dieser Löser die künstliche Dissipation im subsonischen Bereich und erhält hydrostatische Gleichgewichte. Der Löser ist besonders stabil, da er eine diskrete Entropieungleichung erfüllt, die Positivität von Dichte und interner Energie bewahrt und Schachbrettmuster unterdrückt. Das zweite Verfahren löst die idealen MHD Gleichungen und kombiniert verschiedene Ansätze, um die einzelnen numerischen Herausforderungen zu bewältigen. Für einen effizienten Umgang mit niedrigen Machzahlen wird eine Variante des HLLD Lösers mit künstlich niedriger Dissipation sowie eine teilweise implizite Zeitdiskretisierung zur Lockerung der CFL Zeitschrittbeschränkung gewählt. Eine Deviation Well-Balancing Methode wird angewendet, um magnetohydrostatische Gleichgewichte zu bewahren und dadurch das Ausmaß von räumlichen Diskretisierungsfehlern in stark geschichteten Atmosphären zu reduzieren. Das dritte Verfahren verwendet einen IMEX Ansatz, welcher auf einer Aufspaltung der MHD Gleichungen basiert. Das Teilsystem mit langsamen Ausbreitungsgeschwindigkeiten wird durch eine zeit-explizite Godunov-artige Methode diskretisiert, während das Teilsystem mit schnellen Ausbreitungsgeschwindigkiten implizit durch zentrale finite Differenzen diskretisiert wird. Numerische Dissipationsterme und die CFL Zeitschrittbeschränkung der Methode hängen somit nur von den langsamen Wellen des expliziten Teils ab, so dass die Methode besonders für subsonische Regime geeignet ist. Deviation Well-Balancing gewährleistet die Erhaltung a priori bekannter magnetohydrostatischer Gleichgewichte. Die drei Verfahren werden auf numerische Experimente für die kompressiblen Euler und idealen MHD Gleichungen angewendet und zeigen darin ihre Fähigkeit, Strömungen in Regimen mit niedrigen Machzahlen und starker Schichtung auch auf groben diskreten Gittern akkurat zu simulieren. KW - Magnetohydrodynamik KW - Numerische Strömungssimulation KW - Finite-Volumen-Methode KW - relaxation method KW - IMEX scheme KW - low Mach number KW - well-balanced Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-363303 ER -