Diffusive effects in local instabilities of a baroclinic axisymmetric vortex

We present a local stability analysis of an idealized model of the stratified vortices that appear in geophysical settings. The base flow comprises an axisymmetric vortex with background rotation and an out-of-plane stable stratification, and a radial stratification in the thermal wind balance with the out-of-plane momentum gradient. Solving the local stability equations along fluid particle trajectories in the base flow, the dependence of short-wavelength instabilities on the Schmidt number $Sc$ (ratio between momentum and mass diffusivities) is studied, in the presence of curvature effects. In the diffusion-free limit, the well-known symmetric instability is recovered. In the viscous, double-diffusive regime, instability characteristics are shown to depend on three non-dimensional parameters (including $Sc$ ), and two different instabilities are identified: (i) a monotonic instability (same as symmetric instability at $Sc = 1$ ), and (ii) an oscillatory instability (absent at $Sc = 1$ ). Separating the base flow and perturbation characteristics, two each of base flow and perturbation parameters (apart from $Sc$ ) are identified, and the entire parameter space is explored for the aforementioned instabilities. In comparison with $Sc = 1$ , monotonic and oscillatory instabilities are shown to significantly expand the instability region in the space of base flow parameters as $Sc$ moves away from unity. Neutral stability boundaries on the plane of $Sc$ and a modified gradient Richardson number are then identified for both these instabilities. In the absence of curvature effects, our results are shown to be consistent with previous studies based on normal mode analysis, thus establishing that the local stability approach is well suited to capturing symmetric and double-diffusive instabilities. The paper concludes with a discussion of curvature effects, and the likelihood of monotonic and oscillatory instabilities in typical oceanic settings.

Effects of Schmidt number on the short-wavelength instabilities in stratified vortices

We present a local stability analysis to investigate the effects of differential diffusion between momentum and density (quantified by the Schmidt number $Sc$) on the three-dimensional, short-wavelength instabilities in planar vortices with a uniform stable stratification along the vorticity axis. Assuming small diffusion in both momentum and density, but arbitrary values for $Sc$, we present a detailed analytical/numerical analysis for three different classes of base flows: (i) an axisymmetric vortex, (ii) an elliptical vortex and (iii) the flow in the neighbourhood of a hyperbolic stagnation point. While a centrifugally stable axisymmetric vortex remains stable for any $Sc$, it is shown that $Sc$ can have significant effects in a centrifugally unstable axisymmetric vortex: the range of unstable perturbations increases when $Sc$ is taken away from unity, with the extent of increase being larger for $Sc\ll 1$ than for $Sc\gg 1$. Additionally, for $Sc>1$, we report the possibility of oscillatory instability. In an elliptical vortex with a stable stratification, $Sc\neq 1$ is shown to non-trivially influence the three different inviscid instabilities (subharmonic, fundamental and superharmonic) that have been previously reported, and also introduce a new branch of oscillatory instability that is not present at $Sc=1$. The unstable parameter space for the subharmonic (instability IA) and fundamental (instability IB) inviscid instabilities are shown to be significantly increased for $Sc<1$ and $Sc>1$, respectively. Importantly, for sufficiently small and large $Sc$, respectively, the maximum growth rate for instabilities IA and IB occurs away from the inviscid limit. The new oscillatory instability (instability III) is shown to occur only for sufficiently small $Sc<1$, the signature of which is nevertheless present with zero growth rate in the inviscid limit. The Schmidt number is then shown to play no role in the evolution of transverse perturbations on the flow around a hyperbolic stagnation point with a stable stratification. We conclude by discussing the physical length scales associated with the $Sc\neq 1$ instabilities, and their potential relevance in various realistic settings.

Early azimuthal instability during drop impact

When a drop impacts on a liquid surface its bottom is deformed by lubrication pressure and it entraps a thin disc of air, thereby making contact along a ring at a finite distance from the centreline. The outer edge of this contact moves radially at high speed, governed by the impact velocity and bottom radius of the drop. Then at a certain radial location an ejecta sheet emerges from the neck connecting the two liquid masses. Herein, we show the formation of an azimuthal instability at the base of this ejecta, in the sharp corners at the two sides of the ejecta. They promote regular radial vorticity, thereby breaking the axisymmetry of the motions on the finest scales. The azimuthal wavenumber grows with the impact Weber number, based on the bottom curvature of the drop, reaching over 400 streamwise streaks around the periphery. This instability occurs first at Reynolds numbers ($Re$) of ${\sim}7000$, but for larger $Re$ is overtaken by the subsequent axisymmetric vortex shedding and their interactions can form intricate tangles, loops or chains.

Advective balance in pipe-formed vortex rings

Vorticity distributions in axisymmetric vortex rings produced by a piston–pipe apparatus are numerically studied over a range of Reynolds numbers, $Re$, and stroke-to-diameter ratios, $L/D$. It is found that a state of advective balance, such that $\unicode[STIX]{x1D701}\equiv \unicode[STIX]{x1D714}_{\unicode[STIX]{x1D719}}/r\approx F(\unicode[STIX]{x1D713},t)$, is achieved within the region (called the vortex ring bubble) enclosed by the dividing streamline. Here $\unicode[STIX]{x1D701}\equiv \unicode[STIX]{x1D714}_{\unicode[STIX]{x1D719}}/r$ is the ratio of azimuthal vorticity to cylindrical radius, and $\unicode[STIX]{x1D713}$ is the Stokes streamfunction in the frame of the ring. Some, but not all, of the $Re$ dependence in the time evolution of $F(\unicode[STIX]{x1D713},t)$ can be captured by introducing a scaled time $\unicode[STIX]{x1D70F}=\unicode[STIX]{x1D708}t$, where $\unicode[STIX]{x1D708}$ is the kinematic viscosity. When $\unicode[STIX]{x1D708}t/D^{2}\gtrsim 0.02$, the shape of $F(\unicode[STIX]{x1D713})$ is dominated by the linear-in-$\unicode[STIX]{x1D713}$ component, the coefficient of the quadratic term being an order of magnitude smaller. An important feature is that, as the dividing streamline ($\unicode[STIX]{x1D713}=0$) is approached, $F(\unicode[STIX]{x1D713})$ tends to a non-zero intercept which exhibits an extra $Re$ dependence. This and other features are explained by a simple toy model consisting of the one-dimensional cylindrical diffusion equation. The key ingredient in the model responsible for the extra $Re$ dependence is a Robin-type boundary condition, similar to Newton’s law of cooling, that accounts for the edge layer at the dividing streamline.