Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Constellation: The autonomous control and data acquisition system for dynamic experimental setups

The operation of instruments and detectors in laboratory or beamline environments presents a complex challenge, requiring stable operation of multiple concurrent devices, often controlled by separate hardware and software solutions. These environments frequently undergo modifications, such as the inclusion of different auxiliary devices depending on the experiment or facility, adding further complexity. The successful management of such dynamic configurations demands a flexible and robust system capable of controlling data acquisition, monitoring experimental setups, enabling seamless reconfiguration, and integrating new devices with limited effort. This paper presents Constellation, a flexible and network-distributed control and data acquisition software framework tailored to laboratory and beamline environments, that addresses the limitations of existing solutions. The framework is designed with a focus on extensibility, providing a streamlined interface for instrument integration. It supports efficient system setup via network discovery mechanisms, promotes stability through autonomous operational features, and provides comprehensive documentation and supporting tools for operators and application developers such as controllers and logging interfaces. At the core of the architectural design is the autonomy of the individual components, called satellites, which can make independent decisions about their operation and communicate these decisions to other components. This paper introduces the design principles and framework architecture of Constellation, presents the available graphical user interfaces, shares insights from initial successful deployments, and provides an outlook on future developments and applications.

Autonomy↗
DOE OSTI2021

Power grid frequency prediction using spatiotemporal modeling

Understanding power system dynamics is essential for interarea oscillation analysis and the detection of grid instabilities. The FNET/GridEye is a GPS-synchronized wide-area frequency measurement network that provides an accurate picture of the normal real-time operational condition of the power system-dynamics, giving rise to new and intricate spatiotemporal patterns of power loads. We propose to model FNET/GridEye grid frequency data from the U.S. Eastern Interconnection with a spatiotemporal statistical model. We predict the frequency data at locations without observations, a critical need during disruption events where measurement data are inaccessible. Spatial information is accounted for either as neighboring measurements in the form of covariates or with a spatiotemporal correlation model captured by a latent Gaussian field. Finally, the proposed method is useful in estimating power system dynamic response from limited phasor measurements and holds promise for predicting instability that may lead to undesirable effects such as cascading outages.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

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↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2025

Direct statistical simulation of the Lorenz96 system in model reduction approaches

Direct statistical simulation (DSS) of nonlinear dynamical systems bypasses the traditional route of accumulating statistics by lengthy direct numerical simulations by solving the equations that govern the statistics themselves. DSS suffers, however, from the curse of dimensionality as the statistics (such as correlations) generally have higher dimensions than the underlying dynamical variables. Here we investigate two approaches to reduce the dimensionality of DSS, illustrating each method with numerical experiments with the Lorenz96 dynamical system. The forms of DSS chosen here involve approximate closures at second and third order in the equal-time cumulants. We demonstrate significant reduction in computational effort that can be achieved without sacrificing the accuracy of DSS. The methods developed here can be applied to turbulent fluid and magnetohydrodynamical systems. Published by the American Physical Society 2025

Li, Kuan↗
DOE OSTI2023

Constraining chaos: Enforcing dynamical invariants in the training of reservoir computers

Drawing on ergodic theory, we introduce a novel training method for machine learning based forecasting methods for chaotic dynamical systems. The training enforces dynamical invariants—such as the Lyapunov exponent spectrum and the fractal dimension—in the systems of interest, enabling longer and more stable forecasts when operating with limited data. The technique is demonstrated in detail using reservoir computing, a specific kind of recurrent neural network. Finally, results are given for the Lorenz 1996 chaotic dynamical system and a spectral quasi-geostrophic model of the atmosphere, both typical test cases for numerical weather prediction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Automating Discovery of Physics-Informed Neural State Space Models via Learning and Evolution

Recent works exploring deep learning application to dynamical systems modeling have demonstrated that embedding physical priors into neural networks can yield more effective, physically-realistic, and data-efficient models. However, in the absence of complete prior knowledge of a dynamical system's physical characteristics, determining the optimal structure and optimization strategy for these models can be difficult. In this work, we explore methods for discovering neural state space dynamics models for system identification. Starting with a design space of block-oriented state space models and structured linear maps with strong physical priors, we encode these components into a model genome alongside network structure, penalty constraints, and optimization hyperparameters. Demonstrating the overall utility of the design space, we employ an asynchronous genetic search algorithm that alternates between model selection and optimization and obtains accurate physically consistent models of three physical systems: an aerodynamics body, a continuous stirred tank reactor, and a two tank interacting system.

