Search

The process referred to as "semi-convection" in astrophysics and "double-diffusive convection in the diffusive regime" in Earth and planetary sciences occurs in stellar and planetary interiors in regions which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. In this series of papers, we analyze the results of an extensive suite of three-dimensional (3D) numerical simulations of the process, and ultimately propose a new 1D prescription for heat and compositional transport in this regime which can be used in stellar or planetary structure and evolution models. In a preliminary study of the phenomenon, Rosenblum et al. showed that, after saturation of the primary instability, a system can evolve in one of two possible ways: the induced turbulence either remains homogeneous, with very weak transport properties, or transitions into a thermo-compositional staircase where the transport rate is much larger (albeit still smaller than in standard convection). In this paper, we show that this dichotomous behavior is a robust property of semi-convection across a wide region of parameter space. We propose a simple semi-analytical criterion to determine whether layer formation is expected or not, and at what rate it proceeds, as a function of the background stratification and of the diffusion parameters (viscosity, thermal diffusivity, and compositional diffusivity) only. The theoretical criterion matches the outcome of our numerical simulations very adequately in the computationally accessible "planetary" parameter regime and can be extrapolated to the stellar parameter regime. Subsequent papers will address more specifically the question of quantifying transport in the layered case and in the non-layered case.

Double-diffusive convection, often referred to as semi-convection in astrophysics, occurs in thermally and compositionally stratified systems which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. This process has been given relatively little attention so far, and its properties remain poorly constrained. In this paper, we present and analyze a set of three-dimensional simulations of this phenomenon in a Cartesian domain under the Boussinesq approximation. We find that in some cases the double-diffusive convection saturates into a state of homogeneous turbulence, but with turbulent fluxes several orders of magnitude smaller than those expected from direct overturning convection. In other cases, the system rapidly and spontaneously develops closely packed thermo-compositional layers, which later successively merge until a single layer is left. We compare the output of our simulations with an existing theory of layer formation in the oceanographic context and find very good agreement between the model and our results. The thermal and compositional mixing rates increase significantly during layer formation and increase even further with each merger. We find that the heat flux through the staircase is a simple function of the layer height. We conclude by proposing a new approach to studying transport by double-diffusive convection in astrophysics.

We present the first three-dimensional simulations of fingering convection performed at parameter values approaching those relevant for astrophysics. Our simulations reveal the existence of simple asymptotic scaling laws for turbulent heat and compositional transport, which can be straightforwardly extrapolated from our numerically tractable values to the true astrophysical regime. Our investigation also indicates that thermo-compositional "staircases," a key consequence of fingering convection in the ocean, cannot form spontaneously in the fingering regime in stellar interiors. Our proposed empirically determined transport laws thus provide simple prescriptions for mixing by fingering convection in a variety of astrophysical situations, and should, from here on, be used preferentially over older and less accurate parameterizations. They also establish that fingering convection does not provide sufficient extra-mixing to explain observed chemical abundances in red giant branch stars.

Regions of the ocean's thermocline unstable to salt fingering are often observed to host thermohaline staircases, stacks of deep well-mixed convective layers separated by thin stably-stratified interfaces. Decades after their discovery, however, their origin remains controversial. In this paper we use 3D direct numerical simulations to shed light on the problem. We study the evolution of an analogous double-diffusive system, starting from an initial statistically homogeneous fingering state and find that it spontaneously transforms into a layered state. By analysing our results in the light of the mean-field theory developed in Paper I, a clear picture of the sequence of events resulting in the staircase formation emerges. A collective instability of homogeneous fingering convection first excites a field of gravity waves, with a well-defined vertical wavelength. However, the waves saturate early through regular but localized breaking events, and are not directly responsible for the formation of the staircase. Meanwhile, slower-growing, horizontally invariant but vertically quasi-periodic gamma-modes are also excited and grow according to the gamma-instability mechanism. Our results suggest that the nonlinear interaction between these various mean-field modes of instability leads to the selection of one particular gamma-mode as the staircase progenitor. Upon reaching a critical amplitude, this progenitor overturns into a fully-formed staircase. We conclude by extending the results of our simulations to real oceanic parameter values, and find that the progenitor gamma-mode is expected to grow on a timescale of a few hours, and leads to the formation of a thermohaline staircase in about one day with an initial spacing of the order of one to two metres.

Double-diffusive instabilities are often invoked to explain enhanced transport in stably-stratified fluids. The most-studied natural manifestation of this process, fingering convection, commonly occurs in the ocean's thermocline and typically increases diapycnal mixing by two orders of magnitude over molecular diffusion. Fingering convection is also often associated with structures on much larger scales, such as thermohaline intrusions, gravity waves and thermohaline staircases. In this paper, we present an exhaustive study of the phenomenon from small to large scales. We perform the first three-dimensional simulations of the process at realistic values of the heat and salt diffusivities and provide accurate estimates of the induced turbulent transport. Our results are consistent with oceanic field measurements of diapycnal mixing in fingering regions. We then develop a generalized mean-field theory to study the stability of fingering systems to large-scale perturbations, using our calculated turbulent fluxes to parameterize small-scale transport. The theory recovers the intrusive instability, the collective instability, and the gamma-instability as limiting cases. We find that the fastest-growing large-scale mode depends sensitively on the ratio of the background gradients of temperature and salinity (the density ratio). While only intrusive modes exist at high density ratios, the collective and gamma-instabilities dominate the system at the low density ratios where staircases are typically observed. We conclude by discussing our findings in the context of staircase formation theory.