Showing 1–50 of 50 results for author: Brown, B P

The solar dynamo begins near the surface

Authors: Geoffrey M Vasil, Daniel Lecoanet, Kyle Augustson, Keaton J Burns, Jeffrey S Oishi, Benjamin P Brown, Nicholas Brummell, Keith Julien

Abstract: The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pin… ▽ More The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the Sun's outer 5-10%, the "Near-Surface Shear Layer". Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the Magneto-Rotational Instability prominent in accretion-disk theory and observed in laboratory experiments. Together, these two facts prompt the general question: Is it possible that the solar dynamo is a near-surface instability? Here, we report strong affirmative evidence in stark contrast to traditional paradigms focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes. State-of-the-art numerical simulations corroborate these estimates and, strikingly, reproduce hemispherical magnetic current helicity laws. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, impacting Earth's electromagnetic infrastructure. △ Less

Submitted 13 April, 2024; v1 submitted 11 April, 2024; originally announced April 2024.

Comments: In press, Nature

Internally heated and fully compressible convection: flow morphology and scaling laws

Authors: Whitney T. Powers, Evan H. Anders, Benjamin P. Brown

Abstract: In stars and planets natural processes heat convective flows in the bulk of a convective region rather than at hard boundaries. By characterizing how convective dynamics are determined by the strength of an internal heating source we can gain insight into the processes driving astrophysical convection. Internally heated convection has been studied extensively in incompressible fluids, but the effe… ▽ More In stars and planets natural processes heat convective flows in the bulk of a convective region rather than at hard boundaries. By characterizing how convective dynamics are determined by the strength of an internal heating source we can gain insight into the processes driving astrophysical convection. Internally heated convection has been studied extensively in incompressible fluids, but the effects of stratification and compressibility have not been examined in detail. In this work, we study fully compressible convection driven by a spatially uniform heating source in 2D and 3D Cartesian, hydrodynamic simulations. We use a fixed temperature upper boundary condition which results in a system that is internally heated in the bulk and cooled at the top. We find that the flow speed, as measured by the Mach number, and turbulence, as measured by the Reynolds number, can be independently controlled by separately varying the characteristic temperature gradient from internal heating and the diffusivities. 2D simulations at a fixed Mach number (flow speed) demonstrate consistent power at low wavenumber as diffusivities are decreased. We observe convection where the velocity distribution is skewed towards cold, fast downflows, and that the flow speed is related to the length scale and entropy gradient of the upper boundary where the downflows are driven. We additionally find a heat transport scaling law which is consistent with prior incompressible work. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 22 pages, 12 figures, submitted to Phys. Rev. Fluids

Journal ref: Phys. Rev. Fluids 9 (2024), 043501

The photometric variability of massive stars due to gravity waves excited by core convection

Authors: Evan H. Anders, Daniel Lecoanet, Matteo Cantiello, Keaton J. Burns, Benjamin A. Hyatt, Emma Kaufman, Richard H. D. Townsend, Benjamin P. Brown, Geoffrey M. Vasil, Jeffrey S. Oishi, Adam S. Jermyn

Abstract: Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometri… ▽ More Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present the first 3D simulations of massive star convection extending from the star's center to near its surface, with realistic stellar luminosities. Using these simulations, we make the first prediction of photometric variability due to convectively-driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, so the red noise must be generated by an alternative process. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: As accepted for publication in Nature Astronomy except for final editorial revisions. Supplemental materials available online at https://doi.org/10.5281/zenodo.7764997 . We have also sonified our results to make them more accessible, see https://github.com/evanhanders/gmode_variability_paper/blob/main/sound/gmode_sonification.pdf

Journal ref: Nature Astronomy, 2023, 10.1038/s41550-023-02040-7

The Puzzling Structure of Solar Convection: Window into the Dynamo

Authors: Nicholas A. Featherstone, Kyle C. Augustson, Jonathan M. Aurnou, Catherine Blume, Benjamin P. Brown, Nicholas Brummell, Keaton J. Burns, Michael A. Calkins, Maria Camisassa, Mausumi Dikpati, Yuhong Fan, J. R. Fuentes, Gustavo Guerrero, Bradley W. Hindman, Keith Julien, Irina N. Kitiashvili, Lydia Korre, Daniel Lecoanet, Bhishek Manek, Loren Matilsky, Mark Miesch, Nicholas J. Nelson, Jeffrey S. Oishi, Whitney T. Powers, Matthias Rempel , et al. (3 additional authors not shown)

Abstract: The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective… ▽ More The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics 2024-2033. 7 pages; 4 figures

The Stability of Prendergast Magnetic Fields

Authors: Emma Kaufman, Daniel Lecoanet, Evan H. Anders, Benjamin P. Brown, Geoffrey M. Vasil, Jeffrey S. Oishi, Keaton J. Burns

Abstract: Convection in massive main sequence stars generates large scale magnetic fields in their cores which persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components which are expected to stabilize each other. Previous analytic and numerical calculations did not find any eviden… ▽ More Convection in massive main sequence stars generates large scale magnetic fields in their cores which persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components which are expected to stabilize each other. Previous analytic and numerical calculations did not find any evidence for instability of the Prendergast field over short timescales. In this paper, we present numerical simulations which show a long timescale, linear instability of this magnetic field. We find the instability to be robust to changes in boundary conditions and it is not stabilized by strong stable stratification. The instability is a resistive instability, and the growth rate has a power-law dependence on the resistivity, in which the growth rate decreases as the resistivity decreases. We estimate the growth rate of the instability in stars by extrapolating this power-law to stellar values of the resistivity. The instability is sufficiently rapid to destabilize the magnetic field on timescales shorter than the stellar evolution timescale, indicating that the Prendergast field is not a good model to use in studies of magnetic fields in stars. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 9 pages, 8 figures

Moosinesq Convection in the Cores of Moosive Stars

Authors: Evan H. Anders, Evan B. Bauer, Adam S. Jermyn, Samuel J. Van Kooten, Benjamin P. Brown, Eric W. Hester, Mindy Wilkinson, Jared A. Goldberg, Tania Varesano, Daniel Lecoanet

Abstract: Stars with masses $\gtrsim 4 \times 10^{27}M_{\rm{moose}} \approx 1.1 M_\odot$ have core convection zones during their time on the main sequence. In these moosive stars, convection introduces many uncertainties in stellar modeling. In this Letter, we build upon the Boussinesq approximation to present the first-ever simulations of Moosinesq convection, which captures the complex geometric structure… ▽ More Stars with masses $\gtrsim 4 \times 10^{27}M_{\rm{moose}} \approx 1.1 M_\odot$ have core convection zones during their time on the main sequence. In these moosive stars, convection introduces many uncertainties in stellar modeling. In this Letter, we build upon the Boussinesq approximation to present the first-ever simulations of Moosinesq convection, which captures the complex geometric structure of the convection zones of these stars. These flows are bounded in a manner informed by the majestic terrestrial Alces alces (moose) and could have important consequences for the evolution of these stars. We find that Moosinesq convection results in very interesting flow morphologies and rapid heat transfer, and posit this as a mechanism of biomechanical thermoregulation. △ Less

Submitted 30 March, 2022; originally announced April 2022.

Comments: 7 pages, 3 figures, 1 moose, no logic in this place

Stellar convective penetration: parameterized theory and dynamical simulations

