Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lyapunov exponent”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Stochastic stability properties of jump linear systems

Jump linear systems are defined as a family of linear systems with randomly jumping parameters (usually governed by a Markov jump process) and are used to model systems subject to failures or changes in structure. The authors study stochastic stability properties in jump linear systems and the relationship among various moment and sample path stability properties. It is shown that all second moment stability properties are equivalent and are sufficient for almost sure sample path stability, and a testable necessary and sufficient condition for second moment stability is derived. The Lyapunov exponent method for the study of almost sure sample stability is discussed, and a theorem which characterizes the Lyapunov exponents of jump linear systems is presented.

Feng, Xiangbo↗
NASA NTRS1989

One-degree-of-freedom motion induced by modeled vortex shedding

The motion of an elastically supported cylinder forced by a nonlinear, quasi-static, aerodynamic model with the unusual feature of a motion-dependent forcing frequency was studied. Numerical solutions for the motion and the Lyapunov exponents are presented for three forcing amplitudes and two frequencies (1.0 and 1.1 times the Strouhal frequency). Initially, positive Lyapunov exponents occur and the motion can appear chaotic. After thousands of characteristic times, the motion changes to a motion (verified analytically) that is periodic and damped. This periodic, damped motion was not observed experimentally, thus raising questions concerning the modeling.

Yates, L. A.↗
NASA NTRS1993

Large-scale isentropic mixing properties of the Antarctic polar vortex from analyzed winds

Winds derived from analyzed geopotential height fields are used to study quasi-horizontal mixing by the large-scale flow in the lower stratosphere during austral spring. This is the period when the Antarctic ozone hole appears and disappears. Trajectories are computed for large ensembles of particles initially inside and outside the main polar vortex. Mixing and transport are diagnosed through estimates of finite time Lyapunov exponents and Lagrangian dispersion statistics of the tracer trajectories. At 450 K and above prior to the vortex breakdown: Lyapunov exponents are a factor of 2 smaller inside the vortex than outside; diffusion coefficients are an order of magnitude smaller inside than outside the vortex; and the trajectories reveal little exchange of air across the vortex boundary. At lower levels (425 and 400 K) mixing is greater, and there is substantial exchange of air across the vortex boundary. In some years there are large wave events that expel small amounts of vortex air into the mid-latitudes. At the end of the spring season during the vortex breakdown there is rapid mixing of air across the vortex boundary, which is evident in the mixing diagnostics and the tracer trajectories.

Bowman, Kenneth P.↗
DOE OSTI2022

Holography of the photon ring

Space-based next-generation interferometers propose to measure the Lyapunov exponents of the nearly bound geodesics that comprise the photon ring surrounding the black hole M87*. We argue that these classical Lyapunov exponents equal the quantum Ruelle resonances describing the late-time approach to thermal equilibrium of the quantum microstate holographically dual to any Kerr black hole such as M87*. Moreover, we identify 'near-ring regions' in the phase space of fields propagating on Kerr that exhibit critical behavior, including emergent conformal symmetries. These are analogues for sub-extremal Kerr of the much-studied 'near-horizon regions' of (near-)extremal black holes. Here, the emergent conformal symmetries greatly constrain the observational predictions for the fine photon ring substructure around M87* and for quasinormal gravitational-wave ringdowns, as well as any proposal for a quantum holographic dual to the Kerr black hole. More generally, we hope that our identification of several universal features of Kerr spectroscopy provides a useful starting point for a bottom-up approach to holography for astrophysical black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Extension of Liouville Formalism to Postinstability Dynamics

A mathematical formalism has been developed for predicting the postinstability motions of a dynamic system governed by a system of nonlinear equations and subject to initial conditions. Previously, there was no general method for prediction and mathematical modeling of postinstability behaviors (e.g., chaos and turbulence) in such a system. The formalism of nonlinear dynamics does not afford means to discriminate between stable and unstable motions: an additional stability analysis is necessary for such discrimination. However, an additional stability analysis does not suggest any modifications of a mathematical model that would enable the model to describe postinstability motions efficiently. The most important type of instability that necessitates a postinstability description is associated with positive Lyapunov exponents. Such an instability leads to exponential growth of small errors in initial conditions or, equivalently, exponential divergence of neighboring trajectories. The development of the present formalism was undertaken in an effort to remove positive Lyapunov exponents. The means chosen to accomplish this is coupling of the governing dynamical equations with the corresponding Liouville equation that describes the evolution of the flow of error probability. The underlying idea is to suppress the divergences of different trajectories that correspond to different initial conditions, without affecting a target trajectory, which is one that starts with prescribed initial conditions.

Zak, Michail↗
DOE OSTI2025

Fuzzy spheres in stringy matrix models: quantifying chaos in a mixed phase space

