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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2020

CP2K: An Electronic Structure and Molecular Dynamics Software Package - Quickstep: Efficient and Accurate Electronic Structure Calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension. TDK has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 716142). VRR has been supported by the Swiss National Science Foundation in the form of Ambizione grant No. PZ00P2 174227 and RZK by the Natural Sciences and Engineering Research Council of Canada (NSERC) through Discovery Grants (RGPIN-2016-0505). GKS and CJM are supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. UK based work was funded under the embedded CSE programme of the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk), grants eCSE03-011, eCSE06-6, eCSE08-9, eCSE13-17 and the EPSRC (EP/P022235/1) grant “Surface and Interface Toolkit for the Materials Chemistry Community". Computational resources were provided by the Swiss National Supercomputing Centre (CSCS) and Compute Canada. The generous allocation of computing time on the FPGA-based supercomputer “Noctua" at PC2 is kindly acknowledged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2012

Minimal Coordinate Formulation of Contact Dynamics in Operational Space

In recent years, complementarity techniques have been developed for modeling non-smooth contact and collision dynamics problems for multi-link robotic systems. Normally, in this approach, a linear complementarity problem (LCP) is set up using 6n non-minimal coordinates for a system with n links together with all the unilateral constraints and inter-link bilateral constraints on the system. In this paper, we use operational space dynamics to develop a complementarity formulation for contact and collision dynamics that uses minimal coordinates. The use of such non-redundant coordinates results in much smaller size LCP problems and the automatic enforcement of the inter-link bilateral constraints. Furthermore, we exploit operational space low-order computational algorithms to overcome some of the bottlenecks in using minimal coordinates.

contact dynamics↗

NASA NTRS1983

Development of monofilar rotor hub vibration absorber

A design and ground test program was conducted to study the performance of the monofilar absorber for vibration reduction on a four-bladed helicopter. A monofilar is a centrifugal tuned two degree-of-freedom rotor hub absorber that provides force attenuation at two frequencies using the same dynamic mass. Linear and non-linear analyses of the coupled monofilar/airframe system were developed to study tuning and attenuation characteristics. Based on the analysis, a design was fabricated and impact bench tests verified the calculated non-rotating natural frequencies and mode shapes. Performance characteristics were measured using a rotating absorber test facility. These tests showed significant attenuation of fixed-system 4P hub motions due to 3P inplane rotating-system hub forces. In addition, detuning effects of the 3P monofilar modal response were small due to the nonlinearities and tuning pin slippage. However, attenuation of 4P hub motions due to 5P inplane hub forces was poor. The performance of the 5P monofilar modal response was degraded by torsional motion of the dynamic mass relative to the support arm which resulted in binding of the dynamic components. Analytical design studies were performed to evaluate this torsional motion problem. An alternative design is proposed which may alleviate the torsional motion of the dynamic mass.

Duh, J.↗

NASA NTRS2004

State-Dependent Riccati Equation Regulation of Systems with State and Control Nonlinearities

The state-dependent Riccati equations (SDRE) is the basis of a technique for suboptimal feedback control of a nonlinear quadratic regulator (NQR) problem. It is an extension of the Riccati equation used for feedback control of linear problems, with the addition of nonlinearities in the state dynamics of the system resulting in a state-dependent gain matrix as the solution of the equation. In this paper several variations on the SDRE-based method will be considered for the feedback control problem with control nonlinearities. The control nonlinearities may result in complications in the numerical implementation of the control, which the different versions of the SDRE method must try to overcome. The control methods will be applied to three test problems and their resulting performance analyzed.

Beeler, Scott C.↗

NASA NTRS1992

Spacecraft nonlinear control

The feedback linearization technique is applied to the problem of spacecraft attitude control and momentum management with control moment gyros (CMGs). The feedback linearization consists of a coordinate transformation, which transforms the system to a companion form, and a nonlinear feedback control law to cancel the nonlinear dynamics resulting in a linear equivalent model. Pole placement techniques are then used to place the closed-loop poles. The coordinate transformation proposed here evolves from three output functions of relative degree four, three, and two, respectively. The nonlinear feedback control law is presented. Stability in a neighborhood of a controllable torque equilibrium attitude (TEA) is guaranteed and this fact is demonstrated by the simulation results. An investigation of the nonlinear control law shows that singularities exist in the state space outside the neighborhood of the controllable TEA. The nonlinear control law is simplified by a standard linearization technique and it is shown that the linearized nonlinear controller provides a natural way to select control gains for the multiple-input, multiple-output system. Simulation results using the linearized nonlinear controller show good performance relative to the nonlinear controller in the neighborhood of the TEA.

Sheen, Jyh-Jong↗

DOE OSTI2023

