Evolution of a vortex in a strongly stratified shear flow. Part 1. Asymptotic analysis

International audienceIn this paper, we investigate the dynamics of an initially vertical vortex embedded in a shear flow in a stratified fluid by means of a long-wavelength analysis. The main goal is to determine, whether or not, the Kelvin-Helmholtz instability can develop unconditionally as speculated by Lilly (1983). The analysis is performed in the case of the Lamb-Oseen vortex profile and a shear flow uniform in the horizontal and varying sinusoidally along the vertical using the assumptionk z a 0 F h ⌧ 1, wherek z is the vertical wavenumber, a 0 the vortex radius and F h the horizontal Froude number based on the circulation of the vortex. The results show that the vortex axis is not only advected in the direction of the shear flow but also in the perpendicular direction owing to the self-induced motion of the vortex. In addition, internal waves are transiently excited at the beginning, generating an initial non-hydrostatic regime. Their relative e↵ects on the displacements of the vortex axis are weak except initially. The angular velocity of the vortex decays owing to a dynamic e↵ect and viscous e↵ects related to the vertical shear. The former e↵ect is due to the squeezing of isopycnals in the vortex core which implies a decrease of the vertical vorticity to satisfy potential vorticity conservation. In addition, a horizontal velocity field with an azimuthal wavenumber m = 2 is generated meaning that the shape of the vortex becomes slightly elliptical. We further show that the minimum asymptotic Richardson number is bounded min(Ri) > 3.43 whenk z a 0 F h ⌧ 1 and therefore can not go below the critical value 1/4. This is because the growth of the vertical shear of the horizontal velocity of the vortex saturates owing to the decay of its angular velocity and because the squeezing of isopycnals increases the stratification strength. This suggests that the shear instability can not always develop in strongly stratified flows, contrary to the conjecture of Lilly (1983). These predictions will be tested against Direct Numerical Simulations in part 2