Search

The dynamics of mechanical systems such as turbomachinery with multiple blades are often modeled by arrays of periodically driven coupled nonlinear oscillators. It is known that such systems may have multiple stable vibrational modes, and transitions between them may occur under the influence of random factors. A methodology for finding most probable escape paths and estimating the transition rates in the small noise limit is developed and applied to a collection of arrays of coupled monostable oscillators with cubic nonlinearity, small damping, and harmonic external forcing. The methodology is built upon the action plot method (Beri et al. 2005) and relies on the large deviation theory, optimal control theory, and the Floquet theory. The action plot method is promoted to non-autonomous high-dimensional systems, and a method for solving the arising optimization problem with discontinuous objective function restricted to a certain manifold is proposed. The most probable escape paths between stable vibrational modes in arrays of up to five oscillators and the corresponding quasipotential barriers are computed and visualized. The dependence of the quasipotential barrier on the parameters of the system is discussed., Comment: 44 pages, 30 figures, 3 tables

The study of noise-driven transitions occurring rarely on the time-scale of systems modeled by SDEs is of crucial importance for understanding such phenomena as genetic switches in living organisms and magnetization switches of the Earth. For a gradient SDE, the predictions for transition times and paths between its metastable states are done using the potential function. For a nongradient SDE, one needs to decompose its forcing into a gradient of the so-called quasipotential and a rotational component, which cannot be done analytically in general. We propose a methodology for computing the quasipotential for highly dissipative and chaotic systems built on the example of Lorenz'63 with an added stochastic term. It is based on the ordered line integral method, a Dijkstra-like quasipotential solver, and combines 3D computations in whole regions, a dimensional reduction technique, and 2D computations on radial meshes on manifolds or their unions. Our collection of source codes is available on M. Cameron's web page and on GitHub., Comment: 38 pages, 17 figures, 5 3D figures are linked to movies posted on youtube, software is available on GitHub

QPot is an R package for analyzing two-dimensional systems of stochastic differential equations. It provides users with a wide range of tools to simulate, analyze, and visualize the dynamics of these systems. One of QPot's key features is the computation of the quasi-potential, an important tool for studying stochastic systems. Quasi-potentials are particularly useful for comparing the relative stabilities of equilibria in systems with alternative stable states. This paper describes QPot's primary functions, and explains how quasi-potentials can yield insights about the dynamics of stochastic systems. Three worked examples guide users through the application of QPot's functions.