Keith Julien

Rapidly rotating turbulent Rayleigh-Bénard convection

Turbulent Boussinesq convection under the influence of rapid rotation (i.e. with comparable characteristic rotation and convection timescales) is studied. The transition to turbulence proceeds through a relatively simple bifurcation sequence, starting with unstable convection rolls at moderate Rayleigh ( Ra ) and Taylor numbers ( Ta ) and culminating in a state dominated by coherent plume structures at high Ra and Ta . Like non-rotating turbulent convection, the rapidly rotating state exhibits a simple power-law dependence on Ra for all statistical properties of the flow. When the fluid layer is bounded by no-slip surfaces, the convective heat transport ( Nu − 1, where Nu is the Nusselt number) exhibits scaling with Ra 2/7 similar to non-rotating laboratory experiments. When the boundaries are stress free, the heat transport obeys ‘classical’ scaling ( Ra 1/3 ) for a limited range in Ra , then appears to undergo a transition to a different law at Ra ≈ 4 × 10 7 . Important dynamical differences between rotating and non-rotating convection are observed: aside from the (expected) differences in the boundary layers due to Ekman pumping effects, angular momentum conservation forces all plume structures created at flow-convergent sites of the heated and cooled boundaries to spin-up cyclonically; the resulting plume/cyclones undergo strong vortex-vortex interactions which dramatically alter the mean state of the flow and result in a finite background temperature gradient as Ra → ∞, holding Ra / Ta fixed.

BAROCLINIC VORTICITY PRODUCTION IN PROTOPLANETARY DISKS. I. VORTEX FORMATION

The formation of vortices in protoplanetary disks is explored via pseudospectral numerical simulations of an anelastic-gas model. This model is a coupled set of equations for vorticity and temperature in two dimensions that includes baroclinic vorticity production and radiative cooling. Vortex formation is unambiguously shown to be caused by baroclinicity, because (1) these simulations have zero initial perturbation vorticity and a nonzero initial temperature distribution, and (2) turning off the baroclinic term halts vortex formation, as shown by an immediate drop in kinetic energy and vorticity. Vortex strength increases with larger background temperature gradients, warmer background temperatures, larger initial temperature perturbations, higher Reynolds number, and higher resolution. In the simulations presented here, vortices form when the background temperatures are ~200 K and vary radially as r-0.25, the initial vorticity perturbations are zero, the initial temperature perturbations are 5% of the background, and the Reynolds number is 109. A sensitivity study consisting of 74 simulations showed that as resolution and Reynolds number increase, vortices can form with smaller initial temperature perturbations, lower background temperatures, and smaller background temperature gradients. For the parameter ranges of these simulations, the disk is shown to be convectively stable by the Solberg-Hoiland criteria.

BAROCLINIC VORTICITY PRODUCTION IN PROTOPLANETARY DISKS. II. VORTEX GROWTH AND LONGEVITY

The factors affecting vortex growth in convectively stable protoplanetary disks are explored using numerical simulations of a two-dimensional anelastic-gas model that includes baroclinic vorticity production and radiative cooling. The baroclinic feedback, in which anomalous temperature gradients produce vorticity through the baroclinic term and vortices then reinforce these temperature gradients, is found to be an important process in the rate of growth of vortices in the disk. Factors that strengthen the baroclinic feedback include fast radiative cooling, high thermal diffusion, and large radial temperature gradients in the background temperature. When the baroclinic feedback is sufficiently strong, anticyclonic vortices form from initial random perturbations and maintain their strength for the duration of the simulation, for over 600 orbital periods. Based on both simulations and a simple vortex model, we find that the local angular momentum transport due to a single vortex may be inward or outward, depending on its orientation. The global angular momentum transport is highly variable in time and is sometimes negative and sometimes positive. This result is for an anelastic-gas model and does not include shocks that could affect angular momentum transport in a compressible-gas disk.

Statistical and physical balances in low Rossby number Rayleigh–Bénard convection

