Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chaos and regularity”

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

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

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

Computational Study of Chaotic and Ordered Solutions of the Kuramoto-Sivashinsky Equation

We report the results of extensive numerical experiments on the Kuramoto-Sivashinsky equation in the strongly chaotic regime as the viscosity parameter is decreased and increasingly more linearly unstable modes enter the dynamics. General initial conditions are used and evolving states do not assume odd-parity. A large number of numerical experiments are employed in order to obtain quantitative characteristics of the dynamics. We report on different routes to chaos and provide numerical evidence and construction of strange attractors with self-similar characteristics. As the 'viscosity' parameter decreases the dynamics becomes increasingly more complicated and chaotic. In particular it is found that regular behavior in the form of steady state or steady state traveling waves is supported amidst the time-dependent and irregular motions. We show that multimodal steady states emerge and are supported on decreasing windows in parameter space. In addition we invoke a self-similarity property of the equation, to show that these profiles are obtainable from global fixed point attractors of the Kuramoto-Sivashinsky equation at much larger values of the viscosity.

Smyrlis, Yiorgos S.↗
NASA NTRS1999

Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions

We examine the nonlinear stability of the Wisdom-Holman (WH) symplectic mapping applied to the integration of perturbed, highly eccentric (e-0.9) two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories for this class of problems, unless the step size chosen 1s small enough that PeriaPse is always resolved, in which case the method is generically stable. This 'radial orbit instability' persists even for weakly perturbed systems. Using the Stark problem as a fiducial test case, we investigate the dynamical origin of this instability and argue that the numerical chaos results from the overlap of step-size resonances; interestingly, for the Stark-problem many of these resonances appear to be absolutely stable. We similarly examine the robustness of several alternative integration methods: a time-regularized version of the WH mapping suggested by Mikkola; the potential-splitting (PS) method of Duncan, Levison, Lee; and two original methods incorporating approximations based on Stark motion instead of Keplerian motion. The two fixed point problem and a related, more general problem are used to conduct a comparative test of the various methods for several types of motion. Among the algorithms tested, the time-transformed WH mapping is clearly the most efficient and stable method of integrating eccentric, nearly Keplerian orbits in the absence of close encounters. For test particles subject to both high eccentricities and very close encounters, we find an enhanced version of the PS method-incorporating time regularization, force-center switching, and an improved kernel function-to be both economical and highly versatile. We conclude that Stark-based methods are of marginal utility in N-body type integrations. Additional implications for the symplectic integration of N-body systems are discussed.

Rauch, Kevin P.↗
DOE OSTI2020

Solving Coupled Cluster Equations by the Newton Krylov Method

We describe using the Newton Krylov method to solve the coupled cluster equation. The method uses a Krylov iterative method to compute the Newton correction to the approximate coupled cluster amplitude. The multiplication of the Jacobian with a vector, which is required in each step of a Krylov iterative method such as the Generalized Minimum Residual (GMRES) method, is carried out through a finite difference approximation, and requires an additional residual evaluation. The overall cost of the method is determined by the sum of the inner Krylov and outer Newton iterations. We discuss the termination criterion used for the inner iteration and show how to apply pre-conditioners to accelerate convergence. We will also examine the use of regularization technique to improve the stability of convergence and compare the method with the widely used direct inversion of iterative subspace (DIIS) methods through numerical examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Ergodic and nonergodic many-body dynamics in strongly nonlinear lattices

The study of nonlinear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi et al. [Los Alamos Scientific Laboratory Report No. LA-1940, 1955 (unpublished)]. Here, we introduce a family of such systems which consist of chains of N harmonically coupled particles with the nonlinearity introduced by confining the motion of each individual particle to a box or stadium with hard walls. The stadia are arranged on a one-dimensional lattice but they individually do not have to be one dimensional, thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of N . Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh Gordon theory, we do not see any evidence for the kind of nonergodicity famously seen in the work of Fermi et al. Finally, we examine the chain with particles confined to two-dimensional stadia where the individual stadium is already chaotic and find a much more chaotic phase space at small system sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Nonlinear problems in flight dynamics

