Phase transition in time-reversible Navier-Stokes equations

We present a comprehensive study of the statistical features of a three-dimensional (3D) time-reversible truncated Navier-Stokes (RNS) system, wherein the standard viscosity ν is replaced by a fluctuating thermostat that dynamically compensates for fluctuations in the total energy. We analyze the statistical features of the RNS steady states in terms of a non-negative dimensionless control parameter R_{r}, which quantifies the balance between the fluctuations of kinetic energy at the forcing length scale ℓ_{f} and the total energy E_{0}. For small R_{r}, the RNS equations are found to produce "warm" stationary statistics, e.g., characterized by the partial thermalization of the small scales. For large R_{r}, the stationary solutions have features akin to standard hydrodynamic ones: they have compact energy support in k space and are essentially insensitive to the truncation scale k_{max}. The transition between the two statistical regimes is observed to be smooth but rather sharp. Using insights from a diffusion model of turbulence (Leith model), we argue that the transition is in fact akin to a continuous second-order phase transition, where R_{r} indeed behaves as a thermodynamic control parameter, e.g., a temperature. A relevant order parameter can be suitably defined in terms of a (normalized) enstrophy, while the symmetry-breaking parameter h is identified as (one over) the truncation scale k_{max}. We find that the signatures of the phase transition close to the critical point R_{r}^{★} can essentially be deduced from a heuristic mean-field Landau free energy. This point of view allows us to reinterpret the relevant asymptotics in which the dynamical ensemble equivalence conjectured by Gallavotti [Phys. Lett. A 223, 91 (1996)PYLAAG0375-960110.1016/S0375-9601(96)00729-3] could hold true. We argue that Gallavotti's limit is precisely the joint limit R_{r}→[over >]R_{r}^{★} and h→[over >]0, with the overset symbol ">" indicating that those limits are approached from above. The limit therefore relates to the statistical features at the critical point. In this regime, our numerics indicate that the low-order statistics of the 3D RNS are indeed qualitatively similar to those observed in direct numerical simulations of the standard Navier-Stokes equations with viscosity chosen so as to match the average value of the reversible thermostat. This result suggests that Gallavotti's equivalence conjecture could indeed be of relevance to model 3D turbulent statistics, and provides a clear guideline for further numerical investigations involving higher resolutions.