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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1991

Indications of low dimensionality in magnetospheric dynamics

Using three separate but related approaches, the question of whether the dynamic response of the magnetosphere to the solar wind input may be described by a low-order system of equations is examined. First, it is determined that the dimension of the subset (the attractor) in the high-dimensional magnetospheric phase space associated with the westward auroral electrojet (AL) index for some of the data sets compiled by Bargatze et al. (1985) is 4.0 + or - 0.2, seemingly independent of activity level. Second, direct modeling of the magnetosphere, considering the bulk properties of the tail plasma, leads to a system of equations that is similar to those previously reported as a dripping faucet model; here, the focus is specifically on the prediction of a natural frequency in this model. Finally, a peak is identified with the predicted frequency in power spectra of AL computed for intervals with both low and high activity. Peaks at other frequencies also appear in the spectra, and such resonances would be expected for a chaotic nonlinear oscillator. Combining these approaches it is concluded that at least some aspects of magnetospheric dynamics may be meaningfully modeled by low-dimensional sets of equations.

Roberts, D. A.↗
NASA NTRS1992

Simulation of particle mixing by turbulent convective flows on the connection machine

Mixing of particles by chaotic flow fields was simulated on the Connection Machine. Each cell was assigned to the processor, and the coordinates of particles residing on the cell were kept in the local memory of the processor. This approach implies the exchange between the local memories, when a particle moves from one cell to another. Approximately 10 exp 5 particles were injected into a time-dependent flow field obtained by solving the nonlinear system of partial differential equations describing turbulent thermal convection. The flow field was calculated on a CRAY, and data were transferred to a CM-200 through a high-speed HIPPI channel.

Malevsky, A. V.↗
NASA NTRS1990

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 140 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce thus that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism/ - Pauli matrices/. The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗
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.↗
NASA NTRS1991

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 240 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce, thus, that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism (- Pauli matrices). The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗
NASA NTRS1996

Heat and Momentum Transfer Studies in High Reynolds Number Wavy Films at Normal and Reduced Gravity Conditions

We examined the effect of the gas flow on the liquid film when the gas flows in the countercurrent direction in a vertical pipe at normal gravity conditions. The most dramatic effect of the simultaneous flow of gas and liquid in pipes is the greatly increased transport rates of heat, mass, and momentum. In practical situations this enhancement can be a benefit or it can result in serious operational problems. For example, gas-liquid flow always results in substantially higher pressure drop and this is usually undesirable. However, much higher heat transfer coefficients can be expected and this can obviously be of benefit for purposes of design. Unfortunately, designers know so little of the behavior of such two phase systems and as a result these advantages are not utilized. Due to the complexity of the second order boundary model as well as the fact that the pressure variation across the film is small compared to the imposed gas phase pressure, the countercurrent gas flow affect was studied for the standard boundary layer model. A different stream function that can compensate the shear stress affect was developed and this stream function also can predict periodic solutions. The discretized model equations were transformed to a traveling wave coordinate system. A stability analysis of these sets of equations showed the presence of a Hopf bifurcation for certain values of the traveling wave velocity and the shear stress. The Hopf celerity was increased due to the countercurrent shear. For low flow rate the increases of celerity are more than for the high flow rate, which was also observed in experiments. Numerical integration of a traveling wave simplification of the model also predicts the existence of chaotic large amplitude, nonperiodic waves as observed in the experiments. The film thickness was increased by the shear.

Balakotaiah, V.↗
NASA NTRS1992

Inhibition of electron thermal conduction by electromagnetic instabilities

Heat flux inhibition by electromagnetic instabilities in a hot magnetized plasma is investigated. Low-frequency electromagnetic waves become unstable due to anisotropy of the electron distribution function. The chaotic magnetic field thus generated scatters the electrons with a specific effective mean free path. Saturation of the instability due to wave-wave interaction, nonlinear scattering, wave propagation, and collisional damping is considered. The effective mean free path is found self-consistently, using a simple model to estimate saturation level and scattering, and is shown to decrease with the temperature gradient length. The results, limited to the assumptions of the model, are applied to astrophysical systems. For some interstellar clouds the instability is found to be important. Collisional damping stabilizes the plasma, and the heat conduction can be dominated by superthermal electrons.

Levinson, Amir↗
DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗
NASA NTRS1991

Non-Gaussian microwave background fluctuations from nonlinear gravitational effects

Whether the statistics of primordial fluctuations for structure formation are Gaussian or otherwise may be determined if the Cosmic Background Explorer (COBE) Satellite makes a detection of the cosmic microwave-background temperature anisotropy delta T(sub CMB)/T(sub CMB). Non-Gaussian fluctuations may be generated in the chaotic inflationary model if two scalar fields interact nonlinearly with gravity. Theoretical contour maps are calculated for the resulting Sachs-Wolfe temperature fluctuations at large angular scales (greater than 3 degrees). In the long-wavelength approximation, one can confidently determine the nonlinear evolution of quantum noise with gravity during the inflationary epoch because: (1) different spatial points are no longer in causal contact; and (2) quantum gravity corrections are typically small-- it is sufficient to model the system using classical random fields. If the potential for two scalar fields V(phi sub 1, phi sub 2) possesses a sharp feature, then non-Gaussian fluctuations may arise. An explicit model is given where cold spots in delta T(sub CMB)/T(sub CMB) maps are suppressed as compared to the Gaussian case. The fluctuations are essentially scale-invariant.

Salopek, D. S.↗
NASA NTRSDate not supplied

Ocean-Land Teleconnections and Chaotic Atmospheric Variability in Climate Models, with Implications for Real-World Behavior