We consider a truncation of the BMN matrix model to a configuration of two fuzzy spheres, described by two coupled non-linear oscillators dependent on the mass parameter μ. The classical phase diagram of the system generically (μ ≠ 0) contains three equilibrium points: two centers and a center-saddle; as μ → 0 the system exhibits a pitchfork bifurcation. We demonstrate that the system is exactly integrable in quadratures for μ = 0, while for very large values of μ, it approaches another integrable point characterized by two harmonic oscillators. The classical phase space is mixed, containing both integrable islands and chaotic regions, as evidenced by the classical Lyapunov spectrum. At the quantum level, we explore indicators of early and late time chaos. The eigenvalue spacing is best described by a Brody distribution, which interpolates between Poisson and Wigner distributions; it dovetails, at the quantum level, the classical results and reemphasizes the notion that the quantum system is mixed. We also study the spectral form factor and the quantum Lyapunov exponent, as defined by out-of-time-ordered correlators. These two indicators of quantum chaos exhibit weak correlations with the Brody distribution. We speculate that the behavior of the system as μ → 0 dominates the spectral form factor and the quantum Lyapunov exponent, making these indicators of quantum chaos less effective in the context of a mixed phase space.

AdS-CFT correspondence↗
DOE OSTI2022

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Chaos in a spatially-developing plane mixing layer

A spatially-developing plane mixing layer was analyzed for chaotic behavior. A direct numerical simulation of the Navier-Stokes equations in a 2-D domain infinite in y and having inflow-outflow boundary conditions in x was used for data. Spectra, correlation dimension and the largest Lyapunov exponent were computed as functions of downstream distance x. When forced at a single (fundamental) frequency with maximum amplitude, the flow is periodic at the inflow but becomes aperiodic with increasing x. The aperiodic behavior is caused by the presence of a noisy subharmonic caused by the feedback between the necessarily nonphysical inflow and outflow boundary conditions. In order to overshadow this noise the flow was also studied with the same fundamental forcing and added random forcing of amplitude upsilon prime sub R/delta U = 0.01 at the inlet. Results were qualitatively the same in both cases: for small x, spectral peaks were sharp and dimension was nearly 1, but as x increased a narrowband spectral peak grew, spectra decayed exponentially at high frequencies and dimension increased to greater than 3. Based on these results, the flow appears to exhibit deterministic chaos. However, at no location was the largest Lyapunov exponent found to be significantly greater than zero.

Broze, J. G.↗
DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗
DOE OSTI2022

Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters

Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraining chaos: Enforcing dynamical invariants in the training of reservoir computers

Drawing on ergodic theory, we introduce a novel training method for machine learning based forecasting methods for chaotic dynamical systems. The training enforces dynamical invariants—such as the Lyapunov exponent spectrum and the fractal dimension—in the systems of interest, enabling longer and more stable forecasts when operating with limited data. The technique is demonstrated in detail using reservoir computing, a specific kind of recurrent neural network. Finally, results are given for the Lorenz 1996 chaotic dynamical system and a spectral quasi-geostrophic model of the atmosphere, both typical test cases for numerical weather prediction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Towards improved neutral exhaust in the HSX stellarator

To improve our knowledge of neutral exhaust in non-resonant divertors, a particle exhaust study is performed for the Helically Symmetric eXperiment (HSX). The chaotic magnetic field topology in the plasma edge is examined with connection lengths and Lyapunov exponents. Plasma flux tubes are identified and these flux tubes deposit particles onto plasma facing components. This helps determine the areas where neutrals are generated and provides guidance on where to put baffles. As shown using EMC3-EIRENE simulations, the introduction of baffles can help increase the neutral pressure at and away from the strike line locations. Furthermore, it is demonstrated that baffle design cannot only be based on placing structures at locations with small connection lengths.

Baffles↗
DOE OSTI2021

Green’s functions on a renormalized lattice: An improved method for the integer quantum Hall transition

Highlights: • Speed up of numerical calculations by exact renormalization group steps. • Practical approach to handle delicate scaling analysis. • Absence of evidence to prefer marginal corrections over power-law corrections. We introduce a performance-optimized method to simulate localization problems on bipartite tight-binding lattices. It combines an exact renormalization group step to reduce the sparseness of the original problem with the recursive Green’s function method. We apply this framework to investigate the critical behavior of the integer quantum Hall transition of a tight-binding Hamiltonian defined on a simple square lattice. In addition, we employ an improved scaling analysis that includes two irrelevant exponents to characterize the shift of the critical energy as well as the corrections to the dimensionless Lyapunov exponent. We compare our findings with the results of a conventional implementation of the recursive Green’s function method, and we put them into broader perspective in view of recent development in this field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Metal-insulator transition in a random Hubbard model

In this work, we examine the metal-insulator transition in a half-filled Hubbard model of electrons with random and all-to-all hopping and exchange, and an on-site nonrandom repulsion, the Hubbard $\textit{U}$. We argue that recent numerical results of Cha et al. can be understood in terms of a deconfined critical point between a disordered Fermi liquid and an insulating spin glass. We find a deconfined critical point in a previously proposed large $\textit{M}$ theory, which generalizes the SU(2) spin symmetry to SU($\textit{M}$) and obtains exponents for the electron and spin correlators which agree with those of Cha et al. This critical point is expected to have the maximal quantum Lyapunov exponent. We also present a renormalization group analysis, and argue for the presence of an additional metallic spin glass phase at half-filling and small $\textit{U}$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Early-Time Exponential Instabilities in Nonchaotic Quantum Systems