Authors: Evan H. Anders, Adam S. Jermyn, Daniel Lecoanet, Benjamin P. Brown

Abstract: Most stars host convection zones in which heat is transported directly by fluid motion, but the behavior of convective boundaries is not well understood. Here we present 3D numerical simulations which exhibit penetration zones: regions where the entire luminosity \emph{could} be carried by radiation, but where the temperature gradient is approximately adiabatic and convection is present. To parame… ▽ More Most stars host convection zones in which heat is transported directly by fluid motion, but the behavior of convective boundaries is not well understood. Here we present 3D numerical simulations which exhibit penetration zones: regions where the entire luminosity \emph{could} be carried by radiation, but where the temperature gradient is approximately adiabatic and convection is present. To parameterize this effect, we define the "penetration parameter" $\mathcal{P}$ which compares how far the radiative gradient deviates from the adiabatic gradient on either side of the Schwarzschild convective boundary. Following Roxburgh (1989) and Zahn (1991), we construct an energy-based theoretical model in which $\mathcal{P}$ controls the extent of penetration. We test this theory using 3D numerical simulations which employ a simplified Boussinesq model of stellar convection. The convection is driven by internal heating and we use a height-dependent radiative conductivity; this allows us to separately specify $\mathcal{P}$ and the stiffness $\mathcal{S}$ of the radiative-convective boundary. We find significant convective penetration in all simulations. Our simple theory describes the simulations well. Penetration zones can take thousands of overturn times to develop, so long simulations or accelerated evolutionary techniques are required. In stars, we expect $\mathcal{P} \approx 1$ and in this regime our results suggest that convection zones may extend beyond the Schwarzschild boundary by up to $\sim$20-30% of a mixing length. We present a MESA stellar model of the Sun which employs our parameterization of convective penetration as a proof of concept. We discuss prospects for extending these results to more realistic stellar contexts. △ Less

Submitted 11 March, 2022; v1 submitted 21 October, 2021; originally announced October 2021.

Comments: Appears as accepted for publication in ApJ. Supplemental materials available at https://doi.org/10.5281/zenodo.5131118

Journal ref: ApJ 926 169 (2022)

Single-hemisphere dynamos in M-dwarf stars

Authors: Benjamin P. Brown, Jeffrey S. Oishi, Geoffrey M. Vasil, Daniel Lecoanet, Keaton J. Burns

Abstract: M-dwarf stars below a certain mass are convective from their cores to their photospheres. These fully convective objects are extremely numerous, very magnetically active, and the likely hosts of many exoplanets. Here we study, for the first time, dynamo action in simulations of stratified, rotating fully convective M-dwarf stars. Importantly, we use new techniques to capture the correct full ball… ▽ More M-dwarf stars below a certain mass are convective from their cores to their photospheres. These fully convective objects are extremely numerous, very magnetically active, and the likely hosts of many exoplanets. Here we study, for the first time, dynamo action in simulations of stratified, rotating fully convective M-dwarf stars. Importantly, we use new techniques to capture the correct full ball geometry down to the center of the star. We find surprising dynamo states in these systems, with the global-scale mean fields confined strongly to a single hemisphere, in contrast to prior stellar dynamo solutions. These hemispheric-dynamo stars are likely to have profoundly different interactions with their surroundings, with important implications for exoplanet habitability and stellar spindown. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Comments: 7 pages, 3 figures; submitted to ApJL

Convective dynamics with mixed temperature boundary conditions: why thermal relaxation matters and how to accelerate it

Authors: Evan H. Anders, Geoffrey M. Vasil, Benjamin P. Brown, Lydia Korre

Abstract: Astrophysical simulations of convection frequently impose different thermal boundary conditions at the top and the bottom of the domain in an effort to more accurately model natural systems. In this work, we study Rayleigh-Benard convection (RBC) under the Boussinesq approximation. We examine simulations with mixed temperature boundary conditions in which the flux is fixed at the bottom boundary a… ▽ More Astrophysical simulations of convection frequently impose different thermal boundary conditions at the top and the bottom of the domain in an effort to more accurately model natural systems. In this work, we study Rayleigh-Benard convection (RBC) under the Boussinesq approximation. We examine simulations with mixed temperature boundary conditions in which the flux is fixed at the bottom boundary and the temperature is fixed at the top ("FT"). We aim to understand how FT boundaries change the nature of the convective solution compared to the traditional choice of thermal boundaries, in which the temperature is fixed at the top and bottom of the domain ("TT"). We demonstrate that the timescale of thermal relaxation for FT simulations is dependent upon the initial conditions. "Classic" initial conditions which employ a hydrostatically -- and thermally -- balanced linear temperature profile exhibit a long thermal relaxation. This long relaxation is not seen in FT simulations which use a TT simulation's nonlinear state as initial conditions ("TT-to-FT"). In the thermally relaxed, statistically stationary state, the mean behavior of an FT simulation corresponds to an equivalent simulation with TT boundaries, and time- and volume-averaged flow statistics like the Nusselt number and the Peclet number are indistinguishable between FT and TT simulations. FT boundaries are fundamentally asymmetric, and we examine the asymmetries that these boundaries produce in the flow. We find that the fixed-flux boundary produces more extreme temperature events than the fixed-temperature boundary. However, these near-boundary asymmetries do not measurably break the symmetry in the convective interior. We briefly explore rotating RBC to demonstrate that our findings with respect to thermal relaxation carry over to this more complex case, and to show the power of TT-to-FT initial conditions. △ Less

Submitted 12 September, 2020; v1 submitted 28 February, 2020; originally announced March 2020.

Comments: Published open access in Physical Review Fluids. The arXiv v3 appears roughly as published, but the published version is better proofread

Journal ref: Phys. Rev. Fluids 5, 083501 (2020)

The Magnetorotational Instability Prefers Three Dimensions

Authors: Jeffrey S. Oishi, Geoffrey M. Vasil, Morgan Baxter, Andrew Swan, Keaton J. Burns, Daniel Lecoanet, Benjamin P. Brown

Abstract: The magnetorotational instability (MRI) occurs when a weak magnetic field destabilises a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disk theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI unstable, and they take a wide range of profiles, including near-critical. We show that the… ▽ More The magnetorotational instability (MRI) occurs when a weak magnetic field destabilises a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disk theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI unstable, and they take a wide range of profiles, including near-critical. We show that the fastest growing modes of an ideal magnetofluid are three-dimensional provided the shear rate, $S$, is near the two-dimensional onset value, $S_c$. For a Keplerian shear, three-dimensional modes are unstable above $S\approx0.10S_c$, and dominate the two-dimensional modes until $S\approx2.05S_{c}$. These three-dimensional modes dominate for shear profiles relevant to stars and at magnetic Prandtl numbers relevant to liquid-metal laboratory experiments. Significant numbers of rapidly growing three-dimensional modes remain well past $2.05S_{c}$. These finding are significant in three ways. First, weakly nonlinear theory suggests that the MRI saturates by pushing the shear rate to its critical value. This can happen for systems, like stars and laboratory experiments, that can rearrange their angular velocity profiles. Second, the non-normal character and large transient growth of MRI modes should be important whenever three-dimensionality exists. Finally, three-dimensional growth suggests direct dynamo action driven from the linear instability. △ Less

Submitted 2 December, 2019; originally announced December 2019.

Comments: 15 pages, 5 figures. Proceedings of the Royal Society A, accepted

