Source author record

Geoffroy Lesur

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

19 published item(s)

preprint2026arXiv

Magnetic winds in resistive compact binary discs

Dwarf Novae and low-mass X-ray binaries are eruptive binary systems comprised of a Roche-lobe overflowing solar-type star and an accreting compact object. Their recurrence time can be explained by a low-accreting phase, the quiescence, during which the angular momentum transport parameter is inferred to be $α\approx 0.01$ by the Disc Instability Model. Non-magnetics mechanisms, such as spiral wave transport, only achieve angular momentum transport an order of magnitude too low, at best, because these discs are so thin in quiescence. During this phase, the Magneto-rotational Instability is known to be suppressed by the increased resistivity of the weakly ionised plasma. Studying these thin magnetised discs is a numerical challenge because of the wide range of scales to be resolved. Thanks to the new GPU-accelerated code Idefix, we produce global 3D MHD simulations of a very thin disc $(H/R = 0.01)$ for the first time. We explore the possibility that an MHD wind arises and increases accretion in low magnetic Reynolds number $(\mathrm{Rm}\approx100)$ and realistic plasma parameter ($β\approx 1000$) regimes. We observe that the MRI is only quenched in the resistive disc bulk but survives in the disc atmosphere. This drives strong accretion and wind launching. We quantify the efficiency of the resulting wind and measure its global effect on the disc. We explore the effect of the initial disc magnetisation and compare the accretion/ejection regime with and without resistivity.

preprint2020arXiv

Magnetic field transport in compact binaries

Dwarf novae (DNe) and low mass X-ray binaries (LMXBs) show eruptions that are thought to be due to a thermal-viscous instability in their accretion disk. These eruptions provide constraints on angular momentum transport mechanisms. We explore the idea that angular momentum transport could be controlled by the dynamical evolution of the large scale magnetic field. We study the impact of different prescriptions for the magnetic field evolution on the dynamics of the disk. This is a first step in confronting the theory of magnetic field transport with observations. We develop a version of the disk instability model that evolves the density, the temperature and the large scale vertical magnetic flux together. We take into account the accretion driven by turbulence or by a magnetized outflow. To evolve the magnetic flux, we use a toy model with physically motivated prescriptions depending mainly on the local magnetization. We find that allowing magnetic flux to be advected inwards provides the best agreement with DNe lightcurves. This leads to a hybrid configuration with an inner magnetized disk, driven by angular momentum losses to an MHD outflow, sharply transiting to an outer weakly-magnetized turbulent disk, where the eruptions are triggered. The dynamical impact is equivalent to truncating a viscous disk so that it does not extend down to the compact object, with the truncation radius dependent on the magnetic flux and evolving as $\dot{M}^{-2/3}$. Models of DNe and LMXBs lightcurves typically require the outer, viscous disk to be truncated in order to match observations. There is no generic explanation for this truncation. We propose that it is a natural outcome of the presence of large-scale magnetic fields in both DNe and LMXBs, the magnetic flux accumulating towards the center to produce a magnetized disk with a fast accretion timescale.

preprint2016arXiv

On the survival of zombie vortices in protoplanetary discs

Recently it has been proposed that the zombie vortex instability (ZVI) could precipitate hydrodynamical activity and angular momentum transport in unmagnetised regions of protoplanetary discs, also known as "dead zones". In this letter we scrutinise, with high resolution 3D spectral simulations, the onset and survival of this instability in the presence of viscous and thermal physics. First, we find that the ZVI is strongly dependent on the nature of the viscous operator. Although the ZVI is easily obtained with hyper-diffusion, it is difficult to sustain with physical (second order) diffusion operators up to Reynolds numbers as high as 10^7. This sensitivity is probably due to the ZVI's reliance on critical layers, whose characteristic lengthscale, structure, and dynamics are controlled by viscous diffusion. Second, we observe that the ZVI is sensitive to radiative processes, and indeed only operates when the Peclet number is greater than a critical value ~10^4, or when the cooling time is longer than ~10 Omega^{-1}. As a consequence, the ZVI struggles to appear at R>0.3 AU in standard 0.01 M_sun T Tauri disc models, though younger more massive disks provide a more hospitable environment. Together these results question the prevalence of the ZVI in protoplanetary discs.

preprint2016arXiv

