Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Non linear dynamics”

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 307 records · Page 17

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Development of a non-linear simulation for generic hypersonic vehicles - ASUHS1

A nonlinear simulation is developed to model the longitudinal motion of a vehicle in hypersonic flight. The equations of motion pertinent to this study are presented. Analytic expressions for the aerodynamic forces acting on a hypersonic vehicle which were obtained from Newtonian Impact Theory are further developed. The control surface forces are further examined to incorporate vehicle elastic motion. The purpose is to establish feasible equations of motion which combine rigid body, elastic, and aeropropulsive dynamics for use in nonlinear simulations. The software package SIMULINK is used to implement the simulation. Also discussed are issues needing additional attention and potential problems associated with the implementation (with proposed solutions).

Salas, Juan↗

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

Videopanorama Frame Rate Requirements Derived from Visual Discrimination of Deceleration During Simulated Aircraft Landing

In order to determine the required visual frame rate (FR) for minimizing prediction errors with out-the-window video displays at remote/virtual airport towers, thirteen active air traffic controllers viewed high dynamic fidelity simulations of landing aircraft and decided whether aircraft would stop as if to be able to make a turnoff or whether a runway excursion would be expected. The viewing conditions and simulation dynamics replicated visual rates and environments of transport aircraft landing at small commercial airports. The required frame rate was estimated using Bayes inference on prediction errors by linear FRextrapolation of event probabilities conditional on predictions (stop, no-stop). Furthermore estimates were obtained from exponential model fits to the parametric and non-parametric perceptual discriminabilities d' and A (average area under ROC-curves) as dependent on FR. Decision errors are biased towards preference of overshoot and appear due to illusionary increase in speed at low frames rates. Both Bayes and A - extrapolations yield a framerate requirement of 35 < FRmin < 40 Hz. When comparing with published results [12] on shooter game scores the model based d'(FR)-extrapolation exhibits the best agreement and indicates even higher FRmin > 40 Hz for minimizing decision errors. Definitive recommendations require further experiments with FR > 30 Hz.

deceleration perception↗

Dispersion of ensembles of non-interacting particles

The dynamics of an ensemble of noninteracting particles dispersing from a common origin and moving in a common force field with an initial distribution of momenta is analyzed using an approach where the particles are considered as a continuum described by a phase-space distribution function. General solutions are obtained for both the distribution function and the associated spatial density function. The linear case of small departures from circular orbits in an axisymmetric gravitational field is treated along with the specific case of particle dispersion from an object in a circular orbit in the same type of field. Numerical results are presented for the latter case, and consideration is given to the inverse problem of determining the initial time and velocity distribution from knowledge of the ensemble structure at a later time. Explicit results are provided for the case of an ellipsoidal distribution of initial momenta, and a numerical procedure is indicated for treating more general cases.

Heard, W. B.↗

Vibration/Libration Interaction Dynamics During the Orbiter Based Deployment of Flexible Members

Essential features of a general formulation for studying librational dynamics of a large class of spacecraft during deployment of flexible members are reviewed. The formulation is applicable to a variety of missions ranging from deployment of antennas, booms and solar panels to manufacturing of trusses for space platforms using the space shuttle. The governing nonlinear, non-autonomous and coupled equations of motion are extremely difficult to solve even with the help of a computer, not to mention the cost involved. To get some appreciation as to the complex interactions between flexibility, deployment and attitude dynamics as well as to help pursue stability and control analysis, the equations are linearized about their nominal deflected equilibrium configuration. The procedure is applied to the Space Shuttle based deployment of boom and plate-like members. Results suggest substantial influence of the inertia parameter, flexural rigidity of the appendages, orbit eccentricity, deployment velocity, initial conditions, etc. on the system response. The results should prove useful in planning of the Orbiter based experiments aimed at assessing effectiveness of procedures for studying dynamics and control of flexible orbiting members.

Modi, V. J.↗

Novel Optical Diagnostic Techniques for Studying Particle Deposition Upon Large Cylinders in a Sheared Suspension

On a twelve-month voyage to Mars, one astronaut will require at least two tons of potable water and two tons of pure oxygen. Efficient, reliable fluid reclamation is therefore necessary for manned space exploration. Space habitats require a compact, flexible, and robust apparatus capable of solid-fluid mechanical separation over a wide range of fluid and particle densities and particle sizes. In space, centrifugal filtration, where particles suspended in fluid are captured by rotating fixed-fiber mat filters, is a logical candidate for mechanical separation. Non-colloidal particles are deposited on the fibers due to inertial impaction or direct interception. Since rotation rates are easily adjustable, inertial effects are the most practical way to control separation rates for a wide variety of multiphase mixtures in variable gravity environments. Understanding how fluid inertia and differential fluid-particle inertia, characterized by the Reynolds and Stokes numbers, respectively, affect deposition is critical in optimizing filtration in a microgravity environment. This work will develop non-intrusive optical diagnostic techniques for directly visualizing where and when non-colloidal particles deposit upon, or contact, solid surfaces: 'particle proximity sensors'. To model particle deposition upon a single filter fiber, these sensors will be used in ground-based experiments to study particle dynamics as in the vicinity of a large (compared with the particles) cylinder in a simply sheared (i.e., linearly-varying, zero-mean velocity profile) neutrally-buoyant, refractive-index matched solid-liquid suspension.

