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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Flight-Determined Subsonic Longitudinal Stability and Control Derivatives of the F-18 High Angle of Attack Research Vehicle (HARV) with Thrust Vectoring

The subsonic longitudinal stability and control derivatives of the F-18 High Angle of Attack Research Vehicle (HARV) are extracted from dynamic flight data using a maximum likelihood parameter identification technique. The technique uses the linearized aircraft equations of motion in their continuous/discrete form and accounts for state and measurement noise as well as thrust-vectoring effects. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics over the aircraft, particularly at high angles of attack. Thrust vectoring was implemented using electrohydraulically-actuated nozzle postexit vanes and a specialized research flight control system. During maneuvers, a control system feature provided independent aerodynamic control surface inputs and independent thrust-vectoring vane inputs, thereby eliminating correlations between the aircraft states and controls. Substantial variations in control excitation and dynamic response were exhibited for maneuvers conducted at different angles of attack. Opposing vane interactions caused most thrust-vectoring inputs to experience some exhaust plume interference and thus reduced effectiveness. The estimated stability and control derivatives are plotted, and a discussion relates them to predicted values and maneuver quality.

Iliff, Kenneth W.↗
NASA NTRS1989

Manipulation strategies for massive space payloads

Control for the bracing strategy is being examined. It was concluded earlier that trajectory planning must be improved to best achieve the bracing motion. Very interesting results were achieved which enable the inverse dynamics of flexible arms to be calculated for linearized motion in a more efficient manner than previously published. The desired motion of the end point beginning at t=0 and ending at t=t sub f is used to calculate the required torque at the joint. The solution is separated into a causal function that is zero for t is less than 0 and an accusal function which is zero for t is greater than t sub f. A number of alternative end point trajectories were explored in terms of the peak torque required, the amount of anticipatory action, and other issues. The single link case is the immediate subject and an experimental verification of that case is being performed. Modeling with experimental verification of closed chain dynamics continues. Modeling effort has pointed out inaccuracies that result from the choice of numerical techniques used to incorporate the closed chain constraints when modeling our experimental prototype RALF (Robotic Arm Large and Flexible). Results were compared to TREETOPS, a multi body code. The experimental verification work is suggesting new ways to make comparisons with systems having structural linearity and joint and geometric nonlinearity. The generation of inertial forces was studied with a small arm that will damp the large arm's vibration.

Book, Wayne J.↗
DOE OSTI2021

A snapshot review—Fluctuations in quantum materials: from skyrmions to superconductivity

Abstract By measuring a linear response function directly, such as the dynamic susceptibility, one can understand fundamental material properties. However, a fresh perspective can be offered by studying fluctuations. This can be related back to the dynamic susceptibility through the fluctuation–dissipation theorem, which relates the fluctuations in a system to its response, an alternate route to access the physics of a material. Here, we describe a new X-ray tool for material characterization that will offer an opportunity to uncover new physics in quantum materials using this theorem. We provide details of the method and discuss the requisite analysis techniques in order to capitalize on the potential to explore an uncharted region of phase space. This is followed by recent results on a topological chiral magnet, together with a discussion of current work in progress. We provide a perspective on future measurements planned for work in unconventional superconductivity. Graphical abstract We describe a new X-ray tool for material characterization that will offer an opportunity to uncover new physics in quantum materials using coherent, short-pulsed X-rays. We provide details of the method and discuss the requisite analysis techniques in order to capitalize on the potential to explore an uncharted region of phase space. This is followed by recent results on a topological chiral magnet, together with a discussion of current work in progress. We provide a perspective on future measurements planned for work in unconventional superconductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Learning thermodynamic master equations for open quantum systems

The characterization of Hamiltonians and other components of open quantum dynamical systems plays a crucial role in quantum computing and other applications. Scientific machine learning techniques have been applied to this problem in a variety of ways, including by modeling with deep neural networks. However, the majority of mathematical models describing open quantum systems are linear, and the natural nonlinearities in learnable models have not been incorporated using physical principles. We present a data-driven model for open quantum systems that includes learnable, thermodynamically consistent terms. The trained model is interpretable, as it directly estimates the system Hamiltonian and linear components of coupling to the environment. We validate the model on synthetic two and three-level data, as well as experimental two-level data collected from a quantum device at Lawrence Livermore National Laboratory.

Mathematics and Computing↗
NASA NTRS1975

Dynamic analysis of a flexible spacecraft with rotating components. Volume 1: Analytical developments

Analytical procedures and digital computer code are presented for the dynamic analysis of a flexible spacecraft with rotating components. Topics, considered include: (1) nonlinear response in the time domain, and (2) linear response in the frequency domain. The spacecraft is assumed to consist of an assembly of connected rigid or flexible subassemblies. The total system is not restricted to a topological connection arrangement and may be acting under the influence of passive or active control systems and external environments. The analytics and associated digital code provide the user with the capability to establish spacecraft system nonlinear total response for specified initial conditions, linear perturbation response about a calculated or specified nominal motion, general frequency response and graphical display, and spacecraft system stability analysis.

