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 37 records · Page 2

On The Computational Capabilities of Physical Systems: Relationship With Conventional Computer Science - Part 2

In the first of this pair of papers, it was proven that there cannot be a physical computer to which one can properly pose any and all computational tasks concerning the physical universe. It was then further proven that no physical computer C can correctly carry out all computational tasks that can be posed to C. As a particular example, this result means that no physical computer that can, for any physical system external to that computer, take the specification of that external system's state as input and then correctly predict its future state before that future state actually occurs; one cannot build a physical computer that can be assured of correctly "processing information faster than the universe does". These results do not rely on systems that are infinite, and/or non-classical, and/or obey chaotic dynamics. They also hold even if one uses an infinitely fast, infinitely dense computer, with computational powers greater than that of a Turing Machine. This generality is a direct consequence of the fact that a novel definition of computation - "physical computation" - is needed to address the issues considered in these papers, which concern real physical computers. While this novel definition does not fit into the traditional Chomsky hierarchy, the mathematical structure and impossibility results associated with it have parallels in the mathematics of the Chomsky hierarchy. This second paper of the pair presents a preliminary exploration of some of this mathematical structure. Analogues of Chomskian results concerning universal Turing Machines and the Halting theorem are derived, as are results concerning the (im)possibility of certain kinds of error-correcting codes. In addition, an analogue of algorithmic information complexity, "prediction complexity", is elaborated. A task-independent bound is derived on how much the prediction complexity of a computational task can differ for two different reference universal physical computers used to solve that task, a bound similar to the "encoding" bound governing how much the algorithm information complexity of a Turing machine calculation can differ for two reference universal Turing machines. Finally, it is proven that either the Hamiltonian of our universe proscribes a certain type of computation, or prediction complexity is unique (unlike algorithmic information complexity), in that there is one and only version of it that can be applicable throughout our universe.

Wolpert, David H.↗

On the Computational Capabilities of Physical Systems: The Impossibility of Infallible Computation - Part 1

In this first of two papers, strong limits on the accuracy of physical computation are established. First it is proven that there cannot be a physical computer C to which one can pose any and all computational tasks concerning the physical universe. Next it is proven that no physical computer C can correctly carry out any computational task in the subset of such tasks that can be posed to C. This result holds whether the computational tasks concern a system that is physically isolated from C, or instead concern a system that is coupled to C. As a particular example, this result means that there cannot be a physical computer that can, for any physical system external to that computer, take the specification of that external system's state as input and then correctly predict its future state before that future state actually occurs; one cannot build a physical computer that can be assured of correctly 'processing information faster than the universe does'. The results also mean that there cannot exist an infallible, general-purpose observation apparatus, and that there cannot be an infallible, general-purpose control apparatus. These results do not rely on systems that are infinite, and/or non-classical, and/or obey chaotic dynamics. They also hold even if one uses an infinitely fast, infinitely dense computer, with computational powers greater than that of a Turing Machine. This generality is a direct consequence of the fact that a novel definition of computation - a definition of 'physical computation' - is needed to address the issues considered in these papers. While this definition does not fit into the traditional Chomsky hierarchy, the mathematical structure and impossibility results associated with it have parallels in the mathematics of the Chomsky hierarchy. The second in this pair of papers presents a preliminary exploration of some of this mathematical structure, including in particular that of prediction complexity, which is a 'physical computation analogue' of algorithmic information complexity. It is proven in that second paper that either the Hamiltonian of our universe proscribes a certain type of computation, or prediction complexity is unique (unlike algorithmic information complexity), in that there is one and only version of it that can be applicable throughout our universe.

Wolpert, David H.↗

Transposable elements as a molecular evolutionary force

This essay addresses the paradoxes of the complex and highly redundant genomes. The central theses developed are that: (1) the distinctive feature of complex genomes is the existence of epigenetic mechanisms that permit extremely high levels of both tandem and dispersed redundancy; (2) the special contribution of transposable elements is to modularize the genome; and (3) the labilizing forces of recombination and transposition are just barely contained, giving a dynamic genetic system of ever increasing complexity that verges on the chaotic.

Review↗

Dynamic Modeling of Off-Nominal Operation in Advanced Life Support Systems

