176 research outputs found
One-Point Probability Distribution Functions of Supersonic Turbulent Flows in Self-Gravitating Media
Turbulence is essential for understanding the structure and dynamics of
molecular clouds and star-forming regions. There is a need for adequate tools
to describe and characterize the properties of turbulent flows. One-point
probability distribution functions (pdf's) of dynamical variables have been
suggested as appropriate statistical measures and applied to several observed
molecular clouds. However, the interpretation of these data requires comparison
with numerical simulations. To address this issue, SPH simulations of driven
and decaying, supersonic, turbulent flows with and without self-gravity are
presented. In addition, random Gaussian velocity fields are analyzed to
estimate the influence of variance effects. To characterize the flow
properties, the pdf's of the density, of the line-of-sight velocity centroids,
and of the line centroid increments are studied. This is supplemented by a
discussion of the dispersion and the kurtosis of the increment pdf's, as well
as the spatial distribution of velocity increments for small spatial lags. From
the comparison between different models of interstellar turbulence, it follows
that the inclusion of self-gravity leads to better agreement with the observed
pdf's in molecular clouds. The increment pdf's for small spatial lags become
exponential for all considered velocities. However, all the processes
considered here lead to non-Gaussian signatures, differences are only gradual,
and the analyzed pdf's are in addition projection dependent. It appears
therefore very difficult to distinguish between different physical processes on
the basis of pdf's only, which limits their applicability for adequately
characterizing interstellar turbulence.Comment: 38 pages (incl. 17 figures), accepted for publication in ApJ, also
available with full resolution figures at
http://www.strw.leidenuniv.nl/~klessen/Preprint
Gravitational Collapse in Turbulent Molecular Clouds. I. Gasdynamical Turbulence
Observed molecular clouds often appear to have very low star formation
efficiencies and lifetimes an order of magnitude longer than their free-fall
times. Their support is attributed to the random supersonic motions observed in
them. We study the support of molecular clouds against gravitational collapse
by supersonic, gas dynamical turbulence using direct numerical simulation.
Computations with two different algorithms are compared: a particle-based,
Lagrangian method (SPH), and a grid-based, Eulerian, second-order method
(ZEUS). The effects of both algorithm and resolution can be studied with this
method. We find that, under typical molecular cloud conditions, global collapse
can indeed be prevented, but density enhancements caused by strong shocks
nevertheless become gravitationally unstable and collapse into dense cores and,
presumably, stars. The occurance and efficiency of local collapse decreases as
the driving wave length decreases and the driving strength increases. It
appears that local collapse can only be prevented entirely with unrealistically
short wave length driving, but observed core formation rates can be reproduced
with more realistic driving. At high collapse rates, cores are formed on short
time scales in coherent structures with high efficiency, while at low collapse
rates they are scattered randomly throughout the region and exhibit
considerable age spread. We suggest that this naturally explains the observed
distinction between isolated and clustered star formation.Comment: Minor revisions in response to referee, thirteen figures, accepted to
Astrophys.
The structure of self-gravitating clouds
To study the interaction of star-formation and turbulent molecular cloud
structuring, we analyse numerical models and observations of self-gravitating
clouds using the Delta-variance as statistical measure for structural
characteristics. In the models we resolve the transition from purely
hydrodynamic turbulence to gravitational collapse associated with the formation
and mass growth of protostellar cores. We compare models of driven and freely
decaying turbulence with and without magnetic fields. Self-gravitating
supersonic turbulence always produces a density structure that contains most
power on the smallest scales provided by collapsed cores as soon as local
collapse sets in. This is in contrast to non-self-gravitating hydrodynamic
turbulence where the Delta-variance is dominated by large scale structures. To
detect this effect in star-forming regions observations have to resolve the
high density contrast of protostellar cores with respect to their ambient
molecular cloud. Using the 3mm continuum map of a star-forming cluster in
Serpens we show that the dust emission traces the full density evolution. On
the contrary, the density range accessible by molecular line observations is
insufficient for this analysis. Only dust emission and dust extinction
observations are able to to determine the structural parameters of star-forming
clouds following the density evolution during the gravitational collapse.Comment: 12 pages, 9 figures, A&A in pres
Dynamical Expansion of H II Regions from Ultracompact to Compact Sizes in Turbulent, Self-Gravitating Molecular Clouds
The nature of ultracompact H II regions (UCHRs) remains poorly determined. In
particular, they are about an order of magnitude more common than would be
expected if they formed around young massive stars and lasted for one dynamical
time, around 10^4 yr. We here perform three-dimensional numerical simulations
of the expansion of an H II region into self-gravitating, radiatively cooled
gas, both with and without supersonic turbulent flows. In the laminar case, we
find that H II region expansion in a collapsing core produces nearly spherical
shells, even if the ionizing source is off-center in the core. This agrees with
analytic models of blast waves in power-law media. In the turbulent case, we
find that the H II region does not disrupt the central collapsing region, but
rather sweeps up a shell of gas in which further collapse occurs. Although this
does not constitute triggering, as the swept-up gas would eventually have
collapsed anyway, it does expose the collapsing regions to ionizing radiation.
We suggest that these regions of secondary collapse, which will not all
themselves form massive stars, may form the bulk of observed UCHRs. As the
larger shell will take over 10^5 years to complete its evolution, this could
solve the timescale problem. Our suggestion is supported by the ubiquitous
observation of more diffuse emission surrounding UCHRs.Comment: accepted to ApJ, 40 pages, 13 b/w figures, changes from v1 include
analytic prediction of radio luminosity, better description of code testing,
and many minor changes also in response to refere
Simulating Stellar Merger using HPX/Kokkos on A64FX on Supercomputer Fugaku
The increasing availability of machines relying on non-GPU architectures,
such as ARM A64FX in high-performance computing, provides a set of interesting
challenges to application developers. In addition to requiring code portability
across different parallelization schemes, programs targeting these
architectures have to be highly adaptable in terms of compute kernel sizes to
accommodate different execution characteristics for various heterogeneous
workloads. In this paper, we demonstrate an approach to code and performance
portability that is based entirely on established standards in the industry. In
addition to applying Kokkos as an abstraction over the execution of compute
kernels on different heterogeneous execution environments, we show that the use
of standard C++ constructs as exposed by the HPX runtime system enables superb
portability in terms of code and performance based on the real-world Octo-Tiger
astrophysics application. We report our experience with porting Octo-Tiger to
the ARM A64FX architecture provided by Stony Brook's Ookami and Riken's
Supercomputer Fugaku and compare the resulting performance with that achieved
on well established GPU-oriented HPC machines such as ORNL's Summit, NERSC's
Perlmutter and CSCS's Piz Daint systems. Octo-Tiger scaled well on
Supercomputer Fugaku without any major code changes due to the abstraction
levels provided by HPX and Kokkos. Adding vectorization support for ARM's SVE
to Octo-Tiger was trivial thanks to using standard C+
- …