BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jerome Guilet

Jerome Guilet contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-ph.EPastro-ph.GA

Trust snapshot

Quick read

Trust 21 - Emerging
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2015arXiv

Neutrino viscosity and drag: impact on the magnetorotational instability in protoneutron stars

The magnetorotational instability (MRI) is a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. The diffusion of neutrinos trapped in the PNS induces a transport of momentum, which can be modelled as a viscosity on length-scales longer than the neutrino mean free path. This neutrino-viscosity can slow down the growth of MRI modes to such an extent that a minimum initial magnetic field strength of $\gtrsim 10^{12}\, {\rm G}$ is needed for the MRI to grow on a sufficiently short time-scale to potentially affect the explosion. It is uncertain whether the magnetic field of fast rotating progenitor cores is strong enough to yield such an initial magnetic field in PNS. At MRI wavelengths shorter than the neutrino mean free path, on the other hand, neutrino radiation does not act as a viscosity but rather induces a drag on the velocity with a damping rate independent of the wavelength. We perform a linear analysis of the MRI in this regime, and apply our analytical results to the PNS structure from a one-dimensional numerical simulation. We show that in the outer layers of the PNS, the MRI can grow from weak magnetic fields at wavelengths shorter than the neutrino mean free path, while deeper in the PNS MRI growth takes place in the viscous regime and requires a minimum magnetic field strength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical simulations of the magnetorotational instability in protoneutron stars: I. Influence of buoyancy

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The stress-pressure relationship in simulations of MRI-induced turbulence

