134 research outputs found
SPH Simulations with Reconfigurable Hardware Accelerator
We present a novel approach to accelerate astrophysical hydrodynamical
simulations. In astrophysical many-body simulations, GRAPE (GRAvity piPE)
system has been widely used by many researchers. However, in the GRAPE systems,
its function is completely fixed because specially developed LSI is used as a
computing engine. Instead of using such LSI, we are developing a special
purpose computing system using Field Programmable Gate Array (FPGA) chips as
the computing engine. Together with our developed programming system, we have
implemented computing pipelines for the Smoothed Particle Hydrodynamics (SPH)
method on our PROGRAPE-3 system. The SPH pipelines running on PROGRAPE-3 system
have the peak speed of 85 GFLOPS and in a realistic setup, the SPH calculation
using one PROGRAPE-3 board is 5-10 times faster than the calculation on the
host computer. Our results clearly shows for the first time that we can
accelerate the speed of the SPH simulations of a simple astrophysical phenomena
using considerable computing power offered by the hardware.Comment: 27 pages, 13 figures, submitted to PAS
Evolution of Collisionally Merged Massive Stars
We investigate the evolution of collisionally merged stars with mass of ~100 Msun which might be formed in dense star clusters. We assumed that massive stars with several tens Msun collide typically after ~1Myr of the formation of the cluster and performed hydrodynamical simulations of several collision events. Our simulations show that after the collisions, merged stars have extended envelopes and their radii are larger than those in the thermal equilibrium states and that their interiors are He-rich because of the stellar evolution of the progenitor stars. We also found that if the mass-ratio of merging stars is far from unity, the interior of the merger product is not well mixed and the elemental abundance is not homogeneous. We then followed the evolution of these collision products by a one dimensional stellar evolution code. After an initial contraction on the Kelvin-Helmholtz (thermal adjustment) timescale (~10^{3-4} yr), the evolution of the merged stars traces that of single homogeneous stars with corresponding masses and abundances, while the initial contraction phase shows variations which depend on the mass ratio of the merged stars. We infer that, once runaway collisions have set in, subsequent collisions of the merged stars take place before mass loss by stellar winds becomes significant. Hence, stellar mass loss does not inhibit the formation of massive stars with mass of ~1000Msun
Numerical Simulations of Globular Cluster Formation
We examine various physical processes associated with the formation of
globular clusters by using the three-dimensional Smoothed Particle
Hydrodynamics (SPH) code. Our code includes radiative cooling of gases, star
formation, energy feedback from stars including stellar winds and supernovae,
and chemical enrichment by stars. We assume that, in the collapsing galaxy,
isothermal cold clouds form through thermal condensations and become
proto-globular clouds. We calculate the size of proto-globular clouds by
solving the linearized equations for perturbation. We compute the evolution of
the inner region of the proto-cloud with our SPH code for various initial
radius and initial composition of gases. When the initial gases contain no
heavy elements, the evolution of proto-clouds sensitively depends on the
initial radius. For a smaller initial radius, the initial star burst is so
intense that the subsequent star formation occurs in the central regions to
form a dense star cluster as massive as the globular cluster. When the initial
gases contain some heavy elements, the metallicity of gases affects the
evolution and the final stellar mass. If the initial radius of the
proto-globular clouds was relatively large, the formation of a star cluster as
massive as the globular clusters requires the initial metallicity as high as
[Fe/H] . The self-enrichment of heavy elements in the star cluster
does not occur in all cases.Comment: Accpeted for publication in the ApJ. Correctiong errors in Table
Nucleosynthesis in Type II Supernovae
Presupernova evolution and explosive nucleosynthesis in massive stars for
main-sequence masses from 13 to 70 are calculated. We
examine the dependence of the supernova yields on the stellar mass,
^{12}C(\alpha, \gamma) ^{16}O} rate, and explosion energy. The supernova
yields integrated over the initial mass function are compared with the solar
abundances.Comment: 1 Page Latex source, 10 PostScript figures, to appear in Nuclear
Physics A, Vol. A616 (1997
PyCOOL - a Cosmological Object-Oriented Lattice code written in Python
There are a number of different phenomena in the early universe that have to
be studied numerically with lattice simulations. This paper presents a graphics
processing unit (GPU) accelerated Python program called PyCOOL that solves the
evolution of scalar fields in a lattice with very precise symplectic
integrators. The program has been written with the intention to hit a sweet
spot of speed, accuracy and user friendliness. This has been achieved by using
the Python language with the PyCUDA interface to make a program that is easy to
adapt to different scalar field models. In this paper we derive the symplectic
dynamics that govern the evolution of the system and then present the
implementation of the program in Python and PyCUDA. The functionality of the
program is tested in a chaotic inflation preheating model, a single field
oscillon case and in a supersymmetric curvaton model which leads to Q-ball
production. We have also compared the performance of a consumer graphics card
to a professional Tesla compute card in these simulations. We find that the
program is not only accurate but also very fast. To further increase the
usefulness of the program we have equipped it with numerous post-processing
functions that provide useful information about the cosmological model. These
include various spectra and statistics of the fields. The program can be
additionally used to calculate the generated curvature perturbation. The
program is publicly available under GNU General Public License at
https://github.com/jtksai/PyCOOL . Some additional information can be found
from http://www.physics.utu.fi/tiedostot/theory/particlecosmology/pycool/ .Comment: 23 pages, 12 figures; some typos correcte
Nucleosynthesis in Type Ia Supernovae
Among the major uncertainties involved in the Chandrasekhar mass models for
Type Ia supernovae are the companion star of the accreting white dwarf (or the
accretion rate that determines the carbon ignition density) and the flame speed
after ignition. We present nucleosynthesis results from relatively slow
deflagration (1.5 - 3 % of the sound speed) to constrain the rate of accretion
from the companion star. Because of electron capture, a significant amount of
neutron-rich species such as ^{54}Cr, ^{50}Ti, ^{58}Fe, ^{62}Ni, etc. are
synthesized in the central region. To avoid the too large ratios of
^{54}Cr/^{56}Fe and ^{50}Ti/^{56}Fe, the central density of the white dwarf at
thermonuclear runaway must be as low as \ltsim 2 \e9 \gmc. Such a low central
density can be realized by the accretion as fast as \dot M \gtsim 1 \times
10^{-7} M_\odot yr^{-1}. These rapidly accreting white dwarfs might correspond
to the super-soft X-ray sources.Comment: 10 page LaTeX, 7 PostScript figures, to appear in Nuclear Physics A,
Vol. A621 (1997
- …