Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system 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 1,117 records · Page 62

Numerical methods for the analysis of sampled-data systems and for the computation of system zeros

MARSYAS is a computer aided control system design package for the simulation and analysis of open loop and closed loop dynamic systems. Outlined here is the numerical and theoretical basis of the MARSYAS functions developed during the summer of 1991. In particular, the numerical computation of the matrix exponential e(exp A) = I + A + A(exp 2)/2! + A(exp 3)/3! + ... and the numerical computation of the finite system zeros of a dynamic system (continuous time or discrete time) are presented.

Hodel, A. Scottedward↗

Resonances in a Chaotic Attractor Crisis of the Lorenz Flow

Local bifurcations of stationary points and limit cycles have successfully been characterized in terms of the critical exponents of these solutions. Lyapunov exponents and their associated covariant Lyapunov vectors have been proposed as tools for supporting the understanding of critical transitions in chaotic dynamical systems. However, it is in general not clear how the statistical properties of dynamical systems change across a boundary crisis during which a chaotic attractor collides with a saddle. This behavior is investigated here for a boundary crisis in the Lorenz flow, for which neither the Lyapunov exponents nor the covariant Lyapunov vectors provide a criterion for the crisis. Instead, the convergence of the time evolution of probability densities to the invariant measure, governed by the semigroup of transfer operators, is expected to slow down at the approach of the crisis. Such convergence is described by the eigenvalues of the generator of this semigroup, which can be divided into two families, referred to as the stable and unstable Ruelle–Pollicott resonances, respectively. The former describes the convergence of densities to the attractor (or escape from a repeller) and is estimated from many short time series sampling the state space. The latter is responsible for the decay of correlations, or mixing, and can be estimated from a long times series, invoking ergodicity. It is found numerically for the Lorenz flow that the stable resonances do approach the imaginary axis during the crisis, as is indicative of the loss of global stability of the attractor. On the other hand, the unstable resonances, and a fortiori the decay of correlations, do not flag the proximity of the crisis, thus questioning the usual design of early warning indicators of boundary crises of chaotic attractors and the applicability of response theory close to such crises.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Krylov vector methods for model reduction and control of flexible structures

Krylov vectors and the concept of parameter matching are combined here to develop model-reduction algorithms for structural dynamics systems. The method is derived for a structural dynamics system described by a second-order matrix differential equation. The reduced models are shown to have a promising application in the control of flexible structures. It can eliminate control and observation spillovers while requiring only the dynamic spillover terms to be considered. A model-order reduction example and a flexible structure control example are provided to show the efficacy of the method.

Su, Tzu-Jeng↗

On Information Transfer-Based Characterization of Power System Stability

In this study, we present a novel approach to identify the generators and states responsible for the small-signal stability of power networks. To this end, the newly developed notion of information transfer between the states of a dynamical system is used. In particular, using the concept of information transfer, which characterizes influence between the various states of a dynamical system, we identify the generators and states which are responsible for causing instability of the power network. While characterizing influence from state to state, information transfer can also describe influence from state to modes thereby generalizing the well-known notion of participation factor while at the same time overcoming some of the limitations of the participation factor. The developed framework is applied to reproduce known results for the three bus system, identifying the various causes of instabilities, and is extended to IEEE 39 bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Validation of the IMCE RFDA System for Dynamic Elastic Moduli Measurements

The focus of this report is to describe the new system that has been acquired for the Low Activation Material Development and Analysis Laboratory (LAMDA) to perform the measurement of the dynamic Young’s and shear moduli, and to show that the results of the new system are more easily analyzed and more accurate than the current systems. The current system is the GrindoSonic MK5 Impulse Excitation and Vibration Test System from J.W. Lemmens, Inc., and the new system is the Resonant Frequency and Damping Analyser (RFDA) from Integrated Material Control Engineering (IMCE). A comparison between the results from both systems was performed with nine specimen shapes and five graphite grades. All of the results showed agreement of the moduli values within 1-2%. The only outliers were a 20% disagreement for one specimen shape and a 5-10% disagreement for a second specimen shape, but it was shown that this occurred because the GrindoSonic system was reporting the frequency for the flexural harmonic peak rather than the fundamental torsional peak.

36 MATERIALS SCIENCE↗

Converting DYNAMO simulations to Powersim Studio simulations

DYNAMO is a computer program for building and running 'continuous' simulation models. It was developed by the Industrial Dynamics Group at the Massachusetts Institute of Technology for simulating dynamic feedback models of business, economic, and social systems. The history of the system dynamics method since 1957 includes many classic models built in DYANMO. It was not until the late 1980s that software was built to take advantage of the rise of personal computers and graphical user interfaces that DYNAMO was supplanted. There is much learning and insight to be gained from examining the DYANMO models and their accompanying research papers. We believe that it is a worthwhile exercise to convert DYNAMO models to more recent software packages. We have made an attempt to make it easier to turn these models into a more current system dynamics software language, Powersim © Studio produced by Powersim AS 2 of Bergen, Norway. This guide shows how to convert DYNAMO syntax into Studio syntax.

