Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex 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 199 records · Page 11

Saddle point approaches in strong field physics and generation of attosecond pulses

Attoscience is the emerging field that accesses the fastest electronic processes occurring at the atomic and molecular length scales with attosecond (1 as = 10 –18 s) time resolution having wide ranging physical, chemical, material science and biological applications. The quintessential and one of the most fundamental processes in this domain is the generation of phase locked XUV attosecond pulses. The theoretical approach to understand the process incorporates a fully quantum or semi classical or relativistic description of coherent charge dynamics in intense ultrashort electromagnetic fields driving a quantum system (an atom, a molecule, solid band gap materials or surface plasmas). Modelling of such physical and dynamical systems in science and also in many other branches often leads to equations represented in terms of complex multi-dimensional integrals. These integrals can often be solved using the stationary phase approximation, which leads to a series of equations identifying the points in the multi-dimensional space, having most significant contributions in their evaluation. These points are usually indicated as saddle points. The description of the dynamics of quantum mechanical or relativistic systems that results from such an approach enables near to classical physics intuitive perceptions of the processes under investigation. Thus, the saddle point methods are very powerful and valuable general theoretical tools to obtain asymptotic expressions of such solutions and help also to gain physical insights on the underlying phenomena. Such techniques developed in the past have been adapted to study the emission of as pulses by different physical systems and have been widely employed in calculating and estimating the response of matter to intense electromagnetic pulses on ultrafast time scales. Here we provide an extensive disposition of the saddle point approaches unifying their ubiquitous applications within the domain of attoscience valid for simple atomic to more complex condensed matter systems undergoing ultrafast dynamics and present current trends and advancements in the field. In this review we would delineate the methodology, present a synthesis of seminal works and describe the state of the art applications. Finally we also address ultrashort time dynamics of novel materials that have gained much attention recently, namely lower dimensional material systems and micro-plasma systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Battery Energy Storage Scenario Analyses Using the Lithium-Ion Battery Resource Assessment (LIBRA) Model

Meeting aggressive carbon emission goals will entail widespread deployment of renewable sources of electricity. Because these sources are variable, there is a need to develop scalable energy storage technologies. The U.S. Department of Energy is supporting efforts to increase U.S. manufacturing and recycling capabilities for LIBs and to decrease costs of stationary storage batteries. Many factors influence the domestic manufacturing and cost of stationary storage batteries, including availability of critical raw materials (lithium, cobalt, and nickel), competition from various demand sectors (consumer electronics, vehicles, and battery energy storage), resource recovery (recycling), government policies, and learning in the industry, among other factors. Understanding how these factors interact and identifying synergies and bottlenecks is important for developing effective strategies for the LIB stationary energy storage system. We developed the Lithium-Ion Battery Resource Analysis (LIBRA) model as a tool to help stakeholders better understand the following types of questions: What are the roles of R&D, industrial learning, and scaling of demand in lowering the cost of battery energy storage system production? How do the intersections between the EV and stationary storage sectors affect the battery supply chain? For various stationary storage and EV penetration scenarios, what volumes of critical materials might be required and what role can resource recovery play? What does expected demand for both EVs and stationary storage portend for mineral resources and overall mineral scarcity? The LIBRA model is developed using a System dynamics (SD) modeling approach to represent interactions across the segments of the battery materials supply chain. System dynamics models can capture the complex interactions and feedback between the various system components that influence supply and demand. The LIBRA model is comprised of several interacting modules that represent specific portions of the LIB supply chain. The model tracks the buildout of the domestic LIB industry over time (2020 - 2050) and in the context of competing demands for raw materials, recycling, and markets for LIBs. The LIBRA model represents major systemic feedback loops and delays across the supply chain. This report provides a complete documentation for the LIBRA model, including model assumptions, data, scenario analysis results, and sensitivity analysis of the model's input space.

25 ENERGY STORAGE↗

An anomaly detection framework for dynamic systems using a Bayesian hierarchical framework

