Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chaotic 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic “active turbulence.” Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Overall, our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hybrid gels via bulk interfacial complexation of supramolecular polymers and polyelectrolytes

Hierarchical self-assembly leading to organized supramolecular structures across multiple length scales has been of great recent interest. Earlier work from our laboratory reported the complexation of peptide amphiphile (PA) supramolecular polymers with oppositely charged polyelectrolytes into a single solid membrane at a macroscopic interface. We report here the formation of bulk gels with many internal interfaces between the covalent and supramolecular polymer components formed by the rapid chaotic mixing of solutions, one containing negatively charged PA nanofibers and the other the positively charged biopolymer chitosan. We found that formation of a contact layer at the interface of the solutions locks the formation of hydrogels with lamellar microstructure. The nanofiber morphology of the supramolecular polymer is essential to this process since gels do not form when solutions of supramolecular assemblies form spherical micelles. Here, we found that rheological properties of the gels can be tuned by changing the relative amounts of each component. Furthermore, both positively and negatively charged proteins are easily encapsulated within the contact layer of the gel, which provides an interesting biomedical function for these systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An SYK-inspired model with density–density interactions: Spectral & wave function statistics, Green’s function and phase diagram

The Sachdev–Ye–Kitaev (SYK) model is a rare example of a strongly-interacting system that is analytically tractable. Tractability arises because the model is largely structureless by design and therefore artificial: while the interaction is restricted to two-body terms, interaction matrix elements are “randomized” and therefore the corresponding interaction operator does not commute with the local density. Unlike conventional density–density-type interactions, the SYK-interaction is, in this sense, not integrable. We here investigate a variant of the (complex) SYK model, which restores this integrability. It features a randomized single-body term and a density–density-type interaction. We present numerical investigations suggesting that the model exhibits two integrable phases separated by several intermediate phases including a chaotic one. The chaotic phase carries several characteristic SYK-signatures including in the spectral statistics and the frequency scaling of the Green’s function and therefore should be adiabatically connected to the non-Fermi liquid phase of the original SYK model. Thus, our model Hamiltonian provides a bridge from the SYK-model towards microscopic realism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shaker-structure interaction modeling and analysis for nonlinear force appropriation testing

Nonlinear force appropriation is an extension of its linear counterpart where sinusoidal excitation is applied to a structure with a modal shaker and phase quadrature is achieved between the excitation and response. While a standard practice in modal testing, modal shaker excitation has the potential to alter the dynamics of the structure under test. Previous studies have been conducted to address several concerns, but this work specifically focuses on a shaker-structure interaction phenomenon which arises during the force appropriation testing of a nonlinear structure. Under pure-tone sinusoidal forcing, a nonlinear structure may respond not only at the fundamental harmonic but also potentially at sub- or superharmonics, or it can even produce aperiodic and chaotic motion in certain cases. Shaker-structure interaction occurs when the response physically pushes back against the shaker attachment, producing non-fundamental harmonic content in the force measured by the load cell, even for pure tone voltage input to the shaker. This work develops a model to replicate these physics and investigates their influence on the response of a nonlinear normal mode of the structure. Experimental evidence is first provided that demonstrates the generation of harmonic content in the measured load cell force during a force appropriation test. This interaction is replicated by developing an electromechanical model of a modal shaker attached to a nonlinear, three-mass dynamical system. Several simulated experiments are conducted both with and without the shaker model in order to identify which effects are specifically due to the presence of the shaker. Finally, the results of these simulations are then compared to the undamped nonlinear normal modes of the structure under test to evaluate the influence of shaker-structure interaction on the identified system’s dynamics.

42 ENGINEERING↗
DOE OSTI2022

Cretaceous sequence stratigraphy of the northern Baltimore Canyon Trough: Implications for basin evolution and carbon storage