97 MATHEMATICS AND COMPUTING↗

Statistical Properties for Flows with Unbounded Roof Function, Including the Lorenz Attractor

For geometric Lorenz attractors (including the classical Lorenz attractor) we obtain a greatly simplified proof of the central limit theorem which applies also to the more general class of codimension two singular hyperbolic attractors. We also obtain the functional central limit theorem and moment estimates, as well as iterated versions of these results. A consequence is deterministic homogenisation (convergence to a stochastic differential equation) for fast-slow dynamical systems whenever the fast dynamics is singularly hyperbolic of codimension two.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Detecting Faults By Use Of Hidden Markov Models

Frequency of false alarms reduced. Faults in complicated dynamic system (e.g., antenna-aiming system, telecommunication network, or human heart) detected automatically by method of automated, continuous monitoring. Obtains time-series data by sampling multiple sensor outputs at discrete intervals of t and processes data via algorithm determining whether system in normal or faulty state. Algorithm implements, among other things, hidden first-order temporal Markov model of states of system. Mathematical model of dynamics of system not needed. Present method is "prior" method mentioned in "Improved Hidden-Markov-Model Method of Detecting Faults" (NPO-18982).

Smyth, Padhraic J.↗

Attitude Design for the LADEE Mission

The Lunar Atmosphere and Dust Environment Explorer (LADEE) satellite successfully completed its 148-day science investigation in a low-altitude, near-equatorial lunar orbit on April 18, 2014. The LADEE spacecraft was built, managed and operated by NASA's Ames Research Center (ARC). The Mission Operations Center (MOC) was located at Ames and was responsible for activity planning, command sequencing, trajectory and attitude design, orbit determination, and spacecraft operations. The Science Operations Center (SOC) was located at Goddard Space Flight Center and was responsible for science planning, data archiving and distribution. This paper details attitude design and operations support for the LADEE mission. LADEE's attitude design was shaped by a wide range of instrument pointing requirements that necessitated regular excursions from the baseline one revolution per orbit "Ram" attitude. Such attitude excursions were constrained by a number of flight rules levied to protect instruments from the Sun, avoid geometries that would result in simultaneous occlusion of LADEE's two star tracker heads, and maintain the spacecraft within its thermal and power operating limits. To satisfy LADEE's many attitude requirements and constraints, a set of rules and conventions was adopted to manage the complexity of this design challenge and facilitate the automation of ground software that generated pointing commands spanning multiple days of operations at a time. The resulting LADEE Flight Dynamics System (FDS) that was developed used Visual Basic scripts that generated instructions to AGI's Satellite Tool Kit (STK) in order to derive quaternion commands at regular intervals that satisfied LADEE's pointing requirements. These scripts relied heavily on the powerful "align and constrain" capability of STK's attitude module to construct LADEE's attitude profiles and the slews to get there. A description of the scripts and the attitude modeling they embodied is provided. One particular challenge analysts faced was in the design of LADEE maneuver attitudes. A flight rule requiring pre-maneuver verification of in-flight maneuver conditions by ground operators prior to burn execution resulted in the need to accommodate long periods in the maneuver attitude. This in turn complicated efforts to satisfy star tracker interference and communication constraints in lunar orbit. In response to this challenge, a graphical method was developed and used to survey candidate rotation angles about the thrust vector. This survey method is described and an example of its use on a particular LADEE maneuver is discussed. Finally, the software and methodology used to satisfy LADEE's attitude requirements are also discussed in the context of LADEE's overall activity planning effort. In particular, the way in which strategic schedules of instrument and engineering activities were translated into actual attitude profiles at the tactical level, then converted into precise quaternion commands to achieve those pointing goals is explained. In order to reduce the risk of time-consuming re-planning efforts, this process included the generation of long-term projections of constraint violation predictions for individual attitude profiles that could be used to establish keep-out time-frames for particular attitude profiles. The challenges experienced and overall efficacy of both the overall LADEE ground system and the attitude components of the Flight Dynamics System in meeting LADEE's varied pointing requirements are discussed.

LADEE↗

Cooperative storage of shared files in a parallel computing system with dynamic block size

Improved techniques are provided for parallel writing of data to a shared object in a parallel computing system. A method is provided for storing data generated by a plurality of parallel processes to a shared object in a parallel computing system. The method is performed by at least one of the processes and comprises: dynamically determining a block size for storing the data; exchanging a determined amount of the data with at least one additional process to achieve a block of the data having the dynamically determined block size; and writing the block of the data having the dynamically determined block size to a file system. The determined block size comprises, e.g., a total amount of the data to be stored divided by the number of parallel processes. The file system comprises, for example, a log structured virtual parallel file system, such as a Parallel Log-Structured File System (PLFS).