A comprehensive framework is proposed for the description and analysis of nonlinear problems in flight dynamics. Emphasis is placed on the aerodynamic component as the major source of nonlinearities in the flight dynamic system. Four aerodynamic flows are examined to illustrate the richness and regularity of the flow structures and the nature of the flow structures and the nature of the resulting nonlinear aerodynamic forces and moments. A framework to facilitate the study of the aerodynamic system is proposed having parallel observational and mathematical components. The observational component, structure is described in the language of topology. Changes in flow structure are described via bifurcation theory. Chaos or turbulence is related to the analogous chaotic behavior of nonlinear dynamical systems characterized by the existence of strange attractors having fractal dimensionality. Scales of the flow are considered in the light of ideas from group theory. Several one and two degree of freedom dynamical systems with various mathematical models of the nonlinear aerodynamic forces and moments are examined to illustrate the resulting types of dynamical behavior. The mathematical ideas that proved useful in the description of fluid flows are shown to be similarly useful in the description of flight dynamic behavior.

Chapman, G. T.↗
NASA NTRS2002

The Global Geophysical Fluids Center of IERS (and its Special Bureau for Mantle)

The Global Geophysical Fluids Center (GGFC) was established by the International Earth Rotation Service (IERS) on IERS's 10th anniversary day January 1, 1998, in an effort to expand IERS's services to the scientific community. Under the GGFC, eight Special Bureaus (SB) have been selected, each to be responsible for research and data service activities pertaining to mass transports and related geophysical processes in specific components of the Earth system, or "global geophysical fluids," including the atmosphere, oceans, solid Earth, core, and geophysical processes of gravity, loading, tides and hydrological cycles. GGFC and the SBs have the responsibility of supporting, facilitating, and providing services to the worldwide research community, in areas related to the variations in Earth rotation, gravity field and geocenter that are caused by mass transport in the global geophysical fluids. These minute variations have been observed by various space geodetic techniques, as effective remote sensing tools, with ever increasing precision/accuracy and temporal/spatial resolution. The GGFC and SBs have organized dedicated workshops and special sessions at international conferences, published articles, and held regular business meetings. The SBs also maintain individual website for data services and information exchanges. See URL . In particular, the SB for Mantle focuses on large-scale mass redistributions that occur in the mantle in association with various dynamic processes, including seismic activities, the post-glacial rebound, and mantle convections.

Chao, Benjamin F.↗
NASA NTRS2019

Approximate Entropy and Sample Entropy: A Comprehensive Tutorial

Approximate Entropy and Sample Entropy are two algorithms for determining the regularity of series of data based on the existence of patterns. Despite their similarities, the theoretical ideas behind those techniques are different but usually ignored. This paper aims to be a complete guideline of the theory and application of the algorithms, intended to explain their characteristics in detail to researchers from different fields. While initially developed for physiological applications, both algorithms have been used in other fields such as medicine, telecommunications, economics or Earth sciences. In this paper, we explain the theoretical aspects involving Information Theory and Chaos Theory, provide simple source codes for their computation, and illustrate the techniques with a step by step example of how to use the algorithms properly. This paper is not intended to be an exhaustive review of all previous applications of the algorithms but rather a comprehensive tutorial where no previous knowledge isrequired to understand the methodology.

Delgado-Bonal, Alfonso↗
DOE OSTI2021

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lie-algebraic Kähler sigma models with U(1) isotropy

We discuss various questions that emerge in connection with the Lie-algebraic deformation of the sigma model in two dimensions. First, we supersymmetrize the original model endowing it with the minimal and extended supersymmetries. Then we derive the general hypercurrent anomaly in both cases. In the latter case this anomaly is one-loop but is somewhat different from the standard expressions one can find in the literature because the target manifold is nonsymmetric. We also show how to introduce the twisted masses and the term, and study the Bogomol’nyi–Prasad–Sommerfield equation for instantons, in particular the value of the topological charge. Then we demonstrate that the second loop in the function of the nonsupersymmetric Lie-algebraic sigma model is due to an infrared effect. To this end we use a supersymmetric regularization. We also conjecture that the above statement is valid for higher loops too, similar to the parallel phenomenon in four-dimensional super-Yang-Mills. In the second part of the paper we develop a special dimensional reduction—namely, starting from the two-dimensional Lie-algebraic model we arrive at a quasi-exactly solvable quantum-mechanical problem of the Lamé type.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Time-dependent relaxation of observables in complex quantum systems

