433 research outputs found
Smooth-Particle Phase Stability with density and density-gradient potentials
Stable fluid and solid particle phases are essential to the simulation of
continuum fluids and solids using Smooth Particle Applied Mechanics. We show
that density-dependent potentials, such as Phi=(1/2)Sum (rho-rho_0)^2, along
with their corresponding constitutive relations, provide a simple means for
characterizing fluids and that a special stabilization potential, Phi=(1/2)Sum
(delrho)^2, not only stabilizes crystalline solid phases (or meshes) but also
provides a surface tension which is missing in the usual
density-dependent-potential approach. We illustrate these ideas for
two-dimensional square, triangular, and hexagonal lattices.Comment: 10 pages, 5 figure
Relativistic Hydrodynamic Flows Using Spatial and Temporal Adaptive Structured Mesh Refinement
Astrophysical relativistic flow problems require high resolution
three-dimensional numerical simulations. In this paper, we describe a new
parallel three-dimensional code for simulations of special relativistic
hydrodynamics (SRHD) using both spatially and temporally structured adaptive
mesh refinement (AMR). We used the method of lines to discretize the SRHD
equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme
for time integration. For spatial reconstruction, we have implemented piecewise
linear method (PLM), piecewise parabolic method (PPM), third order convex
essentially non-oscillatory (CENO) and third and fifth order weighted
essentially non-oscillatory (WENO) schemes. Flux is computed using either
direct flux reconstruction or approximate Riemann solvers including HLL,
modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR
part of the code is built on top of the cosmological Eulerian AMR code {\sl
enzo}. We discuss the coupling of the AMR framework with the relativistic
solvers. Via various test problems, we emphasize the importance of resolution
studies in relativistic flow simulations because extremely high resolution is
required especially when shear flows are present in the problem. We also
present the results of two 3d simulations of astrophysical jets: AGN jets and
GRB jets. Resolution study of those two cases further highlights the need of
high resolutions to calculate accurately relativistic flow problems.Comment: 14 pages, 23 figures. A section on 3D GRB jet simulation added.
Accepted by ApJ
Tidal exposure or microhabitats: what determines sandy-beach nematode zonation? A case study of a macrotidal ridge-and-runnel sandy beach in Belgium
Lately, across-shore zonation has been found to be more important in structuring the nematode community of a tropical macrotidal sandy beach than microhabitat heterogeneity. To evaluate whether this zonation pattern applies to a temperate beach, a macrotidal ridge-and-runnels sandy beach in the North Sea was studied. We investigated whether a similar zonation occurs in sandbar and runnel microhabitats, and whether the runnels harbour a different community from the subtidal. Our results indicate that nematode communities from runnel and sandbar habitats are significantly different. In addition, horizontal zonation patterns for nematode communities differ between both habitats. Nematode assemblages from sandbars are divided to lower, middle and upper beach while upper and middle runnels cluster together. The subtidal and upper runnels showed dissimilar nematode assemblages, although runnels showed the same dominant species (Daptonema normandicum), which increases its abundance towards the upper runnels. This study illustrates the importance of microhabitat heterogeneity, which resulted in different zonation patterns across the sandy beach examined. The divergent zonation between sandbars and runnels in the macrotidal temperate sandy beach, compared with the pattern observed for a subtropical sandy beach with similar morphodynamics, indicates that generalizations about nematode distribution patterns should be made with caution
Solving One Dimensional Scalar Conservation Laws by Particle Management
We present a meshfree numerical solver for scalar conservation laws in one
space dimension. Points representing the solution are moved according to their
characteristic velocities. Particle interaction is resolved by purely local
particle management. Since no global remeshing is required, shocks stay sharp
and propagate at the correct speed, while rarefaction waves are created where
appropriate. The method is TVD, entropy decreasing, exactly conservative, and
has no numerical dissipation. Difficulties involving transonic points do not
occur, however inflection points of the flux function pose a slight challenge,
which can be overcome by a special treatment. Away from shocks the method is
second order accurate, while shocks are resolved with first order accuracy. A
postprocessing step can recover the second order accuracy. The method is
compared to CLAWPACK in test cases and is found to yield an increase in
accuracy for comparable resolutions.Comment: 15 pages, 6 figures. Submitted to proceedings of the Fourth
International Workshop Meshfree Methods for Partial Differential Equation
The Chemical Compositions of the Type II Cepheids -- The BL Her and W Vir Variables
Abundance analyses from high-resolution optical spectra are presented for 19
Type II Cepheids in the Galactic field. The sample includes both short-period
(BL Her) and long-period (W Vir) stars. This is the first extensive abundance
analysis of these variables. The C, N, and O abundances with similar spreads
for the BL Her and W Vir show evidence for an atmosphere contaminated with
-process and CN-cycling products. A notable anomaly of the BL Her
stars is an overabundance of Na by a factor of about five relative to their
presumed initial abundances. This overabundance is not seen in the W Vir stars.
The abundance anomalies running from mild to extreme in W Vir stars but not
seen in the BL Her stars are attributed to dust-gas separation that provides an
atmosphere deficient in elements of high condensation temperature, notably Al,
Ca, Sc, Ti, and -process elements. Such anomalies have previously been seen
among RV Tau stars which represent a long-period extension of the variability
enjoyed by the Type II Cepheids. Comments are offered on how the contrasting
abundance anomalies of BL Her and W Vir stars may be explained in terms of the
stars' evolution from the blue horizontal branch.Comment: 41 pages including 11 figures and 4 tables; Accepted for publication
in Ap
Effects of a Supermassive Black Hole Binary on a Nuclear Gas Disk
We study influence of a galactic central supermassive black hole (SMBH)
binary on gas dynamics and star formation activity in a nuclear gas disk by
making three-dimensional Tree+SPH simulations. Due to orbital motions of SMBHs,
there are various resonances between gas motion and the SMBH binary motion. We
have shown that these resonances create some characteristic structures of gas
in the nuclear gas disk, for examples, gas elongated or filament structures,
formation of gaseous spiral arms, and small gas disks around SMBHs. In these
gaseous dense regions, active star formations are induced. As the result, many
star burst regions are formed in the nuclear region.Comment: 19 pages, 11 figures, accepted for publication in Ap
SPH Simulations of Counterrotating Disk Formation in Spiral Galaxies
We present the results of Smoothed Particle Hydrodynamics (SPH) simulations
of the formation of a massive counterrotating disk in a spiral galaxy. The
current study revisits and extends (with SPH) previous work carried out with
sticky particle gas dynamics, in which adiabatic gas infall and a retrograde
gas-rich dwarf merger were tested as the two most likely processes for
producing such a counterrotating disk. We report on experiments with a cold
primary similar to our Galaxy, as well as a hot, compact primary modeled after
NGC 4138. We have also conducted numerical experiments with varying amounts of
prograde gas in the primary disk, and an alternative infall model (a spherical
shell with retrograde angular momentum). The structure of the resulting
counterrotating disks is dramatically different with SPH. The disks we produce
are considerably thinner than the primary disks and those produced with sticky
particles. The time-scales for counterrotating disk formation are shorter with
SPH because the gas loses kinetic energy and angular momentum more rapidly.
Spiral structure is evident in most of the disks, but an exponential radial
profile is not a natural byproduct of these processes. The infalling gas shells
that we tested produce counterrotating bulges and rings rather than disks. The
presence of a considerable amount of preexisting prograde gas in the primary
causes, at least in the absence of star formation, a rapid inflow of gas to the
center and a subsequent hole in the counterrotating disk. In general, our SPH
experiments yield stronger evidence to suggest that the accretion of massive
counterrotating disks drives the evolution of the host galaxies towards earlier
(S0/Sa) Hubble types.Comment: To appear in ApJ. 20 pages LaTex 2-column with 3 tables, 23 figures
(GIF) available at this site. Complete gzipped postscript preprint with
embedded figures available from http://tarkus.pha.jhu.edu/~thakar/cr3.html (3
Mb
The mass function
We present the mass functions for different mass estimators for a range of
cosmological models. We pay particular attention to how universal the mass
function is, and how it depends on the cosmology, halo identification and mass
estimator chosen. We investigate quantitatively how well we can relate observed
masses to theoretical mass functions.Comment: 14 pages, 12 figures, to appear in ApJ
The Structure of Isothermal, Self-gravitating Gas Spheres for Softened Gravity
A theory for the structure of isothermal, self-gravitating gas spheres in
pressure equilibrium in a softened gravitational field is developed. The one
parameter spline softening proposed by Hernquist & Katz (1989) is used. We show
that the addition of this extra scale parameter implies that the set of
equilibrium solutions constitute a one-parameter family, rather than the one
and only one isothermal sphere solution for Newtonian gravity. We demonstrate
the perhaps somewhat surprising result that for any finite choice of softening
length and temperature, it is possible to deposit an arbitrarily large mass of
gas in pressure equilibrium and with a non-singular density distribution inside
of r_0 for any r_0 > 0. The theoretical predictions of our models are compared
with the properties of the small, massive, quasi-isothermal gas clumps which
typically form in numerical Tree-SPH simulations of 'passive' galaxy formation
of Milky Way sized galaxies. We find reasonable agreement despite the neglect
of rotational support in the models. We comment on whether the hydrodynamical
resolution in our numerical simulation of galaxy formation is sufficient, and
finally we conclude that one should be cautious, when comparing results of
numerical simulations involving gravitational softening and hydrodynamical
smoothing, with reality.Comment: 22 pages Latex + 12 figure
Tidal disruption of dark matter halos around proto-globular clusters
Tidal disruption of dark matter halos around proto-globular clusters in a
halo of a small galaxy is studied in the context of the hierarchical clustering
scenario by using semi-cosmological N-body/SPH simulations assuming the
standard cold dark matter model (). Our analysis on formation and
evolution of the galaxy and its substructures archives until . In such
a high-redshift universe, the Einstein-de Sitter universe is still a good
approximation for a recently favored -dominated universe, and then our
results does not depend on the choice of cosmology. In order to resolve small
gravitationally-bound clumps around galaxies and consider radiative cooling
below , we adopt a fine mass resolution (m_{\rm SPH} = 1.12 \times
10^3 \Msun). Because of the cooling, each clump immediately forms a
`core-halo' structure which consists of a baryonic core and a dark matter halo.
The tidal force from the host galaxy mainly strips the dark matter halo from
clumps and, as a result, theses clumps get dominated by baryons. Once a clump
is captured by the host halo, its mass drastically decreases each pericenter
passage. At , more than half of the clumps become baryon dominated
systems (baryon mass/total mass ). Our results support the tidal
evolution scenario of the formation of globular clusters and baryon dominated
dwarf galaxies in the context of the cold dark matter universe.Comment: 9page, 13 figures. Accepted for publication in ApJ. A high-resolution
PDF of the paper can be obtained from http://th.nao.ac.jp/~takayuki/ApJ05
- âŠ