Abstract
The butterfly effect is today commonly identified with the sensitive dependence of deterministic chaotic systems upon initial conditions. However, this is only one facet of the notion of unpredictability pioneered by Lorenz, who actually predicted that multiscale fluid flows could spontaneously lose their deterministic nature and become intrinsically random. This effect, which is radically different from chaos, have remained out of reach for detailed physical observations. Here we show that this scenario is inherent to the elementary Kelvin–Helmholtz hydrodynamical instability of an initially singular shear layer. We moreover provide evidence that the resulting macroscopic flow displays universal statistical properties that are triggered by, but independent of specific physical properties at microscales. This spontaneous stochasticity is interpreted as an Eulerian counterpart to Richardson’s relative dispersion of Lagrangian particles, giving substance to the intrinsic nature of randomness in turbulence.
Introduction
In his seminal works on unpredictability of atmospheric motion, Edward N. Lorenz^{1,2} formulated a visionary conjecture on the sensitive dependence of deterministic systems upon initial perturbations or errors. He distinguished between two kinds of unstable behaviors: Either the mismatch between two replica of the same system can be made arbitrarily small by sufficiently reducing their initial discrepancy or, alternatively, the two systems reach diverging states, no matter how small they initially differ from each other. The first type of behavior arises in chaotic systems, wherein small perturbations are amplified exponentially in time. This phenomenon is now celebrated as the “butterfly effect”^{3} and has been widely used to link unpredictability and complexity^{4}.
The second behavior corresponds to a much more intrinsic materialization of unpredictability. In Lorenz’s own words: “certain formally deterministic fluid systems which possess many scales of motion are observationally indistinguishable from indeterministic systems; specifically, that two states of the system differing initially by a small observational error will evolve into two states differing as greatly as randomly chosen states of the system within a finitetime interval, which cannot be lengthened by reducing the amplitude of the initial error”^{2}. According to Lorenz, this second type of unpredictability pertains to fluid flows with a sufficient amount of kinetic energy at small scales, meaning that the velocity field should have a singular spatial structure. As pointed out by Palmer et al.^{5}, the dynamical formulation of such fluid systems must be illposed: solutions do not depend continuously on initial conditions, therefore, allowing a finitetime separation between initially undistinguishable systems. In spite of being supported by phenomenological arguments^{6,7,8,9,10,11,12,13}, a clear physical evidence is still lacking whether such a behavior indeed arises in genuine fluid flows, hereby making them infinitely less predictable than chaos.
Physically, this scenario could be the relevant one to address fundamental instabilities in classical fluid mechanics. An example is the Kelvin–Helmholtz (KH) instability, which describes the growth of a shear layer from an initially discontinuous velocity profile^{14}. This instability is an important constituent in dynamical phenomena^{15} ranging from the microworld of quantum fluids^{16,17} to macroscopic motions in Earth’s atmosphere and oceans^{18}, extending further to astrophysical problems of supernovae^{19} and interstellar clouds^{20}. While formally deterministic, the singular shear layer problem is however illposed in ideal fluid mechanics^{21}. The dominant viewpoint generally attributed to the French mathematician J. Hadamard is that such problems are physically meaningless, unless suitably regularized, so that solutions become uniquely determined and exhibit continuous dependence upon a class of initial data^{22,23}. This viewpoint has motivated numerous studies employing analytic perturbation of the vortex sheet initial datum, that carefully analyze finitetime solutions of such deterministic and wellposed problems^{24,25,26,27}; it also motivates the use of KH instability as a testbed numerical problem^{28}, and suggests the deterministic approach to unveil the physics of shear layer flows^{29,30}. The Hadamard viewpoint, however, appears quite at odds with Lorenz’s intuition. It also fails to account for the essential footprint of KH flows: their visually striking largescale features, that allow for identifying easily occurrences of KH instabilities across the physical communities, regardless of the detailed triggering mechanisms. This is the fundamental and longstanding observation that certain macroscopic statistical features of freely evolving shear flows are apparently only mildly dependent upon initial conditions^{31,32}.
Here we show that the classical Kelvin–Helmholtz singular shear layer becomes in fact wellposed and physically relevant when formulated in a probabilistic sense. Such a change from a deterministic to a probabilistic viewpoint reconciles Hadamardʼs to Lorenzʼs. To unveil the intrinsically random nature of the shear layer dynamics, we employ a probabilistic approach inspired from earlier studies connected to the theory of spontaneous stochasticity^{33,34,35,36,37,38} describing unpredictability of trajectories of Lagrangian particles advected by rough (nonsmooth) deterministic velocity fields, a key notion beneath the modern view on turbulent mixing^{39,40,41,42,43}. This approach makes use of microscale regularizing mechanisms, such as viscous dissipation and thermal fluctuations, to select relevant physical solutions at the macroscale. Inferring from this connection, our main result is to reveal the intrinsic spacetime randomness of the shear layer growth, characterized by the existence of a universal nontrivial probability distribution of the velocity field. Therefore, despite unpredictability, the resulting ideal flows possess welldefined statistical properties at finite times. These properties are triggered by, but not sensitive to the nature of infinitesimal microscale details.
Results
The singular shear layer
The classical mathematical formulation of the KH instability refers to perturbations of the idealized interface between two parallel streams with the velocity difference U. We will focus on twodimensional incompressible formulation, as it is more accessible for accurate numerical simulations; in applications, this formulation refers to largescale motions in the atmosphere and ocean^{44,45}. Such streams are defined by the constant horizontal velocities: \(U(x,y)=U/2\) in the upper halfplane \(y\, > \, 0\) and \(u(x,y)=U/2\) in the lower halfplane \(y\, < \, 0\); the vertical speed \(v(x,y)\) is everywhere zero. This flow can be seen as a vortex sheet localized on the horizontal axis \(y=0\), where the vorticity field \(\omega = \partial_x v  \partial_y u\) has the Diracdelta singularity, \(\omega=U\delta(y)\). This is a steady solution to the inviscid incompressible Euler equation. For the ideal flow, the wellknown linear theory^{14} predicts that a small perturbation with wavelength l has the exponential temporal growth \(\sim \exp (\pi Ut/l)\). The growth rate becomes infinite when the perturbation scale vanishes, implying the explosive breakdown of linear theory and signaling the illposedness of the underlying equations^{46}, as required for Lorenz’s intrinsic unpredictability. In order to select physically relevant solutions, we employ simultaneously two different smallscale mechanisms.
The first mechanism refers to a vanishingly small perturbation, which can be viewed as the effect of a small noise induced, e.g., by microscopic fluctuations^{7}. Specifically, we consider initial conditions for the vortex sheet with infinitesimal smallscale perturbations of the form
where \(\varepsilon\) is a small parameter and \(\eta(x)\) is a perturbation profile generated by white noise. The second mechanism desingularizes the equations of motion at vanishingly small scales, such that the resulting regularized equations can be used to evolve the flow. In order to explore the physical properties of the regularized dynamics, we focus on two fundamentally distinct formulations: (a) the viscous or hyperviscous (Navier–Stokes) regularization, which is controlled by the small viscosity parameter \(\nu \, > \, 0\), and (b) the pointvortex (Birkhoff–Rott) approximation, where the regularization is determined by the number of vortices \(N_b\) and controlled by the small parameter \(\nu \propto {N}_{b}^{1}\). Here, the first case employs a natural dissipative mechanism, and the flow is considered in the socalled Eulerian formulation. For deterministic initial data, adding a vanishingly small viscosity is known to select weak solutions of the 2D Euler equations but is not enough to ensure uniqueness^{47,48}. Differently, the second case follows the Lagrangian approach in fluid mechanics, also relevant for superfluids^{49,50}, in which one follows point vortices advected by the induced velocity field.
Recent numerical studies^{32} of the chaotic pointvortex system with finite \(N_b\) indicate a type of statistical universality of the nonlinear dynamics with respect to the initial data over a wide range of timescales. While finite amplitude of any random initial perturbation trivially entails statistical description, robustness of the macroscopic features with respect to the perturbation size points towards a deeper phenomenon. In this work we unveil the origin of dynamical universality, showing that for vanishingly small values of viscosity and noise, the KH instability develops into a universal and spontaneously stochastic flow.
By “spontaneous stochasticity”, we mean that such a flow is intrinsically random. “Universality” signifies that the statistical properties of the flow are independent of the type of regularization, except at very small scales and times. The mathematical formulation of this phenomenon corresponds to the vanishing viscosity and noise, \(\nu \to 0\) and \(\varepsilon \to 0\), where the two limits are taken simultaneously in a suitable fashion, such that physically meaningful norms of the initial disturbance (such as the kinetic energy) vanish; see the Methods section for more details.
This limiting procedure grants the originally illposed inviscid problem with a wellposed statistical interpretation where the asymptotic flow is unambiguously defined as a nontrivial stochastic process. Mathematically^{51,52}, one expects that each realization of this limiting process solves the incompressible Euler equations in a weak sense for the initial ideal vortex sheet, \(\omega = U\delta(y)\). Physically, this implies that infinitesimal effects of microscale regularization and noise do not select a unique deterministic KH flow, but rather a welldefined statistical solution, with universal nonstationary probability distribution.
The spontaneous butterfly effect
In order to reveal how the intrinsic randomness in the vortex layer emerges at small times, we first focus on the Navier–Stokes dynamics, with the viscous parameter \(\nu\) here also specifying the noise amplitude as \(\varepsilon \propto \nu\). We estimate the random component of the flow by measuring how two solutions, \((u,v)\) and \((u^{\prime} ,v^{\prime} )\), which differ initially by a very small random perturbation of size \(\varepsilon\), split in time. This is done by introducing the separation energy
where \(L\) is a prescribed (macroscopic) observation scale and the average \(\langle \cdot \rangle_{x,\eta}\) is over both the xdirection and the statistical realizations of the initial noise. With this choice of normalization, \({\mathcal{E}}\) varies between 0 and 1, and those two extremal values are reached when the two solutions, respectively, coincide or fully decorrelate over the observational window. The results are presented in Fig. 1, where the separation energy is plotted for the three different sets of parameters \(\nu\) and \(\varepsilon\).
Initially the flows are very close and diverge exponentially as seen in the inset of Fig. 1. This exponential growth is governed by a positive Lyapunov exponent, which is the distinctive feature of the usual butterfly effect^{3}. However, it is apparent that the growth rate depends strongly on the viscous parameter \(\nu\). In fact, the main part of Fig. 1 demonstrates that the exponential growth is a transient behavior. The transient region is shifted to smaller and smaller times in the combined limit of infinitesimal viscosity and initial perturbation, \(\nu \to 0\) and \(\varepsilon \to 0\). At larger times, the separation energy reaches the asymptotics \({\mathcal{E}}\approx 0.14t\), which is independent of the initial separation and has the order of the total energy within the flow. This behavior quantifies spontaneous stochasticity for infinitesimal viscosity and initial disturbance: solutions which are initially undistinguishable and deterministic become distinct and random at finite times. It is analogous to Richardson’s relative dispersion in fully developed turbulent flow^{53}, where the distance between two Lagrangian fluid elements grows algebraically as \(\sim \!t^{3/2}\), independently of how close they initially are. The explosive separation of Fig. 1 can be seen as the Eulerian counterpart to this phenomenon as it pertains to the full description of the flow. This type of randomness is spontaneous since it builds up in infinitesimal time and requires only infinitesimal perturbations to be revealed.
Time development of the vortex layer
We will now argue that while the growth of the vortex layer is an illposed dynamics from the deterministic viewpoint, it becomes wellposed in a probabilistic sense, meaning that macroscopic statistical features of the flow are independent of the noise and of the type of regularization. Fig. 2 shows a typical vorticity distribution for numerical simulations with the two regularization methods; see the Methods section for details of numerical implementations. Panels (a) and (c) correspond to the viscous regularization based on direct numerical simulations of the Navier–Stokes equations in a square domain with periodic boundary conditions. Panels (b) and (d) correspond to simulations of the Birkhoff–Rott dynamics. Visually, the vorticity distributions obtained for each type of regularization look akin to each other on macroscopic scales. It is only upon zooming into the microscopic structure, that the distinction between the viscous and the pointvortex regularization becomes apparent, as the first is continuous and the second is discrete.
On macroscopic scales, the development of the shear layer consists in a cascade process of collisions and the subsequent mergers of smaller vortex blobs into larger ones; see Fig. 2. Here, the blobs are small concentrated regions of high vorticity, surrounded by halos of lowervorticity streaks^{28,31,54,55}. In the collision process, a part of the vorticity is scattered away from the blobs, and it may be absorbed by the same or other blobs at later times. Measurements suggests that the vorticity is divided into two approximately equal parts corresponding to the blobs and to the background flow, while the enstrophy is mostly carried by the blobs.
We now verify numerically that the vortex layer can be qualified as being universal and spontaneously stochastic, namely that the increase of numerical resolution along with decreasing viscosity and noise yield a unique stochastic solution at finite times. By the scaling theory^{56}, one may expect that statistical properties of the macroscopic dynamics depend only on the two independent quantities, namely, the velocity jump\(U\) and the mixing length \(\ell(t)\). The latter describes the width of the vortex layer, and can be conveniently defined as
where \({\bf{R}} = (x,y)\). The average \(\langle \cdot\rangle_{x,\eta}\) is over both the xdirection and the statistical realizations of the initial noise. For the Birkhoff–Rott regularization, this quantity represents the standarddeviation of pointvortices across the vortex layer.
From dimensional analysis, we expect that the mixing layer grows linearly in time: \(\ell = \alpha U t\), where \(\alpha\) is a dimensionless coefficient. Figure 3 presents the numerical results shown in log–log scale, which verify this asymptotic linear scaling law. The vortex layer is formed after a short transient time, and this transient time decreases as we decrease the regularization parameter and noise. The subsequent evolution yields the universal prefactor \(\alpha\simeq 0.029\) for both Navier–Stokes and Birkhoff–Rott regularizations; see the inset in Fig. 3.
Statistical universality
We now study multiscale statistical properties of the flow. We start by analyzing the evolution of the normalized vorticity profile,
obtained from the vorticity distribution averaged with respect to the xdirection and the statistical realizations. In the Birkhoff–Rott case, \(p_1\) relates to the onepoint distribution of pointvortices along the ydirection so that\(\int p_1(y,t) \, dy=\ell(t)\). Numerical vorticity profiles are shown in Fig. 4. Panel (a) highlights the selfsimilarity of this onepoint statistical quantity, that is \({p}_{1}(y,t)\approx {P}_{1}\left(y/\ell (t)\right)\). Panel (b) confirms that the asymptotic \(P_1\), corresponding to vanishing amplitudes of the initial perturbations, is identical for both the Navier–Stokes and Birkhoff–Rott regularizations.
As a secondorder statistical observable, we introduce the isotropic twopoint correlator
written in dimensionless form. It is constructed from the covariance of the vorticity field between the points \({\bf{R}}=(x,y)\) and \({\bf{R}} + {\boldsymbol{\rho}}=(x+\rho_x,y+\rho_y)\) at a given distance \(\left\Vert {\boldsymbol{\rho }}\right\Vert =r\). The function \(p_2(r,t)\) is strongly dominated by the relative positions of concentrated blobs of vorticity (see Fig. 2), for which the integrated product of vorticities in Eq. (4) is maximal.
Results of numerical simulations are shown in Fig. 5, where Panel (a) displays the asymptotic selfsimilar relation \({p}_{2}(r,t)\approx {P}_{2}\left(r/\ell (t)\right)\). Panel (b) demonstrates that the asymptotic form of \(p_2\) is universal: at a given finite time, it collapses onto a single universal function \(P_2\) in the limit of vanishing regularization and noise for both Navier–Stokes and Birkhoff–Rott formulations. One can see that the profiles for different regularizations are distinct only at small scales, which reflects the different smallscale structures visible in the bottom plots of Fig. 2. The function \(P_2\) has a pronounced maximum around \(r=1.4\ell(t)\) characterizing a typical distance between nearest vortex blobs in Fig. 2. For very large distances \(r \gg \ell(t)\), it converges to the limiting value \(P_2 \to 2\). This value can be computed by substituting the vortex sheet expression \(\omega=U\delta(y)\) into Eq. (4), which approximates the vortex layer at large scales. It is remarkable that a constant state develops also at small distances \(r \ll \ell(t)\), with the universal value estimated as \(P_2 \to 2.4\); see Fig. 5. This value characterizes the uniform statistical distribution of the vorticity at distances, which are small compared to the mixing length, but still larger than the regularization scale. In our view, this constant asymptotic value reflects the selfsimilar nature of the reconnection process, where vortex blobs constantly attract each other, before merging into larger blobs.
Discussion
Opposed to the exponential growth of observational errors in chaotic systems, now known as the butterfly effect, Lorenz anticipated that finitetime randomness can emerge in some formally deterministic multiscale fluid systems irrespective of how small the initial uncertainty^{2}. In our work, we demonstrated that this scenario is relevant for the flows generated by the Kelvin–Helmholtz instability of a singular shear layer. Being illposed for the ideal fluid mechanics, this system triggers a spontaneous stochastic process for the vorticity field, when considered in the combined limit of vanishing microscale regularization and vanishing random perturbation. In particular, considering hyperviscous Navier–Stokes vs. pointvortex Birkhoff–Rott formulations, we showed that this limiting stochastic process is universal.
The recently developed theory for particles advected by nonsmooth deterministic velocity fields^{33,34,35,36,37} demonstrate similar spontaneously stochastic properties when infinitesimal noise is introduced: fluid parcels separate in finite time no matter how close they are initially. Such an approach proved to be relevant for particles in turbulent flows^{57,58,59,60,61}, where it is referred to as the Richardson superdiffusion^{53}. Within this modern perception, Lorenz’s scenario of intrinsic randomness can be seen as the Eulerian counterpart of spontaneous stochasticity applied to the flow fields rather than to individual particles^{11}. In our point of view, this extension is rather natural, because the Eulerian description of fluid dynamics involves the transport of active quantities, such as momentum or vorticity. A comprehensive mathematical framework of this effect dealing with statistical solutions to fluid equations^{62,63}, however, is still to be developed^{51}.
Whether the finitetime unpredictability is caused by a small though finite “butterfly” as in chaos or by an arbitrarily small perturbation as in spontaneous stochasticity, the distinction may be illusive for an unprepared observer. However, these two phenomena represent fundamentally different physical and mathematical mechanisms that bridge the deterministic classical world to genuinely random systems. As such, spontaneous stochasticity could then possibly provide the salient mechanism that generates intrinsic randomness in natural phenomena. It provides a conceptual framework that justifies the relevance of stochastic modeling for apparently deterministic fluid dynamical systems^{64}, as well as for systematic statistical mechanics approaches^{32}. In addition to potential applications in fluid mechanics, which of course include the open problem of developed turbulence^{65} or the physics of boundary layers^{43,66}, the relevance of the spontaneous stochasticity phenomenon may extend to other nonlinear field theories featuring multiscale dynamics and singularities, e.g., wave turbulence^{67}, nonlinear optics^{68}, and astrophysics^{69,70}.
Methods
Navier–Stokes regularization
Our simulations were performed using the opensource GHOST parallel solver^{71,72}, that employs a pseudospectral scheme with standard 2/3 dealiasing on spatial grids with \({N}_{a}^{2}\) points, and secondorder Runge–Kutta scheme in time. The timestep is fixed and characterized by Courant number \(\sim 0.2\). Specifically, the Navier–Stokes equations in vorticity form
are integrated within a twodimensional periodic domain \((x,y) \in [0,2\pi]^2\) Results presented in the paper are obtained using a hyperviscous linear dissipation defined through its Fourierspace representation as \({\widehat{D}}_{\nu }=\nu {k}^{2}{(k/{k}_{\max })}^{6}\widehat{\omega }\), where the regularization wavenumber \({k}_{\max }={N}_{a}/3\) is used; this is a common tool to localize viscous effects at small scales. Simulations with the standard viscous dissipation \(D_\nu = \nu \Delta \omega\) were also carried out employing up to 16, 384^{2} grid points: they lead to the same conclusions though for a shorter interval of scales due to numerical limitations.
To comply with periodicity, the initial conditions Eq. (1) are implemented along two parallel lines: \(y=0\) and \(y=\pi\) so that the total vorticity is zero. The noise \(\eta(x)\) is generated by a centered uniform stochastic process with spatial correlations \(\left\langle \eta (x)\eta ({x}^{\prime})\right\rangle =\delta (x{x}^{\prime})/3\). Simulations are stopped at sufficiently small times to avoid interaction between the layers and nonlocal effects of a finite domain. Then, the two resulting vortex layers are analyzed separately in the respective domains.
We performed nine simulations for each of the three parameter sets from Table 1. These parameters follow the rule \({k}_{\max }\propto {N}_{a}\), \(\nu \propto {N}_{a}^{3/2}\) and \(\varepsilon \propto {N}_{a}^{3/2}\), so that \(\nu \to 0\) and \(\varepsilon\to 0\) simultaneously with the increasing resolution \(N_a\). To avoid spurious Gibbs oscillations of rough initial conditions, the initial state is diffused by hyperviscosity until the largestwavenumber \({k}_{\max }\) has negligeable Fourier contribution, comparable to the machine precision. This takes a short initialization time \({t}_{d}\propto {N}_{a}^{1/2}\), vanishing with increasing resolution. After this initial smoothing, the energy and enstrophy (squared norm of the vorticity) of the initial perturbation in expression Eq. (1) scale as \({\mathcal{E}}({0}^{+}) \sim {N}_{a}^{3}{\mathrm{log}}\,{N}_{a}\) and \({\mathcal{Z}}({0}^{+}) \sim {N}_{a}^{1}\): Both quantities vanish upon increasing resolution. The initial (regularized) conditions therefore converge to the straight (unperturbed) vortex sheet in both energy and enstrophy metrics. We observe that our choice of scaling guarantees that the most unstable linear mode in our simulation is \(\propto N_a\) upon increasing resolution.
Finally, please note that choosing simultaneous scaling rules for \(\varepsilon \to 0\) and \(\nu \to 0\) is dictated by the fact that these limits in principle do not commute. In particular, letting \(\varepsilon \to 0\) before \(\nu \to 0\), one should expect the process to converge towards the deterministic trivial baseflow.
Birkhoff–Rott regularization
In this case, we use a domain that is 1periodic in the xdirection and unlimited in the ydirection. The initial condition Eq. (1) is approximated by a periodic discrete row of pointvortices located on the line \(y=0\) at positions \(x_n(0)=n/N_b, n=1,\cdots , {N_{b}}\), and carrying the initial vorticity \({\omega }_{n}=\frac{U}{{N}_{b}}\delta (x{x}_{n})\delta (y)(1+\varepsilon \ {N}_{b}^{1/2}\ {\eta }_{n})\), where the \(\eta_n\)’s are independent Gaussian variables with unit variance. Using the complex position variables \(z_n=x_n+iy_n\), Biot–Savart law prescribes the advection of each pointvortex as^{54,73}
where the star subscript denotes complex conjugation. We performed simulations of the dynamics Eq. (6) with \(U=1,\ \varepsilon=10^{5}\), and \(N_b=2^{15}, 2^{16}\) and \(2^{17}\) vortices using the fourthorder Runge–Kutta scheme. The timestep \(\Delta t\) is chosen adaptively to satisfy \(\Delta t\,<\,0.1\ \mathop{\min }_{n\ne m}\left\{ {z}_{n}{z}_{m} / {\dot{z}}_{n}{\dot{z}}_{m} \right\}\). This sort of Courant–Friedrichs–Lewy upperbound is enough to ensure stability and accuracy of the numerical scheme. Indeed, with such a choice, the kinetic energy, defined as
is conserved with a highenough accuracy. At the end of the simulations, the relative error in energy remains below 10^{−4}.
Qualitatively similar CFL criteria were employed in previous simulations of pointvortex systems in similar setups^{32,74}. While this CFL criterion in principle does not preclude divergence of dynamical trajectories due to numerical noise, it is enough to ensure convergence in a statistical sense. Upon increasing resolutions, the spacetime distributions of pointvortices indeed converge towards a nontrivial measure. In mathematical terms, the observed statistical universality reflects convergence in distribution^{75}.
Finitesize effects
Formally, the scaling theory for the singular vortex sheet addresses flows evolving within an unbounded domain. In our numerical simulations, this domain has a finite size in the xdirection, taken as \(L=2\pi\) for the Navier–Stokes and \(L=1\) for the Birkhoff–Rott models. In all figures, we displayed the results in dimensionless form using \(t_*=L/U\) as the time unit, for which the mixing layer invades roughly 3% of the computational domain.
Data availability
The data that support the findings of this study are available from the corresponding author on request.
Code availability
The simulation and postprocessing codes that have been used to produce the results of this study are available from the corresponding author on request.
References
 1.
