We investigate the vorticity-preserving properties of the compressible, second-order residual-based scheme, "RBV2". The scheme has been extensively tested on hydrodynamical problems, and has been shown to exhibit remarkably accurate results on the propagation of inviscid compressible vortices, airfoil-vortex interactions on a curvilinear mesh, vortex mergers in an astrophysical accretion disk, and the establishment of a two-dimensional inverse cascade in high-resolution turbulent simulations. Here, we demonstrate that RBV2 sustains the analytic solution for a one-dimensional shear flow. We assess the fidelity by which the algorithm maintains a skewed shear flow, and present convergence tests to quantify the magnitude of the expected numerical dispersion. We propose an adjustment to the dissipation in the algorithm that retains the vorticity-preserving qualities, and accurately incorporates external body forces, and demonstrate that it indefinitely maintains a steady-state hydrostatic equilibrium between a generic acceleration and a density gradient. We present a novel numerical assessment of vorticity preservation for discrete-wavenumber vortical modes up to the Nyquist wavenumber. We apply this assessment to RBV2 in order to quantify the extent to which the scheme preserves vorticity for the full Euler equations. We find that RBV2 perfectly preserves vorticity for modes with symmetric wavenumbers, i.e., $k_x=k_y$, and that the error increases with asymmetry. We simulate the dynamical interaction of vortices in a protoplanetary disk to demonstrate the utility of the updated scheme for rendering astrophysical flows replete with vortices and turbulence. We conclude that RBV2 is a competitive treatment for evolving vorticity-dominated astrophysical flows with minimal dissipation.

A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to as Z-modes. The associated critical layer is the location where the Doppler shifted frequency of a distant Rossby wave equals the local Brunt-Vaisala frequency. The minimum required Rossby number for instability, Ro= 0.2, is confirmed for parameter values reported in the literature. It is also found that the shear layer supports the instability in the limit where stratification vanishes. The zombie vortex instability is examined in a jet model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity forcing due to unstable Z-modes is shown to result in the creation of a jet flow at the critical layer emerging as the result of the competition between the vertical lifting of perturbation radial vorticity and the radial transport of perturbation vertical vorticity. We find that the picture of this instability leading to a form of nonlinearly driven self-replicating pattern of creation and destruction is warranted: a parent jet spawns a growing child jet at associated critical layers. A mature child jet creates a next generation of child jets at associated critical layers while simultaneously destroying its parent jet via the Rossby wave instability.

The Goldreich and Ward (1973) (axisymmetric) gravitational instability of a razor thin particle layer occurs when the Toomre parameter $Q_T \equiv c_p \Omega_0 / \pi G \Sigma_p < 1$ ($c_p$ being the particle dispersion velocity). Ward(1976,2000) extended this analysis by adding the effect of gas drag upon particles and found that even when $Q_T > 1$, sufficiently long waves were always unstable. Youdin (2005a,b) carried out a detailed analysis and showed that the instability allows chondrule-sized ($\sim 1 $ mm) particles to undergo radial clumping with reasonable growth times even in the presence of a moderate amount of turbulent stirring. The analysis of Youdin includes the role of turbulence in setting the thickness of the dust layer and in creating a turbulent particle pressure in the momentum equation. However, he ignores the effect of turbulent mass diffusivity on the disturbance wave. Here we show that including this effect reduces the growth-rate significantly, by an amount that depends on the level of turbulence, and reduces the maximum intensity of turbulence the instability can withstand by 1 to 3 orders of magnitude. The instability is viable only when turbulence is extremely weak and the solid to gas surface density of the particle layer is considerably enhanced over minimum-mass-nebula values. A simple mechanistic explanation of the instability shows how the azimuthal component of drag promotes instability while the radial component hinders it. A gravito-diffusive overstability is also possible but never realized in the nebula models.

Several mechanisms have been identified that create dense particle clumps in the solar nebula. The present work is concerned with the gravitational collapse of such clumps, idealized as being spherically symmetric. Calculations using the two-fluid model are performed (almost) up to the time when a central density singularity forms. The end result of the study is a parametrization for this time, in order that it may be compared with timescales for various disruptive effects to which clumps may be subject. An important effect is that as the clump compresses, it also compresses the gas due to drag. This increases gas pressure which retards particle collapse and leads to oscillation in the size and density of the clump. The ratio of gravitational force to gas pressure gives a two-phase Jeans parameter, $J_t$, which is the classical Jeans parameter with the sound speed replaced by an the wave speed in a coupled two-fluid medium. Its use makes the results insensitive to the initial density ratio of particles to gas as a separate parameter. An ordinary differential equation model is developed which takes the form of two coupled non-linear oscillators and reproduces key features of the simulations. Finally, a parametric study of the time to collapse is performed and a formula (fit to the simulations) is developed. In the incompressible limit $J_t \to 0$, collapse time equals sedimentation time. As $J_t$ increases, the collapse time decreases roughly linearly with $J_t$ until $J_t \gtrsim 0.4$ when it becomes approximately equal to the dynamical time.

We outline a scenario which traces a direct path from freely-floating nebula particles to the first 10-100km-sized bodies in the terrestrial planet region, producing planetesimals which have properties matching those of primitive meteorite parent bodies. We call this "primary accretion". The scenario draws on elements of previous work, and introduces a new critical threshold for planetesimal formation. We presume the nebula to be weakly turbulent, which leads to dense concentrations of aerodynamically size-sorted particles having properties like those observed in chondrites. The fractional volume of the nebula occupied by these dense zones or clumps obeys a probability distribution as a function of their density, and the densest concentrations have particle mass density 100 times that of the gas. However, even these densest clumps are prevented by gas pressure from undergoing gravitational instability in the traditional sense (on a dynamical timescale). While in this state of arrested development, they are susceptible to disruption by the ram pressure of the differentially orbiting nebula gas. However, self-gravity can preserve sufficiently large and dense clumps from ram pressure disruption, allowing their entrained particles to sediment gently but inexorably towards their centers, producing 10-100 km "sandpile" planetesimals. Localized radial pressure fluctuations in the nebula, and interactions between differentially moving dense clumps, will also play a role that must be allowed for in future studies. The scenario is readily extended from meteorite parent bodies to primary accretion throughout the solar system.