Researcher profile

Mark S. Miesch

Mark S. Miesch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
23works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

23 published item(s)

preprint2016arXiv

A Three-Dimensional Babcock-Leighton Solar Dynamo Model: Initial Results with Axisymmetric Flows

The main objective of this paper is to introduce the STABLE (Surface flux Transport And Babcock-LEighton) solar dynamo model. STABLE is a 3D Babcock-Leighton/Flux Transport dynamo model in which the source of poloidal field is the explicit emergence, distortion, and dispersal of bipolar magnetic regions (BMRs). Here we describe the STABLE model in more detail than we have previously and we verify it by reproducing a 2D mean-field benchmark. We also present some representative dynamo simulations, focusing on the special case of kinematic magnetic induction and axisymmetric flow fields. Not all solutions are supercritical; it can be a challenge for the BL mechanism to sustain the dynamo when the turbulent diffusion near the surface is $\geq 10^{12}$ cm$^2$ s$^{-1}$. However, if BMRs are sufficiently large, deep, and numerous, then sustained, cyclic, dynamo solutions can be found that exhibit solar-like features. Furthermore, we find that the shearing of radial magnetic flux by the surface differential rotation can account for most of the net toroidal flux generation in each hemisphere, as has been recently argued for the Sun by Cameron & Schussler (2015).

preprint2016arXiv

Magnetic Helicity Reversals in a Cyclic Convective Dynamo

We investigate the role of magnetic helicity in promoting cyclic magnetic activity in a global, 3D, magnetohydrodynamic (MHD) simulation of a convective dynamo. This simulation is characterized by coherent bands of toroidal field that exist within the convection zone, with opposite polarities in the northern and southern hemispheres. Throughout most of the cycle, the magnetic helicity in these bands is negative in the northern hemisphere and positive in the southern hemisphere. However, during the declining phase of each cycle, this hemispheric rule reverses. We attribute this to a global restructuring of the magnetic topology that is induced by the interaction of the bands across the equator. This band interaction appears to be ultimately responsible for, or at least associated with, the decay and subsequent reversal of both the toroidal bands and the polar fields. We briefly discuss the implications of these results within the context of solar observations, which also show some potential evidence for toroidal band interactions and helicity reversals.

preprint2016arXiv

Velocity Amplitudes in Global Convection Simulations: The Role of the Prandtl Number and Near-Surface Driving

Several lines of evidence suggest that the velocity amplitude in global simulations of solar convection, U, may be systematically over-estimated. Motivated by these recent results, we explore the factors that determine U and we consider how these might scale to solar parameter regimes. To this end, we decrease the thermal diffusivity $κ$ along two paths in parameter space. If the kinematic viscosity $ν$ is decreased proportionally with $κ$ (fixing the Prandtl number $P_r = ν/κ$), we find that U increases but asymptotes toward a constant value, as found by Featherstone & Hindman (2016). However, if $ν$ is held fixed while decreasing $κ$ (increasing $P_r$), we find that U systematically decreases. We attribute this to an enhancement of the thermal content of downflow plumes, which allows them to carry the solar luminosity with slower flow speeds. We contrast this with the case of Rayleigh-Benard convection which is not subject to this luminosity constraint. This dramatic difference in behavior for the two paths in parameter space (fixed $P_r$ or fixed $ν$) persists whether the heat transport by unresolved, near-surface convection is modeled as a thermal conduction or as a fixed flux. The results suggest that if solar convection can operate in a high-$P_r$ regime, then this might effectively limit the velocity amplitude. Small-scale magnetism is a possible source of enhanced viscosity that may serve to achieve this high-$P_r$ regime.

preprint2015arXiv

A Compressible High-Order Unstructured Spectral Difference Code for Stratified Convection in Rotating Spherical Shells

We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are treated as rotating spherical shells filled with stratified gas. The hydrodynamic equations are discretized by a robust and efficient high-order Spectral Difference Method (SDM) on unstructured meshes. The computational stencil of the spectral difference method is compact and advantageous for parallel processing. CHORUS demonstrates excellent parallel performance for all test cases reported in this paper, scaling up to 12,000 cores on the Yellowstone High-Performance Computing cluster at NCAR. The code is verified by defining two benchmark cases for global convection in Jupiter and the Sun. CHORUS results are compared with results from the ASH code and good agreement is found. The CHORUS code creates new opportunities for simulating such varied phenomena as multi-scale solar convection, core convection, and convection in rapidly-rotating, oblate stars.

