BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ewald Mueller

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEgr-qcastro-phastro-ph.COphysics.flu-dyn

Catalog footprint

What is connected

17works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

17 published item(s)

preprint2020arXiv

Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations

The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. 1D overshooting descriptions introduce a free parameter ($f_{ov}$) that is in general not well constrained from observations. Especially in small central convective regions the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional hydrodynamic simulations can be used to study the size of the overshooting region and the involved mixing processes. Here we show how one can calibrate an overshooting parameter by performing 2D Maestro simulations of Zero-Age-Main-Sequence stars ranging from $1.3$ to $3.5 M_\odot$. The simulations cover the convective cores of the stars and a large fraction of the surrounding radiative envelope. We follow the convective flow for at least 20 convective turnover times, while the longest simulation covers 430 turnover time scales. This allows us to study how the mixing as well as the convective boundary evolve with time, and how the resulting entrainment can be interpreted in terms of overshooting parameters. We find that increasing the overshooting parameter $f_{ov}$ beyond a certain value in the initial model of our simulations, changes the mixing behaviour completely. This result can be used to put limits on the overshooting parameter. We find $0.010 < f_{ov} < 0.017$ to be in good agreement with our simulations of a $3.5 M_\odot$ mass star. We also identify a diffusive mixing component due to internal gravity waves (IGW) that is active throughout the convectively stable layer, but likely overestimated in our simulations. Furthermore, applying our calibration method to simulations of less massive stars suggests a need for a mass-dependent overshooting description where the mixing in terms of the pressure scale height is reduced for small convective cores.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Magnetic Field Amplification During the Common Envelope Phase

During the common envelope (CE) phase, a giant star in a binary system overflows its Roche lobe and unstable mass transfer leads to a spiral-in of the companion, resulting in a close binary system or in a merger of the stellar cores. Dynamo processes during the CE phase have been proposed as a mechanism to generate magnetic fields that are important for forming magnetic white dwarfs (MWDs) and for shaping planetary nebulae. Here, we present the first magnetohydrodynamics simulations of the dynamical spiral-in during a CE phase. We find that magnetic fields are strongly amplified in the accretion stream around the $1M_\odot$ companion as it spirals into the envelope of a $2M_\odot$ RG. This leads to field strengths of 10 to 100 kG throughout the envelope after 120 d. The magnetic field amplification is consistent with being driven by the magnetorotational instability. The field strengths reached in our simulation make the magnetic field interesting for diagnostic purposes, but they are dynamically irrelevant. They are also too small to explain the formation of the highest fields found in MWDs, but may be relevant for luminous red novae, and detecting magnetic fields in these events would support the scenario as proposed here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutrino viscosity and drag: impact on the magnetorotational instability in protoneutron stars

The magnetorotational instability (MRI) is a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. The diffusion of neutrinos trapped in the PNS induces a transport of momentum, which can be modelled as a viscosity on length-scales longer than the neutrino mean free path. This neutrino-viscosity can slow down the growth of MRI modes to such an extent that a minimum initial magnetic field strength of $\gtrsim 10^{12}\, {\rm G}$ is needed for the MRI to grow on a sufficiently short time-scale to potentially affect the explosion. It is uncertain whether the magnetic field of fast rotating progenitor cores is strong enough to yield such an initial magnetic field in PNS. At MRI wavelengths shorter than the neutrino mean free path, on the other hand, neutrino radiation does not act as a viscosity but rather induces a drag on the velocity with a damping rate independent of the wavelength. We perform a linear analysis of the MRI in this regime, and apply our analytical results to the PNS structure from a one-dimensional numerical simulation. We show that in the outer layers of the PNS, the MRI can grow from weak magnetic fields at wavelengths shorter than the neutrino mean free path, while deeper in the PNS MRI growth takes place in the viscous regime and requires a minimum magnetic field strength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical simulations of the magnetorotational instability in protoneutron stars: I. Influence of buoyancy

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Supernova 1987A: neutrino-driven explosions in three dimensions and light curves

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-Dimensional Simulations of Core-Collapse Supernovae: From Shock Revival to Shock Breakout

We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: comparison with the BSSN formulation in spherical symmetry

We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal", or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preserving outer boundary conditions, nor does it need any modifications of the equations for evolutions of black holes. We perform several tests and compare the performance of the fCCZ4 system, for different choices of certain free parameters, with that of BSSN. Confirming earlier results we find that, for an optimal choice of these parameters, and for neutron-star spacetimes, the violations of the Hamiltonian constraint can be between 1 and 3 orders of magnitude smaller in the fCCZ4 system than in the BSSN formulation. For black-hole spacetimes, on the other hand, any advantages of fCCZ4 over BSSN are less evident.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Relativistic collapse and explosion of rotating supermassive stars with thermonuclear effects

