We consider convection in an internally heated layer of fluid that is bounded below by a perfect insulator and above by a poor conductor. The poorly conducting boundary is modelled by a fixed heat flux. Using solely analytical methods, we find linear and energy stability thresholds for the static state, and we construct a lower bound on the mean temperature that applies to all flows. The linear stability analysis yields a Rayleigh number above which the static state is linearly unstable ($R_L$), and the energy analysis yields a Rayleigh number below which it is globally stable ($R_E$). For various boundary conditions on the velocity, exact expressions for $R_L$ and $R_E$ are found using long-wavelength asymptotics. Each $R_E$ is strictly smaller than the corresponding $R_L$ but is within 1%. The lower bound on the mean temperature is proven for no-slip velocity boundary conditions using the background method. The bound guarantees that the mean temperature of the fluid, relative to that of the top boundary, grows with the heating rate ($H$) no slower than $H^{2/3}$.

This work reviews basic features of both Rayleigh-B\'enard (RB) convection and internally heated (IH) convection, along with findings on IH convection from laboratory experiments and numerical simulations. In the first chapter, six canonical models of convection are described: three configurations of IH convection driven by constant and uniform volumetric heating, and three configurations of RB convection driven by the boundary conditions. The IH models are distinguished by differing pairs of thermal boundary conditions: top and bottom boundaries of equal temperature, an insulating bottom with heat flux fixed at the top, and an insulating bottom with temperature fixed at the top. The RB models also are distinguished by whether temperatures or heat fluxes are fixed at the top and bottom boundaries. Integral quantities important to heat transport are discussed, including the mean fluid temperature, the mean temperature difference between the boundaries, and the mean convective heat transport. Integral relations and bounds are presented, and further bounds are conjectured for the IH cases. The second chapter presents results that can be derived mathematically from the governing equations: linear and nonlinear stability thresholds of static states, and parameter-dependent bounds. Known bounds are on mean temperatures in IH convection and on convective transport in RB convection. The third chapter reviews numerical simulations and laboratory experiments on IH convection, emphasizing quantitative results.

We simulate three-dimensional, horizontally periodic Rayleigh-B\'enard convection between free-slip horizontal plates, rotating about a distant horizontal axis. When both the temperature difference between the plates and the rotation rate are sufficiently large, a strong horizontal wind is generated that is perpendicular to both the rotation vector and the gravity vector. The wind is turbulent, large-scale, and vertically sheared. Horizontal anisotropy, engendered here by rotation, appears necessary for such wind generation. Most of the kinetic energy of the flow resides in the wind, and the vertical turbulent heat flux is much lower on average than when there is no wind.

We report on direct numerical simulations of two-dimensional, horizontally periodic Rayleigh-B\'enard convection, focusing on its ability to drive large-scale horizontal flow that is vertically sheared. For the Prandtl numbers ($Pr$) between 1 and 10 simulated here, this large-scale shear can be induced by raising the Rayleigh number ($Ra$) sufficiently, and we explore the resulting convection for $Ra$ up to $10^{10}$. When present in our simulations, the sheared mean flow accounts for a large fraction of the total kinetic energy, and this fraction tends towards unity as $Ra\to\infty$. The shear helps disperse convective structures, and it reduces vertical heat flux; in parameter regimes where one state with large-scale shear and one without are both stable, the Nusselt number of the state with shear is smaller and grows more slowly with $Ra$. When the large-scale shear is present with $Pr\lesssim2$, the convection undergoes strong global oscillations on long timescales, and heat transport occurs in bursts. Nusselt numbers, time-averaged over these bursts, vary non-monotonically with $Ra$ for $Pr=1$. When the shear is present with $Pr\gtrsim3$, the flow does not burst, and convective heat transport is sustained at all times. Nusselt numbers then grow roughly as powers of $Ra$, but the growth rates are slower than any previously reported for Rayleigh-B\'enard convection without large-scale shear. We find the Nusselt numbers grow proportionally to $Ra^{0.077}$ when $Pr=3$ and to $Ra^{0.19}$ when $Pr=10$. Analogies with tokamak plasmas are described.