genetic algorithms, neural architecture search, ne↗
DOE OSTI2024

Limitations for Quantum Algorithms to Solve Turbulent and Chaotic Systems

We investigate the limitations of quantum computers for solving nonlinear dynamical systems. In particular, we tighten the worst-case bounds of the quantum Carleman linearisation (QCL) algorithm answering one of their open questions. We provide a further significant limitation for any quantum algorithm that aims to output a quantum state that approximates the normalized solution vector. Given a natural choice of coordinates for a dynamical system with one or more positive Lyapunov exponents and solutions that grow sub-exponentially, we prove that any such algorithm has complexity scaling at least exponentially in the integration time. As such, an efficient quantum algorithm for simulating chaotic systems or regimes is likely not possible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation

RTI International, in collaboration with Pacific Northwest National Laboratory (PNNL), Carbon Capture Simulation for Industry Impact (CCSI 2 ), Electricity Power Research Institute (EPRI), and West Virginia University (WVU), successfully completed a joint research effort in developing a cost-effective, resilient, load-following advanced CO 2 capture technology for natural gas power plants. The project’s objective was to develop a CO 2 capture process that maximizes the net present value (NPV) of the electricity sale by minimizing the levelized cost of electricity (LCOE) under dynamic plant loads and high renewable penetration environments. The two key innovations developed in this project were the use of (i) advanced water-lean solvents (WLSs) and (ii) process intensification equipment, such as a rotating packed bed (RPB) absorber and dual-stage flash regeneration. The process’s low CO 2 capture cost is realized through WLSs’ low energy required for solvent regeneration, which lowers the operating cost while RPBs intensify the absorption process and reduces the power plant capital cost. A suite of advanced computational and simulation packages was implemented to guide the process design, validate the dynamic response of the capture plant, evaluate system-wide performance, and maximize the power plant’s profit. The project also engaged with power producers and other stakeholders to ensure its technical relevance and techno-economic viability. The development of this highly disruptive CO 2 capture technology could accelerate the industry adoption and thereby lower the greenhouse gas emissions of the U.S. power sector. Deployment of this technology can increase the reliability and decrease the cost of electricity generation in the U.S. by enabling the use of low-carbon fossil fuels to balance fluctuations of renewable energy availability.

03 NATURAL GAS↗
DOE OSTI2020

Early-Time Exponential Instabilities in Nonchaotic Quantum Systems

The majority of classical dynamical systems are chaotic and exhibit the butterfly effect: a minute change in initial conditions has exponentially large effects later on. But this phenomenon is difficult to reconcile with quantum mechanics. One of the main goals in the field of quantum chaos is to establish a correspondence between the dynamics of classical chaotic systems and their quantum counterparts. In isolated systems in the absence of decoherence, there is such a correspondence in dynamics, but it usually persists only over a short time window, after which quantum interference washes out classical chaos. In this work, we demonstrate that quantum mechanics can also play the opposite role and generate exponential instabilities in classically nonchaotic systems within this early-time window. Our calculations employ the out-of-time-ordered correlator (OTOC)—a diagnostic that reduces to the Lyapunov exponent in the classical limit but is well defined for general quantum systems. We show that certain classically nonchaotic models, such as polygonal billiards, demonstrate a Lyapunov-like exponential growth of the OTOC at early times with Planck’s-constant-dependent rates. This behavior is sharply contrasted with the slow early-time growth of the analog of the OTOC in the systems’ classical counterparts. These results suggest that classical-to-quantum correspondence in dynamics is violated in the OTOC even before quantum interference develops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗
DOE OSTI2022

Variational multiscale reinforcement learning for discovering reduced order closure models of nonlinear spatiotemporal transport systems

Abstract A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Driven superconducting vortex dynamics in systems with twofold anisotropy in the presence of pinning

