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,225 records · Page 68

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗

Stability Assessment and Tuning of an Adaptively Augmented Classical Controller for Launch Vehicle Flight Control

Recently, a robust and practical adaptive control scheme for launch vehicles [ [1] has been introduced. It augments a classical controller with a real-time loop-gain adaptation, and it is therefore called Adaptive Augmentation Control (AAC). The loop-gain will be increased from the nominal design when the tracking error between the (filtered) output and the (filtered) command trajectory is large; whereas it will be decreased when excitation of flex or sloshing modes are detected. There is a need to determine the range and rate of the loop-gain adaptation in order to retain (exponential) stability, which is critical in vehicle operation, and to develop some theoretically based heuristic tuning methods for the adaptive law gain parameters. The classical launch vehicle flight controller design technics are based on gain-scheduling, whereby the launch vehicle dynamics model is linearized at selected operating points along the nominal tracking command trajectory, and Linear Time-Invariant (LTI) controller design techniques are employed to ensure asymptotic stability of the tracking error dynamics, typically by meeting some prescribed Gain Margin (GM) and Phase Margin (PM) specifications. The controller gains at the design points are then scheduled, tuned and sometimes interpolated to achieve good performance and stability robustness under external disturbances (e.g. winds) and structural perturbations (e.g. vehicle modeling errors). While the GM does give a bound for loop-gain variation without losing stability, it is for constant dispersions of the loop-gain because the GM is based on frequency-domain analysis, which is applicable only for LTI systems. The real-time adaptive loop-gain variation of the AAC effectively renders the closed-loop system a time-varying system, for which it is well-known that the LTI system stability criterion is neither necessary nor sufficient when applying to a Linear Time-Varying (LTV) system in a frozen-time fashion. Therefore, a generalized stability metric for time-varying loop=gain perturbations is needed for the AAC.

VanZwieten, Tannen↗

Finite element appendage equations for hybrid coordinate dynamic analysis.

The increasingly common practice of idealizing a spacecraft as a collection of interconnected rigid bodies to some of which are attached linearly elastic flexible appendages leads to equations of motion expressed in terms of a combination of discrete coordinates describing the arbitrary rotational motions of the rigid bodies and distributed or modal coordinates describing the small, time-varying deformations of the appendages: such a formulation is said to employ a hybrid system of coordinates. In the present paper the existing literature is extended to provide hybrid coordinate equations of motion for a finite element model of a flexible appendage attached to a rigid base undergoing unrestricted motions and some of the advantages of the finite element approach are noted. Transformations to the modal coordinates appropriate for the general case are provided.

Likins, P. W.↗

Numerical methods for multibody systems

This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented.

Glowinski, Roland↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

An application of high authority/low authority control and positivity

Control Dynamics Company (CDy), in conjunction with NASA Marshall Space Flight Center (MSFC), has supported the U.S. Air Force Wright Aeronautical Laboratory (AFWAL) in conducting an investigation of the implementation of several DOD controls techniques. These techniques are to provide vibration suppression and precise attitude control for flexible space structures. AFWAL issued a contract to Control Dynamics to perform this work under the Active Control Technique Evaluation for Spacecraft (ACES) Program. The High Authority Control/Low Authority Control (HAC/LAC) and Positivity controls techniques, which were cultivated under the DARPA Active Control of Space Structures (ACOSS) Program, were applied to a structural model of the NASA/MSFC Ground Test Facility ACES configuration. The control systems design were accomplished and linear post-analyses of the closed-loop systems are provided. The control system designs take into account effects of sampling and delay in the control computer. Nonlinear simulation runs were used to verify the control system designs and implementations in the facility control computers. Finally, test results are given to verify operations of the control systems in the test facility.

Seltzer, S. M.↗

Dynamic loading of spur gears with linear or parabolic tooth profile modification

A computer simulation was conducted to investigate the effects of both linear and parabolic tooth profile modification on the dynamic response of low-contact-ratio spur gears. The effect of the total amount of modification and the length of the modification zone were studied at various loads and speeds to find the optimal profile modification for minimal dynamic loading. Design charts consisting of normalized maximum dynamic load curves were generated for gear systems operated at various loads and with different tooth profile modification. An optimum profile modification can be determined from these design charts to minimize the dynamic loads of spur gear systems.

Lin, Hsiang Hsi↗

Recursive Newton-Euler formulation of manipulator dynamics

A recursive Newton-Euler procedure is presented for the formulation and solution of manipulator dynamical equations. The procedure includes rotational and translational joints and a topological tree. This model was verified analytically using a planar two-link manipulator. Also, the model was tested numerically against the Walker-Orin model using the Shuttle Remote Manipulator System data. The hinge accelerations obtained from both models were identical. The computational requirements of the model vary linearly with the number of joints. The computational efficiency of this method exceeds that of Walker-Orin methods. This procedure may be viewed as a considerable generalization of Armstrong's method. A six-by-six formulation is adopted which enhances both the computational efficiency and simplicity of the model.

Nasser, M. G.↗

Dynamic loading of spur gears with linear or parabolic tooth profile modifications

A computer simulation was conducted to investigate the effects of both linear and parabolic tooth profile modification on the dynamic response of low-contact-ratio spur gears. The effect of the total amount of modification and the length of the modification zone were studied at various loads and speeds to find the optimal profile modification for minimal dynamic loading. Design charts consisting of normalized maximum dynamic load curves were generated for gear systems operated at various loads and with different tooth profile modification. An optimum profile modification can be determined from these design charts to minimize the dynamic loads of spur gear systems.

Lin, H. H.↗

Hybrid state-space self-tuning control of uncertain linear systems

The paper presents a hybrid state-space self-tuner using a new dual-rate sampling scheme for digital adaptive control of continuous-time uncertain linear systems. A state-space-based recursive least-squares algorithm, together with a variable forgetting factor, is used for direct estimations of both the equivalent discrete-time uncertain linear system parameters and the associated discrete-time state of a continuous-time uncertain linear system from the sampled input and output data. An analogue optimal regional pole-placement design method is used for designing an optimal observer-based analogue controller. A suboptimal observer-based digital controller is then designed from the designed analogue controller using digital redesign technique. To enhance the robustness of parameter identification and state estimation algorithms, a dynamic bound for a class of uncertain bilinear parameters and a fast-rate digital controller are developed at each fast-sampling period. Also, to accommodate computation loads and computation delay for developing the advanced hybrid self-tuner, the designed analogue controller and observer gains are both updated at each slow-sampling period. This control technique has been successfully applied to benchmark control problems.

Shieh, L. S.↗

Diminishing sign anomaly and scaling behavior of the mixed-state Hall resistivity in Tl2Ba2Ca2Cu3O10 films containing columnar defects

The issues of sign reversal of the Hall voltage and scaling between longitudinal (rho(xx)) and Hall (rho(xy)) resistivities are studied in Tl2Ba2Ca2Cu3O10 films in which the vortex dynamics is drastically changed by flux pinning at heavy-ion-irradiation-induced linear defects. While the sign anomaly diminishes with increasing defect concentration, the power law rho(xy) is approximately equal to rho(xx) exp beta, beta = 1.85 +/- 0.1, holds even after irradiation. This result shows that the scaling is a universal feature of the mixed state in this system. The sign anomaly, on the other hand, is not consistent with a model that invokes pinning-induced backflow in the vortex core as the mechanism for this effect.

Budhani, R. C.↗

Geometry of thin liquid sheet flows

Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Chubb, Donald L.↗

A Multi-Resolution Nonlinear Mapping Technique for Design and Analysis Application

This report describes a nonlinear mapping technique where the unknown static or dynamic system is approximated by a sum of dimensionally increasing functions (one-dimensional curves, two-dimensional surfaces, etc.). These lower dimensional functions are synthesized from a set of multi-resolution basis functions, where the resolutions specify the level of details at which the nonlinear system is approximated. The basis functions also cause the parameter estimation step to become linear. This feature is taken advantage of to derive a systematic procedure to determine and eliminate basis functions that are less significant for the particular system under identification. The number of unknown parameters that must be estimated is thus reduced and compact models obtained. The lower dimensional functions (identified curves and surfaces) permit a kind of "visualization" into the complexity of the nonlinearity itself.

Phan, Minh Q.↗

A Multi-Resolution Nonlinear Mapping Technique for Design and Analysis Applications

This report describes a nonlinear mapping technique where the unknown static or dynamic system is approximated by a sum of dimensionally increasing functions (one-dimensional curves, two-dimensional surfaces, etc.). These lower dimensional functions are synthesized from a set of multi-resolution basis functions, where the resolutions specify the level of details at which the nonlinear system is approximated. The basis functions also cause the parameter estimation step to become linear. This feature is taken advantage of to derive a systematic procedure to determine and eliminate basis functions that are less significant for the particular system under identification. The number of unknown parameters that must be estimated is thus reduced and compact models obtained. The lower dimensional functions (identified curves and surfaces) permit a kind of "visualization" into the complexity of the nonlinearity itself.

Phan, Minh Q.↗

Adjusting a Full Body Model to Mitigate Inverse Kinematics Artifacts in Opensim

BACKGROUND: In support of Vibration Isolation and Stabilization (VIS) system development for Human Health Countermeasures (HHC) exercise systems in space, such as the European Enhanced Exploration Exercise Device (E4D) [1], dynamic quantities required to model the response of a proposed VIS while considering the effect of VIS motion on the forces between the human and VIS platform were obtained using motion capture data [2]. On occasion, large-amplitude oscillatory spikes were found in the subject’s linear and angular momentum derivatives, affecting analyses that depend on forces and moments derived from motion capture. The purpose of this investigation was to identify causes of these artifacts and techniques for their resolution. METHODS AND RESULTS: To obtain the required human dynamic quantities to drive the VIS simulation, motion capture data was collected containing recorded trajectories of passive retroreflective markers on body landmarks of an exercising subject. Since the full body Rajagopal model [3] was originally used to enhance gait analysis, upper body joints did not require large Ranges of Motion (ROM). We thus modified the Rajagopal model [4, 5] to allow it to be used for upper body intensive exercises like those common to the E4D. OpenSim Inverse Kinematics (IK) [6] was performed using these scaled subject models to generate the joint angles throughout the exercise while minimizing marker error. At times, the arms were observed to ‘snap’ from one configuration to another, causing spike artifacts. Following IK, a custom OpenSim plugin [7] was used to determine the required dynamic quantities including the rates of change of the linear and angular momenta of the human. Since motion capture is recorded at a larger time step than required by the VIS simulation, the human center of mass location was fit with splines and a second derivative taken to obtain the momentum derivative, allowing the VIS simulation to maintain conservation of momentum when appropriate. In a few cases during this stage, artifacts much larger than the expected noise of the second derivatives were introduced. Investigation of cases containing artifacts revealed several modifications that could be made to the OpenSim model to improve IK results. Since OpenSim models use Euler angles and rotation sequences, ‘gimbal lock’ would be encountered in the arms when raised 90 degrees to the side (e.g., T-pose, some hang clean exercise, etc.). This was resolved by reorienting the horizontal axes at the shoulder joint by 45 degrees, placing ‘gimbal lock’ outside common arm ROM, with the arm ROMs adjusted following this change. Elbow and wrist ROMs could also be adjusted to allow realistic motion while at the same time limiting the likelihood of unrealistic orientations. On occasion, the arms flipped backwards when raised above the head. This was prevented by using medial elbow markers in scaling and IK. When medial markers were not available, the acromial joint location in the unscaled model was shifted before model scaling to better align the arm with available markers. Lastly, artifacts which became apparent after taking the second derivatives of spline-fit data were found to occur in cases when the pelvis rotation limit prevented the full range of motion of an exercise. These issues were resolved by unclamping the pelvis rotation limit. Through this investigation, an understanding of conditions leading to IK artifacts was acquired allowing the automation of artifact detection. These artifact detection and mitigation techniques can be applied toward modeling of upper body motions in aerospace and other fields for improved IK results.

C A Bell↗

A parametric LQ approach to multiobjective control system design

The synthesis of a constant parameter output feedback control law of constrained structure is set in a multiple-objective linear quadratic regulator (MOLQR) framework. The use of intuitive objective functions such as model-following ability and closed-loop trajectory sensitivity allows multiple-objective decision-making techniques, such as the surrogate worth tradeoff method, to be applied. For the continuous-time deterministic problem with an infinite time horizon, dynamic compensators as well as static output feedback controllers can be synthesized using a descent Anderson-Moore algorithm modified to impose linear equality constraints on the feedback gains by moving in feasible directions. Results of three different examples are presented, including a unique reformulation of the sensitivity reduction problem.

Kyr, Doug↗

The Sun's Impact on Climate

We provide an overview of the impact of the Sun on the Earth atmosphere and climate system, focused on heating of Earth's atmosphere and oceans. We emphasize the importance of the spectral measurements of SIM and SOLSTICE- that we must know how solar variations are distributed over ultraviolet, visible, and infrared wavelengths, since these have separate characteristic influences on Earth's ozone layer, clouds, and upper layers of the oceans. Emphasis is also given to understanding both direct and indirect influences of the Sun on the Earth, which involve feedbacks between Earth's stratosphere, troposphere, and oceans, each with unique time scales, dynamics, chemistry, and biology, interacting non-linearly. Especially crucial is the role of all three phases of water on Earth, water vapor being the primary greenhouse gas in the atmosphere, the importance of trace gases such as CO2 arising from their absorption in the "water vapor window" at 800 - 1250/cm (12.5 to 8 microns). Melting of polar ice is one major response to the post-industrial global warming, enhanced due to "ice-albedo" feedback. Finally, water in liquid form has a major influence due to cloud albedo feedback, and also due to the oceans' absorption of solar radiation, particularly at visible wavelengths, through the visible "liquid water window" that allows penetration of visible light deep into the mixed layer, while nearby ultraviolet and infrared wavelengths do not penetrate past the upper centimeter ocean surface skin layer. A large fraction of solar energy absorbed by the oceans goes into the latent heat of evaporation. Thus the solar heating of the atmosphere-ocean system is strongly coupled through the water cycle of evaporation, cloud formation, precipitation, surface runoff and ice formation, to Earth's energy budget and climate, each different climate component responding to variations in different solar spectral bands, at ultraviolet, visible and infrared wavelengths.

Cahalan, Robert↗