Entropy Rain: Dilution and Compression of Thermals in Stratified Domains

Authors: Evan H. Anders, Daniel Lecoanet, Benjamin P. Brown

Abstract: Large-scale convective flows called giant cells were once thought to transport the Sun's luminosity in the solar convection zone, but recent observations have called their existence into question. In place of large-scale flows, some authors have suggested the solar luminosity may instead be transported by small droplets of rapidly falling, low entropy fluid. This "entropy rain" could propagate as… ▽ More Large-scale convective flows called giant cells were once thought to transport the Sun's luminosity in the solar convection zone, but recent observations have called their existence into question. In place of large-scale flows, some authors have suggested the solar luminosity may instead be transported by small droplets of rapidly falling, low entropy fluid. This "entropy rain" could propagate as dense vortex rings, analogous to rising buoyant thermals in the Earth's atmosphere. In this work, we develop an analytical theory describing the evolution of dense, negatively buoyant thermals. We verify the theory with 2D cylindrical and 3D cartesian simulations of laminar, axisymmetric thermals in highly stratified atmospheres. Our results show that dense thermals fall in two categories: a stalling regime in which the droplets slow down and expand, and a falling regime in which the droplets accelerate and shrink as they propagate downwards. We estimate that solar downflows are in the falling regime and maintain their entropy perturbation against diffusion until they reach the base of the convection zone. This suggests that entropy rain may be an effective nonlocal mechanism for transporting the solar luminosity. △ Less

Submitted 16 October, 2019; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: As published in ApJ. 12 pages, including 6 figures, 2 tables, and 3 appendices

Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods

Authors: Keaton J. Burns, Geoffrey M. Vasil, Jeffrey S. Oishi, Daniel Lecoanet, Benjamin P. Brown

Abstract: Numerical solutions of partial differential equations enable a broad range of scientific research. The Dedalus Project is a flexible, open-source, parallelized computational framework for solving general partial differential equations using spectral methods. Dedalus translates plain-text strings describing partial differential equations into efficient solvers. This paper details the numerical meth… ▽ More Numerical solutions of partial differential equations enable a broad range of scientific research. The Dedalus Project is a flexible, open-source, parallelized computational framework for solving general partial differential equations using spectral methods. Dedalus translates plain-text strings describing partial differential equations into efficient solvers. This paper details the numerical method that enables this translation, describes the design and implementation of the codebase, and illustrates its capabilities with a variety of example problems. The numerical method is a first-order generalized tau formulation that discretizes equations into banded matrices. This method is implemented with an object-oriented design. Classes for spectral bases and domains manage the discretization and automatic parallel distribution of variables. Discretized fields and mathematical operators are symbolically manipulated with a basic computer algebra system. Initial value, boundary value, and eigenvalue problems are efficiently solved using high-performance linear algebra, transform, and parallel communication libraries. Custom analysis outputs can also be specified in plain text and stored in self-describing portable formats. The performance of the code is evaluated with a parallel scaling benchmark and a comparison to a finite-volume code. The features and flexibility of the codebase are illustrated by solving several examples: the nonlinear Schrodinger equation on a graph, a supersonic magnetohydrodynamic vortex, quasigeostrophic flow, Stokes flow in a cylindrical annulus, normal modes of a radiative atmosphere, and diamagnetic levitation. The Dedalus code and the example problems are available online at http://dedalus-project.org/. △ Less

Submitted 2 December, 2019; v1 submitted 24 May, 2019; originally announced May 2019.

Comments: 40 pages, 18 figures. Accepted to Physical Review Research

Journal ref: Phys. Rev. Research 2, 023068 (2020)

The "Sphered Cube": A New Method for the Solution of Partial Differential Equations in Cubical Geometry

Authors: Keaton J. Burns, Daniel Lecoanet, Geoffrey M. Vasil, Jeffrey S. Oishi, Benjamin P. Brown

Abstract: A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere. By choosing to embed the cube inside of the sphere, one obtains a discretization that is free of any sharp edges or corners. Taking full advantage of the simpl… ▽ More A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere. By choosing to embed the cube inside of the sphere, one obtains a discretization that is free of any sharp edges or corners. Taking full advantage of the simple geometry of the sphere, spectral bases based on spin-weighted spherical harmonics and Jacobi polynomials, which properly capture the regularity of scalar, vector and tensor components in spherical coordinates, can be applied to obtain moderately efficient and accurate numerical solutions of partial differential equations in the cube. This technique demonstrates the advantages of these bases over other methods for solving PDEs in spherical coordinates. We present results for a test case of incompressible hydrodynamics in cubical geometry: Rayleigh-Bénard convection with fully Dirichlet boundary conditions. Analysis of the simulations provides what is, to our knowledge, the first result on the scaling of the heat flux with the thermal forcing for this type of convection in a cube in a sphere. △ Less

Submitted 29 March, 2019; originally announced March 2019.

Comments: 10 pages, 3 figures, 1 cube, 1 sphere

Breezing through the space environment of Barnard's Star b

Authors: Julián D. Alvarado-Gómez, Cecilia Garraffo, Jeremy J. Drake, Benjamin P. Brown, Jeffrey S. Oishi, Sofia P. Moschou, Ofer Cohen

Abstract: A physically realistic stellar wind model based on Alfvén wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distribu… ▽ More A physically realistic stellar wind model based on Alfvén wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distributions of three different stars that constitute admissible magnetic proxies of this object. Under these considerations, Barnard's Star b experiences less intense wind pressure than the much more close-in planet Proxima~b and the planets of the TRAPPIST-1 system. The milder wind conditions are more a result of its much greater orbital distance rather than in differences in the surface magnetic field strengths of the host stars. The dynamic pressure experienced by the planet is comparable to present-day Earth values, but it can undergo variations by factors of several during current sheet crossings in each orbit. The magnetospause standoff distance would be $\sim$\,$20 - 40$\,\% smaller than that of the Earth for an equivalent planetary magnetic field strength. △ Less

Submitted 29 March, 2019; v1 submitted 1 January, 2019; originally announced January 2019.

Comments: J. D. Alvarado-Gómez & C. Garraffo contributed equally to this work. 8 Pages, 4 Figures. Accepted to The Astrophysical Journal Letters

Predicting the Rossby number in convective experiments

Authors: Evan H. Anders, Cathryn M. Manduca, Benjamin P. Brown, Jeffrey S. Oishi, Geoffrey M. Vasil

Abstract: The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of… ▽ More The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of common inputs to numerical models. The Predictive Rossby number can be specified independent of Rayleigh number, allowing suites of numerical solutions to separate the degree of rotational constraint from the strength of the driving of convection. We examine the scaling of convective transport in terms of the Nusselt number and the degree of turbulence in terms of the Reynolds number of the flow, and we find scaling laws nearly identical to those in nonrotational convection at low Rossby number when the Predictive Rossby number is held constant. Finally, we describe the boundary layers as a function of increasing turbulence at constant Rossby number. △ Less

Submitted 15 January, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: As accepted in ApJ. 9 pages (including references). 4 figures, 2 tables

Journal ref: The Astrophysical Journal. Volume 872, Number 2 (2019)

Accelerated evolution of convective simulations

Authors: Evan H. Anders, Benjamin P. Brown, Jeffrey S. Oishi

