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.
Global Modeling and Assimilation Office (GMAO) is close to having a working 4DVAR system and has developed a linearized version of GEOS-5.This talk outlines a series of improvements made to the linearized dynamics, physics and trajectory.Of particular interest is the development of linearized cloud microphysics, which provides the framework for 'all-sky' data assimilation.
We consider problem of dynamics, control, and uncertainty quantification for quadcopter. We use the 6DOF model of quadcopter dynamics, linear quadratic regulator and linear quadratic Gaussian control of quadcopter in the presence of dynamical disturbances, measurement noise, hidden dynamical variables, dashing GPS signal, and wind gusts to predict quadcopter trajectory. We identify key sources of uncertainties and report on progress in development of a system that estimates the probability of safety-critical events using a set of algorithms based on the trajectory predictions.
The spontaneous activity of 40 otolith afferents and 44 canal afferents was recorded in 4 alert, intact squirrel monkeys. Polarization vectors and response properties of otolith afferents were determined during static re-orientations relative to gravity and during Earth-horizontal, sinusoidal, linear oscillations. Canal afferents were tested for sensitivity to linear accelerations. For regular otolith afferents, a significant correlation between upright discharge rate and sensitivity to dynamic acceleration in the horizontal plane was observed. This correlation was not present in irregular units. The sensitivity of otolith afferents to both static tilts and dynamic linear acceleration was much greater in irregularly discharging units than in regularly discharging units. The spontaneous activity and static and dynamic response properties of regularly discharging otolith afferents were similar to those reported in barbiturate-anesthetized squirrel monkeys. Irregular afferents also had similar dynamic response properties when compared to anesthetized monkeys. However, this sample of irregular afferents in alert animals had higher resting discharge rates and greater sensitivity to static tilts. The majority of otolith polarization vectors were oriented near the horizontal in the plane of the utricular maculae; however, directions of maximum sensitivity were different during dynamic and static testing. Canal afferents were not sensitive to static tilts or linear oscillations of the head.
The problem of identifying linear dynamical systems is studied by considering structural and deterministic properties of linear systems that have an impact on stochastic identification algorithms. In particular considered is parametrization of linear systems so that there is a unique solution and all systems in appropriate class can be represented. It is assumed that a parametrization of system matrices has been established from a priori knowledge of the system, and the question is considered of when the unknown parameters of this system can be identified from input/output observations. It is assumed that the transfer function can be asymptotically identified, and the conditions are derived for the local, global and partial identifiability of the parametrization. Then it is shown that, with the right formulation, identifiability in the presence of feedback can be treated in the same way. Similarly the identifiability of parametrizations of systems driven by unobserved white noise is considered using the results from the theory of spectral factorization.
Linearized dynamic equations and mechanical model for propellant sloshing in Lunar Excursion Module Descent Propulsion System /LEMDPS/
The topics considered are related to linear control system design, adaptive control, failure detection, control under failure, system reliability, and large-scale systems and decentralized control. It is pointed out that the design of a linear feedback control system which regulates a process about a desirable set point or steady-state condition in the presence of disturbances is a very important problem. The linearized dynamics of the process are used for design purposes. The typical linear-quadratic design involving the solution of the optimal control problem of a linear time-invariant system with respect to a quadratic performance criterion is considered along with gain reduction theorems and the multivariable phase margin theorem. The stumbling block in many adaptive design methodologies is associated with the amount of real time computation which is necessary. Attention is also given to the desperate need to develop good theories for large-scale systems, the beginning of a microprocessor revolution, the translation of the Wiener-Hopf theory into the time domain, and advances made in dynamic team theory, dynamic stochastic games, and finite memory stochastic control.
Smoothing equations for a linear continuous dynamic system with linear discrete measurements, derived from the discrete results of Rauch, Tung, and Striebel (1965), (R-T-S), are used to extend, through recursive updating, the previously published results of Bryson and Frazier (1963), (B-F), and yield a modified Bryson and Frazier, (M-B-F), algorithm. A comparison of the (M-B-F) and (R-T-S) algorithms leads to the conclusion that the former is to be preferred because it entails less computation, less storage, and less instability. It is felt that the presented (M-B-F) smoothing algorithm is a practical mechanization and should be of value in smoothing discretely observed dynamic linear systems.
A new approach for computing open-loop time-optimal control inputs for uncertain linear dynamical systems is developed. In particular, the single-axis, rest-to-rest maneuvering problem of flexible spacecraft in the presence of uncertainty in model parameters is considered. Robustified time-optimal control inputs are obtained by solving a parameter optimization problem subject to robustness constraints. A simple dynamical system with a rigid-body mode and one flexible mode is used to illustrate the concept.
The Kalman filter is the optimal filter in the presence of known gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions. Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz model as well as more realistic models of the means and atmosphere. A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter situations to allow for correct update of the ensemble members. The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to be quite puzzling in that results state estimates are worse than for their filter analogue. In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use the Lorenz model to test and compare the behavior of a variety of implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.
The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.
Multiloop linear-quadratic state-feedback (LQSF) regulators are shown to be robust against a variety of large dynamical linear time-invariant and memoryless nonlinear time-varying variations in open-loop dynamics. The results are interpreted in terms of the classical concepts of gain and phase margin, thus strengthening the link between classical and modern feedback theory.
The most medial of the deep cerebellar nuclei, the fastigial nucleus (FN), receives sensory vestibular information and direct inhibition from the cerebellar vermis. We investigated the signal processing in the primate FN by recording single-unit activities during translational motion, rotational motion, and eye movements. Firing rate modulation during horizontal plane translation in the absence of eye movements was observed in all non-eye-movement-sensitive cells and 26% of the pursuit eye-movement-sensitive neurons in the caudal FN. Many non-eye-movement-sensitive cells recorded in the rostral FN of three fascicularis monkeys exhibited convergence of signals from both the otolith organs and the semicircular canals. At low frequencies of translation, the majority of these rostral FN cells changed their firing rates in phase with head velocity rather than linear acceleration. As frequency increased, FN vestibular neurons exhibited a wide range of response dynamics with most cells being characterized by increasing phase leads as a function of frequency. Unlike cells in the vestibular nuclei, none of the rostral FN cells responded to rotational motion alone, without simultaneously exhibiting sensitivity to translational motion. Modulation during earth-horizontal axis rotation was observed in more than half (77%) of the neurons, although with smaller gains than during translation. In contrast, only 47% of the cells changed their firing rates during earth-vertical axis rotations in the absence of a dynamic linear acceleration stimulus. These response properties suggest that the rostral FN represents a main processing center of otolith-driven information for inertial motion detection and spatial orientation.
A NASA laboratory for evaluating linear detector arrays for remote sensing applications from the visible to shortwave infrared is described. Nominal requirements are presented on the critical parameters of signal-to-noise ratio, modulation transfer function, response linearity, dynamic range, and detector-to-detector response uniformity. Results of tests on a custom developed multi-spectral linear array from RCA and two commercial arrays from Texas Instruments and Reticon for the visible and near-infrared are presented. The custom array was clearly superior in signal-to-noise ratio and dynamic range to its off-the-shelf counterparts. Allthe arrays were found to have somewhat less than desirable modulation transfer functions, but excellent response linearity.
Coupling of the Reynolds-averaged Navier-Stokes equations, rigid-body dynamics, and a pitch attitude control law is demonstrated in two- and three-dimensions. The application problem was the separation of a canard-controlled store from an open-flow rectangular cavity bay at a freestream Mach number of 1.2. The transient flowfield was computed using a diagonal scheme in an overset mesh framework, with the resultant aerodynamic loads used as the forcing functions in the nonlinear dynamics equations. The proportional and rate gyro sensitivities were computed a priori using pole placement techniques for the linearized dynamical equations. These fixed gain values were used in the controller for the nonlinear simulation. Reasonable comparison between the full and linearized equations for a perturbed two-dimensional missile was found. Also in two-dimensions, a controlled store was found to possess improved separation characteristics over a canard-fixed store. In three-dimensions, trajectory comparisons with wind-tunnel data for the canard-fixed case will be made. In addition, it will be determined if a canard-controlled store is an effective means of improving cavity store separation characteristics.
The paper deals with the problem of system identifiability for a linear dynamical system. Two theorems are given relating the sufficient condition of system identifiability for certain linear structures to the total number of inputs and outputs. The principle of least squares, that seeks the minimization of a cost function is employed to carry out the system identification process. To illustrate the concept of the paper, a structural model of a beam with point masses is examined. Parameter identification methods are studied and a random search technique is introduced.
The problem addressed is that of obtaining reduced-order component models for use in simulating the dynamics of a multibody system. In certain cases, nonlinear system models may be constructed using linear dynamic models for each component, but allowing large angle motion between components. Without some form of model reduction, system models constructed in this manner may be too large for use in control system design and simulation trades. This paper analyzes one method of component model reduction that allows systems level requirements (e.g., capturing the effect of body 1 reaction wheel noise on body 2 camera pointing) to aid in the selection of the reduced-order component models. Briefly stated, important modes are selected at the system level and projected onto the components, and reduced-order components are then assembled into a reduced-order system model that retains the projected modes.
Invertibility of linear time invariant dynamical control systems, discussing L-inverses, sequential circuits, inherent integration, reproducibility, pointwise and functional inverses
An analytical and computerized study of the steady state and transient response of a phosphoric acid fuel cell (PAFC) system was completed. Parametric studies and sensitivity analyses of the PAFC system's operation were accomplished. Four non-linear dynamic models of the fuel cell stack, reformer, shift converters, and heat exchangers were developed based on nonhomogeneous non-linear partial differential equations, which include the material, component, energy balance, and electrochemical kinetic features. Due to a lack of experimental data for the dynamic response of the components only the steady state results were compared with data from other sources, indicating reasonably good agreement. A steady state simulation of the entire system was developed using, nonlinear ordinary differential equations. The finite difference method and trial-and-error procedures were used to obtain a solution. Using the model, a PAFC system, that was developed under NASA Grant, NCC3-17, was improved through the optimization of the heat exchanger network. Three types of cooling configurations for cell plates were evaluated to obtain the best current density and temperature distributions. The steady state solutions were used as the initial conditions in the dynamic model. The transient response of a simplified PAFC system, which included all of the major components, subjected to a load change was obtained. Due to the length of the computation time for the transient response calculations, analysis on a real-time computer was not possible. A simulation of the real-time calculations was developed on a batch type computer. The transient response characteristics are needed for the optimization of the design and control of the whole PAFC system. All of the models, procedures and simulations were programmed in Fortran and run on IBM 370 computers at Cleveland State University and the NASA Lewis Research Center.