Source author record

Andreas Weyhausen

Andreas Weyhausen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

gr-qc physics.comp-ph

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2015arXiv

A Pseudospectral Method for Gravitational Wave Collapse

We present a new pseudospectral code, bamps, for numerical relativity written with the evolution of collapsing gravitational waves in mind. We employ the first order generalized harmonic gauge formulation. The relevant theory is reviewed and the numerical method is critically examined and specialized for the task at hand. In particular we investigate formulation parameters, gauge and constraint preserving boundary conditions well-suited to non-vanishing gauge source functions. Different types of axisymmetric twist-free moment of time symmetry gravitational wave initial data are discussed. A treatment of the axisymmetric apparent horizon condition is presented with careful attention to regularity on axis. Our apparent horizon finder is then evaluated in a number of test cases. Moving on to evolutions, we investigate modifications to the generalized harmonic gauge constraint damping scheme to improve conservation in the strong field regime. We demonstrate strong-scaling of our pseudospectral penalty code. We employ the Cartoon method to efficiently evolve axisymmetric data in our 3+1 dimensional code. We perform test evolutions of Schwarzschild perturbed by gravitational waves and by gauge pulses, both to demonstrate the use of our blackhole excision scheme and for comparison with earlier results. Finally numerical evolutions of supercritical Brill waves are presented to demonstrate durability of the excision scheme for the dynamical formation of a blackhole.

preprint2015arXiv

Solving 3D relativistic hydrodynamical problems with WENO discontinuous Galerkin methods

Discontinuous Galerkin (DG) methods coupled to WENO algorithms allow high order convergence for smooth problems and for the simulation of discontinuities and shocks. In this work, we investigate WENO-DG algorithms in the context of numerical general relativity, in particular for general relativistic hydrodynamics. We implement the standard WENO method at different orders, a compact (simple) WENO scheme, as well as an alternative subcell evolution algorithm. To evaluate the performance of the different numerical schemes, we study non-relativistic, special relativistic, and general relativistic testbeds. We present the first three-dimensional simulations of general relativistic hydrodynamics, albeit for a fixed spacetime background, within the framework of WENO-DG methods. The most important testbed is a single TOV-star in three dimensions, showing that long term stable simulations of single isolated neutron stars can be obtained with WENO-DG methods.

preprint2013arXiv

Collapse of Nonlinear Gravitational Waves in Moving-Puncture Coordinates

We study numerical evolutions of nonlinear gravitational waves in moving-puncture coordinates. We adopt two different types of initial data -- Brill and Teukolsky waves -- and evolve them with two independent codes producing consistent results. We find that Brill data fail to produce long-term evolutions for common choices of coordinates and parameters, unless the initial amplitude is small, while Teukolsky wave initial data lead to stable evolutions, at least for amplitudes sufficiently far from criticality. The critical amplitude separates initial data whose evolutions leave behind flat space from those that lead to a black hole. For the latter we follow the interaction of the wave, the formation of a horizon, and the settling down into a time-independent trumpet geometry. We explore the differences between Brill and Teukolsky data and show that for less common choices of the parameters -- in particular negative amplitudes -- Brill data can be evolved with moving-puncture coordinates, and behave similarly to Teukolsky waves.

preprint2012arXiv

Constraint damping for the Z4c formulation of general relativity

One possibility for avoiding constraint violation in numerical relativity simulations adopting free-evolution schemes is to modify the continuum evolution equations so that constraint violations are damped away. Gundlach et. al. demonstrated that such a scheme damps low amplitude, high frequency constraint violating modes exponentially for the Z4 formulation of General Relativity. Here we analyze the effect of the damping scheme in numerical applications on a conformal decomposition of Z4. After reproducing the theoretically predicted damping rates of constraint violations in the linear regime, we explore numerical solutions not covered by the theoretical analysis. In particular we examine the effect of the damping scheme on low-frequency and on high-amplitude perturbations of flat spacetime as well and on the long-term dynamics of puncture and compact star initial data in the context of spherical symmetry. We find that the damping scheme is effective provided that the constraint violation is resolved on the numerical grid. On grid noise the combination of artificial dissipation and damping helps to suppress constraint violations. We find that care must be taken in choosing the damping parameter in simulations of puncture black holes. Otherwise the damping scheme can cause undesirable growth of the constraints, and even qualitatively incorrect evolutions. In the numerical evolution of a compact static star we find that the choice of the damping parameter is even more delicate, but may lead to a small decrease of constraint violation. For a large range of values it results in unphysical behavior.