BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Freddy Bouchet

Freddy Bouchet appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechphysics.flu-dynphysics.ao-phnlin.CDmath-phmath.MPphysics.comp-phphysics.data-anphysics.class-phphysics.plasm-phastro-ph.HEComputationMachine Learningmath.DSmath.NA

Catalog footprint

What is connected

36works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

Coupling rare event algorithms with data-based learned committor functions using the analogue Markov chain

Rare events play a crucial role in many physics, chemistry, and biology phenomena, when they change the structure of the system, for instance in the case of multistability, or when they have a huge impact. Rare event algorithms have been devised to simulate them efficiently, avoiding the computation of long periods of typical fluctuations. We consider here the family of splitting or cloning algorithms, which are versatile and specifically suited for far-from-equilibrium dynamics. To be efficient, these algorithms need to use a smart score function during the selection stage. Committor functions are the optimal score functions. In this work we propose a new approach, based on the analogue Markov chain, for a data-based learning of approximate committor functions. We demonstrate that such learned committor functions are extremely efficient score functions when used with the Adaptive Multilevel Splitting algorithm. We illustrate our approach for a gradient dynamics in a three-well potential, and for the Charney-DeVore model, which is a paradigmatic toy model of multistability for atmospheric dynamics. For these two dynamics, we show that having observed a few transitions is enough to have a very efficient data-based score function for the rare event algorithm. This new approach is promising for use for complex dynamics: the rare events can be simulated with a minimal prior knowledge and the results are much more precise than those obtained with a user-designed score function.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deep Learning-based Extreme Heatwave Forecast

Because of the impact of extreme heat waves and heat domes on society and biodiversity, their study is a key challenge. We specifically study long-lasting extreme heat waves, which are among the most important for climate impacts. Physics driven weather forecast systems or climate models can be used to forecast their occurrence or predict their probability. The present work explores the use of deep learning architectures, trained using outputs of a climate model, as an alternative strategy to forecast the occurrence of extreme long-lasting heatwaves. This new approach will be useful for several key scientific goals which include the study of climate model statistics, building a quantitative proxy for resampling rare events in climate models, study the impact of climate change, and should eventually be useful for forecasting. Fulfilling these important goals implies addressing issues such as class-size imbalance that is intrinsically associated with rare event prediction, assessing the potential benefits of transfer learning to address the nested nature of extreme events (naturally included in less extreme ones). We train a Convolutional Neural Network, using 1000 years of climate model outputs, with large-class undersampling and transfer learning. From the observed snapshots of the surface temperature and the 500 hPa geopotential height fields, the trained network achieves significant performance in forecasting the occurrence of long-lasting extreme heatwaves. We are able to predict them at three different levels of intensity, and as early as 15 days ahead of the start of the event (30 days ahead of the end of the event).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Instantons and the path to intermittency in turbulent flows

Processes leading to anomalous fluctuations in turbulent flows, referred to as intermittency, are still challenging. We consider cascade trajectories through scales as realizations of a stochastic Langevin process for which multiplicative noise is an intrinsic feature of the turbulent state. The trajectories are conditioned on their entropy exchange. Such selected trajectories concentrate around an optimal path, called instanton, which is the minimum of an effective action. The action is derived from the Langevin equation, estimated from measured data. In particular instantons with negative entropy pinpoint the trajectories responsible for the emergence of non-Gaussian statistics at small-scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Path large deviations for the kinetic theory of weak turbulence

We consider a generic Hamiltonian system of nonlinear interacting waves with 3-wave interactions. In the kinetic regime of wave turbulence, which assumes weak nonlinearity and large system size, the relevant observable associated with the wave amplitude is the empirical spectral density that appears as the natural precursor of the spectral density, or spectrum, for finite system size. Following classical derivations of the Peierls equation for the moment generating function of the wave amplitudes in the kinetic regime, we propose a large deviation estimate for the dynamics of the empirical spectral density, where the number of admissible wavenumbers, which is proportional to the volume of the system, appears as the natural large deviation parameter. The large deviation stochastic Hamiltonian that quantifies the minus of the log probability of a trajectory is computed within the kinetic regime which assumes the Random Phase approximation for weak nonlinearity. We compare this Hamiltonian with the one for a system of modes interacting in a mean-field way with the empirical spectrum. Its relationship with the Random Phase and Amplitude approximation is discussed. Moreover, for the specific case when no forces and dissipation are present, a few fundamental properties of the large deviation dynamics are investigated. We show that the latter conserves total energy and momentum, as expected for a 3-wave interacting systems. In addition, we compute the equilibrium quasipotential and check that global detailed balance is satisfied at the large deviation level. Finally, we discuss briefly some physical applications of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamical large deviations for homogeneous systems with long range interactions and the Balescu--Guernsey--Lenard equation

