Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the resulting differential rotation is markedly different between the two sets. The fixed-entropy-gradient simulations have strong differential rotation contrast and isocontours tilted along radial lines (in good agreement with the Sun's interior rotation revealed by helioseismology), whereas the fixed-entropy simulations have weaker contrast and contours tilted in the opposite sense. We examine in detail the force balances in our models and find that the poleward transport of heat by Busse columns drives a thermal wind responsible for the different rotation profiles. We conclude that the Sun's strong differential rotation along radial lines may result from the solar emissivity being invariant with latitude (which is similar to the fixed-entropy-gradient condition in our models) and the poleward transport of heat by Busse columns. In future work on convection in the solar context, we strongly advise modelers to use a fixed-gradient outer boundary condition.

In recent years, observers have found that the fraction of M-stars demonstrating significant magnetic activity transitions sharply from roughly $10\%$ for main-sequence stars earlier (more massive) than spectral type M3.5 (0.35 M$_\odot$) to nearly $90\%$ for stars later than M3.5. The latter are typically rotating quite rapidly, suggesting differing spin-down histories. Tantalizingly, it is also later than M3.5 at which main-sequence stars become fully convective, and may no longer contain a tachocline, a layer of rotational shear revealed by helioseismology to separate the radiative zone (RZ) and convection zone (CZ). We turn here to the more massive M-stars to study the impact such a layer may have on their internal dynamics. Using the spherical 3D MHD simulation code Rayleigh, we compare the properties of convective dynamos generated within quickly rotating (1, 2, and 4 $\Omega_\odot$) M2 (0.4 M$_\odot$) stars, with the computational domain either terminating at the base of the convection zone or permitting overshoot into the underlying stable region. We find that a tachocline is not necessary for the organization of strong toroidal wreaths of magnetism in these stars, though its presence can increase the coupling of mean field amplitudes to the stellar rotation rate. Additionally, in stars that undergo periodic cycles, we find that the presence of a tachocline tends to make the cycles both longer and more regular than they would have otherwise been. Finally, we find that the tachocline helps to enhance the surface poloidal fields and organize them into larger spatial scales, both of which provide favorable conditions for more rapid angular momentum loss through a magnetized stellar wind.

The calling card of solar magnetism is the sunspot cycle, during which sunspots regularly reverse their polarity sense every 11 years. However, a number of more complicated time-dependent behaviors have also been identified. In particular, there are temporal modulations associated with active longitudes and hemispheric asymmetry, when sunspots appear at certain solar longitudes or else in one hemisphere preferentially. So far, a direct link between between this asymmetric temporal behavior and the underlying solar dynamo has remained elusive. In this work, we present results from global, 3D magnetohydrodynamic (MHD) simulations, which for the first time display both behavior reminiscent of the sunspot cycle (regular polarity reversals and equatorward migration of internal magnetic field) and asymmetric, irregular behavior that in the simulations we interpret as active longitudes and hemispheric asymmetry. The simulations are thus bistable, in that the turbulent convection can stably support two distinct flavors of magnetism at different times, in superposition, or with smooth transitions from one state to the other. We discuss this new family of dynamo models in the context of the extensive observations of the Sun's surface magnetic field with the Solar and Heliospheric Observatory (SOHO) and the Solar Dynamics Observatory (SDO), as well as earlier observations of sunspot number and synoptic maps. We suggest that the solar dynamo itself may be bistable in nature, exhibiting two types of temporal behavior in the magnetic field.

Rotational scaling relationships are examined for the degree of equipartition between magnetic and kinetic energies in stellar convection zones. These scaling relationships are approached from two paradigms, with first a glance at scaling relationship built upon an energy-balance argument and second a look at a force-based scaling. The latter implies a transition between a nearly-constant inertial scaling when in the asymptotic limit of minimal diffusion and magnetostrophy, whereas the former implies a weaker scaling with convective Rossby number. Both scaling relationships are then compared to a suite of 3D convective dynamo simulations with a wide variety of domain geometries, stratifications, and range of convective Rossby numbers.

The advent of helioseismology has determined in detail the average rotation rate of the Sun as a function of radius and latitude. These data immediately reveal two striking boundary layers of shear in the solar convection zone (CZ): a tachocline at the base, where the differential rotation of the CZ transitions to solid-body rotation in the radiative zone, and a 35-Mm-thick near-surface shear layer (NSSL) at the top, where the rotation rate slows by about 5% with increasing radius. Though asteroseismology cannot probe the differential rotation of distant stars to the same level of detail that helioseismology can achieve for the Sun, it is possible that many cool stars with outer convective envelopes possess similar differential rotation characteristics, including both a tachocline and a NSSL. Here we present the results of 3D global hydrodynamic simulations of spherical-shell convection for a Sun-like star at different levels of density contrast across the shell. The simulations with high stratification possess characteristics of near-surface shear, especially at low latitudes. We discuss in detail the dynamical balance of torques giving rise to the NSSL in our models and interpret what these balances imply for the real Sun. We further discuss the dynamical causes that may serve to wipe out near-surface shear at high latitudes, and conclude by offering some theories as to how this problem might be tackled in future work.

The dynamical origins of the Sun's tachocline and near-surface shear layer (NSSL) are still not well understood. We have attempted to self-consistently reproduce a NSSL in numerical simulations of a solar-like convection zone by increasing the density contrast across rotating, 3D spherical shells. We explore the hypothesis that high density contrast leads to near-surface shear by creating a rotationally unconstrained layer of fast flows near the outer surface. Although our high-contrast models do have near-surface shear, it is confined primarily to low latitudes (between $\pm15^\circ$). Two distinct types of flow structures maintain the shear dynamically: rotationally $\textit{constrained}$ Busse columns aligned with the rotation axis and fast, rotationally $\textit{unconstrained}$ downflow plumes that deplete angular momentum from the outer fluid layers. The plumes form at all latitudes, and in fact are more efficient at transporting angular momentum inward at high latitudes. The presence of Busse columns at low latitudes thus appears essential to creating near-surface shear in our models. We conclude that a solar-like NSSL is unobtainable from a rotationally unconstrained outer fluid layer alone. In numerical models, the shear is eliminated through the advection of angular momentum by the meridional circulation. Therefore, a detailed understanding how the solar meridional circulation is dynamically achieved will be necessary to elucidate the origin of the Sun's NSSL.

M-type stars are quickly stepping into the forefront as some of the best candidates in searches for habitable Earth-like exoplanets, and yet many M-dwarfs exhibit extraordinary flaring events which would bombard otherwise habitable planets with ionizing radiation. In recent years, observers have found that the fraction of M-stars demonstrating significant magnetic activity transitions sharply from roughly $10\%$ for main-sequence stars earlier (more massive) than spectral type M3.5 (0.35 M$_\odot$) to nearly $90\%$ for stars later than M3.5. Suggestively, it is also later than M3.5 at which main-sequence stars become fully convective, and may no longer contain a tachocline. Using the spherical 3D MHD simulation code Rayleigh, we compare the peak field strengths, topologies, and time dependencies of convective dynamos generated within a quickly rotating (2 $\Omega_\odot$) M2 (0.4 M$_\odot$) star, with the computational domain either terminating at the base of the convection zone or including the tachocline. We find that while both models generate strong ($\sim$10kG), wreathlike toroidal fields exhibiting polarity reversals, the tachocline model provided a further reservoir for the toroidal field, which slowed the average reversal period from 100 rotations to more than 220 rotations and increased the spectral power of the low-order modes of the near-surface radial field by a factor of 4.

Global 3D simulations of solar giant-cell convection have provided significant insight into the processes which yield the Sun's observed differential rotation and cyclic dynamo action. However, as we move to higher resolution simulations a variety of codes have encountered what has been termed the convection conundrum. As these simulations increase in resolution and hence the level of turbulence achieved, they tend to produce weak or even anti-solar differential rotation patterns associated with a weak rotational influence (high Rossby number) due to large convective velocities. One potential culprit for this convection conundrum is the upper boundary condition applied in most simulations which is generally impenetrable. Here we present an alternative stochastic plume boundary condition which imposes small-scale convective plumes designed to mimic near-surface convective downflows, thus allowing convection to carry the majority of the outward solar energy flux up to and through our simulated upper boundary. The use of a plume boundary condition leads to significant changes in the convective driving realized in the simulated domain and thus to the convective energy transport, the dominant scale of the convective enthalpy flux, and the relative strength of the strongest downflows, the downflow network, and the convective upflows. These changes are present even far from the upper boundary layer. Additionally, we demonstrate that in spite of significant changes, giant cell morphology in the convective patterns is still achieved with self-organization of the imposed boundary plumes into downflow lanes, cellular patterns, and even rotationally-aligned banana cells in equatorial regions. This plume boundary presents an alternative pathway for 3D global convection simulations where driving is non-local and may provide a new approach towards addressing the convection conundrum.

We seek to characterize how the change of global rotation rate influences the overall dynamics and large scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the latter profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number $R_{of} = \tilde{\omega}/2\Omega_* \sim \tilde{v}/2\Omega_* R_*$, which we calibrate based on our 3-D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map ($R_{of} \gtrsim 1.5-2$) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

This paper provides a brief look at dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Two simple models are examined, where one that assumes magnetostrophy and another that includes the effects of inertia. These models are then compared to a suite of convective dynamo simulations of the convective core of a main-sequence B-type star and applied to its later evolutionary stages.

The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through 3D global simulations of convective core dynamos operating within a young $\Msun{10}$ B-type star, using the anelastic spherical harmonic (ASH) code. These simulations capture the inner 65\% of this star by radius, encompassing the convective nuclear-burning core (about 23\% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05\% to 16\% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces.

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.

The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulations relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the $\alpha$-$\Omega$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified.

Using a new implementation of ring-diagram helioseismology, we ascertain the strength and spatial scale of convective flows throughout the near-surface shear layer. Our ring-diagram technique employs highly overlapped analysis regions and an efficient method of 3D inversion to measure convective motions with a resolution that ranges from $3 \ \mathrm{Mm}$ at the surface to $80 \ \mathrm{Mm}$ at the base of the layer. We find the rms horizontal flow speed to peak at $427 \ \mathrm{m \ s^{-1}}$ at the photosphere and fall to a minimum of $124 \ \mathrm{m \ s^{-1}}$ between $20 \ \mathrm{Mm}$ and $30 \ \mathrm{Mm}$. From the velocity amplitude and the dominant horizontal scales seen at each depth, we infer the level of rotational influence on convection to be low near the surface, but transition to a significant level at the base of the near-surface shear layer with a Rossby number varying between 2.2 to as low as 0.1.

The global-scale dynamo action achieved in a simulation of a Sun-like star rotating at thrice the solar rate is assessed. The 3-D MHD Anelastic Spherical Harmonic (ASH) code, augmented with a viscosity minimization scheme, is employed to capture convection and dynamo processes in this G-type star. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. It is found that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals occurring roughly every 6.2 years. The magnetic energy also rises and falls with a regular period. The magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the processes leading to polarity reversals are more complex, appearing to arise from the interaction of convection with the mean toroidal fields. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation, and potentially to a nonlinear dynamo wave. This simulation also enters a grand minimum lasting roughly 20~years, after which the dynamo recovers its regular polarity cycles.

Inferences of sub-surface flow velocities using local domain ring-diagram helioseismology depend on measuring the frequency splittings of oscillation modes seen in acoustic power spectra. Current methods for making these measurements utilize maximum-likelihood fitting techniques to match a model of modal power to the spectra. The model typically describes a single oscillation mode, and each mode in a given power spectrum is fit independently. We present a new method that produces measurements with greater reliability and accuracy by fitting multiple modes simultaneously. We demonstrate how this method permits measurements of sub-surface flows deeper into the Sun while providing higher uniformity in data coverage and velocity response closer to the limb of the solar disk. While the previous fitting method performs better for some measurements of low-phase-speed modes, we find this new method to be particularly useful for high phase-speed modes and small spatial areas.

We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.

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.

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$\Omega_\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.

We examine the role of small-scale granulation in helping to drive supergranulation and even larger scales of convection. The granulation is modeled as localized cooling events introduced at the upper boundary of a 3-D simulation of compressible convection in a rotating spherical shell segment. With a sufficient number of stochastic cooling events compared to uniform cooling, we find that supergranular scales are realized, along with a differential rotation that becomes increasingly solar-like.

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.

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 $\Omega$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $\alpha$-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]

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.

As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance of the flows established in rotating spherical shell segments that capture much of the NSSL with the Curved Spherical Segment (CSS) code using two different diffusion schemes. The CSS code is a new massively parallel modeling tool capable of simulating 3-D compressible MHD convection with a realistic solar stratification in rotating spherical shell segments.

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 $\alpha$-effect explains the regeneration of poloidal field.

Main-sequence massive stars possess convective cores that likely harbor strong dynamo action. To assess the role of core convection in building magnetic fields within these stars, we employ the 3-D anelastic spherical harmonic (ASH) code to model turbulent dynamics within a 10 solar mass main-sequence (MS) B-type star rotating at 4 times the solar rate. We find that strong (900 kG) magnetic fields arise within the turbulence of the core and penetrate into the stably stratified radiative zone. These fields exhibit complex, time-dependent behavior including reversals in magnetic polarity and shifts between which hemisphere dominates the total magnetic energy.

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.

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.

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.

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.

Using high-resolution ring analysis (HRRA) we deduce subsurface flows within magnetic active regions and within quiet sun. With this procedure we are capable of measuring flows with a horizontal spatial resolution of 2 degrees in heliographic angle (or roughly 20 Mm). From the resulting flow fields we deduce mean inflow rates into active regions, mean circulation speeds around active regions, and probability density functions (PDFs) of properties of the flow field. These analyses indicate that active regions have a zonal velocity that exceeds that of quiet sun at the same latitude by 20 m/s, yet active regions advect poleward at the same rate as quiet sun. We also find that almost all active regions possess a mean inflow (20-30 m/s) and a cyclonic circulation (~5 m/s) at their peripheries, whereas their cores, where the sunspots are located, are zones of strong anticyclonic outflow (~50 m/s). From the PDFs, we find that active regions modify the structure of convection with a scale greater than that of supergranulation. Instead of possessing an asymmetry between inflows and outflows (with a larger percentage of the surface occupied by outflows), as is seen in quiet sun, active regions possess symmetric distributions.

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.

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.

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.

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.

The global scales of solar convection are studied through three-dimensional simulations of compressible convection carried out in spherical shells of rotating fluid which extend from the base of the convection zone to within 15 Mm of the photosphere. Such modelling at the highest spatial resolution to date allows study of distinctly turbulent convection, revealing that coherent downflow structures associated with giant cells continue to play a significant role in maintaining the strong differential rotation that is achieved. These giant cells at lower latitudes exhibit prograde propagation relative to the mean zonal flow, or differential rotation, that they establish, and retrograde propagation of more isotropic structures with vortical character at mid and high latitudes. The interstices of the downflow networks often possess strong and compact cyclonic flows. The evolving giant-cell downflow systems can be partly masked by the intense smaller scales of convection driven closer to the surface, yet they are likely to be detectable with the helioseismic probing that is now becoming available. Indeed, the meandering streams and varying cellular subsurface flows revealed by helioseismology must be sampling contributions from the giant cells, yet it is difficult to separate out these signals from those attributed to the faster horizontal flows of supergranulation. To aid in such detection, we use our simulations to describe how the properties of giant cells may be expected to vary with depth, how their patterns evolve in time, and analyze the statistical features of correlations within these complex flow fields.

Core convection and dynamo activity deep within rotating A-type stars of 2 solar masses are studied with 3--D nonlinear simulations. Our modeling considers the inner 30% by radius of such stars, thus capturing within a spherical domain the convective core and a modest portion of the surrounding radiative envelope. The MHD equations are solved using the ASH code to examine turbulent flows and magnetic fields, both of which exhibit intricate time dependence. By introducing small seed magnetic fields into our progenitor hydrodynamic models rotating at one and four times the solar rate, we assess here how the vigorous convection can amplify those fields and sustain them against ohmic decay. Dynamo action is indeed realized, ultimately yielding magnetic fields that are in energy equipartion with the flow. Such magnetism reduces the differential rotation obtained in the progenitors, partly by Maxwell stresses that transport angular momentum poleward and oppose the Reynolds stresses in the latitudinal balance. In contrast, in the radial direction we find that the Maxwell and Reynolds stresses may act together to transport angular momentum. The central columns of slow rotation established in the progenitors are weakened, with the differential rotation waxing and waning in strength as the simulations evolve. We assess the morphology of the flows and magnetic fields, their complex temporal variations, and the manner in which dynamo action is sustained. Differential rotation and helical convection are both found to play roles in giving rise to the magnetic fields. The magnetism is dominated by strong fluctuating fields throughout the core, with the axisymmetric (mean) fields there relatively weak.

The operation of the solar global dynamo appears to involve many dynamical elements. Self-consistent MHD simulations which realistically incorporate all of these processes are not yet computationally feasible, though some elements can now be studied with reasonable fidelity. Here we consider the manner in which turbulent compressible convection within the bulk of the solar convection zone can generate large-scale magnetic fields through dynamo action. We accomplish this through a series of three-dimensional numerical simulations of MHD convection within rotating spherical shells using our ASH code on massively parallel supercomputers. Since differential rotation is a key ingredient in all dynamo models, we also examine here the nature of the rotation profiles that can be sustained within the deep convection zone as strong magnetic fields are built and maintained. We find that the convection is able to maintain a solar-like angular velocity profile despite the influence of Maxwell stresses which tend to oppose Reynolds stresses and thus reduce the latitudinal angular velocity contrast throughout the convection zone. The dynamo-generated magnetic fields exhibit a complex structure and evolution, with radial fields concentrated in downflow lanes and toroidal fields organized into twisted ribbons which are extended in longitude and which achieve field strengths of up to 5000 G. The flows and fields exhibit substantial kinetic and magnetic helicity although systematic hemispherical patterns are only apparent in the former. Fluctuating fields dominate the magnetic energy and account for most of the back-reaction on the flow via Lorentz forces. Mean fields are relatively weak and do not exhibit systematic latitudinal propagation or periodic polarity reversals as in the sun. This may be attributed to the absence of a tachocline.

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.

We present the results of 3--D nonlinear simulations of magnetic dynamo action by core convection within A-type stars of 2 solar masses, at a range of rotation rates. We consider the inner 30% by radius of such stars, with the spherical domain thereby encompassing the convective core and a portion of the surrounding radiative envelope. The compressible Navier-Stokes equations, subject to the anelastic approximation, are solved to examine highly nonlinear flows that span multiple scale heights, exhibit intricate time dependence, and admit magnetic dynamo action. Small initial seed magnetic fields are found to be amplified greatly by the convective and zonal flows. The central columns of strikingly slow rotation found in some of our progenitor hydrodynamic simulations continue to be realized in some simulations to a lesser degree, with such differential rotation arising from the redistribution of angular momentum by the nonlinear convection and magnetic fields. We assess the properties of the magnetic fields thus generated, the extent of convective penetration, the magnitude of the differential rotation, and the excitation of gravity waves within the radiative envelope.

We present the results of 3--D simulations of core convection within A-type stars of 2 solar masses, at a range of rotation rates. We consider the inner 30% by radius of such stars, thereby encompassing the convective core and some of the surrounding radiative envelope. We utilize our anelastic spherical harmonic (ASH) code, which solves the compressible Navier-Stokes equations in the anelastic approximation, to examine highly nonlinear flows that can span multiple scale heights. The cores of these stars are found to rotate differentially, with central cylindrical regions of strikingly slow rotation achieved in our simulations of stars whose convective Rossby number (R_{oc}) is less than unity. Such differential rotation results from the redistribution of angular momentum by the nonlinear convection that strongly senses the overall rotation of the star. Penetrative convective motions extend into the overlying radiative zone, yielding a prolate shape (aligned with the rotation axis) to the central region in which nearly adiabatic stratification is achieved. This is further surrounded by a region of overshooting motions, the extent of which is greater at the equator than at the poles, yielding an overall spherical shape to the domain experiencing at least some convective mixing. We assess the overshooting achieved as the stability of the radiative exterior is varied, and the weak circulations that result in that exterior. The convective plumes serve to excite gravity waves in the radiative envelope, ranging from localized ripples of many scales to some remarkable global resonances.

We discuss recent advances made in modelling the complex magnetohydrodynamics of the Sun using our anelastic spherical harmonics (ASH) code. We have conducted extensive 3--D simulations of compressible convection in rotating spherical shells with and without magnetic fields, to study the coupling between global-scale convection and rotation in seeking to understand how the solar differential rotation is established and maintained. Such simulations capable of studying fairly turbulent convection have been enabled by massively parallel supercomputers. The resulting convection within domains that capture a good fraction of the bulk of the convection zone is highly time dependent and intricate, and is dominated by intermittent upflows and networks of strong downflows. A high degree of coherent structures involving downflowing plumes can be embedded in otherwise chaotic flow fields. These vortical structures play a significant role in yielding Reynolds stresses that serve to redistribute angular momentum, leading to differential rotation profiles with pole to equator contrasts of about 30% in angular velocity $\Omega$ and some constancy along radial line at mid latitudes, thereby making good contact with deductions from helioseismology. When a magnetic field is introduced, a dynamo regime can be found that does not destroy the strong differential rotation achieved in pure hydrodynamics cases. The magnetic fields are found to concentrate around the downflowing networks and to have significant north-south asymmetry and helicity.

The advent of massively parallel supercomputing has begun to permit explicit 3--D simulations of turbulent convection occurring within the cores of early-type main sequence stars. Such studies should complement the stellar structure and evolution efforts that have so far largely employed 1--D nonlocal mixing length descriptions for the transport, mixing and overshooting achieved by core convection. We have turned to A-type stars as representative of many of the dynamical challenges raised by core convection within rotating stars. The differential rotation and meridional circulations achieved deep within the star by the convection, the likelihood of sustained magnetic dynamo action there, and the bringing of fresh fuel into the core by overshooting motions, thereby influencing main sequence lifetimes, all constitute interesting dynamical questions that require detailed modelling of global-scale convection. Using our anelastic spherical harmonic (ASH) code tested on the solar differential rotation problem, we have conducted a series of 3--D spherical domain simulations that deal with a simplified description of the central regions of rotating A-type stars, i.e a convectively unstable core is surrounded by a stable radiative envelope. A sequence of 3--D simulations are used to assess the properties of the convection (its global patterns, differential rotation, meridional circulations, extent and latitudinal variation of the overshooting) as transitions are made between laminar and turbulent states by changing the effective diffusivities, rotation rates, and subadiabaticity of the radiative exterior. We report on the properties deduced from these models for both the extent of penetration and the profile of rotation sustained by the convection.

The differential rotation of the sun, as deduced from helioseismology, exhibits a prominent radial shear layer near the top of the convection zone wherein negative radial gradients of angular velocity are evident in the low- and mid-latitude regions spanning the outer 5% of the solar radius. Supergranulation and related scales of turbulent convection are likely to play a significant role in the maintenance of such radial gradients, and may influence dynamics on a global scale in ways that are not yet understood. To investigate such dynamics, we have constructed a series of three-dimensional numerical simulations of turbulent compressible convection within spherical shells, dealing with shallow domains to make such modeling computationally tractable. These simulations are the first models of solar convection in a spherical geometry that can explicitly resolve both the largest dynamical scales of the system (of order the solar radius) as well as smaller-scale convective overturning motions comparable in size to solar supergranulation (20--40 Mm). We find that convection within these simulations spans a large range of horizontal scales, and that the radial angular velocity gradient in these models is typically negative, especially in low- and mid-latitude regions. Analyses of the angular momentum transport indicates that such gradients are maintained by Reynolds stresses associated with the convection, transporting angular momentum inward to balance the outward transport achieved by viscous diffusion and large-scale flows in the meridional plane. We suggest that similar mechanisms associated with smaller-scale convection in the sun may contribute to the maintenance of the observed radial shear layer located immediately below the solar photosphere.

The intense turbulence present in the solar convection zone is a major challenge to both theory and simulation as one tries to understand the origins of the striking differential rotation profile with radius and latitude that has been revealed by helioseismology. The differential rotation must be an essential element in the operation of the solar magnetic dynamo and its cycles of activity, yet there are many aspects of the interplay between convection, rotation and magnetic fields that are still unclear. We have here carried out a series of 3--D numerical simulations of turbulent convection within deep spherical shells using our anelastic spherical harmonic (ASH) code on massively parallel supercomputers. These studies of the global dynamics of the solar convection zone concentrate on how the differential rotation and meridional circulation are established. We have analyzed the transport of angular momentum in establishing such differential rotation, and clarified the roles played by Reynolds stresses and the meridional circulation in this process. We have found that the Reynolds stresses are crucial in transporting angular momentum toward the equator. The effects of baroclinicity (thermal wind) have been found to have a modest role in the resulting mean zonal flows. The simulations have produced differential rotation profiles within the bulk of the convection zone that make reasonable contact with ones inferred from helioseismic inversions, namely possessing a fast equator, an angular velocity difference of about 30% from equator to pole, and some constancy along radial lines at mid-latitudes.

Helioseismic observations have detected small temporal variations of the rotation rate below the solar surface corresponding to the so-called `torsional oscillations' known from Doppler measurements of the surface. These appear as bands of slower and faster than average rotation moving equatorward. Here we establish, using complementary helioseismic observations over four years from the GONG network and from the MDI instrument on board SOHO, that the banded flows are not merely a near-surface phenomenon: rather they extend downward at least 60 Mm (some 8% of the total solar radius) and thus are evident over a significant fraction of the nearly 200 Mm depth of the solar convection zone.