Yoda, M.↗

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass and momentum perturbations

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 11-12 July 1981 CCOPE case study indicated two episodes of coherent waves. While geostrophic adjustment, shearing instability, and terrain were all implicated separately or in combination as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to define the genesis processes from observations alone. The first part of this paper, 3D Numerical Modeling Studies of Terrain-Induced Mass/Momentum Perturbations, employs a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed gravity wave episode. The meso-beta scale numerical model is used to study various simulations of the role of multiple geostrophic adjustment processes in focusing a region for gravity wave genesis. The second part of this paper, Linear Theory and Theoretical Modeling, investigates the response of non-resting rotating homogeneous and continuously stratified Boussinesq models of the terrestrial atmosphere to temporally impulsive and uniformly propagating three-dimensional localized zonal momentum sources representative of midlatitude jet streaks. The methods of linear perturbation theory applied to the potential vorticity (PV) and wave field equations are used to study the geostrophic adjustment dynamics. The total zonal and meridional wind perturbations are separated into geostrophic and ageostrophic components in order to define and follow the evolution of both the primary and secondary mesocirculations accompanying midlatitude jetogenesis forced by geostrophic adjustment processes. This problem is addressed to help fill the gap in understanding the dynamics and structure of mesoscale inertia-gravity waves forced by geostrophic adjustment processes in simple two-dimensional quiescent current systems and those produced by mesoscale numerical models simulating the orographic and diabatic perturbation of three-dimensional quasi-geostrophically balanced synoptic scale jet streaks associated with complex baroclinic severe storm producing environments.

Lin, Yuh-Lang↗

Droplet Deformation in an Extensional Flow: The Role of Surfactant Physical Chemistry

Surfactant-induced Marangoni effects strongly alter the stresses exerted along fluid particle interfaces. In low gravity processes, these stresses can dictate the system behavior. The dependence of Marangoni effects on surfactant physical chemistry is not understood, severely impacting our ability to predict and control fluid particle flows. A droplet in an extensional flow allows the controlled study of stretching and deforming interfaces. The deformations of the drop allow both Marangoni stresses, which resist tangential shear, and Marangoni elasticities, which resist surface dilatation, to develop. This flow presents an ideal model system for studying these effects. Prior surfactant-related work in this flow considered a linear dependence of the surface tension on the surface concentration, valid only at dilute surface concentrations, or a non-linear framework at concentrations sufficiently dilute that the linear approximation was valid. The linear framework becomes inadequate for several reasons. The finite dimensions of surfactant molecules must be taken into account with a model that includes surfaces saturation. Nonideal interactions between adsorbed surfactant molecules alter the partitioning of surfactant between the bulk and the interface, the dynamics of surfactant adsorptive/desorptive exchange, and the sensitivity of the surface tension to adsorbed surfactant. For example, cohesion between hydrocarbon chains favors strong adsorption. Cohesion also slows the rate of desorption from interfaces, and decreases the sensitivity of the surface tension to adsorbed surfactant. Strong cohesive interactions result in first order surface phase changes with a plateau in the surface tension vs surface concentration. Within this surface concentration range, the surface tension is decoupled from surface concentration gradients. We are engaged in the study of the role of surfactant physical chemistry in determining the Marangoni stresses on a drop in an extensional flow in a numerical and experimental program. Using surfactants whose dynamics and equilibrium behavior have been characterized in our laboratory, drop deformation will be studied in ground-based experiment. In an accompanying numerical study, predictive drop deformations will be determined based on the isotherm and equation of state determined in our laboratory. This work will improve our abilities to predict and control all fluid particle flows.

Stebe, Kathleen J.↗

Turbomachinery Application of Lagrangian Dynamics to the Motion of Continuous Discrete Rotors

The stability/instability condition of a turbine rotor with axisymmetric supports is determined in the presence of gyroscopic loads and rub-induced destabilizing forces. A modal representation of the turbine engine is used, with one mode in each of the vertical and horizontal planes. The use of non-spinning rotor modes permits an explicit treatment of gyroscopic effects. The two linearized modal equations of motion of a rotor with axisymmetric supports are reduced to a single equation in a complex variable. The resulting eigenvalues yield explicit expressions at the stability boundary, for the whirl frequency as well as the required damping in the presence of the available rub-induced destabilization. Conversely, the allowable destabilization in the presence of the available damping is also given.

Source record↗

Non-linear guidance laws for automatic orbital rendezvous

In the automatic assembly of space stations or the automatic rendezvous and docking of a space-shuttle vehicle with a space station, it is of interest to determine suitable guidance laws for orbital rendezvous which are amenable to on-board implementation. In this paper, various forms of nonlinear guidance laws for automatic orbital rendezvous of two bodies represented by point masses in a central Newtonian gravitational field are derived. A basic requirement is that the approach be smooth and nonoscillatory. The global dynamic behavior of each of the resulting feedback systems is determined analytically and/or by computer simulation.

