TY - THES A1 - Müller, Melanie J. I. T1 - Bidirectional transport by molecular motors T1 - Bidirektionaler Transport durch molekulare Motoren N2 - In biological cells, the long-range intracellular traffic is powered by molecular motors which transport various cargos along microtubule filaments. The microtubules possess an intrinsic direction, having a 'plus' and a 'minus' end. Some molecular motors such as cytoplasmic dynein walk to the minus end, while others such as conventional kinesin walk to the plus end. Cells typically have an isopolar microtubule network. This is most pronounced in neuronal axons or fungal hyphae. In these long and thin tubular protrusions, the microtubules are arranged parallel to the tube axis with the minus ends pointing to the cell body and the plus ends pointing to the tip. In such a tubular compartment, transport by only one motor type leads to 'motor traffic jams'. Kinesin-driven cargos accumulate at the tip, while dynein-driven cargos accumulate near the cell body. We identify the relevant length scales and characterize the jamming behaviour in these tube geometries by using both Monte Carlo simulations and analytical calculations. A possible solution to this jamming problem is to transport cargos with a team of plus and a team of minus motors simultaneously, so that they can travel bidirectionally, as observed in cells. The presumably simplest mechanism for such bidirectional transport is provided by a 'tug-of-war' between the two motor teams which is governed by mechanical motor interactions only. We develop a stochastic tug-of-war model and study it with numerical and analytical calculations. We find a surprisingly complex cooperative motility behaviour. We compare our results to the available experimental data, which we reproduce qualitatively and quantitatively. N2 - In biologischen Zellen transportieren molekulare Motoren verschiedenste Frachtteilchen entlang von Mikrotubuli-Filamenten. Die Mikrotubuli-Filamente besitzen eine intrinsische Richtung: sie haben ein "Plus-" und ein "Minus-"Ende. Einige molekulare Motoren wie Dynein laufen zum Minus-Ende, während andere wie Kinesin zum Plus-Ende laufen. Zellen haben typischerweise ein isopolares Mikrotubuli-Netzwerk. Dies ist besonders ausgeprägt in neuronalen Axonen oder Pilz-Hyphen. In diesen langen röhrenförmigen Ausstülpungen liegen die Mikrotubuli parallel zur Achse mit dem Minus-Ende zum Zellkörper und dem Plus-Ende zur Zellspitze gerichtet. In einer solchen Röhre führt Transport durch nur einen Motor-Typ zu "Motor-Staus". Kinesin-getriebene Frachten akkumulieren an der Spitze, während Dynein-getriebene Frachten am Zellkörper akkumulieren. Wir identifizieren die relevanten Längenskalen und charakterisieren das Stauverhalten in diesen Röhrengeometrien mit Hilfe von Monte-Carlo-Simulationen und analytischen Rechnungen. Eine mögliche Lösung für das Stauproblem ist der Transport mit einem Team von Plus- und einem Team von Minus-Motoren gleichzeitig, so dass die Fracht sich in beide Richtungen bewegen kann. Dies wird in Zellen tatsächlich beobachtet. Der einfachste Mechanismus für solchen bidirektionalen Transport ist ein "Tauziehen" zwischen den beiden Motor-Teams, das nur mit mechanischer Interaktion funktioniert. Wir entwickeln ein stochastisches Tauzieh-Modell, das wir mit numerischen und analytischen Rechnungen untersuchen. Es ergibt sich ein erstaunlich komplexes Motilitätsverhalten. Wir vergleichen unsere Resultate mit den vorhandenen experimentellen Daten, die wir qualitativ und quantitativ reproduzieren. KW - molekulare Motoren KW - bidirektionaler intrazellulärer Transport KW - Tauziehen KW - stochastische Prozesse KW - kooperative Phänomene KW - molecular motors KW - bidirectional intracellular transport KW - tug-of-war KW - stochastic processes KW - cooperative phenomena Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-18715 ER - TY - THES A1 - Löffler, Frank T1 - Numerical simulations of neutron star - black hole mergers T1 - Computersimulationen von Zusammenstößen eines schwarzen Lochs und eines Neutronensterns N2 - Collisions of black holes and neutron stars, named mixed binaries in the following, are interesting because of at least two reasons. Firstly, it is expected that they emit a large amount of energy as gravitational waves, which could be measured by new detectors. The form of those waves is expected to carry information about the internal structure of such systems. Secondly, collisions of such objects are the prime suspects of short gamma ray bursts. The exact mechanism for the energy emission is unknown so far. In the past, Newtonian theory of gravitation and modifications to it were often used for numerical simulations of collisions of mixed binary systems. However, near to such objects, the gravitational forces are so strong, that the use of General Relativity is necessary for accurate predictions. There are a lot of problems in general relativistic simulations. However, systems of two neutron stars and systems of two black holes have been studies extensively in the past and a lot of those problems have been solved. One of the remaining problems so far has been the use of hydrodynamic on excision boundaries. Inside excision regions, no evolution is carried out. Such regions are often used inside black holes to circumvent instabilities of the numerical methods near the singularity. Methods to handle hydrodynamics at such boundaries have been described and tests are shown in this work. One important test and the first application of those methods has been the simulation of a collapsing neutron star to a black hole. The success of these simulations and in particular the performance of the excision methods was an important step towards simulations of mixed binaries. Initial data are necessary for every numerical simulation. However, the creation of such initial data for general relativistic situations is in general very complicated. In this work it is shown how to obtain initial data for mixed binary systems using an already existing method for initial data of two black holes. These initial data have been used for evolutions of such systems and problems encountered are discussed in this work. One of the problems are instabilities due to different methods, which could be solved by dissipation of appropriate strength. Another problem is the expected drift of the black hole towards the neutron star. It is shown, that this can be solved by using special gauge conditions, which prevent the black hole from moving on the computational grid. The methods and simulations shown in this work are only the starting step for a much more detailed study of mixed binary system. Better methods, models and simulations with higher resolution and even better gauge conditions will be focus of future work. It is expected that such detailed studies can give information about the emitted gravitational waves, which is important in view of the newly built gravitational wave detectors. In addition, these simulations could give insight into the processes responsible for short gamma ray bursts. N2 - Zusammenstöße eines schwarzen Lochs und eines Neutronensterns, im Folgenden "gemischte Zusammenstöße" genannt, sind aus wenigstens zwei Gründen interessant. Erstens wird erwartet, dass dabei große Mengen Energie als Gravitationswellen freigesetzt werden und diese mit neuen Detektoren gemessen werden können. Die Form dieser Wellen verrät viel über die Beschaffenheit eines solchen Systems und stellt neben elektromagnetischen Wellen eine wichtige Informationsquelle dar. Zweitens sind Zusammenstöße von kompakten Objekten wie Neutronensternen und schwarze Löchern sehr wahrscheinlich die Ursache sogenannter kurzer Gammastrahlungsblitze. Deren genauer Mechanismus für die Umwandlung der gewaltigen Energiemengen, die bei diesen Blitzen ausgesandt werden, ist jedoch bisher unbekannt. Computersimulationen von Zusammenstößen eines gemischten Systems wurden bisher oft unter Benutzung der Newtonschen Gravitationstheorie, bzw. Korrekturen dazu, durchgeführt. In der Nähe so kompakte Objekte wie schwarzer Löcher oder Neutronensterne ist jedoch die Gravitationswirkung so stark, dass Näherungen wie die erwähnten Korrekturen der Newtonschen Gravitationstheorie zu ungenau sind. Eine Benutzung der allgemeinen Relativitätstheorie ist daher für dieses Problem unumgänglich. Die Probleme allgemein-relativistischer Simulationen sind vielfältig. Jedoch wurden Binärsysteme zweier schwarzer Löcher und zweier Neutronensterne schon eingehend untersucht, und so viele Probleme, die auch Simulationen gemischter Systeme betreffen, gelöst. Eins der bisher ausstehenden Probleme war die Behandlung der Hydrodynamik an Ausschneiderändern; Rändern zu Gebieten, die in der Zeitentwicklung der Simulation ignoriert werden. Solche Ränder werden zum Beispiel innerhalb eines schwarzen Lochs benutzt, um Instabilitäten des Programms in der Nähe der Singularität zu vermeiden. Methoden, solche Ränder zu behandeln wurden in der Arbeit entwickelt, getestet und gezeigt, dass sie verlässlich arbeiten. Ein wichtiger Test für diese Methoden, der gleichzeitig der Gewinnung neuer Erkenntnisse diente, war deren Anwendung auf Simulationen von zu schwarzen Löchern kollabierenden, rotierenden Sternen. Der Erfolg, diese Simulationen ohne Probleme mit den erwähnten Methoden durchzuführen, war ein wichtiger Schritt zu Simulationen gemischter Binärsysteme. Für Computersimulationen sind Anfangsdaten notwendig, die das gewünschte Problem beschreiben. Die Erstellung solcher Anfangsdaten ist jedoch unter Benutzung der allgemeinen Relativitätstheorie ausser in Spezialfällen sehr komplex. Wir zeigen, wie man einen schon vorhandenen Algorithmus für Anfangsdaten für zwei schwarze Löcher ändern kann, um Anfangsdaten für ein gemischtes Binärsystem zu erhalten. Diese Anfangsdaten wurden für Simulationen eines gemischten Binärsystems benutzt. Während dieser Simulationen traten mehrere Probleme auf. Zwei dieser Probleme waren numerische Instabilitäten unterschiedlicher Herkunft. Beide konnten jedoch mit angepasst starker Dissipation (der künstliche Entnahme von hochfrequenter Energie aus dem System) unterdrückt werden. Ein weiteres Problem war die erwartete Bewegung des schwarzen Lochs in Richtung des Neutronensterns. Da ein Teil des Simulationsgebietes innerhalb des schwarzen Lochs ausgeschnitten wird und das verwendete Programm bewegte Ausschneidegebiete nicht behandeln kann, darf sich das schwarze Loch jedoch auf dem Gitter kaum bewegen. Wir haben dieses Problem durch eine an das Problem angepasste Eichbedingung gelöst, die auf Bewegungen des scheinbaren Horizons reagiert und die Position des schwarzen Lochs auf diese Weise nahezu konstant hält. Die Methoden und Simulationen dieser Arbeit sind nur der Anfangspunkt einer ausführlichen Studie von Binärsystemen eines schwarzen Lochs und eines Neutronensterns. Bessere Methoden, Modelle und Simulationen mit höherer Auflösung und besser an das System angepassten Koordinaten werden Mittelpunkt zukünftiger Arbeit sein. Es wird erwartet, dass solche detailierten Studien Erkenntnisse über die abgestrahlten Gravitationswellen liefern, die gerade in Hinblick auf die neuen Gravitationswellendetektoren wichtig sind. Weiterhin könnten diese Simulationen dabei helfen, die Prozesse, die kurze Gammastrahlungsblitze hervorrufen, und über die im Moment kaum etwas bekannt ist, aufzuklären. KW - Relativistische Astrophysik KW - Neutronensterne KW - schwarze Löcher KW - neutron stars KW - black holes Y1 - 2005 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-7743 ER - TY - THES A1 - Miteva, Rositsa Stoycheva T1 - Electron acceleration at localized wave structures in the solar corona T1 - Elektronenbeschleunigung an lokalen Wellenstrukturen in der Sonnenkorona N2 - Our dynamic Sun manifests its activity by different phenomena: from the 11-year cyclic sunspot pattern to the unpredictable and violent explosions in the case of solar flares. During flares, a huge amount of the stored magnetic energy is suddenly released and a substantial part of this energy is carried by the energetic electrons, considered to be the source of the nonthermal radio and X-ray radiation. One of the most important and still open question in solar physics is how the electrons are accelerated up to high energies within (the observed in the radio emission) short time scales. Because the acceleration site is extremely small in spatial extent as well (compared to the solar radius), the electron acceleration is regarded as a local process. The search for localized wave structures in the solar corona that are able to accelerate electrons together with the theoretical and numerical description of the conditions and requirements for this process, is the aim of the dissertation. Two models of electron acceleration in the solar corona are proposed in the dissertation: I. Electron acceleration due to the solar jet interaction with the background coronal plasma (the jet--plasma interaction) A jet is formed when the newly reconnected and highly curved magnetic field lines are relaxed by shooting plasma away from the reconnection site. Such jets, as observed in soft X-rays with the Yohkoh satellite, are spatially and temporally associated with beams of nonthermal electrons (in terms of the so-called type III metric radio bursts) propagating through the corona. A model that attempts to give an explanation for such observational facts is developed here. Initially, the interaction of such jets with the background plasma leads to an (ion-acoustic) instability associated with growing of electrostatic fluctuations in time for certain range of the jet initial velocity. During this process, any test electron that happen to feel this electrostatic wave field is drawn to co-move with the wave, gaining energy from it. When the jet speed has a value greater or lower than the one, required by the instability range, such wave excitation cannot be sustained and the process of electron energization (acceleration and/or heating) ceases. Hence, the electrons can propagate further in the corona and be detected as type III radio burst, for example. II. Electron acceleration due to attached whistler waves in the upstream region of coronal shocks (the electron--whistler--shock interaction) Coronal shocks are also able to accelerate electrons, as observed by the so-called type II metric radio bursts (the radio signature of a shock wave in the corona). From in-situ observations in space, e.g., at shocks related to co-rotating interaction regions, it is known that nonthermal electrons are produced preferably at shocks with attached whistler wave packets in their upstream regions. Motivated by these observations and assuming that the physical processes at shocks are the same in the corona as in the interplanetary medium, a new model of electron acceleration at coronal shocks is presented in the dissertation, where the electrons are accelerated by their interaction with such whistlers. The protons inflowing toward the shock are reflected there by nearly conserving their magnetic moment, so that they get a substantial velocity gain in the case of a quasi-perpendicular shock geometry, i.e, the angle between the shock normal and the upstream magnetic field is in the range 50--80 degrees. The so-accelerated protons are able to excite whistler waves in a certain frequency range in the upstream region. When these whistlers (comprising the localized wave structure in this case) are formed, only the incoming electrons are now able to interact resonantly with them. But only a part of these electrons fulfill the the electron--whistler wave resonance condition. Due to such resonant interaction (i.e., of these electrons with the whistlers), the electrons are accelerated in the electric and magnetic wave field within just several whistler periods. While gaining energy from the whistler wave field, the electrons reach the shock front and, subsequently, a major part of them are reflected back into the upstream region, since the shock accompanied with a jump of the magnetic field acts as a magnetic mirror. Co-moving with the whistlers now, the reflected electrons are out of resonance and hence can propagate undisturbed into the far upstream region, where they are detected in terms of type II metric radio bursts. In summary, the kinetic energy of protons is transfered into electrons by the action of localized wave structures in both cases, i.e., at jets outflowing from the magnetic reconnection site and at shock waves in the corona. N2 - Die Sonne ist ein aktiver Stern, was sich nicht nur in den allseits bekannten Sonnenflecken, sondern auch in Flares manifestiert. Während Flares wird eine große Menge gespeicherter, magnetischer Energie in einer kurzen Zeit von einigen Sekunden bis zu wenigen Stunden in der Sonnenkorona freigesetzt. Dabei werden u.a. energiereiche Elektronen erzeugt, die ihrerseits nichtthermische Radio- und Röntgenstrahlung, wie sie z.B. am Observatorium für solare Radioastronomie des Astrophysikalischen Instituts Potsdam (AIP) in Tremsdorf und durch den NASA-Satelliten RHESSI beobachtet werden, erzeugen. Da diese Elektronen einen beträchtlichen Anteil der beim Flare freigesetzten Energie tragen, ist die Frage, wie Elektronen in kurzer Zeit auf hohe Energien in der Sonnenkorona beschleunigt werden, von generellem astrophysikalischen Interesse, da solche Prozesse auch in anderen Sternatmosphären und kosmischen Objekten, wie z.B. Supernova-Überresten, stattfinden. In der vorliegenden Dissertation wird die Elektronenbeschleunigung an lokalen Wellenstrukturen im Plasma der Sonnenkorona untersucht. Solche Wellen treten in der Umgebung der magnetischen Rekonnektion, die als ein wichtiger Auslöser von Flares angesehen wird, und in der Nähe von Stoßwellen, die infolge von Flares erzeugt werden, auf. Generell werden die Elektronen als Testteilchen behandelt. Sie werden durch ihre Wechselwirkung mit den elektrischen und magnetischen Feldern, die mit den Plasmawellen verbunden sind, beschleunigt. Infolge der magnetischen Rekonnektion als Grundlage des Flares werden starke Plasmaströmungen (sogenannte Jets) erzeugt. Solche Jets werden im Licht der weichen Röntgenstrahlung, wie z.B. durch den japanischen Satelliten YOHKOH, beobachtet. Mit solchen Jets sind solare Typ III Radiobursts als Signaturen von energiereichen Elektronenstrahlen in der Sonnenkorona verbunden. Durch die Wechselwirkung eines Jets mit dem umgebenden Plasma werden lokal elektrische Felder erzeugt, die ihrerseits Elektronen beschleunigen können. Dieses hier vorgestellte Szenarium kann sehr gut die Röntgen- und Radiobeobachtungen von Jets und den damit verbundenen Elektronenstrahlen erklären. An koronalen Stoßwellen, die infolge Flares entstehen, werden Elektronen beschleunigt, deren Signatur man in der solaren Radiostrahlung in Form von sogenannten Typ II Bursts beobachten kann. Stoßwellen in kosmischen Plasmen können mit Whistlerwellen (ein spezieller Typ von Plasmawellen) verbunden sein. In der vorliegenden Arbeit wird ein Szenarium vorgestellt, das aufzeigt, wie solche Whistlerwellen an koronalen Stoßwellen erzeugt werden und durch ihre resonante Wechselwirkung mit den Elektronen dieselben beschleunigen. Dieser Prozess ist effizienter als bisher vorgeschlagene Mechanismen und kann deshalb auch auf andere Stoßwellen im Kosmos, wie z.B. an Supernova-Überresten, zur Erklärung der dort erzeugten Radio- und Röntgenstrahlung dienen. KW - Elektronenbeschleunigung KW - Sonnenkorona KW - Jets KW - Stoßwellen KW - Nichtlineare Wellen KW - Electron acceleration KW - Solar corona KW - Jets KW - Shock waves KW - Nonlinear waves Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-14775 ER - TY - THES A1 - Nickelt-Czycykowski, Iliya Peter T1 - Aktive Regionen der Sonnenoberfläche und ihre zeitliche Variation in zweidimensionaler Spektro-Polarimetrie T1 - Active regions on the solar surface and their temporal variation in two-dimensional spectropolarimetry N2 - Die Arbeit beschreibt die Analyse von Beobachtungen zweier Sonnenflecken in zweidimensionaler Spektro-Polarimetrie. Die Daten wurden mit dem Fabry-Pérot-Interferometer der Universität Göttingen am Vakuum-Turm-Teleskop auf Teneriffa erfasst. Von der aktiven Region NOAA 9516 wurde der volle Stokes-Vektor des polarisierten Lichts in der Absorptionslinie bei 630,249 nm in Einzelaufnahmen beobachtet, und von der aktiven Region NOAA 9036 wurde bei 617,3 nm Wellenlänge eine 90-minütige Zeitserie des zirkular polarisierten Lichts aufgezeichnet. Aus den reduzierten Daten werden Ergebniswerte für Intensität, Geschwindigkeit in Beobachtungsrichtung, magnetische Feldstärke sowie verschiedene weitere Plasmaparameter abgeleitet. Mehrere Ansätze zur Inversion solarer Modellatmosphären werden angewendet und verglichen. Die teilweise erheblichen Fehlereinflüsse werden ausführlich diskutiert. Das Frequenzverhalten der Ergebnisse und Abhängigkeiten nach Ort und Zeit werden mit Hilfe der Fourier- und Wavelet-Transformation weiter analysiert. Als Resultat lässt sich die Existenz eines hochfrequenten Bandes für Geschwindigkeitsoszillationen mit einer zentralen Frequenz von 75 Sekunden (13 mHz) bestätigen. In größeren photosphärischen Höhen von etwa 500 km entstammt die Mehrheit der damit zusammenhängenden Schockwellen den dunklen Anteilen der Granulen, im Unterschied zu anderen Frequenzbereichen. Die 75-Sekunden-Oszillationen werden ebenfalls in der aktiven Region beobachtet, vor allem in der Lichtbrücke. In den identifizierten Bändern oszillatorischer Power der Geschwindigkeit sind in einer dunklen, penumbralen Struktur sowie in der Lichtbrücke ausgeprägte Strukturen erkennbar, die sich mit einer Horizontalgeschwindigkeit von 5-8 km/s in die ruhige Sonne bewegen. Diese zeigen einen deutlichen Anstieg der Power, vor allem im 5-Minuten-Band, und stehen möglicherweise in Zusammenhang mit dem Phänomen der „Evershed-clouds“. Eingeschränkt durch ein sehr geringes Signal-Rausch-Verhältnis und hohe Fehlereinflüsse werden auch Magnetfeldvariationen mit einer Periode von sechs Minuten am Übergang von Umbra zu Penumbra in der Nähe einer Lichtbrücke beobachtet. Um die beschriebenen Resultate zu erzielen, wurden bestehende Visualisierungsverfahren der Frequenzanalyse verbessert oder neu entwickelt, insbesondere für Ergebnisse der Wavelet-Transformation. N2 - The publication describes the analysis of two sunspot observations in two-dimensional spectropolarimetry. The data was obtained with the Fabry-Pérot-interferometer of Göttingen University at the German Vacuum Tower Telescope on Tenerife. Of the active region NOAA 9516 the full Stokes vector of polarised light was observed in the absorption line at 630.249 nm in single scans. A ninety minute time series of circular polarised light of the active region NOAA 9036 was observed at 617.3 nm wavelength. From the reduced data results for intensity, line-of-sight velocity, magnetic field strength as well as several other plasma parameters are inferred. Different approaches to solar atmosphere model inversion are applied and compared. The significant influence of errors is discussed in detail. The frequency behaviour of the results and spatial and temporal dependencies are further analysed by Fourier and wavelet transformation. As a result the existence of a high frequency band of velocity oscillations with a central frequency of about 75-seconds (13 mHz) can be confirmed. In greater heights of about 500 km the majority of the corresponding shock waves are derived from darker parts of the granules in contrast to the dominant five-minute-oscillations. 75-second-oscillations can also be observed in the active region, especially in the light bridge. In the identified bands of oscillatory velocity power, distinct structures become visible in a penumbral dark structure as well as in the light bridge that move into the quiet sun with a horizontal speed of 5-8 km/s. They show an increase in power, mostly 5-minute-band, and may be related to the Evershed cloud phenomenon. Under the constraint of a very low signal-to-noise ratio and high error influence, magnetic field variations of a 6-minute period are also observed in an umbral-penumbral transition area close to a light bridge. To derive these results, existing visualisation methods for frequency analysis where improved or newly developed, especially so for wavelet transform results. KW - Sonne: Oszillationen KW - Sonne: Sonnenflecken KW - Instrumente: Polarimeter KW - Instrumente: Spektrographen KW - Frequenzanalyse KW - Sun: oszillations KW - Sun: sunspots KW - instrumentation: polarimeters KW - instrumentation: spectrographs KW - frequency analysis Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-25524 ER - TY - THES A1 - Kellermann, Thorsten T1 - Accurate numerical relativity simulations of non-vacuumspace-times in two dimensions and applications to critical collapse T1 - Exakte numerisch relativistische Simulationen der Nicht-Vakuum-Raum-Zeit in zwei Dimensionen und deren Anwendung zu Problemen des kritischen Kollaps N2 - This Thesis puts its focus on the physics of neutron stars and its description with methods of numerical relativity. In the first step, a new numerical framework the Whisky2D code will be developed, which solves the relativistic equations of hydrodynamics in axisymmetry. Therefore we consider an improved formulation of the conserved form of these equations. The second part will use the new code to investigate the critical behaviour of two colliding neutron stars. Considering the analogy to phase transitions in statistical physics, we will investigate the evolution of the entropy of the neutron stars during the whole process. A better understanding of the evolution of thermodynamical quantities, like the entropy in critical process, should provide deeper understanding of thermodynamics in relativity. More specifically, we have written the Whisky2D code, which solves the general-relativistic hydrodynamics equations in a flux-conservative form and in cylindrical coordinates. This of course brings in 1/r singular terms, where r is the radial cylindrical coordinate, which must be dealt with appropriately. In the above-referenced works, the flux operator is expanded and the 1/r terms, not containing derivatives, are moved to the right-hand-side of the equation (the source term), so that the left hand side assumes a form identical to the one of the three-dimensional (3D) Cartesian formulation. We call this the standard formulation. Another possibility is not to split the flux operator and to redefine the conserved variables, via a multiplication by r. We call this the new formulation. The new equations are solved with the same methods as in the Cartesian case. From a mathematical point of view, one would not expect differences between the two ways of writing the differential operator, but, of course, a difference is present at the numerical level. Our tests show that the new formulation yields results with a global truncation error which is one or more orders of magnitude smaller than those of alternative and commonly used formulations. The second part of the Thesis uses the new code for investigations of critical phenomena in general relativity. In particular, we consider the head-on-collision of two neutron stars in a region of the parameter space where two final states a new stable neutron star or a black hole, lay close to each other. In 1993, Choptuik considered one-parameter families of solutions, S[P], of the Einstein-Klein-Gordon equations for a massless scalar field in spherical symmetry, such that for every P > P⋆, S[P] contains a black hole and for every P < P⋆, S[P] is a solution not containing singularities. He studied numerically the behavior of S[P] as P → P⋆ and found that the critical solution, S[P⋆], is universal, in the sense that it is approached by all nearly-critical solutions regardless of the particular family of initial data considered. All these phenomena have the common property that, as P approaches P⋆, S[P] approaches a universal solution S[P⋆] and that all the physical quantities of S[P] depend only on |P − P⋆|. The first study of critical phenomena concerning the head-on collision of NSs was carried out by Jin and Suen in 2007. In particular, they considered a series of families of equal-mass NSs, modeled with an ideal-gas EOS, boosted towards each other and varied the mass of the stars, their separation, velocity and the polytropic index in the EOS. In this way they could observe a critical phenomenon of type I near the threshold of black-hole formation, with the putative solution being a nonlinearly oscillating star. In a successive work, they performed similar simulations but considering the head-on collision of Gaussian distributions of matter. Also in this case they found the appearance of type-I critical behaviour, but also performed a perturbative analysis of the initial distributions of matter and of the merged object. Because of the considerable difference found in the eigenfrequencies in the two cases, they concluded that the critical solution does not represent a system near equilibrium and in particular not a perturbed Tolmann-Oppenheimer-Volkoff (TOV) solution. In this Thesis we study the dynamics of the head-on collision of two equal-mass NSs using a setup which is as similar as possible to the one considered above. While we confirm that the merged object exhibits a type-I critical behaviour, we also argue against the conclusion that the critical solution cannot be described in terms of equilibrium solution. Indeed, we show that, in analogy with what is found in, the critical solution is effectively a perturbed unstable solution of the TOV equations. Our analysis also considers fine-structure of the scaling relation of type-I critical phenomena and we show that it exhibits oscillations in a similar way to the one studied in the context of scalar-field critical collapse. N2 - Diese Arbeit legt seinen Schwerpunkt auf die Physik von Neutronensternen und deren Beschreibung mit Methoden der numerischen Relativitätstheorie. Im ersten Schritt wird eine neue numerische Umgebung, der Whisky2D Code entwickelt, dieser löst die relativistischen Gleichungen der Hydrodynamik in Axialymmetrie. Hierzu betrachten wir eine verbesserte Formulierung der sog. "flux conserved formulation" der Gleichungen. Im zweiten Teil wird der neue Code verwendet , um das kritische Verhalten zweier kollidierenden Neutronensternen zu untersuchen. In Anbetracht der Analogie, um Übergänge in der statistischen Physik Phase werden wir die Entwicklung der Entropie der Neutronensterne während des gesamten Prozesses betrachten. Ein besseres Verständnis der Evolution von thermodynamischen Größen, wie der Entropie in kritischer Prozess, sollte zu einem tieferen Verständnis der relativistischen Thermodynamik führen. Der Whisky2D Code, zur Lösung Gleichungen relativistischer Hydrodynamik wurde in einer „flux conserved form“ und in zylindrischen Koordinaten geschrieben. Hierdurch entstehen 1 / r singuläre Terme, wobei r der ist, die entsprechend behandelt werden müssen. In früheren Arbeiten, wird der Operator expandiert und die 1 / r spezifisch Therme auf die rechte Seite geschrieben, so dass die linke Seite eine Form annimmt, die identisch ist mit der kartesischen Formulierung. Wir nennen dies die Standard-Formulierung. Eine andere Möglichkeit ist, die Terme nicht zu expandieren, den und den 1/r Term in die Gleichung hinein zu ziehen. Wir nennen dies die Neue-Formulierung. Die neuen Gleichungen werden mit den gleichen Verfahren wie im kartesischen Fall gelöst. Aus mathematischer Sicht ist keine Unterschiede zwischen den beiden Formulierungen zu erwarten, erst die numerische Sicht zeigt die Unterschiede auf. Versuche zeigen, dass die Neue-Formulierung numerische Fehler um mehrere Größenordnungen reduziert. Der zweite Teil der Dissertation verwendet den neuen Code für die Untersuchung kritischer Phänomene in der allgemeinen Relativitätstheorie. Insbesondere betrachten wir die Kopf-auf-Kollision zweier Neutronensterne in einem Bereich des Parameter Raums, deren zwei mögliche Endzustände entweder einen neuen stabilen Neutronenstern oder ein Schwarzes Loch darstellen. Im Jahr 1993, betrachtete Choptuik Ein-Parameter-Familien von Lösungen, S [P], der Einstein-Klein-Gordon-Gleichung für ein masseloses Skalarfeld in sphärischer Symmetrie, so dass für jedes P> P ⋆, S[P] ein Schwarzes Loch enthalten ist und jedes P