We evaluate the Cretaceous stratigraphy and carbon sequestration potential of the northern Baltimore Canyon Trough (NBCT) using >10,000 km of multi-channel seismic profiles integrated with geophysical logs, biostratigraphy, and lithology from 29 offshore wells. We identify and map six sequences resolved primarily at the stage level. Accommodation was dominated by thermal and non-thermal subsidence, though sequence boundaries correlate with regional and global sea-level changes, and the record is modified by igneous intrusion, active faulting, and changes in sediment supply and sources. Our stratigraphic maps illustrate a primary southern (central Appalachian) Early Cretaceous source that migrated northward during the Aptian and Albian. During the Cenomanian, sedimentation rates in the NBCT increased and depocenters shifted northward and landward. We show that deposition occurred in three phases: (1) earlier Cretaceous paleoenvironments were primarily terrestrial indicated by variable amplitude, chaotic seismic facies, serrated gamma logs, and heterolithic sandstones and mudstones with terrestrial microfossils; (2) the Albian to Cenomanian was dominated by deltaic paleoenvironments indicated by blocky, funnel-shaped, gamma-ray logs and clinoforms characterized by continuous high-amplitude seismic reflections with well-defined terminations; and (3) the Cenomanian and younger was marine shelf, inferred from mudstoneprone lithologies, peak gamma-ray values in well logs, and foraminiferal evidence. Long-term transgression and maximum water depths at the Cenomanian/Turonian boundary correlative with Ocean Anoxic Event 2 were followed by a regression and relative sea-level fall. We show that porous and permeable sandstones of three Aptian to Cenomanian highstand systems tracts are high-volume reservoirs for supercritical CO 2 storage that are confined by overlying deep water mudstones.

58 GEOSCIENCES↗
DOE OSTI2023

Consequences and benefits of utilizing continuous vibro-impact representations in constrained pipeline conveying fluid systems

Here, the effectiveness of continuous vibro-impact forcing representations for the cantilevered pipe that conveys fluid is explored and analyzed. The previously accepted forcing model utilizing a smoothened trilinear spring is estimated using three continuous forcing representations, namely, polynomial, rational polynomial, and hyperbolic tangent. The accuracy of the estimated forcing functions is investigated and analyzed by calculating the root mean square error, and bifurcation diagrams are generated and compared to the nominal system. Additionally, the dynamic response of the system is further characterized using Poincare maps, power spectra, and basins of attraction. Once all continuous forcing representations are analyzed and compared to the nominal system, the computational cost of each method is examined, and further limitations of the hyperbolic tangent method are discovered. It is proved that the hyperbolic tangent forcing representation most accurately captures the dynamic response of the pipeline, and the least accurate representation is the rational polynomial representation. Additionally, considerable computational cost is saved when employing the hyperbolic tangent representation compared to the discontinuous representation.

42 ENGINEERING↗
DOE OSTI2022

Scaling CFB risers: Maintaining microstructure dynamics

Over the years there have been numerous approaches to develop a scaleup methodology for circulating fluidized beds. Many of these past approaches have relied on macroscopic definitions of similitude to define dimensionless groups such as Reynolds number, Froude number, density ratios, and dimensional ratios involving different length scales. These approaches have demonstrated various levels of success, but often fail when used outside the bounds of experimental datasets. The issue with these methods is that they rely on macroscopic similitude with the hopes that those conditions would provide microscopic similitude as it does in single phase systems where the Buckingham Pi techniques have given reasonable success. Two of the main reasons for failure are 1) changing fluidization regimes and 2) changing material properties that result in the material falling into a different Geldart group. Both are certain to change the microstructure (local instantaneous volume solids fraction) and changes to the microstructure affect the interphase transport properties of the reactor system, namely the mass transfer and heat transfer. Here, this work presents a new approach towards scaling based upon maintaining similitude at the microscale. It reviews a novel dimensionless flow regime map developed from analysis of the deterministic chaos parameters and higher order moments to ensure microstructure similitude. The work then examines four cases from the literature to assess the scaling approach used: (1) same particle-different scale risers for Geldart Group B particles, (2) Glicksman scaling of atmospheric circulating fluidized bed combustor, (3) same particle-different scale risers for Geldart Group A particles and (4) Glicksman scaling for pressurized gasification. General findings are that none of the scaling approaches produced directly applicable microstructure similar behavior based upon the chaotic and higher moment analysis. The positive take away is that the approach can define scale conditions that will provide similar microstructure that can be us to predict transport behavior.

