Universitätsbibliothek

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.

In this thesis, the broad band emission, especially in the gamma-ray and radio band, of the active galaxy IC 310 located in the Perseus cluster of galaxies was investigated. The main experimental methods were Cherenkov astronomy using the MAGIC telescopes and high resolution very long baseline interferometry (VLBI) at radio frequencies (MOJAVE, EVN). Additionally, data of the object in different energy bands were studied and a multi-wavelength campaign has been organized and conducted. During the campaign, an exceptional bright gamma-ray flare at TeV energies was found with the MAGIC telescopes. The results were compared to theoretical acceleration and emission models for explaining the high energy radiation of active galactic nuclei. Many open questions regarding the particle acceleration to very high energies in the jets of active galactic nuclei, the particle content of the jets, or how the jets are launched, were addressed in this thesis by investigating the variability of IC 310 in the very high energy band. It is argued that IC310 was originally mis-classified as a head-tail radio galaxy. Instead, it shows a variability behavior in the radio, X-ray, and gamma-ray band similar to the one found for blazars. These are active galactic nuclei that are characterized by flux variability in all observed energy bands and at all observed time scales. They are viewed at a small angle between the jet axis and the line-of-sight. Thus, strong relativistic beaming influences the variability properties of blazars. Observations of IC 310 with the European VLBI Network helped to find limits for the angle between the jet axis and the line-of-sight, namely 10 deg - 20 deg. This places IC 310 at the borderline between radio galaxies (larger angles) and blazars (smaller angles). During the gamma-ray outburst detected at the beginning of the multi-wavelength campaign, flux variability as short as minutes was measured. The spectrum during the flare can be described by a simple power-law function over two orders of magnitude in energy up to ~10 TeV. Compared to previous observations, no significant variability of the spectral shape was found. Together with the constraint on the viewing angle, this challenges the currently accepted models for particle acceleration at shock waves in the jets. Alternative models, such as stars moving through the jets, mini-jets in the jet caused, e.g., by reconnection events, or gap acceleration in a pulsar-like magnetosphere around the black hole were investigated. It was found that only the latter can explain all observational findings, which at least suggests that it could even be worthwhile to reconsider published investigations of AGN with this new knowledge in mind. The first multi-wavelength campaign was successfully been conducted in 2012/2013, including ground-based as well as space-based telescopes in the radio, optical, ultraviolet, X-ray, and gamma-ray energy range. No pronounced variability was found after the TeV flare in any energy band. The X-ray data showed a slightly harder spectrum when the emission was brighter. The long-term radio light curve indicated a flickering flux variability, but no strong hint for a new jet component was found from VLBI images of the radio jet. In any case, further analysis of the existing multi-wavelength data as well as complimentary measurements could provide further exciting insights, e.g., about the broad band spectral energy distribution. Overall, it can be stated that IC 310 is a key object for research of active galactic nuclei in the high-energy band due to its proximity and its peculiar properties regarding flux variability and spectral behavior. Such objects are ideally suited for studying particle acceleration, jet formation, and other physical effects and models which are far from being fully understood.

Die Herkunft hochenergetischer solarer Teilchen konnte in den vergangenen Jahren eindeutig auf Schockbeschleunigung an koronalen Masseauswürfen zurückgeführt werden. Durch resonante Interaktionen zwischen Wellen und Teilchen werden zum einen geladene Teilchen unter Veränderung ihrer Energie gestreut, zum anderen wird die Dynamik der Plasmawellen in solchen Beschleunigungsregionen durch diese Prozesse von selbstgenerierten Wellenmoden maßgeblich beeinflusst. Mittels numerischer Modellierungen wurden im Rahmen dieser Arbeit die grundlegenden physikalischen Regimes der Turbulenz und des Teilchentransports beschrieben. Die Simulation der Plasmadynamik bedient sich der Methodik der Magnetohydrodynamik, wohingegen kinetische Einzelteilchen durch die elementaren Bewegungsgleichungen der Elektrodynamik berechnet werden. Es konnten die Turbulenztheorien von Goldreich und Sridhar unter heliosphärischen Bedingungen bei drei solaren Radien bestätigt werden. Vor allem zeigten sich Hinweise für das Erreichen der kritischen Balance, einem Schlüsselparameter dieser Theorien. Weiterhin werden Ergebnisse der dynamischen Entwicklung angeregter Wellenmoden präsentiert, in denen die Bedeutsamkeit für die gesamte Turbulenz gezeigt werden konnte. Als zentraler Prozess bei hohen Energien hat sich das wave-steepening herausgestellt, das als effizienter Energietransportmechanismus in paralleler Richtung zum Hintergrundmagnetfeld identifiziert wurde und somit turbulente Strukturen bei hohen parallelen Wellenzahlen erklärt, deren Entstehung das Goldreich-Sridhar Modell nicht beschreiben kann. Darüber hinaus wurden grundlegende Erkenntnisse über die quasilineare Theorie des Teilchentransports erzielt. Im Speziellen konnte ein tieferes Verständnis für die Interpretation der Diffusionskoeffizienten von Welle-Teilchen Wechselwirkungen erlangt werden. Simulationen zur Streuung an angeregten Wellenmoden zeigten erstmals komplexe resonante Strukturen die im Rahmen analytischer Modelle nicht mehr adäquat beschrieben werden können.