We establish a large deviation principle for time dependent trajectories (paths) of the empirical density of $N$ particles with long range interactions, for homogeneous systems. This result extends the classical kinetic theory that leads to the Balescu--Guernsey--Lenard kinetic equation, by the explicit computation of the probability of typical and large fluctuations. The large deviation principle for the paths of the empirical density is obtained through explicit computations of a large deviation Hamiltonian. This Hamiltonian encodes all the cumulants for the fluctuations of the empirical density, after time averaging of the fast fluctuations. It satisfies a time reversal symmetry, related to the detailed balance for the stochastic process of the empirical density. This explains in a very simple way the increase of the macrostate entropy for the most probable states, while the stochastic process is time reversible, and describes the complete stochastic process at the level of large deviations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computing return times or return periods with rare event algorithms

The average time between two occurrences of the same event, referred to as its return time (or return period), is a useful statistical concept for practical applications. For instance insurances or public agency may be interested by the return time of a 10m flood of the Seine river in Paris. However, due to their scarcity, reliably estimating return times for rare events is very difficult using either observational data or direct numerical simulations. For rare events, an estimator for return times can be built from the extrema of the observable on trajectory blocks. Here, we show that this estimator can be improved to remain accurate for return times of the order of the block size. More importantly, we show that this approach can be generalised to estimate return times from numerical algorithms specifically designed to sample rare events. So far those algorithms often compute probabilities, rather than return times. The approach we propose provides a computationally extremely efficient way to estimate numerically the return times of rare events for a dynamical system, gaining several orders of magnitude of computational costs. We illustrate the method on two kinds of observables, instantaneous and time-averaged, using two different rare event algorithms, for a simple stochastic process, the Ornstein-Uhlenbeck process. As an example of realistic applications to complex systems, we finally discuss extreme values of the drag on an object in a turbulent flow.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Instantons for the destabilization of the inner Solar System

For rare events, path probabilities often concentrate close to a predictable path, called instanton. First developed in statistical physics and field theory, instantons are action minimizers in a path integral representation. For chaotic deterministic systems, where no such action is known, shall we expect path probabilities to concentrate close to an instanton? We address this question for the dynamics of the terrestrial bodies of the Solar System. It is known that the destabilization of the inner Solar System might occur with a low probability, within a few hundred million years, or billion years, through a resonance between the motions of Mercury and Jupiter perihelia. In a simple deterministic model of Mercury dynamics, we show that the first exit time of such a resonance can be computed. We predict the related instanton and demonstrate that path probabilities actually concentrate close to this instanton, for events which occur within a few hundred million years. We discuss the possible implications for the actual Solar System.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Is the Boltzmann equation reversible? A large deviation perspective on the irreversibility paradox

We consider the kinetic theory of dilute gases in the Boltzmann--Grad limit. We propose a new perspective based on a large deviation estimate for the probability of the empirical distribution dynamics. Assuming Boltzmann molecular chaos hypothesis (Stosszahlansatz), we derive a large deviation rate function, or action, that describes the stochastic process for the empirical distribution. The quasipotential for this action is the negative of the entropy, as should be expected. While the Boltzmann equation appears as the most probable evolution, corresponding to a law of large numbers, the action describes a genuine reversible stochastic process for the empirical distribution, in agreement with the microscopic reversibility. As a consequence, this large deviation perspective gives the expected meaning to the Boltzmann equation and explains its irreversibility as the natural consequence of limiting the physical description to the most probable evolution. More interestingly, it also quantifies the probability of any dynamical evolution departing from solutions of the Boltzmann equation. This picture is fully compatible with the heuristic classical view of irreversibility, but makes it much more precise in various ways. We also explain that this large deviation action provides a natural gradient structure for the Boltzmann equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical study of extreme mechanical force exerted by a turbulent flow on a bluff body by direct and rare-event sampling techniques