Quantum dynamics simulations of the 2D spectroscopy for exciton polaritons

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Nonlinear Architectures Can Alter the Dynamics of Polymer–Nanoparticle Composites

Polymer nanocomposites exhibit remarkable physical properties that are attractive for many applications. These systems have been so far investigated by using linear polymer chains; the role of polymer matrix architecture on local dynamics, bulk rheology and nanoparticle motion remains unexplored. Here, using quasielastic neutron scattering, bulk rheology and x-ray photon correlation spectroscopy, we investigated nanocomposites with spherical silica nanoparticles well dispersed in poly (ethylene oxide) matrices having different architectures (specifically linear, stars and hyperbranched). The results reveal that the mechanical reinforcement of the nanocomposites with the nonlinear polymers can be altered by orders of magnitude with respect to the conventional nanocomposite with the linear polymer. Additionally, the polymer compactness and inter-penetrability are found to play crucial roles in determining their bulk rheology. At microscopic level, average segmental dynamics is remarkably slowed down by the attractive NPs in the matrices of high degree of branching whereas no significant effect is observed in linear polymer matrix at the same NP loading. In addition, the nanoscale dynamics of particles in the compact nonlinear matrices exhibit strong decoupling from the bulk viscoelasticity, allowing their fast relaxation even at ≈ 30% by volume. These results provide an experimental evidence that macromolecular architecture is a powerful new tool for tuning the bulk rheological properties as well as the nanoscale dynamics of PNCs without the need for changing polymer molecular weight, nanoparticle size, shape, loading or dispersion state.

36 MATERIALS SCIENCE↗

NASA NTRS2011

Chance-Constrained Guidance With Non-Convex Constraints

Missions to small bodies, such as comets or asteroids, require autonomous guidance for descent to these small bodies. Such guidance is made challenging by uncertainty in the position and velocity of the spacecraft, as well as the uncertainty in the gravitational field around the small body. In addition, the requirement to avoid collision with the asteroid represents a non-convex constraint that means finding the optimal guidance trajectory, in general, is intractable. In this innovation, a new approach is proposed for chance-constrained optimal guidance with non-convex constraints. Chance-constrained guidance takes into account uncertainty so that the probability of collision is below a specified threshold. In this approach, a new bounding method has been developed to obtain a set of decomposed chance constraints that is a sufficient condition of the original chance constraint. The decomposition of the chance constraint enables its efficient evaluation, as well as the application of the branch and bound method. Branch and bound enables non-convex problems to be solved efficiently to global optimality. Considering the problem of finite-horizon robust optimal control of dynamic systems under Gaussian-distributed stochastic uncertainty, with state and control constraints, a discrete-time, continuous-state linear dynamics model is assumed. Gaussian-distributed stochastic uncertainty is a more natural model for exogenous disturbances such as wind gusts and turbulence than the previously studied set-bounded models. However, with stochastic uncertainty, it is often impossible to guarantee that state constraints are satisfied, because there is typically a non-zero probability of having a disturbance that is large enough to push the state out of the feasible region. An effective framework to address robustness with stochastic uncertainty is optimization with chance constraints. These require that the probability of violating the state constraints (i.e., the probability of failure) is below a user-specified bound known as the risk bound. An example problem is to drive a car to a destination as fast as possible while limiting the probability of an accident to 10(exp -7). This framework allows users to trade conservatism against performance by choosing the risk bound. The more risk the user accepts, the better performance they can expect.

FROM↗

NASA NTRS1999

X-33 Attitude Control Using the XRS-2200 Linear Aerospike Engine

The Vehicle Control Systems Team at Marshall Space Flight Center, Structures and Dynamics Laboratory, Guidance and Control Systems Division is designing, under a cooperative agreement with Lockheed Martin Skunkworks, the Ascent, Transition, and Entry flight attitude control systems for the X-33 experimental vehicle. Test flights, while suborbital, will achieve sufficient altitudes and Mach numbers to test Single Stage To Orbit, Reusable Launch Vehicle technologies. Ascent flight control phase, the focus of this paper, begins at liftoff and ends at linear aerospike main engine cutoff (MECO). The X-33 attitude control system design is confronted by a myriad of design challenges: a short design cycle, the X-33 incremental test philosophy, the concurrent design philosophy chosen for the X-33 program, and the fact that the attitude control system design is, as usual, closely linked to many other subsystems and must deal with constraints and requirements from these subsystems. Additionally, however, and of special interest, the use of the linear aerospike engine is a departure from the gimbaled engines traditionally used for thrust vector control (TVC) in launch vehicles and poses certain design challenges. This paper discusses the unique problem of designing the X-33 attitude control system with the linear aerospike engine, requirements development, modeling and analyses that verify the design.

Hall, Charles E.↗

NASA NTRS1978