Lorenz, E. N. Deterministic nonperiodic flow. J. Atmos. Sci. 20, 130–141 (1963).
 2.
Lorenz, E. N. The predictability of a flow which possesses many scales of motion. Tellus 21, 289–307 (1969).
 3.
Gleick, J. Chaos: Making a new science (Open Road Media, New York, 2011).
 4.
Boffetta, G., Cencini, M., Falcioni, M. & Vulpiani, A. Predictability: a way to characterize complexity. Phys. Rep. 356, 367–474 (2002).
 5.
Palmer, T. N., Döring, A. & Seregin, G. The real butterfly effect. Nonlinearity 27, R123 (2014).
 6.
Leith, C. E. & Kraichnan, R. H. Predictability of turbulent flows. J. Atmos. Sci 29, 1041–1058 (1972).
 7.
Ruelle, D. Microscopic fluctuations and turbulence. Phys. Lett. A 72, 81–82 (1979).
 8.
Eyink, G. L. Turbulence noise. J. Stat. Phys. 83, 955–1019 (1996).
 9.
Palmer, T. N. Predicting uncertainty in forecasts of weather and climate. Rep. Prog. Phys. 63, 71 (2000).
 10.
Boffetta, G. & Musacchio, S. Predictability of the inverse energy cascade in 2D turbulence. Phys. Fluids 13, 1060–1062 (2001).
 11.