This study investigates, by means of numerical simulations, extreme mechanical force exerted by a turbulent flow impinging on a bluff body, and examines the relevance of two distinct rare-event algorithms to efficiently sample these events. The drag experienced by a square obstacle placed in a turbulent channel flow (in two dimensions) is taken as a representative case study. Direct sampling shows that extreme fluctuations are closely related to the presence of a strong vortex blocked in the near wake of the obstacle. This vortex is responsible for a significant pressure drop between the forebody and the base of the obstacle, thus yielding a very high value of the drag. Two algorithms are then considered to speed up the sampling of such flow scenarii, namely the AMS and the GKTL algorithms. The general idea behind these algorithms is to replace a long simulation by a set of much shorter ones, running in parallel, with dynamics that are replicated or pruned, according to some specific rules designed to sample large-amplitude events more frequently. These algorithms have been shown to be relevant for a wide range of problems in statistical physics, computer science, biochemistry. The present study is the first application to a fluid-structure interaction problem. Practical evidence is given that the fast sweeping time of turbulent fluid structures past the obstacle has a strong influence on the efficiency of the rare-event algorithm. While the AMS algorithm does not yield significant run-time savings as compared to direct sampling, the GKTL algorithm appears to be effective to sample very efficiently extreme fluctuations of the time-averaged drag and estimate related statistics such as return times.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Atmospheric bistability and abrupt transitions to superrotation: wave-jet resonance and Hadley cell feedbacks

Strong eastward jets at the equator have been observed in many planetary atmospheres and simulated in numerical models of varying complexity. However, the nature of the transition from a conventional state of the general circulation, with easterlies or weak westerlies in the tropics, to such a superrotating state remains unclear. Is it abrupt or continuous? This question may have far-reaching consequences, as it may provide a mechanism for abrupt climate change in a planetary atmosphere, both through the loss of stability of the conventional circulation and through potential noise-induced transitions in the bistability range. We study two feedbacks which may lead to bistability between a conventional and a superrotating state: the Hadley cell feedback and a wave-jet resonance feedback. We delineate the regime of applicability of these two mechanisms in a simple model of zonal acceleration budget at the equator. Then, we show using numerical simulations of the axisymmetric primitive equations that the wave-jet resonance feedback indeed leads to robust bistability, while the bistability governed by the Hadley cell feedback, although observed in our numerical simulations, is much more fragile in a multilevel model.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Computation of extremes values of time averaged observables in climate models with large deviation techniques

One of the goals of climate science is to characterize the statistics of extreme and potentially dangerous events in the present and future climate. Extreme events like heat waves, droughts, or floods due to persisting rains are characterized by large anomalies of the time average of an observable over a long time. The framework of Donsker-Varadhan large deviation theory could therefore be useful for their analysis. In this paper we discuss how concepts and numerical algorithms developed in relation with large deviation theory can be applied to study extreme, rare fluctuations of time averages of surface temperatures at regional scale with comprehensive numerical climate models. We study the convergence of large deviation functions for the time averaged European surface temperature obtained with direct numerical simulation of the climate model Plasim, and discuss their climate implications. We show how using a rare event algorithm can improve the efficiency of the computation of the large deviation rate functions. We discuss the relevance of the large deviation asymptotics for applications, and we show how rare event algorithms can be used also to improve the statistics of events on time scales shorter than the one needed for reaching the large deviation asymptotics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equilibrium statistical mechanics and energy partition for the shallow water model

The aim of this paper is to use large deviation theory in order to compute the entropy of macrostates for the microcanonical measure of the shallow water system. The main prediction of this full statistical mechanics computation is the energy partition between a large scale vortical flow and small scale fluctuations related to inertia-gravity waves. We introduce for that purpose a discretized model of the continuous shallow water system, and compute the corresponding statistical equilibria. We argue that microcanonical equilibrium states of the discretized model in the continuous limit are equilibrium states of the actual shallow water system. We show that the presence of small scale fluctuations selects a subclass of equilibria among the states that were previously computed by phenomenological approaches that were neglecting such fluctuations. In the limit of weak height fluctuations, the equilibrium state can be interpreted as two subsystems in thermal contact: one subsystem corresponds to the large scale vortical flow, the other subsystem corresponds to small scale height and velocity fluctuations. It is shown that either a non-zero circulation or rotation and bottom topography are required to sustain a non-zero large scale flow at equilibrium. Explicit computation of the equilibria and their energy partition is presented in the quasi-geostrophic limit for the energy-enstrophy ensemble. The possible role of small scale dissipation and shocks is discussed. A geophysical application to the Zapiola anticyclone is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Perturbative calculation of quasi-potential in non-equilibrium diffusions: a mean-field example