preprint2015arXiv

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

preprint2015arXiv

Meridional Circulation in Solar and Stellar Convection Zones

We present a series of 3-D nonlinear simulations of solar-like convection, carried out using the Anelastic Spherical Harmonic (ASH) code, that are designed to isolate those processes that drive and shape meridional circulations within stellar convection zones. These simulations have been constructed so as to span the transition between solar-like differential rotation (fast equator/slow poles) and ``anti-solar' differential rotation (slow equator/fast poles). Solar-like states of differential rotation, arising when convection is rotationally constrained, are characterized by a very different convective Reynolds stress than anti-solar regimes, wherein convection only weakly senses the Coriolis force. We find that the angular momentum transport by convective Reynolds stress plays a central role in establishing the meridional flow profiles in these simulations. We find that the transition from single-celled to multi-celled meridional circulation profiles in strong and weak regimes of rotational constraint is linked to a change in the convective Reynolds stress, a clear demonstration of gyroscopic pumping. Latitudinal thermal variations differ between these different regimes, with those in the solar-like regime conspiring to suppress a single cell of meridional circulation, whereas the cool poles and warm equator established in the anti-solar states tend to promote single-celled circulations. Though the convective angular momentum transport becomes radially inward at mid-latitudes in anti-solar regimes, it is the meridional circulation that is primarily responsible for establishing a rapidly-rotating pole. We conclude with a discussion of how these results relate to the Sun, and suggest that the Sun may lie near the transition between rapidly-rotating and slowly-rotating regimes.

preprint2014arXiv

A 3D Babcock-Leighton Solar Dynamo Model

We present a 3D kinematic solar dynamo model in which poloidal field is generated by the emergence and dispersal of tilted sunspot pairs (more generally Bipolar Magnetic Regions, or BMRs). The axisymmetric component of this model functions similarly to previous 2D Babcock-Leighton (BL) dynamo models that employ a double-ring prescription for poloidal field generation but we generalize this prescription into a 3D flux emergence algorithm that places BMRs on the surface in response to the dynamo-generated toroidal field. In this way, the model can be regarded as a unification of BL dynamo models (2D in radius/latitude) and surface flux transport models (2D in latitude/longitude) into a more self-consistent framework that captures the full 3D structure of the evolving magnetic field. The model reproduces some basic features of the solar cycle including an 11-yr periodicity, equatorward migration of toroidal flux in the deep convection zone, and poleward propagation of poloidal flux at the surface. The poleward-propagating surface flux originates as trailing flux in BMRs, migrates poleward in multiple non-axisymmetric streams (made axisymmetric by differential rotation and turbulent diffusion), and eventually reverses the polar field, thus sustaining the dynamo. In this letter we briefly describe the model, initial results, and future plans.

preprint2014arXiv

Differential rotation in main-sequence solar-like stars: Qualitative inference from asteroseismic data

Understanding differential rotation of Sun-like stars is of great importance for insight into the angular momentum transport in these stars. One means of gaining such information is that of asteroseismology. By a forward modeling approach we analyze in a qualitative manner the impact of different differential rotation profiles on the splittings of p-mode oscillation frequencies. The optimum modes for inference on differential rotation are identified along with the best value of the stellar inclination angle. We find that in general it is not likely that asteroseismology can be used to make an unambiguous distinction between a rotation profile such as, e.g., a conical Sun-like profile and a cylindrical profile. In addition, it seems unlikely that asteroseismology of Sun-like stars will result in inferences on the radial profile of the differential rotation, such as can be done for, e.g., red giants. At best one could possibly obtain the sign of the radial differential rotation gradient. Measurements of the extent of the latitudinal differential from frequency splitting are, however, more promising. One very interesting aspect that could likely be tested from frequency splittings is whether the differential rotation is solar-like or anti-solar-like in nature, in the sense that a solar-like profile has an equator rotating faster than the poles.

preprint2014arXiv

Generating buoyant magnetic flux ropes in solar-like convective dynamos