Mailybaev, A. A. Spontaneously stochastic solutions in onedimensional inviscid systems. Nonlinearity 29, 2238 (2016).
 12.
Mailybaev, A. A. Toward analytic theory of the RayleighTaylor instability: lessons from a toy model. Nonlinearity 30, 2466 (2017).
 13.
Biferale, L., Boffetta, G., Mailybaev, A. A. & Scagliarini, A. RayleighTaylor turbulence with singular nonuniform initial conditions. Phys. Rev. Fluids 3, 092601 (2018).
 14.
Drazin, P. G. Introduction to Hydrodynamic Stability, Vol. 32 (Cambridge University Press, Cambridge, 2002).
 15.
Matsuoka, C. KelvinHelmholtz Instability and Rollup. Scholarpedia 9, 11821 (2014).
 16.
Blaauwgeers, R. et al. Shear flow and KelvinHelmholtz instability in superfluids. Phys. Rev. Lett. 89, 155301 (2002).
 17.
Takeuchi, H., Suzuki, N., Kasamatsu, K., Saito, H. & Tsubota, M. Quantum KelvinHelmholtz instability in phaseseparated twocomponent BoseEinstein condensates. Phys. Rev. B 81, 094517 (2010).
 18.
Smyth, W. D. & Moum, J. N. Ocean mixing by KelvinHelmholtz instability. Oceanography 25, 140–149 (2012).
 19.
