213 research outputs found
Mean-field and direct numerical simulations of magnetic flux concentrations from vertical field
Strongly stratified hydromagnetic turbulence has previously been found to
produce magnetic flux concentrations if the domain is large enough compared
with the size of turbulent eddies. Mean-field simulations (MFS) using
parameterizations of the Reynolds and Maxwell stresses show a negative
effective magnetic pressure instability and have been able to reproduce many
aspects of direct numerical simulations (DNS) regarding the growth rate of this
large-scale instability, shape of the resulting magnetic structures, and their
height as a function of magnetic field strength. Unlike the case of an imposed
horizontal field, for a vertical one, magnetic flux concentrations of
equipartition strength with the turbulence can be reached. This results in
magnetic spots that are reminiscent of sunspots. Here we want to find out under
what conditions magnetic flux concentrations with vertical field occur and what
their internal structure is. We use a combination of MFS, DNS, and implicit
large-eddy simulations to characterize the resulting magnetic flux
concentrations in forced isothermal turbulence with an imposed vertical
magnetic field. We confirm earlier results that in the kinematic stage of the
large-scale instability the horizontal wavelength of structures is about 10
times the density scale height. At later times, even larger structures are
being produced in a fashion similar to inverse spectral transfer in helically
driven turbulence. Using turbulence simulations, we find that magnetic flux
concentrations occur for different values of the Mach number between 0.1 and
0.7. DNS and MFS show magnetic flux tubes with mean-field energies comparable
to the turbulent kinetic energy. The resulting vertical magnetic flux tubes are
being confined by downflows along the tubes and corresponding inflow from the
sides, which keep the field concentrated.Comment: 16 pages, 22 figures, Astron. Astrophys., in pres
New scaling for the alpha effect in slowly rotating turbulence
Using simulations of slowly rotating stratified turbulence, we show that the
alpha effect responsible for the generation of astrophysical magnetic fields is
proportional to the logarithmic gradient of kinetic energy density rather than
that of momentum, as was previously thought. This result is in agreement with a
new analytic theory developed in this paper for large Reynolds numbers. Thus,
the contribution of density stratification is less important than that of
turbulent velocity. The alpha effect and other turbulent transport coefficients
are determined by means of the test-field method. In addition to forced
turbulence, we also investigate supernova-driven turbulence and stellar
convection. In some cases (intermediate rotation rate for forced turbulence,
convection with intermediate temperature stratification, and supernova-driven
turbulence) we find that the contribution of density stratification might be
even less important than suggested by the analytic theory.Comment: 10 pages, 9 figures, revised version, Astrophys. J., in pres
Linear analysis of the vertical shear instability: outstanding issues and improved solutions
This equipment is funded by BIS National E-Infrastructure capital grant ST/K000373/1 and STFC DiRAC Operations grant ST/K0003259/1. DiRAC is part of the National E-Infrastructure
A Field-length based refinement criterion for adaptive mesh simulations of the interstellar medium
Adequate modelling of the multiphase interstellar medium requires optically
thin radiative cooling, comprising an inherent thermal instability. The size of
the occurring condensation and evaporation interfaces is determined by the
so-called Field-length, which gives the dimension at which the instability is
significantly damped by thermal conduction. Our aim is to study the relevance
of conduction scale effects in the numerical modelling of a bistable medium and
check the applicability of conventional and alternative adaptive mesh
techniques. The low physical value of the thermal conduction within the ISM
defines a multiscale problem, hence promoting the use of adaptive meshes. We
here introduce a new refinement strategy that applies the Field condition by
Koyama & Inutsuka as a refinement criterion. The described method is very
similar to the Jeans criterion for gravitational instability by Truelove and
efficiently allows to trace the unstable gas situated at the thermal
interfaces. We present test computations that demonstrate the greater accuracy
of the newly proposed refinement criterion in comparison to refinement based on
the local density gradient. Apart from its usefulness as a refinement trigger,
we do not find evidence in favour of the Field criterion as a prerequisite for
numerical stability.Comment: 5 pages, 5 figures, accepted for publication in A&
Low-mass planet migration in three-dimensional wind-driven inviscid discs: a negative corotation torque
We present simulations of low-mass planet–disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems
Global Hydromagnetic Simulations of Protoplanetary Disks with Stellar Irradiation and Simplified Thermochemistry
Outflows driven by large-scale magnetic fields likely play an important role
in the evolution and dispersal of protoplanetary disks, and in setting the
conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD
model of these outflows by incorporating radiative transfer and simplified
thermochemistry, with the twin aims of exploring how heating influences wind
launching, and illustrating how such models can be tested through observations
of diagnostic spectral lines. Our model disks launch magnetocentrifugal
outflows primarily through magnetic tension forces, so the mass-loss rate
increases only moderately when thermochemical effects are switched on. For
typical field strengths, thermochemical and irradiation heating are more
important than magnetic dissipation. We furthermore find that the entrained
vertical magnetic flux diffuses out of the disk on secular timescales as a
result of non-ideal MHD. Through post-processing line radiative transfer, we
demonstrate that spectral line intensities and moment-1 maps of atomic oxygen,
the HCN molecule, and other species show potentially observable differences
between a model with a magnetically driven outflow and one with a weaker,
photoevaporative outflow. In particular, the line shapes and velocity
asymmetries in the moment-1 maps could enable the identification of outflows
emanating from the disk surface.Comment: 35 pages, 20 figures, accepted for publication in Ap
Direct simulations of a supernova-driven galactic dynamo
Supernovae are known to be the dominant energy source for driving turbulence
in the interstellar medium. Yet, their effect on magnetic field amplification
in spiral galaxies is still poorly understood. Previous analytical models,
based on the evolution of isolated, non-interacting supernova remnants,
predicted a dominant vertical pumping that would render dynamo action
improbable. In the present work, we address the issue of vertical transport,
which is thought to be the key process that inhibits dynamo action in the
galactic context. We aim to demonstrate that supernova driving is a powerful
mechanism to amplify galactic magnetic fields. We conduct direct numerical
simulations in the framework of resistive magnetohydrodynamics. Our local box
model of the interstellar medium comprises optically-thin radiative cooling, an
external gravitational potential, and background shear. Dynamo coefficients for
mean-field models are measured by means of passive test fields. Our simulations
show that supernova-driven turbulence in conjunction with shear leads to an
exponential amplification of the mean magnetic field. We found turbulent
pumping to be directed inward and approximately balanced by a galactic wind.Comment: minor changes, 4 pages, 4 figures, accepted for publication in A&A
Shearingbox-implementation for the central-upwind, constraint-transport MHD-code NIRVANA
We describe the implementation of the shearingbox approach into the
Godunov-type central-upwind/constraint-transport magnetohydrodynamics code
NIRVANA. This will allow for applications which require sheared-periodic
boundary conditions as typically used in local Cartesian simulations of
differentially rotating systems. We present the algorithm in detail and discuss
necessary modifications in the numerical fluxes in order to preserve conserved
quantities and to fulfill other analytical constraints as good as seem feasible
within the numerical scheme. We check the source terms which come with the
shearingbox formulation by investigating the conservation of the epicyclic mode
energy. We also perform more realistic simulations of the magneto-rotational
instability with initial zero-net-flux vertical magnetic field and compare the
obtained stresses and energetics with previous non-conservative results
exploring the same parameter regime.Comment: 17 pages, 4 figures, to be published in CP
Galaxies in box: A simulated view of the interstellar medium
We review progress in the development of physically realistic three
dimensional simulated models of the galaxy.We consider the scales from star
forming molecular clouds to the full spiral disc. Models are computed using
hydrodynamic (HD) or magnetohydrodynamic (MHD) equations and may include cosmic
ray or tracer particles. The range of dynamical scales between the full galaxy
structure and the turbulent scales of supernova (SN) explosions and even cloud
collapse to form stars, make it impossible with current computing tools and
resources to resolve all of these in one model. We therefore consider a
hierarchy of models and how they can be related to enhance our understanding of
the complete galaxy.Comment: Chapter in Large Scale Magnetic Fields in the Univers
- …