-
Viscous Dissipation and Dynamics in Simulations of Rotating, Stratified Plane-layer Convection
Authors:
Simon R. W. Lance,
Laura K. Currie,
Matthew K. Browning
Abstract:
Convection in stars and planets must be maintained against viscous and Ohmic dissipation. Here, we present the first systematic investigation of viscous dissipation in simulations of rotating, density-stratified plane layers of convection. Our simulations consider an anelastic ideal gas, and employ the open-source code Dedalus. We demonstrate that when the convection is sufficiently vigorous, the…
▽ More
Convection in stars and planets must be maintained against viscous and Ohmic dissipation. Here, we present the first systematic investigation of viscous dissipation in simulations of rotating, density-stratified plane layers of convection. Our simulations consider an anelastic ideal gas, and employ the open-source code Dedalus. We demonstrate that when the convection is sufficiently vigorous, the integrated dissipative heating tends towards a value that is independent of viscosity or thermal diffusivity, but depends on the imposed luminosity and the stratification. We show that knowledge of the dissipation provides a bound on the magnitude of the kinetic energy flux in the convection zone. In our non-rotating cases with simple flow fields, much of the dissipation occurs near the highest possible temperatures, and the kinetic energy flux approaches this bound. In the rotating cases, although the total integrated dissipation is similar, it is much more uniformly distributed (and locally balanced by work against the stratification), with a consequently smaller kinetic energy flux. The heat transport in our rotating simulations is in good agreement with results previously obtained for 3D Boussinesq convection, and approaches the predictions of diffusion-free theory.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
Simulations of solar and stellar dynamos and their theoretical interpretation
Authors:
Petri J. Käpylä,
Matthew K. Browning,
Allan Sacha Brun,
Gustavo Guerrero,
Jörn Warnecke
Abstract:
We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scal…
▽ More
We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scale dynamo. We continue to studies that model the Sun at different ages and to studies of stars of different masses and evolutionary stages. Both simulations and observations indicate that rotation, measured by the Rossby number which is the ratio of rotation period and convective turnover time, is a key ingredient in setting the overall level and characteristics of magnetic activity. Finally, efforts to understand global 3D simulations in terms of mean-field dynamo theory are discussed.
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
Convection with Misaligned Gravity and Rotation: Simulations and Rotating Mixing Length Theory
Authors:
Laura K. Currie,
Adrian J. Barker,
Yoram Lithwick,
Matthew K. Browning
Abstract:
We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes,…
▽ More
We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.
△ Less
Submitted 6 February, 2020;
originally announced February 2020.
-
High-Energy Photon and Particle Effects onExoplanet Atmospheres and Habitability
Authors:
Jeremy J. Drake,
Julián D. Alvarado-Gómez,
Vladimir Airapetian,
P. Wilson Cauley,
Costanza Argiroffi,
Matthew K. Browning,
Damian J. Christian,
Ofer Cohen,
Lia Corrales,
William Danchi,
Miguel de Val-Borro,
Chuanfei Dong,
William Forman,
Kevin France,
Elena Gallo,
Katherine Garcia-Sage,
Cecilia Garraffo,
Dawn M. Gelino,
Guillaume Gronoff,
H. Moritz Günther,
Graham M. Harper,
Raphaëlle D. Haywood,
Margarita Karovska,
Vinay Kashyap,
Joel Kastner
, et al. (19 additional authors not shown)
Abstract:
It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current inst…
▽ More
It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.
△ Less
Submitted 28 March, 2019;
originally announced March 2019.
-
The Radius and Entropy of a Magnetized, Rotating Fully-convective Star: Analysis With Depth-dependent Mixing Length Theories
Authors:
Lewis G. Ireland,
Matthew K. Browning
Abstract:
Some low-mass stars appear to have larger radii than predicted by standard 1D structure models; prior work has suggested that inefficient convective heat transport, due to rotation and/or magnetism, may ultimately be responsible. We examine this issue using 1D stellar models constructed using Modules for Experiments in Stellar Astrophysics (MESA). First, we consider standard models that do not exp…
▽ More
Some low-mass stars appear to have larger radii than predicted by standard 1D structure models; prior work has suggested that inefficient convective heat transport, due to rotation and/or magnetism, may ultimately be responsible. We examine this issue using 1D stellar models constructed using Modules for Experiments in Stellar Astrophysics (MESA). First, we consider standard models that do not explicitly include rotational/magnetic effects, with convective inhibition modeled by decreasing a depth-independent mixing length theory (MLT) parameter $α_{\text{MLT}}$ (following Cox et al. 1981; Chabrier et al. 2007). We provide formulae linking changes in $α_{\text{MLT}}$ to changes in the interior specific entropy, and hence to the stellar radius. Next, we modify the MLT formulation in MESA to mimic explicitly the influence of rotation and magnetism, using formulations suggested by Stevenson (1979) and MacDonald & Mullan (2014) respectively. We find rapid rotation in these models has a negligible impact on stellar structure, primarily because a star's adiabat, and hence its radius, is predominantly affected by layers near the surface; convection is rapid and largely uninfluenced by rotation there. Magnetic fields, if they influenced convective transport in the manner described by MacDonald & Mullan (2014), could lead to more noticeable radius inflation. Finally, we show that these non-standard effects on stellar structure can be fabricated using a depth-dependent $α_{\text{MLT}}$: a non-magnetic, non-rotating model can be produced that is virtually indistinguishable from one that explicitly parameterizes rotation and/or magnetism using the two formulations above. We provide formulae linking the radially-variable $α_{\text{MLT}}$ to these putative MLT reformulations.
△ Less
Submitted 7 March, 2018;
originally announced March 2018.
-
The magnitude of viscous dissipation in strongly stratified two-dimensional convection
Authors:
Laura K. Currie,
Matthew K. Browning
Abstract:
Convection in astrophysical systems must be maintained against dissipation. Although the effects of dissipation are often assumed to be negligible, theory suggests that in strongly stratified convecting fluids, the dissipative heating rate can exceed the luminosity carried by convection. Here we explore this possibility using a series of numerical simulations. We consider two-dimensional numerical…
▽ More
Convection in astrophysical systems must be maintained against dissipation. Although the effects of dissipation are often assumed to be negligible, theory suggests that in strongly stratified convecting fluids, the dissipative heating rate can exceed the luminosity carried by convection. Here we explore this possibility using a series of numerical simulations. We consider two-dimensional numerical models of hydrodynamic convection in a Cartesian layer under the anelastic approximation and demonstrate that the dissipative heating rate can indeed exceed the imposed luminosity. We establish a theoretical expression for the ratio of the dissipative heating rate to the luminosity emerging at the upper boundary, in terms only of the depth of the layer and the thermal scale height. In particular, we show that this ratio is independent of the diffusivities and confirm this with a series of numerical simulations. Our results suggest that dissipative heating may significantly alter the internal dynamics of stars and planets.
△ Less
Submitted 27 July, 2017;
originally announced July 2017.
-
Results from a set of three-dimensional numerical experiments of a hot Jupiter atmosphere
Authors:
N. J. Mayne,
F. Debras,
I. Baraffe,
John Thuburn,
David S. Amundsen,
David M. Acreman,
Chris Smith,
Matthew K. Browning,
James Manners,
Nigel Wood
Abstract:
We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom bou…
▽ More
We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow.
△ Less
Submitted 3 April, 2017;
originally announced April 2017.
-
The Suppression and Promotion of Magnetic Flux Emergence in Fully Convective Stars
Authors:
Maria A. Weber,
Matthew K. Browning,
Suzannah Boardman,
Joshua Clarke,
Samuel Pugsley,
Edward Townsend
Abstract:
Evidence of surface magnetism is now observed on an increasing number of cool stars. The detailed manner by which dynamo-generated magnetic fields giving rise to starspots traverse the convection zone still remains unclear. Some insight into this flux emergence mechanism has been gained by assuming bundles of magnetic field can be represented by idealized thin flux tubes (TFTs). Weber & Browning (…
▽ More
Evidence of surface magnetism is now observed on an increasing number of cool stars. The detailed manner by which dynamo-generated magnetic fields giving rise to starspots traverse the convection zone still remains unclear. Some insight into this flux emergence mechanism has been gained by assuming bundles of magnetic field can be represented by idealized thin flux tubes (TFTs). Weber & Browning (2016) have recently investigated how individual flux tubes might evolve in a 0.3 solar-mass M dwarf by effectively embedding TFTs in time-dependent flows representative of a fully convective star. We expand upon this work by initiating flux tubes at various depths in the upper 50-75% of the star in order to sample the differing convective flow pattern and differential rotation across this region. Specifically, we comment on the role of differential rotation and time-varying flows in both the suppression and promotion of the magnetic flux emergence process.
△ Less
Submitted 15 March, 2017;
originally announced March 2017.
-
Modeling the Rise of Fibril Magnetic Fields in Fully Convective Stars
Authors:
Maria A. Weber,
Matthew K. Browning
Abstract:
Many fully convective stars exhibit a wide variety of surface magnetism, including starspots and chromospheric activity. The manner by which bundles of magnetic field traverse portions of the convection zone to emerge at the stellar surface is not especially well understood. In the Solar context, some insight into this process has been gleaned by regarding the magnetism as consisting partly of ide…
▽ More
Many fully convective stars exhibit a wide variety of surface magnetism, including starspots and chromospheric activity. The manner by which bundles of magnetic field traverse portions of the convection zone to emerge at the stellar surface is not especially well understood. In the Solar context, some insight into this process has been gleaned by regarding the magnetism as consisting partly of idealized thin flux tubes (TFT). Here, we present the results of a large set of TFT simulations in a rotating spherical domain of convective flows representative of a 0.3 solar-mass, main-sequence star. This is the first study to investigate how individual flux tubes in such a star might rise under the combined influence of buoyancy, convection, and differential rotation. A time-dependent hydrodynamic convective flow field, taken from separate 3D simulations calculated with the anelastic equations, impacts the flux tube as it rises. Convective motions modulate the shape of the initially buoyant flux ring, promoting localized rising loops. Flux tubes in fully convective stars have a tendency to rise nearly parallel to the rotation axis. However, the presence of strong differential rotation allows some initially low latitude flux tubes of moderate strength to develop rising loops that emerge in the near-equatorial region. Magnetic pumping suppresses the global rise of the flux tube most efficiently in the deeper interior and at lower latitudes. The results of these simulations aim to provide a link between dynamo-generated magnetic fields, fluid motions, and observations of starspots for fully convective stars.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
Theoretical limits on magnetic field strengths in low-mass stars
Authors:
Matthew K. Browning,
Maria A. Weber,
Gilles Chabrier,
Angela P. Massey
Abstract:
Observations have suggested that some low-mass stars have larger radii than predicted by 1-D structure models. Some theoretical models have invoked very strong interior magnetic fields (of order 1 MG or more) as a possible cause of such large radii. Whether fields of that strength could in principle by generated by dynamo action in these objects is unclear, and we do not address the matter directl…
▽ More
Observations have suggested that some low-mass stars have larger radii than predicted by 1-D structure models. Some theoretical models have invoked very strong interior magnetic fields (of order 1 MG or more) as a possible cause of such large radii. Whether fields of that strength could in principle by generated by dynamo action in these objects is unclear, and we do not address the matter directly. Instead, we examine whether such fields could remain in the interior of a low mass object for a significant time, and whether they would have any other obvious signatures. First, we estimate timescales for the loss of strong fields by magnetic buoyancy instabilities. We consider a range of field strengths and simple morphologies, including both idealized flux tubes and smooth layers of field. We confirm some of our analytical estimates using thin flux tube magnetohydrodynamic (MHD) simulations of the rise of buoyant fields in a fully-convective M-dwarf. Separately, we consider the Ohmic dissipation of such fields. We find that dissipation provides a complementary constraint to buoyancy: while small-scale, fibril fields might be regenerated faster than they rise, the dissipative heating associated with such fields would in some cases greatly exceed the luminosity of the star. We show how these constraints combine to yield limits on the internal field strength and morphology in low-mass stars. In particular, we find that for stars of 0.3 solar masses, no fields in flux tubes stronger than about 800 kG are simultaneously consistent with both constraints.
△ Less
Submitted 17 December, 2015;
originally announced December 2015.
-
Helicity inversion in spherical convection as a means for equatorward dynamo wave propagation
Authors:
Lúcia D. V. Duarte,
Johannes Wicht,
Matthew K. Browning,
Thomas Gastine
Abstract:
We discuss here a purely hydrodynamical mechanism to invert the sign of the kinetic helicity, which plays a key role in determining the direction of propagation of cyclical magnetism in most models of dynamo action by rotating convection. Such propagation provides a prominent, and puzzling constraint on dynamo models. In the Sun, active regions emerge first at mid-latitudes, then appear nearer the…
▽ More
We discuss here a purely hydrodynamical mechanism to invert the sign of the kinetic helicity, which plays a key role in determining the direction of propagation of cyclical magnetism in most models of dynamo action by rotating convection. Such propagation provides a prominent, and puzzling constraint on dynamo models. In the Sun, active regions emerge first at mid-latitudes, then appear nearer the equator over the course of a cycle, but most previous global-scale dynamo simulations have exhibited poleward propagation (if they were cyclical at all). Here, we highlight some simulations in which the direction of propagation of dynamo waves is altered primarily by an inversion of the kinetic helicity throughout much of the interior, rather than by changes in the differential rotation. This tends to occur in cases with a low Prandtl number and internal heating, in regions where the local density gradient is relatively small. We analyse how this inversion arises, and contrast it to the case of convection that is either highly columnar (i.e., rapidly rotating) or locally very stratified; in both of those situations, the typical profile of kinetic helicity (negative throughout most of the northern hemisphere) instead prevails.
△ Less
Submitted 3 December, 2015; v1 submitted 18 November, 2015;
originally announced November 2015.
-
The UM, a fully--compressible, non--hydrostatic, deep atmosphere GCM, applied to hot Jupiters
Authors:
N. J. Mayne,
I. Baraffe,
David M. Acreman,
Chris Smith,
Matthew K. Browning,
David Skålid Amundsen,
Nigel Wood,
John Thuburn,
David R. Jackson
Abstract:
We are adapting the Global Circulation Model (GCM) of the UK Met Office, the so--called Unified Model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric model…
▽ More
We are adapting the Global Circulation Model (GCM) of the UK Met Office, the so--called Unified Model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment). Within the same numerical scheme ENDGame supports solution to the dynamical equations under varying degrees of simplification. We present results from a simple, shallow (in atmospheric domain) hot Jupiter model (SHJ), and a more realistic (with a deeper atmosphere) HD 209458b test case. For both test cases we find that the large--scale, time--averaged (over the 1200 days prescribed test period), dynamical state of the atmosphere is relatively insensitive to the level of simplification of the dynamical equations. However, problems exist when attempting to reproduce the results for these test cases derived from other models. For the SHJ case the lower (and upper) boundary intersects the dominant dynamical features of the atmosphere meaning the results are heavily dependent on the boundary conditions. For the HD 209458b test case, when using the more complete dynamical models, the atmosphere is still clearly evolving after 1200 days, and in a transient state. Solving the complete (deep atmosphere and non--hydrostatic) dynamical equations allows exchange between the vertical and horizontal momentum of the atmosphere, via Coriolis and metric terms. Subsequently, interaction between the upper atmosphere and the deeper more slowly evolving (radiatively inactive) atmosphere significantly alters the results, and acts over timescales longer than 1200 days.
△ Less
Submitted 8 November, 2013;
originally announced November 2013.
-
Magnetic Cycles in a Convective Dynamo Simulation of a Young Solar-type Star
Authors:
Benjamin P. Brown,
Mark S. Miesch,
Matthew K. Browning,
Allen Sacha Brun,
Juri Toomre
Abstract:
Young solar-type stars rotate rapidly and many are magnetically active; some undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the convective envelope of a solar-type star rotating at 5 times the current solar rotation rate…
▽ More
Young solar-type stars rotate rapidly and many are magnetically active; some undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the convective envelope of a solar-type star rotating at 5 times the current solar rotation rate. Striking global-scale magnetic wreaths appear in the midst of the turbulent convection zone and show rich time-dependence. The dynamo exhibits cyclic activity and undergoes quasi-periodic polarity reversals where both the global-scale poloidal and toroidal fields change in sense on a roughly 1500 day time scale. These magnetic activity patterns emerge spontaneously from the turbulent flow and are more organized temporally and spatially than those realized in our previous simulations of the solar dynamo. We assess in detail the competing processes of magnetic field creation and destruction within our simulations that contribute to the global-scale reversals. We find that the mean toroidal fields are built primarily through an $Ω$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $α$-effect. During a reversal the magnetic wreaths propagate towards the polar regions, and this appears to arise from a poleward propagating dynamo wave. The primary response in the convective flows involves the axisymmetric differential rotation which shows variations associated with the poleward propagating magnetic wreaths. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent poleward propagating magnetic fields deep below the solar surface. [abridged]
△ Less
Submitted 9 February, 2011;
originally announced February 2011.
-
Global-scale Magnetism (and Cycles) in Dynamo Simulations of Stellar Convection Zones
Authors:
Benjamin P Brown,
Matthew K Browning,
Allan Sacha Brun,
Mark S Miesch,
Juri Toomre
Abstract:
Young solar-type stars rotate rapidly and are very magnetically active. The magnetic fields at their surfaces likely originate in their convective envelopes where convection and rotation can drive strong dynamo action. Here we explore simulations of global-scale stellar convection in rapidly rotating suns using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these…
▽ More
Young solar-type stars rotate rapidly and are very magnetically active. The magnetic fields at their surfaces likely originate in their convective envelopes where convection and rotation can drive strong dynamo action. Here we explore simulations of global-scale stellar convection in rapidly rotating suns using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. We explore one case rotates five times faster than the Sun in detail. This dynamo simulation, called case D5, has repeated quasi-cyclic reversals of global-scale polarity. We compare this case D5 to the broader family of simulations we have been able to explore and discuss how future simulations and observations can advance our understanding of stellar dynamos and magnetism.
△ Less
Submitted 30 December, 2010;
originally announced January 2011.
-
Persistent Magnetic Wreaths in a Rapidly Rotating Sun
Authors:
Benjamin P. Brown,
Matthew K. Browning,
Allen Sacha Brun,
Mark S. Miesch,
Juri Toomre
Abstract:
When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magneti…
▽ More
When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for a system rotating at three times the current solar rotation rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in this system. Striking wreaths of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a surprise, for it has been widely believed that such magnetic structures should be disrupted by magnetic buoyancy or turbulent pumping. Thus, many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. We examine how these persistent magnetic wreaths are maintained by dynamo processes and explore whether a classical mean-field $α$-effect explains the regeneration of poloidal field.
△ Less
Submitted 12 November, 2010;
originally announced November 2010.
-
Global-scale wreath-building dynamos in stellar convection zones
Authors:
Benjamin P Brown,
Matthew K Browning,
Allan S Brun,
Mark S Miesch,
Juri Toomre
Abstract:
When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity…
▽ More
When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial or temporal organization. Wreath-building dynamos are present to some degree at all rotation rates, but are most evident in the more rapidly rotating simulations.
△ Less
Submitted 1 November, 2010;
originally announced November 2010.
-
Global magnetic cycles in rapidly rotating younger suns
Authors:
Nicholas J. Nelson,
Benjamin P. Brown,
Matthew K. Browning,
Allan Sacha Brun,
Mark S. Miesch,
Juri Toomre
Abstract:
Observations of sun-like stars rotating faster than our current sun tend to exhibit increased magnetic activity as well as magnetic cycles spanning multiple years. Using global simulations in spherical shells to study the coupling of large-scale convection, rotation, and magnetism in a younger sun, we have probed effects of rotation on stellar dynamos and the nature of magnetic cycles. Major 3-D M…
▽ More
Observations of sun-like stars rotating faster than our current sun tend to exhibit increased magnetic activity as well as magnetic cycles spanning multiple years. Using global simulations in spherical shells to study the coupling of large-scale convection, rotation, and magnetism in a younger sun, we have probed effects of rotation on stellar dynamos and the nature of magnetic cycles. Major 3-D MHD simulations carried out at three times the current solar rotation rate reveal hydromagnetic dynamo action that yields wreaths of strong toroidal magnetic field at low latitudes, often with opposite polarity in the two hemispheres. Our recent simulations have explored behavior in systems with considerably lower diffusivities, achieved with sub-grid scale models including a dynamic Smagorinsky treatment of unresolved turbulence. The lower diffusion promotes the generation of magnetic wreaths that undergo prominent temporal variations in field strength, exhibiting global magnetic cycles that involve polarity reversals. In our least diffusive simulation, we find that magnetic buoyancy coupled with advection by convective giant cells can lead to the rise of coherent loops of magnetic field toward the top of the simulated domain.
△ Less
Submitted 28 October, 2010;
originally announced October 2010.
-
Magnetic Cycles and Meridional Circulation in Global Models of Solar Convection
Authors:
Mark S. Miesch,
Benjamin P. Brown,
Matthew K. Browning,
Allan Sacha Brun,
Juri Toomre
Abstract:
We review recent insights into the dynamics of the solar convection zone obtained from global numerical simulations, focusing on two recent developments in particular. The first is quasi-cyclic magnetic activity in a long-duration dynamo simulation. Although mean fields comprise only a few percent of the total magnetic energy they exhibit remarkable order, with multiple polarity reversals and syst…
▽ More
We review recent insights into the dynamics of the solar convection zone obtained from global numerical simulations, focusing on two recent developments in particular. The first is quasi-cyclic magnetic activity in a long-duration dynamo simulation. Although mean fields comprise only a few percent of the total magnetic energy they exhibit remarkable order, with multiple polarity reversals and systematic variability on time scales of 6-15 years. The second development concerns the maintenance of the meridional circulation. Recent high-resolution simulations have captured the subtle nonlinear dynamical balances with more fidelity than previous, more laminar models, yielding more coherent circulation patterns. These patterns are dominated by a single cell in each hemisphere, with poleward and equatorward flow in the upper and lower convection zone respectively. We briefly address the implications of and future of these modeling efforts.
△ Less
Submitted 2 October, 2010; v1 submitted 30 September, 2010;
originally announced September 2010.
-
Wreathes of Magnetism in Rapidly Rotating Suns
Authors:
Benjamin P. Brown,
Matthew K. Browning,
Mark S. Miesch,
Allan Sacha Brun,
Juri Toomre
Abstract:
When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magne…
▽ More
When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for two systems, rotating at three and five times the current solar rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in these systems. Striking wreathes of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a great surprise, for many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. Some dynamos achieved in these rapidly rotating states build persistent global-scale fields which maintain amplitude and polarity for thousands of days. In the case at five times the solar rate, the dynamo can undergo cycles of activity, with fields varying in strength and even changing polarity. As the magnetic fields wax and wane in strength, the primary response in the convective flows involves the axisymmetric differential rotation, which begins to vary on similar time scales. Bands of relatively fast and slow fluid propagate toward the poles on time scales of roughly 500 days. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent a response to poleward propagating magnetic field deep below the solar surface.
△ Less
Submitted 12 June, 2009;
originally announced June 2009.
-
Three-Dimensional Simulations of Solar and Stellar Dynamos: The Influence of a Tachocline
Authors:
M. S. Miesch,
M. K. Browning,
A. S. Brun,
J. Toomre,
B. P. Brown
Abstract:
We review recent advances in modeling global-scale convection and dynamo processes with the Anelastic Spherical Harmonic (ASH) code. In particular, we have recently achieved the first global-scale solar convection simulations that exhibit turbulent pumping of magnetic flux into a simulated tachocline and the subsequent organization and amplification of toroidal field structures by rotational she…
▽ More
We review recent advances in modeling global-scale convection and dynamo processes with the Anelastic Spherical Harmonic (ASH) code. In particular, we have recently achieved the first global-scale solar convection simulations that exhibit turbulent pumping of magnetic flux into a simulated tachocline and the subsequent organization and amplification of toroidal field structures by rotational shear. The presence of a tachocline not only promotes the generation of mean toroidal flux, but it also enhances and stabilizes the mean poloidal field throughout the convection zone, promoting dipolar structure with less frequent polarity reversals. The magnetic field generated by a convective dynamo with a tachocline and overshoot region is also more helical overall, with a sign reversal in the northern and southern hemispheres. Toroidal tachocline fields exhibit little indication of magnetic buoyancy instabilities but may be undergoing magneto-shear instabilities.
△ Less
Submitted 18 November, 2008;
originally announced November 2008.
-
Rapidly Rotating Suns and Active Nests of Convection
Authors:
Benjamin P. Brown,
Matthew K. Browning,
Allan Sacha Brun,
Mark S. Miesch,
Juri Toomre
Abstract:
In the solar convection zone, rotation couples with intensely turbulent convection to drive a strong differential rotation and achieve complex magnetic dynamo action. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the global-scale patterns of convection in s…
▽ More
In the solar convection zone, rotation couples with intensely turbulent convection to drive a strong differential rotation and achieve complex magnetic dynamo action. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the global-scale patterns of convection in such stars and the flows of differential rotation and meridional circulation which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of localized convection emerges. Convection near the equator in these systems is dominated by one or two nests in longitude of locally enhanced convection, with quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation Delta Omega grows with more rapid rotation while the relative shear Delta Omega/Omega_0 decreases. In contrast, at more rapid rotation the meridional circulations decrease in energy and peak velocities and break into multiple cells of circulation in both radius and latitude.
△ Less
Submitted 12 August, 2008;
originally announced August 2008.
-
Strong Dynamo Action in Rapidly Rotating Suns
Authors:
Benjamin P. Brown,
Matthew K. Browning,
Allan Sacha Brun,
Mark S. Miesch,
Nicholas J. Nelson,
Juri Toomre
Abstract:
Stellar dynamos are driven by complex couplings between rotation and turbulent convection, which drive global-scale flows and build and rebuild stellar magnetic fields. When stars like our sun are young, they rotate much more rapidly than the current solar rate. Observations generally indicate that more rapid rotation is correlated with stronger magnetic activity and perhaps more effective dynam…
▽ More
Stellar dynamos are driven by complex couplings between rotation and turbulent convection, which drive global-scale flows and build and rebuild stellar magnetic fields. When stars like our sun are young, they rotate much more rapidly than the current solar rate. Observations generally indicate that more rapid rotation is correlated with stronger magnetic activity and perhaps more effective dynamo action. Here we examine the effects of more rapid rotation on dynamo action in a star like our sun. We find that vigorous dynamo action is realized, with magnetic field generated throughout the bulk of the convection zone. These simulations do not possess a penetrative tachocline of shear where global-scale fields are thought to be organized in our sun, but despite this we find strikingly ordered fields, much like sea-snakes of toroidal field, which are organized on global scales. We believe this to be a novel finding.
△ Less
Submitted 10 January, 2008;
originally announced January 2008.
-
Rapid Rotation, Active Nests of Convection and Global-scale Flows in Solar-like Stars
Authors:
Benjamin P. Brown,
Matthew K. Browning,
Allan Sacha Brun,
Mark S. Miesch,
Juri Toomre
Abstract:
In the solar convection zone, rotation couples with intensely turbulent convection to build global-scale flows of differential rotation and meridional circulation. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the patterns of convection in such stars and th…
▽ More
In the solar convection zone, rotation couples with intensely turbulent convection to build global-scale flows of differential rotation and meridional circulation. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the patterns of convection in such stars and the global-scale flows which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of spatially localized convection emerges. Convection near the equator in these systems is dominated by one or two patches of locally enhanced convection, with nearly quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation, as measured by latitudinal angular velocity contrast, Delta_Omega, increases with more rapid rotation while the relative shear, Delta_Omega/Omega, decreases. In contrast, at more rapid rotation the meridional circulations decrease in both energy and peak velocities and break into multiple cells of circulation in both radius and latitude.
△ Less
Submitted 10 January, 2008;
originally announced January 2008.
-
Simulations of core convection in rotating A-type stars: Magnetic dynamo action
Authors:
Allan Sacha Brun,
Matthew K. Browning,
Juri Toomre
Abstract:
Core convection and dynamo activity deep within rotating A-type stars of 2 solar masses are studied with 3--D nonlinear simulations. Our modeling considers the inner 30% by radius of such stars, thus capturing within a spherical domain the convective core and a modest portion of the surrounding radiative envelope. The MHD equations are solved using the ASH code to examine turbulent flows and mag…
▽ More
Core convection and dynamo activity deep within rotating A-type stars of 2 solar masses are studied with 3--D nonlinear simulations. Our modeling considers the inner 30% by radius of such stars, thus capturing within a spherical domain the convective core and a modest portion of the surrounding radiative envelope. The MHD equations are solved using the ASH code to examine turbulent flows and magnetic fields, both of which exhibit intricate time dependence. By introducing small seed magnetic fields into our progenitor hydrodynamic models rotating at one and four times the solar rate, we assess here how the vigorous convection can amplify those fields and sustain them against ohmic decay. Dynamo action is indeed realized, ultimately yielding magnetic fields that are in energy equipartion with the flow. Such magnetism reduces the differential rotation obtained in the progenitors, partly by Maxwell stresses that transport angular momentum poleward and oppose the Reynolds stresses in the latitudinal balance. In contrast, in the radial direction we find that the Maxwell and Reynolds stresses may act together to transport angular momentum. The central columns of slow rotation established in the progenitors are weakened, with the differential rotation waxing and waning in strength as the simulations evolve. We assess the morphology of the flows and magnetic fields, their complex temporal variations, and the manner in which dynamo action is sustained. Differential rotation and helical convection are both found to play roles in giving rise to the magnetic fields. The magnetism is dominated by strong fluctuating fields throughout the core, with the axisymmetric (mean) fields there relatively weak.
△ Less
Submitted 3 October, 2006;
originally announced October 2006.
-
Inferring physical conditions in interstellar clouds of H_2
Authors:
Matthew K. Browning,
Jason Tumlinson,
J. Michael Shull
Abstract:
We have developed a code that models the formation, destruction, radiative transfer, and vibrational/rotational excitation of H_2 in a detailed fashion. We discuss how such codes, together with FUSE observations of H_2 in diffuse and translucent lines of sight, may be used to infer various physical parameters. We illustrate the effects of changes in the major physical parameters (UV radiation fi…
▽ More
We have developed a code that models the formation, destruction, radiative transfer, and vibrational/rotational excitation of H_2 in a detailed fashion. We discuss how such codes, together with FUSE observations of H_2 in diffuse and translucent lines of sight, may be used to infer various physical parameters. We illustrate the effects of changes in the major physical parameters (UV radiation field, gas density, metallicity), and we point out the extent to which changes in one parameter may be mirrored by changes in another. We provide an analytic formula for the molecular fraction, f_H2, as a function of cloud column density, radiation field, and grain formation rate of H_2. Some diffuse and translucent lines of sight may be concatenations of multiple distinct clouds viewed together. Such situations can give rise to observables that agree with the data, complicating the problem of uniquely identifying one set of physical parameters with a line of sight. Finally, we illustrate the application of our code to an ensemble of data, such as the FUSE survey of H_2 in the Large and Small Magellanic Clouds (LMC/SMC), in order to constrain the elevated UV radiation field intensity and reduced grain formation rate of H_2 in those low- metallicity environments.
△ Less
Submitted 19 September, 2002;
originally announced September 2002.
-
A FUSE Survey of Interstellar Molecular Hydrogen in the Small and Large Magellanic Clouds
Authors:
Jason Tumlinson,
J. Michael Shull,
Brian L. Rachford,
Matthew K. Browning,
Theodore P. Snow,
Alex W. Fullerton,
Edward B. Jenkins,
Blair D. Savage,
Paul A. Crowther,
H. Warren Moos,
Kenneth R. Sembach,
George Sonneborn,
Donald G. York
Abstract:
We describe a moderate-resolution FUSE survey of H2 along 70 sight lines to the Small and Large Magellanic Clouds, using hot stars as background sources. FUSE spectra of 67% of observed Magellanic Cloud sources (52% of LMC and 92% of SMC) exhibit absorption lines from the H2 Lyman and Werner bands between 912 and 1120 A. Our survey is sensitive to N(H2) >= 10^14 cm^-2; the highest column densiti…
▽ More
We describe a moderate-resolution FUSE survey of H2 along 70 sight lines to the Small and Large Magellanic Clouds, using hot stars as background sources. FUSE spectra of 67% of observed Magellanic Cloud sources (52% of LMC and 92% of SMC) exhibit absorption lines from the H2 Lyman and Werner bands between 912 and 1120 A. Our survey is sensitive to N(H2) >= 10^14 cm^-2; the highest column densities are log N(H2) = 19.9 in the LMC and 20.6 in the SMC. We find reduced H2 abundances in the Magellanic Clouds relative to the Milky Way, with average molecular fractions <f_H2> = 0.010 (+0.005, -0.002) for the SMC and <f_H2> = 0.012 (+0.006, -0.003) for the LMC, compared with <f_H2> = 0.095 for the Galactic disk over a similar range of reddening. The dominant uncertainty in this measurement results from the systematic differences between 21 cm radio emission and Lya in pencil-beam sight lines as measures of N(HI). These results imply that the diffuse H2 masses of the LMC and SMC are 8 x 10^6 Msun and 2 x 10^6 Msun, respectively, 2% and 0.5% of the H I masses derived from 21 cm emission measurements. The LMC and SMC abundance patterns can be reproduced in ensembles of model clouds with a reduced H2 formation rate coefficient, R ~ 3 x 10^-18 cm^3 s^-1, and incident radiation fields ranging from 10 - 100 times the Galactic mean value. We find that these high-radiation, low-formation-rate models can also explain the enhanced N(4)/N(2) and N(5)/N(3) rotational excitation ratios in the Clouds. We use H2 column densities in low rotational states (J = 0 and 1) to derive a mean kinetic and/or rotational temperature <T_01> = 82 +/- 21 K for clouds with N(H2) >= 10^16 cm^-2, similar to Galactic gas. We discuss the implications of this work for theories of star formation in low-metallicity environments. [Abstract abridged]
△ Less
Submitted 10 October, 2001;
originally announced October 2001.