BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Oliver Gressel

Oliver Gressel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.GAastro-ph.HEastro-ph.SRastro-phastro-ph.CO

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

Finite-time response of dynamo mean-field effects in magnetorotational turbulence

Accretion disc turbulence along with its effect on large-scale magnetic fields plays an important role in understanding disc evolution in general, and the launching of astrophysical jets in particular. Motivated by enabling a comprehensive sub-grid description for global long-term simulations of accretions discs, we aim to further characterize the transport coefficients emerging in local simulations of magnetorotational disc turbulence. For the current investigation, we leverage a time-dependent version of the test-field method, which is sensitive to the turbulent electromotive force (EMF) generated as a response to a set of pulsating background fields. We obtain Fourier spectra of the transport coefficients as a function of oscillation frequency. These are well approximated by a simple response function, describing a finite-time build-up of the EMF as a result of a time-variable mean magnetic field. For intermediate timescales (i.e., slightly above the orbital frequency), we observe a significant phase lag of the EMF compared to the causing field. Augmented with our previous result on a non-local closure relation in space, and incorporated into a suitable mean-field description that we briefly sketch out here, the new framework will allow to drop the restrictive assumption of scale separation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Role of Disk Winds in the Evolution and Dispersal of Protoplanetary Disks

The assembly and architecture of planetary systems strongly depend on the physical processes governing the evolution and dispersal of protoplanetary disks. Since Protostars and Planets VI, new observations and theoretical insights favor disk winds as being one of those key processes. This chapter provides a comprehensive review of recent observations probing outflowing gas launched over a range of disk radii for a wide range of evolutionary stages, enabling an empirical understanding of how winds evolve. In parallel, we review theoretical advancements in both magnetohydrodynamic and photoevaporative disk wind models and identify predictions that can be confronted with observations. By linking theory and observations we critically assess the role of disk winds in the evolution and dispersal of protoplanetary disks. Finally, we explore the impact of disk winds on planet formation and evolution and highlight theoretical work, observations, and critical tests for future progress.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the spatial and temporal non-locality of dynamo mean-field effects in supersonic interstellar turbulence

The interstellar medium of the Milky Way and nearby disk galaxies harbours large-scale coherent magnetic fields of Microgauss strength, that can be explained via the action of a mean-field dynamo. As in our previous work, we aim to quantify dynamo effects that are self-consistently emerging in realistic direct magnetohydrodynamic simulations, but we generalise our approach to the case of a non-local (non-instantaneous) closure relation, described by a convolution integral in space (time). To this end, we leverage our comprehensive simulation framework for the supernova-regulated turbulent multi-phase interstellar medium. By introducing spatially (temporally) modulated mean fields, we extend the previously used test-field method to the spectral realm -- providing the Fourier representation of the convolution kernels. The resulting spectra of the dynamo mean-field coefficients that we obtain broadly match expectations and allow to rigorously constrain the degree of scale separation in the Galactic dynamo. A surprising result is found for the diamagnetic pumping term, which increases in amplitude when going to smaller scales. Our results amount to the most comprehensive description of dynamo mean-field effects in the Galactic context to date. Surveying the relevant parameter space and quenching behaviour, this will ultimately enable the development of assumption-free sub-grid prescriptions for otherwise unresolved global galaxy simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Turbulent transport coefficients in galactic dynamo simulations using singular value decomposition

Coherent magnetic fields in disc galaxies are thought to be generated by a large-scale (or mean-field) dynamo operating in their interstellar medium. A key driver of mean magnetic field growth is the turbulent electromotive force (EMF), which represents the influence of correlated small-scale (or fluctuating) velocity and magnetic fields on the mean field. The EMF is usually expressed as a linear expansion in the mean magnetic field and its derivatives, with the dynamo tensors as expansion coefficients. Here, we adopt the singular value decomposition (SVD) method to directly measure these turbulent transport coefficients in a simulation of the turbulent interstellar medium that realizes a large-scale dynamo. Specifically, the SVD is used to least-square fit the time series data of the EMF with that of the mean field and its derivatives, to determine these coefficients. We demonstrate that the spatial profiles of the EMF reconstructed from the SVD coefficients match well with that taken directly from the simulation. Also, as a direct test, we use the coefficients to simulate a 1-D mean-field dynamo model and find an overall similarity in the evolution of the mean magnetic field between the dynamo model and the direct simulation. We also compare the results with those which arise using simple regression and the ones obtained previously using the test-field (TF) method, to find reasonable qualitative agreement. Overall, the SVD method provides an effective post-processing tool to determine turbulent transport coefficients from simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Misaligned Disks in the Binary Protostar IRS 43