In stochastic systems with weak noise, the logarithm of the stationary distribution becomes proportional to a large deviation rate function called the quasi-potential. The quasi-potential, and its characterization through a variational problem, lies at the core of the Freidlin-Wentzell large deviations theory%.~\cite{freidlin1984}.In many interacting particle systems, the particle density is described by fluctuating hydrodynamics governed by Macroscopic Fluctuation Theory%, ~\cite{bertini2014},which formally fits within Freidlin-Wentzell's framework with a weak noise proportional to $1/\sqrt{N}$, where $N$ is the number of particles. The quasi-potential then appears as a natural generalization of the equilibrium free energy to non-equilibrium particle systems. A key physical and practical issue is to actually compute quasi-potentials from their variational characterization for non-equilibrium systems for which detailed balance does not hold. We discuss how to perform such a computation perturbatively in an external parameter $λ$, starting from a known quasi-potential for $λ=0$. In a general setup, explicit iterative formulae for all terms of the power-series expansion of the quasi-potential are given for the first time. The key point is a proof of solvability conditions that assure the existence of the perturbation expansion to all orders. We apply the perturbative approach to diffusive particles interacting through a mean-field potential. For such systems, the variational characterization of the quasi-potential was proven by Dawson and Gartner%. ~\cite{dawson1987,dawson1987b}. Our perturbative analysis provides new explicit results about the quasi-potential and about fluctuations of one-particle observables in a simple example of mean field diffusions: the Shinomoto-Kuramoto model of coupled rotators%. ~\cite{shinomoto1986}. This is one of few systems for which non-equilibrium free energies can be computed and analyzed in an effective way, at least perturbatively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Population dynamics method with a multi-canonical feedback control

We discuss the Giardinà-Kurchan-Peliti population dynamics method for evaluating large deviations of time averaged quantities in Markov processes [Phys. Rev. Lett. \textbf{96}, 120603 (2006)]. This method exhibits systematic errors which can be large in some circumstances, particularly for systems with weak noise, with many degrees of freedom, or close to dynamical phase transitions. We show how these errors can be mitigated by introducing control forces within the algorithm. These forces are determined by an iteration-and-feedback scheme, inspired by multicanonical methods in equilibrium sampling. We demonstrate substantially improved results in a simple model and we discuss potential applications to more complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rare event computation in deterministic chaotic systems using genealogical particle analysis

In this paper we address the use of rare event computation techniques to estimate small over-threshold probabilities of observables in determin-istic dynamical systems. We demonstrate that the genealogical particle analysis algorithms can be successfully applied to a toy model of atmospheric dynamics, the Lorenz '96 model. We furthermore use the Ornstein-Uhlenbeck system to illustrate a number of implementation issues. We also show how a time-dependent objective function based on the fluctuation path to a high threshold can greatly improve the performance of the estimator compared to a fixed-in-time objective function.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computing transition rates for the 1-D stochastic Ginzburg--Landau--Allen--Cahn equation for finite-amplitude noise with a rare event algorithm

In this paper we compute and analyse the transition rates and duration of reactive trajectories of the stochastic 1-D Allen-Cahn equations for both the Freidlin-Wentzell regime (weak noise or temperature limit) and finite-amplitude white noise, as well as for small and large domain. We demonstrate that extremely rare reactive trajectories corresponding to direct transitions between two metastable states are efficiently computed using an algorithm called adaptive multilevel splitting. This algorithm is dedicated to the computation of rare events and is able to provide ensembles of reactive trajectories in a very efficient way. In the small noise limit, our numerical results are in agreement with large-deviation predictions such as instanton-like solutions, mean first passages and escape probabilities. We show that the duration of reactive trajectories follows a Gumbel distribution like for one degree of freedom systems. Moreover, the mean duration growths logarithmically with the inverse temperature. The prefactor given by the potential curvature grows exponentially with size. The main novelty of our work is that we also perform an analysis of reactive trajectories for large noises and large domains. In this case, we show that the position of the reactive front is essentially a random walk. This time, the mean duration grows linearly with the inverse temperature and quadratically with the size. Using a phenomenological description of the system, we are able to calculate the transition rate, although the dynamics is described by neither Freidlin--Wentzell or Eyring--Kramers type of results. Numerical results confirm our analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalisation of the Eyring-Kramers transition rate formula to irreversible diffusion processes

In the small noise regime, the average transition time between metastable states of a reversible diffusion process is described at the logarithmic scale by Arrhenius' law. The Eyring-Kramers formula classically provides a subexponential prefactor to this large deviation estimate. For irreversible diffusion processes, the equivalent of Arrhenius' law is given by the Freidlin-Wentzell theory. In this paper, we compute the associated prefactor and thereby generalise the Eyring-Kramers formula to irreversible diffusion processes. In our formula, the role of the potential is played by Freidlin-Wentzell's quasipotential, and a correction depending on the non-Gibbsianness of the system along the instanton is highlighted. Our analysis relies on a WKB analysis of the quasistationary distribution of the process in metastable regions, and on a probabilistic study of the process in the neighbourhood of saddle-points of the quasipotential.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large Deviations in Fast-Slow Systems