Cold dark matter constitutes a basic tenet of modern cosmology, essential for our understanding of structure formation in the Universe. Since its first discovery by means of spectroscopic observations of the dynamics of the Coma cluster some 80 years ago, mounting evidence of its gravitational pull and its impact on the geometry of space-time has build up across a wide range of scales, from galaxies to the entire Hubble flow. The apparent lack of electromagnetic coupling and independent measurements of the energy density of baryonic matter from the primordial abundances of light elements show the non-baryonic nature of dark matter, and its clustering properties prove that it is cold, i.e. that it has a temperature lower than its mass during the time of radiation-matter equality. A generic particle candidate for cold dark matter are weakly interacting massive particles at the electroweak symmetry-breaking scale, such as the neutralinos in R-parity conserving supersymmetry. Such particles would naturally freeze-out with a cosmologically relevant relic density at early times in the expanding Universe. Subsequent clustering of matter would recover annihilation interactions between the dark matter particles to some extent and thus lead to potentially observable high-energy emission from the decaying unstable secondaries produced in annihilation events. The spectra of the secondaries would permit a determination of the mass and annihilation cross section, which are crucial for the microphysical identification of the dark matter. This the central motivation for indirect dark matter searches. However, presently neither the indirect searches, nor the complementary direct searches based on the detection of elastic scattering events, nor the production of candidate particles in collider experiments, has yet provided unequivocal evidence for dark matter. This does not come as a surprise, since the dark matter particles interact only through weak interactions and therefore the corresponding secondary emission must be extremely faint. It turns out that even for the strongest mass concentrations in the Universe, the dark matter annihilation signal is expected to not exceed the level of competing astrophysical sources. Thus, the discrimination of the putative dark matter annihilation signal from the signals of the astrophysical inventory has become crucial for indirect search strategies. In this thesis, a novel search strategy will be developed and exemplified in which target selection across a wide range of masses, astrophysical background estimation, and multiwavelength signatures play the key role. It turns out that the uncertainties regarding the halo profile and the boost due to surviving substructure are bigger for halos at the lower end of the observed mass scales, i.e. in the regime of dwarf galaxies and below, while astrophysical backgrounds tend to become more severe for massive dark matter halos such as clusters of galaxies. By contrast, the uncertainties due to unknown details of particle physics are invariant under changes of the halo mass. Therefore, the different scaling behaviors can be employed to significantly cut down on the uncertainties in observations of different targets covering a major part of the involved mass scales. This strategical approach was implemented in the scientific program carried out with the MAGIC telescope system. Observations of dwarf galaxies and the Virgo- and Perseus clusters of galaxies have been carried out and, at the time of writing, result in some of the most stringent constraints on weakly interacting massive particles from indirect searches. Here, the low-threshold design of the MAGIC telescope system plays a crucial role, since the bulk of the high-energy photons, produced with a high multiplicity during the fragmentation of unstable dark matter annihilation products, are emitted at energies well below the dark matter mass scale. The upper limits severely constrain less generic, but more prolific scenarios characterized by extraordinarily high annihilation efficiencies.

One key scientific program of the MAGIC telescope project is the discovery and detection of blazars. They constitute the most prominent extragalactic source class in the very high energy (VHE) Gamma-ray regime with 29 out of 34 known objects (as of April 2010). Therefore a major part of the available observation time was spent in the last years on high-frequency peaked blazars. The selection criteria were chosen to increase the detection probability. As the X-ray flux is believed to be correlated to the VHE Gamma-ray flux, only X-ray selected sources with a flux F(X) > 2 μJy at 1 keV were considered. To avoid strong attenuation of the Gamma-rays in the extragalactic infrared background, the redshift was restricted to values between z < 0.15 and z < 0.4, depending on the declination of the objects. The latter determines the zenith distance during culmination which should not exceed 30° (for z < 0.4) and 45° (for z < 0.15), respectively. Between August 2005 and April 2009, a sample of 24 X-ray selected high-frequency peaked blazars has been observed with the MAGIC telescope. Three of them were detected including 1ES 1218+304 being the first high-frequency peaked BL Lacertae object (HBL) to be discovered with MAGIC in VHE Gamma-rays. One previously detected object was not confirmed as VHE emitter in this campaign by MAGIC. A set of 20 blazars previously not detected will be treated more closely in this work. In this campaign, during almost four years ~ 450 hrs or ~ 22% of the available observation time for extragalactic objects were dedicated to investigate the baseline emission of blazars and their broadband spectral properties in this emission state. For the sample of 20 objects in a redshift range of 0.018 < z < 0.361 integral flux upper limits in the VHE range on the 99.7% confidence level (corresponding to 3 standard deviations) were calculated resulting in values between 2.9% and 14.7% of the integral flux of the Crab Nebula. As the distribution of significances of the individual objects shows a clear shift to positive values, a stacking method was applied to the sample. For the whole set of 20 objects, an excess of Gamma-rays was found with a significance of 4.5 standard deviations in 349.5 hours of effective exposure time. For the first time a signal stacking in the VHE regime turned out to be successful. The measured integral flux from the cumulative signal corresponds to 1.4% of the Crab Nebula flux above 150 GeV with a spectral index α = −3.15±0.57. None of the objects showed any significant variability during the observation time and therefore the detected signal can be interpreted as the baseline emission of these objects. For the individual objects lower limits on the broad-band spectral indices αX−Gamma between the X-ray range at 1 keV and the VHE Gamma-ray regime at 200 GeV were calculated. The majority of objects show a spectral behaviour as expected from the source class of HBLs: The energy output in the VHE regime is in general lower than in X-rays. For the stacked blazar sample the broad-band spectral index was calculated to αX−Gamma = 1.09, confirming the result found for the individual objects. Another evidence for the revelation of the baseline emission is the broad-band spectral energy distribution (SED) comprising archival as well as contemporaneous multi-wavelength data from the radio to the VHE band. The SEDs of known VHE Gamma-ray sources in low flux states matches well the SED of the stacked blazar sample.