System failures, off-nominal operation, or unexpected interruptions in processing capability can cause unanticipated instabilities in Advanced Life Support (ALS) systems, even long after they are repaired. Much current modeling assumes ALS systems are static and linear, but ALS systems are actually dynamic and nonlinear, especially when failures and off nominal operation are considered. Modeling and simulation provide a way to study the stability and time behavior of nonlinear dynamic ALS systems under changed system configurations or operational scenarios. The dynamic behavior of a nonlinear system can be fully explored only by computer simulation over the full range of inputs and initial conditions. Previous simulations of BIO-Plex in SIMULINK, a toolbox of Matlab, were extended to model the off-nominal operation and long-term dynamics of partially closed physical/chemical and bioregenerative life support systems. System nonlinearity has many interesting potential consequences. Different equilibrium points may be reached for different initial conditions. The system stability can depend on the exact system inputs and initial conditions. The system may oscillate or even in rare cases behave chaotically. Temporary internal hardware failures or external perturbations in ALS systems can lead to dynamic instability and total ALS system failure. Appropriate control techniques can restore reliable operation and minimize the effects of dynamic instabilities due to anomalies or perturbations in a life support system.

Jones, Harry↗

A Non-Intrusive Algorithm for Sensitivity Analysis of Chaotic Flow Simulations

We demonstrate a novel algorithm for computing the sensitivity of statistics in chaotic flow simulations to parameter perturbations. The algorithm is non-intrusive but requires exposing an interface. Based on the principle of shadowing in dynamical systems, this algorithm is designed to reduce the effect of the sampling error in computing sensitivity of statistics in chaotic simulations. We compare the effectiveness of this method to that of the conventional finite difference method.

Blonigan, Patrick J.↗

The applicability of Lyapunov characteristic numbers in the study of the stability of satellite orbits

It is pointed out that the Lyapunov Characteristic Numbers constitute a new tool for determining stability of trajectories of dynamical systems, or, even more generally, of solutions of systems of ordinary differential equations. In contrast with the characteristic exponents, which apply only to periodic solutions, the Lyapunov Characteristic Numbers apply to arbitrary nonperiodic solutions as well. A description is presented of the numerical experiments which have been made in order to investigate the practical value of the Lyapunov Characteristic Number and the Kolmogorov Entropy for the purpose of estimating the stability of trajectories and/or numerical integration methods in celestial mechanics. It is found that the Lyapunov Characteristic Numbers are extremely useful for the classification of the solutions of nonintegrable dynamical systems, especially in order to distinguish between quasi-periodic and chaotic solutions. However, the Lyapunov Characteristics Numbers do not appear to be useful for the purpose of evaluating numerical integration methods.

Broucke, R.↗

Solar system dynamics

The rotational dynamics of irregularly shaped satellites and the origin of Kirkwood Gaps are discussed. The chaotic tumbling of Hyperion and the anomalously low eccentricity of Deimos are examined. The Digital Orrery is used to explore the phase space of the ellipic restricted three body problem near the principal commensurabilities (2/1, 5/2, 3/1, and 3/2). The results for the 3/1 commensurability are in close agreement with those found earlier with the algebraic mapping method. Large chaotic zones are associated with the 3/1, 2/1 and 5/2 resonances, where there are gaps in the distribution of asteroids. The region near the 3/2 resonance, where the Hilda group of asteroids is located, is largely devoid of chaotic behavior. Thus, there is a qualitative agreement between the character of the motion and the distribution of asteroids.

Wisdom, Jack↗

Variations of the Milankovitch frequencies in time

