Stephan Stellmach

Boundary layer control of rotating convection systems.

Turbulent rotating convection controls many observed features of stars and planets, such as magnetic fields, atmospheric jets and emitted heat flux patterns. It has long been argued that the influence of rotation on turbulent convection dynamics is governed by the ratio of the relevant global-scale forces: the Coriolis force and the buoyancy force. Here, however, we present results from laboratory and numerical experiments which exhibit transitions between rotationally dominated and non-rotating behaviour that are not determined by this global force balance. Instead, the transition is controlled by the relative thicknesses of the thermal (non-rotating) and Ekman (rotating) boundary layers. We formulate a predictive description of the transition between the two regimes on the basis of the competition between these two boundary layers. This transition scaling theory unifies the disparate results of an extensive array of previous experiments, and is broadly applicable to natural convection systems.

TURBULENT MIXING AND LAYER FORMATION IN DOUBLE-DIFFUSIVE CONVECTION: THREE-DIMENSIONAL NUMERICAL SIMULATIONS AND THEORY

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 Schwarzchild 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.

Approaching the asymptotic regime of rapidly rotating convection: Boundary layers versus interior dynamics

Rapidly rotating Rayleigh-Bénard convection is studied by combining results from direct numerical simulations (DNS), laboratory experiments, and asymptotic modeling. The asymptotic theory is shown to provide a good description of the bulk dynamics at low, but finite Rossby number. However, large deviations from the asymptotically predicted heat transfer scaling are found, with laboratory experiments and DNS consistently yielding much larger Nusselt numbers than expected. These deviations are traced down to dynamically active Ekman boundary layers, which are shown to play an integral part in controlling heat transfer even for Ekman numbers as small as 10^{-7}. By adding an analytical parametrization of the Ekman transport to simulations using stress-free boundary conditions, we demonstrate that the heat transfer jumps from values broadly compatible with the asymptotic theory to states of strongly increased heat transfer, in good quantitative agreement with no-slip DNS and compatible with the experimental data. Finally, similarly to nonrotating convection, we find no single scaling behavior, but instead that multiple well-defined dynamical regimes exist in rapidly rotating convection systems.

Numerical Models of the Geodynamo: From Fundamental Cartesian Models to 3D Simulations of Field Reversals

The series Advances in Geophysical and Environmental Mechanics and Mathematics (AGEM2) reports new developments on all topics of geophysical, environmental and astrophysical mechanics and mathematical research and teaching quickly and informally, but with high quality and the explicit aim to summarize and communicate current knowledge in an accessible way. Books published in this series are conceived as bridging material between upper level textbooks and the forefront of research to serve the following purposes: to be a compact and up-to-date source of reference on a well defined topic; to serve as an accessible introduction to a specialized field to advanced students and non-specialized researchers from related areas; to be a source of advanced teaching material for specialized seminars, courses and schools. Both monographs and multi-author volumes will be considered for publication. Edited volumes should, however, consist of a very limited number of contributions only. Proceedings will not be considered for the series. Volumes published in AGEM2 are disseminated both in print and in electronic formats, the electronic archive being available at springerlink.com. The series content is indexed, abstracted and referenced in many abstracting and information services, bibliographic networks, subscription agencies, library networks, and consortia. Proposals should be sent to the Editor or directly to the managing editor at Springer.

Chemical Transport and Spontaneous Layer Formation in Fingering Convection in Astrophysics

A region of a star that is stable to convection according to the Ledoux criterion may nevertheless undergo additional mixing if the mean molecular weight increases with radius. This process is called fingering (thermohaline) convection and may account for some of the unexplained mixing in stars such as those that have been polluted by planetary infall and those burning 3He. We propose a new model for mixing by fingering convection in the parameter regime relevant for stellar (and planetary) interiors. Our theory is based on physical principles and supported by three-dimensional direct numerical simulations. We also discuss the possibility of formation of thermocompositional staircases in fingering regions, and their role in enhancing mixing. Finally, we provide a simple algorithm to implement this theory in one-dimensional stellar codes, such as KEPLER and MESA.

Dynamics of fingering convection I: Small-scale fluxes and large-scale instabilities