Wang, C.Y. & Chevalier, R. A. Instabilities and clumping in type IA supernova remnants. Astrophys. J. 549, 1119 (2001).
 20.
Vietri, M., Ferrara, A. & Miniati, F. The survival of interstellar clouds against KelvinHelmholtz instabilities. Astrophys. J. 483, 262 (1997).
 21.
Székelyhidi, L. Weak solutions to the incompressible Euler equations with vortex sheet initial data. C. R. Math. 349, 1063–1066 (2011).
 22.
Hadamard, J. Sur les problèmes aux dérivées partielles et leur signification physique (Princeton university bulletin, 1902).
 23.
Sobolev, S. Partial Differential Equations of Mathematical Physics (Elsevier, Amsterdam, 2016).
 24.
Caflisch, R., Orellana, O. & Siegel, M. Singularity formulation and illposedness for vortex sheets. SIAM J. Math. Anal. 20, 293–307 (1989).
 25.
Duchon, J. & Robert, R. Global vortex sheet solutions of Euler equations in the plane. J. Diff. Eq. 73, 215–224 (1988).
 26.
Cowley, S., Baker, G. & Tanveer, S. On the formation of moore curvature singularities in vortex sheets. J. Fluid Mech. 378, 233–267 (1999).
 27.
Sulem, C., Sulem, P., Bardos, C. & Frisch, U. Finite time analyticity for the two and three dimensional kelvinhelmholtz instability. Com. Math. Phys 80, 485–516 (1981).
 28.
