Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chaotic dynamics”

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 343 records · Page 19

Nonlinear dynamics and quantum chaos of a family of kicked p -spin models

Herein we introduce kicked p-spin models describing a family of transverse Ising-like models for an ensemble of spin-1/2 particles with all-to-all p-body interaction terms occurring periodically in time as delta-kicks. This is the natural generalization of the well-studied quantum kicked top (p = 2) [Haake, Kus', and Scharf, Z. Phys. B 65, 381 (1987)]. We fully characterize the classical nonlinear dynamics of these models, including the transition to global Hamiltonian chaos. The classical analysis allows us to build a classification for this family of models, distinguishing between p = 2 and p > 2, and between models with odd and even p's. Quantum chaos in these models is characterized in both kinematic and dynamic signatures. For the latter, we show numerically that the growth rate of the out-of-time-order correlator is dictated by the classical Lyapunov exponent. Finally, we argue that the classification of these models constructed in the classical system applies to the quantum system as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics of two fluids under periodic acceleration

The evolution of the interface between two fluids confined in a rectangular cavity is investigated numerically to predict its transient behavior. These computations address mixing characteristics of fluids under microgravity conditions, particularly g-jitter conditions resulting from aircraft vibration or crew motion, and have applications in solution crystal growth. The two-dimensional formulation employs the Boussinesq approximation and treats the mixing of two fluids as an initial value problem with a prescribed concentration field. For fluid mixing of practical applications inside a cavity, it is shown that nonlinear convective transport can dominate over both viscous and molecular diffusion. However, viscous diffusion can become important for low Reynolds number or certain cavity sizes. A stacking phenomenon is shown to occur for aspect ratios (Ar) approximately between 5 to 10. For aspect ratios in the neighborhood of 0.1 to 0.2, inner cells evolve at the interface with increasing time. For a square cavity of Ar = 1, chaotic mixing of the fluid appears to occur for a Stokes-Reynolds number greater than about 5.

Duval, Walter M. B.↗

2.1.3.404 - WEC Array Power Management and Output Simulation Tool

An array of wave energy converter (WEC) devices has variations in power output due to the chaotic nature of the waves. Eliminating or mitigating the power fluctuations is important for reducing the integration impacts of WEC plants in both distribution and transmission grids, and in standalone isolated power systems. Reduced variability of WEC-generated power in combination with energy storage or power management control at each WEC and at the array level will help increasing hosting capacity of distribution feeders for this type of variable renewable generation, and minimization of electric losses. To reduce risks and risk perception, gain key stakeholder acceptance, and enable developers to design effective and compatible energy plants with arrays of WEC devices, it requires the use of modeling tools to simulate the resource environment, device dynamics, utility power system response, as well as the development of an interface for an array controller. The project will create a publicly accessible numerical modeling framework to empower the wave energy sector to design projects of various scales (kW-100s MW), which are optimized on a plant performance basis and are compatible with different power systems and wave conditions. The framework will integrate with WEC-Sim, a wave environment model (SWAN-FUNWAVE), as well as established relevant electrical analysis tools, to model the grid system and interconnection, and optimize power output and power management for the WEC array.

POWER TRANSMISSION AND DISTRIBUTION,TIDAL AND WAVE↗

Active Fluids Form System-Spanning Filamentary Networks

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

active nematics↗

Section on Observed Impacts on El Nino

Agricultural applications of El Nino forecasts are already underway in some countries and need to be evaluated or re-evaluated. For example, in Peru, El Nino forecasts have been incorporated into national planning for the agricultural sector, and areas planted with rice and cotton (cotton being the more drought-tolerant crop) are adjusted accordingly. How well are this and other such programs working? Such evaluations will contribute to the governmental and intergovernmental institutions, including the Inter-American Institute for Global Change Research and the US National Ocean and Atmospheric Agency that are fostering programs to aid the effective use of forecasts. As El Nino climate forecasting grows out of the research mode into operational mode, the research focus shifts to include the design of appropriate modes of utilization. Awareness of and sensitivity to the costs of prediction errors also grow. For example, one major forecasting model failed to predict the very large El Nino event of 1997, when Pacific sea-surface temperatures were the highest on record. Although simple correlations between El Nino events and crop yields may be suggestive, more sophisticated work is needed to understand the subtleties of the interplay among the global climate system, regional climate patterns, and local agricultural systems. Honesty about the limitations of an forecast is essential, especially when human livelihoods are at stake. An end-to-end analysis links tools and expertise from the full sequence of ENSO cause-and-effect processes. Representatives from many disciplines are needed to achieve insights, e.g, oceanographers and atmospheric scientists who predict El Nino events, climatologists who drive global climate models with sea-surface temperature predictions, agronomists who translate regional climate connections in to crop yield forecasts, and economists who analyze market adjustments to the vagaries of climate and determine the value of climate forecasts. Methods include historical studies to understand past patterns and to test hindcasts of the prediction tools, crop modeling, spatial analysis and remote sensing. This research involves expanding, deepening, and applying the understanding of physical climate to the fields of agronomy and social science; and the reciprocal understanding of crop growth and farm economics to climatology. Delivery of a regional climate forecast with no information about how the climate forecast was derived limits its effectiveness. Explanation of a region's major climate driving forces helps to place a seasonal forecast in context. Then, a useful approach is to show historical responses to previous El Nino events, and projections, with uncertainty intervals, of crop response from dynamic process crop growth models. Regional ID forecasts should be updated with real-time weather conditions. Since every El Nino event is different, it is important to track, report and advise on each new event as it unfolds. The stability of human enterprises depends on understanding both the potentialities and the limits of predictability. Farmers rely on past experience to anticipate and respond to fluctuations in the biophysical systems on which their livelihoods depend. Now scientists are improving their ability to predict some major elements of climate variability. The improvements in the reliability of El Nino forecasts are encouraging, but seasonal forecasts for agriculture are not, and will probably never be completely infallible, due to the chaotic nature of the climate system. Uncertainties proliferate as we extend beyond Pacific sea-surface temperatures to climate teleconnections and agricultural outcomes. The goal of this research is to shed as a clear light as possible on these inherent uncertainties and thus to contribute to the development of appropriate responses to El Nino and other seasonal forecasts for a range of stakeholders, which, ultimately, includes food consumers everywhere.

Rosenzweig, Cynthia↗

Effective Power-Law Dependence of Lyapunov Exponents on the Central Mass in Galaxies

Using both numerical and analytical approaches, we demonstrate the existence of an effective power-law relation L alpha m(sup p) between themean Lyapunov exponent L of stellar orbits chaotically scattered by a supermassive black hole (BH) in the centre of a galaxy and the mass parameter m, i.e. ratio of the mass of the BH over the mass of the galaxy. The exponent p is found numerically to obtain values in the range p approximately equals 0.3-0.5. We propose a theoretical interpretation of these exponents, based on estimates of local 'stretching numbers', i.e. local Lyapunov exponents at successive transits of the orbits through the BH's sphere of influence. We thus predict p = 2/3 − q with q approximately equaling 0.1-0.2. Our basic model refers to elliptical galaxy models with a central core. However, we find numerically that an effective power-law scaling of L with m holds also in models with central cusp, beyond a mass scale up to which chaos is dominated by the influence of the cusp itself. We finally show numerically that an analogous law exists also in disc galaxies with rotating bars. In the latter case, chaotic scattering by the BH affects mainly populations of thick tube-like orbits surrounding some low-order branches of the x(sub 1) family of periodic orbits, as well as its bifurcations at low-order resonances, mainly the inner Lindblad resonance and the 4/1 resonance. Implications of the correlations between L and m to determining the rate of secular evolution of galaxies are discussed.

chaos aEuro" galaxies:↗

Spectral anomalies and broken symmetries in maximally chaotic quantum maps

Spectral statistics such as the level spacing statistics and spectral form factor (SFF) are widely expected to accurately identify “ergodicity,” including the presence of underlying macroscopic symmetries, in generic quantum systems ranging from quantized chaotic maps to interacting many-body systems. By studying various quantizations of maximally chaotic maps that break a discrete classical symmetry upon quantization, we demonstrate that this approach can be misleading and fail to detect macroscopic symmetries. Notably, the same classical map can exhibit signatures of different random matrix symmetry classes in short-range spectral statistics depending on the quantization. While the long-range spectral statistics encoded in the early time ramp of the SFF are more robust and correctly identify macroscopic symmetries in several common quantizations, we also demonstrate analytically and numerically that the presence of Berry-like phases in the quantization leads to spectral anomalies, which break this correspondence. Finally, we provide numerical evidence that long-range spectral rigidity remains directly correlated with ergodicity in the quantum dynamical sense of visiting a complete orthonormal basis.

Shou, Laura [Univ. of Maryland, College Park, MD (↗

MHD-induced beta limits in the Large Helical Device

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Influence of Planetary Mass on the Dynamical Lifetime of Planetary Systems

Recent numerical and analytic studies of planetary orbits have demonstrated the importance of resonances and chaos in destabilizing planetary systems. Newton's "clockwork" description of regular, predictable planetary orbits has been replaced by a view in which many systems can have long but finite lifetimes. This new knowledge has altered our perceptions of the later stages of planetary growth and of the stability of planetary systems. Stability criteria are inexact and time dependent. Most previous studies have focused on the effects in initial planetary orbits on the stability of the system. We are conducting an investigation which focuses on the dependence of stability criteria on planetary mass. Synthetic systems are created by increasing the masses of the planets in our Solar System or of the moons of a particular planet; these systems are then integrated until orbit crossing occurs. We have found that over some ranges, the time until orbit crossing varies to a good approximation as a power clothe factor by which the masses of the secondaries arc increased; some scatter occurs as a consequence of vie chaotic nature of orbital evolution. The slope of this power law varies substantially from system to system, and for moons it is mildly dependent on the inclusion of the planet's quadrupole moment in the gravitational potential.

Lissauer, J. J.↗

Reversible to irreversible transitions for ac driven skyrmions on periodic substrates

Abstract Using atomistic simulations, we investigate the dynamical behavior of magnetic skyrmions in dimer and trimer molecular crystal arrangements, as well as bipartite lattices at 3/2 and 5/2 fillings, under ac driving over a square array of anisotropy defects. For low ac amplitudes, at all fillings reversible motion appears in which the skyrmions return to their original positions at the end of each ac drive cycle and the diffusion is zero. We also identify two distinct irreversible regimes. The first is a translating regime in which the skyrmions form channels of flow in opposing directions and translate by one substrate lattice constant per ac drive cycle. The translating state appears in the dimer and trimer arrangements, and produces pronounced peaks in the diffusivity in the direction perpendicular to the external drive. For larger ac amplitudes, we find chaotic irreversible motion in which the skyrmions can randomly exchange places with each other over time, producing long-time diffusive behavior both parallel and perpendicular to the ac driving direction.

36 MATERIALS SCIENCE↗

Multiple grid problems on concurrent-processing computers

Three computer codes were studied which make use of concurrent processing computer architectures in computational fluid dynamics (CFD). The three parallel codes were tested on a two processor multiple-instruction/multiple-data (MIMD) facility at NASA Ames Research Center, and are suggested for efficient parallel computations. The first code is a well-known program which makes use of the Beam and Warming, implicit, approximate factored algorithm. This study demonstrates the parallelism found in a well-known scheme and it achieved speedups exceeding 1.9 on the two processor MIMD test facility. The second code studied made use of an embedded grid scheme which is used to solve problems having complex geometries. The particular application for this study considered an airfoil/flap geometry in an incompressible flow. The scheme eliminates some of the inherent difficulties found in adapting approximate factorization techniques onto MIMD machines and allows the use of chaotic relaxation and asynchronous iteration techniques. The third code studied is an application of overset grids to a supersonic blunt body problem. The code addresses the difficulties encountered when using embedded grids on a compressible, and therefore nonlinear, problem. The complex numerical boundary system associated with overset grids is discussed and several boundary schemes are suggested. A boundary scheme based on the method of characteristics achieved the best results.

Eberhardt, D. S.↗

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗

Atmospheric Dynamics of the Outer Planets

The giant planets-Jupiter, Saturn, Uranus, and Neptune-are fluid objects. The winds are powered by absorbed sunlight, as on earth, and by internal heat left over from planetary formation. The main constituents of the atmospheres are hydrogen and helium. The clouds are made of ammonia, hydrogen sulphide, and water. All four giant planets are banded, with multiple zonal jet streams. Even Uranus, whose spin axis is tipped by 98deg relative to the orbit axis, shows latitudinal banding and zonal jets. Equator-to-pole temperature differences are close to zero. Wind speeds are larger than on earth and do not decrease with distance from the sun. Although the power/area at Neptune is only 1/20 that at Jupiter, the winds at Neptune are three times stronger. Stable vortices like the Great Red Spot of Jupiter and similar spots on Neptune come in all size ranges and exhibit a variety of behaviours including merging, orbiting, filament ejection, and oscillating in both shape and position. At least at cloud-top levels, 90% of the long-lived vortices are anticyclonic and sit in anticyclonic shear zones. Features in the cyclonic zones tend to be chaotic, with lifetimes of several days or less. These mesoscale eddies tend to have lightning in them, which suggests that they get their energy from moist convection. The rate of conversion of eddy kinetic energy into kinetic energy of the zonal jets is more than 10% of the power/area radiated by Jupiter. This fraction is more than an order of magnitude larger than on earth. Several lines of evidence now indicate that the winds at cloud-top levels are the surface manifestation of deep-rooted motions that extend into the interior and are presumably driven by internal heat.

Ingersoll, Andrew P.↗

Computing Sensitivities in Evolutionary Systems: A Real-time Reduced Order Modeling Strategy

We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.

Reduced order model↗

Tensor-based flow reconstruction from optimally located sensor measurements

Reconstructing high-resolution flow fields from sparse measurements is a major challenge in fluid dynamics. Existing methods often vectorize the flow by stacking different spatial directions on top of each other, hence confounding the information encoded in different dimensions. Here, we introduce a tensor-based sensor placement and flow reconstruction method which retains and exploits the inherent multidimensionality of the flow. We derive estimates for the flow reconstruction error, storage requirements and computational cost of our method. We show, with examples, that our tensor-based method is significantly more accurate than similar vectorized methods. Furthermore, the variance of the error is smaller when using our tensor-based method. While the computational cost of our method is comparable to similar vectorized methods, it reduces the storage cost by several orders of magnitude. The reduced storage cost becomes even more pronounced as the dimension of the flow increases. We demonstrate the efficacy of our method on three examples: a chaotic Kolmogorov flow, in situ and satellite measurements of the global sea surface temperature and three-dimensional unsteady simulated flow around a marine research vessel.

Mechanics↗

Universal Anomalous Diffusion of Quantized Vortices in Ultraquantum Turbulence

In classical viscous fluids, turbulent eddies are known to be responsible for the rapid spreading of embedded particles. However, in inviscid quantum fluids where the turbulence is induced by a chaotic tangle of quantized vortices, dispersion of the particles can be achieved via a nonclassical mechanism, i.e., their binding to the evolving vortices. However, knowledge on how the vortices diffuse and spread in quantum-fluid turbulence is very limited, especially for the so-called ultraquantum turbulence (UQT) generated by a random tangle of vortices. Here we report a systematic numerical study of the apparent diffusion of vortices in UQT in superfluid helium-4 using the full Biot-Savart simulation. We reveal that the vortices in the superfluid exhibit a universal anomalous diffusion (superdiffusion) at small times, which transits to normal diffusion at large times. This behavior is found to be the result of a generic scaling property of the vortex velocity. Our simulation at finite temperatures also nicely reproduces recent experimental observations. Lastly, the knowledge obtained from this study may form the base for understanding turbulent transport and universal vortex dynamics in various quantum fluids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Designing capture trajectories to unstable periodic orbits around Europa

The hostile environment of third body perturbations restricts a mission designer's ability to find well-behaved reproducible capture trajectories when dealing with limited control authority as is typical with low-thrust missions. The approach outlined in this paper confronts this shortcoming by utilizing dynamical systems theory and an extensive preexisting database of Restricted Three Body Problem (RTBP) periodic orbits. The stable manifolds of unstable periodic orbits are utilized to attract a spacecraft towards Europa. By selecting an appropriate periodic orbit, a mission designer can control important characteristics of the captured state including stability, minimum altitudes, characteristic inclinations, and characteristic radii among others. Several free parameters are optimized in the non-trivial mapping from the RTBP to a more realistic model. Although the ephemeris capture orbit is ballistic by design, low-thrust is used to target the state that leads to the capture orbit, control the spacecraft after arriving on the unstable quasi-periodic orbit, and begin the spiral down towards the science orbit. The approach allows a mission designer to directly target fuel efficient captures at Europa in an ephemeris model. Furthermore, it provides structure and controllability to the design of capture trajectories that reside in a chaotic environment.

low-thrust↗