Abstract We consider time-dependent relaxation of observables in quantum systems of chaotic and regular type. Using statistical arguments and exact numerical solutions we show that the spread of the initial wave function in the Hilbert space and the main characteristics of evolution of observables have certain generic features. The study compares examples of regular dynamics, a completely chaotic case of the Gaussian orthogonal ensemble, a bosonic system with random interactions, and a fully realistic case of the time evolution of various initial non-stationary states in the nuclear shell model. In the case of the Gaussian orthogonal ensemble we show that the survival probability obtained analytically also fully defines the relaxation timescale of observables. This is not the case in general. Using the realistic nuclear shell model and the quadrupole moment as an observable we demonstrate that the relaxation time is significantly longer than defined by the survival probability of the initial state. The full analysis does not show the presence of an analog of the Lyapunov exponent characteristic for examples of classical chaos.

Volya, Alexander (ORCID:0000000217656466)↗
NASA NTRS1989

Subtropical middle atmosphere dynamics observed by the Chung Li radar

The Chung Li Radar (24.91 N; 121.24 E) has been operating since 1986. A five beam observational configuration was used on a regular basis to study the various dynamics processes in the atmosphere-lower stratosphere height region. Due to its geographical location, the annual Typhoon and Mei-Yu seasons provide good opportunities to study the various interesting dynamic processes such as instabilities, generation of gravity waves, wave mean field interaction, etc. Three dimensional air motions due to these fronts are presented. Special cases of gravity wave generation, propagation and their effects on the turbulent layers are discussed.

Fu, I. J.↗
NASA NTRS1984

Origin and Evolution of the Uranian and Neptunian Satellites: Some Dynamical Considerations

The satellite system of Neptune is so irregular that some formation mechanism or subsequent dynamical evolution that sets it apart from the other, more regular, satellite systems is obviously indicated. McKinnon argued that satellite capture is the most likely possibility and has shown that tidal circularization of Triton's presumably highly eccentric initial orbit probably resulted in melting of the satellite's interior. The satellite system of Uranus, although somewhat bland, also has a number of special features that indicate an interesting dynamical history. These include the anomalously high orbital inclination of Miranda and the probable coexistence of small satellites and narrow rings inside the plant's Roche limit. The possibility that orbital evolution due to tidal dissipation is involved in both of these phenomena is discussed. Other topics discussed are: the origin of rings; the formation of coorbital satellites; the lack of stable orbit-orbit resonances in the Uranian satellite system; and chaos, tidal heating and the shapes of Miranda and Ariel.

Dermott, S. F.↗
DOE OSTI2023

Insights into ENSO Diversity from an Intermediate Coupled Model. Part II: Role of Nonlinear Dynamics and Stochastic Forcing

Abstract In this study, we investigate how a single leading linear El Niño–Southern Oscillation (ENSO) mode, as studied in Part I, leads to the irregular coexistence of central Pacific (CP) and eastern Pacific (EP) ENSO, a phenomenon known as ENSO spatiotemporal diversity. This diversity is fundamentally generated by deterministic nonlinear pathways to chaos via the period-doubling route and, more prevailingly, the subharmonic resonance route with the presence of a seasonally varying basic state. When residing in the weakly nonlinear regime, the coupled system sustains a weak periodic oscillation with a mixed CP/EP pattern as captured by the linear ENSO mode. With a stronger nonlinearity effect, the ENSO behavior experiences a period-doubling bifurcation. The single ENSO orbit splits into coexisting CP-like and EP-like ENSO orbits. A sequence of period-doubling bifurcation results in an aperiodic oscillation featuring irregular CP and EP ENSO occurrences. The overlapping of subharmonic resonances between ENSO and the seasonal cycle allows this ENSO irregularity and diversity to be more readily excited. In the strongly nonlinear regime, the coupled system is dominated by regular EP ENSO. The deterministic ENSO spatiotemporal diversity is thus confined to a relatively narrow range corresponding to a moderately unstable ENSO mode. Stochastic forcing broadens this range and allows ENSO diversity to occur when the ENSO mode is weakly subcritical. A close relationship among a weakened mean zonal temperature gradient, stronger ENSO activity, and more (fewer) occurrences of EP (CP) ENSO is noted, indicating that ENSO–mean state interaction may yield ENSO regime modulations on the multidecadal time scale.

Geng, Licheng↗
NASA NTRS2010

Airborne Dust, "The Good Guy or the Bad Guy": How Much do We Know?