Complex systems are susceptible to many types of anomalies, faults, and abnormal behavior caused by a variety of off-nominal conditions that may ultimately result in major failures or catastrophic events. Early and accurate detection of these anomalies using system inputs and outputs collected from sensors and smart devices has become a challenging problem and an active area of research in many application domains. In this article, we present a new Bayesian hierarchical framework that is able to model the relationship between system inputs (sensor measurements) and outputs (response variables) without imposing strong distributional/parametric assumptions while using only a subset of training samples and sensor attributes. Then, an optimal cost-sensitive anomaly detection framework is proposed to determine whether a sample is an anomalous one taking into consideration the trade-off between misclassification errors and detection rates. The model can be used for both supervised and unsupervised settings depending on the availability of data regarding the behavior of the system under anomaly conditions. The unsupervised model is particularly useful when it is prohibitive to identify in advance the anomalies that a system may present and where no data are available regarding the behavior of the system under anomaly conditions. A Bayesian hierarchical setting is used to structure the proposed framework and help with accommodating uncertainty, imposing interpretability, and controlling the sparsity and complexity of the proposed anomaly detection framework. A Markov chain Monte Carlo algorithm is also developed for model training using past data. Furthermore, the numerical experiments conducted using a simulated data set and a wind turbine data set demonstrate the successful application of the proposed work for system response modeling and anomaly detection.

17 WIND ENERGY↗

Adaptive Space-Time Methods for Large Scale Optimal Design

When modeling complex physical systems with advanced dynamics, such as shocks and singularities, many classic methods for solving partial differential equations can return inaccurate or unusable results. One way to resolve these complex dynamics is through r-adaptive refinement methods, in which a fixed number of mesh points are shifted to areas of high interest. The mesh refinement map can be found through the solution of the Monge-Ampére equation, a highly nonlinear partial differential equation. Due to its nonlinearity, the numerical solution of the Monge-Ampére equation is nontrivial and has previously required computationally expensive methods. In this report, we detail our novel optimization-based, multigrid-enabled solver for a low-order finite element approximation of the Monge-Ampére equation. This fast and scalable solver makes r-adaptive meshing more readily available for problems related to large-scale optimal design. Beyond mesh adaptivity, our report discusses additional applications where our fast solver for the Monge-Ampére equation could be easily applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Embracing fine-root system complexity in terrestrial ecosystem modeling

Projecting the dynamics and functioning of the biosphere requires a holistic consideration of whole-ecosystem processes. However, biases toward leaf, canopy, and soil modeling since the 1970s have constantly left fine-root systems being rudimentarily treated. As accelerated empirical advances in the last two decades establish clearly functional differentiation conferred by the hierarchical structure of fine-root orders and associations with mycorrhizal fungi, a need emerges to embrace this complexity to bridge the data-model gap in still extremely uncertain models. In this work, we propose a three-pool structure comprising transport and absorptive fine roots with mycorrhizal fungi (TAM) to model vertically resolved fine-root systems across organizational and spatial–temporal scales. Emerging from a conceptual shift away from arbitrary homogenization, TAM builds upon theoretical and empirical foundations as an effective and efficient approximation that balances realism and simplicity. A proof-of-concept demonstration of TAM in a big-leaf model both conservatively and radically shows robust impacts of differentiation within fine-root systems on simulating carbon cycling in temperate forests. Theoretical and quantitative support warrants exploiting its rich potentials across ecosystems and models to confront uncertainties and challenges for a predictive understanding of the biosphere. Echoing a broad trend of embracing ecological complexity in integrative ecosystem modeling, TAM may offer a consistent framework where modelers and empiricists can work together toward this grand goal.

59 BASIC BIOLOGICAL SCIENCES↗

Role of Interactions and Correlations on Collective Dynamics of Molecular Motors Along Parallel Filaments

Cytoskeletal motors known as motor proteins are molecules that drive cellular transport along several parallel cytoskeletal filaments and support many biological processes. Experimental evidences suggest that they interact with the nearest molecules of their filament while performing any mechanical work. These interactions modify the microscopic level properties of motor proteins. In this work, a new version of two-channel totally asymmetric simple exclusion process, that incorporates the intra-channel interactions in a thermodynamically consistent way, is proposed. As the existing approaches for multi-channel systems deviate from analyzing the combined effect of inter and intra-channel interactions, a new approach known as modified vertical cluster mean field is developed. The approach along with Monte Carlo simulations successfully encounters some correlations and computes the complex dynamic properties of the system. Role of symmetry of interactions and inter-channel coupling is observed on the phase diagrams, maximal particle current and its corresponding optimal interaction strength. Surprisingly, for all values of coupling rate and most of the interaction splittings, the optimal interaction strength corresponding to maximal current belongs to the case of weak repulsive interactions. Moreover, for weak interaction splittings and with an increase in the coupling rate, the optimal interaction strength tends towards the known experimental results. The effect of coupling as well as interaction energy is also measured for correlations. They are found to be short-range and weaker for repulsive and weak attractive interactions while they are long-range and stronger for large attractions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Off-policy deep reinforcement learning with automatic entropy adjustment for adaptive online grid emergency control

