BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stephan Stellmach

Stephan Stellmach appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPphysics.flu-dynphysics.geo-ph

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2016arXiv

A nonlinear model for rotationally constrained convection with Ekman pumping

It is a well established result of linear theory that the influence of differing mechanical boundary conditions, i.e., stress-free or no-slip, on the primary instability in rotating convection becomes asymptotically small in the limit of rapid rotation. This is accounted for by the diminishing impact of the viscous stresses exerted within Ekman boundary layers and the associated vertical momentum transport by Ekman pumping. By contrast, in the nonlinear regime recent experiments and supporting simulations are now providing evidence that the efficiency of heat transport remains strongly influenced by Ekman pumping in the rapidly rotating limit. In this paper, a reduced model is developed for the case of 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 predominately 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 for spatially resolving 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 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 a heat transport that diverges to infinity. This layer arrests this blowup resulting in finite heat transport at a significantly enhanced value.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of Ekman pumping on quasi-geostrophic Rayleigh-Benard 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new model for mixing by double-diffusive convection (semi-convection). II. The transport of heat and composition through layers

Regions of stellar and planetary interiors that are unstable according to the Schwarzschild criterion, but stable according to the Ledoux criterion, are subject to a form of oscillatory double-diffusive (ODD) convection often called "semi-convection". In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and composition by ODD convection, and ultimately propose a new 1D prescription that can be used in stellar and planetary structure and evolution models. The first paper in this series demonstrated that under certain conditions ODD convection spontaneously transitions from an initially homogeneously turbulent state into a staircase of convective layers, which results in a substantial increase in the transport of heat and composition. Here, we present simulations of ODD convection in this layered regime, we describe the dynamical behavior of the layers, and we derive empirical scaling laws for the transport through layered convection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A New Model for Mixing By Double-Diffusive Convection (Semi-Convection). III. Thermal and Compositional Transport Through Non-Layered ODDC

Oscillatory double-diffusive convection (ODDC) (also known as semi- convection) refers to a type of double diffusive instability that occurs in regions of planetary and stellar interiors which have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. (2011) first showed that ODDC can either spontaneously form layers, which significantly enhance the transport of heat and chemical species compared to mi- croscopic transport, or remain in a state dominated by large scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series have focused on identifying under what condi- tions layers form (Mirouh et al. 2012), and quantifying transport through layered systems (Wood et al. 2013). Here we proceed to characterize transport through systems that are unstable to the ODDC instability, but do not undergo spon- taneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both 2D and 3D numerical simulations and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport by ODDC.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anelastic Versus Fully Compressible Turbulent Rayleigh-Bénard Convection

Numerical simulations of turbulent Rayleigh-Bénard convection in an ideal gas, using either the anelastic approximation or the fully compressible equations, are compared. Theoretically, the anelastic approximation is expected to hold in weakly superadiabatic systems with $ε= ΔT / T_r \ll 1$, where $ΔT$ denotes the superadiabatic temperature drop over the convective layer and $T_r$ the bottom temperature. Using direct numerical simulations, a systematic comparison of anelastic and fully compressible convection is carried out. With decreasing superadiabaticity $ε$, the fully compressible results are found to converge linearly to the anelastic solution with larger density contrasts generally improving the match. We conclude that in many solar and planetary applications, where the superadiabaticity is expected to be vanishingly small, results obtained with the anelastic approximation are in fact more accurate than fully compressible computations, which typically fail to reach small $ε$ for numerical reasons. On the other hand, if the astrophysical system studied contains $ε\sim O(1)$ regions, such as the solar photosphere, fully compressible simulations have the advantage of capturing the full physics. Interestingly, even in weakly superadiabatic regions, like the bulk of the solar convection zone, the errors introduced by using artificially large values for $ε$ for efficiency reasons remain moderate. If quantitative errors of the order of $10\%$ are acceptable in such low $ε$ regions, our work suggests that fully compressible simulations can indeed be computationally more efficient than their anelastic counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Double-diffusive mixing in stellar interiors in the presence of horizontal gradients

We have identified an important source of mixing in stellar radiation zones, that would arise whenever two conditions are satisfied: (1) the presence of an inverse vertical compositional gradient, and (2) the presence of density-compensating horizontal gradients of temperature (alternatively, entropy) and composition. The former can be caused naturally by any off-center burning process, by atomic diffusion, or by surface accretion. The latter could be caused by rotation, tides, meridional flows, etc. The linear instability and its nonlinear development have been well-studied in the oceanographic context. It is known to drive the formation of stacks of fingering layers separated by diffusive interfaces, called intrusions. Using 3D numerical simulations of the process in the astrophysically-relevant region of parameter space, we find similar results, and demonstrate that the material transport in the intrusive regime can be highly enhanced compared with pure diffusion, even in systems which would otherwise be stable to fingering (thermohaline) convection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The compressional beta effect: a source of zonal winds in planets?

Giant planets like Jupiter and Saturn feature strong zonal wind patterns on their surfaces. Although several different mechanisms that may drive these jets have been proposed over the last decades, the origin of the zonal winds is still unclear. Here, we explore the possibility that the interplay of planetary rotation with the compression and expansion of the convecting fluid can drive multiple deep zonal jets by a compressional Rhines-type mechanism, as originally proposed by Ingersoll and Pollard (1982). In a certain limit, this deep mechanism is shown to be mathematically analogous to the classical Rhines mechanism possibly operating at cloud level. Jets are predicted to occur on a compressional Rhines length $l_R = (2 Ω\langle H_ρ^{-1} \rangle v_{jet}^{-1} )^{-1/2}$, where $Ω$ is the angular velocity, $\langle H_ρ^{-1} \rangle$ is the mean inverse density scale height and $v_{jet}$ is the typical jet velocity. Two-dimensional numerical simulations using the anelastic approximation reveal that this mechanism robustly generates jets of the predicted width, and that it typically dominates the dynamics in systems deeper than $O(l_R)$. Potential vorticity staircases are observed to form spontaneously and are typically accompanied by unstably stratified buoyancy staircases. The mechanism only operates at large rotation rates, exceeding those typically reached in three-dimensional simulations of deep convection in spherical shells. Applied to Jupiter and Saturn, the compressional Rhines scaling reasonably fits the available observations. Interestingly, even weak vertical density variations such as those in the Earth core can give rise to a large number of jets, leading to fundamentally different flow structures than predicted by the Boussinesq models typically used in this context.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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 helium-3. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new model for mixing by double-diffusive convection (semi-convection): I. The conditions for layer formation

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 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. (2011) 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 numerically accessible "planetary" parameter regime, and can easily 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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 and numerical experimental results on 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 Lambda-effect, with an amplitude uniquely related to the convective heat flux. The model is then extended to the anelastic case as well as the fully compressible 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerically determined transport laws for fingering ("thermohaline") convection in astrophysics

We present the first three-dimensional simulations of fingering convection performed in a parameter regime close to the one relevant for astrophysics, and reveal the existence of simple asymptotic scaling laws for turbulent heat and compositional transport. These laws can straightforwardly be extrapolated to the true astrophysical regime. Our investigation also indicates that thermocompositional "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 RGB stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Turbulent mixing and layer formation in double-diffusive convection: 3D 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.

0 citations0 reviewsNo code linkedOpen access