The sensitivity of the amplitudes and frequencies in the development of the Earth's orbital and rotational elements involved in the astronomical theory of paleoclimates (eccentricity, obliquity, and climate precession), to the Earth-Moon distance and consequently to the length of the day and to the dynamical ellipticity of the Earth has been discussed for the last billions of years. The shortening of the Earth-Moon distance and of the length of the day, as well as the lengthening of the dynamical ellipticity of the Earth back in time induce a shortening of the fundamental astronomical periods for precession and obliquity. At the same time, the amplitudes of the different terms in the development of the obliquity are undergoing a relative enlargement of about 50 percent at 2 x 10(exp 9) yr BP but the independent term is increasing very weakly (less than 0.1 percent). In other words, the value of the obliquity, which lies within a range of 21.7 to 24.9 deg over the Quarternary was restricted to a range of 22.5 to 24.1 deg at 2 x 10(exp 9) yr BP. On the other hand, the amplitudes in the development of the climatic precession do not change. Moreover, these changes in the frequencies and amplitudes for both obliquity and climatic precession are larger for longer period terms. Finally, the periods in the eccentricity development are not influenced by the variation of the lunar distance. But the motion of the solar system, especially of the inner planets, was shown to be chaotic. It means that it is impossible to compute the exact motion of the planets over more than about 100 Myr, and the fundamental frequencies of the systems are not fixed quantities, but are slowly varying with time. As long as we consider the most important terms, the maximum deviation from the present-day value of the 19-kyr precessional period due to the chaotic motion of the solar system only does not reach more than a few tens of years around 80 Myr BP. Therefore the shortening of the obliquity and climatic precession periods is mostly driven by the change in the lunar distance and the consequent variations in the dynamical ellipticity of the Earth's angular speed. At first sight, the deviation in the period for the eccentricity can be neglected, as the chaotic behavior of the solar system implies a relative change of the main periods by less than 0.2 percent, 1.4 percent, and 1.9 percent respectively, this maximum change being achieved around 80 Myr BP. This implies, in particular, that the eccentricity periods for Quarternary climate studies may be considered more or less constant for pre-Quaternay times and equal to their Quaternary values.

Loutre, Marie-France↗

Some recent results of nonlinear acoustics in combustion chambers

The chief purpose of this paper is to provide a brief review of some recent work on nonlinear acoustics applicable to combustion instabilities, largely in the framework of an approximate analysis. Some connections will be made with the modern theory of nonlinear dynamical systems, including very recent and incomplete attempts by others to assess the possible chaotic behavior observed in laboratory tests.

Culick, F. E. C.↗

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗

Probable Disastrous Consequences of Collision Between Unknown Small (100 m) Asteroids with Known (Approximately 1 km) Near Earth Orbiting (NEO) Asteroids

The long-term stability of the Solar System is not well understood. Ironically its stability is taken for granted even though our knowledge of all the constituents [comets, asteroids. (The Asteroid Belt between Mars and Jupiter, Trojan Asteroids, Kuiper belt, Ort Cloud), planetoids, planets, moons, etc], and its long-term dynamics cannot be easily computed. At best one might say that the solar system is chaotic, but much of the time it seems to exists near a quasi-stationary state. An asteroid that passes near the Earth regularly returns with clock-like precision. Taking into account every known detail of its path through the solar system, its orbit is calculated forward thousands of years with no untoward calamity on the horizon. And then one day, this passive visitor slams into the Earth during a sunny afternoon picnic! Can this happen? Unfortunately, this is a real possibility in the ordinary history of the solar system. In fact our knowledge of the solar system in the small is sketchy, as will be pointed out. Events, which lie outside our awareness, can precipitate disasters that we may perceive when it's too late to launch effective counter measures. In this work, one such scenario is described and the direct consequences for the Earth are calculated.

Smalley, Larry↗

Mission analysis for the EM-1 CubeSats EQUULEUS and OMOTENASHI

EQUULEUS is a Lunar L2 orbiter and a 6-Unit CubeSat by JAXA and the University of Tokyo. OMOTENASHI is a 6-Unit CubeSat by JAXA, the world’s smallest Lunar lander. EQUULEUS and OMOTENASHI are among the 13 secondary payloads selected by NASA to be launched with Exploration Mission-1 in 2019. Despite their limited size and cost, EQUULEUS and OMOTENASHI are challenging missions, especially in terms of trajectory design and control. EQUULEUS exploits the Earth-Sun-Moon chaotic dynamics and enter a libration point orbit around the L2 of the Earth-Moon system, using a new water propulsion system with low thrust and little propellant. This “Orbit Control Experiment” is one of the main objectives of the mission. OMOTENASHI executes a semi-hard landing that requires breaking the spacecraft to a stop just a few-hundred meters above the Moon’s surface. Both missions present new and unique challenges, where the design of the nominal trajectory is mainly driven by the constrains on orbital control capabilities, and operational and robustness considerations. This paper presents the current baselines, and give an overview of the new techniques developed for their design.