This study introduces a unique approach to estimating a fundamental but unmeasurable property of Nature: ρo 2 , the fraction of the variance of a given meteorological variable that is determined by boundary forcing (specifically, SSTs) as opposed to chaotic atmospheric dynamics. Here, we focus on continental monthly-averaged near-surface air temperature, or T2M. Rather than determining the fractions of SST-explained T2M variance found for a specific atmospheric general circulation model (AGCM) and then assuming that these fractions also represent Nature, we perform an analysis that entails: (i) evaluating a specific AGCM’s ability to capture observed variability with an infinite number of ensemble members, (ii) quantifying the degree to which different modeling systems agree on SST-land teleconnections, and (iii) combining these results into our novel estimates of ρo 2 . The derived ρo 2 fields suggest, with caveats, that AGCMs tend to overestimate the fraction of explained T2M variance (the signal-to-total variance ratio) in the Tropics.

Randal D Koster↗
NASA NTRS1989

On the Nature of Navier-stokes Turbulence

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The solutions, separate exponentially with time, having a positive Liapunov exponent. Thus the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number is reduced. Several qualitatively different flows are noted. These are, fully chaotic, complex period, weakly chaotic, simple periodic, and fixed-point. Of these, only the fully chaotic flows are classified as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, is both time-dependent and aperiodic. Turbulent solutions are obtained in which energy cascades from large to small-scale motions. In general, the spectral energy transfer takes place between wavenumber bands that are considerably separated. The special transfer can occur either as a result of nonlinear turbulence self-interaction or by interaction of turbulence with mean gradients. Turbulent systems are compared with those studied in kinetic theory. The two types of systems are fundamentally different (continuous and dissipative as opposed to discrete and conservative), but there are similarities. For instance, both are nonlinear and show sensitive dependence on initial conditions. Also, the turbulent and molecular stress tensors are identical if the macroscopic velocities for the turbulent stress are replaced by molecular velocities.

Deissler, Robert G.↗
DOE OSTI2021

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2018

The Chaotic Long-Term X-Ray Variability of 4U 1705-44

The low-mass X-ray binary 4U1705-44 exhibits dramatic long-term X-ray time variability with a time-scale of several hundred days. The All-SkyMonitor (ASM) aboard the Rossi X-ray Timing Explorer (RXTE) and the Japanese Monitor of All-sky X-ray Image (MAXI) aboard the International Space Station together have continuously observed the source from 1995 December through 2014 May. The combined ASM-MAXI data provide a continuous time series over 50 times the length of the time-scale of interest. Topological analysis can help us identify 'fingerprints' in the phase space of a system unique to its equations of motion. The Birman-Williams theorem postulates that if such fingerprints are the same between two systems, then their equations of motion must be closely related. The phase-space embedding of the source light curve shows a strong resemblance to the double-welled non-linear Duffing oscillator. We explore a range of parameters for which the Duffing oscillator closely mirrors the time evolution of 4U1705-44. We extract low period, unstable periodic orbits from the 4U1705-44 and Duffing time series and compare their topological information. The Duffing and 4U1705-44 topological properties are identical, providing strong evidence that they share the same underlying template. This suggests that we can look to the Duffing equation to help guide the development of a physical model to describe the long-term X-ray variability of this and other similarly behaved X-ray binary systems.

Rossi X-ray Timing Explorer (RXTE)↗
DOE OSTI2020

Scalable Asynchronous Domain Decomposition Solvers

We discuss how parallel implementations of linear iterative solvers generally alternate between phases of data exchange and phases of local computation. Increasingly large problem sizes and more heterogeneous compute architectures make load balancing and the design of low latency network interconnects that are able to satisfy the communication requirements of linear solvers very challenging tasks. In particular, global communication patterns such as inner products become increasingly limiting at scale. We explore the use of asynchronous communication based on one-sided Message Passing Interface primitives in the context of domain decomposition solvers. In particular, a scalable asynchronous two-level Schwarz method is presented. We discuss practical issues encountered in the development of a scalable solver and show experimental results obtained on a state-of-the-art supercomputer system that illustrate the benefits of asynchronous solvers in load balanced as well as load imbalanced scenarios. Using the novel method, we can observe speedups of up to four times over its classical synchronous equivalent.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Temporal studies of compact galactic X-ray sources

Advances in X-ray astronomy due to temporal studies are discussed. Attention is given to X-ray temporal variability categorized into periodic, quasi-periodic, and chaotic. Approximately 95% of all known X-ray sources exhibit chaotic variability, in which there is no apparent pattern or periodicities. The orbits of seven binary X-ray sources are presented, along with empirical findings of neutron star masses. In addition, the interior structure of a supergiant star, HD 77581, is investigated by attempting to measure its apsidal motion. Graphs are provided for the power density spectrum of 4U 1626-67, GX301-2 Doppler delay data, and empirically determined error contours for the mass and radius of four companions stars in X-ray binary systems.

Rappaport, S.↗
NASA NTRS1988

The dynamics of coorbital satellite systems

The dynamical behavior of N coorbital satellites moving with the same average mean motion around a primary has been studied both analytically and by numerical integrations for N from 2 to 9 satellites. Simplified dynamical equations are used to determine the different stationary configurations and their local stability against infinitesimal perturbations. The ring of equally spaced identical satellites is found to be locally unstable for N of six or less, while for N from two to eight there exists another stable compact stationary configuration with separations of 60 deg or less between adjacent satellites. For N of seven or more the equally spaced configuration becomes locally stable, and for N of nine or more it is the only stationary configuration. The motion becomes chaotic for large-amplitude perturbations. The chaotic motion fills a restricted region of phase space whose outer boundary is determined by the maximum velocity curve.

Salo, H.↗