Turbulence-induced melting of a nonequilibrium vortex crystal in a forced thin fluid film

Abstract

To develop an understanding of recent experiments on the turbulence-induced melting of a periodic array of vortices in a thin fluid film, we perform a direct numerical simulation of the two-dimensional Navier-Stokes equations forced such that, at low Reynolds numbers, the steady state of the film is a square lattice of vortices. We find that, as we increase the Reynolds number, this lattice undergoes a series of nonequilibrium phase transitions, first to a crystal with a different reciprocal lattice and then to a sequence of crystals that oscillate in time. Initially the temporal oscillations are periodic; this periodic behaviour becomes more and more complicated, with increasing Reynolds number, until the film enters a spatially disordered nonequilibrium statistical steady that is turbulent. We study this sequence of transitions by using fluid-dynamics measures, such as the Okubo-Weiss parameter that distinguishes between vortical and extensional regions in the flow, ideas from nonlinear dynamics, e.g., \Poincare maps, and theoretical methods that have been developed to study the melting of an equilibrium crystal or the freezing of a liquid and which lead to a natural set of order parameters for the crystalline phases and spatial autocorrelation functions that characterise short- and long-range order in the turbulent and crystalline phases, respectively.