Lecoanet, D. et al. A validated nonlinear KelvinHelmholtz benchmark for numerical hydrodynamics. Mon. Not. R. Astron. Soc 455, 4274–4288 (2015).
 29.
Corcos, G. in The Role of Coherent Structures in Modelling Turbulence and Mixing (Springer, New York, 1981) pp. 10–40.
 30.
Corcos, G. & Sherman, F. The mixing layer: deterministic models of a turbulent flow. part 1. introduction and the twodimensional flow. J. Fluid Mech. 139, 29–65 (1984).
 31.
Sommeria, J., Staquet, C. & Robert, R. Final equilibrium state of a twodimensional shear layer. J. Fluid Mech. 233, 661–689 (1991).
 32.
Suryanarayanan, S., Narasimha, R. & Dass, N. Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics. Phys. Rev. E 89, 013009 (2014).
 33.
VandenEijnden, E. Generalized flows, intrinsic stochasticity, and turbulent transport. Proc. Natl Acad. Sci. USA 97, 8200–8205 (2000).
 34.
Falkovich, G., Gawȩdzki, K. & Vergassola, M. Particles and fields in fluid turbulence. Rev. Mod. Phys. 73, 913 (2001).
 35.
Gawedzki, K. in Intermittency in Turbulent Flows (ed Vassilicos, J. C.) (Cambridge University Press, Cambridge, 2001) pp. 86–104.
 36.