We present results of general relativistic simulations of collapsing supermassive stars with and without rotation using the two-dimensional general relativistic numerical code Nada, which solves the Einstein equations written in the BSSN formalism and the general relativistic hydrodynamics equations with high resolution shock capturing schemes. These numerical simulations use an equation of state which includes effects of gas pressure, and in a tabulated form those associated with radiation and the electron-positron pairs. We also take into account the effect of thermonuclear energy released by hydrogen and helium burning. We find that objects with a mass of 5x10^{5} solar mass and an initial metallicity greater than Z_{CNO}~0.007 do explode if non-rotating, while the threshold metallicity for an explosion is reduced to Z_{CNO}~0.001 for objects uniformly rotating. The critical initial metallicity for a thermonuclear explosion increases for stars with mass ~10^{6} solar mass. For those stars that do not explode we follow the evolution beyond the phase of black hole formation. We compute the neutrino energy loss rates due to several processes that may be relevant during the gravitational collapse of these objects. The peak luminosities of neutrinos and antineutrinos of all flavors for models collapsing to a BH are ~10^{55} erg/s. The total radiated energy in neutrinos varies between ~10^{56} ergs for models collapsing to a BH, and ~10^{45}-10^{46} ergs for models exploding.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Core-Collapse Supernovae: Explosion Dynamics, Neutrinos and Gravitational Waves

The quest for the supernova explosion mechanism has been one of the outstanding challenges in computational astrophysics for several decades. Simulations have now progressed to a stage at which the solution appears close and neutrino and gravitational wave signals from self-consistent explosion models are becoming available. Here we focus one of the recent advances in supernova modeling, the inclusion of general relativity in multi-dimensional neutrino hydrodynamics simulations, and present the latest simulation results for an 11.2 and a 15 solar mass progenitor. We also mention 3D effects as another aspect in supernova physics awaiting further, more thorough investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rayleigh-Taylor finger instability mixing in hydrodynamic shell convection models

Mixing processes in stars driven by composition gradients as a result of the Rayleigh-Taylor instability are not anticipated. They are supported only by hydrodynamic studies of stellar convection. We find that such mixing occurs below the bottom edge of convection zones in our multidimensional hydrodynamic shell convection models. It operates at interfaces created by off-center nuclear burning, where less dense gas with higher mean molecular weight is located above denser gas with a lower mean molecular weight. We discuss the mixing under various conditions with hydrodynamic convection models based on stellar evolutionary calculations of the core helium flash in a 1.25 Msun star, the core carbon flash in a 9.3 Msun star, and of oxygen burning shell in a star with a mass of 23 Msun. We simulate the hydrodynamic behavior of shell convection during various phases of stellar evolution with the Eulerian hydrodynamics code HERAKLES in two and three spatial dimensions. Initial models for this purpose are obtained by state-of-the-art stellar evolutionary codes, namely GARSTEC, STAREVOL, and TYCHO for the core helium flash, core carbon flash, and oxygen shell burning, respectively. Most of our analysis is performed for two-dimensional hydrodynamic models of shell convection during the core helium flash at its peak covering approximately 250 convective turnover timescales. The mixing manifests itself in the form of overdense and cold fingers enriched with matter of higher mean molecular weight, originating from density fluctuations at the lower boundary of the convective shell, and "shooting" down into the core. They result from the Rayleigh-Taylor instability at the lower convection zone boundary due to a negative mean molecular weight gradient. They do not appear when the mean molecular weight gradient is positive.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The core helium flash revisited III. From Pop I to Pop III stars