Abstract: High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in… ▽ More High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in a simulation and using its characteristics to adjust the mean thermodynamic profile within the domain towards its evolved state. We study Rayleigh-Bénard convection as a test case for AE. Evolved flow properties of AE solutions are measured to be within a few percent of solutions which are reached through standard evolution (SE) over a full thermal timescale. At the highest values of the Rayleigh number at which we compare SE and AE, we find that AE solutions require roughly an order of magnitude fewer computing hours to evolve than SE solutions. △ Less

Submitted 21 August, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: As published in PRFluids. 14 pages total with 6 figures, 2 tables, and 2 appendices

Journal ref: Phys. Rev. Fluids 3, 083502 (2018)

Tensor calculus in spherical coordinates using Jacobi polynomials, Part-II: Implementation and Examples

Authors: Daniel Lecoanet, Geoffrey M. Vasil, Keaton J. Burns, Benjamin P. Brown, Jeffrey S. Oishi

Abstract: We present a simulation code which can solve broad ranges of partial differential equations in a full sphere. The code expands tensorial variables in a spectral series of spin-weighted spherical harmonics in the angular directions and a scaled Jacobi polynomial basis in the radial direction, as described in Part-I. Nonlinear terms are calculated by transforming from the coefficients in the spectra… ▽ More We present a simulation code which can solve broad ranges of partial differential equations in a full sphere. The code expands tensorial variables in a spectral series of spin-weighted spherical harmonics in the angular directions and a scaled Jacobi polynomial basis in the radial direction, as described in Part-I. Nonlinear terms are calculated by transforming from the coefficients in the spectral series to the value of each quantity on the physical grid, where it is easy to calculate products and perform other local operations. The expansion makes it straightforward to solve equations in tensor form (i.e., without decomposition into scalars). We propose and study several unit tests which demonstrate the code can accurately solve linear problems, implement boundary conditions, and transform between spectral and physical space. We then run a series of benchmark problems proposed in Marti et al (2014), implementing the hydrodynamic and magnetohydrodynamic equations. We are able to calculate more accurate solutions than reported in Marti et al 2014 by running at higher spatial resolution and using a higher-order timestepping scheme. We find the rotating convection and convective dynamo benchmark problems depend sensitively on details of timestepping and data analysis. We also demonstrate that in low resolution simulations of the dynamo problem, small changes in a numerical scheme can lead to large changes in the solution. To aid future comparison to these benchmarks, we include the source code used to generate the data, as well as the data and analysis scripts used to generate the figures. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: Submitted to JCP

Convective dynamics and disequilibrium chemistry in the atmospheres of giant planets and brown dwarfs

Authors: Baylee Bordwell, Benjamin P. Brown, Jeffrey S. Oishi

Abstract: Disequilibrium chemical processes have a large effect upon the spectra of substellar objects. To study these effects, dynamical disequilibrium has been parameterized using the quench and eddy diffusion approximations, but little work has been done to explore how these approximations perform under realistic planetary conditions in different dynamical regimes. As a first step in addressing this prob… ▽ More Disequilibrium chemical processes have a large effect upon the spectra of substellar objects. To study these effects, dynamical disequilibrium has been parameterized using the quench and eddy diffusion approximations, but little work has been done to explore how these approximations perform under realistic planetary conditions in different dynamical regimes. As a first step in addressing this problem, we study the localized, small scale convective dynamics of planetary atmospheres by direct numerical simulation of fully compressible hydrodynamics with reactive tracers using the Dedalus code. Using polytropically-stratified, plane parallel atmospheres in 2- and 3-D, we explore the quenching behavior of different abstract chemical species as a function of the dynamical conditions of the atmosphere as parameterized by the Rayleigh number. We find that in both 2- and 3-D, chemical species quench deeper than would be predicted based on simple mixing length arguments. Instead, it is necessary to employ length scales based on the chemical equilibrium profile of the reacting species in order to predict quench points and perform chemical kinetics modeling in 1-D. Based on the results of our simulations, we provide a new length scale, derived from the chemical scale height, which can be used to perform these calculations. This length scale is simple to calculate from known chemical data and makes reasonable predictions for our dynamical simulations. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: 12 pages, 9 figures; Published in the Astrophysical Journal

Journal ref: Baylee Bordwell et al 2018 ApJ 854 8

Perspectives on Reproducibility and Sustainability of Open-Source Scientific Software from Seven Years of the Dedalus Project

Authors: Jeffrey S. Oishi, Benjamin P. Brown, Keaton J. Burns, Daniel Lecoanet, Geoffrey M. Vasil

Abstract: As the Science Mission Directorate contemplates establishing an open code policy, we consider it timely to share our experiences as the developers of the open-source partial differential equation solver Dedalus. Dedalus is a flexible framework for solving partial differential equations. Its development team primarily uses it for studying stellar and planetary astrophysics. Dedalus was developed or… ▽ More As the Science Mission Directorate contemplates establishing an open code policy, we consider it timely to share our experiences as the developers of the open-source partial differential equation solver Dedalus. Dedalus is a flexible framework for solving partial differential equations. Its development team primarily uses it for studying stellar and planetary astrophysics. Dedalus was developed originally for astrophysical fluid dynamics (AFD), though it has found a much broader user base, including applied mathematicians, plasma physicists, and oceanographers. Here, we will focus on issues related to open-source software from the perspective of AFD. We use the term AFD with the understanding that astrophysics simulations are inherently multi-physics: fluid dynamics coupled with some combination of gravitational dynamics, radiation transfer, relativity, and magnetic fields. In practice, a few well-known open-source simulation packages represent a large fraction of published work in the field. However, we will argue that an open-code policy should encompass not just these large simulation codes, but also the input files and analysis scripts. It is our interest that NASA adopt an open-code policy because without it, reproducibility in computational science is needlessly hampered. △ Less

Submitted 24 January, 2018; originally announced January 2018.

Comments: A white paper submitted to the National Academies of Sciences, Engineering, and Medicine's Best Practices for a Future Open Code Policy for NASA Space Science

Convective heat transport in stratified atmospheres at low and high Mach number

Authors: Evan H. Anders, Benjamin P. Brown

Abstract: We study fully compressible convection in the context of plane-parallel, polytropically stratified atmospheres. We perform a suite of 2D and 3D simulations in which we vary the initial superadiabaticity ($ε$) and the Rayleigh number (Ra) while fixing the initial density stratification, aspect ratio, and Prandtl number. The evolved heat transport, quantified by the Nusselt number (Nu), follows scal… ▽ More We study fully compressible convection in the context of plane-parallel, polytropically stratified atmospheres. We perform a suite of 2D and 3D simulations in which we vary the initial superadiabaticity ($ε$) and the Rayleigh number (Ra) while fixing the initial density stratification, aspect ratio, and Prandtl number. The evolved heat transport, quantified by the Nusselt number (Nu), follows scaling relationships similar to those found in the well-studied, incompressible Rayleigh-Bénard problem. This scaling holds up in both 2D and 3D and is not appreciably affected by the magnitude of $ε$. △ Less

Submitted 29 August, 2017; v1 submitted 20 November, 2016; originally announced November 2016.

Comments: 8 pages, 4 figures, as published in PRFluids

Journal ref: Phys. Rev. Fluids 2, 083501 (2017)

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

Authors: Daniel Lecoanet, Josiah Schwab, Eliot Quataert, Lars Bildsten, F. X. Timmes, Keaton J. Burns, Geoffrey M. Vasil, Jeffrey S. Oishi, Benjamin P. Brown

