BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Toby S. Wood

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPphysics.flu-dynastro-ph.HEastro-ph.GAcond-mat.supr-connucl-thphysics.ao-phphysics.geo-phphysics.plasm-ph

Catalog footprint

What is connected

14works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

14 published item(s)

preprint2023arXiv

Shear-driven magnetic buoyancy in the solar tachocline: The mean electromotive force due to rotation

The leading theoretical paradigm for the Sun's magnetic cycle is an $αω$-dynamo process, in which a combination of differential rotation and turbulent, helical flows produces a large-scale magnetic field that reverses every 11 years. Most $αω$ solar dynamo models rely on differential rotation in the solar tachocline to generate a strong toroidal field. The most problematic part of such models is then the production of the large-scale poloidal field, via a process known as the $α$-effect. Whilst this is usually attributed to small-scale convective motions under the influence of rotation, the efficiency of this regenerative process has been called into question by some numerical simulations. Motivated by likely conditions within the tachocline, the aim of this paper is to investigate an alternative mechanism for the poloidal field regeneration, namely the magnetic buoyancy instability in a shear-generated, rotating magnetic layer. We use a local, fully compressible model in which an imposed vertical shear winds up an initially vertical magnetic field. The field ultimately becomes buoyantly unstable, and we measure the resulting mean electromotive force (EMF). For sufficiently rapid rotation, we find that a significant component of the mean EMF is aligned with the direction of the mean magnetic field, which is the characteristic feature of the classical $αω$-dynamo model. Our results therefore suggest that magnetic buoyancy could contribute directly to the generation of large-scale poloidal field in the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Three-Dimensional Magnetothermal Simulations of Magnetar Outbursts

The defining trait of magnetars, the most strongly magnetized neutron stars (NSs), is their transient activity in the X/$γ$-bands. In particular, many of them undergo phases of enhanced emission, the so-called outbursts, during which the luminosity rises by a factor $\sim$10$-$1000 in a few hours to then decay over months/years. Outbursts often exhibit a thermal spectrum, associated with the appearance of hotter regions on the surface of the star, which subsequently change in shape and cool down. Here we simulate the unfolding of a sudden, localized heat injection in the external crust of a NS with a 3D magneto-thermal evolution code, finding that this can reproduce the main features of magnetar outbursts. A full 3D treatment allows us to study for the first time the inherently asymmetric hot-spots which appear on the surface of the star as the result of the injection and to follow the evolution of their temperature and shape. We investigate the effects produced by different physical conditions in the heated region, highlighting in particular how the geometry of the magnetic field plays a key role in determining the properties of the event.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Superconducting phases in a two-component microscale model of neutron star cores

We identify the possible ground states for a mixture of two superfluid condensates (one neutral, the other electrically charged) using a phenomenological Ginzburg-Landau model. While this framework is applicable to any interacting condensed-matter mixture of a charged and a neutral component, we focus on nuclear matter in neutron star cores, where proton and neutron condensates are coupled via non-dissipative entrainment. We employ the Skyrme interaction to determine the neutron star's equilibrium composition, and hence obtain realistic coefficients for our Ginzburg-Landau model at each depth within the star's core. We then use the Ginzburg-Landau model to determine the ground state in the presence of a magnetic field. In this way, we obtain superconducting phase diagrams for six representative Skyrme models, revealing the microphysical magnetic flux distribution throughout the neutron star core. The phase diagrams are rather complex and the locations of most of the phase transitions can only be determined through numerical calculations. Nonetheless, we find that for all equations of state considered in this work, much of the outer core exhibits type-1.5 superconductivity, rather than type-II superconductivity as is generally assumed. For local magnetic field strengths $\lesssim 10^{14} \, {\rm G}$, the magnetic flux is distributed inhomogeneously, with bundles of magnetic fluxtubes separated by flux-free Meissner regions. We provide an approximate criterion to determine the transition between this type-1.5 phase and the type-I region in the inner core.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the saturation mechanism of the fluctuation dynamo at ${\text{Pr}_\mathrm{M}} \ge 1$

