Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Parameter Estimation for a Hybrid Adaptive Flight Controller

This paper expands on the hybrid control architecture developed at the NASA Ames Research Center by addressing issues related to indirect adaptation using the recursive least squares (RLS) algorithm. Specifically, the hybrid control architecture is an adaptive flight controller that features both direct and indirect adaptation techniques. This paper will focus almost exclusively on the modifications necessary to achieve quality indirect adaptive control. Additionally this paper will present results that, using a full non -linear aircraft model, demonstrate the effectiveness of the hybrid control architecture given drastic changes in an aircraft s dynamics. Throughout the development of this topic, a thorough discussion of the RLS algorithm as a system identification technique will be provided along with results from seven well-known modifications to the popular RLS algorithm.

Campbell, Stefan F.↗
NASA NTRS2003

Self-Supervised Dynamical Systems

Some progress has been made in a continuing effort to develop mathematical models of the behaviors of multi-agent systems known in biology, economics, and sociology (e.g., systems ranging from single or a few biomolecules to many interacting higher organisms). Living systems can be characterized by nonlinear evolution of probability distributions over different possible choices of the next steps in their motions. One of the main challenges in mathematical modeling of living systems is to distinguish between random walks of purely physical origin (for instance, Brownian motions) and those of biological origin. Following a line of reasoning from prior research, it has been assumed, in the present development, that a biological random walk can be represented by a nonlinear mathematical model that represents coupled mental and motor dynamics incorporating the psychological concept of reflection or self-image. The nonlinear dynamics impart the lifelike ability to behave in ways and to exhibit patterns that depart from thermodynamic equilibrium. Reflection or self-image has traditionally been recognized as a basic element of intelligence. The nonlinear mathematical models of the present development are denoted self-supervised dynamical systems. They include (1) equations of classical dynamics, including random components caused by uncertainties in initial conditions and by Langevin forces, coupled with (2) the corresponding Liouville or Fokker-Planck equations that describe the evolutions of probability densities that represent the uncertainties. The coupling is effected by fictitious information-based forces, denoted supervising forces, composed of probability densities and functionals thereof. The equations of classical mechanics represent motor dynamics that is, dynamics in the traditional sense, signifying Newton s equations of motion. The evolution of the probability densities represents mental dynamics or self-image. Then the interaction between the physical and metal aspects of a monad is implemented by feedback from mental to motor dynamics, as represented by the aforementioned fictitious forces. This feedback is what makes the evolution of probability densities nonlinear. The deviation from linear evolution can be characterized, in a sense, as an expression of free will. It has been demonstrated that probability densities can approach prescribed attractors while exhibiting such patterns as shock waves, solitons, and chaos in probability space. The concept of self-supervised dynamical systems has been considered for application to diverse phenomena, including information-based neural networks, cooperation, competition, deception, games, and control of chaos. In addition, a formal similarity between the mathematical structures of self-supervised dynamical systems and of quantum-mechanical systems has been investigated.

Zak, Michail↗
DOE OSTI2023

Bi-Level Linear Programming Model for Automatic Load Shedding: A Distributed Wide-Area Measurement System-based Solution

Load shedding is currently implemented as a two-step based approach. In the first step, manual load shedding is taken place, were system operators, using estimates, inform distribution utilities of predicted stressful conditions. Information provided include the potential use of energy reserves, as well as load shedding amount. In a second step, automatic load shedding is done. The latter is realized using protection relays. While considering frequency variation, pre-defined values of load to be shed and correspondent number of stages for such to be realized are transformed into relay settings. Under-frequency protection relays use only local measurements towards decision making, thus operate in a decentralized architecture. Decision making is done in milliseconds plus breaker time. While this approach has provided much system reliability, considering the new smart grid paradigm, where system dynamics are much faster due to increasing renewable resources penetration, in some operating conditions it will generate sub-optimal solutions, such as islanding. Phasor measurement units provide a source of information which can be useful for this problem. Centralized architecture-based solutions for automatic load shedding, as present in the state-of-the-art, require though total processing times which are not acceptable for real-life implementation. In this work, considering the above, a bi-level linear programming model is presented. The model is implemented considering a distributed architecture while leveraging phasor measurement units data. The upper-level model estimates the current system state. Results of this model are embedded in a lower-level model, which decision variables are the location and load value to be shed. Easy-to-implement model, built-on the classic weighted least squares solution, highlight potential aspects towards real-life applications.

Bretas, Arturo Suman↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