Abstract: It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much l… ▽ More It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution. △ Less

Submitted 25 September, 2016; v1 submitted 29 March, 2016; originally announced March 2016.

Comments: Accepted to ApJ

M Dwarf Flare Continuum Variations on One-Second Timescales: Calibrating and Modeling of ULTRACAM Flare Color Indices

Authors: Adam F. Kowalski, Mihalis Mathioudakis, Suzanne L. Hawley, John P. Wisniewski, Vik S. Dhillon, Tom R. Marsh, Eric J. Hilton, Benjamin P. Brown

Abstract: We present a large dataset of high cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the NUV and optical continuum spectral evolution on timescales of ~1 second. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously-obtained l… ▽ More We present a large dataset of high cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the NUV and optical continuum spectral evolution on timescales of ~1 second. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously-obtained low-resolution spectra during two moderate-sized flares in the dM4.5e star YZ CMi. By avoiding the spectral complexity within the broadband Johnson filters, the ULTRACAM filters are shown to characterize bona-fide continuum emission in the NUV, blue, and red wavelength regimes. The NUV/blue flux ratio in flares is equivalent to a Balmer jump ratio, and the blue/red flux ratio provides an estimate for the color temperature of the optical continuum emission. We present a new "color-color" relationship for these continuum flux ratios at the peaks of the flares. Using the RADYN and RH codes, we interpret the ULTRACAM filter emission using the dominant emission processes from a radiative-hydrodynamic flare model with a high nonthermal electron beam flux, which explains a hot, T~10,000 K, color temperature at blue-to-red optical wavelengths and a small Balmer jump ratio as are observed in moderate-sized and large flares alike. We also discuss the high time-resolution, high signal-to-noise continuum color variations observed in YZ CMi during a giant flare, which increased the NUV flux from this star by over a factor of 100. △ Less

Submitted 15 February, 2016; originally announced February 2016.

Comments: 88 pages, 15 figures, 6 tables, accepted for publication in the Astrophysical Journal, data available online at the Zenodo repository

Tensor calculus in polar coordinates using Jacobi polynomials

Authors: Geoffrey M. Vasil, Keaton J. Burns, Daniel Lecoanet, Sheehan Olver, Benjamin P. Brown, Jeffrey S. Oishi

Abstract: Spectral methods are an efficient way to solve partial differential equations on domains possessing certain symmetries. The utility of a method depends strongly on the choice of spectral basis. In this paper we describe a set of bases built out of Jacobi polynomials, and associated operators for solving scalar, vector, and tensor partial differential equations in polar coordinates on a unit disk.… ▽ More Spectral methods are an efficient way to solve partial differential equations on domains possessing certain symmetries. The utility of a method depends strongly on the choice of spectral basis. In this paper we describe a set of bases built out of Jacobi polynomials, and associated operators for solving scalar, vector, and tensor partial differential equations in polar coordinates on a unit disk. By construction, the bases satisfy regularity conditions at r=0 for any tensorial field. The coordinate singularity in a disk is a prototypical case for many coordinate singularities. The work presented here extends to other geometries. The operators represent covariant derivatives, multiplication by azimuthally symmetric functions, and the tensorial relationship between fields. These arise naturally from relations between classical orthogonal polynomials, and form a Heisenberg algebra. Other past work uses more specific polynomial bases for solving equations in polar coordinates. The main innovation in this paper is to use a larger set of possible bases to achieve maximum bandedness of linear operations. We provide a series of applications of the methods, illustrating their ease-of-use and accuracy. △ Less

Submitted 15 August, 2016; v1 submitted 25 September, 2015; originally announced September 2015.

Comments: 48 pages, 8 figures. Accepted for publication in the Journal of Computational Physics

A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics

Authors: Daniel Lecoanet, Michael McCourt, Eliot Quataert, Keaton J. Burns, Geoffrey M. Vasil, Jeffrey S. Oishi, Benjamin P. Brown, James M. Stone, Ryan M. O'Leary

Abstract: The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack… ▽ More The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form. △ Less

Submitted 27 December, 2021; v1 submitted 11 September, 2015; originally announced September 2015.

Comments: Reference solution snapshots can be found at https://doi.org/10.5281/zenodo.5803759

Numerical Simulations of Internal Wave Generation by Convection in Water

Authors: Daniel Lecoanet, Michael Le Bars, Keaton J. Burns, Geoffrey M. Vasil, Benjamin P. Brown, Eliot Quataert, Jeffrey S. Oishi

Abstract: Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory… ▽ More Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. In order to isolate the physical mechanism exciting internal waves, we use the novel spectral code Dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations over excite high frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal wave generation by convection. △ Less

Submitted 17 June, 2015; v1 submitted 9 December, 2014; originally announced December 2014.

Comments: Accepted to PRE

Conduction in low Mach number flows: Part I Linear & weakly nonlinear regimes

Authors: Daniel Lecoanet, Benjamin P. Brown, Ellen G. Zweibel, Keaton J. Burns, Jeffrey S. Oishi, Geoffrey M. Vasil

Abstract: Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law - the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion - is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entro… ▽ More Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law - the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion - is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entropy. This paper investigates the differences between temperature and entropy diffusion for both linear internal gravity waves and weakly nonlinear convection. In addition to simulating the two thermal conduction models with the fully compressible Navier-Stokes equations, we also study their effects in the reduced, "sound-proof" anelastic and pseudo-incompressible equations. We find that in the linear and weakly nonlinear regimes, temperature and entropy diffusion give quantitatively similar results, although there are some larger errors in the pseudo-incompressible equations with temperature diffusion due to inaccuracies in the equation of state. Extrapolating our weakly nonlinear results, we speculate that differences between temperature and entropy diffusion might become more important for strongly turbulent convection. △ Less

Submitted 20 October, 2014; originally announced October 2014.

Comments: 17 pages, accepted to ApJ

Properties of 42 Solar-type Kepler Targets from the Asteroseismic Modeling Portal

Authors: T. S. Metcalfe, O. L. Creevey, G. Dogan, S. Mathur, H. Xu, T. R. Bedding, W. J. Chaplin, J. Christensen-Dalsgaard, C. Karoff, R. Trampedach, O. Benomar, B. P. Brown, D. L. Buzasi, T. L. Campante, Z. Celik, M. S. Cunha, G. R. Davies, S. Deheuvels, A. Derekas, M. P. Di Mauro, R. A. Garcia, J. A. Guzik, R. Howe, K. B. MacGregor, A. Mazumdar , et al. (17 additional authors not shown)

Abstract: Recently the number of main-sequence and subgiant stars exhibiting solar-like oscillations that are resolved into individual mode frequencies has increased dramatically. While only a few such data sets were available for detailed modeling just a decade ago, the Kepler mission has produced suitable observations for hundreds of new targets. This rapid expansion in observational capacity has been acc… ▽ More Recently the number of main-sequence and subgiant stars exhibiting solar-like oscillations that are resolved into individual mode frequencies has increased dramatically. While only a few such data sets were available for detailed modeling just a decade ago, the Kepler mission has produced suitable observations for hundreds of new targets. This rapid expansion in observational capacity has been accompanied by a shift in analysis and modeling strategies to yield uniform sets of derived stellar properties more quickly and easily. We use previously published asteroseismic and spectroscopic data sets to provide a uniform analysis of 42 solar-type Kepler targets from the Asteroseismic Modeling Portal (AMP). We find that fitting the individual frequencies typically doubles the precision of the asteroseismic radius, mass and age compared to grid-based modeling of the global oscillation properties, and improves the precision of the radius and mass by about a factor of three over empirical scaling relations. We demonstrate the utility of the derived properties with several applications. △ Less