On approximating hereditary dynamics by systems of ordinary differential equations

The paper deals with methods of obtaining approximate solutions to linear retarded functional differential equations (hereditary systems). The basic notion is to project the infinite dimensional space of initial functions for the hereditary system onto a finite dimensional subspace. Within this framework, two particular schemes are discussed. The first uses well-known piecewise constant approximations, while the second is a new method based on piecewise linear approximating functions. Numerical results are given.

Cliff, E. M.↗

NASA NTRS1983

Structural optimization with dynamic behavior constraints

The minimum weight optimum design of damped linearly elastic structural systems subjected to periodic loading with behavior constraints on maximum deflections and side constraints on design variables is addressed. Attention is focused on the two major impediments to an optimal solution: (1) the time parametric nature of the behavior constraints; and (2) the severe nonconvexity of the design space. A solution method based on upper bound approximations for the behavior constraints and an innovative mathematical programming scheme for seeking the optimal frequency subspace is set forth. Numerical results for several test problems illustrate the effectiveness of the method reported.

Mills-Curran, W. C.↗

NASA NTRS1984

Non-Gaussian approach for parametric random vibration of non-linear structures

The dynamic response of a nonlinear, single degree of freedom structural system subjected to a physically white noise parametric excitation is investigated. The Ito stochastic calculus is employed to derive a general differential equation for the moments of the response coordinates. The differential equations of moments of any order are found to be coupled with higher order moments. A non-Gaussian closure scheme is developed to truncate the moment equations up to fourth order. The statistical of the stationary response are computed numerically and compared with analytical solutions predicted by a Gaussian closure scheme and the stochastic averaging method. It is found that the computed results exhibit the jump phenomenon which is typical of the characteristics of deterministic nonlinear systems. In addition, the numerical algorithm leads to multiple solutions all of which give positive mean squares. However, two of these solutions are found to violate the properties of high order moments. One solution preserves the moments properties and demonstrates that the system achieves a stationary response.

Ibrahim, R. A.↗

NASA NTRS1994

Attitude stability of a spinning spacecraft during appendage deployment/retraction

The work presented is motivated by the need for a national satellite rescue policy, not the ad hoc policy now in place. In studying different approaches for a national policy, the issue of capture and stabilization of a tumbling spacecraft must be addressed. For a rescue mission involving a tumbling spacecraft, it may be advantageous to have a rescue vehicle which is compact and 'rigid' during the rendezvous/capture phase. After capture, passive stabilization techniques could be utilized as an efficient means of detumbling the resulting system (i.e., both the rescue vehicle and captures spacecraft). Since the rescue vehicle is initially compact and 'rigid,' significant passive stabilization through energy dissipation can only be achieved through the deployment of flexible appendages. Once stabilization is accomplished, retraction of the appendages before maneuvering the system to its final destination may also prove advantageous. It is therefore of paramount interest that we study the effect of appendage deployment/retraction on the attitude stability of a spacecraft. Particular interest should be paid to appendage retraction, since if this process is destabilizing, passive stabilization as proposed may not be useful. Over the past three decades, it has been an 'on-again-off-again affair' with the problem of spacecraft appendage deployment. In most instances, these studies have been numerical simulations of specific spacecraft configurations for which there were specific concerns. The primary focus of these studies was the behavior of the appendage during deployment; the effects of appendage retraction was considered only in one of these studies. What is missing in the literature is a thorough study of the effects of appendage deployment/retraction on the attitude stability of a spacecraft. This paper presents a rigorous analysis of the stability of a spinning spacecraft during the deployment or the retraction of an appendage. The analysis is simplified such that meaningful insights into the problem can be inferred; it is not overly simplified such that critical dynamical behavior is neglected. The system is analyzed assuming that the spacecraft hub is rigid. The appendage deployment mechanism is modeled as a point mass on a massless rod whose length undergoes prescribed changes. Simplified flexibility effects of the appendage are included. The system is examined for stability by linearizing the equations in terms of small deviations from steady, noninterfering coning motion. Routh's procedure for analyzing small deviations from steady motion in dynamical systems is utilized in the analysis. The system of equations are nondimensionalized to facilitate parametric studies. The results are presented in terms of a reduced number of nondimensional parameters so that some general conclusions may be drawn. Verification of the linear analysis is presented through numerical simulations of the complete nonlinear, nonautonomous, coupled equations.

Fitz-Coy, Norman↗

DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗

NASA NTRSDate not supplied

Aeroservoelastic Control Law Development for the Integrated Adaptive Wing Technology Maturation Wind-Tunnel Test