Our Sun exhibits strong convective dynamo action which results in magnetic flux bundles emerging through the stellar surface as magnetic spots. Global-scale dynamo action is believed to generate large-scale magnetic structures in the deep solar interior through the interplay of convection, rotation, and shear. Portions of these large-scale magnetic structures are then believed to rise through the convective layer, forming magnetic loops which then pierce the photosphere as sunspot pairs. Previous global simulations of 3D MHD convection in rotating spherical shells have demonstrated mechanisms whereby large-scale magnetic wreaths can be generated in the bulk of the convection zone. Our recent simulations have achieved sufficiently high levels of turbulence to permit portions of these wreaths to become magnetically buoyant and rise through the simulated convective layer through a combination of magnetic buoyancy and advection by convective giant cells. These buoyant magnetic loops are created in the bulk of the convective layer as strong Lorentz force feedback in the cores of the magnetic wreaths dampen small-scale convective motions, permitting the amplification of local magnetic energies to over 100 times the local kinetic energy. While the magnetic wreaths are largely generated the shearing of axisymmetric poloidal magnetic fields by axisymmetric rotational shear (the $Ω$-effect), the loops are amplified to their peak field strengths before beginning to rise by non-axisymmetric processes. This further extends and enhances a new paradigm for the generation of emergent magnetic flux bundles, which we term turbulence-enabled magnetic buoyancy.

preprint2013arXiv

A Theory on the Convective Origins of Active Longitudes on Solar-like Stars

Using a thin flux tube model in a rotating spherical shell of turbulent, solar-like convective flows, we find that the distribution of emerging flux tubes in our simulation is inhomogeneous in longitude, with properties similar to those of active longitudes on the Sun and other solar-like stars. The large-scale pattern of flux emergence our simulations produce exhibits preferred longitudinal modes of low order, drift with respect to a fixed reference system, and alignment across the Equator at low latitudes between 15 degrees. We suggest that these active-longitude-like emergence patterns are the result of columnar, rotationally aligned giant cells present in our convection simulation at low latitudes. If giant convecting cells exist in the bulk of the solar convection zone, this phenomenon, along with differential rotation, could in part provide an explanation for the behavior of active longitudes.

preprint2012arXiv

Buoyant Magnetic Loops Generated by Global Convective Dynamo Action

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.

preprint2012arXiv

Comparing Simulations of Rising Flux Tubes Through the Solar Convection Zone with Observations of Solar Active Regions: Constraining the Dynamo Field Strength

We study how active-region-scale flux tubes rise buoyantly from the base of the convection zone to near the solar surface by embedding a thin flux tube model in a rotating spherical shell of solar-like turbulent convection. These toroidal flux tubes that we simulate range in magnetic field strength from 15 kG to 100 kG at initial latitudes of 1 degree to 40 degrees in both hemispheres. This article expands upon Weber, Fan, and Miesch (Astrophys. J., 741, 11, 2011) (Article 1) with the inclusion of tubes with magnetic flux of 10^20 Mx and 10^21 Mx, and more simulations of the previously investigated case of 10^22 Mx, sampling more convective flows than the previous article, greatly improving statistics. Observed properties of active regions are compared to properties of the simulated emerging flux tubes, including: the tilt of active regions in accordance with Joy's Law as in Article 1, and in addition the scatter of tilt angles about the Joy's Law trend, the most commonly occurring tilt angle, the rotation rate of the emerging loops with respect to the surrounding plasma, and the nature of the magnetic field at the flux tube apex. We discuss how these diagnostic properties constrain the initial field strength of the active region flux tubes at the bottom of the solar convection zone, and suggest that flux tubes of initial magnetic field strengths of \geq 40 kG are good candidates for the progenitors of large (10^21 Mx to 10^22 Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 10^22 Mx. With the addition of more magnetic flux values and more simulations, we find that for all magnetic field strengths, the emerging tubes show a positive Joy's Law trend, and that this trend does not show a statistically significant dependence on the magnetic flux.

preprint2012arXiv

Convective Babcock-Leighton Dynamo Models

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.

preprint2012arXiv

Magnetic Wreaths and Cycles in Convective Dynamos

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.

preprint2012arXiv

On the Amplitude of Convective Velocities in the Deep Solar Interior

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone based only on the observed properties of the differential rotation and meridional circulation together with simple and robust dynamical balances obtained from the fundamental MHD equations. The linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation across isosurfaces of specific angular momentum is linked to the angular momentum transport by the convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 30 m s$^{-1}$ in the upper CZ ($r \sim 0.95 R$) and at least 8 m s$^{-1}$ in the lower CZ ($r \sim 0.75 R$) in order to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a probe of the rotational influence, we are further able to show that the characteristic length scale of deep convective motions must be no smaller than 5.5--30 Mm. These results are compatible with convection models but suggest that the efficiency of the turbulent transport assumed in advection-dominated flux-transport dynamo models is generally not consistent with the mean flows they employ.

preprint2011arXiv

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

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.

preprint2011arXiv

Gyroscopic Pumping in the Solar Near-Surface Shear Layer

We use global and local helioseismic inversions to explore the prevailing dynamical balances in the solar Near-Surface Shear Layer (NSSL). The differential rotation and meridional circulation are intimately linked, with a common origin in the turbulent stresses of the upper solar convection zone. The existence and structure of the NSSL cannot be attributed to the conservation of angular momentum by solar surface convection, as is often supposed. Rather, the turbulent angular momentum transport accounts for the poleward meridional flow while the often overlooked meridional force balance maintains the mid-latitude rotational shear. We suggest that the base of the NSSL is marked by a transition from baroclinic to turbulent stresses in the meridional plane which suppress Coriolis-induced circulations that would otherwise establish a cylindrical rotation profile. The turbulent angular momentum transport must be non-diffusive and directed radially inward. Inferred mean flows are consistent with the idea that turbulent convection tends to mix angular momentum but only if the mixing efficiency is inhomogeneous and/or anisotropic. The latitudinal and longitudinal components of the estimated turbulent transport are comparable in amplitude and about an order of magnitude larger than the vertical component. We estimate that it requires 2--4% of the solar luminosity to maintain the solar NSSL against the inertia of the mean flow. Most of this energy is associated with the turbulent transport of angular momentum out of the layer, with a spin-down time scale of $\sim$ 600 days. We also address implications of these results for numerical modeling of the NSSL.

preprint2011arXiv

Magnetic Cycles in a Convective Dynamo Simulation of a Young Solar-type Star

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]