The presence of magnetic fields in many astrophysical objects is due to dynamo action, whereby a part of the kinetic energy is converted into magnetic energy. A turbulent dynamo that produces magnetic field structures on the same scale as the turbulent flow is known as the fluctuation dynamo. We use numerical simulations to explore the nonlinear, statistically steady state of the fluctuation dynamo in driven turbulence. We demonstrate that as the magnetic field growth saturates, its amplification and diffusion are both affected by the back-reaction of the Lorentz force upon the flow. The amplification of the magnetic field is reduced due to stronger alignment between the velocity field, magnetic field, and electric current density. Furthermore, we confirm that the amplification decreases due to a weaker stretching of the magnetic field lines. The enhancement in diffusion relative to the field line stretching is quantified by a decrease in the computed local value of the magnetic Reynolds number. Using the Minkowski functionals, we quantify the shape of the magnetic structures produced by the dynamo as magnetic filaments and ribbons in both kinematic and saturated dynamos and derive the scalings of the typical length, width, and thickness of the magnetic structures with the magnetic Reynolds number. We show that all three of these magnetic length scales increase as the dynamo saturates. The magnetic intermittency, strong in the kinematic dynamo (where the magnetic field strength grows exponentially) persists in the statistically steady state, but intense magnetic filaments and ribbons are more volume-filling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new model for mixing by double-diffusive convection (semi-convection). II. The transport of heat and composition through layers

Regions of stellar and planetary interiors that are unstable according to the Schwarzschild criterion, but stable according to the Ledoux criterion, are subject to a form of oscillatory double-diffusive (ODD) convection often called "semi-convection". In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and composition by ODD convection, and ultimately propose a new 1D prescription that can be used in stellar and planetary structure and evolution models. The first paper in this series demonstrated that under certain conditions ODD convection spontaneously transitions from an initially homogeneously turbulent state into a staircase of convective layers, which results in a substantial increase in the transport of heat and composition. Here, we present simulations of ODD convection in this layered regime, we describe the dynamical behavior of the layers, and we derive empirical scaling laws for the transport through layered convection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Transport by meridional circulations in solar-type stars

Transport by meridional flows has significant consequences for stellar evolution, but is difficult to capture in global-scale numerical simulations because of the wide range of timescales involved. Stellar evolution models therefore usually adopt parameterizations for such transport based on idealized laminar or mean-field models. Unfortunately, recent attempts to model this transport in global simulations have produced results that are not consistent with any of these idealized models. In an effort to explain the discrepancies between global simulations and idealized models, we here use three-dimensional local Cartesian simulations of compressible convection to study the efficiency of transport by meridional flows below a convection zone in several parameter regimes of relevance to the Sun and solar-type stars. In these local simulations we are able to establish the correct ordering of dynamical timescales, although the separation of the timescales remains unrealistic. We find that, even though the generation of internal waves by convective overshoot produces a high degree of time dependence in the meridional flow field, the mean flow has the qualitative behavior predicted by laminar, "balanced" models. In particular, we observe a progressive deepening, or "burrowing", of the mean circulation if the local Eddington-Sweet timescale is shorter than the viscous diffusion timescale. Such burrowing is a robust prediction of laminar models in this parameter regime, but has never been observed in any previous numerical simulation. We argue that previous simulations therefore underestimate the transport by meridional flows.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of the solar tachocline III: Numerical solutions of the Gough and McIntyre model

We present the first numerical simulations of the solar interior to exhibit a tachocline consistent with the Gough and McIntyre (1998) model. We find nonlinear, axisymmetric, steady-state numerical solutions in which: (1) a large-scale primordial field is confined within the radiation zone by downwelling meridional flows that are gyroscopically pumped in the convection zone (2) the radiation zone is in almost-uniform rotation, with a rotation rate consistent with observations (3) the bulk of the tachocline is magnetic free, in thermal-wind balance and in thermal equilibrium and (4) the interaction between the field and the flows takes place within a very thin magnetic boundary layer, the tachopause, located at the bottom of the tachocline. We show that the thickness of the tachocline scales with the amplitude of the meridional flows exactly as predicted by Gough and McIntyre. We also determine the parameter conditions under which such solutions can be obtained, and provide a simple explanation for the failure of previous numerical attempts at reproducing the Gough and McIntyre model. Finally, we discuss the implications of our findings for future numerical models of the solar interior, and for future observations of the Sun and other stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of metallic walls on dynamos generated by laminar boundary-driven flow in a spherical domain