Rapidly rotating Rayleigh–Bénard convection is studied using an asymptotically reduced equation set valid in the limit of low Rossby numbers. Four distinct dynamical regimes are identified: a disordered cellular regime near threshold, a regime of weakly interacting convective Taylor columns at larger Rayleigh numbers, followed for yet larger Rayleigh numbers by a breakdown of the convective Taylor columns into a disordered plume regime characterized by reduced efficiency and finally by geostrophic turbulence. The transitions are quantified by examining the properties of the horizontally and temporally averaged temperature and thermal dissipation rate. The maximum of the thermal dissipation rate is used to define the width of the thermal boundary layer. In contrast to the non-rotating Rayleigh–Bénard convection, the temperature drop across this layer decreases monotonically with increasing Rayleigh number and does not saturate. The breakdown of the convective Taylor column regime is attributed to the onset of convective instability of the thermal boundary layer and confirmed using the explicit linear stability analysis. Horizontal spectra of the vorticity, vertical velocity and temperature fluctuations are computed and their evolution with time is elucidated. A large-scale barotropic mode evolves from random initial conditions on an extremely long time scale and leads to continued evolution of the nominally saturated Nusselt number and its variance over very long times. The results are used to provide insights into the dynamics of rapidly rotating convection outside the asymptotic regime described by the reduced equations.

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.

Langmuir–Submesoscale Interactions: Descriptive Analysis of Multiscale Frontal Spindown Simulations

AbstractThe interactions between boundary layer turbulence, including Langmuir turbulence, and submesoscale processes in the oceanic mixed layer are described using large-eddy simulations of the spindown of a temperature front in the presence of submesoscale eddies, winds, and waves. The simulations solve the surface-wave-averaged Boussinesq equations with Stokes drift wave forcing at a resolution that is sufficiently fine to capture small-scale Langmuir turbulence. A simulation without Stokes drift forcing is also performed for comparison. Spatial and spectral properties of temperature, velocity, and vorticity fields are described, and these fields are scale decomposed in order to examine multiscale fluxes of momentum and buoyancy. Buoyancy flux results indicate that Langmuir turbulence counters the restratifying effects of submesoscale eddies, leading to small-scale vertical transport and mixing that is 4 times greater than in the simulations without Stokes drift forcing. The observed fluxes are also sh...

Reduced models for fluid flows with strong constraints

The presence of a dominant balance in the equations for fluid flow can be exploited to derive an asymptotically exact but simpler set of governing equations. These permit semianalytical and/or numerical explorations of parameter regimes that would otherwise be inaccessible to direct numerical simulation. The derivation of the resulting reduced models is illustrated here for (i) rapidly rotating convection in a plane layer, (ii) convection in a strong magnetic field, and (iii) the magnetorotational instability in accretion disks and the results used to extend our understanding of these systems in the strongly nonlinear regime.

Plumes in rotating convection. Part 1. Ensemble statistics and dynamical balances