Submitted 29 September, 2014; v1 submitted 14 February, 2014; originally announced February 2014.

Comments: 12 emulateapj pages, 9 figures, 1 online-only extended figure, 1 table, ApJS accepted (typo corrected in Eq.8)

Journal ref: Astrophys.J.Supp.214:27,2014

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

Authors: Geoffrey M. Vasil, Daniel Lecoanet, Benjamin P. Brown, Toby S. Wood, Ellen G. Zweibel

Abstract: 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. Ge… ▽ More 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. △ Less

Submitted 28 February, 2013; originally announced March 2013.

Comments: 23 pages, submitted to ApJ

Buoyant Magnetic Loops Generated by Global Convective Dynamo Action

Authors: Nicholas J. Nelson, Benjamin P. Brown, A. Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: Our global 3D simulations of convection and dynamo action in a Sun-like star reveal that persistent wreaths of strong magnetism can be built within the bulk of the convention zone. Here we examine the characteristics of buoyant magnetic structures that are self-consistently created by dynamo action and turbulent convective motions in a simulation with solar stratification but rotating at three tim… ▽ More Our global 3D simulations of convection and dynamo action in a Sun-like star reveal that persistent wreaths of strong magnetism can be built within the bulk of the convention zone. Here we examine the characteristics of buoyant magnetic structures that are self-consistently created by dynamo action and turbulent convective motions in a simulation with solar stratification but rotating at three times the current solar rate. These buoyant loops originate within sections of the magnetic wreaths in which turbulent flows amplify the fields to much larger values than is possible through laminar processes. These amplified portions can rise through the convective layer by a combination of magnetic buoyancy and advection by convective giant cells, forming buoyant loops. We measure statistical trends in the polarity, twist, and tilt of these loops. Loops are shown to preferentially arise in longitudinal patches somewhat reminiscent of active longitudes in the Sun, although broader in extent. We show that the strength of the axisymmetric toroidal field is not a good predictor of the production rate for buoyant loops or the amount of magnetic flux in the loops that are produced. △ Less

Submitted 21 December, 2012; originally announced December 2012.

Comments: Solar Physics format, 22 pages, accepted by Solar Physics for topical issue related NSO Workshop #26: Solar Origins of Space Weather and Space Climate

Magnetic Activity Cycles in the Exoplanet Host Star epsilon Eridani

Authors: T. S. Metcalfe, A. P. Buccino, B. P. Brown, S. Mathur, D. R. Soderblom, T. J. Henry, P. J. D. Mauas, R. Petrucci, J. C. Hall, S. Basu

Abstract: The active K2 dwarf epsilon Eri has been extensively characterized, both as a young solar analog and more recently as an exoplanet host star. As one of the nearest and brightest stars in the sky, it provides an unparalleled opportunity to constrain stellar dynamo theory beyond the Sun. We confirm and document the 3 year magnetic activity cycle in epsilon Eri originally reported by Hatzes and cowor… ▽ More The active K2 dwarf epsilon Eri has been extensively characterized, both as a young solar analog and more recently as an exoplanet host star. As one of the nearest and brightest stars in the sky, it provides an unparalleled opportunity to constrain stellar dynamo theory beyond the Sun. We confirm and document the 3 year magnetic activity cycle in epsilon Eri originally reported by Hatzes and coworkers, and we examine the archival data from previous observations spanning 45 years. The data show coexisting 3 year and 13 year periods leading into a broad activity minimum that resembles a Maunder minimum-like state, followed by the resurgence of a coherent 3 year cycle. The nearly continuous activity record suggests the simultaneous operation of two stellar dynamos with cycle periods of 2.95+/-0.03 years and 12.7+/-0.3 years, which by analogy with the solar case suggests a revised identification of the dynamo mechanisms that are responsible for the so-called "active" and "inactive" sequences as proposed by Bohm-Vitense. Finally, based on the observed properties of epsilon Eri we argue that the rotational history of the Sun is what makes it an outlier in the context of magnetic cycles observed in other stars (as also suggested by its Li depletion), and that a Jovian-mass companion cannot be the universal explanation for the solar peculiarities. △ Less

Submitted 18 December, 2012; originally announced December 2012.

Comments: 6 pages, 3 figures, 1 table, ApJ Letters (accepted)

Journal ref: Astrophys.J.763:L26,2013

Magnetic Wreaths and Cycles in Convective Dynamos

Authors: Nicholas J. Nelson, Benjamin P. Brown, A. Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: Solar-type stars exhibit a rich variety of magnetic activity. Seeking to explore the convective origins of this activity, we have carried out a series of global 3D magnetohydrodynamic (MHD) simulations with the anelastic spherical harmonic (ASH) code. Here we report on the dynamo mechanisms achieved as the effects of artificial diffusion are systematically decreased. The simulations are carried ou… ▽ More Solar-type stars exhibit a rich variety of magnetic activity. Seeking to explore the convective origins of this activity, we have carried out a series of global 3D magnetohydrodynamic (MHD) simulations with the anelastic spherical harmonic (ASH) code. Here we report on the dynamo mechanisms achieved as the effects of artificial diffusion are systematically decreased. The simulations are carried out at a nominal rotation rate of three times the solar value (3$Ω_\odot$), but similar dynamics may also apply to the Sun. Our previous simulations demonstrated that convective dynamos can build persistent toroidal flux structures (magnetic wreaths) in the midst of a turbulent convection zone and that high rotation rates promote the cyclic reversal of these wreaths. Here we demonstrate that magnetic cycles can also be achieved by reducing the diffusion, thus increasing the Reynolds and magnetic Reynolds numbers. In these more turbulent models, diffusive processes no longer play a significant role in the key dynamical balances that establish and maintain the differential rotation and magnetic wreaths. Magnetic reversals are attributed to an imbalance in the poloidal magnetic induction by convective motions that is stabilized at higher diffusion levels. Additionally, the enhanced levels of turbulence lead to greater intermittency in the toroidal magnetic wreaths, promoting the generation of buoyant magnetic loops that rise from the deep interior to the upper regions of our simulated domain. The implications of such turbulence-induced magnetic buoyancy for solar and stellar flux emergence are also discussed. △ Less

Submitted 13 November, 2012; originally announced November 2012.

Comments: 21 pages, 16 figures, accepted for publication in ApJ

Optimized boundary driven flows for dynamos in a sphere

Authors: I. V. Khalzov, B. P. Brown, C. M. Cooper, D. B. Weisberg, C. B. Forest

