Search for: All records

We study the dynamics of salt fingers in the regime of slow salinity diffusion (small inverse Lewis number) and strong stratification (large density ratio), focusing on regimes relevant to Earth’s oceans. Using three-dimensional direct numerical simulations in periodic domains, we show that salt fingers exhibit rich, multiscale dynamics in this regime, with vertically elongated fingers that are twisted into helical shapes at large scales by mean flows and disrupted at small scales by isotropic eddies. We use a multiscale asymptotic analysis to motivate a reduced set of partial differential equations that filters internal gravity waves and removes inertia from all parts of the momentum equation except for the Reynolds stress that drives the helical mean flow. When simulated numerically, the reduced equations capture the same dynamics and fluxes as the full equations in the appropriate regime. The reduced equations enforce zero helicity in all fluctuations about the mean flow, implying that the symmetry-breaking helical flow is generated spontaneously by strictly non-helical fluctuations. 

Abstract Polluted white dwarfs (WDs) with small surface convection zones deposit significant concentrations of heavy elements to the underlying radiative interior, presumably driving thermohaline convection. Current models of polluted WDs frequently fail to account for this effect, although its inclusion can increase the inferred accretion rate by orders of magnitude. It has been argued that this instability cannot be treated as a continuous mixing process and thus should not be considered in these models. In this work, we study 3D simulations of a thermohaline-unstable layer propagating into an underlying stable region, approximating the polluted WD scenario. We find that although thermohaline convection works to reduce driving gradients somewhat, the front continues to propagate and the system remains unstable. Importantly, the turbulent flux of metals broadly dominates over the diffusive flux in quantitative agreement with existing mixing prescriptions implemented in some stellar evolution models (except slightly below the boundary of the propagating front, where recent prescriptions neglect overshoot-like effects). Thus, our results broadly support polluted WD models that include thermohaline mixing in their estimates of the settling rate. 

Abstract Fingering convection (also known as thermohaline convection) is a process that drives the vertical transport of chemical elements in regions of stellar radiative zones where the mean molecular weight increases with radius. Recently, Harrington & Garaud used three-dimensional direct numerical simulations (DNS) to show that a vertical magnetic field can dramatically enhance the rate of chemical mixing by fingering convection. Furthermore, they proposed a so-called “parasitic saturation” theory to model this process. Here, we test their model over a broad range of parameter space, using a suite of DNS of magnetized fingering convection, varying the magnetic Prandtl number, magnetic field strength, and composition gradient. We find that the rate of chemical mixing measured in the simulations is not always predicted accurately by their existing model, in particular when the magnetic diffusivity is large. We then present an extension of the Harrington & Garaud model which resolves this issue. When applied to stellar parameters, it recovers the results of Harrington & Garaud except in the limit where fingering convection becomes marginally stable, where the new model is preferred. We discuss the implications of our findings for stellar structure and evolution. 

Abstract Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we propose a new framework to facilitate direct comparison: using carbon-to-nitrogen measurements from the Sloan Digital Sky Survey-IV APOGEE survey as a probe of mixing and a fluid parameter known as the reduced density ratio from one-dimensional stellar evolution programs, we compare the observed amount of extra mixing on the upper giant branch to predicted trends from three-dimensional fluid dynamics simulations. Using this method, we are able to empirically constrain how mixing efficiency should vary with the reduced density ratio. We find the observed amount of extra mixing is strongly correlated with the reduced density ratio and that trends between reduced density ratio and fundamental stellar parameters are robust across choices for modeling prescription. We show that stars with available mixing data tend to have relatively low density ratios, which should inform the regimes selected for future simulation efforts. Finally, we show that there is increased mixing at low reduced density ratios, which is consistent with current hydrodynamical models of thermohaline mixing. The introduction of this framework sets a new standard for theoretical modeling efforts, as validation for not only the amount of extra mixing, but trends between the degree of extra mixing and fundamental stellar parameters is now possible. 

Abstract Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions that are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild unstable but Ledoux stable due to a composition gradient. Over many convective overturn timescales, the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists.