89 research outputs found
Accuracy requirements to test the applicability of the random cascade model to supersonic turbulence
A model, which is widely used for inertial rang statistics of supersonic
turbulence in the context of molecular clouds and star formation, expresses
(measurable) relative scaling exponents Z_p of two-point velocity statistics as
a function of two parameters, beta and Delta. The model relates them to the
dimension D of the most dissipative structures, D=3-Delta/(1-beta). While this
description has proved most successful for incompressible turbulence
(beta=Delta=2/3, and D=1), its applicability in the highly compressible regime
remains debated. For this regime, theoretical arguments suggest D=2 and
Delta=2/3, or Delta=1. Best estimates based on 3D periodic box simulations of
supersonic isothermal turbulence yield Delta=0.71 and D=1.9, with uncertainty
ranges of Delta in [0.67, 0.78] and D in [2.04,1.60]. With these 5-10\%
uncertainty ranges just marginally including the theoretical values of
Delta=2/3 and D=2, doubts remain whether the model indeed applies and, if it
applies, for what values of beta and Delta. We use a Monte Carlo approach to
mimic actual simulation data and examine what factors are most relevant for the
fit quality. We estimate that 0.1% (0.05%) accurate Z_p, with p=1...5, should
allow for 2% (1%) accurate estimates of beta and Delta in the highly
compressible regime, but not in the mildly compressible regime. We argue that
simulation-based Z_p with such accuracy are within reach of today's computer
resources. If this kind of data does not allow for the expected high quality
fit of beta and Delta, then this may indicate the inapplicability of the model
for the simulation data. In fact, other models than the one we examine here
have been suggested.Comment: 8 pages, 8 figures, accepted by Astronomy and Astrophysic
Supersonic turbulence in 3D isothermal flow collision
Colliding supersonic bulk flows shape observable properties and internal
physics of various astrophysical objects, like O-star winds, molecular clouds,
galactic sheets, binaries, or gamma-ray bursts. Using numerical simulations, we
show that the bulk flows leave a clear imprint on the collision zone, its mean
properties and the turbulence it naturally develops. Our model setup consists
of 3D head-on colliding isothermal hydrodynamical flows with Mach numbers
between 2 and 43. Simulation results are in line with expectations from
self-similarity: root mean square Mach numbers (Mrms) scale linearly with
upstream Mach numbers, mean densities remain limited to a few times the
upstream density. The density PDF is not log-normal. The turbulence is
inhomogeneous: weaker in the zone center than close to the confining shocks. It
is anisotropic: while Mrms is generally supersonic, Mrms transverse to the
upstream flow is always subsonic. We argue that uniform, isothermal, head-on
colliding flows generally disfavor isotropic, supersonic turbulence. The
anisotropy carries over to other quantities like the density variance - Mach
number relation. Structure functions differ depending on whether they are
computed along a line-of-sight perpendicular or parallel to the upstream flow.
We suggest that such line-of-sight effects should be kept in mind when
interpreting turbulence characteristics derived from observations.Comment: 20 pages, 14 figures, 4 tables, accepted by Astronomy and
Astrophysic
Relativistic magnetic reconnection in collisionless ion-electron plasmas explored with particle-in-cell simulations
Magnetic reconnection is a leading mechanism for magnetic energy conversion
and high-energy non-thermal particle production in a variety of high-energy
astrophysical objects, including ones with relativistic ion-electron plasmas
(e.g., microquasars or AGNs) - a regime where first principle studies are
scarce. We present 2D particle-in-cell (PIC) simulations of low
ion-electron plasmas under relativistic conditions, i.e., with inflow magnetic
energy exceeding the plasma rest-mass energy. We identify outstanding
properties: (i) For relativistic inflow magnetizations (here ), the reconnection outflows are dominated by thermal agitation instead of
bulk kinetic energy. (ii) At large inflow electron magnetization (), the reconnection electric field is sustained more by bulk inertia than by
thermal inertia. It challenges the thermal-inertia-paradigm and its
implications. (iii) The inflows feature sharp transitions at the entrance of
the diffusion zones. These are not shocks but results from particle ballistic
motions, all bouncing at the same location, provided that the thermal velocity
in the inflow is far smaller than the inflow E cross B bulk velocity. (iv)
Island centers are magnetically isolated from the rest of the flow, and can
present a density depletion at their center. (v) The reconnection rates are
slightly larger than in non-relativistic studies. They are best normalized by
the inflow relativistic Alfv\'en speed projected in the outflow direction,
which then leads to rates in a close range (0.14-0.25) thus allowing for an
easy estimation of the reconnection electric field.Comment: Submitted to A&
The energetics of relativistic magnetic reconnection: ion-electron repartition and particle distribution hardness
Collisionless magnetic reconnection is a prime candidate to account for
flare-like or steady emission, outflow launching, or plasma heating, in a
variety of high-energy astrophysical objects, including ones with relativistic
ion-electron plasmas. But the fate of the initial magnetic energy in a
reconnection event remains poorly known: what is the amount given to kinetic
energy, the ion/electron repartition, and the hardness of the particle
distributions? We explore these questions with 2D particle-in-cell simulations
of ion-electron plasmas. We find that 45 to 75% of the total initial magnetic
energy ends up in kinetic energy, this fraction increasing with the inflow
magnetization. Depending on the guide field strength, ions get from 30 to 60%
of the total kinetic energy. Particles can be separated into two populations
that only weakly mix: (i) particles initially in the current sheet, heated by
its initial tearing and subsequent contraction of the islands; and (ii)
particles from the background plasma that primarily gain energy via the
reconnection electric field when passing near the X-point. Particles (ii) tend
to form a power-law with an index , that
depends mostly on the inflow Alfv\'en speed and magnetization
of species , with for electrons to for increasing .
The highest particle Lorentz factor, for ions or electrons, increases roughly
linearly with time for all the relativistic simulations. This is faster, and
the spectra can be harder, than for collisionless shock acceleration. We
discuss applications to microquasar and AGN coronae, to extragalactic jets, and
to radio lobes. We point out situations where effects such as Compton drag or
pair creation are important.Comment: 15 pages, submitted to A&
- …