Recent high angular resolution ($\sim$0.2") ALMA observations of the 1.1 mm continuum and of HCO+ J=3-2 and HCN J=3-2 gas towards the binary protostar IRS 43 reveal multiple Keplerian disks which are significantly misaligned ($\gt$ 60$^\circ$), both in inclination and position angle and also with respect to the binary orbital plane. Each stellar component has an associated circumstellar disk while the binary is surrounded by a circumbinary disk. Together with archival VLA measurements of the stellar positions over 25 years, and assuming a circular orbit, we use our continuum measurements to determine the binary separation, a = 74 $\pm$ 4 AU, and its inclination, i $\lt$ 30$^\circ$. The misalignment in this system suggests that turbulence has likely played a major role in the formation of IRS 43.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Characterizing the mean-field dynamo in turbulent accretion disks

The formation and evolution of a wide class of astrophysical objects is governed by turbulent, magnetized accretion disks. Understanding their secular dynamics is of primary importance. Apart from enabling mass accretion via the transport of angular momentum, the turbulence affects the long-term evolution of the embedded magnetic flux, which in turn regulates the efficiency of the transport. In this paper, we take a comprehensive next step towards an effective mean-field model for turbulent astrophysical disks by systematically studying the key properties of magnetorotational turbulence in vertically-stratified, isothermal shearing boxes. This allows us to infer emergent properties of the ensuing chaotic flow as a function of the shear parameter as well as the amount of net-vertical flux. Using the test-field method, we furthermore characterize the mean-field dynamo coefficients that describe the long-term evolution of large-scale fields. We simultaneously infer the vertical shape and the spectral scale dependence of these closure coefficients, with the latter describing non-local contributions to the turbulent electromotive force. Based on this, we infer a scale-separation ratio of about ten for the large-scale dynamo. We finally determine scaling properties of the mean-field dynamo coefficients. The relevant component of the dynamo α effect is found to scale linearly with the shear rate, as is the corresponding turbulent diffusion, η. Together, these scalings allow us to predict, in a quantitative manner, the cycle period of the well-known butterfly diagram. This lends new support to the importance of the αΩ mechanism in determining the evolution of large-scale magnetic fields in turbulent accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Global simulations of protoplanetary disks with ohmic resistivity and ambipolar diffusion

Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Do magnetic fields influence gas rotation in galaxies?

We aim to estimate the contribution of the radial component of the Lorentz force to the gas rotation in several types of galaxies. Using typical parameters for the exponential scale of synchrotron emission and the scale length of HI gas, under the assumption of equipartition between the energies of cosmic rays and total magnetic fields, we derive the Lorentz force and compare it to the gravitational force in the radial component of the momentum equation. We distinguish the different contributions between the large-scale and the small-scale turbulent fields by Reynolds averaging. We compare these findings with a dynamical dynamo model. We find a possible reduction of circular gas velocity in the very outer parts and an increase inside a radius of four times the synchrotron scale length. Sufficiently localized radial reversals of the magnetic field may cause characteristic modulations in the gas rotation curve with typical amplitudes of 10-20 km/s. It is unlikely that the magnetic field contributes to the flat rotation in the outer parts of galaxies. If anything, it will \emph{impede} the gravitationally supported rotation, demanding for an even higher halo mass to explain the observed rotation profile. We speculate that this may have consequences for ram pressure stripping and the truncation of the stellar disc.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamo effects in magnetorotational turbulence with finite thermal diffusivity

We investigate the saturation level of hydromagnetic turbulence driven by the magnetorotational instability in the case of vanishing net flux. Motivated by a recent paper of Bodo, Cattaneo, Mignone, & Rossi, we here focus on the case of a non-isothermal equation of state with constant thermal diffusivity. The central aim of the paper is to complement the previous result with closure parameters for mean-field dynamo models, and to test the hypothesis that the dynamo is affected by the mode of heat transport. We perform computer simulations of local shearing-box models of stratified accretion disks with approximate treatment of radiative heat transport, which is modeled via thermal conduction. We study the effect of varying the (constant) thermal diffusivity, and apply different vertical boundary conditions. In the case of impenetrable vertical boundaries, we confirm the transition from mainly conductive to mainly convective vertical heat transport below a critical thermal diffusivity. This transition is however much less dramatic when more natural outflow boundary conditions are applied. Similarly, the enhancement of magnetic activity in this case is less pronounced. Nevertheless, heating via turbulent dissipation determines the thermodynamic structure of accretion disks, and clearly affects the properties of the related dynamo. This effect may however have been overestimated in previous work, and a careful study of the role played by boundaries will be required.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global hydromagnetic simulations of a planet embedded in a dead zone: gap opening, gas accretion and formation of a protoplanetary jet

We present global hydrodynamic and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (i) The accretion flow in the Hill sphere is intrinsically 3D for hydrodynamic and MHD models. Net inflow toward the planet is dominated by high latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (ii) Gap opening magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere and a CPD whose structure displays high levels of variability. (iii) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch a highly collimated, one-sided jet is observed. (iv) The CPD's surface density $\sim30{\rm\,g\,cm^{-2}}$, small enough for significant ionization and turbulence to develop. (v) The accretion rate onto the planet in the MHD simulation reaches a steady value $8 \times 10^{-3} {\rm M_\oplus yr^{-1}}$, and is similar in the viscous hydrodynamic runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Toward a hybrid dynamo model for the Milky Way

(Abridged) Based on the rapidly increasing all-sky data of Faraday rotation measures and polarised synchrotron radiation, the Milky Way's magnetic field is now modelled with an unprecedented level of detail and complexity. We aim to complement this heuristic approach with a physically motivated, quantitative Galactic dynamo model -- a model that moreover allows for the evolution of the system as a whole, instead of just solving the induction equation for a fixed static disc. Building on the framework of mean-field magnetohydrodynamics and extending it to the realm of a hybrid evolution, we perform three-dimensional global simulations of the Galactic disc. Closure coefficients embodying the mean-field dynamo are calibrated against resolved box simulations of supernova-driven interstellar turbulence. The emerging dynamo solutions comprise a mixture of the dominant axisymmetric S0 mode, with even parity, and a subdominant A0 mode, with odd parity. Notably, such a superposition of modes creates a strong localised vertical field on one side of the Galactic disc. We moreover find significant radial pitch angles, which decay with radius -- explained by flaring of the disc. In accordance with previous work, magnetic instabilities appear to be restricted to the less-stirred outer Galactic disc. Their main effect is to create strong fields at large radii such that the radial scale length of the magnetic field increases from 4 kpc (for the case of a mean-field dynamo alone) to about 10 kpc in the hybrid models. There remain aspects (e.g., spiral arms, X-shaped halo fields, fluctuating fields) that are not captured by the current model and that will require further development towards a fully dynamical evolution. Nevertheless, the work presented demonstrates that a hybrid modelling of the Galactic dynamo is feasible and can serve as a foundation for future efforts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dead zones as safe-havens for planetesimals: influence of disc mass and external magnetic field

(Abridged) Planetesimals embedded in a protoplanetary disc are stirred by gravitational torques exerted by density fluctuations in the surrounding turbulence. In particular, planetesimals in a disc supporting fully developed magneto-rotational turbulence are readily excited to velocity dispersions above the threshold for catastrophic disruption, halting planet formation. We aim to examine the stirring of planetesimals lying instead in a magnetically-decoupled midplane dead zone, stirred only by spiral density waves propagating out of the disc's magnetically-coupled turbulent surface layers. We extend previous studies to include a wider range of disc models, and explore the effects of varying the disc column density and external magnetic field strength. [...] The strength of the stirring is found to be independent of the gas surface density, which is contrary to the increase with disc mass expected from a simple linear wave picture. The discrepancy arises from the shearing out of density waves as they propagate into the dead zone, resulting in density structures near the midplane that exert weaker stochastic torques on average. We provide a simple analytic fit to our numerically obtained torque amplitudes that accounts for this effect. The stirring on the other hand depends sensitively on the net vertical magnetic flux, up to a saturation level above which magnetic forces dominate in the turbulent layers. For the majority of our models, the equilibrium planetesimal velocity dispersions lie between the thresholds for disrupting strong and weak aggregates, suggesting that collision outcomes will depend on material properties. However, discs with relatively weak magnetic fields yield reduced stirring, and their dead zones provide safe-havens even for the weakest planetesimals against collisional destruction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Linear and nonlinear evolution of the vertical shear instability in accretion discs

(Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as a function of the power-law indices of the temperature or entropy, the thermal relaxation time of the fluid, and the viscosity. We present a stability analysis of the problem that we use to interpret the simulation results. We find that disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the thermodynamic response of the fluid is isothermal, or the thermal evolution time is comparable to or shorter than the local dynamical time scale. These discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear' linear instability. Development of the instability involves excitation of vertical breathing and corrugation modes in the disc, with the corrugation modes in particular being a feature of the nonlinear saturated state. Instability operates when the dimensionless disc kinematic viscosity nu < 10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates a quasi-turbulent flow, and the Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by the end of the simulation. The vertical shear instability in disc models which include realistic thermal physics has yet to be examined. Should it occur, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties, and observational appearance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the magnetic quenching of mean-field effects in supersonic interstellar turbulence

The emergence of large-scale magnetic fields observed in the diffuse interstellar medium is explained by a turbulent dynamo. The underlying transport coefficients have previously been extracted from numerical simulations. So far, this was restricted to the kinematic regime, but we aim to extend our analysis into the realm of dynamically important fields. This marks an important step on which derived mean-field models rely to explain observed equipartition strength fields. As in previous work, we diagnose turbulent transport coefficients by means of the test-field method. We derive quenching functions for the dynamo α effect, diamagnetic pumping, and turbulent diffusivity, which are compared with theoretical expectations. At late times, we observe the suppression of the vertical wind. Because this potentially affects the removal of small-scale magnetic helicity, new concerns arise about circumventing constraints imposed by the conservation of magnetic helicity at high magnetic Reynolds numbers. While present results cannot safely rule out this possibility, the issue only becomes important at late stages and is absent when the dynamo is quenched by the wind itself.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The role of star formation for the galactic dynamo

Magnetic field amplification by a fast dynamo is seen in local box simulations of SN-driven ISM turbulence, where the self-consistent emergence of large-scale fields agrees very well with its mean-field description. We accordingly derive scaling laws of the turbulent transport coef- ficients in dependence of the SN rate, density and rotation. These provide the input for global simulations of regular magnetic fields in galaxies within a mean-field MHD framework. Using a Kennicutt-Schmidt relation between the star formation (SF) rate and midplane density, we can reduce the number of free parameters in our global models. We consequently present dynamo models for different rotation curves and radial density distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the dynamics of planetesimals embedded in turbulent protoplanetary discs with dead zones

(abridged) Accretion in protoplanetary discs is thought to be driven by [...] turbulence via the magnetorotational instability (MRI). Recent work has shown that a planetesimal swarm embedded in a fully turbulent disc is subject to strong excitation of the velocity dispersion, leading to collisional destruction of bodies with radii R_p < 100 km. Significant diffusion of planetesimal semimajor axes also arises, leading to large-scale spreading of the planetesimal population throughout the inner regions of the protoplanetary disc, in apparent contradiction of constraints provided by the distribution of asteroids within the asteroid belt. In this paper, we examine the dynamics of planetesimals embedded in vertically stratified turbulent discs, with and without dead zones. Our main aims are to examine the turbulent excitation of the velocity dispersion, and the radial diffusion, of planetesimals in these discs. We employ three dimensional MHD simulations [...], along with an equilibrium chemistry model [...] We find that planetesimals in fully turbulent discs develop large random velocities that will lead to collisional destruction/erosion for bodies with sizes below 100 km, and undergo radial diffusion on a scale \sim 2.5 au over a 5 Myr disc life time. But planetesimals in a dead zone experience a much reduced excitation of their random velocities, and equilibrium velocity dispersions lie between the disruption thresholds for weak and strong aggregates for sizes R_p < 100 km. We also find that radial diffusion occurs over a much reduced length scale \sim 0.25 au over the disc life time, this being consistent with solar system constraints. We conclude that planetesimal growth via mutual collisions between smaller bodies cannot occur in a fully turbulent disc. By contrast, a dead zone may provide a safe haven in which km-sized planetesimals can avoid mutual destruction through collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Supernova-driven interstellar turbulence and the galactic dynamo

The fractal shape and multi-component nature of the interstellar medium together with its vast range of dynamical scales provides one of the great challenges in theoretical and numerical astrophysics. Here we will review recent progress in the direct modelling of interstellar hydromagnetic turbulence, focusing on the role of energy injection by supernova explosions. The implications for dynamo theory will be discussed in the context of the mean-field approach. Results obtained with the test field-method are confronted with analytical predictions and estimates from quasilinear theory. The simulation results enforce the classical understanding of a turbulent Galactic dynamo and, more importantly, yield new quantitative insights. The derived scaling relations enable confident global mean-field modelling.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A mean-field approach to the propagation of field patterns in stratified magneto rotational turbulence

Local shearing box simulations of stratified magneto rotational turbulence invariably exhibit cyclic field patterns which propagate away from the disc midplane. A common explanation for this is magnetic buoyancy. The recent analysis by Shi et al. however shows that the flow is buoyantly stable below one disc scale height H, necessitating an alternative explanation in this region. We here conduct and analyse direct numerical simulations to explain the observed behaviour by means of a mean-field description. Apart from the mean radial and azimuthal field, we monitor the small-scale current helicity, which we propose as a key indicator for saturation. Reconstructing the horizontally averaged field, we demonstrate that the problem can be reduced to a one-dimensional induction equation. By means of the so-called test field method, we then determine the underlying closure parameters. Our analysis shows that, apart from a possible direct MRI dynamo, two distinct indirect dynamo mechanisms operate in the disc. This resolves the issue of the "wrong" sign of the MRI dynamo effect. Finally, we use the obtained closure parameters to run a dynamically quenched dynamo model. This model approximately recovers the observed field patterns in the mean fields. Moreover, the model reproduces the prevailing parity and the distinct phase pattern in the small-scale current helicity. The latter property might open a potential route to understand the saturation of MRI induced turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

MRI channel flows in vertically-stratified models of accretion disks

Simulations of the magnetorotational instability (MRI) in 'unstratified' shearing boxes exhibit powerful coherent flows, whereby the fluid vertically splits into countermoving planar jets or `channels'. Channel flows correspond to certain axisymmetric linear MRI modes, and their preponderance follows from the remarkable fact that they are approximate nonlinear solutions of the MHD equations in the limit of weak magnetic fields. We show in this paper, analytically and with one-dimensional numerical simulations, that this property is also shared by certain axisymmetric MRI modes in vertically-stratified shearing boxes. These channel flows rapidly capture significant amounts of magnetic and kinetic energy, and thus are vulnerable to secondary shear instabilities. We examine these parasites in the vertically stratified context, and estimate the maximum amplitudes that channels attain before they are destroyed. These estimates suggest that a dominant channel flow will usually drive the disk's magnetic field to thermal strengths. The prominence of these flows and their destruction place enormous demands on simulations, but channels in their initial stages also offer a useful check on numerical codes. These benchmarks are especially valuable given the increasing interest in the saturation of the stratified MRI. Lastly we speculate on the potential connection between 'run-away' channel flows and outburst behaviour in protostellar and dwarf nova disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Supernova-driven Turbulence and Magnetic Field Amplification in Disk Galaxies

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. Analytical models based on the uncorrelated-ensemble approach predicted that any created field will be expelled from the disk before a significant amplification can occur. By means of direct simulations of supernova-driven turbulence, we demonstrate that this is not the case. Accounting for vertical stratification and galactic differential rotation, we find an exponential amplification of the mean field on timescales of 100Myr. The self-consistent numerical verification of such a "fast dynamo" is highly beneficial in explaining the observed strong magnetic fields in young galaxies. We, furthermore, highlight the importance of rotation in the generation of helicity by showing that a similar mechanism based on Cartesian shear does not lead to a sustained amplification of the mean magnetic field. This finding impressively confirms the classical picture of a dynamo based on cyclonic turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Magnetic field amplification by SN-driven interstellar turbulence

Within the interstellar medium, supernovae are thought to be the prevailing agents in driving turbulence. Until recently, their effects on magnetic field amplification in disk galaxies remained uncertain. Analytical models based on the uncorrelated-ensemble approach predicted that any created field would be expelled from the disk before it could be amplified significantly. By means of direct simulations of supernova-driven turbulence, we demonstrate that this is not the case. Accounting for galactic differential rotation and vertical stratification, we find an exponential amplification of the mean field on timescales of several hundred million years. We especially highlight the importance of rotation in the generation of helicity by showing that a similar mechanism based on Cartesian shear does not lead to a sustained amplification of the mean magnetic field.

0 citations0 reviewsNo code linkedOpen access