There is considerable interest in hydrodynamic instabilities in dead zones of protoplanetary disks as a mechanism for driving angular momentum transport and as a source of particle-trapping vortices to mix chondrules and incubate planetesimal formation. We present simulations with a pseudo-spectral anelastic code and with the compressible code Athena, showing that stably stratified flows in a shearing, rotating box are violently unstable and produce space-filling, sustained turbulence dominated by large vortices with Rossby numbers of order ∼0.2 − 0.3. This Zombie Vortex Instability (ZVI) is observed in both codes and is triggered by Kolmogorov turbulence with Mach numbers less than ∼0.01. It is a common view that if a given constant density flow is stable, then stable vertical stratification should make the flow even more stable.

A previously unknown instability creates space-filling lattices of 3D vortices in linearly stable, rotating, stratified shear flows. The instability starts from an easily excited critical layer. The layer intensifies by drawing energy from the background shear and rolls up into vortices that excite new critical layers and vortices. The vortices self-similarly replicate to create lattices of turbulent vortices. The vortices persist for all time. This self-replication occurs in stratified Couette flows and in the dead zones of protoplanetary disks where it can destabilize Keplerian flows.

In the traditional hybrid RANS/LES approaches for the simulation of wall-bounded fluid turbulence, such as detached-eddy simulation (DES), the whole flow domain is divided into an inner layer and an outer layer. Typically the Reynolds-averaged Navier–Stokes (RANS) equations are used for the inner layer, while large-eddy simulation (LES) is used for the outer layer. The transition from the inner-layer solution to the outer-layer solution is often problematic due to the lack of small-scale dynamics in the RANS region. In this paper, we propose to simulate the whole flow domain by large-eddy simulation while enforcing a Reynolds-stress constraint on the subgrid-scale (SGS) stress model in the inner layer. Both the algebraic eddy-viscosity model and the one-equation Spalart–Allmaras (SA) model have been used to constrain the Reynolds stress in the inner layer. In this way, we improve the LES methodology by allowing the mean flow of the inner layer to satisfy the RANS solution while small-scale dynamics is included. We validate the Reynolds-stress-constrained large-eddy simulation (RSC-LES) model by simulating three-dimensional turbulent channel flow and flow past a circular cylinder. Our model is able to predict mean velocity, turbulent stress and skin-friction coefficients more accurately in turbulent channel flow and to estimate the pressure coefficient after separation more precisely in flow past a circular cylinder compared with the pure dynamic Smagorinsky model (DSM) and DES using the same grid resolution. Furthermore, the computational cost of the RSC-LES is almost the same as that of DES.

In Zombie Vortex Instability (ZVI), perturbations excite critical layers in stratified, rotating shear flow (as in protoplanetary disks), causing them to generate vortex layers, which roll-up into anticyclonic zombie vortices and cyclonic vortex sheets. The process is self-sustaining as zombie vortices perturb new critical layers, spawning a next generation of zombie vortices. Here, we focus on two issues: the minimum threshold of perturbations that trigger self-sustaining vortex generation, and the properties of the late-time zombie turbulence on large and small scales. The critical parameter that determines whether ZVI is triggered is the magnitude of the vorticity on the small scales (and not velocity); the minimum Rossby number needed for instability is Ro crit ∼ 0.2 for β ≡ N/Ω = 2, where N is the Brunt-Väisälä frequency. While the threshold is set by vorticity, it is useful to infer a criterion on the Mach number; for Kolmogorov noise, the critical Mach number scales with Reynolds number: M a crit ∼ Ro crit Re −1/2 . In protoplanetary disks, this is M a crit ∼ 10 −6 . On large scales, zombie turbulence is characterized by anticyclones and cyclonic sheets with typical Rossby number ∼0.3. The spacing of the cyclonic sheets and anticyclones appears to have a "memory" of the spacing of the critical layers. On the small scales, zombie turbulence has no memory of the initial conditions and has a Kolmogorov-like energy spectrum. While our earlier work was in the limit of uniform stratification, we have demonstrated that ZVI works for non-uniform Brunt-Väisälä frequency profiles that may be found in protoplanetary disks.

The Zombie Vortex Instability (ZVI) occurs in the dead zones of protoplanetary disks (PPDs) where perturbations excite baroclinic critical layers, generating "zombie" vortices and turbulence. In this work, we investigate ZVI with nonuniform vertical stratification; while ZVI is triggered in the stratified regions away from the midplane, the subsequent turbulence propagates into and fills the midplane. ZVI turbulence alters the background Keplerian shear flow, creating a steady-state zonal flow. Intermittency is observed, where the flow cycles through near-laminar phases of zonal flow punctuated by chaotic bursts of new vortices. ZVI persists in the presence of radiative damping, as long as the thermal relaxation timescale is more than a few orbital periods. We refute the premature claim by Lesur & Latter (2016) that radiative damping inhibits ZVI for disk radii r 0.3 au. Their conclusions were based on unrealistically short cooling times using opacities with virtually no grain growth. We explore different grain growth and vertical settling scenarios, and find that the gas and dust in off-midplane regions are not necessarily in local thermodynamic equilibrium (LTE) with each other. In such cases, thermal relaxation timescales can be orders of magnitude longer than the optically thin cooling times assuming LTE because of the finite time for energy to be exchanged between gas and dust grains via collisions. With minimal amounts of grain growth and dust settling, the off-midplane regions of disks are susceptible to ZVI and much of the planet-forming regions can be filled with zombie vortices and turbulence.1 The Reynolds number is the ratio of the rate of change of momentum via advection to the magnitude of the viscous force: Re ≡ |(v · ∇)v|/|ν∇ 2 v| ≈ U L/ν ∼ H/ mf p , where L is a characteristic length which we take to be equal to the gas pressure scale height H, U is a characteristic velocity which we take to be equal to H times the orbital frequency Ω, and where ν is the molecular viscosity and ν ∼ cs mf p , where cs is the sound speed and mf p is the mean free path between collisions of gas molecules.