Degenerate ignition of helium in low-mass stars at the end of the red giant branch phase leads to dynamic convection in their helium cores. One-dimensional (1D) stellar modeling of this intrinsically multi-dimensional dynamic event is likely to be inadequate. Previous hydrodynamic simulations imply that the single convection zone in the helium core of metal-rich Pop I stars grows during the flash on a dynamic timescale. This may lead to hydrogen injection into the core, and a double convection zone structure as known from one-dimensional core helium flash simulations of low-mass Pop III stars. We perform hydrodynamic simulations of the core helium flash in two and three dimensions to better constrain the nature of these events. To this end we study the hydrodynamics of convection within the helium cores of a 1.25 \Msun metal-rich Pop I star (Z=0.02), and a 0.85 \Msun metal-free Pop III star (Z=0) near the peak of the flash. These models possess single and double convection zones, respectively. We use 1D stellar models of the core helium flash computed with state-of-the-art stellar evolution codes as initial models for our multidimensional hydrodynamic study, and simulate the evolution of these models with the Riemann solver based hydrodynamics code Herakles which integrates the Euler equations coupled with source terms corresponding to gravity and nuclear burning. The hydrodynamic simulation of the Pop I model involving a single convection zone covers 27 hours of stellar evolution, while the first hydrodynamic simulations of a double convection zone, in the Pop III model, span 1.8 hours of stellar life. We find differences between the predictions of mixing length theory and our hydrodynamic simulations. The simulation of the single convection zone in the Pop I model shows a strong growth of the size of the convection zone due to turbulent entrainment. Hence we predict that for the Pop I model a hydrogen injection phase (i.e. hydrogen injection into the helium core) will commence after about 23 days, which should eventually lead to a double convection zone structure known from 1D stellar modeling of low-mass Pop III stars. Our two and three-dimensional hydrodynamic simulations of the double (Pop III) convection zone model show that the velocity field in the convection zones is different from that predicted by stellar evolutionary calculations. The simulations suggest that the double convection zone decays quickly, the flow eventually being dominated by internal gravity waves.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Fragmentation of general relativistic quasi-toroidal polytropes

We investigate the role of rotational instabilities in the context of black hole formation in relativistic stars. In addition to the standard scenario - an axially symmetric dynamical instability forming a horizon at the star's center - the recently found low-$T/|W|$ instabilities are shown to lead to fragmentation and off-center horizon formation in differentially rotating stars. This process might be an alternative pathway to produce SMBHs from supermassive stars with inefficient angular momentum transport.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Black hole formation through fragmentation of toroidal polytropes

We investigate new paths to black hole formation by considering the general relativistic evolution of a differentially rotating polytrope with toroidal shape. We find that this polytrope is unstable to nonaxisymmetric modes, which leads to a fragmentation into self-gravitating, collapsing components. In the case of one such fragment, we apply a simplified adaptive mesh refinement technique to follow the evolution to the formation of an apparent horizon centered on the fragment. This is the first study of the one-armed instability in full general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Numerical hydrodynamics in special relativity

This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discussed. An evaluation of the various numerical methods is given and future developments are analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The First Second of a Type-II Supernova: Convection, Accretion, and Shock Propagation

One- and two-dimensional hydrodynamical simulations of neutrino-driven supernova explosions are discussed. The simulations cover the phase between the stagnation of the prompt shock and about one second after core bounce. Systematic variation of the neutrino fluxes from the neutrino sphere shows that the explosion energy, explosion time scale, initial mass of the protoneutron star, and explosive nucleosynthesis of iron-group elements depend sensitively on the strength of the neutrino heating during the first few 100 ms after shock formation. Convective overturn in the neutrino-heated region behind the shock is a crucial help for the explosion only in a narrow window of neutrino luminosities. Here powerful explosions can be obtained only in the multi-dimensional case. For higher core-neutrino fluxes also spherically symmetrical models yield energetic explosions, while for lower luminosities even with convection no strong explosions occur.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

The Formation of Disk Galaxies in a Cosmological Context: Structure and Kinematics

We present results concerning the internal structure and kinematics of disk galaxies formed in cosmologically motivated simulations. The calculations include dark matter, gas dynamics, radiative cooling, star formation, supernova feedback and metal enrichment. The initial model is a rigidly rotating overdense sphere with a mass of about 8 10^11 Msol which is perturbed by small scale fluctuations according to a biased CDM power spectrum. Converging, Jeans unstable and rapidly cooling regions are allowed to form stars. Via supernovae, metal enriched gas is returned to the interstellar medium. {}From these initial conditions a galaxy forms which shows the main properties of spiral galaxies: a rotationally supported exponential disk which consists of young stars with about solar metallicity, a slowly rotating halo of old metal poor stars, a bulge of old metal rich stars and a slowly rotating extended halo of dark matter. Bulge, stellar and dark halo are supported by an anisotropic velocity dispersion and have a de Vaucouleurs surface density profile. The flattening of the dark and stellar halo is too large to be explained by rotation only. Whether the flattening of the bulge is caused by an anisotropic velocity dispersion or by its rotation cannot be answered, because of the limited numerical resolution due to gravitational softening. The velocity dispersion and the thickness of the stellar disk increase with the age of the stars. Considering only the young stellar component, the disk is cold (sigma=20 km/sec) and thin (z <1 kpc). The dynamical formation process ends after about 4\,Gyr, when

0 citations0 reviewsNo code linkedOpen access