Showing 1–50 of 52 results for author: Hollerbach, R

Tidal dissipation in stably stratified and semi-convective regions of rotating giant planets: incorporating Coriolis forces

Authors: Christina M. Pontin, Adrian J. Barker, Rainer Hollerbach

Abstract: We study how stably stratified or semi-convective layers alter tidal dissipation rates associated with the generation of inertial, gravito-inertial, interfacial and surface gravity waves in rotating giant planets. We explore scenarios in which stable (non-convective) layers contribute to the high rates of tidal dissipation observed for Jupiter and Saturn in our solar system. Our model is an ideali… ▽ More We study how stably stratified or semi-convective layers alter tidal dissipation rates associated with the generation of inertial, gravito-inertial, interfacial and surface gravity waves in rotating giant planets. We explore scenarios in which stable (non-convective) layers contribute to the high rates of tidal dissipation observed for Jupiter and Saturn in our solar system. Our model is an idealised spherical Boussinesq system incorporating Coriolis forces to study effects of stable stratification and semi-convective layers on tidal dissipation. Our detailed numerical calculations consider realistic tidal forcing and compute the resulting viscous and thermal dissipation rates. The presence of an extended stably stratified fluid core significantly enhances tidal wave excitation of both inertial waves (due to rotation) in the convective envelope and gravito-inertial waves in the dilute core. We show that a sufficiently strongly stratified fluid core enhances inertial wave dissipation in a convective envelope much like a solid core does. We demonstrate that efficient tidal dissipation rates (and associated tidal quality factors $Q'$) -- sufficient to explain the observed migration rates of Saturn's moons -- are predicted at the frequencies of the orbiting moons due to the excitation of inertial or gravito-inertial waves in our models with stable layers (without requiring resonance-locking). Stable layers could also be important for tidal evolution of hot and warm Jupiters, and hot Neptunes, providing efficient tidal circularisation rates. Future work should study more sophisticated planetary models that also account for magnetism and differential rotation, as well as the interaction of inertial waves with turbulent convection. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 26 pages, 18 figures, 1 table, accepted for publication in ApJ (6th Nov 2023)

Three-dimensional magneto-thermal evolution of off-centred dipole magnetic field configurations in neutron stars

Authors: Andrei P. Igoshev, Rainer Hollerbach, Toby Wood

Abstract: Off-centred dipole configurations have been suggested to explain different phenomena in neutron stars, such as natal kicks, irregularities in polarisation of radio pulsars and properties of X-ray emission from millisecond pulsars. Here for the first time we model magneto-thermal evolution of neutron stars with crust-confined magnetic fields and off-centred dipole moments. We find that the dipole s… ▽ More Off-centred dipole configurations have been suggested to explain different phenomena in neutron stars, such as natal kicks, irregularities in polarisation of radio pulsars and properties of X-ray emission from millisecond pulsars. Here for the first time we model magneto-thermal evolution of neutron stars with crust-confined magnetic fields and off-centred dipole moments. We find that the dipole shift decays with time if the initial configuration has no toroidal magnetic field. The decay timescale is inversely proportional to magnetic field. The octupole moment decreases much faster than the quadrupole. Alternatively, if the initial condition includes strong dipolar toroidal magnetic field, the external poloidal magnetic field evolves from centred dipole to off-centred dipole. The surface thermal maps are very different for configurations with weak $B = 10^{13}$ G and strong $B = 10^{14}$ G magnetic fields. In the former case, the magnetic equator is cold while in the latter case it is hot. We model lightcurves and spectra of our magneto-thermal configurations. We found that in the case of cold equator, the pulsed fraction is small (below a few percent in most cases) and spectra are well described with a single blackbody. Under the same conditions models with stronger magnetic fields produce lightcurves with pulsed fraction of tens of percent. Their spectra are significantly better described with two blackbodies. Overall, the magnetic field strength has a more significant effect on bulk thermal emission of neutron stars than does the field geometry. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 23 pages, accepted to MNRAS on 3 August 2023

Tidal dissipation due to the elliptical instability and turbulent viscosity in convection zones in rotating giant planets and stars

Authors: Nils B. de Vries, Adrian J. Barker, Rainer Hollerbach

Abstract: Tidal dissipation in star-planet systems can occur through various mechanisms, among which is the elliptical instability. This acts on elliptically deformed equilibrium tidal flows in rotating fluid planets and stars, and excites inertial waves in convective regions if the dimensionless tidal amplitude ($ε$) is sufficiently large. We study its interaction with turbulent convection, and attempt to… ▽ More Tidal dissipation in star-planet systems can occur through various mechanisms, among which is the elliptical instability. This acts on elliptically deformed equilibrium tidal flows in rotating fluid planets and stars, and excites inertial waves in convective regions if the dimensionless tidal amplitude ($ε$) is sufficiently large. We study its interaction with turbulent convection, and attempt to constrain the contributions of both elliptical instability and convection to tidal dissipation. For this, we perform an extensive suite of Cartesian hydrodynamical simulations of rotating Rayleigh-Bénard convection in a small patch of a planet. We find that tidal dissipation resulting from the elliptical instability, when it operates, is consistent with $ε^3$, as in prior simulations without convection. Convective motions also act as an effective viscosity on large-scale tidal flows, resulting in continuous tidal dissipation (scaling as $ε^2$). We derive scaling laws for the effective viscosity using (rotating) mixing-length theory, and find that they predict the turbulent quantities found in our simulations very well. In addition, we examine the reduction of the effective viscosity for fast tides, which we observe to scale with tidal frequency ($ω$) as $ω^{-2}$. We evaluate our scaling laws using interior models of Hot Jupiters computed with MESA. We conclude that rotation reduces convective length scales, velocities and effective viscosities (though not in the fast tides regime). We estimate that elliptical instability is efficient for the shortest-period Hot Jupiters, and that effective viscosity of turbulent convection is negligible in giant planets compared with inertial waves. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 23 pages, 15 figures, 2 tables; accepted for publication in MNRAS

Tidal dissipation in stratified and semi-convective regions of giant planets

Authors: Christina M. Pontin, Adrian J. Barker, Rainer Hollerbach

Abstract: We study how stably stratified or semi-convective layers alter the tidal dissipation rates associated with the generation of internal waves in planetary interiors. We consider if these layers could contribute to the high rates of tidal dissipation observed for Jupiter and Saturn in our solar system. We use an idealised global spherical Boussinesq model to study the influence of stable stratificati… ▽ More We study how stably stratified or semi-convective layers alter the tidal dissipation rates associated with the generation of internal waves in planetary interiors. We consider if these layers could contribute to the high rates of tidal dissipation observed for Jupiter and Saturn in our solar system. We use an idealised global spherical Boussinesq model to study the influence of stable stratification and semi-convective layers on tidal dissipation rates. We carry out analytical and numerical calculations considering realistic tidal forcing and measure how the viscous and thermal dissipation rates depend on the parameters relating to the internal stratification profile. We find that the strongly frequency-dependent tidal dissipation rate is highly dependent on the parameters relating to the stable stratification, with strong resonant peaks that align with the internal modes of the system. The locations and sizes of these resonances depend on the form and parameters of the stratification, which we explore both analytically and numerically. Our results suggest that stable stratification can significantly enhance the tidal dissipation in particular frequency ranges. Analytical calculations in the low frequency regime give us scaling laws for the key parameters, including the tidal quality factor $Q'$ due to internal gravity waves. Stably stratified layers can significantly contribute to tidal dissipation in solar and extrasolar giant planets, and we estimate substantial tidal evolution for hot Neptunes. Further investigation is needed to robustly quantify the significance of the contribution in realistic interior models, and to consider the contribution of inertial waves. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 29 pages, 17 figures, accepted for publication in ApJ (13th April 2023)

The interactions of the elliptical instability and convection

Authors: N. B. de Vries, A. J. Barker, R. Hollerbach

Abstract: The elliptical instability is an instability of elliptical streamlines, which can be excited by large-scale tidal flows in rotating fluid bodies, and excites inertial waves if the dimensionless tidal amplitude ($ε$) is sufficiently large. It operates in convection zones but its interactions with turbulent convection have not been studied in this context. We perform an extensive suite of Cartesian… ▽ More The elliptical instability is an instability of elliptical streamlines, which can be excited by large-scale tidal flows in rotating fluid bodies, and excites inertial waves if the dimensionless tidal amplitude ($ε$) is sufficiently large. It operates in convection zones but its interactions with turbulent convection have not been studied in this context. We perform an extensive suite of Cartesian hydrodynamical simulations in wide boxes to explore the interactions of the elliptical instability and Rayleigh-Bénard convection. We find that geostrophic vortices generated by the elliptical instability dominate the flow, with energies far exceeding those of the inertial waves. Furthermore, we find that the elliptical instability can operate with convection, but it is suppressed for sufficiently strong convection, primarily by convectively-driven large-scale vortices. We examine the flow in Fourier space, allowing us to determine the energetically dominant frequencies and wavenumbers. We find that power primarily concentrates in geostrophic vortices, in wavenumbers that are convectively unstable, and along the inertial wave dispersion relation, even in non-elliptically deformed convective flows. Examining linear growth rates on a convective background, we find that convective large-scale vortices suppress the elliptical instability in the same way as the geostrophic vortices created by the elliptical instability itself. Finally, convective motions act as an effective viscosity on large-scale tidal flows, providing a sustained energy transfer (scaling as $ε^2$). Furthermore, we find that the energy transfer resulting from bursts of elliptical instability, when it operates, is consistent with the $ε^3$ scaling found in prior work. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 24 pages, 22 figures, 1 table, accepted for publication in Physics of Fluids

Three-dimensional numerical simulations of ambipolar diffusion in NS cores in the one-fluid approximation: instability of poloidal magnetic field

Authors: Andrei P. Igoshev, Rainer Hollerbach

Abstract: We numerically model evolution of magnetic fields inside a neutron star under the influence of ambipolar diffusion in the weak-coupling mode in the one-fluid MHD approximation. Our simulations are three-dimensional and performed in spherical coordinates. Our model covers the neutron star core and includes crust where the magnetic field decay is due to Ohmic decay. We discover an instability of pol… ▽ More We numerically model evolution of magnetic fields inside a neutron star under the influence of ambipolar diffusion in the weak-coupling mode in the one-fluid MHD approximation. Our simulations are three-dimensional and performed in spherical coordinates. Our model covers the neutron star core and includes crust where the magnetic field decay is due to Ohmic decay. We discover an instability of poloidal magnetic field under the influence of ambipolar diffusion. This instability develops in the neutron star core and grows on a timescale of 0.2 dimensionless times, reaching saturation by 2 dimensionless times. The instability leads to formation of azimuthal magnetic field with azimuthal wavenumber $m=14$ (at the moment of saturation) which keeps merging and reaches $m=4$ by 16 dimensionless times. Over the course of our simulations (16 dimensionless times) the surface dipolar magnetic field decays, reaching 20 percent of its original value and keeps decaying. The decay timescale for the total magnetic energy is six dimensionless times. The ambipolar diffusion induces electric currents in the crust where these currents dissipate efficiently. Strong electric currents in the crust lead to heating, which could correspond to luminosities of $\approx 10^{29}$ erg s$^{-1}$ during hundreds of Myrs for an initial magnetic field of $10^{14}$ G. Ambipolar diffusion leads to formation of small-scale magnetic fields at the neutron star surface. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS; 26 pages

Evolution of neutron star magnetic fields

Authors: Andrei P. Igoshev, Sergei B. Popov, Rainer Hollerbach

Abstract: Neutron stars are natural physical laboratories allowing us to study a plethora of phenomena in extreme conditions. In particular, these compact objects can have very strong magnetic fields with non-trivial origin and evolution. In many respects its magnetic field determines the appearance of a neutron star. Thus, understanding the field properties is important for interpretation of observational… ▽ More Neutron stars are natural physical laboratories allowing us to study a plethora of phenomena in extreme conditions. In particular, these compact objects can have very strong magnetic fields with non-trivial origin and evolution. In many respects its magnetic field determines the appearance of a neutron star. Thus, understanding the field properties is important for interpretation of observational data. Complementing this, observations of diverse kinds of neutron stars enable us to probe parameters of electro-dynamical processes at scales unavailable in terrestrial laboratories. In this review we first briefly describe theoretical models of formation and evolution of magnetic field of neutron stars, paying special attention to field decay processes. Then we present important observational results related to field properties of different types of compact objects: magnetars, cooling neutron stars, radio pulsars, sources in binary systems. After that, we discuss which observations can shed light on obscure characteristics of neutron star magnetic fields and their behaviour. We end the review with a subjective list of open problems. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: Review submitted to Universe; 39 pages, 8 figures

3D Magneto-thermal Simulations of Tangled Crustal Magnetic Field in Central Compact Objects

Authors: Andrei P. Igoshev, Konstantinos N. Gourgouliatos, Rainer Hollerbach, Toby S. Wood

Abstract: Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this ca… ▽ More Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magneto-thermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns which consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of $2.5-10\times 10^{47}$ erg have temperatures of hot regions that reach $\approx 0.2$ keV, to be compared with the bulk temperature of $\approx 0.1$ keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of central compact objects. The hot spots produce periodic modulations in light curve with typical amplitudes of $\leq 9-11$ %. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some central compact objects. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 12 pages, 25 figures, accepted for publication ApJ on January 19, 2021

Sleeping beasts: strong toroidal magnetic field in quiescent magnetars explains their large pulsed fraction

Authors: Andrei P. Igoshev, Rainer Hollerbach, Toby Wood, Konstantinos N. Gourgouliatos

Abstract: Magnetars are neutron stars (NSs) with extreme magnetic fields of strength $5 \times 10^{13}$ - $10^{15}$ G. They exhibit transient, highly energetic events, such as short X-ray flashes, bursts and giant flares, all of which are powered by their enormous magnetic energy. Quiescent magnetars have X-ray luminosities between $10^{29}$ and $10^{35}$ erg/s, and are further classified as either persiste… ▽ More Magnetars are neutron stars (NSs) with extreme magnetic fields of strength $5 \times 10^{13}$ - $10^{15}$ G. They exhibit transient, highly energetic events, such as short X-ray flashes, bursts and giant flares, all of which are powered by their enormous magnetic energy. Quiescent magnetars have X-ray luminosities between $10^{29}$ and $10^{35}$ erg/s, and are further classified as either persistent or transient magnetars. Their X-ray emission is modulated with the rotational period of the NS, with a typical relative amplitude (so-called pulsed fraction) between 10-58 per cent, implying that the surface temperature is significantly non-uniform despite the high thermal conductivity of the star's crust. Here, we present the first 3D magneto-thermal MHD simulations of magnetars with strong toroidal magnetic fields. We show that these models, combined with ray propagation in curved space-time, accurately describe the light-curves of most transient magnetars in quiescence and allow us to further constrain their rotational orientation. We find that the presence of a strong toroidal magnetic field explains the observed asymmetry in the surface temperature, and is the main cause of the strong modulation of thermal X-ray emission in quiescence. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: Original version of manuscript submitted to the Nature Astronomy. 11 pages, 6 figures

Powering Central Compact Objects with a Tangled Crustal Magnetic Field

Authors: Konstantinos N. Gourgouliatos, Rainer Hollerbach, Andrei P. Igoshev

Abstract: Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by t… ▽ More Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of $10^{14}$ G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately $10^{10}$ G, as observed in CCOs. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 9 pages, Accepted by MNRAS

Wave propagation in semi-convective regions of giant planets

Authors: Christina M. Pontin, Adrian J. Barker, Rainer Hollerbach, Quentin André, Stéphane Mathis

Abstract: Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilise ordinary convection and produce a stably-stratified layer near the core of these planets. This region could consist of semi-convective layers with a staircase-like density profile, which have multiple convective zones separated by thi… ▽ More Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilise ordinary convection and produce a stably-stratified layer near the core of these planets. This region could consist of semi-convective layers with a staircase-like density profile, which have multiple convective zones separated by thin stably-stratified interfaces, as a result of double-diffusive convection. These layers could have important effects on wave propagation and tidal dissipation that have not been fully explored. We analyse the effects of these layers on the propagation and transmission of internal waves within giant planets, extending prior work in a local Cartesian model. We adopt a simplified global Boussinesq planetary model in which we explore the internal waves in a non-rotating spherical body. We begin by studying the free modes of a region containing semi-convective layers. We then analyse the transmission of internal waves through such a region. The free modes depend strongly on the staircase properties, and consist of modes with both internal and interfacial gravity wave-like behaviour. We determine the frequency shifts of these waves as a function of the number of steps to explore their potential to probe planetary internal structures. We also find that wave transmission is strongly affected by the presence of a staircase. Very large-wavelength waves are transmitted efficiently, but small-scale waves are only transmitted if they are resonant with one of the free modes. The effective size of the core is therefore larger for non-resonant modes. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: 19 pages, 14 figures, accepted for publication in MNRAS on 2nd March 2020

Two types of axisymmetric helical magnetorotational instability in rotating flows with positive shear

Authors: G. Mamatsashvili, F. Stefani, R. Hollerbach, G. Rüdiger

Abstract: We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu… ▽ More We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${\rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${\rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period. △ Less

Submitted 21 September, 2019; v1 submitted 31 October, 2018; originally announced October 2018.

Comments: 15 pages, 13 figures, accepted for publication in Phys. Rev. Fluids

Journal ref: Phys. Rev. Fluids 4, 103905 (2019)

Magnetic fields driven by tidal mixing in radiative stars

Authors: Jérémie Vidal, David Cébron, Nathanaël Schaeffer, Rainer Hollerbach

Abstract: Stellar magnetism plays an important role in stellar evolution theory. Approximatively 10% of observed main sequence (MS) and pre-main-sequence (PMS) radiative stars exhibit surface magnetic fields above the detection limit, raising the question of their origin. These stars host outer radiative envelopes, which are stably stratified. Therefore, they are assumed to be motionless in standard models… ▽ More Stellar magnetism plays an important role in stellar evolution theory. Approximatively 10% of observed main sequence (MS) and pre-main-sequence (PMS) radiative stars exhibit surface magnetic fields above the detection limit, raising the question of their origin. These stars host outer radiative envelopes, which are stably stratified. Therefore, they are assumed to be motionless in standard models of stellar structure and evolution. We focus on rapidly rotating, radiative stars which may be prone to the tidal instability, due to an orbital companion. Using direct numerical simulations in a sphere, we study the interplay between a stable stratification and the tidal instability, and assess its dynamo capability. We show that the tidal instability is triggered regardless of the strength of the stratification (Brunt-V{ä}is{ä}l{ä} frequency). Furthermore, the tidal instability can lead to both mixing and self-induced magnetic fields in stably stratified layers (provided that the Brunt-V{ä}is{ä}l{ä} frequency does not exceed the stellar spin rate in the simulations too much). The application to stars suggests that the resulting magnetic fields could be observable at the stellar surfaces. Indeed, we expect magnetic field strengths up to several Gauss.Consequently, tidally driven dynamos should be considered as a (complementary) dynamo mechanism, possibly operating in radiative MS and PMS stars hosting orbital companions. In particular, tidally driven dynamos may explain the observed magnetism of tidally deformed and rapidly rotating Vega-like stars. △ Less

Submitted 12 January, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Journal ref: Monthly Notices of the Royal Astronomical Society, Oxford University Press (OUP): Policy P - Oxford Open Option A, A Para{î}tre, 000, pp.1 - 18

Magnetic Axis Drift and Magnetic Spot Formation in Neutron Stars with Toroidal Fields

Authors: Konstantinos N. Gourgouliatos, Rainer Hollerbach

Abstract: We explore magnetic field configurations that lead to the formation of magnetic spots on the surface of neutron stars, and to the displacement of the magnetic dipole axis. We find that a toroidally dominated magnetic field is essential for the generation of a single spot with a strong magnetic field. Once a spot forms, it survives for several million years, even after the total magnetic field has… ▽ More We explore magnetic field configurations that lead to the formation of magnetic spots on the surface of neutron stars, and to the displacement of the magnetic dipole axis. We find that a toroidally dominated magnetic field is essential for the generation of a single spot with a strong magnetic field. Once a spot forms, it survives for several million years, even after the total magnetic field has decayed significantly. We find that the dipole axis is not stationary with respect to the neutron star's surface and does not in general coincide with the location of the magnetic spot. This is due to non-axisymmetric instabilities of the toroidal field that displace the poloidal dipole axis at rates that may reach $0.4^{\circ}$ per century. A misaligned poloidal dipole axis with the toroidal field leads to more significant displacement of the dipole axis than the fully aligned case. Finally we discuss the evolution of neutron stars with such magnetic fields on the $P-\dot{P}$ diagram and the observational implications. We find that neutron stars spend a very short time before they cross the Death-Line of the $P-\dot{P}$ diagram, compared to their characteristic ages. Moreover, the maximum intensity of their surface magnetic field is substantially higher than the dipole component of the field. We argue that SGR 0418+5729 could be an example of this type of behaviour, having a weak dipole field, yet hosting a magnetic spot responsible for its magnetar behaviour. The evolution on the pulse profile and braking index of the Crab pulsar, which are attributed to an increase of its obliquity, are compatible with the anticipated drift of the magnetic axis. △ Less

Submitted 13 December, 2017; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 12 pages, 8 figures, Accepted by the Astrophysical Journal

Journal ref: 2018 ApJ 852 21

Transport properties of the Azimuthal Magnetorotational Instability

Authors: Anna Guseva, Ashley P. Willis, Rainer Hollerbach, Marc Avila

Abstract: The magnetorotational instability (MRI) is thought to be a powerful source of turbulence in Keplerian accretion disks. Motivated by recent laboratory experiments, we study the MRI driven by an azimuthal magnetic field in an electrically conducting fluid sheared between two concentric rotating cylinders. By adjusting the rotation rates of the cylinders, we approximate angular velocity profiles… ▽ More The magnetorotational instability (MRI) is thought to be a powerful source of turbulence in Keplerian accretion disks. Motivated by recent laboratory experiments, we study the MRI driven by an azimuthal magnetic field in an electrically conducting fluid sheared between two concentric rotating cylinders. By adjusting the rotation rates of the cylinders, we approximate angular velocity profiles $ω\propto r^{q}$. We perform direct numerical simulations of a steep profile close to the Rayleigh line $q \gtrsim -2 $ and a quasi-Keplerian profile $q \approx -3/2$ and cover wide ranges of Reynolds ($Re\le 4\cdot10^4$) and magnetic Prandtl numbers ($0\le Pm \le 1$). In the quasi-Keplerian case, the onset of instability depends on the magnetic Reynolds number, with $Rm_c \approx 50$, and angular momentum transport scales as $\sqrt{Pm} Re^2$ in the turbulent regime. The ratio of Maxwell to Reynolds stresses is set by $Rm$. At the onset of instability both stresses have similar magnitude, whereas the Reynolds stress vanishes or becomes even negative as $Rm$ increases. For the profile close to the Rayleigh line, the instability shares these properties as long as $Pm\gtrsim0.1$, but exhibits a markedly different character if $Pm\rightarrow 0$, where the onset of instability is governed by the Reynolds number, with $Re_c \approx 1250$, transport is via Reynolds stresses and scales as $Re^2$. At intermediate $Pm=0.01$ we observe a continuous transition from one regime to the other, with a crossover at $Rm=\mathcal{O}(100)$. Our results give a comprehensive picture of angular momentum transport of the MRI with an imposed azimuthal field. △ Less

Submitted 5 October, 2017; v1 submitted 10 May, 2017; originally announced May 2017.

Journal ref: Astrophys. J., 849:92, 2017

Stability and instability of hydromagnetic Taylor-Couette flows

Authors: Günther Rüdiger, Marcus Gellert, Rainer Hollerbach, Manfred Schultz, Frank Stefani

Abstract: Decades ago S. Lundquist, S. Chandrasekhar, P.H. Roberts and R. J.~Tayler first posed questions about the stability of Taylor-Couette flows of conducting material under the influence of large-scale magnetic fields. These and many new questions can now be answered numerically where the nonlinear simulations even provide the instability-induced values of several transport coefficients. The cylindric… ▽ More Decades ago S. Lundquist, S. Chandrasekhar, P.H. Roberts and R. J.~Tayler first posed questions about the stability of Taylor-Couette flows of conducting material under the influence of large-scale magnetic fields. These and many new questions can now be answered numerically where the nonlinear simulations even provide the instability-induced values of several transport coefficients. The cylindrical containers are axially unbounded and penetrated by magnetic background fields with axial and/or azimuthal components. The influence of the magnetic Prandtl number $Pm$ on the onset of the instabilities is shown to be substantial. The potential flow subject to axial fields becomes unstable against axisymmetric perturbations for a certain supercritical value of the averaged Reynolds number $\overline{Rm}=\sqrt{Re\cdot Rm}$ (with $Re$ the Reynolds number of rotation, $Rm$ its magnetic Reynolds number). Rotation profiles as flat as the quasi-Keplerian rotation law scale similarly but only for $Pm\gg 1$ while for $Pm\ll 1$ the instability instead sets in for supercritical $Rm$ at an optimal value of the magnetic field. Among the considered instabilities of azimuthal fields, those of the Chandrasekhar-type, where the background field and the background flow have identical radial profiles, are particularly interesting. They are unstable against nonaxisymmetric perturbations if at least one of the diffusivities is non-zero. For $Pm\ll 1$ the onset of the instability scales with $Re$ while it scales with $\overline{Rm}$ for $Pm\gg 1$. - Even superrotation can be destabilized by azimuthal and current-free magnetic fields; this recently discovered nonaxisymmetric instability is of a double-diffusive character, thus excluding $Pm= 1$. It scales with $Re$ for $Pm\to 0$ and with $Rm$ for $Pm\to \infty$. △ Less

Submitted 5 April, 2018; v1 submitted 29 March, 2017; originally announced March 2017.

Comments: 109 pages, lots of figures

Resistive Tearing Instability in Electron-MHD: Application to Neutron Star Crusts

Authors: Konstantinos N. Gourgouliatos, Rainer Hollerbach

Abstract: We study a resistive tearing instability developing in a system evolving through the combined effect of Hall drift in the Electron-MHD limit and Ohmic dissipation. We explore first the exponential growth of the instability in the linear case and we find the fastest growing mode, the corresponding eigenvalues and dispersion relation. The instability growth rate scales as… ▽ More We study a resistive tearing instability developing in a system evolving through the combined effect of Hall drift in the Electron-MHD limit and Ohmic dissipation. We explore first the exponential growth of the instability in the linear case and we find the fastest growing mode, the corresponding eigenvalues and dispersion relation. The instability growth rate scales as $γ\propto B^{2/3} σ^{-1/3}$ where $B$ is the magnetic field and $σ$ the electrical conductivity. We confirm the development of the tearing resistive instability in the fully non-linear case, in a plane parallel configuration where the magnetic field polarity reverses, through simulations of systems initiating in Hall equilibrium with some superimposed perturbation. Following a transient phase, during which there is some minor rearrangement of the magnetic field, the perturbation grows exponentially. Once the instability is fully developed the magnetic field forms the characteristic islands and X-type reconnection points, where Ohmic decay is enhanced. We discuss the implications of this instability for the local magnetic field evolution in neutron stars' crusts, proposing that it can contribute to heating near the surface of the star, as suggested by models of magnetar post-burst cooling. In particular, we find that a current sheet a few meters thick, covering as little as $1\%$ of the total surface can provide $10^{42}~$erg in thermal energy within a few days. We briefly discuss applications of this instability in other systems where the Hall effect operates such as protoplanetary discs and space plasmas. △ Less

Submitted 8 September, 2016; v1 submitted 26 July, 2016; originally announced July 2016.

Comments: Accepted by MNRAS, 10 pages, 10 Figures

Dynamical model for spindown of solar-type stars

Authors: Aditi Sood, Eun-jin Kim, Rainer Hollerbach

Abstract: Since their formation, stars slow down their rotation rates by the removal of angular momentum from their surfaces, e.g. via stellar winds. Despite the complexity of the processes involved, a traditional model, where the removal of angular momentum loss by magnetic fields is prescribed, has provided a useful framework to understand observational relations between stellar rotation and age and magne… ▽ More Since their formation, stars slow down their rotation rates by the removal of angular momentum from their surfaces, e.g. via stellar winds. Despite the complexity of the processes involved, a traditional model, where the removal of angular momentum loss by magnetic fields is prescribed, has provided a useful framework to understand observational relations between stellar rotation and age and magnetic field strength. Here, a spindown model is proposed where loss of angular momentum by magnetic fields is evolved dynamically, instead of being kinematically prescribed. To this end, we evolve the stellar rotation and magnetic field simultaneously over stellar evolution time by extending our previous work on a dynamo model which incorporates the nonlinear feedback mechanisms on rotation and magnetic fields. Our extended model reproduces key observations and explains the presence of the two branches of (fast and slow rotating) stars which have different relations between rotation rate $Ω$ vs. time (age), magnetic field strength $|B|$ vs. rotation rate, and frequency of magnetic field $ω_{cyc}$ vs. rotation rate. For fast rotating stars we find: (i) an exponential spindown $Ω\propto e^{-1.35t}$, with $t$ measured in Gyrs, (ii) $|B|$ saturates for higher rotation rate, (iii) $ω_{cyc} \propto Ω^{0.85}$. For slow rotating stars we obtain: (i) a power law spindown $Ω\propto t^{-0.52}$, (ii) $|B|$ scales almost linearly with rotation rate, (iii) $ω_{cyc} \propto Ω^{1.16}$. The results obtained are in good agreement with observations. The Vaughan-Preston gap is consistently explained in our model by the shortest spindown timescale in this transition from fast to slow rotators. Our results highlight the importance of self-regulation of magnetic fields and rotation by direct and indirect interactions involving nonlinear feedback in stellar evolution. △ Less

Submitted 18 September, 2016; v1 submitted 23 May, 2016; originally announced May 2016.

Comments: This paper is accepted for publication in ApJ

Suppression of a kinematic dynamo by large shear

Authors: Aditi Sood, Rainer Hollerbach, Eun-jin Kim

Abstract: We numerically solve the magnetic induction equation in a spherical shell geometry, with a kinematically prescribed axisymmetric flow that consists of a superposition of a small-scale helical flow and a large-scale shear flow. The small-scale flow is chosen to be a local analog of the classical Roberts cells, consisting of strongly helical vortex rolls. The large-scale flow is a shearing motion in… ▽ More We numerically solve the magnetic induction equation in a spherical shell geometry, with a kinematically prescribed axisymmetric flow that consists of a superposition of a small-scale helical flow and a large-scale shear flow. The small-scale flow is chosen to be a local analog of the classical Roberts cells, consisting of strongly helical vortex rolls. The large-scale flow is a shearing motion in either the radial or the latitudinal directions. In the absence of large-scale shear, the small-scale flow is an efficient dynamo, in agreement with previous results. Adding increasingly large shear flows strongly suppresses the dynamo efficiency, indicating that shear is not always a favourable ingredient in dynamo action. △ Less

Submitted 18 September, 2016; v1 submitted 2 May, 2016; originally announced May 2016.

Comments: This paper is accepted for publication in Journal of Physics A and in press

Magnetic field evolution in magnetar crusts through three dimensional simulations

Authors: Konstantinos N. Gourgouliatos, Toby Wood, Rainer Hollerbach

Abstract: Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we presen… ▽ More Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we present detailed modelling of the evolution of the magnetic field in the crust of a neutron star through 3-D simulations. We find that, in the plausible scenario of equipartition of energy between global-scale poloidal and toroidal magnetic components, magnetic instabilities transfer energy to non-axisymmetric, kilometre-sized magnetic features, in which the local field strength can greatly exceed that of the global-scale field. These intense small-scale magnetic features can induce high energy bursts through local crust yielding, and the localised enhancement of Ohmic heating can power the star's persistent emission. Thus, the observed diversity in magnetar behaviour can be explained with mixed poloidal-toroidal fields of comparable energies. △ Less

Submitted 5 April, 2016; originally announced April 2016.

Comments: Published by PNAS

Nonaxisymmetric MHD instabilities of Chandrasekhar states in Taylor-Couette geometry

Authors: M. Gellert, G. Rüdiger, M. Schultz, A. Guseva, R. Hollerbach

Abstract: We consider axially periodic Taylor-Couette geometry with insulating boundary conditions. The imposed basic states are so-called Chandrasekhar states, where the azimuthal flow $U_φ$ and magnetic field $B_φ$ have the same radial profiles. Mainly three particular profiles are considered: the Rayleigh limit, quasi-Keplerian, and solid-body rotation. In each case we begin by computing linear instabili… ▽ More We consider axially periodic Taylor-Couette geometry with insulating boundary conditions. The imposed basic states are so-called Chandrasekhar states, where the azimuthal flow $U_φ$ and magnetic field $B_φ$ have the same radial profiles. Mainly three particular profiles are considered: the Rayleigh limit, quasi-Keplerian, and solid-body rotation. In each case we begin by computing linear instability curves and their dependence on the magnetic Prandtl number Pm. For the azimuthal wavenumber m=1 modes, the instability curves always scale with the Reynolds number and the Hartmann number. For sufficiently small Pm these modes therefore only become unstable for magnetic Mach numbers less than unity, and are thus not relevant for most astrophysical applications. However, modes with m>10 can behave very differently. For sufficiently flat profiles, they scale with the magnetic Reynolds number and the Lundquist number, thereby allowing instability also for the large magnetic Mach numbers of astrophysical objects. We further compute fully nonlinear, three-dimensional equilibration of these instabilities, and investigate how the energy is distributed among the azimuthal (m) and axial (k) wavenumbers. In comparison spectra become steeper for large m, reflecting the smoothing action of shear. On the other hand kinetic and magnetic energy spectra exhibit similar behavior: if several azimuthal modes are already linearly unstable they are relatively flat, but for the rigidly rotating case where m=1 is the only unstable mode they are so steep that neither Kolmogorov nor Iroshnikov-Kraichnan spectra fit the results. The total magnetic energy exceeds the kinetic energy only for large magnetic Reynolds numbers Rm>100. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: 12 pages, 14 figures, submitted to ApJ

Generalised quasi-linear approximation of the HMRI

Authors: Adam Child, Rainer Hollerbach, Brad Marston, Steven Tobias

Abstract: Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-… ▽ More Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in describing the statistics of the HMRI even when relatively few large-scale modes are kept in the formalism. We conclude that DSS based on GQL (GCE2) will be significantly more accurate than that based on QL (CE2). △ Less

Submitted 19 February, 2016; originally announced February 2016.

Comments: 17 pages, 8 figures

Magnetar Activity via the Density-Shear Instability in Hall-MHD

Authors: Konstantinos N. Gourgouliatos, Todor Kondic, Maxim Lyutikov, Rainer Hollerbach

Abstract: We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing… ▽ More We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing electron number density and magnetic field, a plane-parallel analogue of an azimuthal or meridional magnetic field in the crust of a magnetar. The growth rate of the instability depends on the Hall properties of the field (magnetic field intensity, electron number density and the corresponding scale-heights), while being insensitive to weak resistivity. Since the Hall effect is the driving process for the evolution of the crustal magnetic field of magnetars, we argue that this instability is critical for systems containing strong meridional or azimuthal fields. We find that this process mediates the formation of localised structures with much stronger magnetic field than the average, which can lead to magnetar activity and accelerate the dissipation of the field and consequently the production of Ohmic heating. Assuming a $5\times10^{14}$G magnetic field at the base of crust, we anticipate that magnetic field as strong as $10^{15}$G will easily develop in regions of typical size of a few $10^{2}$ meters, containing magnetic energy of $10^{43}$erg, sufficient to power magnetar bursts. These active regions are more likely to appear in the magnetic equator where the tangential magnetic field is stronger. △ Less

Submitted 27 July, 2015; originally announced July 2015.

Comments: Accepted by MNRAS Letters, 3 Figures, 1 Table

Three-dimensional simulations of the magnetic stress in a neutron star crust

Authors: Toby S. Wood, Rainer Hollerbach

Abstract: We present the first fully self-consistent three-dimensional model of a neutron star's magnetic field, generated by electric currents in the star's crust via the Hall effect. We find that the global-scale field converges to a Hall-attractor state, as seen in recent axisymmetric models, but that small-scale features in the magnetic field survive even on much longer timescales. These small-scale fea… ▽ More We present the first fully self-consistent three-dimensional model of a neutron star's magnetic field, generated by electric currents in the star's crust via the Hall effect. We find that the global-scale field converges to a Hall-attractor state, as seen in recent axisymmetric models, but that small-scale features in the magnetic field survive even on much longer timescales. These small-scale features propagate toward the dipole equator, where the crustal electric currents organize themselves into a strong equatorial jet. By calculating the distribution of magnetic stresses in the crust, we predict that neutron stars with fields stronger than $10^{14}$G can still be subject to starquakes more than $10^5$yr after their formation. △ Less

Submitted 12 June, 2015; v1 submitted 21 January, 2015; originally announced January 2015.

Comments: Accepted version, prior to copy editing

Journal ref: Phys. Rev. Lett. 114, 191101 (2015)

Tidally driven dynamos in a rotating sphere

Authors: David Cébron, Rainer Hollerbach

Abstract: Large-scale planetary or stellar magnetic fields generated by a dynamo effect are mostly attributed to flows forced by buoyancy forces in electrically conducting fluid layers. However, these large-scale fields may also be controlled by tides, as previously suggested for the star $τ$-boo, Mars or the Early Moon. By simulating a small local patch of a rotating fluid, \cite{Barker2014} have recently… ▽ More Large-scale planetary or stellar magnetic fields generated by a dynamo effect are mostly attributed to flows forced by buoyancy forces in electrically conducting fluid layers. However, these large-scale fields may also be controlled by tides, as previously suggested for the star $τ$-boo, Mars or the Early Moon. By simulating a small local patch of a rotating fluid, \cite{Barker2014} have recently shown that tides can drive small-scale dynamos by exciting a hydrodynamic instability, the so-called elliptical (or tidal) instability. By performing global magnetohydrodynamic simulations of a rotating spherical fluid body, we investigate if this instability can also drive the observed large-scale magnetic fields. We are thus interested by the dynamo threshold and the generated magnetic field in order to test if such a mechanism is relevant for planets and stars. Rather than solving the problem in a geometry deformed by tides, we consider a spherical fluid body and add a body force to mimic the tidal deformation in the bulk of the fluid. This allows us to use an efficient spectral code to solve the magnetohydrodynamic problem. We first compare the hydrodynamic results with theoretical asymptotic results, and numerical results obtained in a truely deformed ellipsoid, which confirms the presence of the elliptical instability. We then perform magnetohydrodynamic simulations, and investigate the dynamo capability of the flow. Kinematic and self-consistent dynamos are finally simulated, showing that the elliptical instability is capable of generating dipole dominated large-scale magnetic field in global simulations of a fluid rotating sphere. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: Astrophysical Journal Letters In press, (accepted) (2014) (accepted)

The effects of rotation on a double-diffusive layer in a rotating spherical shell

Authors: Patrick Blies, Friedrich Kupka, Florian Zaussinger, Rainer Hollerbach

Abstract: So far, numerical studies of double-diffusive layering in turbulent convective flows have neglected the effects of rotation. We undertake a first step into that direction by investigating how Coriolis forces affect a double-diffusive layer inside a rotating spherical shell. For this purpose we have run simulations in a parameter regime where these layers are expected to form and successively incre… ▽ More So far, numerical studies of double-diffusive layering in turbulent convective flows have neglected the effects of rotation. We undertake a first step into that direction by investigating how Coriolis forces affect a double-diffusive layer inside a rotating spherical shell. For this purpose we have run simulations in a parameter regime where these layers are expected to form and successively increased the rate of rotation with the result that fast rotation is found to have a similar stabilising effect on the overall convective flux as an increase of the stability ratio $R_ρ$ has in a non-rotating setup. We have also studied to what extent the regimes of rotational constraints suggested by King, Stellmach, and Buffett (2013) for rotation in the case of Rayleigh-Bénard convection are influenced by double-diffusive convection: their classification could also be applicable to the case of double-diffusive convection in a spherical shell if it is extended to be also a function of the stability ratio $R_ρ$. Furthermore, we examined the ratio of saline and thermal Nusselt numbers and compared our results with models of Spruit (2013), Rosenblum et al. (2011) and Wood, Garaud, and Stellmach (2013). We find our data to be fitted best by Spruit's model. Our result that fast rotation further decreases the convective transport, which is already lowered by double-diffusive convection, could play a major role for e.g. the modeling of the interior of some rapidly rotating giant planets, as gaseous giant planets have recently been proposed to be influenced by double-diffusive convection. △ Less

Submitted 24 April, 2014; originally announced April 2014.

Comments: 31 pages, 15 figures, 3 tables

Density-shear instability in electron MHD

Authors: Toby S Wood, Rainer Hollerbach, Maxim Lyutikov

Abstract: We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic r… ▽ More We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic reconnection regions on scales smaller than the ion skin depth, and in neutron star crusts. We demonstrate that the so-called Hall drift instability, previously argued to be relevant in neutron star crusts, is a resistive tearing instability rather than an instability of the Hall term itself. We argue that the density-shear instability is of greater significance in neutron stars than the tearing instability, because it generally has a faster growth-rate and is less sensitive to geometry and boundary conditions. We prove that, for uniform electron density, EMHD is "at least as stable" as regular, incompressible MHD, in the sense that any field configuration that is stable in MHD is also stable in EMHD. We present a connection between the density-shear instability in EMHD and the magneto-buoyancy instability in anelastic MHD. △ Less

Submitted 11 July, 2014; v1 submitted 8 April, 2014; originally announced April 2014.

Comments: 1 Figure

Journal ref: Physics of Plasmas, 2014, 21, 052110

Experimental evidence for non-axisymmetric magnetorotational instability in a rotating liquid metal exposed to an azimuthal magnetic field

Authors: Martin Seilmayer, Vladimir Galindo, Gunter Gerbeth, Thomas Gundrum, Frank Stefani, Marcus Gellert, Günther Rüdiger, Manfred Schultz, Rainer Hollerbach

Abstract: The azimuthal version of the magnetorotational instability (MRI) is a non-axisymmetric instability of a hydrodynamically stable differentially rotating flow under the influence of a purely or predominantly azimuthal magnetic field. It may be of considerable importance for destabilizing accretion disks, and plays a central role in the concept of the MRI dynamo. We report the results of a liquid met… ▽ More The azimuthal version of the magnetorotational instability (MRI) is a non-axisymmetric instability of a hydrodynamically stable differentially rotating flow under the influence of a purely or predominantly azimuthal magnetic field. It may be of considerable importance for destabilizing accretion disks, and plays a central role in the concept of the MRI dynamo. We report the results of a liquid metal Taylor-Couette experiment that shows the occurrence of an azimuthal MRI in the expected range of Hartmann numbers. △ Less

Submitted 25 June, 2014; v1 submitted 28 November, 2013; originally announced November 2013.

Comments: 5 pages, 5 figures, accepted at Physical Review Letters

Astrophysical and experimental implications from the magnetorotational instability of toroidal fields

Authors: G. Ruediger, M. Gellert, M. Schultz, R. Hollerbach, F. Stefani

Abstract: The interaction of differential rotation and toroidal fields that are current-free in the gap between two corotating axially unbounded cylinders is considered. It is shown that nonaxisymmetric perturbations are unstable if the rotation rate and Alfvén frequency of the field are of the same order, almost independent of the magnetic Prandtl number Pm. For the very steep rotation law Ω\propto R^{-2}… ▽ More The interaction of differential rotation and toroidal fields that are current-free in the gap between two corotating axially unbounded cylinders is considered. It is shown that nonaxisymmetric perturbations are unstable if the rotation rate and Alfvén frequency of the field are of the same order, almost independent of the magnetic Prandtl number Pm. For the very steep rotation law Ω\propto R^{-2} (the Rayleigh limit) and for small Pm the threshold values of rotation and field for this Azimuthal MagnetoRotational Instability (AMRI) scale with the ordinary Reynolds number and the Hartmann number, resp. A laboratory experiment with liquid metals like sodium or gallium in a Taylor-Couette container has been designed on the basis of this finding. For fluids with more flat rotation laws the Reynolds number and the Hartmann number are no longer typical quantities for the instability. For the weakly nonlinear system the numerical values of the kinetic energy and the magnetic energy are derived for magnetic Prandtl numbers \leq 1. We find that the magnetic energy grows monotonically with the magnetic Reynolds number Rm, while the kinetic energy grows with Rm/\sqrt{Pm}. The resulting turbulent Schmidt number, as the ratio of the `eddy' viscosity and the diffusion coefficient of a passive scalar (such as lithium) is of order 20 for Pm=1, but for small Pm it drops to order unity. Hence, in a stellar core with fossil fields and steep rotation law the transport of angular momentum by AMRI is always accompanied by an intense mixing of the plasma, until the rotation becomes rigid. △ Less

Submitted 6 November, 2013; v1 submitted 19 March, 2013; originally announced March 2013.

Comments: 8 pages, 11 figures

Rapidly rotating plane layer convection with zonal flow

Authors: Robert J. Teed, Chris A. Jones, Rainer Hollerbach

Abstract: The onset of convection in a rapidly rotating layer in which a thermal wind is present is studied. Diffusive effects are included. The main motivation is from convection in planetary interiors, where thermal winds are expected due to temperature variations on the core-mantle boundary. The system admits both convective instability and baroclinic instability. We find a smooth transition between the… ▽ More The onset of convection in a rapidly rotating layer in which a thermal wind is present is studied. Diffusive effects are included. The main motivation is from convection in planetary interiors, where thermal winds are expected due to temperature variations on the core-mantle boundary. The system admits both convective instability and baroclinic instability. We find a smooth transition between the two types of modes, and investigate where the transition region between the two types of instability occurs in parameter space. The thermal wind helps to destabilise the convective modes. Baroclinic instability can occur when the applied vertical temperature gradient is stable, and the critical Rayleigh number is then negative. Long wavelength modes are the first to become unstable. Asymptotic analysis is possible for the transition region and also for long wavelength instabilities, and the results agree well with our numerical solutions. We also investigate how the instabilities in this system relate to the classical baroclinic instability in the Eady problem. We conclude by noting that baroclinic instabilities in the Earth's core arising from heterogeneity in the lower mantle could possibly drive a dynamo even if the Earth's core were stably stratified and so not convecting. △ Less

Submitted 1 November, 2011; originally announced November 2011.

Comments: 20 pages, 7 figures

Journal ref: Geophys. Astrophys. Fluid Dynam. 104 (2010) 457-480

The shear-Hall instability in newborn neutron stars

Authors: Todor Kondić, Günther Rüdiger, Rainer Hollerbach

Abstract: Aims. In the first few minutes of a newborn neutron star's life the Hall effect and differential rotation may both be important. We demonstrate that these two ingredients are sufficient for generating a 'shear-Hall instability' and for studying its excitation conditions, growth rates, and characteristic magnetic field patterns. Methods. We numerically solve the induction equation in a spherical… ▽ More Aims. In the first few minutes of a newborn neutron star's life the Hall effect and differential rotation may both be important. We demonstrate that these two ingredients are sufficient for generating a 'shear-Hall instability' and for studying its excitation conditions, growth rates, and characteristic magnetic field patterns. Methods. We numerically solve the induction equation in a spherical shell, with a kinematically prescribed differential rotation profile Ω(s), where s is the cylindrical radius. The Hall term is linearized about an imposed uniform axial field. The linear stability of individual azimuthal modes, both axisymmetric and non-axisymmetric, is then investigated. Results. For the shear-Hall instability to occur, the axial field must be parallel to the rotation axis if Ω(s) decreases outward, whereas if Ω(s) increases outward it must be anti-parallel. The instability draws its energy from the differential rotation, and occurs on the short rotational timescale rather than on the much longer Hall timescale. It operates most efficiently if the Hall time is comparable to the diffusion time. Depending on the precise field strengths B0, either axisymmetric or non-axisymmetric modes may be the most unstable. Conclusions. Even if the differential rotation in newborn neutron stars is quenched within minutes, the shear-Hall instability may nevertheless amplify any seed magnetic fields by many orders of magnitude. △ Less

Submitted 18 October, 2011; originally announced October 2011.

Comments: 4 pages, 3 figures

Helicity and alpha-effect by current-driven instabilities of helical magnetic fields

Authors: M. Gellert, G. Ruediger, R. Hollerbach

Abstract: Helical magnetic background fields with adjustable pitch angle are imposed on a conducting fluid in a differentially rotating cylindrical container. The small-scale kinetic and current helicities are calculated for various field geometries, and shown to have the opposite sign as the helicity of the large-scale field. These helicities and also the corresponding $α$-effect scale with the current hel… ▽ More Helical magnetic background fields with adjustable pitch angle are imposed on a conducting fluid in a differentially rotating cylindrical container. The small-scale kinetic and current helicities are calculated for various field geometries, and shown to have the opposite sign as the helicity of the large-scale field. These helicities and also the corresponding $α$-effect scale with the current helicity of the background field. The $α$-tensor is highly anisotropic as the components $α_{φφ}$ and $α_{zz}$ have opposite signs. The amplitudes of the azimuthal $α$-effect computed with the cylindrical 3D MHD code are so small that the operation of an $αΩ$ dynamo on the basis of the current-driven, kink-type instabilities of toroidal fields is highly questionable. In any case the low value of the $α$-effect would lead to very long growth times of a dynamo in the radiation zone of the Sun and early-type stars of the order of mega-years. △ Less

Submitted 4 January, 2011; v1 submitted 3 January, 2011; originally announced January 2011.

Comments: 6 pages, 7 figures, submitted to MNRAS

Kinematic dynamo in spherical Couette flow

Authors: Xing Wei, Andrew Jackson, Rainer Hollerbach

Abstract: We investigate numerically kinematic dynamos driven by flow of electrically conducting fluid in the shell between two concentric differentially rotating spheres, a configuration normally referred to as spherical Couette flow. We compare between axisymmetric (2D) and fully three dimensional flows, between low and high global rotation rates, between prograde and retrograde differential rotations, be… ▽ More We investigate numerically kinematic dynamos driven by flow of electrically conducting fluid in the shell between two concentric differentially rotating spheres, a configuration normally referred to as spherical Couette flow. We compare between axisymmetric (2D) and fully three dimensional flows, between low and high global rotation rates, between prograde and retrograde differential rotations, between weak and strong nonlinear inertial forces, between insulating and conducting boundaries, and between two aspect ratios. The main results are as follows. Azimuthally drifting Rossby waves arising from the destabilisation of the Stewartson shear layer are crucial to dynamo action. Differential rotation and helical Rossby waves combine to contribute to the spherical Couette dynamo. At a slow global rotation rate, the direction of differential rotation plays an important role in the dynamo because of different patterns of Rossby waves in prograde and retrograde flows. At a rapid global rotation rate, stronger flow supercriticality (namely the difference between the differential rotation rate of the flow and its critical value for the onset of nonaxisymmetric instability) facilitates the onset of dynamo action. A conducting magnetic boundary condition and a larger aspect ratio both favour dynamo action. △ Less

Submitted 6 August, 2012; v1 submitted 16 November, 2010; originally announced November 2010.

Instabilities of Shercliff and Stewartson layers in spherical Couette flow

Authors: Xing Wei, Rainer Hollerbach

Abstract: We explore numerically the flow induced in a spherical shell by differentially rotating the inner and outer spheres. The fluid is also taken to be electrically conducting (in the low magnetic Reynolds number limit), and a magnetic field is imposed parallel to the axis of rotation. If the outer sphere is stationary, the magnetic field induces a Shercliffe layer on the tangent cylinder, the cylinder… ▽ More We explore numerically the flow induced in a spherical shell by differentially rotating the inner and outer spheres. The fluid is also taken to be electrically conducting (in the low magnetic Reynolds number limit), and a magnetic field is imposed parallel to the axis of rotation. If the outer sphere is stationary, the magnetic field induces a Shercliffe layer on the tangent cylinder, the cylinder just touching the inner sphere and parallel to the field. If the magnetic field is absent, but a strong overall rotation is present, Coriolis effects induce a Stewartson layer on the tangent cylinder. The non-axisymmetric instabilities of both types of layer separately have been studied before; here we consider the two cases side by side, as well as the mixed case, and investigate how magnetic and rotational effects interact. We find that if the differential rotation and the overall rotation are in the same direction, the overall rotation may have a destabilizing influence, whereas if the differential rotation and the overall rotation are in the opposite direction, the overall rotation always has a stabilizing influence. △ Less

Submitted 2 August, 2010; originally announced August 2010.

Magnetic spherical Couette flow in linear combinations of axial and dipolar fields

Authors: Xing Wei, Rainer Hollerbach

Abstract: We present axisymmetric numerical calculations of the fluid flow induced in a spherical shell with inner sphere rotating and outer sphere stationary. A magnetic field is also imposed, consisting of particular linear combinations of axial and dipolar fields, chosen to make $B_r=0$ at either the outer sphere, or the inner, or in between. This leads to the formation of Shercliff shear layers at these… ▽ More We present axisymmetric numerical calculations of the fluid flow induced in a spherical shell with inner sphere rotating and outer sphere stationary. A magnetic field is also imposed, consisting of particular linear combinations of axial and dipolar fields, chosen to make $B_r=0$ at either the outer sphere, or the inner, or in between. This leads to the formation of Shercliff shear layers at these particular locations. We then consider the effect of increasingly large inertial effects, and show that an outer Shercliff layer is eventually de-stabilized, an inner Shercliff layer appears to remain stable, and an in-between Shercliff layer is almost completely disrupted even before the onset of time-dependence, which does eventually occur though. △ Less

Submitted 2 August, 2010; originally announced August 2010.

Dissipative Taylor-Couette flows under the influence of helical magnetic fields

Authors: G. Ruediger, M. Gellert, M. Schultz, R. Hollerbach

Abstract: The linear stability of MHD Taylor-Couette flows in axially unbounded cylinders is considered, for magnetic Prandtl number unity. Magnetic fields varying from purely axial to purely azimuthal are imposed, with a general helical field parameterized by β=B_φ/B_z. We map out the transition from the standard MRI for β=0 to the nonaxisymmetric Azimuthal MagnetoRotational Instability (AMRI) for β\to \… ▽ More The linear stability of MHD Taylor-Couette flows in axially unbounded cylinders is considered, for magnetic Prandtl number unity. Magnetic fields varying from purely axial to purely azimuthal are imposed, with a general helical field parameterized by β=B_φ/B_z. We map out the transition from the standard MRI for β=0 to the nonaxisymmetric Azimuthal MagnetoRotational Instability (AMRI) for β\to \infty. For finite β, positive and negative wave numbers m, corresponding to right and left spirals, are no longer identical. The transition from β=0 to β\to\infty includes all the possible forms of MRI with axisymmetric and nonaxisymmetric modes. For the nonaxisymmetric modes, the most unstable mode spirals in the opposite direction to the background field. The standard (β=0) MRI is axisymmetric for weak fields (including the instability with the lowest Reynolds number) but is nonaxisymmetric for stronger fields. If the azimuthal field is due in part to an axial current flowing through the fluid itself (and not just along the central axis), then it is also unstable to the nonaxisymmetric Tayler instability, which is most effective without rotation. For large βthis instability has wavenumber m=1, whereas for β\simeq 1 m=2 is most unstable. The most unstable mode spirals in the same direction as the background field. △ Less

Submitted 12 January, 2010; originally announced January 2010.

Comments: 9 pages, 11 figures

Hall cascades versus instabilities in neutron star magnetic fields

Authors: C. J. Wareing, R. Hollerbach

Abstract: The Hall effect is an important nonlinear mechanism affecting the evolution of magnetic fields in neutron stars. Studies of the governing equation, both theoretical and numerical, have shown that the Hall effect proceeds in a turbulent cascade of energy from large to small scales. We investigate the small-scale Hall instability conjectured to exist from the linear stability analysis of Rheinhard… ▽ More The Hall effect is an important nonlinear mechanism affecting the evolution of magnetic fields in neutron stars. Studies of the governing equation, both theoretical and numerical, have shown that the Hall effect proceeds in a turbulent cascade of energy from large to small scales. We investigate the small-scale Hall instability conjectured to exist from the linear stability analysis of Rheinhardt and Geppert. Identical linear stability analyses are performed to find a suitable background field to model Rheinhardt and Geppert's ideas. The nonlinear evolution of this field is then modelled using a three-dimensional pseudospectral numerical MHD code. Combined with the background field, energy was injected at the ten specific eigenmodes with the greatest positive eigenvalues as inferred by the linear stability analysis. Energy is transferred to different scales in the system, but not into small scales to any extent that could be interpreted as a Hall instability. Any instabilities are overwhelmed by a late-onset turbulent Hall cascade, initially avoided by the choice of background field, but soon generated by nonlinear interactions between the growing eigenmodes. The Hall cascade is shown here, and by several authors elsewhere, to be the dominant mechanism in this system. △ Less

Submitted 14 December, 2009; originally announced December 2009.

Comments: 7 pages, 3 figures, accepted to Astronomy & Astrophysics

Nonlinear Evolution of Axisymmetric Twisted Flux Tubes in the Solar Tachocline

Authors: R. Hollerbach, P. S. Cally

Abstract: We numerically study the evolution of magnetic fields and fluid flows in a thin spherical shell. We take the initial field to be a latitudinally confined, predominantly toroidal flux tube. For purely toroidal, untwisted flux tubes, we recover previously known radial-shredding instabilities, and show further that in the nonlinear regime these instabilities can very effectively destroy the origina… ▽ More We numerically study the evolution of magnetic fields and fluid flows in a thin spherical shell. We take the initial field to be a latitudinally confined, predominantly toroidal flux tube. For purely toroidal, untwisted flux tubes, we recover previously known radial-shredding instabilities, and show further that in the nonlinear regime these instabilities can very effectively destroy the original field. For twisted flux tubes, including also a poloidal component, there are several possibilities, including the suppression of the radial-shredding instability, but also a more directly induced evolution, brought about because twisted flux tubes in general are not equilibrium solutions of the governing equations. △ Less

Submitted 29 December, 2009; v1 submitted 21 September, 2009; originally announced September 2009.

Comments: final version as published in journal

Journal ref: Solar Physics 260, 251-260, 2009

Helical magnetorotational instability in a Taylor-Couette flow with strongly reduced Ekman pumping

Authors: Frank Stefani, Gunter Gerbeth, Thomas Gundrum, Rainer Hollerbach, Jānis Priede, Günther Rüdiger, Jacek Szklarski

Abstract: The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion discs. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interf… ▽ More The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion discs. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interfered with a significant Ekman pumping driven by solid end-caps that confined the instability only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split end-caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper and in better agreement with numerical predictions. By analyzing various parameter dependencies we conclude that the observed HMRI represents a self-sustained global instability rather than a noise-sustained convective one. △ Less

Submitted 7 April, 2009; originally announced April 2009.

Comments: 30 pages, 22 figures, submitted to PRE

Journal ref: Phys.Rev.E80:066303,2009

Liquid metal experiments on the helical magnetorotational instability

Authors: Frank Stefani, Gunter Gerbeth, Thomas Gundrum, Jacek Szklarski, Guenther Ruediger, Rainer Hollerbach

Abstract: The magnetorotational instability (MRI) plays an essential role in the formation of stars and black holes. By destabilizing hydrodynamically stable Keplerian flows, the MRI triggers turbulence and enables outward transport of angular momentum in accretion discs. We present the results of a liquid metal Taylor-Couette experiment under the influence of helical magnetic fields that show typical fea… ▽ More The magnetorotational instability (MRI) plays an essential role in the formation of stars and black holes. By destabilizing hydrodynamically stable Keplerian flows, the MRI triggers turbulence and enables outward transport of angular momentum in accretion discs. We present the results of a liquid metal Taylor-Couette experiment under the influence of helical magnetic fields that show typical features of MRI at Reynolds numbers of the order 1000 and Hartmann numbers of the order 10. Particular focus is laid on an improved experiment in which split end caps are used to minimize the Ekman pumping. △ Less

Submitted 19 December, 2008; originally announced December 2008.

Comments: 7 pages, 4 figures, submitted to Magnetohydrodynamics

Comment on "Helical MRI in magnetized Taylor-Couette flow"

Authors: G. Ruediger, R. Hollerbach

Abstract: Liu et al. [Phys. Rev. E 74, 056302 (2006)] have presented a WKB analysis of the helical magnetorotational instability (HMRI), and claim that it does not exist for Keplerian rotation profiles. We show that if radial boundary conditions are included, the HMRI can exist even for rotation profiles as flat as Keplerian, provided only that at least one of the boundaries is sufficiently conducting. Liu et al. [Phys. Rev. E 74, 056302 (2006)] have presented a WKB analysis of the helical magnetorotational instability (HMRI), and claim that it does not exist for Keplerian rotation profiles. We show that if radial boundary conditions are included, the HMRI can exist even for rotation profiles as flat as Keplerian, provided only that at least one of the boundaries is sufficiently conducting. △ Less

Submitted 10 May, 2007; v1 submitted 7 May, 2007; originally announced May 2007.

Comments: 2 pages, 2 figures, submitted to Phys. Rev. E (comments)

Destabilisation of hydrodynamically stable rotation laws by azimuthal magnetic fields

Authors: G. Ruediger, R. Hollerbach, M. Schultz, D. Elstner

Abstract: We consider the effect of toroidal magnetic fields on hydrodynamically stable Taylor-Couette differential rotation flows. For current-free magnetic fields a nonaxisymmetric m=1 magnetorotational instability arises when the magnetic Reynolds number exceeds O (100). We then consider how this `azimuthal magnetorotational instability' (AMRI) is modified if the magnetic field is not current-free, but… ▽ More We consider the effect of toroidal magnetic fields on hydrodynamically stable Taylor-Couette differential rotation flows. For current-free magnetic fields a nonaxisymmetric m=1 magnetorotational instability arises when the magnetic Reynolds number exceeds O (100). We then consider how this `azimuthal magnetorotational instability' (AMRI) is modified if the magnetic field is not current-free, but also has an associated electric current throughout the fluid. This gives rise to current-driven Tayler instabilities (TI) that exist even without any differential rotation at all. The interaction of the AMRI and the TI is then considered when both electric currents and differential rotation are present simultaneously. The magnetic Prandtl number Pm turns out to be crucial in this case. Large Pm have a destabilizing influence, and lead to a smooth transition between the AMRI and the TI. In contrast, small Pm have a stabilizing influence, with a broad stable zone separating the AMRI and the TI. In this region the differential rotation is acting to stabilize the Tayler instabilities, with possible astrophysical applications (Ap stars). The growth rates of both the AMRI and the TI are largely independent of Pm, with the TI acting on the timescale of a single rotation period, and the AMRI slightly slower, but still on the basic rotational timescale. The azimuthal drift timescale is ~20 rotations, and may thus be a (flip-flop) timescale of stellar activity between the rotation period and the diffusion time. △ Less

Submitted 16 January, 2007; originally announced January 2007.

Comments: 7 pages, 8 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.377:1481-1487,2007

Experiments on the magnetorotational instability in helical magnetic fields

Authors: Frank Stefani, Thomas Gundrum, Gunter Gerbeth, Günther Rüdiger, Jacek Szklarski, Rainer Hollerbach

Abstract: The magnetorotational instability (MRI) plays a key role in the formation of stars and black holes, by enabling outward angular momentum transport in accretion disks. The use of combined axial and azimuthal magnetic fields allows the investigation of this effect in liquid metal flows at moderate Reynolds and Hartmann numbers. A variety of experimental results is presented showing evidence for th… ▽ More The magnetorotational instability (MRI) plays a key role in the formation of stars and black holes, by enabling outward angular momentum transport in accretion disks. The use of combined axial and azimuthal magnetic fields allows the investigation of this effect in liquid metal flows at moderate Reynolds and Hartmann numbers. A variety of experimental results is presented showing evidence for the occurrence of the MRI in a Taylor-Couette flow using the liquid metal alloy GaInSn. △ Less

Submitted 2 January, 2007; originally announced January 2007.

Comments: 15 pages, 13 figures, submitted to New Journal of Physics

Journal ref: NewJ.Phys.9:295,2007

A non-axisymmetric magnetorotational instability of a purely toroidal magnetic field

Authors: R. Hollerbach, G. Rudiger, M. Schultz, D. Elstner

Abstract: We consider the flow of an electrically conducting fluid between differentially rotating cylinders, in the presence of an externally imposed toroidal field B_0 (r_i/r) e_phi. It is known that the classical, axisymmetric magnetorotational instability does not exist for such a purely toroidal imposed field. We show here that a non-axisymmetric magnetorotational instability does exist, having prope… ▽ More We consider the flow of an electrically conducting fluid between differentially rotating cylinders, in the presence of an externally imposed toroidal field B_0 (r_i/r) e_phi. It is known that the classical, axisymmetric magnetorotational instability does not exist for such a purely toroidal imposed field. We show here that a non-axisymmetric magnetorotational instability does exist, having properties very similar to the axisymmetric magnetorotational instability in the presence of an axial field. △ Less

Submitted 16 November, 2006; originally announced November 2006.

Theory of current instability experiments in magnetic Taylor-Couette flows

Authors: G. Ruediger, M. Schultz, D. Shalybkov, R. Hollerbach

Abstract: We consider the linear stability of dissipative MHD Taylor-Couette flow with imposed toroidal magnetic fields. The inner and outer cylinders can be either insulating or conducting; the inner one rotates, the outer one is stationary. The magnetic Prandtl number can be as small as 10-5, approaching realistic liquid-metal values. The magnetic field destabilizes the flow, except for radial profiles… ▽ More We consider the linear stability of dissipative MHD Taylor-Couette flow with imposed toroidal magnetic fields. The inner and outer cylinders can be either insulating or conducting; the inner one rotates, the outer one is stationary. The magnetic Prandtl number can be as small as 10-5, approaching realistic liquid-metal values. The magnetic field destabilizes the flow, except for radial profiles of B$_φ$(R) close to the current-free solution. The profile with B$_{in}$=B$_{out}$ (the most uniform field) is considered in detail. For weak fields the TC-flow is stabilized, until for moderately strong fields the m=1 azimuthal mode dramatically destabilizes the flow again. There is thus a maximum value for the critical Reynolds number. For sufficiently strong fields (as measured by the Hartmann number) the toroidal field is always unstable, even for Re=0. The electric currents needed to generate the required toroidal fields in laboratory experiments are a few kA if liquid sodium is used, somewhat more if gallium is used. Weaker currents are needed for wider gaps, so a wide-gap apparatus could succeed even with gallium. The critical Reynolds numbers are only somewhat larger than the nonmagnetic values, so such an experiment would require only modest rotation rates. △ Less

Submitted 17 November, 2006; v1 submitted 15 November, 2006; originally announced November 2006.

Comments: 8 pages, 4 figures, submitted to Phys. Rev. E

The traveling wave MRI in cylindrical Taylor-Couette flow: comparing wavelengths and speeds in theory and experiment

Authors: Guenther Ruediger, Rainer Hollerbach, Frank Stefani, Thomas Gundrum, Gunter Gerbeth, Robert Rosner

Abstract: We study experimentally the flow of a liquid metal confined between differentially rotating cylinders, in the presence of externally imposed axial and azimuthal magnetic fields. For increasingly large azimuthal fields a wave-like disturbance arises, traveling along the axis of the cylinders. The wavelengths and speeds of these structures, as well as the field strengths and rotation rates at whic… ▽ More We study experimentally the flow of a liquid metal confined between differentially rotating cylinders, in the presence of externally imposed axial and azimuthal magnetic fields. For increasingly large azimuthal fields a wave-like disturbance arises, traveling along the axis of the cylinders. The wavelengths and speeds of these structures, as well as the field strengths and rotation rates at which they arise, are broadly consistent with theoretical predictions of such a traveling wave magnetorotational instability. △ Less

Submitted 18 September, 2006; v1 submitted 4 July, 2006; originally announced July 2006.

Comments: 5 pages, 3 figures, accepted by ApJ

Journal ref: Astrophys.J.649:L145-L148,2006

Experimental evidence for magnetorotational instability in a helical magnetic field

Authors: Frank Stefani, Thomas Gundrum, Gunter Gerbeth, Günther Rüdiger, Manfred Schultz, Jacek Szklarski, Rainer Hollerbach

Abstract: A recent paper [R. Hollerbach and G. Rudiger, Phys. Rev. Lett. 95, 124501 (2005)] has shown that the threshold for the onset of the magnetorotational instability (MRI) in a Taylor-Couette flow is dramatically reduced if both axial and azimuthal magnetic fields are imposed. In agreement with this prediction, we present results of a Taylor-Couette experiment with the liquid metal alloy GaInSn, sho… ▽ More A recent paper [R. Hollerbach and G. Rudiger, Phys. Rev. Lett. 95, 124501 (2005)] has shown that the threshold for the onset of the magnetorotational instability (MRI) in a Taylor-Couette flow is dramatically reduced if both axial and azimuthal magnetic fields are imposed. In agreement with this prediction, we present results of a Taylor-Couette experiment with the liquid metal alloy GaInSn, showing evidence for the existence of the MRI at Reynolds numbers of order 1000 and Hartmann numbers of order 10. △ Less

Submitted 20 June, 2006; originally announced June 2006.

Comments: 4 pages, 4 figures

Journal ref: Phys.Rev.Lett. 97 (2006) 184502

New type of magneto-rotational instability in cylindrical Taylor-Couette flow

Authors: R. Hollerbach, G. Ruediger

Abstract: We study the stability of cylindrical Taylor-Couette flow in the presence of combined axial and azimuthal magnetic fields, and show that adding an azimuthal field profoundly alters the previous results for purely axial fields. For small magnetic Prandtl numbers Pm, the critical Reynolds number Re_c for the onset of the magneto-rotational instability becomes independent of Pm, whereas for purely… ▽ More We study the stability of cylindrical Taylor-Couette flow in the presence of combined axial and azimuthal magnetic fields, and show that adding an azimuthal field profoundly alters the previous results for purely axial fields. For small magnetic Prandtl numbers Pm, the critical Reynolds number Re_c for the onset of the magneto-rotational instability becomes independent of Pm, whereas for purely axial fields it scales as Pm^{-1}. For typical liquid metals, Re_c is then reduced by several orders of magnitude, enough that this new design should succeed in realizing this instability in the laboratory. △ Less

Submitted 1 August, 2005; originally announced August 2005.

Journal ref: Phys.Rev.Lett.95:124501,2005

The stability of MHD Taylor-Couette flow with current-free spiral magnetic fields between conducting cylinders

Authors: G. Ruediger, R. Hollerbach, M. Schultz, D. A. Shalybkov

Abstract: We study the magnetorotational instability in cylindrical Taylor-Couette flow, with the (vertically unbounded) cylinders taken to be perfect conductors, and with externally imposed spiral magnetic fields. The azimuthal component of this field is generated by an axial current inside the inner cylinder, and may be slightly stronger than the axial field. We obtain an instability beyond the Rayleigh… ▽ More We study the magnetorotational instability in cylindrical Taylor-Couette flow, with the (vertically unbounded) cylinders taken to be perfect conductors, and with externally imposed spiral magnetic fields. The azimuthal component of this field is generated by an axial current inside the inner cylinder, and may be slightly stronger than the axial field. We obtain an instability beyond the Rayleigh line, for Reynolds numbers of order 1000 and Hartmann numbers of order 10, and independent of the (small) magnetic Prandtl number. For experiments with R_out = 2 R_in = 10 cm and Omega_out = 0.27 Omega_in, the instability appears for liquid sodium for axial fields of ~20 Gauss and axial currents of \~1200 A. For gallium the numbers are ~50 Gauss and ~3200 A. The vertical cell size is about twice the cell size known for nonmagnetic experiments. △ Less

Submitted 6 July, 2005; originally announced July 2005.

Comments: 5 pages, 7 figures, submitted to Astronomische Nachrichten

End-effects in rapidly rotating cylindrical Taylor-Couette flow

Authors: Rainer Hollerbach, Alexandre Fournier

Abstract: We present numerical simulations of the flow in a rapidly rotating cylindrical annulus. We show that at the rotation rates relevant to the magneto-rotational instability, the flow is strongly constrained by the Taylor-Proudman theorem. As a result, it is controlled almost entirely by the end-plates. We then consider two possible options for minimizing these end-effects, namely (i) simply taking… ▽ More We present numerical simulations of the flow in a rapidly rotating cylindrical annulus. We show that at the rotation rates relevant to the magneto-rotational instability, the flow is strongly constrained by the Taylor-Proudman theorem. As a result, it is controlled almost entirely by the end-plates. We then consider two possible options for minimizing these end-effects, namely (i) simply taking a very long cylinder, and (ii) splitting the end-plates into a series of differentially rotating rings. Regarding option (i), we show that the cylinder would have to be hundreds of times as long as it is wide before end-effects become unimportant in the interior. Since this is clearly not feasible, we turn to option (ii), and show that in order to obtain a smooth angular velocity profile, the end-plates would have to be split into around ten rings. If the end-plates are split into fewer rings, perhaps 3-5, the angular velocity profile will not be smooth, but will instead consist of a series of Stewartson layers at the boundaries from one ring to the next. We suggest therefore that the instabilities one obtains in this system will be the familiar Kelvin-Helmholtz instabilities of these Stewartson layers, rather than the magneto-rotational instability. At best, one might hope to obtain the MRI superimposed on these Kelvin-Helmholtz modes. Any subsequent interpretation of results is thus likely to be quite complicated. △ Less

Submitted 3 June, 2005; originally announced June 2005.

Journal ref: AIP Conf.Proc.733:114-121,2004