LeJan, Y. & Raimond, O. Integration of Brownian vector fields. Ann. Probab. 30, 826–873 (2002).
 37.
Kupiainen, A. Nondeterministic dynamics and turbulent transport. Ann. Henri Poincaré 4, 713–726 (2003).
 38.
Eyink, G. et al. Fluxfreezing breakdown in highconductivity magnetohydrodynamic turbulence. Nature 497, 466 (2013).
 39.
Bernard, D., Gawedzki, K. & Kupiainen, A. Slow modes in passive advection. J. Stat. Phys. 90, 519–569 (1998).
 40.
Shraiman, B. I. & Siggia, E. D. Scalar turbulence. Nature 405, 639 (2000).
 41.
Iyer, K. P., Schumacher, J., Sreenivasan, K. R. & Yeung, P. K. Steep cliffs and saturated exponents in threedimensional scalar turbulence. Phys. Rev. Lett. 121, 264501 (2018).
 42.
Dombre, T. Zero modes, instantons and fluctuations in the kraichnan model. Theory Exp. 2018, 113202 (2018).
 43.
Drivas, T. D. & Eyink, G. L. A lagrangian fluctuationdissipation relation for scalar turbulence. part i. flows with no bounding walls. J. Fluid Mech. 829, 153–189 (2017).
 44.