The incidence of rare events in fast-slow systems is investigated via analysis of the large deviation principle (LDP) that characterizes the likelihood and pathway of large fluctuations of the slow variables away from their mean behavior -- such fluctuations are rare on short timescales but become ubiquitous eventually. This LDP involves an Hamilton-Jacobi equation whose Hamiltonian is related to the leading eigenvalue of the generator of the fast process, and is typically non-quadratic in the momenta -- in other words, the LDP for the slow variables in fast-slow systems is different in general from that of any stochastic differential equation (SDE) one would write for the slow variables alone. It is shown here that the eigenvalue problem for the Hamiltonian can be reduced to a simpler algebraic equation for this Hamiltonian for a specific class of systems in which the fast variables satisfy a linear equation whose coefficients depend nonlinearly on the slow variables, and the fast variables enter quadratically the equation for the slow variables. These results are illustrated via examples, inspired by kinetic theories of turbulent flows and plasma, in which the quasipotential characterizing the long time behavior of the system is calculated and shown again to be different from that of an SDE.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computation of rare transitions in the barotropic quasi-geostrophic equations

We investigate the theoretical and numerical computation of rare transitions in simple geophysical turbulent models. We consider the barotropic quasi-geostrophic and two-dimensional Navier-Stokes equations in regimes where bistability between two coexisting large-scale attractors exist. By means of large deviations and instanton theory with the use of an Onsager-Machlup path integral formalism for the transition probability, we show how one can directly compute the most probable transition path between two coexisting attractors analytically in an equilibrium (Langevin) framework and numerically otherwise. We adapt a class of numerical optimization algorithms known as minimum action methods to simple geophysical turbulent models. We show, that by numerically minimizing an appropriate action functional, in a large deviation limit, one can predict the most likely transition path for a rare transition between two states. By considering examples where theoretical predictions can be made, we show that the minimum action method successfully predicts the most likely transition path. Finally, we discuss the application and extension of such numerical optimization schemes to compute rare transitions observed in direct numerical simulations, experiments and to other, more complex, turbulent systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic theory of jet dynamics in the stochastic barotropic and 2D Navier-Stokes equations

We discuss the dynamics of zonal (or unidirectional) jets for barotropic flows forced by Gaussian stochastic fields with white in time correlation functions. This problem contains the stochastic dynamics of 2D Navier-Stokes equation as a special case. We consider the limit of weak forces and dissipation, when there is a time scale separation between the inertial time scale (fast) and the spin-up or spin-down time (large) needed to reach an average energy balance. In this limit, we show that an adiabatic reduction (or stochastic averaging) of the dynamics can be performed. We then obtain a kinetic equation that describes the slow evolution of zonal jets over a very long time scale, where the effect of non-zonal turbulence has been integrated out. The main theoretical difficulty, achieved in this work, is to analyze the stationary distribution of a Lyapunov equation that describes quasi-Gaussian fluctuations around each zonal jet, in the inertial limit. This is necessary to prove that there is no ultraviolet divergence at leading order in such a way that the asymptotic expansion is self-consistent. We obtain at leading order a Fokker--Planck equation, associated to a stochastic kinetic equation, that describes the slow jet dynamics. Its deterministic part is related to well known phenomenological theories (for instance Stochastic Structural Stability Theory) and to quasi-linear approximations, whereas the stochastic part allows to go beyond the computation of the most probable zonal jet. We argue that the effect of the stochastic part may be of huge importance when, as for instance in the proximity of phase transitions, more than one attractor of the dynamics is present.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Statistical mechanics of the 3D axi-symmetric Euler equations in a Taylor-Couette geometry

In the present paper, microcanonical measures for the dynamics of three dimensional (3D) axially symmetric turbulent flows with swirl in a Taylor-Couette geometry are defined, using an analogy with a long-range lattice model. We compute the relevant physical quantities and argue that two kinds of equilibrium regimes exist, depending on the value of the total kinetic energy. For low energies, the equilibrium flow consists of a purely swirling flow whose toroidal profile depends on the radial coordinate only. For high energies, the typical toroidal field is uniform, while the typical poloidal field is organized into either a single vertical jet or a large scale dipole, and exhibits infinite fluctuations. This unusual phase diagram comes from the poloidal fluctuations not being bounded for the axi-symmetric Euler dynamics, even though the latter conserve infinitely many ''Casimir invariants''. This shows that 3D axially symmetric flows can be considered as intermediate between 2D and 3D flows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A complete theory of low-energy phase diagrams for two-dimensional turbulence steady states and equilibria

