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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Hydrodynamic Instability and Enhanced Transport in Protoplanetary Nebulae

The nature of turbulence (and the enhanced transport it provides) is a key element to comprehend the dynamics, physics and chemistry of the protoplanetary nebulae and consequently the planet formation process. Early accretion disk models postulated the turbulent transport through the well-known "alpha-viscosity" model, introduced by Shakura and Sunyaev in 1973. Since then, the nature of the turbulence in disks has been a subject of investigation. In 1991, the rediscovery by Balbus and Hawley of Chandrasekhar's linear instability in magnetized disks was a breakthrough in the discipline. Unfortunately, the mechanisms leading to turbulence in non-magnetized disks, such as protoplanetary nebulae, remain poorly understood. We will present results from laboratory experiments along with analytical arguments showing that, despite skepticism in the Astrophysical community, differential rotation may indeed be sufficient to trigger and sustain turbulence. We will also propose an alternative viscosity prescription derived from both experiments and analysis.

Richard, Denis T.↗
NASA NTRS2009

Three-dimensional Finite Element Formulation and Scalable Domain Decomposition for High Fidelity Rotor Dynamic Analysis

This paper has two objectives. The first objective is to formulate a 3-dimensional Finite Element Model for the dynamic analysis of helicopter rotor blades. The second objective is to implement and analyze a dual-primal iterative substructuring based Krylov solver, that is parallel and scalable, for the solution of the 3-D FEM analysis. The numerical and parallel scalability of the solver is studied using two prototype problems - one for ideal hover (symmetric) and one for a transient forward flight (non-symmetric) - both carried out on up to 48 processors. In both hover and forward flight conditions, a perfect linear speed-up is observed, for a given problem size, up to the point of substructure optimality. Substructure optimality and the linear parallel speed-up range are both shown to depend on the problem size as well as on the selection of the coarse problem. With a larger problem size, linear speed-up is restored up to the new substructure optimality. The solver also scales with problem size - even though this conclusion is premature given the small prototype grids considered in this study.

Datta, Anubhav↗
NASA NTRS1996

Relating the Hadamard Variance to MCS Kalman Filter Clock Estimation

The Global Positioning System (GPS) Master Control Station (MCS) currently makes significant use of the Allan Variance. This two-sample variance equation has proven excellent as a handy, understandable tool, both for time domain analysis of GPS cesium frequency standards, and for fine tuning the MCS's state estimation of these atomic clocks. The Allan Variance does not explicitly converge for the nose types of alpha less than or equal to minus 3 and can be greatly affected by frequency drift. Because GPS rubidium frequency standards exhibit non-trivial aging and aging noise characteristics, the basic Allan Variance analysis must be augmented in order to (a) compensate for a dynamic frequency drift, and (b) characterize two additional noise types, specifically alpha = minus 3, and alpha = minus 4. As the GPS program progresses, we will utilize a larger percentage of rubidium frequency standards than ever before. Hence, GPS rubidium clock characterization will require more attention than ever before. The three sample variance, commonly referred to as a renormalized Hadamard Variance, is unaffected by linear frequency drift, converges for alpha is greater than minus 5, and thus has utility for modeling noise in GPS rubidium frequency standards. This paper demonstrates the potential of Hadamard Variance analysis in GPS operations, and presents an equation that relates the Hadamard Variance to the MCS's Kalman filter process noises.

Hutsell, Steven T.↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗
NASA NTRS2005

Cycle-averaged dynamics of a periodically driven, closed-loop circulation model

Time-varying elastance models have been used extensively in the past to simulate the pulsatile nature of cardiovascular waveforms. Frequently, however, one is interested in dynamics that occur over longer time scales, in which case a detailed simulation of each cardiac contraction becomes computationally burdensome. In this paper, we apply circuit-averaging techniques to a periodically driven, closed-loop, three-compartment recirculation model. The resultant cycle-averaged model is linear and time invariant, and greatly reduces the computational burden. It is also amenable to systematic order reduction methods that lead to further efficiencies. Despite its simplicity, the averaged model captures the dynamics relevant to the representation of a range of cardiovascular reflex mechanisms. c2004 Elsevier Ltd. All rights reserved.