Electric overloading conditions and contingencies put modern power systems at risk of voltage collapse and blackouts. Load shedding is crucial to maintain voltage stability for grid emergency control. However, the rule- or model-based schemes rely on accurate dynamic system models and face considerable challenges in adapting to various operating conditions and uncertain event occurrences. Here, to address these issues, this paper proposes a novel deep reinforcement learning (DRL)-based voltage stability control algorithm with automatic entropy adjustment (AEA) for grid emergency control. Various dynamic network components for complex system operations are modeled to construct the DRL environment. An off-policy soft actor-critic architecture is developed to maximize the expected reward and policy entropy simultaneously. The AEA mechanism is proposed to facilitate the policy maximum entropy procedure, and the proposed method can automatically provide effective discrete and continuous actions against various fault scenarios. Our approach accomplishes high sampling efficiency, scalability, and auto-adaptivity of the control policies under high uncertainties. Comparative studies with the existing DRL-based control methods in IEEE benchmarks indicate salient performance improvement of the proposed method for dynamic system emergency control.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Chemical Imaging and Dynamical Studies of Reactivity and Emergent Behavior in Complex Interfacial Systems. Final Technical Report

This research program explored the efficacy of using molecular-level manipulation, imaging and scanning tunneling spectroscopy in conjunction with supersonic molecular beam gas-surface scattering to significantly enhance our understanding of chemical processes occurring on well-characterized interfaces. One program focus was on the spatially-resolved emergent behavior of complex reaction systems as a function of the local geometry and density of adsorbate-substrate systems under reaction conditions. Another focus was on elucidating the emergent electronic and related reactivity characteristics of intentionally constructed single and multicomponent atom- and nanoparticle-based materials. We also examined emergent chirality and self-organization in adsorbed molecular systems where collective interactions between adsorbates and the supporting interface lead to spatial symmetry breaking. In many of these studies we combined the advantages of scanning tunneling (STM) and atomic force (AFM) imaging, scanning tunneling local electronic spectroscopy (STS), and reactive supersonic molecular beams to elucidate precise details of interfacial reactivity that had not been observed by more traditional surface science methods. Using these methods, it was possible to examine, for example, the differential reactivity of molecules adsorbed at different bonding sites in conjunction with how reactivity is modified by the local configuration of nearby adsorbates. At the core of this effort was the goal of significantly extending our understanding of interfacial atomic-scale interactions to create, with intent, molecular assemblies and materials with advanced chemical and physical properties. This ambitious program addressed several key topics in DOE Grand Challenge Science, including emergent chemical and physical properties in condensed phase systems, novel uses of chemical imaging, and the development of advanced reactivity concepts in combustion and catalysis including carbon management. These activities directly benefitted national science objectives in the areas of chemical energy production and advanced materials development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing complexity and dynamics in epidemics: geographical barriers and facilitators of foot-and-mouth disease dissemination

Introduction: Physical and non-physical processes that occur in nature may influence biological processes, such as dissemination of infectious diseases. However, such processes may be hard to detect when they are complex systems. Because complexity is a dynamic and non-linear interaction among numerous elements and structural levels in which specific effects are not necessarily linked to any one specific element, cause-effect connections are rarely or poorly observed. Methods: To test this hypothesis, the complex and dynamic properties of geo-biological data were explored with high-resolution epidemiological data collected in the 2001 Uruguayan foot-and-mouth disease (FMD) epizootic that mainly affected cattle. County-level data on cases, farm density, road density, river density, and the ratio of road (or river) length/county perimeter were analyzed with an open-ended procedure that identified geographical clustering in the first 11 epidemic weeks. Two questions were asked: (i) do geo-referenced epidemiologic data display complex properties? and (ii) can such properties facilitate or prevent disease dissemination? Results: Emergent patterns were detected when complex data structures were analyzed, which were not observed when variables were assessed individually. Complex properties–including data circularity–were demonstrated. The emergent patterns helped identify 11 counties as ‘disseminators’ or ‘facilitators’ (F) and 264 counties as ‘barriers’ (B) of epidemic spread. In the early epidemic phase, F and B counties differed in terms of road density and FMD case density. Focusing on non-biological, geographical data, a second analysis indicated that complex relationships may identify B-like counties even before epidemics occur. Discussion: Geographical barriers and/or promoters of disease dispersal may precede the introduction of emerging pathogens. If corroborated, the analysis of geo-referenced complexity may support anticipatory epidemiological policies.