For the 2D Euler equations and related models of geophysical flows, minima of energy--Casimir variational problems are stable steady states of the equations (Arnol'd theorems). The same variational problems also describe sets of statistical equilibria of the equations. In this paper, we make use of Lyapunov--Schmidt reduction in order to study the bifurcation diagrams for these variational problems, in the limit of small energy or, equivalently, of small departure from quadratic Casimir functionals. We show a generic occurrence of phase transitions, either continuous or discontinuous. We derive the type of phase transitions for any domain geometry and any model analogous to the 2D Euler equations. The bifurcations depend crucially on a_4, the quartic coefficient in the Taylor expansion of the Casimir functional around its minima. Note that a_4 can be related to the fourth moment of the vorticity in the statistical mechanics framework. A tricritical point (bifurcation from a continuous to a discontinuous phase transition) often occurs when a_4 changes sign. The bifurcations depend also on possible constraints on the variational problems (circulation, energy). These results show that the analytical results obtained with quadratic Casimir functionals by several authors are non-generic (not robust to a small change in the parameters).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinetic theory for non-equilibrium stationary states in long-range interacting systems

We study long-range interacting systems perturbed by external stochastic forces. Unlike the case of short-range systems, where stochastic forces usually act locally on each particle, here we consider perturbations by external stochastic fields. The system reaches stationary states where external forces balance dissipation on average. These states do not respect detailed balance and support non-vanishing fluxes of conserved quantities. We generalize the kinetic theory of isolated long-range systems to describe the dynamics of this non-equilibrium problem. The kinetic equation that we obtain applies to plasmas, self-gravitating systems, and to a broad class of other systems. Our theoretical results hold for homogeneous states, but may also be generalized to apply to inhomogeneous states. We obtain an excellent agreement between our theoretical predictions and numerical simulations. We discuss possible applications to describe non-equilibrium phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinetic theory of nonequilibrium stochastic long-range systems: Phase transition and bistability

We study long-range interacting systems driven by external stochastic forces that act collectively on all the particles constituting the system. Such a scenario is frequently encountered in the context of plasmas, self-gravitating systems, two-dimensional turbulence, and also in a broad class of other systems. Under the effect of stochastic driving, the system reaches a stationary state where external forces balance dissipation on average. These states have the invariant probability that does not respect detailed balance, and are characterized by non-vanishing currents of conserved quantities. In order to analyze spatially homogeneous stationary states, we develop a kinetic approach that generalizes the one known for deterministic long-range systems; we obtain a very good agreement between predictions from kinetic theory and extensive numerical simulations. Our approach may also be generalized to describe spatially inhomogeneous stationary states. We also report on numerical simulations exhibiting a first-order nonequilibrium phase transition from homogeneous to inhomogeneous states. Close to the phase transition, the system shows bistable behavior between the two states, with a mean residence time that diverges as an exponential in the inverse of the strength of the external stochastic forces, in the limit of low values of such forces.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-classical large deviations for a noisy system with non-isolated attractors

We study the large deviations of a simple noise-perturbed dynamical system having continuous sets of steady states, which mimick those found in some partial differential equations related, for example, to turbulence problems. The system is a two-dimensional nonlinear Langevin equation involving a dissipative, non-potential force, which has the essential effect of creating a line of stable fixed points (attracting line) touching a line of unstable fixed points (repelling line). Using different analytical and numerical techniques, we show that the stationary distribution of this system satisfies in the low-noise limit a large deviation principle containing two competing terms: i) a classical but sub-dominant large deviation term, which can be derived from the Freidlin-Wentzell theory of large deviations by studying the fluctuation paths or instantons of the system near the attracting line, and ii) a dominant large deviation term, which does not follow from the Freidlin-Wentzell theory, as it is related to fluctuation paths of zero action, referred to as sub-instantons, emanating from the repelling line. We discuss the nature of these sub-instantons, and show how they arise from the connection between the attracting and repelling lines. We also discuss in a more general way how we expect these to arise in more general stochastic systems having connected sets of stable and unstable fixed points, and how they should determine the large deviation properties of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sampling microcanonical measures of the 2D Euler equations through Creutz's algorithm: a phase transition from disorder to order when energy is increased

