Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chaotic dynamical systems”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

McMillan map and nonlinear Twiss parameters

In this article we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings originally introduced by Edwin McMillan; the second one is also known as canonical McMillan map. Both of them are simplest symmetric McMillan maps with only one intrinsic parameter, the trace of the Jacobian at the fixed point. While these dynamical systems have numerous of applications and are used in many areas of math and physics, some of their dynamical properties have not been described yet. We fulfill the gap and provide complete description of all stable trajectories including parametrization of invariant curves, Pioncaré rotation numbers and canonical action-angle variables. In the second part we relate these maps with general chaotic map in McMillan-Turaev form. We show that McMillan sextupole and octupole mappings are first order approximations of dynamics around the fixed point, in a similar way as linear map and quadratic invariant (Courant-Snyder invariant in accelerator physics) is the zeroth order approximation (known as linearization). Finally we suggest the new formalism of nonlinear Twiss parameters which incorporate dependence of rotation number as a function of amplitude, in contrast to e.g. betatron phase advance used in accelerator physics which is independent of amplitude. Specifically in application to accelerator physics this new formalism is capable of predicting dynamical aperture around 1-st, 2-nd, 3-rd and 4-th order resonances for flat beams, which is critical for beam injection/extraction.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. Here, in this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region—a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

active nematics↗
DOE OSTI2022

Resilient Hierarchical Networked Control Systems: Secure Controls for Critical Locations and at Edge

Integration of information and communication technology (ICT) offers new opportunities in improving the management and operation of critical infrastructures such as power systems as it allows connection of different sensors and control components via a communication network, leading to the so-called networked control systems (NCS). However, the use of open and pervasive ICT such as the Internet or wireless communication technologies comes at a price of making NCS vulnerable to cyber intrusions/attacks which may cause physical damage. Here, this chapter presents control algorithms to ensure resilient and safe operation of NCS under unknown cyberattacks. Specifically, a variant of dynamic watermarking strategies is presented by embedding encoding/decoding components of chaotic signals into the NCS for secure control for critical locations where the measurement/control signals are transmitted to/from the control center via a communication network. In addition, resilient cooperative control algorithms are discussed to ensure safe operation at edge of the NCS which consists of a large number of distributed controllable devices. Several numerical examples are provided to illustrate the proposed control strategies.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
NASA NTRS2005

Stochastic Representation of Chaos using Terminal Attractors

A nonlinear version of the Liouville equation based upon terminal attractors is proposed for describing post-instability motions of dynamical systems with exponential divergence of trajectories such as those leading to chaos and turbulence. As a result, the post-instability motions are represented by expectations, variances, and higher moments of the state variables as functions of time. The proposed approach can be applied to conservative chaos, and in particular, to n-bodies problem, as well as to dissipative systems, and in particular, to chaotic attractors and turbulence.

stochastic processes↗
NASA NTRS2009

Unsteady Full Annulus Simulations of a Transonic Axial Compressor Stage

Two recent research endeavors in turbomachinery at NASA Glenn Research Center have focused on compression system stall inception and compression system aerothermodynamic performance. Physical experiment and computational research are ongoing in support of these research objectives. TURBO, an unsteady, three-dimensional, Navier-Stokes computational fluid dynamics code commissioned and developed by NASA, has been utilized, enhanced, and validated in support of these endeavors. In the research which follows, TURBO is shown to accurately capture compression system flow range-from choke to stall inception-and also to accurately calculate fundamental aerothermodynamic performance parameters. Rigorous full-annulus calculations are performed to validate TURBO s ability to simulate the unstable, unsteady, chaotic stall inception process; as part of these efforts, full-annulus calculations are also performed at a condition approaching choke to further document TURBO s capabilities to compute aerothermodynamic performance data and support a NASA code assessment effort.

Herrick, Gregory P.↗
NASA NTRS1988

Dynamics of coorbital satellite rings