We present a numerical study of dynamo action in a conducting fluid encased in a metallic spherical shell. Motions in the fluid are driven by differential rotation of the outer metallic shell, which we refer to as "the wall". The two hemispheres of the wall are held in counter-rotation, producing a steady, axisymmetric interior flow consisting of differential rotation and a two-cell meridional circulation with radial inflow in the equatorial plane. From previous studies, this type of flow is known to maintain a stationary equatorial dipole by dynamo action if the magnetic Reynolds number is larger than about 300 and if the outer boundary is electrically insulating. We vary independently the thickness, electrical conductivity, and magnetic permeability of the wall to determine their effect on the dynamo action. The main results are: (a) Increasing the conductivity of the wall hinders the dynamo by allowing eddy currents within the wall, which are induced by the relative motion of the equatorial dipole field and the wall. This processes can be viewed as a skin effect or, equivalently, as the tearing apart of the dipole by the differential rotation of the wall, to which the field lines are anchored by high conductivity. (b) Increasing the magnetic permeability of the wall favors dynamo action by constraining the magnetic field lines in the fluid to be normal to the wall, thereby decoupling the fluid from any induction in the wall. (c) Decreasing the wall thickness limits the amplitude of the eddy currents, and is therefore favorable for dynamo action, provided that the wall is thinner than the skin depth. We explicitly demonstrate these effects of the wall properties on the dynamo field by deriving an effective boundary condition in the limit of vanishing wall thickness.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy conservation and gravity waves in sound-proof treatments of stellar interiors: Part II Lagrangian constrained analysis

The speed of sound greatly exceeds typical flow velocities in many stellar and planetary interiors. To follow the slow evolution of subsonic motions, various sound-proof models attempt to remove fast acoustic waves whilst retaining stratified convection and buoyancy dynamics. In astrophysics, anelastic models typically receive the most attention in the class of sound-filtered stratified models. Generally, anelastic models remain valid in nearly adiabatically stratified regions like stellar convection zones, but may break down in strongly sub-adiabatic, stably stratified layers common in stellar radiative zones. However, studying stellar rotation, circulation, and dynamos requires understanding the complex coupling between convection and radiative zones, and this requires robust equations valid in both regimes. Here we extend the analysis of equation sets begun in Brown Vasil & Zweibel 2012, which studied anelastic models, to two types of pseudo-incompressible models. This class of models has received attention in atmospheric applications, and more recently in studies of white-dwarf supernovae progenitors. We demonstrate that one model conserves energy but the other does not. We use Lagrangian variational methods to extend the energy conserving model to a general equation of state, and dub the resulting equation set the Generalized Pseudo-Incompressible (GPI) model. We show that the GPI equations suitably capture low frequency phenomena in both convection and radiative zones in stars and other stratified systems, and we provide recommendations for converting low-Mach number codes to this equation set.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new model for mixing by double-diffusive convection (semi-convection): I. The conditions for layer formation

The process referred to as "semi-convection" in astrophysics and "double-diffusive convection in the diffusive regime" in Earth and planetary sciences, occurs in stellar and planetary interiors in regions which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. In this series of papers, we analyze the results of an extensive suite of 3D numerical simulations of the process, and ultimately propose a new 1D prescription for heat and compositional transport in this regime which can be used in stellar or planetary structure and evolution models. In a preliminary study of the phenomenon, Rosenblum et al. (2011) showed that, after saturation of the primary instability, a system can evolve in one of two possible ways: the induced turbulence either remains homogeneous, with very weak transport properties, or transitions into a thermo-compositional staircase where the transport rate is much larger (albeit still smaller than in standard convection). In this paper, we show that this dichotomous behavior is a robust property of semi-convection across a wide region of parameter space. We propose a simple semi-analytical criterion to determine whether layer formation is expected or not, and at what rate it proceeds, as a function of the background stratification and of the diffusion parameters (viscosity, thermal diffusivity and compositional diffusivity) only. The theoretical criterion matches the outcome of our numerical simulations very adequately in the numerically accessible "planetary" parameter regime, and can easily be extrapolated to the stellar parameter regime. Subsequent papers will address more specifically the question of quantifying transport in the layered case and in the non-layered case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Confinement of the Sun's interior magnetic field, with implications for lithium burning