We determine how MRI-turbulent stresses depend on gas pressure via a suite of unstratified shearing box simulations. Earlier numerical work reported only a very weak dependence at best, results that call into question the canonical alpha-disk model and the thermal stability results that follow from it. Our simulations, in contrast, exhibit a stronger relationship, and show that previous work was box-size limited: turbulent `eddies' were artificially restricted by the numerical domain rather than by the scale height. Zero-net-flux runs without physical diffusion coefficients yield a stress proportional to $P^{0.5}$, where P is pressure. The stresses are also proportional to the grid length and hence remain numerically unconverged. The same runs with physical diffusivities, however, give a result closer to an alpha-disk: the stress is proportional to $P^{0.9}$. Net-flux simulations without explicit diffusion exhibit stresses proportional to $P^{0.5}$, but stronger imposed fields weaken this correlation. In summary, compressibility is important for the saturation of the MRI, but the exact stress-pressure relationship is difficult to ascertain in local simulations because of numerical convergence issues and the influence of any imposed flux. As a consequence, the interpretation of thermal stability behaviour in local simulations is a problematic enterprise.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global evolution of the magnetic field in a thin disc and its consequences for protoplanetary systems

The strength and structure of the large-scale magnetic field in protoplanetary discs are still unknown, although they could have important consequences for the dynamics and evolution of the disc. Using a mean-field approach in which we model the effects of turbulence through enhanced diffusion coefficients, we study the time-evolution of the large-scale poloidal magnetic field in a global model of a thin accretion disc, with particular attention to protoplanetary discs. With the transport coefficients usually assumed, the magnetic field strength does not significantly increase radially inwards, leading to a relatively weak magnetic field in the inner part of the disc. We show that with more realistic transport coefficients that take into account the vertical structure of the disc and the back-reaction of the magnetic field on the flow as obtained by Guilet & Ogilvie (2012), the magnetic field can significantly increase radially inwards. The magnetic-field profile adjusts to reach an equilibrium value of the plasma $β$ parameter (the ratio of midplane thermal pressure to magnetic pressure) in the inner part of the disc. This value of $β$ depends strongly on the aspect ratio of the disc and on the turbulent magnetic Prandtl number, and lies in the range $10^4-10^7$ for protoplanetary discs. Such a magnetic field is expected to affect significantly the dynamics of protoplanetary discs by increasing the strength of MHD turbulence and launching an outflow. We discuss the implications of our results for the evolution of protoplanetary discs and for the formation of powerful jets as observed in T-Tauri star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Type I planet migration in weakly magnetised laminar discs

[Abridged] The migration of low mass planets has been studied in hydrodynamical disc models for more than three decades, but the impact of a magnetic field in the protoplanetary disc is less known. When the disc's magnetic field is strong enough to prevent horseshoe motion, the corotation torque is replaced by a torque arising from magnetic resonances. For weak enough magnetic fields, horseshoe motion and a corotation torque exist, and recent turbulent MHD simulations have reported the existence of a new component of the corotation torque in the presence of a mean toroidal field. The aim of this paper is to investigate the physical origin and the properties of this new corotation torque. We performed MHD simulations of a low mass planet embedded in a 2D laminar disc threaded by a weak toroidal magnetic field, with the effects of turbulence modelled by a viscosity and a resistivity. We confirm that the interaction between the magnetic field and the horseshoe motion results in an additional corotation torque on the planet, which we dub the MHD torque excess. It is caused by the accumulation of the magnetic field along the downstream separatrices of the planet's horseshoe region, which gives rise to an azimuthally asymmetric underdense region at that location. The properties of the MHD torque excess are characterised by varying the slope of the density, temperature and magnetic field profiles, as well as the diffusion coefficients and the strength of the magnetic field. The sign of the MHD torque excess depends on the density and temperature gradients only, and is positive for profiles expected in protoplanetary discs. Its magnitude is in turn mainly determined by the strength of the magnetic field and the turbulent resistivity. The MHD torque excess can be strong enough to reverse migration, even when the magnetic pressure is less than one percent of the thermal pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transport of magnetic flux and the vertical structure of accretion discs: I. Uniform diffusion coefficients

The evolution of a large-scale poloidal magnetic field in accretion discs is an important problem because of its role in the launching of jets and winds and in determining the intensity of turbulence. In this paper, we develop a formalism to calculate the transport magnetic flux in a thin accretion disc, thus determining its evolution on a viscous/resistive timescale. The governing equations are derived by performing an asymptotic expansion in the limit of a thin disc, in the regime where the magnetic field is dominated by its vertical component. Turbulent viscosity and resistivity are included, with an arbitrary vertical profile that can be adjusted to mimic the vertical structure of the turbulence. At a given radius and time, the rates of transport of mass and magnetic flux are determined by a one-dimensional problem in the vertical direction, in which the radial gradients of various quantities appear as source terms. We solve this problem to obtain the transport rates and the vertical structure of the disc. The present paper is then restricted to the idealised case of uniform diffusion coefficients, while a companion paper will study more realistic vertical profiles of these coefficients. We show the advection of weak magnetic fields to be significantly faster than the advection of mass, contrary to what a crude vertical averaging might suggest. This results from the larger radial velocities away from the mid-plane, which barely affect the mass accretion owing to the low density in these regions but do affect the advection of magnetic flux. Possible consequences of this larger accretion velocity include a potentially interesting time-dependence with the magnetic flux distribution evolving faster than the mass distribution. If the disc is not too thin, this fast advection may also partially solve the long-standing problem of too efficient diffusion of an inclined magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics of an Alfven surface in core collapse supernovae

We investigate the dynamics of an Alfven surface (where the Alfven speed equals the advection velocity) in the context of core collapse supernovae during the phase of accretion on the proto-neutron star. Such a surface should exist even for weak magnetic fields because the advection velocity decreases to zero at the center of the collapsing core. In this decelerated flow, Alfven waves created by the standing accretion shock instability (SASI) or convection accumulate and amplify while approaching the Alfven surface. We study this amplification using one dimensional MHD simulations with explicit physical dissipation. In the linear regime, the amplification continues until the Alfven wavelength becomes as small as the dissipative scale. A pressure feedback that increases the pressure in the upstream flow is created via a non linear coupling. We derive analytic formulae for the maximum amplification and the non linear coupling and check them with numerical simulations to a very good accuracy. We also characterize the non linear saturation of this amplification when compression effects become important, leading to either a change of the velocity gradient, or a steepening of the Alfven wave. Applying these results to core collapse supernovae shows that the amplification can be fast enough to affect the dynamics, if the magnetic field is strong enough for the Alfven surface to lie in the region of strong velocity gradient just above the neutrinosphere. This requires the presence of a strong magnetic field in the progenitor star, which would correspond to the formation of a magnetar under the assumption of magnetic flux conservation. An extrapolation of our analytic formula (taking into account the nonlinear saturation) suggests that the Alfven wave could reach an amplitude of B ~ 10^15 G, and that the pressure feedback could significantly contribute to the pressure below the shock.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the linear growth mechanism driving the stationary accretion shock instability

During stellar core collapse, which eventually leads to a supernovae explosion, the stalled shock is unstable due to the standing accretion shock instability (SASI). This instability induces large-scale non spherical oscillations of the shock, which have crucial consequences on the dynamics and the geometry of the explosion. While the existence of this instability has been firmly established, its physical origin remains somewhat uncertain. Two mechanisms have indeed been proposed to explain its linear growth. The first is an advective-acoustic cycle, where the instability results from the interplay between advected perturbations (entropy and vorticity) and an acoustic wave. The second mechanism is purely acoustic and assumes that the shock is able to amplify trapped acoustic waves. Several arguments favouring the advective-acoustic cycle have already been proposed, however none was entirely conclusive for realistic flow parameters. In this article we give two new arguments which unambiguously show that the instability is not purely acoustic, and should be attributed to the advective-acoustic cycle. First, we extract a radial propagation timescale by comparing the frequencies of several unstable harmonics that differ only by their radial structure. The extracted time matches the advective-acoustic time but strongly disagrees with a purely acoustic interpretation. Second, we present a method to compute purely acoustic modes, by artificially removing advected perturbations below the shock. All these purely acoustic modes are found to be stable, showing that the advected wave is essential to the instability mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The saturation of SASI by parasitic instabilities

The Standing Accretion Shock Instability (SASI) is commonly believed to be responsible for large amplitude dipolar oscillations of the stalled shock during core collapse, potentially leading to an asymmetric supernovae explosion. The degree of asymmetry depends on the amplitude of SASI, which nonlinear saturation mechanism has never been elucidated. We investigate the role of parasitic instabilities as a possible cause of nonlinear SASI saturation. As the shock oscillations create both vorticity and entropy gradients, we show that both Kelvin-Helmholtz and Rayleigh-Taylor types of instabilities are able to grow on a SASI mode if its amplitude is large enough. We obtain simple estimates of their growth rates, taking into account the effects of advection and entropy stratification. In the context of the advective-acoustic cycle, we use numerical simulations to demonstrate how the acoustic feedback can be decreased if a parasitic instability distorts the advected structure. The amplitude of the shock deformation is estimated analytically in this scenario. When applied to the set up of Fernandez & Thompson (2009a), this saturation mechanism is able to explain the dramatic decrease of the SASI power when both the nuclear dissociation energy and the cooling rate are varied. Our results open new perspectives for anticipating the effect, on the SASI amplitude, of the physical ingredients involved in the modeling of the collapsing star.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Toward a magnetohydrodynamic theory of the stationary accretion shock: toy model of the advective-acoustic cycle in a magnetized flow

The effect of a magnetic field on the linear phase of the advective-acoustic instability is investigated, as a first step toward a magnetohydrodynamic (MHD) theory of the stationary accretion shock instability taking place during stellar core collapse. We study a toy model where the flow behind a planar stationary accretion shock is adiabatically decelerated by an external potential. Two magnetic field geometries are considered: parallel or perpendicular to the shock. The entropy-vorticity wave, which is simply advected in the unmagnetized limit, separates into five different waves: the entropy perturbations are advected, while the vorticity can propagate along the field lines through two Alfven waves and two slow magnetosonic waves. The two cycles existing in the unmagnetized limit, advective-acoustic and purely acoustic, are replaced by up to six distinct MHD cycles. The phase differences among the cycles play an important role in determining the total cycle efficiency and hence the growth rate. Oscillations in the growth rate as a function of the magnetic field strength are due to this varying phase shift. A vertical magnetic field hardly affects the cycle efficiency in the regime of super-Alfvenic accretion that is considered. In contrast, we find that a horizontal magnetic field strongly increases the efficiencies of the vorticity cycles that bend the field lines, resulting in a significant increase of the growth rate if the different cycles are in phase. These magnetic effects are significant for large-scale modes if the Alfven velocity is a sizable fraction of the flow velocity.

0 citations0 reviewsNo code linkedOpen access