97 MATHEMATICS AND COMPUTING↗

SU-D-BRF-07: Ultrasound and Fluoroscopy Based Intraoperative Image-Guidance System for Dynamic Dosimetry in Prostate Brachytherapy

Purpose: Prostate brachytherapy is a common treatment method for low-risk prostate cancer patients. Intraoperative treatment planning is known to improve the treatment procedure and the outcome. The current limitation of intraoperative treatment planning is the inability to localize the seeds in relation to the prostate. We developed an image-guidance system to fulfill this need to achieve intraoperative dynamic dosimetry in prostate brachytherapy. Methods: Our system is based on standard imaging equipments available in the operating room, including the transrectal ultrasound (TRUS) and the mobile C-arm. A simple fiducial is added to compute the C-arm pose. Three fluoroscopic images and an ultrasound volume of the seeds and the prostate are acquired and processed by four image processing algorithms: seed segmentation, fiducial detection with pose estimation, seed reconstruction, and seeds-to-TRUS registration. The updated seed positions allow the physician to assess the quality of implantation and dynamically adjust the treatment plan during the course of surgery to achieve improved exit dosimetry. Results: The system was tested on 10 phantoms and 37 patients. Seed segmentation resulted in a 1% false negative and 2% false positive rates. Fiducial detection with pose estimation resulted in a detection rate of 98%. Seed reconstruction had a mean reconstruction error of 0.4 mm. Seeds-to-TRUS registration had a mean registration error of 1.3 mm. The total processing time from image acquisition to registration was approximately 1 minute. Conclusion: We present an image-guidance system for intraoperative dynamic dosimetry in prostate brachytherapy. Using standard imaging equipments and a simple fiducial, our system can be easily adopted in any clinics. Robust image processing algorithms enable accurate and fast computation of the delivered dose. Especially, the system enables detection of possible hot/cold spots during the surgery, allowing the physician to address these in the operating room rather than requiring additional therapy afterwards. This research was supported in part by the National Institutes of Health/National Cancer Institute (NIH/NCI) under grant 2R44CA099374 and grant 1R01CA151395, and in part by the Department of Defense (DoD) under grant W81XWH-05-1-0407.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A deep learning enabler for nonintrusive reduced order modeling of fluid flows

Here in this paper, we introduce a modular deep neural network (DNN) framework for data-driven reduced order modeling of dynamical systems relevant to fluid flows. We propose various DNN architectures which numerically predict evolution of dynamical systems by learning from either using discrete state or slope information of the system. Our approach has been demonstrated using both residual formula and backward difference scheme formulas. However, it can be easily generalized into many different numerical schemes as well. We give a demonstration of our framework for three examples: (i) Kraichnan-Orszag system, an illustrative coupled nonlinear ordinary differential equation, (ii) Lorenz system exhibiting chaotic behavior, and (iii) a nonintrusive model order reduction framework for the two-dimensional Boussinesq equations with a differentially heated cavity flow setup at various Rayleigh numbers. Using only snapshots of state variables at discrete time instances, our data-driven approach can be considered truly nonintrusive since any prior information about the underlying governing equations is not required for generating the reduced order model. Our a posteriori analysis shows that the proposed data-driven approach is remarkably accurate and can be used as a robust predictive tool for nonintrusive model order reduction of complex fluid flows.

42 ENGINEERING↗

Dynamic systems-engineering process - The application of concurrent engineering

A system engineering methodology is described which enables users, particulary NASA and DOD, to accommodate changing needs; incorporate emerging technologies; identify, quantify, and manage system risks; manage evolving functional requirements; track the changing environment; and reduce system life-cycle costs. The approach is a concurrent, dynamic one which starts by constructing a performance model defining the required system functions and the interrelationships. A detailed probabilistic risk assessment of the system elements and their interrelationships is performed, and quantitative analysis of the reliability and maintainability of an engineering system allows its different technical and process failure modes to be identified and their probabilities to be computed. Decision makers can choose technical solutions that maximize an objective function and minimize the probability of failure under resource constraints.

Wiskerchen, Michael J.↗

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↗

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↗

From globally coupled maps to complex-systems biology

Studies of globally coupled maps, introduced as a network of chaotic dynamics, are briefly reviewed with an emphasis on novel concepts therein, which are universal in high-dimensional dynamical systems. They include clustering of synchronized oscillations, hierarchical clustering, chimera of synchronization and desynchronization, partition complexity, prevalence of Milnor attractors, chaotic itinerancy, and collective chaos. The degrees of freedom necessary for high dimensionality are proposed to equal the number in which the combinatorial exceeds the exponential. Future analysis of high-dimensional dynamical systems with regard to complex-systems biology is briefly discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