preprint2011arXiv

The Rise of Active Region Flux Tubes in the Turbulent Solar Convective Envelope

We use a thin flux tube model in a rotating spherical shell of turbulent convective flows to study how active region scale flux tubes rise buoyantly from the bottom of the convection zone to near the solar surface. We investigate toroidal flux tubes at the base of the convection zone with field strengths ranging from 15 kG to 100 kG at initial latitudes ranging from 1 degree to 40 degrees with a total flux of 10^22 Mx. We find that the dynamic evolution of the flux tube changes from convection dominated to magnetic buoyancy dominated as the initial field strength increases from 15 kG to 100 kG. At 100 kG, the development of Omega-shaped rising loops is mainly controlled by the growth of the magnetic buoyancy instability. However, at low field strengths of 15 kG, the development of rising Omega-shaped loops is largely controlled by convective flows, and properties of the emerging loops are significantly changed compared to previous results in the absence of convection. With convection, rise times are drastically reduced (from years to a few months), loops are able to emerge at low latitudes, and tilt angles of emerging loops are consistent with Joy's Law for initial field strengths of greater than or equal to 40 kG. We also examine other asymmetries that develop between the leading and following legs of the emerging loops. Taking all the results together, we find that mid-range field strengths of approximately 40 - 50 kG produce emerging loops that best match the observed properties of solar active regions.

preprint2010arXiv

Global magnetic cycles in rapidly rotating younger suns

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.

preprint2010arXiv

Magnetic Cycles and Meridional Circulation in Global Models of Solar Convection

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.

preprint2010arXiv

Persistent Magnetic Wreaths in a Rapidly Rotating Sun

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.

preprint2006arXiv

Dynamo Action in the Solar Convection Zone and Tachocline: Pumping and Organization of Toroidal Fields

We present the first results from three-dimensional spherical shell simulations of magnetic dynamo action realized by turbulent convection penetrating downward into a tachocline of rotational shear. This permits us to assess several dynamical elements believed to be crucial to the operation of the solar global dynamo, variously involving differential rotation resulting from convection, magnetic pumping, and amplification of fields by stretching within the tachocline. The simulations reveal that strong axisymmetric toroidal magnetic fields (about 3000 G in strength) are realized within the lower stable layer, unlike in the convection zone where fluctuating fields are predominant. The toroidal fields in the stable layer possess a striking persistent antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. The associated mean poloidal magnetic fields there have a clear dipolar geometry, but we have not yet observed any distinctive reversals or latitudinal propagation. The presence of these deep magnetic fields appears to stabilize the sense of mean fields produced by vigorous dynamo action in the bulk of the convection zone.