The 2D Euler equations is the basic example of fluid models for which a microcanical measure can be constructed from first principles. This measure is defined through finite-dimensional approximations and a limiting procedure. Creutz's algorithm is a microcanonical generalization of the Metropolis-Hasting algorithm (to sample Gibbs measures, in the canonical ensemble). We prove that Creutz's algorithm can sample finite-dimensional approximations of the 2D Euler microcanonical measures (incorporating fixed energy and other invariants). This is essential as microcanonical and canonical measures are known to be inequivalent at some values of energy and vorticity distribution. Creutz's algorithm is used to check predictions from the mean-field statistical mechanics theory of the 2D Euler equations (the Robert-Sommeria-Miller theory). We found full agreement with theory. Three different ways to compute the temperature give consistent results. Using Creutz's algorithm, a first-order phase transition never observed previously, and a situation of statistical ensemble inequivalence are found and studied. Strikingly, and contrasting usual statistical mechanics interpretations, this phase transition appears from a disordered phase to an ordered phase (with less symmetries) when energy is increased. We explain this paradox.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Equivalence of the Lagrangian-Averaged Navier-Stokes-α Model and the Rational LES model in Two Dimensions

In the Large Eddy Simulation (LES) framework for modeling a turbulent flow, when the large scale velocity field is defined by low-pass filtering the full velocity field, a Taylor series expansion of the full velocity field in terms of the large scale velocity field leads (at the leading order) to the nonlinear gradient model for the subfilter stresses. Motivated by the fact that while the nonlinear gradient model shows excellent a priori agreement in resolved simulations, the use of this model by itself is problematic, we consider two models that are related, but better behaved: The Rational LES model that uses a sub-diagonal Pade approximation instead of a Taylor series expansion and the Lagrangian Averaged Navier-Stokes-α model that uses a regulariza- tion approach to modeling turbulence. In this article, we show that these two latter models are identical in two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Oceanic rings and jets as statistical equilibrium states

Equilibrium statistical mechanics of two-dimensional flows provides an explanation and a prediction for the self-organization of large scale coherent structures. This theory is applied in this paper to the description of oceanic rings and jets, in the framework of a 1.5 layer quasi-geostrophic model. The theory predicts the spontaneous formation of regions where the potential vorticity is homogenized, with strong and localized jets at their interface. Mesoscale rings are shown to be close to a statistical equilibrium: the theory accounts for their shape, their drift, and their ubiquity in the ocean, independently of the underlying generation mechanism. At basin scale, inertial states presenting mid basin eastward jets (and then different from the classical Fofonoff solution) are described as marginally unstable states. These states are shown to be marginally unstable for the equilibrium statistical theory. In that case, considering a purely inertial limit is a first step toward more comprehensive out of equilibrium studies that would take into account other essential aspects, such as wind forcing.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Solvable phase diagrams and ensemble inequivalence for two-dimensional and geophysical turbulent flows

Using explicit analytical computations, generic occurrence of inequivalence between two or more statistical ensembles is obtained for a large class of equilibrium states of two-dimensional and geophysical turbulent flows. The occurrence of statistical ensemble inequivalence is shown to be related to previously observed phase transitions in the equilibrium flow topology. We find in these turbulent flow equilibria, two mechanisms for the appearance of ensemble equivalences, that were not observed in any physical systems before. These mechanisms are associated respectively with second-order azeotropy (simultaneous appearance of two second-order phase transitions), and with bicritical points (bifurcation from a first-order to two second-order phase transition lines). The important roles of domain geometry, of topography, and of a screening length scale (the Rossby radius of deformation) are discussed. It is found that decreasing the screening length scale (making interactions more local) surprisingly widens the range of parameters associated with ensemble inequivalence. These results are then generalized to a larger class of models, and applied to a complete description of an academic model for inertial oceanic circulation, the Fofonoff flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistical mechanics of two-dimensional and geophysical flows

The theoretical study of the self-organization of two-dimensional and geophysical turbulent flows is addressed based on statistical mechanics methods. This review is a self-contained presentation of classical and recent works on this subject; from the statistical mechanics basis of the theory up to applications to Jupiter's troposphere and ocean vortices and jets. Emphasize has been placed on examples with available analytical treatment in order to favor better understanding of the physics and dynamics. The equilibrium microcanonical measure is built from the Liouville theorem. On this theoretical basis, we predict the output of the long time evolution of complex turbulent flows as statistical equilibria. This is applied to make quantitative models of two-dimensional turbulence, the Great Red Spot and other Jovian vortices, ocean jets like the Gulf-Stream, and ocean vortices. We also present recent results for non-equilibrium situations, for the studies of either the relaxation towards equilibrium or non-equilibrium steady states.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Invariant measures of the 2D Euler and Vlasov equations