Oshima, Kenta↗

Least Squares Shadowing Sensitivity Analysis of Chaotic Flow Around a Two-Dimensional Airfoil

Gradient-based sensitivity analysis has proven to be an enabling technology for many applications, including design of aerospace vehicles. However, conventional sensitivity analysis methods break down when applied to long-time averages of chaotic systems. This breakdown is a serious limitation because many aerospace applications involve physical phenomena that exhibit chaotic dynamics, most notably high-resolution large-eddy and direct numerical simulations of turbulent aerodynamic flows. A recently proposed methodology, Least Squares Shadowing (LSS), avoids this breakdown and advances the state of the art in sensitivity analysis for chaotic flows. The first application of LSS to a chaotic flow simulated with a large-scale computational fluid dynamics solver is presented. The LSS sensitivity computed for this chaotic flow is verified and shown to be accurate, but the computational cost of the current LSS implementation is high.

Blonigan, Patrick J.↗

Determining Chaotic Instabilities in Mechanical Systems

Theoretical developments enable suppression of chaotic structual motions. Theory enables prediction, avoidance, and suppression of chaotic vibrations in structures, especially using dynamic feedback stabilization. In new formulation, motion both repeatable and predictable.

Zak, M. A.↗

Chaos in the sunspot cycle - Analysis and prediction

The variability of solar activity over long time scales, given semiquantitatively by measurements of sunspot numbers, is examined as a nonlinear dynamical system. First, a discussion of the data set used and the techniques utilized to reduce the noise and capture the long-term dynamics inherent in the data is presented. Subsequently, an attractor is reconstructed from the data set using the method of time delays. The reconstructed attractor is then used to determine both the dimension of the underlying system and also the largest Lyapunov exponent, which together indicate that the sunspot cycle is indeed chaotic and also low dimensional. In addition, recent techniques of exploiting chaotic dynamics to provide accurate, short-term predictions are utilized in order to improve upon current forecasting methods and also to place theoretical limits on predictability extent. The results are compared to chaotic solar-dynamo models as a possible physically motivated source of this chaotic behavior.

Mundt, Michael D.↗

The kinematics of buoyancy induced mixing

A model problem to study mixing driven by buoyancy flow fields due to a time dependent body force which consists of a steady and an oscillatory component is considered. Flow fields generated by a time dependent body force, which consists of a steady and an oscillatory component from vibrations or g-jitter, mix two fluids inside a cavity by stretching and folding the interface into various morphological patterns. These patterns show two basic structures depending on the dominant components of the body force. These structures are whirls and tendrils, corresponding to the dominance of either the steady or oscillatory component, respectively. A combination of these structures also occurs in the interaction region. In this region the effect of the steady component opposes the oscillatory component and the result is to smooth the interface, and prevent breakup of the interface. For the parametric range considered, two basic flow regimes occur: convective and chaotic. In the convective regime, the morphological patterns are topological. In the chaotic regime, mixing occurs by a repetitive sequence of bubble formation, necking, and breakup which spreads the concentration field throughout the cavity. The length stretch of the interface increases exponentially with time. The dynamical system characteristics show that in the convective regime, at a point in space, the flow field is oscillatory with elliptical phase space trajectories, and its power spectrum indicates that the flow field responds to the corresponding input frequency. Whereas, in the chaotic regime, the flow field is aperiodic, the phase space trajectories show irregular patterns, and the power spectrum shows a broadband distribution.

Duval, Walter M. B.↗

Low-dimensional chaos in magnetospheric activity from AE time series

The magnetospheric response to the solar-wind input, as represented by the time-series measurements of the auroral electrojet (AE) index, has been examined using phase-space reconstruction techniques. The system was found to behave as a low-dimensional chaotic system with a fractal dimension of 3.6 and has Kolmogorov entropy less than 0.2/min. These indicate that the dynamics of the system can be adequately described by four independent variables, and that the corresponding intrinsic time scale is of the order of 5 min. The relevance of the results to magnetospheric modeling is discussed.

Vassiliadis, D. V.↗