Processes in generating, transporting, and dissipating the airborne dust particles are global phenomena -African dust regularly reaching the Alps; Asian dust seasonally crossing the Pacific into North America, and ultimately the Atlantic into Europe. One of the vital biogeochemical roles dust storms play in Earth's ecosystem is routinely mobilizing mineral dust, as a source of iron, from deserts into oceans for fertilizing the growth of phytoplankton -the basis of the oceanic food chain. Similarly, these dust-laden airs also supply crucial nutrients for the soil of tropical rain forests, the so-called womb of life that hosts 50-90% of the species on Earth. With massive amounts of dust lifted from desert regions and injected into the atmosphere, however, these dust storms often affect daily activities in dramatic ways: pushing grit through windows and doors, forcing people to stay indoors, causing breathing problems, reducing visibility and delaying flights, and by and large creating chaos. Thus, both increasing and decreasing concentrations of doses result in harmful biological effects; so do the airborne dust particles to our Living Earth. Since 1997 NASA has been successfully launching a series of satellites - the Earth Observing System - to intensively study, and gain a better understanding of, the Earth as an integrated system. Through participation in many satellite remote-sensing/retrieval and validation projects over the years, we have gradually developed and refined the SMART (Surface-sensing Measurements for Atmospheric Radiative Transfer) and COMMIT (Chemical, Optical & Microphysical Measurements of In-situ Troposphere) mobile observatories, a suite of surface remote sensing and in-situ instruments that proved to be vital in providing high temporal measurements, which complement the satellite observations. In this talk, we will present SMART-COMMIT which has played key roles, serving as network or supersite, in major international research projects such as the Joint Aerosol Monsoon Experiment (JAMEX), a core element of the Asian Monsoon Years (AMY, 2008-2012). SMART-COMMIT deployments during 2008 AMY/JAMEX were conducted in northwestern China to characterize the properties of dust-laden aerosols. In 2009, SMART-COMMIT also participated in the JAMEX/RAJO-MEGHA (Radiation, Aerosol Joint Observations-Monsoon Experiment in the Gangetic-Himalayan Area; Sanskrit for Dust-Cloud) to study the aerosol properties, solar absorption and the associated atmospheric warming, and the climatic impact of elevated aerosols during the premonsoon season in South Asia. To fully characterize the properties of airborne dust in the field is an important but challenging task. In this seminar, we will present our recent measurements and retrievals of airborne dust properties.

Tsay, Si-Chee↗
NASA NTRS2005

Computation of the Streamfunction and Velocity Potential for Limited and Irregular Domains

An algorithm is proposed for the computation of streamfunction and velocity potential from given horizontal velocity vectors based on solving a minimization problem. To guarantee the uniqueness of the solution and computational reliability of the algorithm, a Tikhonov regularization is applied. The solution implies that the obtained streamfunction and velocity potential have minimal magnitude, while the given velocity vectors can be accurately reconstructed from the computed streamfunction and velocity potential. Because the formulation of the minimization problem allows for circumventing the explicit specification of separate boundary conditions on the streamfunction and velocity potential, the algorithm is easily applicable to irregular domains. By using an advanced minimization algorithm with the use of adjoint techniques, the method is computationally efficient and suitable for problems with large dimensions. An example is presented for coastal oceans to illustrate the practical application of the algorithm.

velocity↗
DOE OSTI2020

Chaos on the hypercube

We analyze the spectral properties of a d-dimensional HyperCubic (HC) lattice model originally introduced by Parisi. The U(1) gauge links of this model give rise to a magnetic flux of constant magnitude φ but random orientation through the faces of the hypercube. The HC model, which also can be written as a model of 2d interacting Majorana fermions, has a spectral flow that is reminiscent of Maldacena-Qi (MQ) model, and its spectrum at φ = 0, actually coincides with the coupling term of the MQ model. As was already shown by Parisi, at leading order in 1/d, the spectral density of this model is given by the density function of the Q-Hermite polynomials, which is also the spectral density of the double-scaled Sachdev-Ye-Kitaev model. Parisi demonstrated this by mapping the moments of the HC model to Q-weighted sums on chord diagrams. We point out that the subleading moments of the HC model can also be mapped to weighted sums on chord diagrams, in a manner that descends from the leading moments. The HC model has a magnetic inversion symmetry that depends on both the magnitude and the orientation of the magnetic flux through the faces of the hypercube. The spectrum for fixed quantum number of this symmetry exhibits a transition from regular spectra at φ = 0 to chaotic spectra with spectral statistics given by the Gaussian Unitary Ensembles (GUE) for larger values of φ. For small magnetic flux, the ground state is gapped and is close to a Thermofield Double (TFD) state.

1/N expansion↗