The NASA Advanced Air Transport Technology (AATT) Project seeks to “explore and develop technologies and concepts for improved energy efficiency and environmental compatibility for fixed wing subsonic transports” [1]. The IAWTM project is a subproject within AATT and is a joint Boeing/NASA effort to, in part, demonstrate the active control systems necessary to employ high-aspect ratio wings for improved aerodynamic efficiency. The aircraft configuration under study for the IAWTM project is a generic jet transport based on the NASA Common Research Model (CRM) whose wing aspect ratio has been increased to 13.5 (as opposed to nominal CRM aspect ratio of 9) and whose development is given in Ref. [2]. A semi-span model of the high-aspect-ratio CRM with 10 trailing edge control surfaces is currently being fabricated and will be delivered to the NASA Langley Transonic Dynamics Tunnel (TDT). There, several tests will be conducted to demonstrate maneuver load alleviation (MLA), gust load alleviation (GLA), active flutter suppression (AFS), and drag reduction, all at transonic conditions. A layout of the wind tunnel model is shown in Fig. 1. The wind-tunnel model is a spar/pod wing design with a relatively rigid fuselage attached to a mounting structure (seen in Fig. 1), which in turn is attached to a load balance behind the tunnel wall. Shims exist that can be attached to force the load to go through the fuselage rather than through the balance so as not to damage the balance in the event of flutter. The spar/pod design is similar to that used in the Boeing truss-braced wing aeroelastic model tested in the TDT in 2013 [3]. There is a central spar running the span of the wing to which the pods are attached that form the aerodynamic shape of the wing. The pods are constructed of ribs for attachment to the spar and skins to create the outer mold line. There are ten active control surfaces that will be used for the test objectives discussed above. Three of these (two outboard ailerons and one inboard aileron) are high-speed electrohydraulic actuators, and the remaining seven are low-speed mini-plane flaps, which have slower servo motors that will be used primarily for drag and maneuver load reductions. Control laws must be developed for each test (MLA, GLA, AFS, and drag reduction) and must be able to work in tandem when two or more control strategies are being tested simultaneously (e.g., when control laws to mitigate maneuver loads and gust loads operate simultaneously). The development of these control laws requires accurate mathematical representations of the aerodynamics, structural dynamics, and control surface dynamics of the system. Taken together, this mathematical representation—which is often in state-space form— is the ASE model. Three low-order ASE models of differing fidelities were developed for the IAWTM test and will be presented here. The differences in fidelity appeared in the aerodynamic modeling—one model assumed purely linear aerodynamics, while the other two models used different techniques to linearize the aerodynamics about a nonlinear mean. Furthermore, simple observer-based full-state feedback controllers for each test (MLA, GLA, AFS, and drag reduction) were designed from these low-order ASE models and implemented in a higher-order, nonlinear computational ASE simulation to assess the effectiveness of each controller. Observers were designed from each of the three lower-order ASE models, and each (linear) observer was then used as an estimate of the nonlinear dynamics of the system. Using the estimated states and the full-state feedback controller, control surface commands were fed back to the higher-order, nonlinear ASE system. The insight gained during this development will be applied to control law development for the physical IAWTM wind-tunnel test.

Josiah M Waite↗

DOE OSTI2022

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1980

Minimum fuel control of a vehicle with a continuously variable transmission

The design and evaluation of a control system for a sedan with a heat engine and a continuously variable transmission, is considered in a effort to minimize fuel consumption and achieve satisfactory dynamic response of vehicle variables as the vehicle is driven over a standard driving cycle. Even though the vehicle system was highly nonlinear, attention was restricted to linear control algorithms which could be easily understood and implemented demonstrated by simulation. Simulation results also revealed that the vehicle could exhibit unexpected dynamic behavior which must be taken into account in any control system design.

Burghart, J. H.↗

NASA NTRS1995

A class of stabilizing controllers for flexible multibody systems

The problem of controlling a class of nonlinear multibody flexible space systems consisting of a flexible central body to which a number of articulated appendages are attached is considered. Collocated actuators and sensors are assumed, and global asymptotic stability of such systems is established under a nonlinear dissipative control law. The stability is shown to be robust to unmodeled dynamics and parametric uncertainties. For a special case in which the attitude motion of the central body is small, the system, although still nonlinear, is shown to be stabilized by linear dissipative control laws. Two types of linear controllers are considered: static dissipative (constant gain) and dynamic dissipative. The static dissipative control law is also shown to provide robust stability in the presence of certain classes of actuator and sensor nonlinearities and actuator dynamics. The results obtained for this special case can also be readily applied for controlling single-body linear flexible space structures. For this case, a synthesis technique for the design of a suboptimal dynamic dissipative controller is also presented. The results obtained in this paper are applicable to a broad class of multibody and single-body systems such as flexible multilink manipulators, multipayload space platforms, and space antennas. The stability proofs use the Lyapunov approach and exploit the inherent passivity of such systems.

Joshi, Suresh M.↗