Astrophysical simulations of convection frequently impose different thermal boundary conditions at the top and the bottom of the domain in an effort to more accurately model natural systems. In this work, we study Rayleigh-Benard convection (RBC) under the Boussinesq approximation. We examine simulations with mixed temperature boundary conditions in which the flux is fixed at the bottom boundary and the temperature is fixed at the top ("FT"). We aim to understand how FT boundaries change the nature of the convective solution compared to the traditional choice of thermal boundaries, in which the temperature is fixed at the top and bottom of the domain ("TT"). We demonstrate that the timescale of thermal relaxation for FT simulations is dependent upon the initial conditions. "Classic" initial conditions which employ a hydrostatically -- and thermally -- balanced linear temperature profile exhibit a long thermal relaxation. This long relaxation is not seen in FT simulations which use a TT simulation's nonlinear state as initial conditions ("TT-to-FT"). In the thermally relaxed, statistically stationary state, the mean behavior of an FT simulation corresponds to an equivalent simulation with TT boundaries, and time- and volume-averaged flow statistics like the Nusselt number and the Peclet number are indistinguishable between FT and TT simulations. FT boundaries are fundamentally asymmetric, and we examine the asymmetries that these boundaries produce in the flow. We find that the fixed-flux boundary produces more extreme temperature events than the fixed-temperature boundary. However, these near-boundary asymmetries do not measurably break the symmetry in the convective interior. We briefly explore rotating RBC to demonstrate that our findings with respect to thermal relaxation carry over to this more complex case, and to show the power of TT-to-FT initial conditions.

Large-scale convective flows called giant cells were once thought to transport the Sun's luminosity in the solar convection zone, but recent observations have called their existence into question. In place of large-scale flows, some authors have suggested the solar luminosity may instead be transported by small droplets of rapidly falling, low entropy fluid. This "entropy rain" could propagate as dense vortex rings, analogous to rising buoyant thermals in the Earth's atmosphere. In this work, we develop an analytical theory describing the evolution of dense, negatively buoyant thermals. We verify the theory with 2D cylindrical and 3D cartesian simulations of laminar, axisymmetric thermals in highly stratified atmospheres. Our results show that dense thermals fall in two categories: a stalling regime in which the droplets slow down and expand, and a falling regime in which the droplets accelerate and shrink as they propagate downwards. We estimate that solar downflows are in the falling regime and maintain their entropy perturbation against diffusion until they reach the base of the convection zone. This suggests that entropy rain may be an effective nonlocal mechanism for transporting the solar luminosity.

The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of common inputs to numerical models. The Predictive Rossby number can be specified independent of Rayleigh number, allowing suites of numerical solutions to separate the degree of rotational constraint from the strength of the driving of convection. We examine the scaling of convective transport in terms of the Nusselt number and the degree of turbulence in terms of the Reynolds number of the flow, and we find scaling laws nearly identical to those in nonrotational convection at low Rossby number when the Predictive Rossby number is held constant. Finally, we describe the boundary layers as a function of increasing turbulence at constant Rossby number.

High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in a simulation and using its characteristics to adjust the mean thermodynamic profile within the domain towards its evolved state. We study Rayleigh-B\'enard convection as a test case for AE. Evolved flow properties of AE solutions are measured to be within a few percent of solutions which are reached through standard evolution (SE) over a full thermal timescale. At the highest values of the Rayleigh number at which we compare SE and AE, we find that AE solutions require roughly an order of magnitude fewer computing hours to evolve than SE solutions.

We study fully compressible convection in the context of plane-parallel, polytropically stratified atmospheres. We perform a suite of 2D and 3D simulations in which we vary the initial superadiabaticity ($\epsilon$) and the Rayleigh number (Ra) while fixing the initial density stratification, aspect ratio, and Prandtl number. The evolved heat transport, quantified by the Nusselt number (Nu), follows scaling relationships similar to those found in the well-studied, incompressible Rayleigh-B\'{e}nard problem. This scaling holds up in both 2D and 3D and is not appreciably affected by the magnitude of $\epsilon$.

The two interferometers of the Laser Interferometry Gravitaional-wave Observatory (LIGO) recently detected gravitational waves from the mergers of binary black hole systems. Accurate calibration of the output of these detectors was crucial for the observation of these events, and the extraction of parameters of the sources. The principal tools used to calibrate the responses of the second-generation (Advanced) LIGO detectors to gravitational waves are systems based on radiation pressure and referred to as Photon Calibrators. These systems, which were completely redesigned for Advanced LIGO, include several significant upgrades that enable them to meet the calibration requirements of second-generation gravitational wave detectors in the new era of gravitational-wave astronomy. We report on the design, implementation, and operation of these Advanced LIGO Photon Calibrators that are currently providing fiducial displacements on the order of $10^{-18}$ m/$\sqrt{\textrm{Hz}}$ with accuracy and precision of better than 1 %.