Atmospheric and oceanic convection often occurs over areas occupied by many localized circulation elements known as plumes . The convective transports therefore may depend not only on the individual elements, but also on the interactions between plumes and the turbulent environment created by other plumes. However, many attempts to understand these plumes focus on individual isolated elements, and the behaviour of an ensemble is not understood. Geophysical convection may be influenced by rotation when the transit time of a convecting element is long compared to an inertial period (for example in deep oceanic convection). Much recent attention has been given to the effect of rotation on individual plumes, but the role of rotation in modifying the behaviour of an ensemble is not fully understood. Here we examine the behaviour of plumes within an ensemble, both with and without rotation, to identify the influence of rotation on ensemble plume dynamics. We identify the coherent structures (plumes) present in numerical solutions of turbulent Rayleigh–Benard convection, a canonical example of a turbulent plume ensemble. We use a conditional sampling compositing technique to extract the typical structure in both non-rotating and rotating solutions. The dynamical balances of these composite plumes are evaluated and compared with entraining plume models. We find many differences between non-rotating and rotating plumes in their transports of mass, buoyancy and momentum. As shown in previous studies, the expansion of the turbulent plume by entrainment of exterior fluid is suppressed by strong rotation. Our most significant new result is quantification of the continuous mixing between the plume and ambient fluid which occurs at high rotation without any net changes in plume volume. This mixing is generated by the plume–plume interactions and acts to reduce the buoyancy anomaly of the plume. By contrast, in the non-rotating case, no such loss of buoyancy by mixing occurs. As a result, the total buoyancy transport by upwardly moving plumes diminishes across the layer in the rotating case, while remaining approximately constant in the non-rotating case. At high values of rotation, the net vertical acceleration is considerably reduced compared to the non-rotating case due to loss of momentum through entrainment and mixing and a decelerating pressure gradient which partially balances the buoyancy-driven acceleration of plumes. As a result of the dilution of buoyancy, the pressure-gradient deceleration and the loss of momentum due to mixing with the environment in the rotating solutions, the conversion of potential energy to kinetic energy is significantly less than that of non-rotating plumes. The combination of efficient lateral mixing and slow vertical movement by the plumes accounts for the unstable mean temperature gradient that occurs in rotating Rayleigh–Benard convection, while the less penetrative convection found at low Rossby number is a consequence of the reduced kinetic energy transport. Within the ensemble of plumes identified by the conditional sampling algorithm, distributions of vertical velocity, buoyancy and vorticity mimic those of the volume as a whole. Plumes cover a small fraction of the total area, yet account for most of the vertical heat flux.

Generalized quasi-geostrophy for spatially anisotropic rotationally constrained flows

Closed reduced equations analogous to the quasi-geostrophic equations are derived in the extratropics for small Rossby numbers and vertical scales that are comparable to or much larger than horizontal scales. On these scales, significant vertical motions are permitted and found to couple to balanced geostrophic dynamics. In the equatorial regions, similar reduced equations are derived for meridional scales much larger than the vertical and zonal scales. These equations are derived by a systematic exploration of different aspect ratios, and Froude and buoyancy numbers, and offer advantages similar to the standard quasi-geostrophic equations for studies of smaller-scale processes and/or of the equatorial regions.

Strongly nonlinear convection cells in a rapidly rotating fluid layer: the tilted f-plane

Investigation of the linear stability problem for rapidly rotating convection on an f-plane has revealed the existence of two distinct scales in the vertical structure of the critical eigenfunctions: a small length scale whose vertical wavenumber kz is comparable with the large horizontal wavenumber k⊥ selected at onset, and a large-scale modulation which forms an envelope on the order of the layer depth d. The small-scale structure in the vertical results from a geostrophic balance imposed by the Taylor–Proudman constraint. This primary balance forces rotational alignment and confines fluid motions to planes perpendicular to the rotation axis. For convective transport in the vertical this constraint must be relaxed. This is achieved by molecular dissipation which allows weak upward (downward) spiralling of hot (cold) fluid elements across the Taylor–Proudman planes and results in a large-scale vertical modulation of the Taylor columns.In the limit of fast rotation (i.e. large Taylor number) a multiple-scales analysis leads to the determination of a critical Rayleigh number as a function of wavenumber, roll orientation and the tilt angle of the f-plane. The corresponding critical eigenfunction represents the core solution; matching to passive Ekman boundary layers is required for a complete solution satisfying boundary conditions.An extension of this analysis, introduced by Bassom & Zhang (1994), is used to describe strongly nonlinear two-dimensional convection, characterized by significant departures of the mean thermal field from its conduction profile. The analysis requires the solution of a nonlinear eigenvalue problem for the Nusselt number (for steady convection) and the Nusselt number and oscillation frequency (for the overstable problem). The solutions of this problem are used to calculate horizontal and vertical heat fluxes, as well as Reynolds stresses, as functions of both the latitude and roll orientation in the horizontal, and these are used to calculate self-consistently north–south and east–west mean flows. These analytical predictions are in good agreement with the results of three-dimensional simulations reported by Hathaway & Somerville (1983).