Boffetta, G. & Ecke, R. E. Twodimensional turbulence. Annu. Rev. Fluid Mech. 44, 427–451 (2012).
 45.
Bouchet, F. & Venaille, A. Statistical mechanics of twodimensional and geophysical flows. Phys. Rep. 515, 227–295 (2012).
 46.
Majda, A. J. & Bertozzi, A. L. Vorticity and Incompressible Flow (Cambridge University Press, Cambridge, 2002).
 47.
Delort, J.M. Existence de nappes de tourbillon en dimension deux. J. Am. Math. Soc. 4, 553–586 (1991).
 48.
Majda, A. J. Remarks on weak solutions for vortex sheets with a distinguished sign. Ind. Univ. Math. J. 42, 921–939 (1993).
 49.
Johnstone, S. et al. Evolution of largescale flow from turbulence in a twodimensional superfluid. Science 364, 1267–1271 (2019).
 50.
Gauthier, G. et al. Giant vortex clusters in a twodimensional quantum fluid. Science 364, 1264–1267 (2019).
 51.
Flandoli, F. Random Perturbation of PDEs and Fluid Dynamic Models: École d’été de Probabilités de SaintFlour XL–2010, Vol. 2015 (Springer Science, Business Media, New York, 2011).
 52.
Drivas, T. D., Mailybaev, A. A., and Raibekas, A., Statistical determinism in nonlipschitz dynamical systems, preprint: http://arXiv.org/abs/2004.03075 (2020).
 53.
