@phdthesis{Klotzky2018,
  author    = {Klotzky, Jens},
  title     = {Well-posedness of a fluid-particle interaction model},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-169009},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2018},
  abstract  = {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{\`e}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.},
  subject      = {Hyperbolische Differentialgleichung},
  language  = {en}
}
@phdthesis{GallegoValencia2017,
  author    = {Gallego Valencia, Juan Pablo},
  title     = {On Runge-Kutta discontinuous Galerkin methods for compressible Euler equations and the ideal magneto-hydrodynamical model},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-148874},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2017},
  abstract  = {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.},
  subject      = {Eulersche Differentialgleichung},
  language  = {en}
}
@phdthesis{Schnuecke2016,
  author    = {Schn{\"u}cke, Gero},
  title     = {Arbitrary Lagrangian-Eulerian Discontinous Galerkin methods for nonlinear time-dependent first order partial differential equations},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-139579},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2016},
  abstract  = {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.},
  subject      = {Galerkin-Methode},
  language  = {en}
}
@phdthesis{Kanbar2023,
  author    = {Kanbar, Farah},
  title     = {Asymptotic and Stationary Preserving Schemes for Kinetic and Hyperbolic Partial Differential Equations},
  edition   = {1. Auflage},
  publisher = {W{\"u}rzburg University Press},
  address   = {W{\"u}rzburg},
  isbn      = {978-3-95826-210-2},
  doi       = {10.25972/WUP-978-3-95826-211-9},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-301903},
  school      = {W{\"u}rzburg University Press},
  pages     = {xiv, 137},
  year      = {2023},
  abstract  = {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.},
  subject      = {Angewandte Mathematik},
  language  = {en}
}
@phdthesis{Zenk2018,
  author    = {Zenk, Markus},
  title     = {On Numerical Methods for Astrophysical Applications},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-162669},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2018},
  abstract  = {Diese Arbeit befasst sich mit der Approximation der L{\"o}sungen von Modellen zur Beschreibung des Str{\"o}mungsverhaltens in Atmosph{\"a}ren. Im Speziellen umfassen die hier behandelten Modelle die kompressiblen Euler Gleichungen der Gasdynamik mit einem Quellterm bez{\"u}glich der Gravitation und die Flachwassergleichungen mit einem nicht konstanten Bodenprofil. Verschiedene Methoden wurden bereits entwickelt um die L{\"o}sungen dieser Gleichungen zu approximieren. Im Speziellen geht diese Arbeit auf die Approximation von L{\"o}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{\"a}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{\"o}glichst hohe Approximationsg{\"u}te zu erreichen. Jedoch sind auch manchmal diese großen Maschinen nicht ausreichend, um die gew{\"u}nschte Qualit{\"a}t zu erreichen. Das Hauptaugenmerk dieser Arbeit ist darauf gerichtet, die Qualit{\"a}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{\"u}rzburg und Prof. R{\"o}pke des Astrophysikalischen Instituts in W{\"u}rzburg entstanden. Das Ziel dieser Kollaboration ist es, Methoden zur Berechnung von stellarer Atmosph{\"a}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{\"u}hrt. Diese erlauben genaue Approximationen von L{\"o}sungen nahe des Gleichgewichts. Die zweite Problemstellung bezieht sich auf die Approximation von Str{\"o}mungen bei kleinen Mach Zahlen. Es ist bekannt, dass L{\"o}sungen der kompressiblen Euler Gleichungen zu L{\"o}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{\"a}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{\"a}tzlich wird gezeigt, dass die Diffusion des hier entwickelten Schemas unabh{\"a}ngig von der Mach Zahl ist. In Kapitel 3 wird ein HLL approximativer Riemann L{\"o}ser f{\"u}r die Approximation der L{\"o}sungen der Flachwassergleichungen mit einem nicht konstanten Bodenprofil angewendet und ein well-balanced Schema entwickelt. Die meisten well-balanced Schemata f{\"u}r die Flachwassergleichungen behandeln nur den Fall eines Fluids im Ruhezustand, die so genannten Lake at Rest L{\"o}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{\"u}hrt um die Vorteile des neuen Verfahrens zu zeigen. In Kapitel 4 wird ein Suliciu Relaxations L{\"o}ser angepasst um die hydrostatischen Gleichgewichte der Euler Gleichungen mit einem Gravitationspotential aufzul{\"o}sen. Die Gleichungen der hydrostatischen Gleichgewichte sind unterbestimmt und lassen deshalb keine Eindeutigen L{\"o}sungen zu. Es wird jedoch gezeigt, dass das neue Schema f{\"u}r eine große Klasse dieser L{\"o}sungen die well-balanced Eigenschaft besitzt. F{\"u}r bestimmte Klassen werden Quadraturformeln zur Approximation des Quellterms entwickelt. Es wird auch gezeigt, dass das Schema robust, d.h. es erh{\"a}lt die Positivit{\"a}t der Masse und Energie, und stabil bez{\"u}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{\"u}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{\"u}r bestimmte Parameter robust ist. Eine Stabilit{\"a}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{\"u}ssen mit kleiner Mach Zahl in durch die Gravitation geschichteten Atmosph{\"a}ren zu untersuchen. Es wird gezeigt, dass das Schema well-balanced ist. Die Robustheit und die Stabilit{\"a}t werden analog zu Kapitel 5 behandelt. Auch hier werden numerische Tests durchgef{\"u}hrt. Es zeigt sich, dass das neu entwickelte Schema in der Lage ist, die Dynamiken besser Aufzul{\"o}sen als vor der Anpassung. Das Kapitel 7 besch{\"a}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{\"o}nnen nicht pr{\"a}sentiert werden. Es wird aufgezeigt, wo sich die Schw{\"a}chen dieses Ansatzes befinden und weiterer Entwicklungsbedarf besteht.},
  subject      = {Str{\"o}mung},
  language  = {en}
}
@phdthesis{Pirner2018,
  author    = {Pirner, Marlies},
  title     = {Kinetic modelling of gas mixtures},
  edition   = {1. Auflage},
  publisher = {W{\"u}rzburg University Press},
  address   = {W{\"u}rzburg},
  isbn      = {978-3-95826-080-1 (Print)},
  doi       = {10.25972/WUP-978-3-95826-081-8},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-161077},
  school      = {W{\"u}rzburg University Press},
  pages     = {xi, 222},
  year      = {2018},
  abstract  = {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.},
  subject      = {Polyatomare Verbindungen},
  language  = {en}
}
@phdthesis{Warnecke2022,
  author    = {Warnecke, Sandra},
  title     = {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},
  edition   = {1. Auflage},
  publisher = {W{\"u}rzburg University Press},
  address   = {W{\"u}rzburg},
  isbn      = {978-3-95826-192-1},
  doi       = {10.25972/WUP-978-3-95826-193-8},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-282378},
  school      = {W{\"u}rzburg University Press},
  pages     = {xiii, 203},
  year      = {2022},
  abstract  = {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.},
  subject      = {Kinetische Gastheorie},
  language  = {en}
}
@phdthesis{Barsukow2018,
  author    = {Barsukow, Wasilij},
  title     = {Low Mach number finite volume methods for the acoustic and Euler equations},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-159965},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2018},
  abstract  = {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.},
  subject      = {Finite-Volumen-Methode},
  language  = {en}
}
@phdthesis{Berberich2021,
  author    = {Berberich, Jonas Philipp},
  title     = {Fluids in Gravitational Fields - Well-Balanced Modifications for Astrophysical Finite-Volume Codes},
  doi       = {10.25972/OPUS-21967},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-219679},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2021},
  abstract  = {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.},
  subject      = {Fluid},
  language  = {en}
}
@phdthesis{Herrmann2021,
  author    = {Herrmann, Marc},
  title     = {The Total Variation on Surfaces and of Surfaces},
  doi       = {10.25972/OPUS-24073},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-240736},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2021},
  abstract  = {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.},
  subject      = {Gestaltoptimierung},
  language  = {en}
}
@phdthesis{Birke2024,
  author    = {Birke, Claudius B.},
  title     = {Low Mach and Well-Balanced Numerical Methods for Compressible Euler and Ideal MHD Equations with Gravity},
  doi       = {10.25972/OPUS-36330},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-363303},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2024},
  abstract  = {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.},
  subject      = {Magnetohydrodynamik},
  language  = {en}
}