The dynamical behavior of a coorbital satellite ring is studied for N = 2-9 satellites, in terms of a simplified dynamical description, in which the motion is reduced to the separation angles between satellites. The number and stability of different kinds of stationary configurations is explored, revealing that equally spaced rings are not stable against small perturbations for N not greater than 6, while for N = 2-8 there exists another, stable compact solution. Integrations of exact equations confirm these results. Moreover, the systems are found to display chaotic characteristics for a certain range of energy. The behavior can be interpreted in terms of maximum velocity curves, defining the allowed region of the motion in the phase space.

Salo, H.↗
NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗
NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗
DOE OSTI2021

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗
NASA NTRS2007

Fuel-Optimal Trajectories in a Planet-Moon Environment Using Multiple Gravity Assists

For low energy spacecraft trajectories such as multi-moon orbiters for the Jupiter system, multiple gravity assists by moons could be used in conjunction with ballistic capture to drastically decrease fuel usage. In this paper, we outline a procedure to obtain a family of zero-fuel multi-moon orbiter trajectories, using a family of Keplerian maps derived by the first author previously. The maps capture well the dynamics of the full equations of motion; the phase space contains a connected chaotic zone where intersections between unstable resonant orbit manifolds provide the template for lanes of fast migration between orbits of different semimajor axes. Patched three body approach is used and the four body problem is broken down into two three-body problems, and the search space is considerably reduced by the use of properties of the Keplerian maps. We also introduce the notion of Switching Region where the perturbations due to the two perturbing moons are of comparable strength, and which separates the domains of applicability of the corresponding two Keplerian maps.

Ross, Shane D.↗
DOE OSTI2023

Maintaining Microstructure – The Path to Successful Technology Maturation in Fluidized Systems

Scaling circulating fluidized bed risers has been a point of contention for nearly a century. There have been numerous attempts to define various methodologies. These have all fallen short of providing a robust approach, primarily due to a lack of maintaining microstructure dynamics – a key factor in maintain interphase heat and mass transfer. Recent work by the author put forward an approach based up maintaining dynamic similarity at different scales by ensuring that the microstructure remained similar by maintaining statistical and chaotic parameters across the scale. That work relied on a dimensionless regime map in which the x-axis was defined by the ratio of the solids flux to the saturation carrying capacity. This latter property was only good within the range of the data and it became evident that extrapolation beyond the limits of the data could induce unrealistic conditions. Therefore, literature was reviewed to develop a better correlation for the saturation carrying capacity. In doing so, a critical diameter concept was developed and applied to the correlations for both Geldart Group A and B materials. The critical diameter for Geldart Group A and B materials is 0.2 m and 0.3 m, respectively. The paper then gives examples on how to use the scaling approach to maintain dynamic similarity across the scales.<br>

Breault, Ronald↗
DOE OSTI2023

Light-driven transitions in quantum paraelectrics

Motivated by recent experiments on pump-induced polar ordering in the quantum paraelectric SrTiO 3 , we study a driven phonon system close to a second-order phase transition. Analyzing its classical dynamics, we find that sufficiently strong driving leads to transitions into polar phases whose structures, determined by the light polarization, are not all accessible in equilibrium. In addition, for certain intensity profiles, we demonstrate the possibility of two-step transitions as a function of fluence. For even stronger field intensities, the possibility of period-doubling and chaotic behavior is demonstrated. Finally we develop a generalized formalism that allows us to consider quantum corrections to the classical dynamics in a systematic fashion. Furthermore, we predict a shift in the critical pump fluence due to quantum fluctuations with a characteristic dependence on the fluence increase rate that should be observable in experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Trotter Errors from Dynamical Structural Instabilities of Floquet Maps in Quantum Simulation