60 APPLIED LIFE SCIENCES↗

Tracking and predicting barges on inland waterways

We present a non-linear, probabilistic prediction model developed and implemented to track spatial location and other navigation characteristics of a barge traveling on the inland waterway system. A pre-filter, to check the validity of the measurements, a non-linear speed estimation process, and a Kalman filter to predict the navigation solution of the barge is developed in this work. Due to the complex dynamics involved in the system, a non-linear stochastic model was developed in state space using system dynamics to represent the process and measurement systems while maintaining the fidelity of an actual system. The algorithm was verified using actual measurements obtained from multiple barges on multiple rivers acquired from different sensors. The results show a reliable and robust prediction algorithm for tracking inland waterway barges.

Randeniya, Duminda I↗

Dynamic granularity of imaging systems

Imaging systems that include a specific source, imaging concept, geometry, and detector have unique properties such as signal-to-noise ratio, dynamic range, spatial resolution, distortions, and contrast. Some of these properties are inherently connected, particularly dynamic range and spatial resolution. It must be emphasized that spatial resolution is not a single number but must be seen in the context of dynamic range and consequently is better described by a function or distribution. We introduce the “dynamic granularity” G dyn as a standardized, objective relation between a detector’s spatial resolution (granularity) and dynamic range for complex imaging systems in a given environment rather than the widely found characterization of detectors such as cameras or films by themselves. We found that this relation can partly be explained through consideration of the signal’s photon statistics, background noise, and detector sensitivity, but a comprehensive description including some unpredictable data such as dust, damages, or an unknown spectral distribution will ultimately have to be based on measurements. Measured dynamic granularities can be objectively used to assess the limits of an imaging system’s performance including all contributing noise sources and to qualify the influence of alternative components within an imaging system. Our article explains the construction criteria to formulate a dynamic granularity and compares measured dynamic granularities for different detectors used in the X-ray backlighting scheme employed at Sandia’s Z-Backlighter facility.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High Performance-Multidisciplinary Simulation for Regional Scale Earthquake Hazard and Risk Assessments

Moderate-to-large magnitude earthquakes represent a significant risk around the world and a major issue across the Department of Energy (DOE) mission interests, ranging from DOE's one-of-a-kind mission critical facilities to national energy systems (nuclear power plants, electric/gas distribution systems, and renewable energy production facilities). As a result, for the DOE there is a strong motivation in promoting and developing a risk-informed and performance-based approach for both the assessment of existing facilities and the design of new structures. DOE National Labs have historically made multidisciplinary contributions to the solution of challenging seismic hazard and risk problems and have been at the forefront for the development of design standards for major facilities. The current approach for estimating the ground motion hazard for critical facilities is based upon a Probabilistic Seismic Hazard Analysis (PSHA), whereby measured ground motion data are collected and homogenized into a probabilistic framework to estimate the hazard for future events. However, recent studies have highlighted the site-specific nature of earthquakes and shown the limitations of a purely empirical approach. Moreover, the design of infrastructure systems has been moving toward a performance-based approach, in which specific performance goals are specified at different levels of excitation and for different limit states. Such an approach requires the availability of robust and reliable computational tools for the nonlinear dynamic analysis of complex structural systems. Ongoing research at Lawrence Berkeley National Laboratory (LBNL) and Lawrence Livermore National Laboratory (LLNL) aims at developing an advanced simulation environment at the exascale for a coupled assessment of earthquake hazard (ground motions) and earthquake risk (structural demands)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Faster network disruption from layered oscillatory dynamics

Nonlinear complex network-coupled systems typically have multiple stable equilibrium states. Following perturbations or due to ambient noise, the system is pushed away from its initial equilibrium, and, depending on the direction and the amplitude of the excursion, it might undergo a transition to another equilibrium. It was recently demonstrated [M. Tyloo, J. Phys. Complex. 3 03LT01 (2022)] that layered complex networks may exhibit amplified fluctuations. Here, I investigate how noise with system-specific correlations impacts the first escape time of nonlinearly coupled oscillators. Interestingly, I show that, not only the strong amplification of the fluctuations is a threat to the good functioning of the network but also the spatial and temporal correlations of the noise along the lowest-lying eigenmodes of the Laplacian matrix. Finally, I analyze first escape times on synthetic networks and compare noise originating from layered dynamics to uncorrelated noise.