Abstract: We perform numerical optimization of the axisymmetric flows in a sphere to minimize the critical magnetic Reynolds number Rm_cr required for dynamo onset. The optimization is done for the class of laminar incompressible flows of von Karman type satisfying the steady-state Navier-Stokes equation. Such flows are determined by equatorially antisymmetric profiles of driving azimuthal (toroidal) veloci… ▽ More We perform numerical optimization of the axisymmetric flows in a sphere to minimize the critical magnetic Reynolds number Rm_cr required for dynamo onset. The optimization is done for the class of laminar incompressible flows of von Karman type satisfying the steady-state Navier-Stokes equation. Such flows are determined by equatorially antisymmetric profiles of driving azimuthal (toroidal) velocity specified at the spherical boundary. The model is relevant to the Madison plasma dynamo experiment (MPDX), whose spherical boundary is capable of differential driving of plasma in the azimuthal direction. We show that the dynamo onset in this system depends strongly on details of the driving velocity profile and the fluid Reynolds number Re. It is found that the overall lowest Rm_cr~200 is achieved at Re~240 for the flow, which is hydrodynamically marginally stable. We also show that the optimized flows can sustain dynamos only in the range Rm_cr<Rm<Rm_cr2, where Rm_cr2 is the second critical magnetic Reynolds number, above which the dynamo is quenched. Samples of the optimized flows and the corresponding dynamo fields are presented. △ Less

Submitted 8 November, 2012; originally announced November 2012.

Journal ref: Physics of Plasmas 19, 112106 (2012)

Energy Conservation and Gravity Waves in Sound-proof Treatments of Stellar Interiors: Part I Anelastic Approximations

Authors: Benjamin P Brown, Geoffrey M Vasil, Ellen G Zweibel

Abstract: Typical flows in stellar interiors are much slower than the speed of sound. To follow the slow evolution of subsonic motions, various sound-proof equations are in wide use, particularly in stellar astrophysical fluid dynamics. These low-Mach number equations include the anelastic equations. Generally, these equations are valid in nearly adiabatically stratified regions like stellar convection zone… ▽ More Typical flows in stellar interiors are much slower than the speed of sound. To follow the slow evolution of subsonic motions, various sound-proof equations are in wide use, particularly in stellar astrophysical fluid dynamics. These low-Mach number equations include the anelastic equations. Generally, these equations are valid in nearly adiabatically stratified regions like stellar convection zones, but may not be valid in the sub-adiabatic, stably stratified stellar radiative interiors. Understanding the coupling between the convection zone and the radiative interior is a problem of crucial interest and may have strong implications for solar and stellar dynamo theories as the interface between the two, called the tachocline in the Sun, plays a crucial role in many solar dynamo theories. Here we study the properties of gravity waves in stably-stratified atmospheres. In particular, we explore how gravity waves are handled in various sound-proof equations. We find that some anelastic treatments fail to conserve energy in stably-stratified atmospheres, instead conserving pseudo-energies that depend on the stratification, and we demonstrate this numerically. One anelastic equation set does conserve energy in all atmospheres and we provide recommendations for converting low-Mach number anelastic codes to this set of equations. △ Less

Submitted 11 July, 2012; originally announced July 2012.

Comments: Accepted for publication in ApJ. 20 pages emulateapj format, 7 figures

Journal ref: 2012, ApJ, 756, 109

Convective Babcock-Leighton Dynamo Models

Authors: Mark S. Miesch, Benjamin P. Brown

Abstract: We present the first global, three-dimensional simulations of solar/stellar convection that take into account the influence of magnetic flux emergence by means of the Babcock-Leighton (BL) mechanism. We have shown that the inclusion of a BL poloidal source term in a convection simulation can promote cyclic activity in an otherwise steady dynamo. Some cycle properties are reminiscent of solar obser… ▽ More We present the first global, three-dimensional simulations of solar/stellar convection that take into account the influence of magnetic flux emergence by means of the Babcock-Leighton (BL) mechanism. We have shown that the inclusion of a BL poloidal source term in a convection simulation can promote cyclic activity in an otherwise steady dynamo. Some cycle properties are reminiscent of solar observations, such as the equatorward propagation of toroidal flux near the base of the convection zone. However, the cycle period in this young sun (rotating three times faster than the solar rate) is very short ($\sim$ 6 months) and it is unclear whether much longer cycles may be achieved within this modeling framework, given the high efficiency of field generation and transport by the convection. Even so, the incorporation of mean-field parameterizations in 3D convection simulations to account for elusive processes such as flux emergence may well prove useful in the future modeling of solar and stellar activity cycles. △ Less

Submitted 12 January, 2012; originally announced January 2012.

Comments: 5 pages, 4 figures, submitted to the Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, 746:L26 (5pp) 2012 February 20

Convection and Differential Rotation Properties of G & K Stars Computed with the ASH Code

Authors: Sean P. Matt, Olivier Do Cao, B. P. Brown, A. S. Brun

Abstract: The stellar luminosity and depth of the convective envelope vary rapidly with mass for G- and K-type main sequence stars. In order to understand how these properties influence the convective turbulence, differential rotation, and meridional circulation, we have carried out 3D dynamical simulations of the interiors of rotating main sequence stars, using the anelastic spherical harmonic (ASH) code.… ▽ More The stellar luminosity and depth of the convective envelope vary rapidly with mass for G- and K-type main sequence stars. In order to understand how these properties influence the convective turbulence, differential rotation, and meridional circulation, we have carried out 3D dynamical simulations of the interiors of rotating main sequence stars, using the anelastic spherical harmonic (ASH) code. The stars in our simulations have masses of 0.5, 0.7, 0.9, and 1.1 M_sun, corresponding to spectral types K7 through G0, and rotate at the same angular speed as the sun. We identify several trends of convection zone properties with stellar mass, exhibited by the simulations. The convective velocities, temperature contrast between up- and down-flows, and meridional circulation velocities all increase with stellar luminosity. As a consequence of the trend in convective velocity, the Rossby number (at a fixed rotation rate) increases and the convective turnover timescales decrease significantly with increasing stellar mass. The 3 lowest mass cases exhibit solar-like differential rotation, in a sense that they show a maximum rotation at the equator and minimum at higher latitudes, but the 1.1 M_sun case exhibits anti-solar rotation. At low mass, the meridional circulation is multi-cellular and aligned with the rotation axis; as the mass increases, the circulation pattern tends toward a unicellular structure covering each hemisphere in the convection zone. △ Less

Submitted 23 November, 2011; originally announced November 2011.

Comments: Accepted for publication in Astronomische Nachrichten; 11 pages, 8 figures

Buoyant Magnetic Loops in a Global Dynamo Simulation of a Young Sun