ExaWind: Predictive Wind Energy Simulations

This presentation describes the ExaWind project and the team's progress in creating a suite of performance-portable codes designed for predictive simulations of wind farms on next-generation exascale-class supercomputers. Such simulations will require the resolution of scales spanning many orders of magnitude, from blade boundary layers to wind farm flow structures. In the U.S., the first exascale systems will be GPU accelerated, and different GPU manufacturers have been chosen for the different systems. At the heart of the ExaWind software is a hybrid-solver approach based on the codes Nalu-Wind and AMR-Wind, which are computational fluid dynamics solvers for the incompressible Navier-Stokes equations. Nalu-Wind is an unstructured-grid code used to resolve wind turbine geometry and blade boundary layers, whereas AMR-Wind is a structured-grid background solver for atmospheric turbulent flow and turbine wake propagation. The models are coupled with overset meshes and global linear systems are approximated through a loose-coupling algorithm. Results will include validation-quality high-fidelity simulations and strong/weak scaling results from the Summit supercomputer.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
NASA NTRS1974

A movable-mass attitude stabilization system for cable-connected artificial-g space stations

The development of an active, momentum-exchange system to be used for attitude stabilization of a class of cable-connected artificial-g space stations is studied. A system which employs a single movable control mass is examined for the control of a space station which has the physical appearance of two cylinders connected axially by cables. The dynamic model for the space station includes its aggregate rigid body rotation and relative torsional rotation between the bodies. A zero torsional stiffness design (one cable) and a maximum torsional stiffness design (eight cables) are examined in various stages of deployment, for selected spin velocities ranging from 4 rpm upwards. A linear, time-invariant, feed-back control system is employed, with gains calculated via a root-specification procedure. The movable mass controller provides critical wobble-damping capability for the crew quarters for all configurations and spin velocity.

Childs, D. W.↗
NASA NTRS1987

Unsteady aerodynamics of blade rows

The requirements placed on an unsteady aerodynamic theory intended for turbomachinery aeroelastic applications are discussed along with a brief description of the various theoretical models that are available to address these requirements. The main emphasis is placed on the description of a linearized inviscid theory which fully accounts for the effects of a nonuniform mean or steady flow on unsteady aerodynamic response. Although this theory has been developed primarily for blade flutter prediction, more general equations are presented which account for unsteady excitations due to incident external aerodynamic disturbances as well as those due to blade motions. The resulting equations consist of a system of three field equations along with conditions imposed at blade, wake and shock surfaces and in the far field. These equations can be solved to determine the fluctuations in all fluid dynamic properties throughout the required solution domain. Example solutions are presented to demonstrate several effects associated with nonuniform steady flows on the linearized unsteady aerodynamic response to prescribed blade motions.

Verdon, J. M.↗
NASA NTRS1971

Statistical analysis of the modal properties of large structural systems.

A theory is developed to predict eigenvalue and eigenvector uncertainty in large dynamic models. The uncertainty is based on physical property uncertainty and should not be confused with numerical roundoff, although the method can be extended to include the latter. The theory, when implemented on a computer, is used to analyze the uncertainties in frequencies and mode shapes based on uncertainties in mass, stiffness, modulus of elasticity, etc. The method incorporates a linear statistical model which is quite adequate for handling property uncertainties of 10% or more. The model is not limited to small systems but uses certain statistical assumptions as well as selective matrix manipulations to keep the size of all matrix operations to within the number of degrees of freedom of the system. Examples are given for two longitudinal vibration problems, and the results are supported by a Monte Carlo analysis.

Collins, J. D.↗
DOE OSTI2022

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales t ∝ μ –1 , with μ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number N of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales τ ≈ µ –1 $\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales τ ≈ μ –1 log(N). Here, we present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1982

Sled Control and Safety System

Computerized system for controlling motion of linear-track accelerator applied to other automated equipment, such as numerically-controlled machine tools and robot manipulators on assembly lines. System controls motions of sled with sine-wave signal created digitally by microprocessor. Dynamic parameters of sled motion are monitored so sled may be stopped safely if malfunction occurs. Sled is capable of sinusoidal accelerations up to 0.5 g with 125-kg load.

Forrest, L. J.↗
NASA NTRS1981

Information and display requirements for aircraft terrain following

