The validity of the anelastic approximation has recently been questioned in the regime of rapidly-rotating compressible convection in low Prandtl number fluids (Calkins et al. 2015). Given the broad usage and the high computational efficiency of sound-proof approaches in this astrophysically relevant regime, this paper clarifies the conditions for a safe application. The potential of the alternative pseudo-incompressible approximation is investigated, which in contrast to the anelastic approximation is shown to never break down for predicting the point of marginal stability. Its accuracy, however, decreases close to the parameters corresponding to the failure of the anelastic approach, which is shown to occur when the sound-crossing time of the domain exceeds a rotation time scale, i.e. for rotational Mach numbers greater than one. Concerning the supercritical case, which is naturally characterised by smaller rotational Mach numbers, we find that the anelastic approximation does not show unphysical behaviour. Growth rates computed with the linearised anelastic equations converge toward the corresponding fully compressible values as the Rayleigh number increases. Likewise, our fully nonlinear turbulent simulations, produced with our fully compressible and anelastic models and carried out in a highly supercritical, rotating, compressible, low Prandtl number regime show good agreement. However, this nonlinear test example is for only a moderately low convective Rossby number of 0.14.

The zonal winds on the surfaces of giant planets vary with latitude. Jupiter and Saturn, for example, have several bands of alternating eastward (prograde) and westward (retrograde) jets relative to the angular velocity of their global magnetic fields. These surface wind profiles are likely manifestations of the variations in depth and latitude of angular velocity deep within the liquid interiors of these planets. Two decades ago it was proposed that this differential rotation could be maintained by vortex stretching of convective fluid columns that span the interiors of these planets from the northern hemisphere surface to the southern hemisphere surface. This now classic mechanism explains the differential rotation seen in laboratory experiments and in computer simulations of, at best, weakly turbulent convection in rotating constant-density fluid spheres. However, these experiments and simulations are poor approximations for the density-stratified strongly-turbulent interiors of giant planets. The long thin global convective columns predicted by the classic geostrophic theory for these planets would likely not develop. Here we propose a much more robust mechanism for maintaining differential rotation in radius based on the local generation of vorticity as rising plumes expand and sinking plumes contract. Our high-resolution two-dimensional computer simulations demonstrate how this mechanism could maintain either prograde or retrograde surface winds in the equatorial region of a giant planet depending on how the density scale height varies with depth.

We present results of nonlinear numerical simulations of gravity wave driven shear flow oscillations in the equatorial plane of the solar radiative interior. These results show that many of the assumptions of quasi-linear theory are not valid. When only two waves are forced (prograde and retrograde) oscillatory mean flow is maintained; but critical layers often form and are dynamically important. When a spectrum of waves is forced, the non-linear wave-wave interactions are dynamically important, often acting to decrease the maintenance of a mean flow. The (in)coherence of such wave-wave interactions must be taken into account when describing wave driven mean flows.

Liu et al. 2008 discuss an important consideration for models of zonal winds deep within giant planets. However, the constraints they propose for the depth of the winds are based on their prescriptions for the internal structures of the magnetic field and zonal winds. The same kinematic analysis applied to other plausible configurations would produce no constraint on the depth to which the winds extend.

We discuss how rotation and binary interactions may be related to the diversity of type Ibc supernovae and long gamma-ray bursts. After presenting recent evolutionary models of massive single and binary stars including rotation, the Tayler-Spruit dynamo and binary interactions, we argue that the nature of SNe Ibc progenitors from binary systems may not significantly differ from that of single star progenitors in terms of rotation, and that most long GRB progenitors may be produced via the quasi-chemically homogeneous evolution at sub-solar metallicity. We also briefly discuss the possible role of magnetic fields generated in the convective core of a massive star for the transport of angular momentum, which is potentially important for future stellar evolution models of supernova and GRB progenitors.