We study the behavior of errors in the quantum simulation of spin systems with long-range multibody interactions resulting from the Trotter-Suzuki decomposition of the time-evolution operator. We identify a regime where the Floquet operator underlying the Trotter decomposition undergoes sharp changes even for small variations in the simulation step size. This results in a time evolution operator that is very different from the dynamics generated by the targeted Hamiltonian, which leads to a proliferation of errors in the quantum simulation. These regions of sharp change in the Floquet operator, referred to as structural instability regions, appear typically at intermediate Trotter step sizes and in the weakly interacting regime, and are thus complementary to recently revealed quantum chaotic regimes of the Trotterized evolution [L. M. Sieberer et al. npj Quantum Inf. 5, 78 (2019); M. Heyl, P. Hauke, and P. Zoller, Sci. Adv. 5, eaau8342 (2019)]. We characterize these structural instability regimes in p-spin models, transverse-field Ising models with all-to-all p-body interactions, and analytically predict their occurrence based on unitary perturbation theory. We further show that the effective Hamiltonian associated with the Trotter decomposition of the unitary time-evolution operator, when the Trotter step size is chosen to be in the structural instability region, is very different from the target Hamiltonian, which explains the large errors that can occur in the simulation in the regions of instability. These results have implications for the reliability of near-term gate-based quantum simulators, and reveal an important interplay between errors and the physical properties of the system being simulated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sensitivity of Grass Fires Burning in Marginal Conditions to Atmospheric Turbulence

Abstract Atmospheric forcing and interactions between the fire and atmosphere are primary drivers of wildland fire behavior. The atmosphere is known to be a chaotic system that, although deterministic, is very sensitive to small perturbations to initial conditions. We assume that as a result of the tight coupling between fire and atmosphere; wildland fire behavior, in turn, should also be sensitive to perturbations in atmospheric initial conditions. Observations suggest that low intensity prescribed fire, in particular, is susceptible to small perturbations in the wind field, which can significantly alter fire spread. Here, we employ a computational fluid dynamics model of coupled fire‐atmosphere interactions to answer the question: How sensitive is fire behavior to small variations in atmospheric turbulence? We perform ensemble simulations of fires in homogenous grass fuels. The only difference between ensemble members is the state of the turbulent atmosphere provided to the model throughout the simulation. The atmospheric state is a function of the initial conditions applied at the start of the simulation and boundary conditions applied throughout the simulation. We find a wide range of outcomes, with area burned ranging from 2,212 to 11,236 m 2 (>400% change), driven primarily by sensitivity to initial conditions, with nonnegligible contributions from boundary condition variability during the initial 30 s of simulation. Our results highlight the need for ensemble simulations, especially when considering fire behavior in marginal burning conditions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

The route to chaos for the Kuramoto-Sivashinsky equation

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609... . Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Papageorgiou, Demetrios T.↗
NASA NTRS1991

The route to chaos for the Kuramoto-Sivashinsky equation

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609 .... Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Papageorgiou, Demetrios T.↗
DOE OSTI2024

Quantum fragmentation in the extended quantum breakdown model

We introduce a one-dimensional (1D) extended quantum breakdown model comprising a fermionic and a spin degree of freedom per site, and featuring a spatially asymmetric breakdown-type interaction between the fermions and spins. Furthermore, our model resembles the breakdown process of particles incident into a cloud chamber with nonzero quantum amplitudes of both exciting and not exciting the local vapor atoms. We analytically show that, in the absence of any magnetic field for the spins, the model exhibits Hilbert space fragmentation within each symmetry sector into exponentially many Krylov subspaces and hence displays nonthermal dynamics. Here, we demonstrate that the fragmentation naturally occurs in an entangled basis and thus provides an example of “quantum fragmentation.” Besides establishing the nature of fragmentation analytically, we also study the long-time behavior of the entanglement entropy and its deviation from the expected Page value as a probe of ergodicity in the system. Upon introducing a magnetic field for the spins, most of the Krylov subspaces merge and the model becomes chaotic. Finally, we study the effects of strong randomness on the system and observe behavior similar to that of many-body localized systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