Non-NASA Center↗
NASA NTRSDate not supplied

Parametric Mechanism Design Through Numerical Optimization and Physics Simulation

Design-Build-Test approaches for developing spaceflight hardware are prohibitively time and cost intensive and often lead to suboptimal mechanism designs. Approaches that couple machine learning and high-fidelity physics simulation could eliminate the need for hardware prototyping and dramatically accelerate the engineering design cycle, ultimately reducing cost. This work presents a modular NASA-developed toolchain to optimize hardware mechanisms in a virtual environment using numerical optimization and multi-body physics simulation. The toolchain enables multi-objective optimization, generates parametric CAD files that can be further post-processed by an end user, and can be expanded to optimize full systems and non-mechanical parameters such as feedback control variables. We demonstrate the toolchain through an independently verifiable design problem that optimizes wheel radius to achieve a desired linear velocity in a rigid-body physics environment when the wheel rotates at a constant angular speed, and then post-process the parametric CAD file of the optimal design generated by the tool before ultimately manufacturing it via 3D printing. We end with a discussion of how the toolchain can incorporate other analysis tools, including finite element analysis, computational fluid dynamics, and granular media simulations.

Optimization↗
NASA NTRSDate not supplied

An Optimization-Based Toolchain for Parametric Mechanism Design

Design-Build-Test approaches for developing spaceflight hardware are prohibitively time and cost intensive and often lead to suboptimal mechanism designs. Approaches that couple machine learning and high-fidelity physics simulation could eliminate the need for hardware prototyping and dramatically accelerate the engineering design cycle, ultimately reducing cost. This work presents a modular NASA-developed toolchain to optimize hardware mechanisms in a virtual environment using numerical optimization and multi-body physics simulation. The toolchain enables multi-objective optimization, generates parametric CAD files that can be further post-processed by an end user, and can be expanded to optimize full systems and non-mechanical parameters such as feedback control variables. We demonstrate the toolchain through an independently verifiable design problem that optimizes wheel radius to achieve a desired linear velocity in a rigid-body physics environment when the wheel rotates at a constant angular speed, and then post-process the parametric CAD file of the optimal design generated by the tool before ultimately manufacturing it via 3D printing. We end with a discussion of how the toolchain can incorporate other analysis tools, including finite element analysis, computational fluid dynamics, and granular media simulations.

Optimization↗
NASA NTRS1993

Research on Earth's rotation and the effect of atmospheric pressure on vertical deformation and sea level variability

The work done under NASA grant NAG5-485 included modelling the deformation of the earth caused by variations in atmospheric pressure. The amount of deformation near coasts is sensitive to the nature of the oceanic response to the pressure. The PSMSL (Permanent Service for Mean Sea Level) data suggest the response is inverted barometer at periods greater than a couple months. Green's functions were constructed to describe the perturbation of the geoid caused by atmospheric and oceanic loading and by the accompanying load-induced deformation. It was found that perturbation of up to 2 cm are possible. Ice mass balance data was used for continental glaciers to look at the glacial contributions to time-dependent changes in polar motion, the lod, the earth's gravitational field, the position of the earth's center-of-mass, and global sea level. It was found that there can be lateral, non-hydrostatic structure inside the fluid core caused by gravitational forcing from the mantle, from the inner core, or from topography at the core/mantle or inner core/outer core boundaries. The nutational and tidal response of a non-hydrostatic earth with a solid inner core was modeled. Monthly, global tide gauge data from PSMSL was used to look at the 18.6-year ocean tide, the 14-month pole tide, the oceanic response to pressure, the linear trend and inter-annual variability in the earth's gravity field, the global sea level rise, and the effects of post glacial rebound. The effects of mantle anelasticity on nutations, earth tides, and tidal variation in the lod was modeled. Results of this model can be used with Crustal Dynamics observations to look at the anelastic dissipation and dispersion at tidal periods. The effects of surface topography on various components of crustal deformation was also modeled, and numerical models were developed of post glacial rebound.