We study how the pattern of thermal convection and differential rotation in the interior of a giant gaseous planet is affected by the presence of a small solid core as a function of the planetary rotation rate. We show, using 2D anelastic, hydrodynamic simulations, that the presence of a small solid core results in significantly different flow structure relative to that of a fully convective interior only if there is little or no planetary rotation.

We present numerical simulations of convective overshoot in a two-dimensional model of the solar equatorial plane. The simulated domain extends from 0.001 R_sun to 0.93 R_sun, spanning both convective and radiative regions. We show that convective penetration leads to a slightly extended, mildly subadiabatic temperature gradient beneath the convection zone below which there is a rapid transition to a strongly subadiabatic region. A slightly higher temperature is maintained in the overshoot region by adiabatic heating from overshooting plumes. This enhanced temperature may partially account for the sound speed discrepancy between the standard solar model and helioseismology. Simulations conducted with tracer particles suggest that a fully mixed region exists down to at least 0.687 R_sun.

We present self-consistent numerical simulations of the sun's convection zone and radiative interior using a two-dimensional model of its equatorial plane. The background reference state is a one-dimensional solar structure model. Turbulent convection in the outer convection zone continually excites gravity waves which propagate throughout the stable radiative interior and deposit their angular momentum. We find that angular velocity variations in the tachocline are driven by angular momentum transported by overshooting convective plumes rather than the nonlinear interaction of waves. The mean flow in the tachocline is time dependent but not oscillatory in direction. Since the forcing in this shallow region can not be described by simple linear waves, it is unlikely that the interaction of such waves is responsible for the solar cycle or the 1.3 year oscillation. However, in the deep radiative interior, the interaction of low amplitude gravity waves, continually excited by the overshooting plumes, is responsible for the angular velocity deviations observed there. Near the center of the model sun the angular velocity deviation is about two orders of magnitude greater than that in the bulk of the radiative region and reverses its direction (prograde to retrograde or vice versa) in the opposite sense of the angular velocity deviations that occur in the tachocline. Our simulations thus demonstrate how angular velocity variations in the solar core are linked to those in the tachocline, which themselves are driven by convective overshooting.

The thermonuclear runaway that culminates in the explosion of a Chandrasekhar mass white dwarf as a Type Ia supernova begins centuries before the star actually explodes. Here, using a 3D anelastic code, we examine numerically the convective flow during the last minute of that runaway, a time that is crucial in determining just where and how often the supernova ignites. We find that the overall convective flow is dipolar, with the higher temperature fluctuations in an outbound flow preferentially on one side of the star. Taken at face value, this suggests an asymmetric ignition that may well persist in the geometry of the final explosion. However, we also find that even a moderate amount of rotation tends to fracture this dipole flow, making ignition over a broader region more likely. Though our calculations lack the resolution to study the flow at astrophysically relevant Rayleigh numbers, we also speculate that the observed dipolar flow will become less organized as the viscosity becomes very small. Motion within the dipole flow shows evidence of turbulence, suggesting that only geometrically large fluctuations (~1 km) will persist to ignite the runaway. We also examine the probability density function for the temperature fluctuations, finding evidence for a Gaussian, rather than exponential distribution, which suggests that ignition sparks may be strongly spatially clustered.

We present numerical simulations of penetrative convection and gravity wave excitation in the Sun. Gravity waves are self-consistently generated by a convective zone overlying a radiative interior. We produce power spectra for gravity waves in the radiative region as well as estimates for the energy flux of gravity waves below the convection zone. We calculate a peak energy flux in waves below the convection zone to be three orders of magnitude smaller than previous estimates for m=1. The simulations show that the linear dispersion relation is a good approximation only deep below the convective-radiative boundary. Both low frequency propagating gravity waves as well as higher frequency standing modes are generated; although we find that convection does not continually drive the standing g-mode frequencies.

After briefly describing the anelastic approximation and Glatzmaier's code, we present results from our preliminary studies of core convection during the hydrogen burning phase of a 15 solar mass star, as well as our most recent results concerning convection in the oxygen shell of a 25 solar mass star.