Bodley, C. S.↗
NASA NTRS1976

A spectral identification technique for adaptive attitude control and pointing of the Space Telescope

The Space Telescope is a 2.4 m class aperture optical telescope having near-diffraction-limited performance. It will be placed into earth orbit by 1980 via the Space Shuttle. The problem considered is how to achieve negligible degradation of the astronomy imaging capability (to 0.005 arc second) due to smearing by pointing motions during observations. Initially, pointing instability sources were identified and a linear stability was used to assess the magnitude of elastic body modes and to design control system compensation regions necessary for subsequent adaptive control. A spectral identification technique for this adaptive attitude control and pointing has been investigated that will alleviate requirements for comprehensive dynamic ground testing. Typical all-digital simulation results describing motions of the telescope line of sight are presented.

Teuber, D. L.↗
NASA NTRS1990

Recursive calculation of geared robot manipulator dynamics

A recursive formulation is presented for the calculation of the inverse and forward dynamics of rigid robot manipulators with gear systems on each joint. The complete effects of the gear ratios and the gyroscopic effects of the spinning motor/gear are included in the recursive formulation. The forward dynamics solution recursively calculates the joint accelerations when given motor torques, and the number of computations grows linearly with the number of links.

Murphy, Steve H.↗
NASA NTRS1991

A comparison of controller designs for an experimental flexible structure

Control systems design and hardware testing are addressed for an experimental structure that displays the characteristics of a typical flexible spacecraft. The results of designing and implementing various control design methodologies are described. The design methodologies under investigation include linear quadratic Gaussian control, static and dynamic dissipative controls, and H-infinity optimal control. Among the three controllers considered, it is shown, through computer simulation and laboratory experiments on the evolutionary structure, that the dynamic dissipative controller gave the best results in terms of vibration suppression and robustness with respect to modeling errors.

Lim, K. B.↗
NASA NTRS1992

Large-amplitude finite element flutter analysis of composite panels in hypersonic flow

Consideration is given to a finite-element approach for determining the nonlinear flutter characteristics of 3D thin laminated composite panels using the full third-order-piston transverse-loading aerodynamic theory. The unsteady hypersonic aerodynamic theory and the von Karman large-deflection-plate theory are used to formulate the aeroelasticity problem. Nonlinear flutter analyses are performed to assess the influence of the higher-order aerodynamic theory on the structure's limit-cycle amplitude and the dynamic pressure of the flow velocity. A solution procedure is presented to solve the nonlinear panel flutter and large-amplitude free-vibration finite-element equations. Nonlinear flutter analyses are performed for different boundary support-conditions and for various system parameters. Linear finite-element flutter for isotropic and composite panels and large-amplitude isotropic panel flutter results are compared with existing classical solutions. The large-amplitude panel flutter results using the full third-order piston aerodynamic theory are presented to assess the influence of the nonlinear aerodynamic theory.

Gray, Carl E., Jr.↗
NASA NTRS2003

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗
NASA NTRS2006

Control Surface Interaction Effects of the Active Aeroelastic Wing Wind Tunnel Model

This paper presents results from testing the Active Aeroelastic Wing wind tunnel model in NASA Langley s Transonic Dynamics Tunnel. The wind tunnel test provided an opportunity to study aeroelastic system behavior under combined control surface deflections, testing for control surface interaction effects. Control surface interactions were observed in both static control surface actuation testing and dynamic control surface oscillation testing. The primary method of evaluating interactions was examination of the goodness of the linear superposition assumptions. Responses produced by independently actuating single control surfaces were combined and compared with those produced by simultaneously actuating and oscillating multiple control surfaces. Adjustments to the data were required to isolate the control surface influences. Using dynamic data, the task increases, as both the amplitude and phase have to be considered in the data corrections. The goodness of static linear superposition was examined and analysis of variance was used to evaluate significant factors influencing that goodness. The dynamic data showed interaction effects in both the aerodynamic measurements and the structural measurements.

Heeg, Jennifer↗
NASA NTRS2024

Revolutionary Vertical Lift Technology Lift + Cruise Concept Vehicle Powertrain Dynamic Model