We discuss invariant measures of partial differential equations such as the 2D Euler or Vlasov equations. For the 2D Euler equations, starting from the Liouville theorem, valid for N-dimensional approximations of the dynamics, we define the microcanonical measure as a limit measure where N goes to infinity. When only the energy and enstrophy invariants are taken into account, we give an explicit computation to prove the following result: the microcanonical measure is actually a Young measure corresponding to the maximization of a mean-field entropy. We explain why this result remains true for more general microcanonical measures, when all the dynamical invariants are taken into account. We give an explicit proof that these microcanonical measures are invariant measures for the dynamics of the 2D Euler equations. We describe a more general set of invariant measures, and discuss briefly their stability and their consequence for the ergodicity of the 2D Euler equations. The extension of these results to the Vlasov equations is also discussed, together with a proof of the uniqueness of statistical equilibria, for Vlasov equations with repulsive convex potentials. Even if we consider, in this paper, invariant measures only for Hamiltonian equations, with no fluxes of conserved quantities, we think this work is an important step towards the description of non-equilibrium invariant measures with fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large time behavior and asymptotic stability of the two-dimensional Euler and linearized Euler equations

We study the asymptotic behavior and the asymptotic stability of the two-dimensional Euler equations and of the two-dimensional linearized Euler equations close to parallel flows. We focus on spectrally stable jet profiles $U(y)$ with stationary streamlines $y_{0}$ such that $U'(y_{0})=0$, a case that has not been studied previously. We describe a new dynamical phenomenon: the depletion of the vorticity at the stationary streamlines. An unexpected consequence, is that the velocity decays for large times with power laws, similarly to what happens in the case of the Orr mechanism for base flows without stationary streamlines. The asymptotic behaviors of velocity and the asymptotic profiles of vorticity are theoretically predicted and compared with direct numerical simulations. We argue on the asymptotic stability of these flow velocities even in the absence of any dissipative mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermodynamics and dynamics of systems with long-range interactions

Thermodynamic and dynamical properties of systems with long-range pairwise interactions (LRI), which decay as $1/r^{d+σ}$ at large distances $r$ in $d$ dimensions, are reviewed. Two broad classes of such systems are discussed. (i) Systems with a slow decay of the interactions, termed "strong" LRI, where the energy is super-extensive. These systems are characterized by unusual properties such as inequivalence of ensembles, negative specific heat, slow decay of correlations, anomalous diffusion and ergodicity breaking. (ii) Systems with faster decay of the interaction potential, where the energy is additive, thus resulting in less dramatic effects. These interactions affect the thermodynamic behavior of systems near phase transitions, where long-range correlations are naturally present. Long-range correlations are often present in systems driven out of equilibrium when the dynamics involves conserved quantities. Steady state properties of driven systems with local dynamics are considered within the framework outlined above.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Ensemble inequivalence, bicritical points and azeotropy for generalized Fofonoff flows

We present a theoretical description for the equilibrium states of a large class of models of two-dimensional and geophysical flows, in arbitrary domains. We account for the existence of ensemble inequivalence and negative specific heat in those models, for the first time using explicit computations. We give exact theoretical computation of a criteria to determine phase transition location and type. Strikingly, this criteria does not depend on the model, but only on the domain geometry. We report the first example of bicritical points and second order azeotropy in the context of systems with long range interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Phase space gaps and ergodicity breaking in systems with long range interactions

We study a generalized isotropic XY-model which includes both two-spin and four-spin mean-field interactions. This model can be solved in the microcanonical ensemble. It is shown that in certain parameter regions the model exhibits gaps in the magnetization at fixed energy, resulting in ergodicity breaking. This phenomenon has previously been reported in anisotropic and discrete spin models. The entropy of the model is calculated and the microcanonical phase diagram is derived, showing the existence of first order phase transitions from the ferromagnetic to a paramagnetic disordered phase. It is found that ergodicity breaking takes place both in the ferromagnetic and the paramagnetic phases. As a consequence, the system can exhibit a stable ferromagnetic phase within the paramagnetic region, and conversely a disordered phase within the magnetically ordered region.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Classification of phase transitions and ensemble inequivalence, in systems with long range interactions

Systems with long range interactions in general are not additive, which can lead to an inequivalence of the microcanonical and canonical ensembles. The microcanonical ensemble may show richer behavior than the canonical one, including negative specific heats and other non-common behaviors. We propose a classification of microcanonical phase transitions, of their link to canonical ones, and of the possible situations of ensemble inequivalence. We discuss previously observed phase transitions and inequivalence in self-gravitating, two-dimensional fluid dynamics and non-neutral plasmas. We note a number of generic situations that have not yet been observed in such systems.

0 citations0 reviewsNo code linkedOpen access