The simplest interior magnetic field B_i that can explain the observed uniform rotation of the Sun's radiative envelope is an axial dipole stabilized by a deep toroidal field. It can explain the uniform rotation only if confined in the polar caps. The field must be prevented from diffusing up into the high-latitude convection zone, whose slower rotation must remain decoupled from the radiative interior. This paper describes new analytical and numerical solutions of the relevant magnetohydrodynamic equations showing that such confinement and decoupling is dynamically possible by means of a laminar "magnetic confinement layer" at the bottom of the tachocline. With realistic values of the microscopic diffusivities, a weak laminar downwelling flow U~10^{-5}cm/s over the poles is enough to enforce exponential decay of B_i with altitude, in a confinement layer only a fraction of a megameter thick. Downwelling in the polar tachocline is implied both by helioseismic observations, combined with elementary dynamics, and by theoretical arguments about the "gyroscopic pumping" that would spread differential rotation downward in the absence of B_i. Our confinement-layer solutions are the first to take account of all the relevant physical effects in a self-consistent mathematical model. The effects include magnetic diffusion, baroclinicity and stable stratification (thermal and compositional), Coriolis effects, and thermal relaxation. We discuss how the confinement layers at each pole might fit into a global dynamical picture of the solar tachocline. That picture, in turn, suggests a new insight into the early Sun and solar lithium depletion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polar confinement of the Sun's interior magnetic field by laminar magnetostrophic flow

The global-scale interior magnetic field needed to account for the Sun's observed differential rotation can be effective only if confined below the convection zone in all latitudes, including the polar caps. Axisymmetric nonlinear MHD solutions are obtained showing that such confinement can be brought about by a very weak downwelling flow U~10^{-5}cm/s over each pole. Such downwelling is consistent with the helioseismic evidence. All three components of the magnetic field decay exponentially with altitude across a thin "magnetic confinement layer" located at the bottom of the tachocline. With realistic parameter values, the thickness of the confinement layer ~10^{-3} of the Sun's radius. Alongside baroclinic effects and stable thermal stratification, the solutions take into account the stable compositional stratification of the helium settling layer, if present as in today's Sun, and the small diffusivity of helium through hydrogen, chi. The small value of chi relative to magnetic diffusivity produces a double boundary-layer structure in which a "helium sublayer" of smaller vertical scale is sandwiched between the top of the helium settling layer and the rest of the confinement layer. Solutions are obtained using both semi-analytical and purely numerical, finite-difference techniques. The confinement-layer flows are magnetostrophic to excellent approximation. More precisely, the principal force balances are between Lorentz, Coriolis, pressure-gradient and buoyancy forces, with relative accelerations and viscous forces negligible. This is despite the kinematic viscosity being somewhat greater than chi. We discuss how the confinement layers at each pole might fit into a global dynamical picture of the solar tachocline. That picture, in turn, suggests a new insight into the early Sun and into the longstanding enigma of solar lithium depletion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Sun's meridional circulation and interior magnetic field

To date, no self-consistent numerical simulation of the solar interior has succeeded in reproducing the observed thinness of the solar tachocline, and the persistence of uniform rotation beneath it. Although it is known that the uniform rotation can be explained by the presence of a global-scale confined magnetic field, numerical simulations have thus far failed to produce any solution where such a field remains confined against outward diffusion. We argue that the problem lies in the choice of parameters for which these numerical simulations have been performed. We construct a simple analytical magneto-hydrodynamic model of the solar interior and identify several distinct parameter regimes. For realistic solar parameter values, our results are in broad agreement with the tachocline model of Gough & McIntyre. In this regime, meridional flows driven at the base of the convection zone are of sufficient amplitude to hold back the interior magnetic field against diffusion. For the parameter values used in existing numerical simulations, on the other hand, we find that meridional flows are significantly weaker and, we argue, unable to confine the interior field. We propose a method for selecting parameter values in future numerical models.

0 citations0 reviewsNo code linkedOpen access