Wahr, John↗
NASA NTRS2016

Improving Fidelity of Launch Vehicle Liftoff Acoustic Simulations

Launch vehicles experience high acoustic loads during ignition and liftoff affected by the interaction of rocket plume generated acoustic waves with launch pad structures. Application of highly parallelized Computational Fluid Dynamics (CFD) analysis tools optimized for application on the NAS computer systems such as the Loci/CHEM program now enable simulation of time-accurate, turbulent, multi-species plume formation and interaction with launch pad geometry and capture the generation of acoustic noise at the source regions in the plume shear layers and impingement regions. These CFD solvers are robust in capturing the acoustic fluctuations, but they are too dissipative to accurately resolve the propagation of the acoustic waves throughout the launch environment domain along the vehicle. A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed to improve such liftoff acoustic environment predictions. The framework combines the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin (DG) solver, Loci/THRUST, developed in the same computational framework. Loci/THRUST employs a low dissipation, high-order, unstructured DG method to accurately propagate acoustic waves away from the source regions across large distances. The DG solver is currently capable of solving up to 4th order solutions for non-linear, conservative acoustic field propagation. Higher order boundary conditions are implemented to accurately model the reflection and refraction of acoustic waves on launch pad components. The DG solver accepts generalized unstructured meshes, enabling efficient application of common mesh generation tools for CHEM and THRUST simulations. The DG solution is coupled with the CFD solution at interface boundaries placed near the CFD acoustic source regions. Both simulations are executed simultaneously with coordinated boundary condition data exchange.

Liever, Peter↗
NASA NTRS1976

The design of digital-adaptive controllers for VTOL aircraft

Design procedures for VTOL automatic control systems have been developed and are presented. Using linear-optimal estimation and control techniques as a starting point, digital-adaptive control laws have been designed for the VALT Research Aircraft, a tandem-rotor helicopter which is equipped for fully automatic flight in terminal area operations. These control laws are designed to interface with velocity-command and attitude-command guidance logic, which could be used in short-haul VTOL operations. Developments reported here include new algorithms for designing non-zero-set-point digital regulators, design procedures for rate-limited systems, and algorithms for dynamic control trim setting.

Stengel, R. F.↗
NASA NTRS2002

Argus: a new instrument for the measurement of the stratospheric dynamical tracers, N2O and CH4

We describe here a new instrument for the simultaneous, in situ measurement of the stratospheric tracer molecules, nitrous oxide (N2O) and methane (CH4). Argus is unique in its small size making it well suited for limited payload atmospheric research platforms. Argus employs second harmonic spectroscopy using tunable lead-salt diode lasers emitting in the mid-infrared. We first explain the Argus design philosophy followed by detailed descriptions of the instrument's optical, mechanical, and thermal sub-systems. Argus employs an in-flight calibration system providing real time calibrations and tightly constrained uncertainty estimates of the returned data. Data analysis is carried out using non-linear least-squares model fits to the acquired second harmonic spectra. A sampling of Argus data acquired on a recent stratospheric research campaign in the Arctic winter is presented.

Methane/analysis↗
NASA NTRS1999

Preview-Based Stable-Inversion for Output Tracking

Stable Inversion techniques can be used to achieve high-accuracy output tracking. However, for nonminimum phase systems, the inverse is non-causal - hence the inverse has to be pre-computed using a pre-specified desired-output trajectory. This requirement for pre-specification of the desired output restricts the use of inversion-based approaches to trajectory planning problems (for nonminimum phase systems). In the present article, it is shown that preview information of the desired output can be used to achieve online inversion-based output tracking of linear systems. The amount of preview-time needed is quantified in terms of the tracking error and the internal dynamics of the system (zeros of the system). The methodology is applied to the online output tracking of a flexible structure and experimental results are presented.

Zou, Qing-Ze↗
NASA NTRS1983

Solar maximum mission: Ground support programs at the Harvard Radio Astronomy Station