Frisch, U. Turbulence: the legacy of AN Kolmogorov (Cambridge University Press, Cambridge, 1995).
 54.
Krasny, R. Desingularization of periodic vortex sheet rollup. J. Comput. Phys. 65, 292–313 (1986).
 55.
Lesieur, M., Staquet, C., LeRoy, P. & Comte, P. The mixing layer and its coherence examined from the point of view of twodimensional turbulence. J. Fluid Mech. 192, 511–534 (1988).
 56.
Barenblatt, G. I. Scaling, Vol. 34 (Cambridge University Press, Cambridge, 2003).
 57.
Jullien, M.C., Paret, J. & Tabeling, P. Richardson pair dispersion in twodimensional turbulence. Phys. Rev. Lett. 82, 2872 (1999).
 58.
Scatamacchia, R., Biferale, L. & Toschi, F. Extreme events in the dispersions of two neighboring particles under the influence of fluid turbulence. Phys. Rev. Lett. 109, 144501 (2012).
 59.
Bitane, R., Homann, H. & Bec, J. Geometry and violent events in turbulent pair dispersion. J. Turbul. 14, 23–45 (2013).
 60.
Bitane, R., Homann, H. & Bec, J. Time scales of turbulent relative dispersion. Phys. Rev. E 86, 045302 (2012).
 61.
Buaria, D., Sawford, B. L. & Yeung, P.K. Characteristics of backward and forward twoparticle relative dispersion in turbulence at different Reynolds numbers. Phys. Fluids 27, 105101 (2015).
 62.
Fjordholm, U., Mishra, S. & Tadmor, E. On the computation of measurevalued solutions. Act. Num 25, 567–679 (2016).
 63.
Fjordholm, U., Lanthaler, S. & Mishra, S. Statistical solutions of hyperbolic conservation laws: foundations. Arch. Rat. Mech. An 226, 809–849 (2017).
 64.
Palmer, T. N. Stochastic weather and climate models. Nat. Rev. Phys 1, 463–471 (2019).
 65.
Brenner, M. P., Hormoz, S. & Pumir, A. Potential singularity mechanism for the Euler equations. Phys. Rev. Fluids 1, 084503 (2016).
 66.
Waidmann, M., Klein, R., Farge, M. & Schneider, K. Energy dissipation caused by boundary layer instability at vanishing viscosity. J. Fluid Mech. 849, 676–717 (2018).
 67.
Nazarenko, S. Wave turbulence, Vol. 825 (Springer Science & Business Media, New York, 2011).
 68.
Kraych, A. E., Agafontsev, D., Randoux, S. & Suret, P. Statistical properties of the nonlinear stage of modulation instability in fiber optics. Phys. Rev. Lett. 123, 093902 (2019).
 69.
Genel, S. et al. A quantification of the butterfly effect in cosmological simulations and implications for galaxy scaling relations. Astrophys. J. 871, 21 (2019).
 70.
Khalatnikov, I. M., Lifshitz, E. M., Khanin, K. M., Shchur, L. N. & Sinai, Y. G. On the stochasticity in relativistic cosmology. J. Stat. Phys. 38, 97–114 (1985).
 71.
Gómez, D. O., Mininni, P. D. & Dmitruk, P. Parallel simulations in turbulent MHD. Phys. Scr. 2005, 123 (2005).
 72.
Mininni, P. D., Rosenberg, D., Reddy, R. & Pouquet, A. A hybrid MPIOpenMP scheme for scalable parallel pseudospectral computations for fluid turbulence. Parallel Comput. 37, 316–326 (2011).
 73.
Aref, H. Point vortex dynamics: a classical mathematics playground. J. Math. Phys. 48, 065401 (2007).
 74.
Aref, H. & Siggia, E. Vortex dynamics of the twodimensional turbulent shear layer. J. Fluid Mech. 100, 705–737 (1980).
 75.
Jacod, J. & Protter, P. Probability Essentials (Springer Science, Business Media, New York, 2012).
Acknowledgements
We thank M. Brachet and P. Mininni for insightful discussions and G. Eyink for sharing his continuously encouraging perspectives on the phenomenon of spontaneous stochasticity. This work was supported by CNPq (grants 303047/20186, 406431/20183) and the BrazilianFrench Network in Mathematics.
Author information
Affiliations
Contributions
S.T., J.B., and A.A.M. contributed equally to the implementation, the analysis, and the redaction of the present study.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Thalabard, S., Bec, J. & Mailybaev, A.A. From the butterfly effect to spontaneous stochasticity in singular shear flows. Commun Phys 3, 122 (2020). https://doi.org/10.1038/s4200502003916
Received:
Accepted:
Published:
Further reading

Peering inside a cough or sneeze to explain enhanced airborne transmission under dry weather
Scientific Reports (2021)

Solvable Intermittent Shell Model of Turbulence
Communications in Mathematical Physics (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.