109 research outputs found
On the plane wave Riemann Problem in Fluid Dynamics
The paper contains a stability analysis of the plane-wave Riemann problem for
the two-dimensional hyperbolic conservation laws for an ideal compressible gas.
It is proved that the contact discontinuity in the plane-wave Riemann problem
is unstable under perturbations. The implications for Godunovs method are
discussed and it is shown that numerical post shock noise can set of a contact
instability. A relation to carbuncle instabilities is established.Comment: 27 pages, 1 figur
Phase appearance or disappearance in two-phase flows
This paper is devoted to the treatment of specific numerical problems which
appear when phase appearance or disappearance occurs in models of two-phase
flows. Such models have crucial importance in many industrial areas such as
nuclear power plant safety studies. In this paper, two outstanding problems are
identified: first, the loss of hyperbolicity of the system when a phase appears
or disappears and second, the lack of positivity of standard shock capturing
schemes such as the Roe scheme. After an asymptotic study of the model, this
paper proposes accurate and robust numerical methods adapted to the simulation
of phase appearance or disappearance. Polynomial solvers are developed to avoid
the use of eigenvectors which are needed in usual shock capturing schemes, and
a method based on an adaptive numerical diffusion is designed to treat the
positivity problems. An alternate method, based on the use of the hyperbolic
tangent function instead of a polynomial, is also considered. Numerical results
are presented which demonstrate the efficiency of the proposed solutions
Pion Interferometry for a Granular Source of Quark-Gluon Plasma Droplets
We examine the two-pion interferometry for a granular source of quark-gluon
plasma droplets. The evolution of the droplets is described by relativistic
hydrodynamics with an equation of state suggested by lattice gauge results.
Pions are assumed to be emitted thermally from the droplets at the freeze-out
configuration characterized by a freeze-out temperature . We find that the
HBT radius decreases if the initial size of the droplets decreases.
On the other hand, depends on the droplet spatial distribution and
is relatively independent of the droplet size. It increases with an increase in
the width of the spatial distribution and the collective-expansion velocity of
the droplets. As a result, the value of can lie close to
for a granular quark-gluon plasma source. The granular model of the emitting
source may provide an explanation to the RHIC HBT puzzle and may lead to a new
insight into the dynamics of the quark-gluon plasma phase transition.Comment: 5 pages, 4 figure
Correlated Gravitational Wave and Neutrino Signals from General-Relativistic Rapidly Rotating Iron Core Collapse
We present results from a new set of 3D general-relativistic hydrodynamic
simulations of rotating iron core collapse. We assume octant symmetry and focus
on axisymmetric collapse, bounce, the early postbounce evolution, and the
associated gravitational wave (GW) and neutrino signals. We employ a
finite-temperature nuclear equation of state, parameterized electron capture in
the collapse phase, and a multi-species neutrino leakage scheme after bounce.
The latter captures the important effects of deleptonization, neutrino cooling
and heating and enables approximate predictions for the neutrino luminosities
in the early evolution after core bounce. We consider 12-solar-mass and
40-solar-mass presupernova models and systematically study the effects of (i)
rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on
dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of
rapidly rotating core collapse is practically independent of progenitor mass
and precollapse structure. Moreover, we show that the effects of neutrino
leakage on the GW signal are strong only in nonrotating or slowly rotating
models in which GW emission is not dominated by inner core dynamics. In rapidly
rotating cores, core bounce of the centrifugally-deformed inner core excites
the fundamental quadrupole pulsation mode of the nascent protoneutron star. The
ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in
the GW signal and correlated strong variations in the rising luminosities of
antineutrino and heavy-lepton neutrinos. We find these features in cores that
collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational
energies sufficient to drive hyper-energetic core-collapse supernova
explosions. Hence, joint GW + neutrino observations of a core collapse event
could deliver strong evidence for or against rapid core rotation. [abridged]Comment: 29 pages, 14 figures. Replaced with version matching published
versio
The impact of viscosity on the morphology of gaseous flows in semidetached binary systems
Results of 3D gas dynamical simulation of mass transfer in binaries are
presented for systems with various values of viscosity. Analysis of obtained
solutions shows that in the systems with low value of viscosity the flow
structure is qualitatively similar to one for systems with high viscosity.
Presented calculations confirm that there is no shock interaction between the
stream from L1 and the forming accretion disk (`hot spot') at any value of
viscosity.Comment: LaTeX, 18 pages, 15 eps-figures, Astron. Reports, in pres
Scalar field induced oscillations of neutron stars and gravitational collapse
We study the interaction of massless scalar fields with self-gravitating
neutron stars by means of fully dynamic numerical simulations of the
Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to
spherical symmetry and the neutron stars are approximated by relativistic
polytropes. Studying the nonlinear dynamics of isolated neutron stars is very
effectively performed within the characteristic formulation of general
relativity, in which the spacetime is foliated by a family of outgoing light
cones. We are able to compactify the entire spacetime on a computational grid
and simultaneously impose natural radiative boundary conditions and extract
accurate radiative signals. We study the transfer of energy from the scalar
field to the fluid star. We find, in particular, that depending on the
compactness of the neutron star model, the scalar wave forces the neutron star
either to oscillate in its radial modes of pulsation or to undergo
gravitational collapse to a black hole on a dynamical timescale. The radiative
signal, read off at future null infinity, shows quasi-normal oscillations
before the setting of a late time power-law tail.Comment: 12 pages, 13 figures, submitted to Phys. Rev.
Simulation techniques for cosmological simulations
Modern cosmological observations allow us to study in great detail the
evolution and history of the large scale structure hierarchy. The fundamental
problem of accurate constraints on the cosmological parameters, within a given
cosmological model, requires precise modelling of the observed structure. In
this paper we briefly review the current most effective techniques of large
scale structure simulations, emphasising both their advantages and
shortcomings. Starting with basics of the direct N-body simulations appropriate
to modelling cold dark matter evolution, we then discuss the direct-sum
technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and
the tree algorithms. Simulations of baryonic matter in the Universe often use
hydrodynamic codes based on both particle methods that discretise mass, and
grid-based methods. We briefly describe Eulerian grid methods, and also some
variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.Comment: 42 pages, 16 figures, accepted for publication in Space Science
Reviews, special issue "Clusters of galaxies: beyond the thermal view",
Editor J.S. Kaastra, Chapter 12; work done by an international team at the
International Space Science Institute (ISSI), Bern, organised by J.S.
Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke
A New Open-Source Code for Spherically-Symmetric Stellar Collapse to Neutron Stars and Black Holes
We present the new open-source spherically-symmetric general-relativistic
(GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR
hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs
radial-gauge, polar-slicing coordinates in which the 3+1 equations simplify
substantially. We discretize the GRHD equations with a finite-volume scheme,
employing piecewise-parabolic reconstruction and an approximate Riemann solver.
GR1D is intended for the simulation of stellar collapse to neutron stars and
black holes and will also serve as a testbed for modeling technology to be
incorporated in multi-D GR codes. Its GRHD part is coupled to various
finite-temperature microphysical equations of state in tabulated form that we
make available with GR1D. An approximate deleptonization scheme for the
collapse phase and a neutrino-leakage/heating scheme for the postbounce epoch
are included and described. We also derive the equations for effective rotation
in 1D and implement them in GR1D. We present an array of standard test
calculations and also show how simple analytic equations of state in
combination with presupernova models from stellar evolutionary calculations can
be used to study qualitative aspects of black hole formation in failing
rotating core-collapse supernovae. In addition, we present a simulation with
microphysical EOS and neutrino leakage/heating of a failing core-collapse
supernova and black hole formation in a presupernova model of a 40 solar mass
zero-age main-sequence star. We find good agreement on the time of black hole
formation (within 20%) and last stable protoneutron star mass (within 10%) with
predictions from simulations with full Boltzmann neutrino radiation
hydrodynamics.Comment: 25 pages, 6 figures, 2 appendices. Accepted for publication to the
Classical and Quantum Gravity special issue for MICRA2009. Code may be
downloaded from http://www.stellarcollapse.org Update: corrected title, small
modifications suggested by the referees, added source term derivation in
appendix
Measuring the angular momentum distribution in core-collapse supernova progenitors with gravitational waves
- …