Self-organisation in protoplanetary disks: global, non-stratified Hall-MHD simulations

Recent observations revealed organised structures in protoplanetary disks, such as axisymmetric rings or horseshoe concen- trations evocative of large-scale vortices. These structures are often interpreted as the result of planet-disc interactions. However, these disks are also known to be unstable to the magneto-rotational instability (MRI) which is believed to be one of the dominant angular momentum transport mechanism in these objects. It is therefore natural to ask if the MRI itself could produce these structures without invoking planets. The nonlinear evolution of the MRI is strongly affected by the low ionisation fraction in protoplanetary disks. The Hall effect in particular, which is dominant in dense and weakly ionised parts of these objects, has been shown to spontaneously drive self- organising flows in shearing box simulations. Here, we investigate the behaviour of global MRI-unstable disc models dominated by the Hall effect and characterise their dynamics. We perform 3D unstratified Hall-MHD simulations of keplerian disks for a broad range of Hall, ohmic and ambipolar Elsasser numbers. We confirm the transition from a turbulent to an organised state as the intensity of the Hall effect is increased. We observe the formation of zonal flows, their number depending on the available magnetic flux and on the intensity of the Hall effect. For intermediate Hall intensity, the flow self-organises into long-lived magnetised vortices. Neither the addition of a toroidal field nor ohmic or ambipolar diffusion drastically change this picture in the range of parameters we have explored. The ability of these structures to trap dust particles in this configuration is demonstrated. We conclude that Hall-MRI driven organisation is a plausible scenario which could explain some of the structures found in recent observations.

preprint2016arXiv

Spiral-driven accretion in protoplanetary discs - II Self-similar solutions

Accretion discs are ubiquitous in the universe and it is a crucial issue to understand how angular momentum and mass are being radially transported in these objects. Here, we study the role played by non-linear spiral patterns within hydrodynamical and non self-gravitating accretion disc assuming that external disturbances such as infall onto the disc may trigger them. To do so, we computed self-similar solutions that describe discs in which a spiral wave propagates. Such solutions present both shocks and critical sonic points that we carefully analyze. For all allowed temperatures and for several spiral shocks, we calculated the wave structure. In particular we inferred the angle of the spiral patern, the stress it exerts on the disc as well as the associated flux of mass and angular momentum as a function of temperature. We quantified the rate of angular momentum transport by means of the dimensionless $α$ parameter. For the thickest disc we considered (corresponding to $h/r$ values of about 1/3), we found values of $α$ as high as $0.1$, and scaling with the temperature $T$ such that $α\propto T^{3/2} \propto (h/r)^3$. The spiral angle scales with the temperature as $\arctan(r/h)$. The existence of these solutions suggests that perturbations occurring at disc outer boundaries, such as for example perturbations due to infall motions, can propagate deep inside the disc and therefore should not be ignored, even when considering small radii.

preprint2015arXiv

Magnetically driven accretion in protoplanetary discs

We characterize magnetically driven accretion at radii between 1 au and 100 au in protoplanetary discs, using a series of local non-ideal magnetohydrodynamic (MHD) simulations. The simulations assume a Minimum Mass Solar Nebula (MMSN) disc that is threaded by a net vertical magnetic field of specified strength. Confirming previous results, we find that the Hall effect has only a modest impact on accretion at 30 au, and essentially none at 100 au. At 1-10 au the Hall effect introduces a pronounced bi-modality in the accretion process, with vertical magnetic fields aligned to the disc rotation supporting a strong laminar Maxwell stress that is absent if the field is anti-aligned. In the anti-aligned case, we instead find evidence for bursts of turbulent stress at 5-10 au, which we tentatively identify with the non-axisymmetric Hall-shear instability. The presence or absence of these bursts depends upon the details of the adopted chemical model, which suggests that appreciable regions of actual protoplanetary discs might lie close to the borderline between laminar and turbulent behaviour. Given the number of important control parameters that have already been identified in MHD models, quantitative predictions for disc structure in terms of only radius and accretion rate appear to be difficult. Instead, we identify robust qualitative tests of magnetically driven accretion. These include the presence of turbulence in the outer disc, independent of the orientation of the vertical magnetic fields, and a Hall-mediated bi-modality in turbulent properties extending from the region of thermal ionization to 10 au.

preprint2015arXiv

