TY - THES A1 - Burkart, Thomas T1 - Der Einfluss des fundamentalen Massenverhältnisses auf die Teilchenbeschleunigung durch Plasmainstabilitäten T1 - The influence of the fundamental mass-ratio on particle acceleration by plasma instabilities N2 - Im Rahmen dieser Arbeit wurde ein dreidimensionaler vollrelativistischer und parallelisierter Particle-in-Cell Code geschrieben, ausführlich getestet und angewandt. Der Code ACRONYM ist variabel einsetzbar und von der Genauigkeit und Stabilität her State-of-the-Art und somit konkurrenzfähig zu den sonstigen in der Astrophysik eingesetzten Codes anderer Gruppen. Die Energie bleibt bis auf einen Fehler von < 0.03% erhalten, die Divergenz des Magnetfeldes bleibt immer unter einem Wert von 10^{-12} und die Skalierung wurde mittlerweile bis zu einem Clustergröße von einigen 10000 CPUs getestet. In dieser Arbeit wurde dann, nach der Entwicklung des Codes, der Einfluss des fundamentalen Massenverhältnisses m_p/m_e auf die Teilchenbeschleunigung durch Plasmainstabilitäten untersucht. Dies ist relevant und wichtig, da in PiC-Simulationen in den allermeisten Fällen nicht mit dem realen Massenverhältnis gerechnet wird, da sonst viel zu viel Rechenleistung benötigt würde, um zu sehen, was mit den Protonen geschieht und was ihr Einfluss auf die leichten Teilchen wie Elektronen und Positronen ist. Zu diesem Zweck wurden Simulationen mit Massenverhältnissen zwischen m_p/m_e = 1.0 und 200.0 durchgeführt. Diese haben alle gemeinsam, dass periodische Randbedingungen verwendet wurden und das zur Verfügung stehende Simulationsgebiet mit jeweils zwei gegeneinander strömenden Plasmapopulationen vollständig gefüllt wurde, um jegliche Art von auftretenden Schocks auszuschließen. Die Rohdaten der einzelnen Simulationen wurden auf vielfältige Art und Weise analysiert, es wurden z.B. Schnitte durch die Teilchenverteilung erstellt, sowie ein- oder zweidimensionale Histogramme und Energieverläufe betrachtet. Dabei haben sich folgende Kernpunkte ergeben: Für Massenverhältnisse bis etwa m_p/m_e = 20 bildet sich die gesamte Zweistrom-Instabilität in nur einer Phase aus, das heißt, es bilden sich von ringförmigen Magnetfeldern umgebene Flussschläuche aus, die dann verschmelzen, bis nur noch zwei übrig sind und alle Teilchen werden über den gesamten Verlauf der Instabilität beschleunigt. Es ist damit zu folgern, dass die unterschiedlich schweren Teilchenspezies Protonen und Elektronen/Positronen durch die relativ nahe beieinander liegenden Massen noch so stark gekoppelt sind, dass sich nur eine Instabilität entwickeln kann. Bei großen Massenverhältnissen (m_p/m_e > 20) ist eine deutliche Trennung in zwei Phasen der Instabilität zu erkennen. Zuerst bilden sich wiederum Flussschläuche aus, diese verschmelzen miteinander (zu zweien oder mehr), bevor der erste Teil der Instabilität abflaut. Anschließend entstehen wieder ringförmige Magnetfelder und Flussschläuche, von denen einer meist deutlich stärker ist als all die anderen, das bedeutet, dass dieser von stärkeren Magnetfeldern umgeben ist und eine höhere Teilchendichte aufweist. Im Rahmen dieser zweigeteilten Instabilität werden die Elektronen und Positronen nur in der ersten Phase signifikant beschleunigt, die deutlich schwereren Protonen gewinnen über den gesamten Zeitraum Energie. Die höchstenergetischen Teilchen erreichen im Ruhesystem der jeweiligen Plasmapopulation Werte um gamma = 250. Man kann daraus für zukünftige Untersuchungen mit Hilfe von Particle-in-Cell Codes den Schluss ziehen, dass Rückschlüsse auf das tatsächliche Verhalten beim realen Massenverhältnis von m_p/m_e = 1836.2 nur aus den Simulationen mit m_p/m_e >> 20 gezogen werden können, da die starke Kopplung der leichten und schweren Teilchen bei kleineren Massenverhältnissen die Ergebnisse sehr stark beeinflusst. Es wurde anhand der gemessenen Zeitpunkte der Instabilitätsmaxima eine Extrapolation durchgeführt, die zeigt, dass die Instabilität beim realen Massenverhältnis etwa bei t = 1400 omega_{pe}^{-1} auftreten würde. Um dies wirklich zu simulieren müsste allerdings mehr als die 1000-fache Anzahl an CPU-Stunden aufgewandt werden. Des weiteren wurde eine Maxwell-Jüttner-Verteilung an die Teilchenverteilungen der einzelnen Simulationen auf dem Höhepunkt der Instabilität gefittet, um sowohl die neue Temperatur des Plasmas als auch die Beschleunigungseffizienz des Prozesses zu berechnen. Die Temperatur erhöht sich demnach durch die Instabilität von etwa 10^8K auf 10^{10} bis 10^{11}K, der Anteil suprathermischer Teilchen beträgt 2 bis 4%. N2 - In this thesis a three-dimensional, fully relativistic and parallelised Particle-in-Cell Code was developed, tested and used for astrophysical purposes. The Code ACRONYM can be used for a variety of different scenarios, it is state-of-the-art in matters of stability and accuracy. After the development the code was used to investigate the influence of the fundamental mass ratio m_p/m_e on particle acceleration by plasma instabilities. This is important, because usually in PiC-simulations the mass ratio used isn't the real one m_p/m_e = 1836.2, because this would take too much CPU-time in order to see what happens to the protons and what is their influence on the lighter particles like electrons and positrons. For this purpose simulations with mass ratios between 1.0 and 200.0 have been performed. They all have in common that periodic boundary conditions were used and that the whole computational domain has been filled with particles that are counterstreaming along the z-direction with gamma approximately 10 each in order to exclude any development of shocks. The resulting main issues are the following: For mass ratios below m_p/m_e approximately 20 the whole instability develops in only one phase, i.e. current filaments surrounded by circular magnetic fields develop and merge together. All particles are accelerated over the whole run, so one can conclude that the different species are still strongly coupled because of the very similar masses of electrons/positrons and the protons and therefore only one instability can arise. For higher mass ratios a distinctive separation of the instability in two phases is observable. First some flux tubes develop and merge until the first phase is over. Afterwards new magnetic fields and flux tubes are arising, where one of them usually is particularly strong compared to the others, i.e. it is surrounded by stronger magnetic fields and holds a much higher particle density. In the context of this split instability, the electrons and positrons are getting accelerated significantly only in the first phase, the much heavier protons gain energy over the whole time. One can therefore conclude for future investigations with PiC codes that informations about the behaviour at the realistic mass ratio of m_p/m_e = 1836.2 can only be gained from the simulations with m_p/m_e >> 20 because of the strong coupling of the light and heavy particles at low mass ratios. An extrapolation to the real mass ratio shows that the peak of the instability would occur approximately seven times later than the runtime of the longest simulation at about t = 1400 omega_{pe}^{-1}, but in order to realize this, at least 1000 times the now used CPU-hours would be necessary. Furthermore the acceleration efficiency for this process was calculated by fitting a Maxwell-Jüttner-Distribution to the particle distribution from the simulations during the peak of the instabilities. The calculated fraction of superthermal particles is in the range of 2 to 4% and the temperatures of the plasma streams rise from 10^8 at the beginning of the simulations to values around 10^10 to 10^11K. KW - Astrophysik KW - Plasmaschwingung KW - Teilchenbeschleunigung KW - Numerisches Modell KW - Particle-in-Cell Code KW - Energiereiches Teilchen KW - Instabilität KW - Plasma KW - Plasmaaufheizung KW - Teilchenbewegung KW - Particle-in-Cell Code Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-56636 ER - TY - THES A1 - Lange, Sebastian T1 - Turbulenz und Teilchentransport in der Heliosphäre - Simulationen von inkompressiblen MHD-Plasmen und Testteilchen - T1 - Turbulence and particle transport within the heliosphere - simulations of incompressible MHD-plasmas and test particles - N2 - 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. N2 - In the past years, the origin of high energetic solar particles could be clearly connected to shock acceleration at coronal mass ejections. Caused by resonant wave-particle interactions, on the one hand, the particles change their energy because of scattering, on the other hand, the dynamics of plasma waves in such acceleration regions are significantly influenced by these processes through self--generated wave modes. In this dissertation, the basic physical regime of turbulence and particle transport were described via numerical modeling. The simulation of the plasma dynamics uses the methodology of magnetohydrodynamics, whereas the kinetic description of single particles is calculated by elementary electrodynamic equations of motion. The common plasma turbulence theories by Goldreich and Sridhar could be confirmed by simulations resembling conditions at three solar radii. Foremost, evidence for the critical balance has been found, which is a key parameter of these theories. Furthermore, results of the dynamic evolution of amplified wavemodes are presented, which are very important for the general turbulence development. In this context, the wave-steepening was identified as a central process, which is an efficient energy transport mechanism in parallel direction to the magnetic background field. This explains turbulent structures at high parallel wavenumbers, which are not described by the Goldreich-Sridhar model. Moreover, a fundamental understanding of the quasilinear theory of particle transport has been achieved. Specifically, more detailed insight into the interpretation of the diffusion coefficients of wave-particle interactions could be obtained. For the first time, simulations of particle scattering at amplified wave modes showed complex resonant structures, which cannot be described by analytical approaches adequately. KW - Heliosphäre KW - Sonnenwind KW - Turbulente Strömung KW - Transportprozess KW - Teilchentransport KW - turbulence KW - high-energy particles KW - heliosphere KW - magnetohydrodynamics KW - Numerisches Modell KW - Turbulenz KW - Energiereiches Teilchen KW - Magnetohydrodynamik Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-74012 ER -