18 research outputs found
Density probability distribution in one-dimensional polytropic gas dynamics
We discuss the generation and statistics of the density fluctuations in
highly compressible polytropic turbulence, based on a simple model and
one-dimensional numerical simulations. Observing that density structures tend
to form in a hierarchical manner, we assume that density fluctuations follow a
random multiplicative process. When the polytropic exponent is equal
to unity, the local Mach number is independent of the density, and our
assumption leads us to expect that the probability density function (PDF) of
the density field is a lognormal. This isothermal case is found to be singular,
with a dispersion which scales like the square turbulent Mach
number , where and is the fluid density.
This leads to much higher fluctuations than those due to shock jump relations.
Extrapolating the model to the case , we find that, as the
Mach number becomes large, the density PDF is expected to asymptotically
approach a power-law regime, at high densities when , and at low
densities when . This effect can be traced back to the fact that the
pressure term in the momentum equation varies exponentially with , thus
opposing the growth of fluctuations on one side of the PDF, while being
negligible on the other side. This also causes the dispersion to
grow more slowly than when . In view of these
results, we suggest that Burgers flow is a singular case not approached by the
high- limit, with a PDF that develops power laws on both sides.Comment: 9 pages + 12 postscript figures. Submitted to Phys. Rev.
Interstellar MHD Turbulence and Star Formation
This chapter reviews the nature of turbulence in the Galactic interstellar
medium (ISM) and its connections to the star formation (SF) process. The ISM is
turbulent, magnetized, self-gravitating, and is subject to heating and cooling
processes that control its thermodynamic behavior. The turbulence in the warm
and hot ionized components of the ISM appears to be trans- or subsonic, and
thus to behave nearly incompressibly. However, the neutral warm and cold
components are highly compressible, as a consequence of both thermal
instability in the atomic gas and of moderately-to-strongly supersonic motions
in the roughly isothermal cold atomic and molecular components. Within this
context, we discuss: i) the production and statistical distribution of
turbulent density fluctuations in both isothermal and polytropic media; ii) the
nature of the clumps produced by thermal instability, noting that, contrary to
classical ideas, they in general accrete mass from their environment; iii) the
density-magnetic field correlation (or lack thereof) in turbulent density
fluctuations, as a consequence of the superposition of the different wave modes
in the turbulent flow; iv) the evolution of the mass-to-magnetic flux ratio
(MFR) in density fluctuations as they are built up by dynamic compressions; v)
the formation of cold, dense clouds aided by thermal instability; vi) the
expectation that star-forming molecular clouds are likely to be undergoing
global gravitational contraction, rather than being near equilibrium, and vii)
the regulation of the star formation rate (SFR) in such gravitationally
contracting clouds by stellar feedback which, rather than keeping the clouds
from collapsing, evaporates and diperses them while they collapse.Comment: 43 pages. Invited chapter for the book "Magnetic Fields in Diffuse
Media", edited by Elisabete de Gouveia dal Pino and Alex Lazarian. Revised as
per referee's recommendation
The chemical enrichment of the ICM from hydrodynamical simulations
The study of the metal enrichment of the intra-cluster and inter-galactic
media (ICM and IGM) represents a direct means to reconstruct the past history
of star formation, the role of feedback processes and the gas-dynamical
processes which determine the evolution of the cosmic baryons. In this paper we
review the approaches that have been followed so far to model the enrichment of
the ICM in a cosmological context. While our presentation will be focused on
the role played by hydrodynamical simulations, we will also discuss other
approaches based on semi-analytical models of galaxy formation, also critically
discussing pros and cons of the different methods. We will first review the
concept of the model of chemical evolution to be implemented in any
chemo-dynamical description. We will emphasise how the predictions of this
model critically depend on the choice of the stellar initial mass function, on
the stellar life-times and on the stellar yields. We will then overview the
comparisons presented so far between X-ray observations of the ICM enrichment
and model predictions. We will show how the most recent chemo-dynamical models
are able to capture the basic features of the observed metal content of the ICM
and its evolution. We will conclude by highlighting the open questions in this
study and the direction of improvements for cosmological chemo-dynamical models
of the next generation.Comment: 25 pages, 11 figures, accepted for publication in Space Science
Reviews, special issue "Clusters of galaxies: beyond the thermal view",
Editor J.S. Kaastra, Chapter 18; 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
Theory and Applications of Non-Relativistic and Relativistic Turbulent Reconnection
Realistic astrophysical environments are turbulent due to the extremely high
Reynolds numbers. Therefore, the theories of reconnection intended for
describing astrophysical reconnection should not ignore the effects of
turbulence on magnetic reconnection. Turbulence is known to change the nature
of many physical processes dramatically and in this review we claim that
magnetic reconnection is not an exception. We stress that not only
astrophysical turbulence is ubiquitous, but also magnetic reconnection itself
induces turbulence. Thus turbulence must be accounted for in any realistic
astrophysical reconnection setup. We argue that due to the similarities of MHD
turbulence in relativistic and non-relativistic cases the theory of magnetic
reconnection developed for the non-relativistic case can be extended to the
relativistic case and we provide numerical simulations that support this
conjecture. We also provide quantitative comparisons of the theoretical
predictions and results of numerical experiments, including the situations when
turbulent reconnection is self-driven, i.e. the turbulence in the system is
generated by the reconnection process itself. We show how turbulent
reconnection entails the violation of magnetic flux freezing, the conclusion
that has really far reaching consequences for many realistically turbulent
astrophysical environments. In addition, we consider observational testing of
turbulent reconnection as well as numerous implications of the theory. The
former includes the Sun and solar wind reconnection, while the latter include
the process of reconnection diffusion induced by turbulent reconnection, the
acceleration of energetic particles, bursts of turbulent reconnection related
to black hole sources as well as gamma ray bursts. Finally, we explain why
turbulent reconnection cannot be explained by turbulent resistivity or derived
through the mean field approach.Comment: 66 pages, 24 figures, a chapter of the book "Magnetic Reconnection -
Concepts and Applications", editors W. Gonzalez, E. N. Parke
Gravitational Collapse and Disk Formation in Magnetized Cores
We discuss the effects of the magnetic field observed in molecular clouds on
the process of star formation, concentrating on the phase of gravitational
collapse of low-mass dense cores, cradles of sunlike stars. We summarize recent
analytic work and numerical simulations showing that a substantial level of
magnetic field diffusion at high densities has to occur in order to form
rotationally supported disks. Furthermore, newly formed accretion disks are
threaded by the magnetic field dragged from the parent core during the
gravitational collapse. These disks are expected to rotate with a sub-Keplerian
speed because they are partially supported by magnetic tension against the
gravity of the central star. We discuss how sub-Keplerian rotation makes it
difficult to eject disk winds and accelerates the process of planet migration.
Moreover, magnetic fields modify the Toomre criterion for gravitational
instability via two opposing effects: magnetic tension and pressure increase
the disk local stability, but sub-Keplerian rotation makes the disk more
unstable. In general, magnetized disks are more stable than their nonmagnetic
counterparts; thus, they can be more massive and less prone to the formation of
giant planets by gravitational instability.Comment: Chapter 16 in "Magnetic Fields in Diffuse Media", Springer-Verlag,
eds. de Gouveia Dal Pino, E., Lazarian, A., Melioli,
The Physics of Star Cluster Formation and Evolution
© 2020 Springer-Verlag. The final publication is available at Springer via https://doi.org/10.1007/s11214-020-00689-4.Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work.Peer reviewe
Mixed model-based additive models for sample extremes
We consider additive models fitting and inference when the response variable is a sample extreme. Non-linear covariate effects are handled using the mixed model representation of penalised splines. A fitting algorithm based on likelihood approximations is derived. The efficacy of the resulting methodology is demonstrated via application to simulated and real data.