The display design procedure for manned vehicle systems, is applied and validated, for a particular scenario. The scenario chosen is that of zero visibility high speed terrain following (V = 466 ft/sec, H = 200 ft) with an A-10 aircraft. The longitudal (linearized) dynamics are considered. The variations in (command path over) terrain pi(t) are modeled as a third order random process. The display design methodology is based on the optimal control model of pilot response, and employs this model in various ways in different phases of the design process. The overall methodology indicates that the design process is intended as a precursor to manned simulation. It provides a rank ordering of candidate displays through a three level process.

Kleinman, D. L.↗
NASA NTRS1984

Differential rotation in a solar-driven quasi-axisymmetric circulation

The concept of a quasi-axisymmetric circulation is used to explore the global scale dynamics of planetary atmospheres. A numerical circulation model applicable to Jupiter is presented, and an analytical study is performed elucidating the conditions leading to differential rotation in an atmosphere which is convectively unstable. A linear system forced by solar differential heating is considered, with nonlinear effects arising from advection being represented in the form of eddy diffusion. An empirical, latitudinal spectrum of the observed zonal wind field on Jupiter is discussed. Numerical solutions are presented which reveal banded wind fields with alternating and equatorial zonal jets and a multicellular Ferrel-Thomson meridional circulation consistent with the observed cloud striations on Jupiter. The vertical derivatives are parameterized to construct a simplified one-layer model.

Mayr, H. G.↗
NASA NTRS1987

Viscous damped space structure for reduced jitter

A technique to provide modal vibration damping in high performance space structures was developed which uses less than one once of incompressible fluid. Up to 50 percent damping can be achieved which can reduce the settling times of the lowest structural mode by as much as 50 to 1. This concept allows the designers to reduce the weight of the structure while improving its dynamic performance. Damping by this technique is purely viscous and has been shown by test to be linear over 5 orders of input magnitude. Amplitudes as low as 0.2 microinch were demonstrated. Damping in the system is independent of stiffness and relatively insensitive to temperature.

Wilson, James F.↗
NASA NTRS1988

Longitudinal long-period dynamics of aerospace craft

Linear analyses are performed to examine the generic aspects of aerospace vehicle longitudinal long-period or trajectory modes of motion. The influence of Mach number, dynamic pressure, thrust-to-drag ratio, and propulsion system thrust laws on the longitudinal trajectory modes is presented in terms of phugoid frequency and damping and height mode stability. The results of these analyses are compared to flying qualities requirements where possible, and potential deficiencies in both the vehicle and the criteria are noted. A preliminary look at possible augmentation schemes to improve potential deficiencies is also presented. Interpretation of the practical consequences of the results is aided by typical time histories. Results indicate that propulsion system characteristics are the dominant influence on the longitudinal long-period flight dynamics of hypersonic aerospace craft. However, straightforward augmentation systems demonstrated the potential to accommodate these influences if the effects are included in the design process. These efforts may be hampered by a lack of design criteria for hypersonic aircraft.

Berry, Donald T.↗
NASA NTRS1952

Flight Investigation of a Mechanical Feel Device in an Irreversible Elevator Control System of a Large Airplane

Report presents the results of measurements of the longitudinal stability and control characteristics of a large airplane using a mechanical feel device in combination with a booster incorporated in the elevator-control system. Tests were made to investigate the feasibility of eliminating the aerodynamic control forces through use of a booster and of providing control-feel forces mechanically. The feel device consisted of a centering spring which restrained the control stick through a linkage which was changed as a function of the dynamic pressure. Provisions were made for trimming and for manual adjustment of the force gradient. The system was designed to approximate the control-force characteristics that would result with a conventional elevator control with linear hinge-moment characteristics.

Brown, B Porter↗
NASA NTRS2017

Understanding the Sun-Earth Libration Point Orbit Formation Flying Challenges For WFIRST and Starshade

In order to fly an occulter in formation with a telescope at the Sun-Earth L2 (SEL2) Libration Point, one must have a detailed understanding of the dy-namics that govern the restricted three body system. For initial purposes, a linear approximation is satisfactory, but operations will require a high-fidelity modeling tool along with strategic targeting methods in order to be successful. This paper focuses on the challenging dynamics of the transfer trajectories to achieve the relative positioning of two spacecraft to fly in formation at SEL2, in our case, the Wide-Field Infrared Survey Telescope (WFIRST) and a proposed Starshade. By modeling the formation transfers using a high fidelity tool, an accurate V approximation can be made to as-sist with the development of the subsystem design required for a WFIRST and Starshade formation flight mission.

Sun-Earth-Moon↗