Abstract We examine the dynamics of superconducting vortices with twofold anisotropic interaction potentials driven over random pinning, and compare the behavior under drives applied along the hard and the soft anisotropy directions. As the driving force increases, the number of topological defects reaches a maximum near the depinning threshold, and then decreases as the vortices form one-dimensional (1D) chains. This coincides with a transition from a pinned nematic to a moving smectic aligned with the soft anisotropy direction. The system is generally more ordered when the drive is applied along the soft direction of the anisotropy. For driving along the hard direction, there is a critical value of the twofold anisotropy above which the system remains aligned with the soft direction. Hysteretic behavior appears upon cycling the driving force, with 1D vortex chains persisting during the decreasing leg below the threshold for chain formation for increasing drive. More anisotropic systems have a greater amount of structural disorder in the moving state. For lower anisotropy, the system forms a moving smectic-A state, while at higher anisotropy, a moving nematic state appears instead.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhancing T c in a composite superconductor/metal bilayer system: A dynamical cluster approximation study

It has been proposed that the superconducting transition temperature T c of an unconventional superconductor with a large pairing scale but strong phase fluctuations can be enhanced by coupling it to a metal. However, the general efficacy of this approach across different parameter regimes remains an open question. Using the dynamical cluster approximation, we study this question in a system composed of an attractive Hubbard layer in the intermediate coupling regime, where the magnitude of the attractive Coulomb interaction |U| is slightly larger than the bandwidth W, hybridized with a noninteracting metallic layer. We find that while the superconducting transition becomes more mean-field-like with increasing interlayer hopping, the superconducting transition temperature T c exhibits a nonmonotonic dependence on the strength of the hybridization t ⊥ . Furthermore, this behavior arises from a reduction of the effective pairing interaction in the correlated layer that outcompetes the growth in the intrinsic pair-field susceptibility induced by the coupling to the metallic layer. We find that the largest T c inferred here for the composite system is comparable to the maximum value currently estimated for the isolated negative-U Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Safety Risk Reliability Model Library

SR2ML is a software package which contains a set of safety and reliability models designed to be interfaced with the INL developed RAVEN code. These models can be employed to perform both static and dynamic system risk analysis and determine risk importance of specific elements of the considered system. Two classes of reliability models have been developed; the first class includes all classical reliability models (Fault-Trees, Event-Trees, Markov models and Reliability Block Diagrams) which have been extended to deal not only with Boolean logic values but also time dependent values. The second class includes several components aging models. Models of these two classes are designed to be included in a RAVEN ensemble model to perform time dependent system reliability analysis (dynamic analysis). Similarly, these models can be interfaced with system analysis codes to determine failure time of systems and evaluate accident progression (static analysis).

Wang, Congjian↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Studying impacts of communication system performance on dynamic stability of networked microgrid

The development of smart grid technologies has resulted in increased interdependence between power and communication systems. Many of the operations in the existing power system rely on a stable and secured communication system. For electrically weak systems and time‐critical applications, this reliance can be even greater, where a small degradation in communication performance can degrade system stability. However, despite inter‐dependencies between power and communication systems, only a few studies have investigated the impacts of communication system performance on power system dynamics. This study investigates the dependencies of power system dynamics operations on a communication system performance. First, a detailed, dynamic networked microgrid model is developed in the GridLAB‐D simulation environment, along with a representative multi‐traffic, multi‐channel, multi‐protocol communication system model, developed in the network simulator (ns‐3). Second, a hierarchical engine for large‐scale infrastructure co‐simulation framework is developed to co‐simulate microgrid dynamics, its communication system, and a microgrid control system. The impact of communication system delays on the dynamic stability of networked microgrids is evaluated for the loss of generation using three use‐cases. While the example use‐cases examine microgrid applications and the impact to resiliency, the framework can be applied to all levels of power system operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Data-driven model construction for industrial asset decision boundary classification

In some embodiments, a system model construction platform may receive, from a system node data store, system node data associated with an industrial asset. The system model construction platform may automatically construct a data-driven, dynamic system model for the industrial asset based on the received system node data. A synthetic attack platform may then inject at least one synthetic attack into the data-driven, dynamic system model to create, for each of a plurality of monitoring nodes, a series of synthetic attack monitoring node values over time that represent simulated attacked operation of the industrial asset. The synthetic attack platform may store, in a synthetic attack space data source, the series of synthetic attack monitoring node values over time that represent simulated attacked operation of the industrial asset. This information may then be used, for example, along with normal operational data to construct a threat detection model for the industrial asset.

97 MATHEMATICS AND COMPUTING↗