Observations of the spectral characteristics of solar radio bursts were made with new dynamic spectrum analyzers of high sensitivity and high reliability, over the frequency range 25-580 MHz. The observations also covered the maximum period of the current solar cycle and the period of international cooperative programs designated as the Solar Maximum Year. Radio data on shock waves generated by solar flares were combined with optical data on coronal transients, taken with equipment on the SMM and other satellites, and then incorporated into computer models for the outward passage of fast-mode MHD shocks through the solar corona. The MHD models are non-linear, time-dependent and for the most recent models, quasi-three-dimensional. They examine the global response of the corona for different types of input pulses (thermal, magnetic, etc.) and for different magnetic topologies (for example, open and closed fields). Data on coronal shocks and high-velocity material ejected from solar flares have been interpreted in terms of a model consisting of three main velocity regimes.

Maxwell, A.↗
NASA NTRS2000

Spatiotemporal processing of linear acceleration: primary afferent and central vestibular neuron responses

Spatiotemporal convergence and two-dimensional (2-D) neural tuning have been proposed as a major neural mechanism in the signal processing of linear acceleration. To examine this hypothesis, we studied the firing properties of primary otolith afferents and central otolith neurons that respond exclusively to horizontal linear accelerations of the head (0.16-10 Hz) in alert rhesus monkeys. Unlike primary afferents, the majority of central otolith neurons exhibited 2-D spatial tuning to linear acceleration. As a result, central otolith dynamics vary as a function of movement direction. During movement along the maximum sensitivity direction, the dynamics of all central otolith neurons differed significantly from those observed for the primary afferent population. Specifically at low frequencies (</=0.5 Hz), the firing rate of the majority of central otolith neurons peaked in phase with linear velocity, in contrast to primary afferents that peaked in phase with linear acceleration. At least three different groups of central response dynamics were described according to the properties observed for motion along the maximum sensitivity direction. "High-pass" neurons exhibited increasing gains and phase values as a function of frequency. "Flat" neurons were characterized by relatively flat gains and constant phase lags (approximately 20-55 degrees ). A few neurons ("low-pass") were characterized by decreasing gain and phase as a function of frequency. The response dynamics of central otolith neurons suggest that the approximately 90 degrees phase lags observed at low frequencies are not the result of a neural integration but rather the effect of nonminimum phase behavior, which could arise at least partly through spatiotemporal convergence. Neither afferent nor central otolith neurons discriminated between gravitational and inertial components of linear acceleration. Thus response sensitivity was indistinguishable during 0.5-Hz pitch oscillations and fore-aft movements. The fact that otolith-only central neurons with "high-pass" filter properties exhibit semicircular canal-like dynamics during head tilts might have important consequences for the conclusions of previous studies of sensory convergence and sensorimotor transformations in central vestibular neurons.

NASA Discipline Neuroscience↗
NASA NTRS2015

Regional Rates of Young US Forest Growth Estimated From Annual Landsat Disturbance History and IKONOS Stereo Imagery

Forests of the Contiguous United States (CONUS) have been found to be a large contributor to the global atmospheric carbon sink. The magnitude and nature of this sink is still uncertain and recent studies have sought to define the dynamics that control its strength and longevity. The Landsat series of satellites has been a vital resource to understand the long-term changes in land cover that can impact ecosystem function and terrestrial carbonstock. We combine annual Landsat forest disturbance history from 1985 to 2011 with single date IKONOS stereoimagery to estimate the change in young forest canopy height and above ground live dry biomass accumulation for selected sites in the CONUS. Our approach follows an approximately linear growth rate following clearing over short intervals and does not estimate the distinct non-linear growth rate over longer intervals.We produced canopy height models by differencing digital surface models estimated from IKONOS stereo pairs with national elevation data (NED). Correlations between height and biomass were established independently using airborne LiDAR, and then applied to the IKONOS-estimated canopy height models. Graphing current biomass against time since disturbance provided biomass accumulation rates. For 20 study sites distributed across five regions of the CONUS, 19 showed statistically significant recovery trends (p is less than 0.001) with canopy growth from 0.26 m yr-1to 0.73 m yr-1. Above ground live dry biomass (AGB) density accumulation ranged from 1.31 t/ha yr-1 to 12.47 t/ha yr-1. Mean forest AGB accumulationwas 6.31 t/ha yr-1 among all sites with significant growth trends. We evaluated the accuracy of our estimates by comparing to field estimated site index curves of growth, airborne LiDAR data, and independent model predictions of C accumulation. Growth estimates found with this approach are consistent with site index curves and total biomass estimates fall within the range of field estimates. This is aviable approach to estimate forest biomass accumulation in regions with clear-cut harvest disturbances.