On the Nature of Magnetic Turbulence in Rotating, Shearing Flows

The local properties of turbulence driven by the magnetorotational instability (MRI) in rotating, shearing flows are studied in the framework of a shearing-box model. Based on numerical simulations, we propose that the MRI-driven turbulence comprises two components: the large-scale shear-aligned strong magnetic field and the small-scale fluctuations resembling magnetohydrodynamic (MHD) turbulence. The energy spectrum of the large-scale component is close to $k^{-2}$, whereas the spectrum of the small-scale component agrees with the spectrum of strong MHD turbulence $k^{-3/2}$. While the spectrum of the fluctuations is universal, the outer-scale characteristics of the turbulence are not; they depend on the parameters of the system, such as the net magnetic flux. However, there is remarkable universality among the allowed turbulent states -- their intensity $v_0$ and their outer scale $λ_0$ satisfy the balance condition $v_0/λ_0\sim \mathrm dΩ/\mathrm d\ln r$, where $\mathrm dΩ/\mathrm d\ln r$ is the local orbital shearing rate of the flow. Finally, we find no sustained dynamo action in the $\mathrm{Pm}=1$ zero net-flux case for Reynolds numbers as high as $45\,000$, casting doubts on the existence of an MRI dynamo in the $\mathrm{Pm}\leq 1$ regime.

preprint2015arXiv

Spiral-driven accretion in protoplanetary discs - I. 2D models

We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 < alpha < 1e-2 for an inflow accretion rate Mout>1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.

preprint2014arXiv

Thanatology in Protoplanetary Discs: the combined influence of Ohmic, Hall, and ambipolar diffusion on dead zones

Protoplanetary discs are poorly ionised due to their low temperatures and high column densities, and are therefore subject to three "non-ideal" magnetohydrodynamic effects: Ohmic dissipation, ambipolar diffusion, and the Hall effect. The existence of magnetically driven turbulence in these discs has been a central question since the discovery of the magnetorotational instability. Early models considered Ohmic diffusion only and led to a scenario of layered accretion, in which a magnetically "dead" zone in the disc midplane is embedded within magnetically "active" surface layers at distances ~1-10 au from the central protostellar object. Recent work has suggested that a combination of Ohmic dissipation and ambipolar diffusion can render both the midplane and surface layers of the disc inactive and that torques due to magnetically driven outflows are required to explain the observed accretion rates. We reassess this picture by performing three-dimensional numerical simulations that include, for the first time, all three non-ideal MHD effects. We find that the Hall effect can generically "revive" dead zones by producing a dominant azimuthal magnetic field and a large-scale Maxwell stress throughout the midplane, provided the angular velocity and magnetic field satisfy Omega.B > 0. The attendant large magnetic pressure modifies the vertical density profile and substantially increases the disc scale height beyond its hydrostatic value. Outflows are produced, but are not necessary to explain accretion rates <10^{-7} Msun/yr. The flow in the disc midplane is essentially laminar, suggesting that dust sedimentation may be efficient. These results demonstrate that, if the MRI is relevant for driving mass accretion in protoplanetary discs, one must include the Hall effect to obtain even qualitatively correct results.

preprint2013arXiv

Local outflows from turbulent accretion disks

The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between MHD turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio beta0 of the thermal and magnetic pressures in the disk midplane equals 10^4. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows.

preprint2013arXiv

Magnetic self-organisation in Hall-dominated magnetorotational turbulence