Authors: Nicholas J. Nelson, Benjamin P. Brown, Allan Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are believed to arise when sections of these magnetic structures become buoyantly unstable and rise from the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic (MHD… ▽ More The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are believed to arise when sections of these magnetic structures become buoyantly unstable and rise from the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic (MHD) simulation in which turbulent convection, stratification, and rotation combine to yield a dynamo that self-consistently generates buoyant magnetic loops. We simulate stellar convection and dynamo action in a spherical shell with solar stratification, but rotating three times faster than the current solar rate. Strong wreaths of toroidal magnetic field are realized by dynamo action in the convection zone. By turning to a dynamic Smagorinsky model for subgrid-scale turbulence, we here attain considerably reduced diffusion in our simulation. This permits the regions of strongest magnetic field in these wreaths to rise toward the top of the convection zone via a combination of magnetic buoyancy instabilities and advection by convective giant cells. Such a global simulation yielding buoyant loops represents a significant step forward in combining numerical models of dynamo action and flux emergence. △ Less

Submitted 23 August, 2011; originally announced August 2011.

Comments: 6 pages, 3 figures, accepted in ApJ Letters

Global-scale Simulations of Stellar Convection and their Observational Predictions

Authors: Benjamin P. Brown

Abstract: Stars on the lower main sequence (F-type through M-type) have substantial convective envelopes beneath their stellar photospheres. Convection in these regions can couple with rotation to build global-scale structures that may be observable by interferometers that can resolve stellar disks. Here I discuss predictions emerging from 3D MHD simulations for solar-type stars with the anelastic spherical… ▽ More Stars on the lower main sequence (F-type through M-type) have substantial convective envelopes beneath their stellar photospheres. Convection in these regions can couple with rotation to build global-scale structures that may be observable by interferometers that can resolve stellar disks. Here I discuss predictions emerging from 3D MHD simulations for solar-type stars with the anelastic spherical harmonic (ASH) code and how these predictions may be observationally tested. The zonal flow of differential rotation is likely the most easily observable signature of dynamics occurring deep within the stellar interior. Generally, we find that rapidly rotating suns have a strong solar-like differential rotation with a prograde equator and retrograde poles while slowly spinning suns may have anti-solar rotation profiles with fast poles and slow equators. The thermal wind balance accompanying the differential rotation may lead to hot and bright poles in the rapid rotators and cooler, darker poles in slow rotators. The convection and differential rotation build global-scale magnetic structures in the bulk of the convection zone, and these wreaths of magnetism may be observable near the stellar surfaces. △ Less

Submitted 29 June, 2011; originally announced June 2011.

Comments: 9 pages, 4 figures. Conference proceedings for MRO 2011 workshop "Resolving the Future of Astronomy with Long-Baseline Interferometry", Socorro, NM, March 28-31, 2011

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.

Comments: 20 pages, 14 figures, emulateapj format; accepted for publication in ApJ. Expanded and published version of sections 5-6 from http://arxiv.org/abs/0906.2407

Journal ref: ApJ, 731 (2011), 69 (19pp)

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.

Comments: 8 pages, 3 figures; to appear in proceedings of Cool Stars 16, Seattle Aug 29 - Sep 2, 2010

Dynamos in Stellar Convection Zones: of Wreaths and Cycles

Authors: Benjamin P Brown

Abstract: We live near a magnetic star whose cycles of activity are driven by dynamo action beneath the surface. In the solar convection zone, rotation couples with plasma motions to build highly organized magnetic fields that erupt at the surface and undergo relatively regular cycles of polarity reversal. Despite our proximity to the Sun, the nature of its dynamo remains elusive, but observations of other… ▽ More We live near a magnetic star whose cycles of activity are driven by dynamo action beneath the surface. In the solar convection zone, rotation couples with plasma motions to build highly organized magnetic fields that erupt at the surface and undergo relatively regular cycles of polarity reversal. Despite our proximity to the Sun, the nature of its dynamo remains elusive, but observations of other solar-type stars show that surface magnetism is a nearly ubiquitous feature. In recent time, numerical simulations of convection and dynamo action have taken tremendous strides forward. Global-scale organization and cyclic magnetism are being achieved by several groups in distinctly different solar and stellar simulations. Here I will talk about advances on the numerical front including wreath-building dynamos which may occupy stellar convection zones. I will discuss the interplay between the new simulations, various classes of mean-field models, and current and upcoming solar and stellar observations. △ Less

Submitted 12 November, 2010; originally announced November 2010.

Comments: 10 pages, 5 figures, a4 format; proceedings for SOHO 24/GONG 2010 conference: "A new era of seismology of the Sun and solar-like stars," Aix-en-Provence, France, June 27-July 4, 2010 (JPCS)

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.

Comments: 17 pages, 9 figures, 1 appendix, emulateapj format; published version of sections 3-4, 7 and appendix from arXiv:0906.2407

Journal ref: 2010, ApJ, 711, 424-438

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.

Comments: 8 pages, 4 figures. To appear in IAU 271: "Astrophysical Dynamics: from Stars to Galaxies"

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.

Comments: 4 pages, 3 figures, from IAU 273: The Physics of Sun and Star Spots

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.

Comments: 8 pages, 2 color figures. To appear in Proc. IAU Symposium 271, "Astrophysical Dynamics, from Stars to Galaxies", ed. A.S. Brun, N.H. Brummell & M.S. Miesch, in press Second version: corrected an error in the caption of Figure 2

Exploring the $P_{cyc}$ vs $P_{rot}$ relation with flux transport dynamo models of solar-like stars

Authors: L. Jouve, B. P. Brown, A. S. Brun

Abstract: Aims: To understand stellar magnetism and to test the validity of the Babcock-Leighton flux transport mean field dynamo models with stellar activity observations Methods: 2-D mean field dynamo models at various rotation rates are computed with the STELEM code to study the sensitivity of the activity cycle period and butterfly diagram to parameter changes and are compared to observational data. T… ▽ More Aims: To understand stellar magnetism and to test the validity of the Babcock-Leighton flux transport mean field dynamo models with stellar activity observations Methods: 2-D mean field dynamo models at various rotation rates are computed with the STELEM code to study the sensitivity of the activity cycle period and butterfly diagram to parameter changes and are compared to observational data. The novelty is that these 2-D mean field dynamo models incorporate scaling laws deduced from 3-D hydrodynamical simulations for the influence of rotation rate on the amplitude and profile of the meridional circulation. These models make also use of observational scaling laws for the variation of differential rotation with rotation rate. Results: We find that Babcock-Leighton flux transport dynamo models are able to reproduce the change in topology of the magnetic field (i.e. toward being more toroidal with increasing rotation rate) but seem to have difficulty reproducing the cycle period vs activity period correlation observed in solar-like stars if a monolithic single cell meridional flow is assumed. It may however be possible to recover the $P_{cyc}$ vs $P_{rot}$ relation with more complex meridional flows, if the profile changes in a particular assumed manner with rotation rate. Conclusions: The Babcock-Leighton flux transport dynamo model based on single cell meridional circulation does not reproduce the $P_{cyc}$ vs $P_{rot}$ relation unless the amplitude of the meridional circulation is assumed to increase with rotation rate which seems to be in contradiction with recent results obtained with 3-D global simulations. △ Less

Submitted 10 November, 2009; originally announced November 2009.

Comments: 12 pages, 8 figures, accepted for publication by A&A 1: AIM, CEA/DSM-CNRS-Univ. Paris 7, IRFU/SAp, France, 2: D.A.M.T.P., Centre for Mathematical Sciences, Univ. of Cambridge, UK, 3: JILA and Department of Astrophysical and Planetary Sciences, Univ. of Colorado, USA

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.

Comments: 22 pages, 16 color figures in emulateapj format. Images heavily degraded, but legible. Submitted to ApJ

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.

Comments: 14 pages, 5 color figures, to appear in Proc. GONG 2008/SOHO XXI Meeting on Solar-Stellar Dynamos as Revealed by Helio and Asteroseismology, held August 15-18, 2008, Boulder, CO, Astronomical Soc. Pac. Conf. Series, volume TBD

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.

Comments: accepted for publication in ApJ. Estimated publication date is December 10, 2008, in volume 689. Images in this version are degraded but legible

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.

Comments: 8 pages, 4 figs. Published in conference proceedings "Unsolved Problems in Stellar Physics", held July 2-6 2007 Cambridge, England

Journal ref: AIP Conf.Proc.948:271-278,2007

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.

Comments: 4 pages, 2 figures. Presentation at AIP 5th Potsdam Thinkshop, June 24-29 2007. Published in Astr. Nachr. Dec 2007

Journal ref: Astron.Nachr.328:1002-1005,2007