42 ENGINEERING↗
DOE OSTI2025

Active-noise-induced dynamic clustering of passive colloidal particles

Active fluids generate spontaneous, often chaotic mesoscale flows. Harnessing these flows to drive soft materials embedded within an active fluid into structures with controlled length scales and lifetimes is a key challenge at the interface between the fields of active matter and nonequilibrium self-assembly. Here, we present a simple and efficient computational approach to model soft materials advected by active fluids, by simulating particles moving in a spatiotemporally correlated noise field. To illustrate our approach, we simulate the dynamical self-organization of repulsive colloids within such an active noise field. The colloids form structures whose sizes and dynamics can be tuned by the correlation time and length of the active fluid, and range from small rotating droplets to clusters with internal flows and system-spanning sizes that vastly exceed the active correlation length. Our results explain how the interplay between active fluid time and length scales and emergent driven assembly can be used to rationally design functional assemblies. More broadly, our approach can be used to efficiently simulate diverse active fluids and other systems with spatiotemporally correlated noise.

Brownian dynamics↗
DOE OSTI2024

Passive Islanding Detection of Inverter-Based Resources in a Noisy Environment

Islanding occurs when a load is energized solely by local generators and can result in frequency and voltage instability, changes in current, and poor power quality. Poor power quality can interrupt industrial operations, damage sensitive electrical equipment, and induce outages upon the resynchronization of the island with the grid. This study proposes an islanding detection method employing a Duffing oscillator to analyze voltage fluctuations at the point of common coupling (PCC) under a high-noise environment. Unlike existing methods, which overlook the noise effect, this paper mitigates noise impact on islanding detection. Power system noise in PCC measurements arises from switching transients, harmonics, grounding issues, voltage sags and swells, electromagnetic interference, and power quality issues that affect islanding detection. Transient events like lightning-induced traveling waves to the PCC can also introduce noise levels exceeding the voltage amplitude by more than seven times, thus disturbing conventional detection techniques. The noise interferes with measurements and increases the nondetection zone (NDZ), causing failed or delayed islanding detection. The Duffing oscillator nonlinear dynamics enable detection capabilities at a high noise level. The proposed method is designed to detect the PCC voltage fluctuations based on the IEEE standard 1547 through the Duffing oscillator. For the voltages beyond the threshold, the Duffing oscillator phase trajectory changes from periodic to chaotic mode and sends an islanded operation command to the inverter. The proposed islanding detection method distinguishes switching transients and faults from an islanded operation. Experimental validation of the method is conducted using a 3.6 kW PV setup.

Energy & Fuels↗
DOE OSTI2022

Reconfigurable training and reservoir computing in an artificial spin-vortex ice via spin-wave fingerprinting

Strongly interacting artificial spin systems are moving beyond mimicking naturally occurring materials to emerge as versatile functional platforms, from reconfigurable magnonics to neuromorphic computing. Typically, artificial spin systems comprise nanomagnets with a single magnetization texture: collinear macrospins or chiral vortices. Here, by tuning nanoarray dimensions we have achieved macrospin–vortex bistability and demonstrated a four-state metamaterial spin system, the ‘artificial spin-vortex ice’ (ASVI). ASVI can host Ising-like macrospins with strong ice-like vertex interactions and weakly coupled vortices with low stray dipolar field. Vortices and macrospins exhibit starkly differing spin-wave spectra with analogue mode amplitude control and mode frequency shifts of Δf = 3.8 GHz. The enhanced bitextural microstate space gives rise to emergent physical memory phenomena, with ratchet-like vortex injection and history-dependent non-linear fading memory when driven through global magnetic field cycles. We employed spin-wave microstate fingerprinting for rapid, scalable readout of vortex and macrospin populations, and leveraged this for spin-wave reservoir computation. ASVI performs non-linear mapping transformations of diverse input and target signals in addition to chaotic time-series forecasting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Vibro-impact analysis and characterization of pipeline conveying fluids with multi-segmented motion-limiting constraints