Wang, P. K. C.↗

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗

Calculation of steady and unsteady airfoil flow fields via the Navier-Stokes equations

A compressible time-dependent procedure for the two-dimensional ensemble averaged Navier-Stokes equations has been applied to the isolated airfoil problem in steady and unsteady flows. The procedure solves the governing equations via the linearized block implicit technique. Turbulence is modeled either via a mixing length or turbulence energy approach. The equations are solved in general non-orthogonal form with no-slip boundary conditions applied at the airfoil surface. Results are presented for airfoils at constant incidence, an airfoil in ramp motion and an airfoil oscillating through a dynamic stall loop. In general, steady converged solutions are obtained within 70 time steps over the range of Mach numbers considered. Comparisons with measured data show good agreement between computation and measurement.

Shamroth, S. J.↗

Dynamic response for thermal control and measurement and fast radiation thermometry

A preliminary evaluation was made by ORNL of a two-color ratio pyrometer (TCRP) for temperature control in the Modular Electromagnetic Levitation (MEL) experiment. A discussion was presented by Eric Spjut at the 1987 NASA Non-Contact Temperature Measurement Workshop (NASA Conf. Publ. 2503, pp. 182-213) in which he described the non-linear characteristics of the time response of TCPs. Researchers replicated his model and results and note that the non-linear response behavior is minimized for small temperature steps at high temperatures. They then used the predicted response in a model for a proportional or integral feedback controller and predicted the control characteristics for heating and cooling a 5-mm diameter sphere of niobium at high (1500 to 2750 K) temperatures. The analysis shows that for a slow (25-ms) time response for a commercial RCRP, overshoots of several hundred kelvins will result from a 100-K decrease in the setpoint, and temperature tracking errors of 14 to 45 K will occur for control temperature ramps of 1000K/s. For a fast (greater than 0.1 ms) time response, the overshoot and ramp response errors are largely eliminated.

Shepard, R. L.↗

On vibrational imperfection sensitivity of Augusti's model structure in the vicinity of a non-linear static state

In this paper we present a closed-form solution for vibrational imperfection sensitivity the effect of small imperfections on the vibrational frequencies of perturbed motion around the static equilibrium state of Augusti's model Structure (a rigid link, pinned at one end to a rigid foundation and supported at the other by a linear extensional spring that retains its horizontality, as the system deflects). We also treat a modified version of that model with attendant slightly different dynamics. It is demonstrated that the vibrational frequencies decreases as the initial imperfections increase.

Elishakoff, Isaac↗

Morphological evidence for local microcircuits in rat vestibular maculae

Previous studies suggested that intramacular, unmyelinated segments of vestibular afferent nerve fibers and their large afferent endings (calyces) on type I hair cells branch. Many of the branches (processes) contain vesicles and are presynaptic to type II hair cells, other processes, intramacular nerve fibers, and calyces. This study used serial section transmission electron microscopy and three-dimensional reconstruction methods to document the origins and distributions of presynaptic processes of afferents in the medial part of the adult rat utricular macula. The ultrastructural research focused on presynaptic processes whose origin and termination could be observed in a single micrograph. Results showed that calyces had 1) vesiculated, spine-like processes that invaginated type I cells and 2) other, elongate processes that ended on type II cells pre- as well as postsynaptically. Intramacular, unmyelinated segments of afferent nerve fibers gave origin to branches that were presynaptic to type II cells, calyces, calyceal processes, and other nerve fibers in the macula. Synapses with type II cells occurred opposite subsynaptic cisternae (C synapses); all other synapses were asymmetric. Vesicles were pleomorphic but were differentially distributed according to process origin. Small, clear-centered vesicles, approximately 40-60 nm in diameter, predominated in processes originating from afferent nerve fibers and basal parts of calyces. Larger vesicles approximately 70-120 nm in diameter having approximately 40-80 nm electron-opaque cores were dominant in processes originating from the necks of calyces. Results are interpreted to indicate the existence of a complex system of intrinsic feedforward (postsynaptic)-feedback (presynaptic) connections in a network of direct and local microcircuits. The morphological findings support the concept that maculae dynamically preprocess linear acceleratory information before its transmission to the central nervous system.

NASA Center ARC↗

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↗

Large-scale galaxy flow from a non-gravitational impulse

A theory is presented describing linear perturbations of an expanding universe containing multiple, independently perturbed, collisionless, gravitationally coupled constituents. Solutions are found in the limit where one initially unperturbed component dominates the total density. The theory is applied to perturbations generated by a nongravitational process in one or more of the light components, as would occur in explosive or radiation-pressure-instability theories of galaxy formation. The apparent dynamical density parameter and correlations between density and velocity amplitude for various populations, are evaluated as a function of cosmic scale factor.

Hogan, Craig J.↗