The Lift + Cruise (L+C) is one of several Revolutionary Vertical Lift Technology (RVLT) concept aircraft identified by NASA to investigate Urban Air Mobility (UAM) requirements. This six-passenger vehicle has a conventional tube-and-wing configuration, but with eight lift rotors (four on each wing) and a pusher rotor on the tail. This arrangement allows it to take off and land vertically as well as fly like an airplane. This report presents a novel dynamic model of the L+C powertrain built using the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS) and the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), both of which are NASA-developed toolboxes, and together are capable of modeling hybrid-electric powertrain architectures. This research effort models the complex behavior of the powertrain at the mechanical timescale (15 to 20 ms). The model consists of a turboshaft engine, generator, direct current bus, energy storage management system, gearboxes, and electric propulsor systems coupled together to create a dynamic model that mimics the operation of the multidomain powertrain. Two control algorithms, designed to equalize the power consumed to the power produced, are compared to assess their practicality: (1) independent proportional-integral (PI) control and (2) linear quadratic integral regulator (LQIR) control. Simulation results highlight the key features present in the model and demonstrate its operation under the two control schemes.

revolutionary vertical lift technology↗
NASA NTRS1989

Computational fluid dynamics on a massively parallel computer

A finite difference code was implemented for the compressible Navier-Stokes equations on the Connection Machine, a massively parallel computer. The code is based on the ARC2D/ARC3D program and uses the implicit factored algorithm of Beam and Warming. The codes uses odd-even elimination to solve linear systems. Timings and computation rates are given for the code, and a comparison is made with a Cray XMP.

Jespersen, Dennis C.↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
NASA NTRS1989

A methodology for formulating a minimal uncertainty model for robust control system design and analysis

In the design and analysis of robust control systems for uncertain plants, the technique of formulating what is termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents the transfer function matrix M(s) of the nominal system, and delta represents an uncertainty matrix acting on M(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unstructured uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, and for real parameter variations the diagonal elements are real. As stated in the literature, this structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the literature addresses methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty. Since have a delta matrix of minimum order would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta model would be useful. A generalized method of obtaining a minimal M-delta structure for systems with real parameter variations is given.

Belcastro, Christine M.↗
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
NASA NTRS2022

How Well do Multisatellite Products Capture the Space-Time Dynamics of Precipitation? Part II: Building an Error Model Through Spectral System Identification

Satellite precipitation products, as all quantitative estimates, come with some inherent degree of uncertainty. To associate a quantitative value of the uncertainty to each individual estimate, error modeling is necessary. Most of the error models proposed so far compute the uncertainty as a function of precipitation intensity only, and only at one specific spatio-temporal scale. We propose a spectral error model which accounts for the neighboring space-time dynamics of precipitation into the uncertainty quantification. Systematic distortions of the precipitation signal and random errors are characterized distinctively in every frequency-wavenumber band in the Fourier domain, to accurately characterize error across scales. The systematic distortions are represented as a deterministic space-time linear filtering term. The random errors are represented as a non-stationary additive noise. The spectral error model is applied to the IMERG multi satellite precipitation product and its parameters are estimated empirically through a system identification approach using the GV-MRMS gauge-radar measurements as reference (“truth”) over the eastern United States. The filtering term is found to be essentially low-pass. While traditional error models attribute most of the error variance to random errors, it is found here that the systematic filtering term explains 48% of the error variance at the native resolution of IMERG. This fact confirms that, at high resolution, filtering effects in satellite precipitation products cannot be ignored, and that the error cannot be represented as a purely random additive or multiplicative term. An important consequence is that precipitation estimates derived from totally different sources shall not be expected to automatically have statistically independent errors.

Precipitation↗
NASA NTRS2013

A Short-Range Distance Sensor with Exceptional Linearity

A sensor has been demonstrated that can measure distance over a total range of about 300 microns to an accuracy of about 0.1 nm (resolution of about 0.01 nm). This represents an exceptionally large dynamic range of operation - over 1,000,000. The sensor is optical in nature, and requires the attachment of a mirror to the object whose distance is being measured. This work resulted from actively developing a white light interferometric system to be used to measure the depths of defects in the Space Shuttle Orbiter windows. The concept was then applied to measuring distance. The concept later expanded to include spectrometer calibration. In summary, broadband (i.e., white) light is launched into a Michelson interferometer, one mirror of which is fixed and one of which is attached to the object whose distance is to be measured. The light emerging from the interferometer has traveled one of two distances: either the distance to the fixed mirror and back, or the distance to the moving mirror and back. These two light beams mix and produce an interference pattern where some wavelengths interfere constructively and some destructively. Sending this light into a spectrometer allows this interference pattern to be analyzed, yielding the net distance difference between the two paths. The unique feature of this distance sensor is its ability to measure accurately distance over a dynamic range of more than one million, the ratio of its range (about 300 microns) to its accuracy (about 0.1 nanometer). Such a large linear operating range is rare and arises here because both amplitude and phase-matching algorithms contribute to the performance. The sensor is limited by the need to attach a mirror of some kind to the object being tracked, and by the fairly small total range, but the exceptional dynamic range should make it of interest.

Simmons, Steven↗