(Received ?? and in revised form ??) Double-diffusive instabilities are often invoked to explain enhanced transport in stablystratified 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 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 -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.

Dynamics of Fingering Convection II: The Formation of Thermohaline Staircases

Regions of the oceans 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.

A model of the entropy flux and Reynolds stress in turbulent convection

We propose a closure model for the transport of entropy and momentum in astrophysical turbulence, intended for application to rotating stellar convective regions. Our closure model is first presented in the Boussinesq formalism, and compared with laboratory experiments and numerical simulations on the Rayleigh-Benard convection and homogeneous Rayleigh-Benard convection. The predicted angular momentum transport properties of the turbulence in the slowly rotating case recover the well-known A-effect, with an amplitude related to the convective heat flux. The model is then extended to the anelastic case. In the special case of spherical symmetry, the predicted radial heat flux is equivalent to that of mixing-length theory. For rotating stars, our model describes the coupled transport of heat and angular momentum, and provides a unified formalism in which to study both differential rotation and thermal inhomogeneities in stellar convection zones.

The effects of Ekman pumping on quasi-geostrophic Rayleigh–Bénard convection

Numerical simulations of 3D, rapidly rotating Rayleigh-Benard convection are performed using an asymptotic quasi-geostrophic model that incorporates the effects of no-slip boundaries through (i) parameterized Ekman pumping boundary conditions, and (ii) a thermal wind boundary layer that regularizes the enhanced thermal fluctuations induced by pumping. The fidelity of the model, obtained by an asymptotic reduction of the Navier-Stokes equations that implicitly enforces a pointwise geostrophic balance, is explored for the first time by comparisons of simulations against the findings of direct numerical simulations and laboratory experiments. Results from these methods have established Ekman pumping as the mechanism responsible for significantly enhancing the vertical heat transport. This asymptotic model demonstrates excellent agreement over a range of thermal forcing for Pr ~1 when compared with results from experiments and DNS at maximal values of their attainable rotation rates, as measured by the Ekman number (E ~ 10^{-7}); good qualitative agreement is achieved for Pr > 1. Similar to studies with stress-free boundaries, four spatially distinct flow morphologies exists. Despite the presence of frictional drag at the upper and/or lower boundaries, a strong non-local inverse cascade of barotropic (i.e., depth-independent) kinetic energy persists in the final regime of geostrophic turbulence and is dominant at large scales. For mixed no-slip/stress-free and no-slip/no-slip boundaries, Ekman friction is found to attenuate the efficiency of the upscale energy transport and, unlike the case of stress-free boundaries, rapidly saturates the barotropic kinetic energy. For no-slip/no-slip boundaries, Ekman friction is strong enough to prevent the development of a coherent dipole vortex condensate. Instead vortex pairs are found to be intermittent, varying in both time and strength.

A nonlinear model for rotationally constrained convection with Ekman pumping

A reduced model is developed for low-Rossby-number convection in a plane layer geometry with no-slip upper and lower boundaries held at fixed temperatures. A complete description of the dynamics requires the existence of three distinct regions within the fluid layer: a geostrophically balanced interior where fluid motions are predominantly aligned with the axis of rotation, Ekman boundary layers immediately adjacent to the bounding plates, and thermal wind layers driven by Ekman pumping in between. The reduced model uses a classical Ekman pumping parameterization to alleviate the need to resolve the Ekman boundary layers. Results are presented for both linear stability theory and a special class of nonlinear solutions described by a single horizontal spatial wavenumber. It is shown that Ekman pumping (which correlates positively with interior convection) allows for significant enhancement in the heat transport relative to that observed in simulations with stress-free boundaries. Without the intermediate thermal wind layer, the nonlinear feedback from Ekman pumping would be able to generate heat transport that diverges to infinity at finite Rayleigh number. This layer arrests this blowup, resulting in finite heat transport at a significantly enhanced value. With increasing buoyancy forcing, the heat transport transitions to a more efficient regime, a transition that is always achieved within the regime of asymptotic validity of the theory, suggesting that this behaviour may be prevalent in geophysical and astrophysical settings. As the rotation rate increases, the slope of the heat transport curve below this transition steepens, a result that is in agreement with observations from laboratory experiments and direct numerical simulations.