-
Tayler-Spruit dynamo in stably stratified rotating fluids: Application to proto-magnetars
Authors:
Paul Barrère,
Jérôme Guilet,
Raphaël Raynaud,
Alexis Reboul-Salze
Abstract:
The formation of highly magnetized young neutron stars, called magnetars, is still a strongly debated question. A promising scenario invokes the amplification of the magnetic field by the Tayler-Spruit dynamo in a proto-neutron star (PNS) spun up by fallback. Barrère et al. 2023 supports this scenario by demonstrating that this dynamo can generate magnetar-like magnetic fields in stably stratified…
▽ More
The formation of highly magnetized young neutron stars, called magnetars, is still a strongly debated question. A promising scenario invokes the amplification of the magnetic field by the Tayler-Spruit dynamo in a proto-neutron star (PNS) spun up by fallback. Barrère et al. 2023 supports this scenario by demonstrating that this dynamo can generate magnetar-like magnetic fields in stably stratified Boussinesq models of a PNS interior. To further investigate the Tayler-Spruit dynamo, we perform 3D-MHD numerical simulations with the MagIC code varying the ratio between the Brunt-Väisälä frequency and the rotation rate. We first demonstrate that a self-sustained dynamo process can be maintained for a Brunt-Väisälä frequency about 4 times higher than the angular rotation frequency. The generated magnetic fields and angular momentum transport follow the analytical scaling laws of Fuller et al. 2019, which confirms our previous results. We also report for the first time the existence of an intermittent Tayler-Spruit dynamo. For a typical PNS Brunt-Väisälä frequency of $10^{3}\,{\rm s}^{-1}$, the axisymmetric toroidal and dipolar magnetic fields range between $1.2\times 10^{15}-2\times 10^{16}\,{\rm G}$ and $1.4\times 10^{13}-3\times 10^{15}\,{\rm G}$, for rotation periods of $1-10\,{\rm ms}$. Thus, our results provide numerical evidence that our scenario can explain the formation of magnetars. As the Tayler-Spruit dynamo is often invoked for the angular momentum transport in stellar radiative zones, our results are also of particular importance in this field and we provide a calibration of the Fuller et al.'s prescription based on our simulations, with a dimensionless normalisation factor of the order of $10^{-2}$.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Coriolis darkening in late-type stars II. Effect of self-sustained magnetic fields in stratified convective envelope
Authors:
C. Pinçon,
L. Petitdemange,
R. Raynaud,
L. J. Garcia,
A. Guseva,
M. Rieutord,
E. Alecian
Abstract:
Modeling the surface brightness distribution of stars is of prime importance to interpret observations. Nevertheless, this remains quite challenging for cool stars as it requires one to model the MHD turbulence that develops in their convective envelope. In Paper I, the effect of the Coriolis acceleration on the surface heat flux has been studied by means of hydrodynamic simulations. In this paper…
▽ More
Modeling the surface brightness distribution of stars is of prime importance to interpret observations. Nevertheless, this remains quite challenging for cool stars as it requires one to model the MHD turbulence that develops in their convective envelope. In Paper I, the effect of the Coriolis acceleration on the surface heat flux has been studied by means of hydrodynamic simulations. In this paper, we aim to investigate the additional effect of dynamo magnetic fields. We focus on an envelope thickness that is representative of either a $\sim0.35~M_\odot$ M dwarf, a young red giant star or a pre-main sequence star. We performed a parametric study using numerical MHD simulations of anelastic convection in thick rotating spherical shells. For each model, we computed the mean surface distribution of the heat flux, and examined the leading-order effect of the magnetic field on the obtained latitudinal luminosity profile. We identify three different regimes. Close to the onset of convection, while the first unstable modes tend to convey heat more efficiently near the equator, magnetic fields are shown to generally enhance the mean heat flux close to the polar regions (and the tangent cylinder). By progressively increasing the Rayleigh number, the development of a prograde equatorial jet was previously shown to make the equator darker when no magnetic field is taken into account. For moderate Rayleigh numbers, magnetic fields can instead inverse the mean pole-equator brightness contrast (which means going from a darker to a brighter equator when a dynamo sets in) and finally induce a similar regime to that found close to the onset of convection. For more turbulent models with larger Rayleigh numbers, magnetic fields alternatively tend to smooth out the brightness contrast. This general behavior is shown to be related to the quenching of the surface differential rotation by magnetic fields.
△ Less
Submitted 13 March, 2024;
originally announced March 2024.
-
Prospects for realtime characterization of core-collapse supernova and neutrino properties
Authors:
Meriem Bendahman,
Isabel Goos,
Joao Coelho,
Matteo Bugli,
Alexis Coleiro,
Sonia El Hedri,
Thierry Foglizzo,
Davide Franco,
Jérôme Guilet,
Antoine Kouchner,
Raphaël Raynaud,
Yahya Tayalati
Abstract:
Core-collapse supernovae (CCSNe) offer extremely valuable insights into the dynamics of galaxies. Neutrino time profiles from CCSNe, in particular, could reveal unique details about collapsing stars and particle behavior in dense environments. However, CCSNe in our galaxy and the Large Magellanic Cloud are rare and only one supernova neutrino observation has been made so far. To maximize the infor…
▽ More
Core-collapse supernovae (CCSNe) offer extremely valuable insights into the dynamics of galaxies. Neutrino time profiles from CCSNe, in particular, could reveal unique details about collapsing stars and particle behavior in dense environments. However, CCSNe in our galaxy and the Large Magellanic Cloud are rare and only one supernova neutrino observation has been made so far. To maximize the information obtained from the next Galactic CCSN, it is essential to combine analyses from multiple neutrino experiments in real time and transmit any relevant information to electromagnetic facilities within minutes. Locating the CCSN, in particular, is challenging, requiring disentangling CCSN localization information from observational features associated with the properties of the supernova progenitor and the physics of the neutrinos. Yet, being able to estimate the progenitor distance from the neutrino signal would be of great help for the optimisation of the electromagnetic follow-up campaign that will start soon after the propagation of the neutrino alert. Existing CCSN distance measurement algorithms based on neutrino observations hence rely on the assumption that neutrino properties can be described by the Standard Model. This paper presents a swift and robust approach to extract CCSN and neutrino physics information, leveraging diverse next-generation neutrino detectors to counteract potential measurement biases from Beyond the Standard Model effects.
△ Less
Submitted 23 November, 2023; v1 submitted 10 November, 2023;
originally announced November 2023.
-
Numerical simulations of the Tayler-Spruit dynamo in proto-magnetars
Authors:
Paul Barrère,
Jérôme Guilet,
Raphaël Raynaud,
Alexis Reboul-Salze
Abstract:
The Tayler-Spruit dynamo is one of the most promising mechanisms proposed to explain angular momentum transport during stellar evolution. Its development in proto-neutron stars spun-up by supernova fallback has also been put forward as a scenario to explain the formation of very magnetized neutron stars called magnetars. Using three-dimensional direct numerical simulations, we model the proto-neut…
▽ More
The Tayler-Spruit dynamo is one of the most promising mechanisms proposed to explain angular momentum transport during stellar evolution. Its development in proto-neutron stars spun-up by supernova fallback has also been put forward as a scenario to explain the formation of very magnetized neutron stars called magnetars. Using three-dimensional direct numerical simulations, we model the proto-neutron star interior as a stably stratified spherical Couette flow with the outer sphere that rotates faster than the inner one. We report the existence of two subcritical dynamo branches driven by the Tayler instability. They differ by their equatorial symmetry (dipolar or hemispherical) and the magnetic field scaling, which is in agreement with different theoretical predictions (by Fuller and Spruit, respectively). The magnetic dipole of the dipolar branch is found to reach intensities compatible with observational constraints on magnetars.
△ Less
Submitted 26 August, 2023; v1 submitted 21 June, 2023;
originally announced June 2023.
-
A new scenario for magnetar formation: Tayler-Spruit dynamo in a proto-neutron star spun up by fallback
Authors:
P. Barrère,
J. Guilet,
A. Reboul-Salze,
R. Raynaud,
H. -T. Janka
Abstract:
Magnetars are isolated young neutron stars characterized by the most intense magnetic fields known in the universe. The origin of their magnetic field is still a challenging question. In situ magnetic field amplification by dynamo action is a promising process to generate ultra-strong magnetic fields in fast-rotating progenitors. However, it is unclear whether the fraction of progenitors harboring…
▽ More
Magnetars are isolated young neutron stars characterized by the most intense magnetic fields known in the universe. The origin of their magnetic field is still a challenging question. In situ magnetic field amplification by dynamo action is a promising process to generate ultra-strong magnetic fields in fast-rotating progenitors. However, it is unclear whether the fraction of progenitors harboring fast core rotation is sufficient to explain the entire magnetar population. To address this point, we propose a new scenario for magnetar formation, in which a slow-rotating proto-neutron star is spun up by the supernova fallback. We argue that this can trigger the development of the Tayler-Spruit dynamo while other dynamo processes are disfavored. Using previous works done on this dynamo and simulations to characterize the fallback, we derive equations modelling the coupled evolution of the proto-neutron star rotation and magnetic field. Their time integration for different fallback masses is successfully compared with analytical estimates of the amplification timescales and saturation value of the magnetic field. We find that the magnetic field is amplified within $20$ to $40$s after the core bounce, and that the radial magnetic field saturates at intensities $10^{14}-10^{15}$G, therefore spanning the full range of magnetar's dipolar magnetic fields. We also compare predictions of two proposed saturation mechanisms showing that magnetar-like magnetic fields can be generated for neutron star spun up to rotation periods $\lesssim8$ms and $\lesssim28$ms, corresponding to fallback masses $\gtrsim4\times10^{-2}{\rm M}_{\odot}$ and $\gtrsim10^{-2}{\rm M}_{\odot}$. Thus, our results suggest that magnetars can be formed from slow-rotating progenitors for fallback masses compatible with recent supernova simulations and leading to plausible initial rotation periods of the proto-neutron star.
△ Less
Submitted 19 October, 2022; v1 submitted 2 June, 2022;
originally announced June 2022.
-
MRI-driven dynamo at very high magnetic Prandtl numbers
Authors:
Jerome Guilet,
Alexis Reboul-Salze,
Raphael Raynaud,
Matteo Bugli,
Basile Gallet
Abstract:
The dynamo driven by the magnetorotational instability (MRI) is believed to play an important role in the dynamics of accretion discs and may also explain the origin of the extreme magnetic fields present in magnetars. Its saturation level is an important open question known to be particularly sensitive to the diffusive processes through the magnetic Prandtl number Pm (the ratio of viscosity to re…
▽ More
The dynamo driven by the magnetorotational instability (MRI) is believed to play an important role in the dynamics of accretion discs and may also explain the origin of the extreme magnetic fields present in magnetars. Its saturation level is an important open question known to be particularly sensitive to the diffusive processes through the magnetic Prandtl number Pm (the ratio of viscosity to resistivity). Despite its relevance to proto-neutron stars and neutron star merger remnants, the numerically challenging regime of high Pm is still largely unknown. Using zero-net flux shearing box simulations in the incompressible approximation, we studied MRI-driven dynamos at unprecedentedly high values of Pm reaching 256. The simulations show that the stress and turbulent energies are proportional to Pm up to moderately high values ($\mathrm{Pm} \sim 50$). At higher Pm, they transition to a new regime consistent with a plateau independent of Pm for $\rm Pm \gtrsim 100$. This trend is independent of the Reynolds number, which may suggest an asymptotic regime where the energy injection and dissipation are independent of the diffusive processes. Interestingly, large values of Pm not only lead to intense small-scale magnetic fields but also to a more efficient dynamo at the largest scales of the box.
△ Less
Submitted 1 September, 2022; v1 submitted 17 May, 2022;
originally announced May 2022.
-
MRI-driven $α-Ω$ dynamos in protoneutron stars
Authors:
Alexis Reboul-Salze,
Jérôme Guilet,
Raphaël Raynaud,
Matteo Bugli
Abstract:
Magnetars are highly magnetized neutron stars that can produce X-ray and soft gamma-ray emissions and that have a dipole of $10^{14}$ G to $10^{15}$ G. A promising mechanism for explaining magnetar formation is magnetic field amplification by the MRI in fast-rotating protoneutron stars (PNS). This scenario is supported by recent global models, which showed that small-scale turbulence can generate…
▽ More
Magnetars are highly magnetized neutron stars that can produce X-ray and soft gamma-ray emissions and that have a dipole of $10^{14}$ G to $10^{15}$ G. A promising mechanism for explaining magnetar formation is magnetic field amplification by the MRI in fast-rotating protoneutron stars (PNS). This scenario is supported by recent global models, which showed that small-scale turbulence can generate a dipole with magnetar-like intensity. However, the impact of buoyancy and density stratification on the efficiency of the MRI in generating a dipole is still unknown. We assess the impact of the density and entropy profiles on the MRI dynamo in a global model of a fast-rotating PNS, which focuses on its outer stratified region that is stable to convection. Using the pseudo-spectral code MagIC, we performed 3D Boussinesq and anelastic MHD simulations in spherical geometry with explicit diffusivities. We performed a parameter study in which we investigate the effect of different approximations and of thermal diffusion. We obtain a self-sustained turbulent MRI-driven dynamo. This confirms most of our previous incompressible results once rescaled for density. The MRI also generates a nondominant equatorial dipole, which represents about 4.3% of the averaged magnetic field strength. Interestingly, in the presence of a density gradient, an axisymmetric magnetic field at large scales oscillates with time, which can be described as a mean-field $αΩ$ dynamo. Buoyancy damps turbulence in the equatorial plane but it has overall a relatively weak influence with a realistic high thermal diffusion. Our results support the ability of the MRI to generate magnetar-like large-scale magnetic fields. They furthermore predict the presence of an $αΩ$ dynamo in the protoneutron star, which could be important to model in-situ magnetic field amplification in core-collapse supernovae. [abridged]
△ Less
Submitted 21 September, 2022; v1 submitted 3 November, 2021;
originally announced November 2021.
-
Gravitational wave signature of proto-neutron star convection: I. MHD numerical simulations
Authors:
Raphaël Raynaud,
Pablo Cerdá-Durán,
Jérôme Guilet
Abstract:
Gravitational waves provide a unique and powerful opportunity to constrain the dynamics in the interior of proto-neutron stars during core collapse supernovae. Convective motions play an important role in generating neutron stars magnetic fields, which could explain magnetar formation in the presence of fast rotation. We compute the gravitational wave emission from proto-neutron star convection an…
▽ More
Gravitational waves provide a unique and powerful opportunity to constrain the dynamics in the interior of proto-neutron stars during core collapse supernovae. Convective motions play an important role in generating neutron stars magnetic fields, which could explain magnetar formation in the presence of fast rotation. We compute the gravitational wave emission from proto-neutron star convection and its associated dynamo, by post-processing three-dimensional MHD simulations of a model restricted to the convective zone in the anelastic approximation. We consider two different proto-neutron star structures representative of early times (with a convective layer) and late times (when the star is almost entirely convective). In the slow rotation regime, the gravitational wave emission follows a broad spectrum peaking at about three times the turnover frequency. In this regime, the inclusion of magnetic fields slightly decreases the amplitude without changing the spectrum significantly compared to a non-magnetised simulation. Fast rotation changes both the amplitude and spectrum dramatically. The amplitude is increased by a factor of up to a few thousands. The spectrum is characterized by several peaks associated to inertial modes, whose frequency scales with the rotation frequency. Using simple physical arguments, we derive scalings that reproduce quantitatively several aspects of these numerical results. We also observe an excess of low-frequency gravitational waves, which appears at the transition to a strong field dynamo characterized by a strong axisymmetric toroidal magnetic field. This signature of dynamo action could be used to constrain the dynamo efficiency in a proto-neutron star with future gravitational wave detections.
△ Less
Submitted 22 October, 2021; v1 submitted 23 March, 2021;
originally announced March 2021.
-
A global model of the magnetorotational instability in protoneutron stars
Authors:
Alexis Reboul-Salze,
Jérôme Guilet,
Raphaël Raynaud,
Matteo Bugli
Abstract:
Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MR…
▽ More
Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MRI at generating a dipole field is still unknown. To answer this question, we study the MRI dynamo in an idealized global model of a fast rotating protoneutron star with differential rotation. We perform 3D incompressible MHD simulations in spherical geometry with explicit diffusivities where the differential rotation is forced at the outer boundary. We vary the initial magnetic field and investigated different magnetic boundary conditions. These simulations were compared to local shearing box simulations. We obtain a self-sustained turbulent MRI-driven dynamo, whose saturated state is independent of the initial magnetic field. The MRI generates a strong turbulent magnetic field of $B \geq 2\times 10^{15}$ G and a non-dominant magnetic dipole, which represents systematically about $5\%$ of the averaged magnetic field strength. Interestingly, this dipole is tilted towards the equatorial plane. We find that local shearing box models can reproduce fairly well several characteristics of global MRI turbulence such as the kinetic and magnetic spectra. The turbulence is nonetheless more vigorous in the local models than in the global ones. Overall, our results support the ability of the MRI to form magnetar-like large-scale magnetic fields. They furthermore predict the presence of a stronger small-scale magnetic field. The resulting magnetic field could be important to power outstanding stellar explosions, such as superluminous supernovae and GRBs.
△ Less
Submitted 27 November, 2020; v1 submitted 7 May, 2020;
originally announced May 2020.
-
Magnetar formation through a convective dynamo in protoneutron stars
Authors:
Raphaël Raynaud,
Jérôme Guilet,
Hans-Thomas Janka,
Thomas Gastine
Abstract:
The release of spin-down energy by a magnetar is a promising scenario to power several classes of extreme explosive transients. However, it lacks a firm basis because magnetar formation still represents a theoretical challenge. Using the first three-dimensional simulations of a convective dynamo based on a protoneutron star interior model, we demonstrate that the required dipolar magnetic field ca…
▽ More
The release of spin-down energy by a magnetar is a promising scenario to power several classes of extreme explosive transients. However, it lacks a firm basis because magnetar formation still represents a theoretical challenge. Using the first three-dimensional simulations of a convective dynamo based on a protoneutron star interior model, we demonstrate that the required dipolar magnetic field can be consistently generated for sufficiently fast rotation rates. The dynamo instability saturates in the magnetostrophic regime with the magnetic energy exceeding the kinetic energy by a factor of up to 10. Our results are compatible with the observational constraints on galactic magnetar field strength and provide strong theoretical support for millisecond protomagnetar models of gamma-ray burst and superluminous supernova central engines.
△ Less
Submitted 14 March, 2020;
originally announced March 2020.
-
Gravity darkening in late-type stars. The Coriolis effect
Authors:
Raphaël Raynaud,
Michel Rieutord,
Ludovic Petitdemange,
Thomas Gastine,
Bertrand Putigny
Abstract:
Recent interferometric data have been used to constrain the brightness distribution at the surface of nearby stars, in particular the so-called gravity darkening that makes fast rotating stars brighter at their poles than at their equator. However, good models of gravity darkening are missing for stars that posses a convective envelope. In order to better understand how rotation affects the heat t…
▽ More
Recent interferometric data have been used to constrain the brightness distribution at the surface of nearby stars, in particular the so-called gravity darkening that makes fast rotating stars brighter at their poles than at their equator. However, good models of gravity darkening are missing for stars that posses a convective envelope. In order to better understand how rotation affects the heat transfer in stellar convective envelopes, we focus on the heat flux distribution in latitude at the outer surface of numerical models. We carry out a systematic parameter study of three-dimensional, direct numerical simulations of anelastic convection in rotating spherical shells. Restricting our investigations to hydrodynamical models with a thermal Prandtl number fixed to unity, we consider both thick and thin (solar-like) shells, and vary the stratification over three orders of magnitude. We measure the heat transfer efficiency in terms of the Nusselt number, defined as the output luminosity normalised by the conductive state luminosity. We report diverse Nusselt number profiles in latitude, ranging from brighter (usually at the onset of convection) to darker equator and uniform profiles. We find that the variations of the surface brightness are mainly controlled by the surface value of the local Rossby number: when the Coriolis force dominates the dynamics, the heat flux is weakened in the equatorial region by the zonal wind and enhanced at the poles by convective motions inside the tangent cylinder. In the presence of a strong background density stratification however, as expected in real stars, the increase of the local Rossby number in the outer layers leads to uniformisation of the surface heat flux distribution.
△ Less
Submitted 8 December, 2017;
originally announced December 2017.
-
Convective dynamo action in a spherical shell: symmetries and modulation
Authors:
Raphaël Raynaud,
Steven M. Tobias
Abstract:
We consider dynamo action driven by three-dimensional rotating anelastic convection in a spherical shell. Motivated by the behaviour of the solar dynamo, we examine the interaction of hydromagnetic modes with different symmetries and demonstrate how complicated interactions between convection, differential rotation and magnetic fields may lead to modulation of the basic cycle. For some parameters,…
▽ More
We consider dynamo action driven by three-dimensional rotating anelastic convection in a spherical shell. Motivated by the behaviour of the solar dynamo, we examine the interaction of hydromagnetic modes with different symmetries and demonstrate how complicated interactions between convection, differential rotation and magnetic fields may lead to modulation of the basic cycle. For some parameters, Type 1 modulation occurs by the transfer of energy between modes of different symmetries with little change in the overall amplitude, for other parameters, the modulation is of Type 2, where the amplitude is significantly affected (leading to grand minima in activity) without significant changes in symmetry. Most importantly, we identify the presence of "supermodulation" in the solutions, where the activity switches chaotically between Type 1 and Type 2 modulation, this is believed to be an important process in solar activity.
△ Less
Submitted 4 July, 2016;
originally announced July 2016.
-
Dipolar dynamos in stratified systems
Authors:
R. Raynaud,
L. Petitdemange,
E. Dormy
Abstract:
Observations of low-mass stars reveal a variety of magnetic field topologies ranging from large-scale, axial dipoles to more complex magnetic fields. At the same time, three-dimensional spherical simulations of convectively driven dynamos reproduce a similar diversity, which is commonly obtained either with Boussinesq models or with more realistic models based on the anelastic approximation, whi…
▽ More
Observations of low-mass stars reveal a variety of magnetic field topologies ranging from large-scale, axial dipoles to more complex magnetic fields. At the same time, three-dimensional spherical simulations of convectively driven dynamos reproduce a similar diversity, which is commonly obtained either with Boussinesq models or with more realistic models based on the anelastic approximation, which take into account the variation of the density with depth throughout the convection zone. Nevertheless, a conclusion from different anelastic studies is that dipolar solutions seem more difficult to obtain as soon as substantial stratifications are considered. In this paper, we aim at clarifying this point by investigating in more detail the influence of the density stratification on dipolar dynamos. To that end, we rely on a systematic parameter study that allows us to clearly follow the evolution of the stability domain of the dipolar branch as the density stratification is increased. The impact of the density stratification both on the dynamo onset and the dipole collapse is discussed and compared to previous Boussinesq results. Furthermore, our study indicates that the loss of the dipolar branch does not ensue from a specific modification of the dynamo mechanisms related to the background stratification, but could instead result from a bias as our observations naturally favour a certain domain in the parameter space characterized by moderate values of the Ekman number, owing to current computational limitations. Moreover, we also show that the critical magnetic Reynolds number of the dipolar branch is scarcely modified by the increase of the density stratification, which provides an important insight into the global understanding of the impact of the density stratification on the stability domain of the dipolar dynamo branch.
△ Less
Submitted 3 March, 2015; v1 submitted 28 February, 2015;
originally announced March 2015.
-
Influence of the mass distribution on the magnetic field topology
Authors:
Raphaël Raynaud,
Ludovic Petitdemange,
Emmanuel Dormy
Abstract:
Context. Three-dimensional spherical dynamo simulations carried out within the framework of the anelastic approximation have revealed that the established distinction between dipolar and multipolar dynamos tends to be less clear than it was in Boussinesq studies. This result was first interpreted as a direct consequence of the existence of a larger number of models with a high equatorial dipole co…
▽ More
Context. Three-dimensional spherical dynamo simulations carried out within the framework of the anelastic approximation have revealed that the established distinction between dipolar and multipolar dynamos tends to be less clear than it was in Boussinesq studies. This result was first interpreted as a direct consequence of the existence of a larger number of models with a high equatorial dipole contribution, together with an intermediate dipole field strength. However, this finding has not been clearly related to specific changes that would have been introduced by the use of the anelastic approximation. Aims. In this paper, we primarily focus on the effects of choosing of a different mass distribution. Indeed, it is likely to have as large consequences as would taking a stratified reference state into account, especially when comparing our results to previous Boussinesq studies. Methods. Our investigation is based on a systematic parameter study of weakly stratified anelastic dynamo models. Results. We show that the tendencies highlighted in previous anelastic dynamo simulations are already present in the Boussinesq limit. Thus they cannot be systematically related to anelastic effects. Actually, a central mass distribution can result in changes in the magnetic field topology that are mainly due to the concentration of convective cells close to the inner sphere.
△ Less
Submitted 8 July, 2014; v1 submitted 18 June, 2014;
originally announced June 2014.
-
Topology and field strength in spherical, anelastic dynamo simulations
Authors:
M. Schrinner,
L. Petitdemange,
R. Raynaud,
E. Dormy
Abstract:
Numerical modelling of convection driven dynamos in the Boussinesq approximation revealed fundamental characteristics of the dynamo-generated magnetic fields and the fluid flow. Because these results were obtained for an incompressible fluid, their validity for gas planets and stars remains to be assessed. A common approach is to take some density stratification into account with the so-called a…
▽ More
Numerical modelling of convection driven dynamos in the Boussinesq approximation revealed fundamental characteristics of the dynamo-generated magnetic fields and the fluid flow. Because these results were obtained for an incompressible fluid, their validity for gas planets and stars remains to be assessed. A common approach is to take some density stratification into account with the so-called anelastic approximation. The validity of previous results obtained in the Boussinesq approximation is tested for anelastic models. We point out and explain specific differences between both types of models, in particular with respect to the field geometry and the field strength, but we also compare scaling laws for the velocity amplitude, the magnetic dissipation time, and the convective heat flux. Our investigation is based on a systematic parameter study of spherical dynamo models in the anelastic approximation. We make use of a recently developed numerical solver and provide results for the test cases of the anelastic dynamo benchmark. The dichotomy of dipolar and multipolar dynamos identified in Boussinesq simulations is also present in our sample of anelastic models. Dipolar models require that the typical length scale of convection is an order of magnitude larger than the Rossby radius. However, the distinction between both classes of models is somewhat less explicit than in previous studies. This is mainly due to two reasons: we found a number of models with a considerable equatorial dipole contribution and an intermediate overall dipole field strength. Furthermore, a large density stratification may hamper the generation of dipole dominated magnetic fields. Previously proposed scaling laws, such as those for the field strength, are similarly applicable to anelastic models. It is not clear, however, if this consistency necessarily implies similar dynamo processes in both settings.
△ Less
Submitted 2 March, 2015; v1 submitted 27 December, 2013;
originally announced December 2013.