Neigh, Christopher S. R.↗
NASA NTRS1985

Modal forced vibration analysis of aerodynamically excited turbosystems

Theoretical aspects of a new capability to determine the vibratory response of turbosystems subjected to aerodynamic excitation are presented. Turbosystems such as advanced turbopropellers with highly swept blades, and axial-flow compressors and turbines can be analyzed using this capability. The capability has been developed and implemented in the April 1984 release of the general purpose finite element program NASTRAN. The dynamic response problem is addressed in terms of the normal modal coordinates of these tuned rotating cyclic structures. Both rigid and flexible hubs/disks are considered. Coriolis and centripetal accelerations, as well as differential stiffness effects are included. Generally non-uniform steady inflow fields and uniform flow fields arbitrarily inclined at small angles with respect to the axis of rotation of the turbosystem are considered sources of aerodynamic excitation. The spatial non-uniformities are considered to be small deviations from a principally uniform inflow. Subsonic and supersonic relative inflows are addressed, with provision for linearly interpolating transonic airloads.

Elchuri, V.↗
NASA NTRS2004

Input/output properties of the lateral vestibular nucleus

This article is a review of work in three species, squirrel monkey, cat, and rat studying the inputs and outputs from the lateral vestibular nucleus (LVN). Different electrophysiological shock paradigms were used to determine the synaptic inputs derived from thick to thin diameter vestibular nerve afferents. Angular and linear mechanical stimulations were used to activate and study the combined and individual contribution of inner ear organs and neck afferents. The spatio-temporal properties of LVN neurons in the decerebrated rat were studied in response to dynamic acceleration inputs using sinusoidal linear translation in the horizontal head plane. Outputs were evaluated using antidromic identification techniques and identified LVN neurons were intracellularly injected with biocytin and their morphology studied.

Efferent Pathways/anatomy & histology/physiology↗
NASA NTRS1993

Active stabilization to prevent surge in centrifugal compression systems

This report documents an experimental and analytical study of the active stabilization of surge in a centrifugal engine. The aims of the research were to extend the operating range of a compressor as far as possible and to establish the theoretical framework for the active stabilization of surge from both an aerodynamic stability and a control theoretic perspective. In particular, much attention was paid to understanding the physical limitations of active stabilization and how they are influenced by control system design parameters. Previously developed linear models of actively stabilized compressors were extended to include such nonlinear phenomena as bounded actuation, bandwidth limits, and robustness criteria. This model was then used to systematically quantify the influence of sensor-actuator selection on system performance. Five different actuation schemes were considered along with four different sensors. Sensor-actuator choice was shown to have a profound effect on the performance of the stabilized compressor. The optimum choice was not unique, but rather shown to be a strong function of some of the non-dimensional parameters which characterize the compression system dynamics. Specifically, the utility of the concepts were shown to depend on the system compliance to inertia ratio ('B' parameter) and the local slope of the compressor speedline. In general, the most effective arrangements are ones in which the actuator is most closely coupled to the compressor, such as a close-coupled bleed valve inlet jet, rather than elsewhere in the flow train, such as a fuel flow modulator. The analytical model was used to explore the influence of control system bandwidth on control effectiveness. The relevant reference frequency was shown to be the compression system's Helmholtz frequency rather than the surge frequency. The analysis shows that control bandwidths of three to ten times the Helmholtz frequency are required for larger increases in the compressor flow range. This has important implications for implementation in gas turbine engines since the Helmholtz frequencies can be over 100 Hz, making actuator design extremely challenging.

Epstein, Alan H.↗