BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Freddy Bouchet

Freddy Bouchet contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechphysics.flu-dynphysics.ao-phnlin.CDphysics.comp-phphysics.data-anphysics.plasm-phastro-ph.HEMachine Learningmath.DS

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 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
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