Previous studies of the cantilevered pipeline conveying fluid system have included motion-limiting constraints in the form of trilinear springs. While this is desirable in experimental scenarios, it may not be representative of real-world applications. Therefore, here, this study focuses on multi-segmented motion-limiting constraints. As this type of motion-limiting constraint has not been investigated with a cantilevered pipeline system, a wide variety of outer and inner constraint stiffness and constraint gap sizes are investigated in this study to gain a comprehensive understanding of how the multi-segmented constraints affect the dynamics of the cantilevered pipeline. In this effort, bifurcation diagrams, phase portraits, Poincare maps, time histories, and power spectra are used to investigate the dynamics of the system, and the fluid flow speeds where dynamic characteristics are considered. In general, it is found that critical flow speeds like when the pipe sticks in the constraints are reduced as the constraint stiffnesses are increased. Additionally, the sticking flow speed occurred at lower flow speeds as the gap sizes of the inner and outer constraints decrease, and a larger constraint offset results in a smaller inner gap size leading to critical behaviors occurring at earlier flow speeds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Secular Resonances in Nonrestricted Hierarchical Triple Systems

In this paper, the averaged Hamiltonian of a nonrestricted hierarchical triple system truncated at the third order is investigated. First, each secular resonant term is studied. For the well-studied secular quadrupole theory, it is analytically reformulated in a different manner in our work. The resonance width is numerically determined and displayed on the √(1−e{sub 1}{sup 2})−√(1−e{sub 2}{sup 2}) plane (also denoted as the e-tilde {sub 1}− e-tilde {sub 2} plane). In terms of the octupole terms, we show that for a near-planar configuration of the system, considerable variations of both the eccentricities of the inner and outer orbits can be generated by a single resonant term. The resonance width for every secular resonant angle from the octupole terms is also numerically determined and displayed on the e-tilde {sub 1}− e-tilde {sub 2} plane. The results show that an orbit flip with a near-perpendicular initial mutual inclination is possible for each secular resonance. By displaying the resonance widths of different resonant terms on the same e-tilde {sub 1}− e-tilde {sub 2} plane, we intuitively show the overlap of different secular resonances. Then, the full averaged Hamiltonian with both quadrupole and octupole terms is investigated using the Poincaré surface of section, with a special focus on the orbit flip. For the cases we exploited, we find that the near-planar flip of the inner orbit can be either regular or chaotic while the outer orbit flip is generally chaotic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Anomalous losses of energetic particles in the presence of an oscillating radial electric field in fusion plasmas

The confinement of energetic particles in nuclear fusion devices is studied in the presence of an oscillating radial electric field and an axisymmetric magnetic equilibrium. Here, it is shown that, despite the poloidal and toroidal symmetries, initially integrable orbits turn into chaotic regions that can potentially intercept the wall of the tokamak, leading to particle losses. It is observed that the losses exhibit algebraic time decay different from the expected exponential decay characteristic of radial diffusive transport. A dynamical explanation of this behaviour is presented, within the continuous time random walk theory. The central point of the analysis is based on the fact that, contrary to the radial displacement, the poloidal angle is not bounded and a proper statistical analysis can therefore be made, showing for the first time that energetic particle transport can be super-diffusive in the poloidal direction and characterised by asymmetric poloidal displacement. The connection between poloidal and radial positions ensured by the conservation of the toroidal canonical momentum, implies that energetic particles spend statistically more time in the inner region of the tokamak than in the outer one, which explains the observed algebraic decay. This indicates that energetic particles might be efficiently slowed down by the thermal population before leaving the system. Also, the asymmetric transport reveals a new possible mechanism of self-generation of momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Poor confinement in stellarators at high energy

High-energy particle resonances can modify particle distributions and even cause significant particle loss. Resonances can be present in any toroidal confinement device and can easily be found numerically. Many stellarators have weak magnetic shear so that large islands and large chaotic regions can be produced by resonant perturbations with small amplitudes. Additionally, while the choice of the field line helicity profile in the plasma can limit the presence of resonances at low particle energy, the resonance location is energy-dependent, and they can move into the plasma at higher energy. If resonances match the toroidal variation of the equilibrium, they can produce wide islands in the phase space of orbits even in the absence of perturbations due to instabilities. These islands increase in size with particle energy and can seriously affect the confinement of high-energy ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