The majority of classical dynamical systems are chaotic and exhibit the butterfly effect: a minute change in initial conditions has exponentially large effects later on. But this phenomenon is difficult to reconcile with quantum mechanics. One of the main goals in the field of quantum chaos is to establish a correspondence between the dynamics of classical chaotic systems and their quantum counterparts. In isolated systems in the absence of decoherence, there is such a correspondence in dynamics, but it usually persists only over a short time window, after which quantum interference washes out classical chaos. In this work, we demonstrate that quantum mechanics can also play the opposite role and generate exponential instabilities in classically nonchaotic systems within this early-time window. Our calculations employ the out-of-time-ordered correlator (OTOC)—a diagnostic that reduces to the Lyapunov exponent in the classical limit but is well defined for general quantum systems. We show that certain classically nonchaotic models, such as polygonal billiards, demonstrate a Lyapunov-like exponential growth of the OTOC at early times with Planck’s-constant-dependent rates. This behavior is sharply contrasted with the slow early-time growth of the analog of the OTOC in the systems’ classical counterparts. These results suggest that classical-to-quantum correspondence in dynamics is violated in the OTOC even before quantum interference develops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Chaos as the cause of randomness in tearing mode onset times in DIII-D ITER baseline scenario plasmas

The first evidence of widespread chaotic plasma dynamics relevant to the onset of disruptive tearing modes is found in DIII-D tokamak discharges. This evidence is obtained by calculating positive maximal Lyapunov exponents from Mirnov signals using the Rosenstein algorithm, implying the chaotic growth and/or rotation of magnetohydrodynamic modes in DIII-D [Luxon, Nucl. Fusion 42, 614 (2002)] plasmas. Such chaotic dynamics offer an explanation for the random tearing mode onset times previously observed in ITER baseline scenario discharges. Further, the Lyapunov timescales associated with this chaos could inhibit reliable prediction of tearing mode onset beyond the angular momentum confinement timescale; this constraint should inform the design of plasma control systems on ITER.

Chaotic dynamics↗
DOE OSTI2022

Observability of fidelity decay at the Lyapunov rate in few-qubit quantum simulations

In certain regimes, the fidelity of quantum states will decay at a rate set by the classical Lyapunov exponent. This serves both as one of the most important examples of the quantum-classical correspondence principle and as an accurate test for the presence of chaos. While detecting this phenomenon is one of the first useful calculations that noisy quantum computers without error correction can perform, a thorough study of the quantum sawtooth map reveals that observing the Lyapunov regime is just beyond the reach of present-day devices. We prove that there are three bounds on the ability of any device to observe the Lyapunov regime and give the first quantitatively accurate description of these bounds: (1) the Fermi golden rule decay rate must be larger than the Lyapunov rate, (2) the quantum dynamics must be diffusive rather than localized, and (3) the initial decay rate must be slow enough for Lyapunov decay to be observable. This last bound, which has not been recognized previously, places a limit on the maximum amount of noise that can be tolerated. The theory implies that an absolute minimum of 6 qubits is required. Recent experiments on IBM-Q and IonQ imply that some combination of a noise reduction by up to 100x per gate and large increases in connectivity and gate parallelization are also necessary. Finally, scaling arguments are given that quantify the ability of future devices to observe the Lyapunov regime based on trade-offs between hardware architecture and performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Local Poincaré algebra from quantum chaos

The local two-dimensional Poincaré algebra near the horizon of an eternal AdS black hole, or in proximity to any bifurcate Killing horizon, is generated by the Killing flow and outward null translations on the horizon. In holography, this local Poincaré algebra is reflected as a pair of unitary flows in the boundary Hilbert space whose generators under modular flow grow and decay exponentially with a maximal Lyapunov exponent. This is a universal feature of many geometric vacua of quantum gravity. To explain this universality, we show that a two-dimensional Poincaré algebra emerges in any quantum system that has von Neumann subalgebras associated with half-infinite modular time intervals (modular future and past subalgebras) in a limit analogous to the near-horizon limit. In ergodic theory, quantum dynamical systems with future or past algebras are called quantum K-systems. The surprising statement is that modular K-systems are always maximally chaotic. Interacting quantum systems in the thermodynamic limit and large N theories above the Hawking-Page phase transition are examples of physical theories with future/past subalgebras. We prove that the existence of (modular) future/past von Neumann subalgebras also implies a second law of (modular) thermodynamics and the exponential decay of (modular) correlators. We generalize our results from the modular flow to any dynamical flow with a positive generator and interpret the positivity condition as quantum detailed balance.

79 ASTRONOMY AND ASTROPHYSICS↗