Discontinuous Galerkin Spectral Element Methods for Astrophysical Flows in Multi-physics Applications

In engineering applications, discontinuous Galerkin methods (DG) have been proven to be a powerful and flexible class of high order methods for problems in computational fluid dynamics. However, the potential benefits of DG for applications in astrophysical contexts is still relatively unexplored in its entirety. To this day, a decent number of studies surveying DG for astrophysical flows have been conducted. But the adoption of DG by the astrophysics community is just beginning to gain traction and integration of DG into established, multi-physics simulation frameworks for comprehensive astrophysical modeling is still lacking. It is our firm believe, that the full potential of novel approaches for numerically solving the fluid equations only shows under the pressure of real-world simulations with all aspects of multi-physics, challenging flow configurations, resolution and runtime constraints, and efficiency metrics on high-performance systems involved. Thus, we see the pressing need to propel DG from the well-trodden path of cataloguing test results under "optimal" laboratory conditions towards the harsh and unforgiving environment of large-scale astrophysics simulations. Consequently, the core of this work is the development and deployment of a robust DG scheme solving the ideal magneto-hydrodynamics equations with multiple species on three-dimensional Cartesian grids with adaptive mesh refinement. We chose to implement DG within the venerable simulation framework FLASH, with a specific focus on multi-physics problems in astrophysics. This entails modifications of the vanilla DG scheme to make it fit seamlessly within FLASH in such a way that all other physics modules can be naturally coupled without additional implementation overhead. A key ingredient is that our DG scheme uses mean value data organized into blocks - the central data structure in FLASH. Having the opportunity to work on mean values, allows us to rely on a rock-solid, monotone Finite Volume (FV) scheme as "backup" whenever the high order DG method fails in cases when the flow gets too harsh. Finding ways to combine the two schemes in a fail-safe manner without loosing primary conservation while still maintaining high order accuracy for smooth, well-resolved flows involves a series of careful considerations, which we document in this thesis. The result of our work is a novel shock capturing scheme - a hybrid between FV and DG - with smooth transitions between low and high order fluxes according to solution smoothness estimators. We present extensive validations and test cases, specifically its interaction with multi-physics modules in FLASH such as (self-)gravity and radiative transfer. We also investigate the benefits and pitfalls of integrating end-to-end entropy stability into our numerical scheme, with special focus on highly compressible turbulent flows and shocks. Our implementation of DG in FLASH allows us to conduct preliminary yet comprehensive astrophysics simulations proving that our new solver is ready for assessments and investigations by the astrophysics community

Numerical Simulations of Shock and Rarefaction Waves Interacting With Interfaces in Compressible Multiphase Flows

Developing a highly accurate numerical framework to study multiphase mixing in high speed flows containing shear layers, shocks, and strong accelerations is critical to many scientific and engineering endeavors. These flows occur across a wide range of scales: from tiny bubbles in human tissue to massive stars collapsing. The lack of understanding of these flows has impeded the success of many engineering applications, our comprehension of astrophysical and planetary formation processes, and the development of biomedical technologies. Controlling mixing between different fluids is central to achieving fusion energy, where mixing is undesirable, and supersonic combustion, where enhanced mixing is important. Iron, found throughout the universe and a necessary component for life, is dispersed through the mixing processes of a dying star. Non-invasive treatments using ultrasound to induce bubble collapse in tissue are being developed to destroy tumors or deliver genes to specific cells. Laboratory experiments of these flows are challenging because the initial conditions and material properties are difficult to control, modern diagnostics are unable to resolve the flow dynamics and conditions, and experiments of these flows are expensive. Numerical simulations can circumvent these difficulties and, therefore, have become a necessary component of any scientific challenge. Advances in the three fields of numerical methods, high performance computing, and multiphase flow modeling are presented: (i) novel numerical methods to capture accurately the multiphase nature of the problem; (ii) modern high performance computing paradigms to resolve the disparate time and length scales of the physical processes; (iii) new insights and models of the dynamics of multiphase flows, including mixing through hydrodynamic instabilities. These studies have direct applications to engineering and biomedical fields such as fuel injection problems, plasma deposition, cancer treatments, and turbomachinery.PhDMechanical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/133458/1/marchdf_1.pd

Influence of A-Posteriori Subcell Limiting on Fault Frequency in Higher-Order DG Schemes

Soft error rates are increasing as modern architectures require increasingly small features at low voltages. Due to the large number of components used in HPC architectures, these are particularly vulnerable to soft errors. Hence, when designing applications that run for long time periods on large machines, algorithmic resilience must be taken into account. In this paper we analyse the inherent resiliency of a-posteriori limiting procedures in the context of the explicit ADER DG hyperbolic PDE solver ExaHyPE. The a-posteriori limiter checks element-local high-order DG solutions for physical admissibility, and can thus be expected to also detect hardware-induced errors. Algorithmically, it can be interpreted as element-local checkpointing and restarting of the solver with a more robust finite volume scheme on a fine subgrid. We show that the limiter indeed increases the resilience of the DG algorithm, detecting and correcting particularly those faults which would otherwise lead to a fatal failure.Comment: 2018 IEEE/ACM 8th Workshop on Fault Tolerance for HPC at eXtreme Scale (FTXS

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

Proceedings of the Third International Mobile Satellite Conference (IMSC 1993)

Satellite-based mobile communications systems provide voice and data communications to users over a vast geographic area. The users may communicate via mobile or hand-held terminals, which may also provide access to terrestrial cellular communications services. While the first and second International Mobile Satellite Conferences (IMSC) mostly concentrated on technical advances, this Third IMSC also focuses on the increasing worldwide commercial activities in Mobile Satellite Services. Because of the large service areas provided by such systems, it is important to consider political and regulatory issues in addition to technical and user requirements issues. Topics covered include: the direct broadcast of audio programming from satellites; spacecraft technology; regulatory and policy considerations; advanced system concepts and analysis; propagation; and user requirements and applications