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.

Convection is the predominant mechanism by which energy and angular momentum are transported in the outer portion of the Sun. The resulting overturning motions are also the primary energy source for the solar magnetic field. An accurate solar dynamo model therefore requires a complete description of the convective motions, but these motions remain poorly understood. Studying stellar convection numerically remains challenging; it occurs within a parameter regime that is extreme by computational standards. The fluid properties of the convection zone are characterized in part by the Prandtl number $\mathrm{Pr}=\nu/\kappa$, where $\nu$ is the kinematic viscosity and $\kappa$ is the thermal diffusion; in stars, $\mathrm{Pr}$ is extremely low, $\mathrm{Pr}\approx 10^{-7}$. The influence of $\mathrm{Pr}$ on the convective motions at the heart of the dynamo is not well understood since most numerical studies are limited to using $\mathrm{Pr}\approx1$. We systematically vary $\mathrm{Pr}$ and the degree of thermal forcing, characterized through a Rayleigh number, to explore its influence on the convective dynamics. For sufficiently large thermal driving, the simulations reach a so-called convective free-fall state where diffusion no longer plays an important role in the interior dynamics. Simulations with a lower $\mathrm{Pr}$ generate faster convective flows and broader ranges of scales for equivalent levels of thermal forcing. Characteristics of the spectral distribution of the velocity remain largely insensitive to changes in $\mathrm{Pr}$. Importantly, we find that $\mathrm{Pr}$ plays a key role in determining when the free-fall regime is reached by controlling the thickness of the thermal boundary layer.

We investigate how rotationally-constrained, deep convection might give rise to supergranulation, the largest distinct spatial scale of convection observed in the solar photosphere. While supergranulation is only weakly influenced by rotation, larger spatial scales of convection sample the deep convection zone and are presumably rotationally influenced. We present numerical results from a series of nonlinear, 3-D simulations of rotating convection and examine the velocity power distribution realized under a range of Rossby numbers. When rotation is present, the convective power distribution possesses a pronounced peak, at characteristic wavenumber $\ell_\mathrm{peak}$, whose value increases as the Rossby number is decreased. This distribution of power contrasts with that realized in non-rotating convection, where power increases monotonically from high to low wavenumbers. We find that spatial scales smaller than $\ell_\mathrm{peak}$ behave in analogy to non-rotating convection. Spatial scales larger than $\ell_\mathrm{peak}$ are rotationally-constrained and possess substantially reduced power relative to the non-rotating system. We argue that the supergranular scale emerges due to a suppression of power on spatial scales larger than $\ell\approx100$ owing to the presence of deep, rotationally-constrained convection. Supergranulation thus represents the largest non-rotationally-constrained mode of solar convection. We conclude that the characteristic spatial scale of supergranulation bounds that of the deep convective motions from above, making supergranulation an indirect measure of the deep-seated dynamics at work in the solar dynamo. Using the spatial scale of supergranulation in conjunction with our numerical results, we estimate an upper bound of 10 m s$^{-1}$ for the Sun's bulk rms convective velocity.

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 $\kappa$ along two paths in parameter space. If the kinematic viscosity $\nu$ is decreased proportionally with $\kappa$ (fixing the Prandtl number $P_r = \nu/\kappa$), we find that U increases but asymptotes toward a constant value, as found by Featherstone & Hindman (2016). However, if $\nu$ is held fixed while decreasing $\kappa$ (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 $\nu$) 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.

Convection plays a central role in the dynamics of any stellar interior, and yet its operation remains largely-hidden from direct observation. As a result, much of our understanding concerning stellar convection necessarily derives from theoretical and computational models. The Sun is, however, exceptional in that regard. The wealth of observational data afforded by its proximity provides a unique testbed for comparing convection models against observations. When such comparisons are carried out, surprising inconsistencies between those models and observations become apparent. Both photospheric and helioseismic measurements suggest that convection simulations may overestimate convective flow speeds on large spatial scales. Moreover, many solar convection simulations have difficulty reproducing the observed solar differential rotation due to this apparent overestimation. We present a series of 3-dimensional (3-D) stellar convection simulations designed to examine how the amplitude and spectral distribution of convective flows are established within a star's interior. While these simulations are non-magnetic and non-rotating in nature, they demonstrate two robust phenomena. When run with sufficiently high Rayleigh number, the integrated kinetic energy of the convection becomes effectively independent of thermal diffusion, but the spectral distribution of that kinetic energy remains sensitive to both of these quantities. A simulation that has converged to a diffusion-independent value of kinetic energy will divide that energy between spatial scales such that low-wavenumber power is overestimated, and high-wavenumber power is underestimated relative to a comparable system possessing higher Rayleigh number. We discuss the implications of these results in light of the current inconsistencies between models and observations.

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.

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.

We describe a 3D finite-difference spherical anelastic MHD (FSAM) code for modeling the subsonic dynamic processes in the solar convective envelope. A comparison of this code with the widely used global spectral anlastic MHD code, ASH (Anelastic Spherical Harmonics), shows that FSAM produces convective flows with statistical properties and mean flows similar to the ASH results. Using FSAM, we first simulate the rotating solar convection in a partial spherical shell domain and obtain a statistically steady, giant-cell convective flow with a solar-like differential rotation. We then insert buoyant toroidal flux tubes near the bottom of the convecting envelope and simulate the rise of the flux tubes in the presence of the giant cell convection. We find that for buoyant flux tubes with an initial field strength of 100 kG, the magnetic buoyancy largely determines the rise of the tubes although strong down flows produce significant undulation and distortion to the shape of the emerging $\Omega$-shaped loops. The convective flows significantly reduce the rise time it takes for the apex of the flux tube to reach the top. For the weakly twisted and untwisted cases we simulated, the apex portion is found to rise nearly radially to the top in about a month, and produce an emerging region (at a depth of about 30 Mm below the photosphere) with an overall tilt angle consistent with the active region tilts, although the emergence pattern is more complex compared to the case without convection. Near the top boundary at a depth of about 30 Mm, the emerging flux shows a retrograde zonal flow of about 345 m/s relative to the mean flow at that latitude.

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.