The magnetorotational instability (MRI) is the most promising mechanism by which angular momentum is efficiently transported outwards in astrophysical discs. However, its application to protoplanetary discs remains problematic. These discs are so poorly ionised that they may not support magnetorotational turbulence in regions referred to as `dead zones'. It has recently been suggested that the Hall effect, a non-ideal magnetohydrodynamic (MHD) effect, could revive these dead zones by enhancing the magnetically active column density by an order of magnitude or more. We investigate this idea by performing local, three-dimensional, resistive Hall-MHD simulations of the MRI in situations where the Hall effect dominates over Ohmic dissipation. As expected from linear stability analysis, we find an exponentially growing instability in regimes otherwise linearly stable in resistive MHD. However, instead of vigorous and sustained magnetorotational turbulence, we find that the MRI saturates by producing large-scale, long-lived, axisymmetric structures in the magnetic and velocity fields. We refer to these structures as zonal fields and zonal flows, respectively. Their emergence causes a steep reduction in turbulent transport by at least two orders of magnitude from extrapolations based upon resistive MHD, a result that calls into question contemporary models of layered accretion. We construct a rigorous mean-field theory to explain this new behaviour and to predict when it should occur. Implications for protoplanetary disc structure and evolution, as well as for theories of planet formation, are briefly discussed.

preprint2012arXiv

Colliding wind binaries and gamma-ray binaries : relativistic version of the RAMSES code

Gamma-ray binaries are colliding wind binaries (CWB) composed of a massive star a non-accreting pulsar with a highly relativistic wind. Particle acceleration at the shocks results in emission going from extended radio emission to the gamma-ray band. The interaction region is expected to show common features with stellar CWB. Performing numerical simulations with the hydrodynamical code RAMSES, we focus on their structure and stability and find that the Kelvin-Helmholtz instability (KHI) can lead to important mixing between the winds and destroy the large scale spiral structure. To investigate the impact of the relativistic nature of the pulsar wind, we extend RAMSES to relativistic hydrodynamics (RHD). Preliminary simulations of the interaction between a pulsar wind and a stellar wind show important similarities with stellar colliding winds with small relativistic corrections.

preprint2012arXiv

Impact of orbital motion on the structure and stability of adiabatic shocks in colliding wind binaries

The collision of winds from massive stars in binaries results in the formation of a double-shock structure with observed signatures from radio to X-rays. We study the structure and stability of the colliding wind region as it turns into a spiral due to orbital motion. We focus on adiabatic winds, where mixing between the two winds is expected to be restricted to the Kelvin-Helmholtz instability (KHI). Mixing of the Wolf-Rayet wind with hydrogen-rich material is important for dust formation in pinwheel nebulae such as WR 104, where the spiral structure has been resolved in infrared. We use the hydrodynamical code RAMSES with an adaptive grid. A wide range of binary systems with different wind velocities and mass loss rates are studied with 2D simulations. A specific 3D simulation is performed to model WR 104. Orbital motion leads to the formation of two distinct spiral arms where the KHI develops differently. We find that the spiral structure is destroyed when there is a large velocity gradient between the winds, unless the collimated wind is much faster. We argue that the KHI plays a major role in maintaining or not the structure. We discuss the consequences for various colliding wind binaries. When stable, there is no straightforward relationship between the spatial step of the spiral, the wind velocities, and the orbital period. Our 3D simulation of WR 104 indicates that the colder, well-mixed trailing arm has more favourable conditions for dust formation than the leading arm. The single-arm infrared spiral follows more closely the mixing map than the density map, suggesting the dust-to-gas ratio may vary between the leading and trailing density spirals. However, the density is much lower than what dust formation models require. Including radiative cooling would lead to higher densities, and also to thin shell instabilities whose impact on the large structure remains unknown.

preprint2012arXiv

Modelling colliding wind binaries with RAMSES, extension to special relativity

We present high resolution simulations with RAMSES of supersonic colliding stellar winds. The collision results in a double shock structure which is subject to different instabilities. The Kelvin-Helmholtz instability (KHI) introduces some mixing and variability. For isothermal winds, the Non-linear Thin Shell Instability violently affects the interaction region. Properly modelling these instabilities requires a high enough resolution and an adapted numerical method, especially when one of the winds strongly dominates the other one. At large scale, orbital motion is expected to turn the shocked zone into a spiral but we find that in some configurations the KHI may disrupt the spiral. A colliding wind structure is also expected in gamma-ray binaries composed of a massive star and a young pulsar which emits a highly relativistic wind. Numerical simulations are necessary to understand the geometry of such systems and should take into account the relativistic nature of the pulsar wind. We implemented a second order Godunov method to solve the equations of relativistic hydrodynamics with RAMSES, including the possibility for Adaptive Mesh Refinement. After a brief overview of our numerical implementation, we will present preliminary simulation that contrast the structure of the flow in gamma-ray and stellar binaries.

preprint2012arXiv

On the interaction between tides and convection

We study the interaction between tides and convection in astrophysical bodies by analysing the effect of a homogeneous oscillatory shear on a fluid flow. This model can be taken to represent the interaction between a large-scale periodic tidal deformation and a smaller-scale convective motion. We first consider analytically the limit in which the shear is of low amplitude and the oscillation period is short compared to the timescales of the unperturbed flow. In this limit there is a viscoelastic response and we obtain expressions for the effective elastic modulus and viscosity coefficient. The effective viscosity is inversely proportional to the square of the oscillation frequency, with a coefficient that can be positive, negative or zero depending on the properties of the unperturbed flow. We also carry out direct numerical simulations of Boussinesq convection in an oscillatory shearing box and measure the time-dependent Reynolds stress. The results indicate that the effective viscosity of turbulent convection falls rapidly as the oscillation frequency is increased, attaining small negative values in the cases we have examined, although significant uncertainties remain because of the turbulent noise. We discuss the implications of this analysis for astrophysical tides.

preprint2012arXiv

The magnetorotational instability as a jet launching mechanism

Magnetorotational turbulence and magnetically driven disc winds are often considered as separate processes. However, realistic astrophysical discs are expected to be subject to both effects, although possibly at different times and locations. We investigate here the potential link between these two phenomena using a mixed numerical and analytical approach. We show in particular that large-scale MRI modes which dominate strongly magnetised discs (plasma beta~10) naturally produce magnetically driven outflows in the nonlinear regime. We show that these outflows share many similarities with local and global disc wind solutions found in the literature. We also investigate the 3D stability of these outflows and show that they are unstable on dynamical timescales. The implications of these results for the transition between a jet-emitting disc and a standard "viscous" disc are discussed.

preprint2010arXiv

Nonlinear energy transfers in accretion discs MRI turbulence. I-Net vertical field case

The magnetorotational instability (MRI) is believed to be responsible for most of the angular momentum transport in accretion discs. However, molecular dissipation processes may drastically change the efficiency of MRI turbulence in realistic astrophysical situations. The physical origin of this dependency is still poorly understood as linear and quasi linear theories fail to explain it. In this paper, we look for the link between molecular dissipation processes and MRI transport of angular momentum in non stratified shearing box simulations including a mean vertical field. We show that magnetic helicity is unimportant in the model we consider. We perform a spectral analysis on the simulations tracking energy exchanges in spectral space when turbulence is fully developed. We find that the energy exchanges are essentially direct (from large to small scale) whereas some non linear interactions appear to be non local in spectral space. We speculate that these non local interactions are responsible for the correlation between turbulent transport and molecular dissipation. We argue that this correlation should then disappear when a significant scale separation is achieved and we discuss several methods by which one can test this hypothesis.

preprint2010arXiv

Supertransient magnetohydrodynamic turbulence in Keplerian shear flows

A subcritical transition to turbulence in magnetized Keplerian shear flows is investigated using a statistical approach. Three-dimensional numerical simulations of the shearing box equations with zero net magnetic flux are employed to determine the transition from decaying to sustained turbulence as a function of the magnetic Reynolds number Rm. The results reveal no clear transition to sustained turbulence as the average lifetime of the transients grows as an exponential function of Rm, in accordance with a type-II supertransient law.

preprint2010arXiv

The Validity of the Super-Particle Approximation during Planetesimal Formation

The formation mechanism of planetesimals in protoplanetary discs is hotly debated. Currently, the favoured model involves the accumulation of meter-sized objects within a turbulent disc, followed by a phase of gravitational instability. At best one can simulate a few million particles numerically as opposed to the several trillion meter-sized particles expected in a real protoplanetary disc. Therefore, single particles are often used as super-particles to represent a distribution of many smaller particles. It is assumed that small scale phenomena do not play a role and particle collisions are not modeled. The super-particle approximation can only be valid in a collisionless or strongly collisional system, however, in many recent numerical simulations this is not the case. In this work we present new results from numerical simulations of planetesimal formation via gravitational instability. A scaled system is studied that does not require the use of super-particles. We find that the scaled particles can be used to model the initial phases of clumping if the properties of the scaled particles are chosen such that all important timescales in the system are equivalent to what is expected in a real protoplanetary disc. Constraints are given for the number of particles needed in order to achieve numerical convergence. We compare this new method to the standard super-particle approach. We find that the super-particle approach produces unreliable results that depend on artifacts such as the gravitational softening in both the requirement for gravitational collapse and the resulting clump statistics. Our results show that short range interactions (collisions) have to be modelled properly.