Researcher profile

Wolfgang Kastaun

Wolfgang Kastaun contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2022arXiv

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) - black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on `conserved' variables which then need to be converted back into the fundamental (`primitive') variables. The corresponding conservative-to-primitive variable recovery procedure, based on root-finding algorithms, constitutes one of the core elements of such GRMHD codes. Recently, a new robust, accurate and efficient recovery scheme called RePrimAnd was introduced, which has demonstrated the ability to always converge to a unique solution. The scheme provides fine-grained error policies to handle invalid states caused by evolution errors, and also provides analytical bounds for the error of all primitive variables. In this work, we describe the technical aspects of implementing the RePrimAnd scheme into the GRMHD code Spritz. To check our implementation as well as to assess the various features of the scheme, we perform a number of GRMHD tests in three dimensions. Our tests, which include critical cases such as a NS collapse to a BH as well as the early evolution (~50 ms) of a Fishbone-Moncrief BH-accrection disk system, show that RePrimAnd is able to support magnetized, low density environments with magnetic-to-fluid pressure ratios as high as 10^4, in situations where the previously used recovery scheme fails.

preprint2021arXiv

Robust Recovery of Primitive Variables in Relativistic Ideal Magnetohydrodynamics

Modern simulation codes for general relativistic ideal magnetohydrodynamics are all facing a long standing technical problem given by the need to recover fundamental variables from those variables that are evolved in time. In the relativistic case, this requires the numerical solution of a system of nonlinear equations. Although several approaches are available, none has proven completely reliable. A recent study comparing different methods showed that all can fail, in particular for the important case of strong magnetization and moderate Lorentz factors. Here, we propose a new robust, efficient, and accurate solution scheme, along with a proof for the existence and uniqueness of a solution, and analytic bounds for the accuracy. Further, the scheme allows us to reliably detect evolution errors leading to unphysical states and automatically applies corrections for typical harmless cases. A reference implementation of the method is made publicly available as a software library. The aim of this library is to improve the reliability of binary neutron star merger simulations, in particular in the investigation of jet formation and magnetically driven winds.

preprint2013arXiv

Implementation of a simplified approach to radiative transfer in general relativity

We describe in detail the implementation of a simplified approach to radiative transfer in general relativity by means of the well-known neutrino leakage scheme (NLS). In particular, we carry out an extensive investigation of the properties and limitations of the NLS for isolated relativistic stars to a level of detail that has not been discussed before in a general-relativistic context. Although the numerous tests considered here are rather idealized, they provide a well-controlled environment in which to understand the relationship between the matter dynamics and the neutrino emission, which is important in order to model the neutrino signals from more complicated scenarios, such as binary neutron-star mergers. When considering nonrotating hot neutron stars we confirm earlier results of one-dimensional simulations, but also present novel results about the equilibrium properties and on how the cooling affects the stability of these configurations. In our idealized but controlled setup, we can then show that deviations from the thermal and weak-interaction equilibrium affect the stability of these models to radial perturbations, leading models that are stable in the absence of radiative losses, to a gravitational collapse to a black hole when neutrinos are instead radiated.

preprint2013arXiv

On the black hole from merging binary neutron stars: how fast can it spin?

The merger of two neutron stars will in general lead to the formation of a torus surrounding a black hole whose rotational energy can be tapped to potentially power a short gamma-ray burst. We have studied the merger of equal-mass binaries with spins aligned with the orbital angular momentum to determine the maximum spin the black hole can reach. Our initial data consists of irrotational binaries to which we add various amounts of rotation to increase the total angular momentum. Although the initial data violates the constraint equations, the use of the constraint-damping CCZ4 formulation yields evolutions with violations smaller than those with irrotational initial data and standard formulations. Interestingly, we find that a limit of $J/M^2 \simeq 0.89$ exists for the dimensionless spin and that any additional angular momentum given to the binary ends up in the torus rather than in the black hole, thus providing another nontrivial example supporting the cosmic censorship hypothesis.

preprint2011arXiv

On the saturation amplitude of the f-mode instability

We investigate strong nonlinear damping effects which occur during high amplitude oscillations of neutron stars, and the gravitational waves they produce. For this, we use a general relativistic nonlinear hydrodynamics code in conjunction with a fixed spacetime (Cowling approximation) and a polytropic equation of state (EOS). Gravitational waves are estimated using the quadrupole formula. Our main interest are l=m=2 f modes subject to the CFS (Chandrasekhar, Friedman, Schutz) instability, but we also investigate axisymmetric and quasiradial modes. We study various models to determine the influence of rotation rate and EOS. We find that axisymmetric oscillations at high amplitudes are predominantly damped by shock formation, while the nonaxisymmetric f modes are mainly damped by wave breaking and, for rapidly rotating models, coupling to nonaxisymmetric inertial modes. From the observed nonlinear damping, we derive upper limits for the saturation amplitude of CFS-unstable f modes. Finally, we estimate that the corresponding gravitational waves for an oscillation amplitude at the upper limit should be detectable with the advanced LIGO and VIRGO interferometers at distances above 10 MPc. This strongly depends on the stellar model, in particular on the mode frequency.