28 research outputs found
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 deterministic dynamical systems. We demonstrate that 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
Relevance of instantons in Burgers turbulence
Instanton calculations are performed in the context of stationary Burgers turbulence to estimate the tails of the probability density function (PDF) of velocity gradients. These results are then compared to those obtained from massive direct numerical simulations (DNS) of the randomly forced Burgers equation. The instanton predictions are shown to agree with the DNS in a wide range of regimes, including those that are far from the limiting cases previously considered in the literature. These results settle the controversy of the relevance of the instanton approach for the prediction of the velocity gradient PDF tail exponents. They also demonstrate the usefulness of the instanton formalism in Burgers turbulence, and suggest that this approach may be applicable in other contexts, such as 2D and 3D turbulence in compressible and incompressible flows
Arclength parametrized Hamilton's equations for the calculation of instantons
A method is presented to compute minimizers (instantons) of action functionals using arclength parametrization of Hamilton's equations. This method can be interpreted as a local variant of the geometric minimum action method introduced to compute minimizers of the Freidlin--Wentzell action functional that arises in the context of large deviation theory for stochastic differential equations. The method is particularly well suited to calculate expectations dominated by noise-induced excursions from deterministically stable fixpoints. Its simplicity and computational efficiency are illustrated here using several examples: a finite-dimensional stochastic dynamical system (an Ornstein--Uhlenbeck model) and two models based on stochastic partial differential equations: the -model and the stochastically driven Burgers equation