A program for interactive computation in linear systems theory Final report
FORTRAN 4 program for analysis and design of dynamic systems
SEARCH · Engineering Papers
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.
FORTRAN 4 program for analysis and design of dynamic systems
Self-reorganization methods for complex linear dynamic systems to compensate for component failures
High speed automated identification and design of dynamic systems, both linear and nonlinear, are discussed. Special emphasis is placed on developing hardware and techniques which are applicable to practical problems. The basic modeling experiment and new results are described. Using the improvements developed successful identification of several systems, including a physical example as well as simulated systems, was obtained. The advantages of parameter signature analysis over signal signature analysis in go-no go testing of operational systems were demonstrated. The feasibility of using these ideas in failure mode prediction in operating systems was also investigated. An improved digital controlled nonlinear function generator was developed, de-bugged, and completely documented.
The problem of estimation of state in nonlinear dynamical systems containing time delays is studied. The plant is specified by a set of nonlinear differential-difference equations. Observations are a nonlinear function of current and/or delayed states. Both contain additive disturbances. The criterion used for the optimal estimates is the integral of the weighted squared error. Using the theory of the calculus of variations, equations are developed for the estimation. They are first expressed in the form of a split boundary value problem, which is then converted to an initial value problem for on-line estimation. The result yields a sequential estimation scheme in which filtered and smoothed estimates are computed in a sequential manner. The applicability of the procedure is demonstrated by a practical example.
Recent advances in nonlinear dynamics demonstrate a remarkable complexity of patterns outside of equilibrium, which are derived from simple basic laws of physics.
A mathematical model is proposed to examine the vibration transmission through rolling element bearings in geared rotor systems. Current bearing models, based on either ideal boundary conditions for the shaft or purely translational stiffness element description, cannot explain how the vibratory motion may be transmitted from the rotating shaft to the casing. This study clarifies this issue qualitatively and quantitatively by developing a comprehensive bearing stiffness matrix of dimension 6 model for the precision rolling element bearings from basic principles. The proposed bearing formulation is extended to analyze the overall geared rotor system dynamics including casing and mounts. The bearing stiffness matrix is included in discrete system models using lumped parameter and/or dynamic finite element techniques. Eigensolution and forced harmonic response due to rotating mass unbalance or kinematic transmission error excitation for a number of examples are computed.
A new mathematical model is proposed to examine the vibration transmission through rolling element bearings in geared rotor systems. Current bearing models, based on either ideal boundary conditions for the shaft or purely translational stiffness element description, cannot explain how the vibratory motion may be transmitted from the rotating shaft to the casing. Experimental results have shown that the casing plate motion is primarily flexural. Here, this issue is clarified qualitatively and quantitatively by developing a comprehensive bearing stiffness matrix of dimension 6 to model precision rolling element bearings using basic principles. The proposed bearing stiffness matrix is partially verified using available analytical and experimental data, and is completely characterized. The study extends the proposed bearing formulation to analyze the overall geared rotor system dynamics including casing and mounts. The bearing stiffness matrix is included in discrete system models using lumped parameter and/or dynamic finite element techniques.
An accurate system-bath model to investigate the quantum dynamics of hydrogen atoms chemisorbed on graphene is presented. The system comprises a hydrogen atom and the carbon atom from graphene that forms the covalent bond, and it is described by a previously developed 4D potential energy surface based on density functional theory ab initio data. The bath describes the rest of the carbon lattice and is obtained from an empirical force field through inversion of a classical equilibrium correlation function describing the hydrogen motion. By construction, model building easily accommodates improvements coming from the use of higher level electronic structure theory for the system. Further, it is well suited to a determination of the system-environment coupling by means of ab initio molecular dynamics. This paper details the system-bath modeling and shows its application to the quantum dynamics of vibrational relaxation of a chemisorbed hydrogen atom, which is here investigated at T = 0 K with the help of the multi-configuration time-dependent Hartree method. Paper II deals with the sticking dynamics.
Active large angle slewing maneuvers of a multi-body flexible dynamic system are investigated. An appropriate state variable transformation and a feedback linearization technique are employed to transform the dynamics of the nonlinear system to a new state that is more amenable to control design procedures. Closed-loop feedback algorithms are implemented to perform slewing maneuvers, while simultaneously suppressing flexural vibrations of the system. Stability of this class of nonlinear systems is also investigated, whereby a sufficient condition for asymptotic stability of the system is established. Numerical examples are presented to demonstrate the proposed active control algorithms.
Advances in machine learning and deep neural networks has enabled complex engineering tasks like image recognition, anomaly detection, regression, and multi-objective optimization, to name but a few. The complexity of the algorithm architecture, e.g., the number of hidden layers in a deep neural network, typically grows with the complexity of the problems they are required to solve, leaving little room for interpreting (or explaining) the path that results in a specific solution. This drawback is particularly relevant for autonomous aerospace and aviation systems, where certifications require a complete understanding of the algorithm behavior in all possible scenarios. Including physics knowledge in such data-driven tools may improve the interpretability of the algorithms, thus enhancing model validation against events with low probability but relevant for system certification. Such events include, for example, spacecraft or aircraft sub-system failures, for which data may not be available in the training phase. This paper investigates a recent physics-informed learning algorithm for identification of system dynamics, and shows how the governing equations of a system can be extracted from data using sparse regression. The learned relationships can be utilized as a surrogate model which, unlike typical data-driven surrogate models, relies on the learned underlying dynamics of the system rather than large number of fitting parameters. The work shows that the algorithm can reconstruct the differential equations underlying the observed dynamics using a single trajectory when no uncertainty is involved. However, the training set size must increase when dealing with stochastic systems, e.g., nonlinear dynamics with random initial conditions.
This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .
The research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of the NASA Microgravity Science and Applications Program is described. The primary focus was on investigating fundamental questions concerning the propagation and extinction of premixed flames in Earth gravity and in microgravity environments. The approach was to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. A combination of one-dimensional and two-dimensional simulations was used to investigate the effects of curvature and dilution on ignition and propagation of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the mechanism leading to cellular instability, and to study the effects of gravity on the transition to cellular structure. A flame in a microgravity environment can be extinguished without external losses, and the mechanism leading to cellular structure is not preferential diffusion but a thermo-diffusive instability. The simulations have also lead to a better understanding of the interactions between buoyancy forces and the processes leading to thermo-diffusive instability.
Conventional Event-Tree (ET) based methodologies are extensively used as tools to perform reliability and safety assessment of complex and critical engineering systems. One of the disadvantages of these methods is that timing/sequencing of events and system dynamics are not explicitly accounted for in the analysis. In order to overcome these limitations several techniques, also know as Dynamic Probabilistic Risk Assessment (DPRA), have been developed. Monte-Carlo (MC) and Dynamic Event Tree (DET) are two of the most widely used D-PRA methodologies to perform safety assessment of Nuclear Power Plants (NPP). In the past two years, the Idaho National Laboratory (INL) has developed its own tool to perform Dynamic PRA: RAVEN (Reactor Analysis and Virtual control ENvironment). RAVEN has been designed to perform two main tasks: 1) control logic driver for the new Thermo-Hydraulic code RELAP-7 and 2) post-processing tool. In the first task, RAVEN acts as a deterministic controller in which the set of control logic laws (user defined) monitors the RELAP-7 simulation and controls the activation of specific systems. Moreover, the control logic infrastructure is used to model stochastic events, such as components failures, and perform uncertainty propagation. Such stochastic modeling is deployed using both MC and DET algorithms. In the second task, RAVEN processes the large amount of data generated by RELAP-7 using data-mining based algorithms. This report focuses on the analysis of dynamic stochastic systems using the newly developed RAVEN DET capability. As an example, a DPRA analysis, using DET, of a simplified pressurized water reactor for a Station Black-Out (SBO) scenario is presented.
Dissipative entanglement generation protocols embrace open quantum system dynamics to generate long-lived entangled states. In this letter, we report that zero-delay anti-bunching in the $g^{(2)}(\tau)$ correlation function is a signature of dissipative entanglement between actively driven quantum emitters coupled to a plasmonic resirvoir. The anti-bunching effect is a universal signature for entangled bi-partite states generated by relaxation to sub-radiant subspaces. We discuss how the anti-bunching can be validated on experimentally feasible timescales by a Hanburry Brown and Twiss interferometer and analyze the entanglement generation robustness with respect to variations in local detunings, system couplings, and driving fields.
While a new type of power take-off (PTO) based on mechanical-motion-rectifier has been developed, its dynamic model is over simplified and cannot properly predict its nonlinear system characteristics. There is a need to develop a high fidelity PTO model to predict the system dynamics. An equivalent circuit would be a convenient tool to model the system behaviour and can easily apply to a circuit-based wave-to-wire (W2W) model. In this paper, a methodology for deriving a circuit-based high-fidelity model through the PTO design is introduced following by the extraction of the parameters through experimental results. Finally, the PTO circuit model is validated by the bench test results under the time domain
Chaotic transport is a subject of paramount importance in a variety of problems in plasma physics, specially those related to anomalous transport and turbulence. On the other hand, a great deal of information on chaotic transport can be obtained from simple dynamical systems like two-dimensional area-preserving (symplectic) maps, where powerful mathematical results like KAM theory are available. In this work we review recent works on transport barriers in area-preserving maps, focusing on systems which do not obey the so-called twist property. For such systems KAM theory no longer holds everywhere and novel dynamical features show up as nonresistive reconnection, shearless curves and shearless bifurcations. After presenting some general features using a standard nontwist mapping, we consider magnetic field line maps for magnetically confined plasmas in tokamaks.
Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.
This paper presents a novel technique for generating a convex system model from an experimentally derived data set which features variance among repeated measurements. The convex system model developed as a test case characterizes the dynamic system losses of a vanadium redox flow battery as a function of the active power output and the battery state of charge. The technique hinges on a pre-cleaning via clustering procedure which precedes the formation of a planar convex hull comprised of triangular simplices. The clustering procedure efficiently reduces the experimental data set while mitigating variance among repeated measurements and removes outliers. Ultimately, the lower evelope of the planar convex hull serves as the desired convex system model. The proposed technique reduces systematic model error which is otherwise present when directly developing a planar convex hull model based on an unreduced data set.