97 MATHEMATICS AND COMPUTING↗

On the modeling and nonlinear dynamics of autonomous Silva-Young type chaotic oscillators with flat power spectrum

The nonlinear dynamics of fourth-order Silva-Young type chaotic oscillators with flat power spectrum recently introduced by Tamaseviciute and collaborators is considered. In this type of oscillators, a pair of semiconductor diodes in an anti-parallel connection acts as the nonlinear component necessary for generating chaotic oscillations. Based on the Shockley diode equation and an appropriate selection of the state variables, a smooth mathematical model (involving hyperbolic sine and cosine functions) is derived for a better description of both the regular and chaotic dynamics of the system. The complex behavior of the oscillator is characterized in terms of its parameters by using time series, bifurcation diagrams, Lyapunov exponents' plots, Poincaré sections, and frequency spectra. It is shown that the onset of chaos is achieved via the classical period-doubling and symmetry restoring crisis scenarios. Some PSPICE simulations of the nonlinear dynamics of the oscillator are presented in order to confirm the ability of the proposed mathematical model to accurately describe/predict both the regular and chaotic behaviors of the oscillator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗

High Fidelity Simulations of a Methane-Oxygen Rotating Detonation Rocket Engine

The rotating detonation combustor (RDC) is an important realization of pressure gain combustion for practical applications. The RDC system is characterized by a highly unsteady flow field, with multiple reflected pressure waves following the strong detonations and the existence of partially-burnt gases entrained within the detonation annulus during the injector recovery process. While experimental efforts have provided macroscopic properties of an RDC with a given fuel-oxidizer mixture and geometry, limited optical access and reduced lifespan of flow-field diagnostic equipment constrain understanding of the mechanisms that lend to wave stability, controllability, and sustainability. To this end, high-fidelity numerical simulations of the Air Force Research Laboratory (AFRL) methane-oxygen rotating detonation rocket engine (RDRE) are performed to provide detailed insight into the detonation physics and anomalous behavior within the combustor. The results highlight the merit of high-fidelity numerical studies in analyzing the combustion processes and injector dynamics of a complex combustion system. The simulations capture the general wave behavior observed within the experimental studies. Three different operating conditions, with varying global equivalence ratio and system mass flow rate, are considered: (1) baseline configuration, (2) high equivalence ratio, and (3) high mass flow rate. The high mass flow rate case results in increased detonation wave strength whereas the high equivalence ratio results in increased deflagrative heat release. The high equivalence ratio case results in chaotic behavior with two pairs of counter-rotating detonation waves, whereas a single pair of co-rotating waves is observed in the baseline and high mass flow rate cases. Further, parasitic combustion of the reactant mixture ahead of the detonation wave diminishes its strength and reduces propagation speed to approximately 50-60\% of the Chapman-Jouguet (CJ) condition. In this work, a detailed analysis of the detonation wave behavior, injector dynamics, and the statistics of combustion parameters are presented.

Prakash, Supraj↗

Experiments and Analyses of Data Transfers Over Wide-Area Dedicated Connections

Dedicated wide-area network connections are increasingly employed in high-performance computing and big data scenarios. One might expect the performance and dynamics of data transfers over such connections to be easy to analyze due to the lack of competing traffic. However, non-linear transport dynamics and end-system complexities (e.g., multi-core hosts and distributed filesystems) can in fact make analysis surprisingly challenging. We present extensive measurements of memory-to-memory and disk-to-disk file transfers over 10 Gbps physical and emulated connections with 0–366 ms round trip times (RTTs). For memory-to-memory transfers, profiles of both TCP and UDT throughput as a function of RTT show concave and convex regions; large buffer sizes and more parallel flows lead to wider concave regions, which are highly desirable. TCP and UDT both also display complex throughput dynamics, as indicated by their Poincare maps and Lyapunov exponents. For disk-to-disk transfers, we determine that high throughput can be achieved via a combination of parallel I/O threads, parallel network threads, and direct I/O mode. Our measurements also show that Lustre filesystems can be mounted over long-haul connections using LNet routers, although challenges remain in jointly optimizing file I/O and transport method parameters to achieve peak throughput.

Rao, Nageswara S.↗