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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Aerodynamic characteristics of the standard dynamics model in coning motion at Mach 0.6

A wind tunnel test was conducted on the Standard Dynamics Model (a simplified generic fighter aircraft shape) undergoing coning motion at Mach 0.6. Six component force and moment data are presented for a range of angle of attack, sideslip, and coning rates. At the relatively low non-dimensional coning rate employed (omega b/2V less than or equal to 0.04), the lateral aerodynamic characteristics generally show a linear variation with coning rate.

Jermey, C.↗
NASA NTRS2010

Shallow and Deep Latent Heating Modes Over Tropical Oceans Observed with TRMM PR Spectral Latent Heating Data

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of its formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the Tropics. Associated with this rainfall is a vast amount of heat, which is known as latent heat. It arises mainly due to the phase change of water vapor condensing into liquid droplets; three-fourths of the total heat energy available to the Earth's atmosphere comes from tropical rainfall. In addition, fresh water provided by tropical rainfall and its variability exerts a large impact upon the structure and motions of the upper ocean layer. Three-dimensional distributions of latent heating estimated from Tropical Rainfall Measuring Mission Precipitation Radar (TRMM PR)utilizing the Spectral Latent Heating (SLH) algorithm are analyzed. Mass-weighted and vertically integrated latent heating averaged over the tropical oceans is estimated as approx.72.6 J/s (approx.2.51 mm/day), and that over tropical land is approx.73.7 J/s (approx.2.55 mm/day), for 30degN-30degS. It is shown that non-drizzle precipitation over tropical and subtropical oceans consists of two dominant modes of rainfall systems, deep systems and congestus. A rough estimate of shallow mode contribution against the total heating is about 46.7 % for the average tropical oceans, which is substantially larger than 23.7 % over tropical land. While cumulus congestus heating linearly correlates with the SST, deep mode is dynamically bounded by large-scale subsidence. It is notable that substantial amount of rain, as large as 2.38 mm day-1 in average, is brought from congestus clouds under the large-scale subsiding circulation. It is also notable that even in the region with SST warmer than 28 oC, large-scale subsidence effectively suppresses the deep convection, remaining the heating by congestus clouds. Our results support that the entrainment of mid-to-lower-tropospheric dry air, which accompanies the large-scale subsidence is the major factor suppressing the deep convection. Therefore, representation of the realistic entrainment is very important for proper reproduction of precipitation distribution and resultant large-scale circulation.

Takayabu, Yukari N.↗
NASA NTRS1998

Efficient and Optimal Attitude Determination Using Recursive Global Positioning System Signal Operations

In this paper, a new and efficient algorithm is developed for attitude determination from Global Positioning System signals. The new algorithm is derived from a generalized nonlinear predictive filter for nonlinear systems. This uses a one time-step ahead approach to propagate a simple kinematics model for attitude determination. The advantages of the new algorithm over previously developed methods include: it provides optimal attitudes even for coplanar baseline configurations; it guarantees convergence even for poor initial conditions; it is a non-iterative algorithm; and it is computationally efficient. These advantages clearly make the new algorithm well suited to on-board applications. The performance of the new algorithm is tested on a dynamic hardware simulator. Results indicate that the new algorithm accurately estimates the attitude of a moving vehicle, and provides robust attitude estimates even when other methods, such as a linearized least-squares approach, fail due to poor initial starting conditions.

Crassidis, John L.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRS2018

Solution of Ordinary Differential Equations in Gradient-Based Multidisciplinary Design Optimization

A gradient-based approach to multidisciplinary design optimization enables efficient scalability to large numbers of design variables. However, the need for derivatives causes difficulties when integrating ordinary differential equations (ODEs) in models. To simplify this, we propose the use of the general linear methods framework, which unifies all Runge-Kutta and linear multistep methods. This approach enables rapid implementation of integration methods without the need to differentiate each one, even in a gradient-based optimization context. We also develop a new parallel time integration algorithm that enables vectorization across time steps. We present a set of benchmarking results using a stiff ODE, a non-stiff nonlinear ODE, and an orbital dynamics ODE, and compare integration methods. In a modular gradient-based multidisciplinary design optimization context, we find that the new parallel time integration algorithm with high-order implicit methods, especially Gauss-Legendre collocation, is the best choice for a broad range of problems.

Hwang, John T.↗
NASA NTRSDate not supplied

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flowfields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperaturelimited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗
NASA NTRS2000

Ordering Unstructured Meshes for Sparse Matrix Computations on Leading Parallel Systems

The ability of computers to solve hitherto intractable problems and simulate complex processes using mathematical models makes them an indispensable part of modern science and engineering. Computer simulations of large-scale realistic applications usually require solving a set of non-linear partial differential equations (PDES) over a finite region. For example, one thrust area in the DOE Grand Challenge projects is to design future accelerators such as the SpaHation Neutron Source (SNS). Our colleagues at SLAC need to model complex RFQ cavities with large aspect ratios. Unstructured grids are currently used to resolve the small features in a large computational domain; dynamic mesh adaptation will be added in the future for additional efficiency. The PDEs for electromagnetics are discretized by the FEM method, which leads to a generalized eigenvalue problem Kx = AMx, where K and M are the stiffness and mass matrices, and are very sparse. In a typical cavity model, the number of degrees of freedom is about one million. For such large eigenproblems, direct solution techniques quickly reach the memory limits. Instead, the most widely-used methods are Krylov subspace methods, such as Lanczos or Jacobi-Davidson. In all the Krylov-based algorithms, sparse matrix-vector multiplication (SPMV) must be performed repeatedly. Therefore, the efficiency of SPMV usually determines the eigensolver speed. SPMV is also one of the most heavily used kernels in large-scale numerical simulations.

Oliker, Leonid↗
NASA NTRSDate not supplied

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flow fields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperature limited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗
NASA NTRS1997

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations (DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those of a spatially evolving jet, a temporal jet problem was solved, using periodicity in the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible appli(,a- tion to active noise suppression. In addition, the data generated can be used to compute, various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved DNS of a turbulent jet flow is produced.

Rao, Ram Mohan↗
NASA NTRS1997

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations(DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those, of a spatially evolving jet, a temporal jet problem was solved, using periodicity ill the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible application to active noise suppression. In addition, the data generated can be used to compute various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved Direct Numerical Simulation (DNS) of a turbulent jet flow is produced.

Rao, Ram Mohan↗
NASA NTRS2004

Proceedings from the 2nd International Symposium on Formation Flying Missions and Technologies

Topics discussed include: The Stellar Imager (SI) "Vision Mission"; First Formation Flying Demonstration Mission Including on Flight Nulling; Formation Flying X-ray Telescope in L2 Orbit; SPECS: The Kilometer-baseline Far-IR Interferometer in NASA's Space Science Roadmap Presentation; A Tight Formation for Along-track SAR Interferometry; Realization of the Solar Power Satellite using the Formation Flying Solar Reflector; SIMBOL-X : Formation Flying for High-Energy Astrophysics; High Precision Optical Metrology for DARWIN; Close Formation Flight of Micro-Satellites for SAR Interferometry; Station-Keeping Requirements for Astronomical Imaging with Constellations of Free-Flying Collectors; Closed-Loop Control of Formation Flying Satellites; Formation Control for the MAXIM Mission; Precision Formation Keeping at L2 Using the Autonomous Formation Flying Sensor; Robust Control of Multiple Spacecraft Formation Flying; Virtual Rigid Body (VRB) Satellite Formation Control: Stable Mode-Switching and Cross-Coupling; Electromagnetic Formation Flight (EMFF) System Design, Mission Capabilities, and Testbed Development; Navigation Algorithms for Formation Flying Missions; Use of Formation Flying Small Satellites Incorporating OISL's in a Tandem Cluster Mission; Semimajor Axis Estimation Strategies; Relative Attitude Determination of Earth Orbiting Formations Using GPS Receivers; Analysis of Formation Flying in Eccentric Orbits Using Linearized Equations of Relative Motion; Conservative Analytical Collision Probabilities for Orbital Formation Flying; Equations of Motion and Stability of Two Spacecraft in Formation at the Earth/Moon Triangular Libration Points; Formations Near the Libration Points: Design Strategies Using Natural and Non-Natural Ares; An Overview of the Formation and Attitude Control System for the Terrestrial Planet Finder Formation Flying Interferometer; GVE-Based Dynamics and Control for Formation Flying Spacecraft; GNC System Design for a New Concept of X-Ray Distributed Telescope; GNC System for the Deployment and Fine Control of the DARWIN Free-Flying Interferometer; Formation Algorithm and Simulation Testbed; and PLATFORM: A Formation Flying, RvD and Robotic Validation Test-bench.

Source record↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗
NASA NTRSDate not supplied

Changes in Oculomotor Behavior and PVT Reaction Time During One Night of Sleep Deprivation

INTRODUCTION: The Psychomotor Vigilance Test (PVT) is a widely used objective measure of sustained attention, alertness, and fatigue. Previous studies consistently show that prolonged wakefulness and disrupted sleep patterns lead to increased lapses of attention and slower reaction times. These findings have critical implications for various professions, including healthcare, aviation, transportation safety, and spaceflight. In this study, we employed linear mixed models (LMM) to investigate hourly-measured 5-minute PVT reaction times during one night of total sleep deprivation. Our goal was to characterize the contribution of homeostatic sleep pressure (time awake) and circadian phase (salivary melatonin levels) on the dynamics of PVT reaction time. MATERIALS & METHODS: Data from twelve human participants were used in the analysis. Participants were healthy non-smokers, aged 18 and 40, with normal sleep habits defined as Pittsburg Sleep Quality Index scores < 5, and Morningness – Eveningness Questionnaire scores > 42 and < 58, respectively. Participants followed a constant-routine protocol and arrived at the NASA sleep laboratory approximately 1 to 2 hours after their habitual waking time. Tympanic temperature was measured every 30 min to provide a real-time estimate of circadian phase. Participants were required to stay awake throughout the laboratory experiment, which continued until their tympanic temperature had returned to baseline levels. This approach ensured that we could capture each participant's circadian trough and recovery, typically occurring between 24 and 26 hours after waking. The PVT and saliva melatonin were measured every hour, starting approximately from 3 hours after awakening. More information about the measurement protocol can be found in the publication by Stone et al. We used PVT reaction times between 100 and 10000 ms to calculate Reciprocal Reaction Times (RRTs). LMM models were employed to assess the contribution of time awake and saliva melatonin levels on changes in RRTs. Two models were created: the first model treated time awake as a fixed effect and the second model considered both the salivary melatonin level and time awake as fixed effects. Participants were included as a random effect in both models. We compared the models using an Analysis of Variance (ANOVA). The models were computed using the 'lmer' function from the 'lme4' package in R, and LMM p-values were determined using Satterthwaite's method. RESULTS: We found that time awake had a statistically significant effect on PVT RRTs during one night of sleep deprivation [χ2 = 2132.4, df = 1, p < 0.001]. In addition, the model including both saliva melatonin levels and time awake was statistically significantly better than one that uses just time awake [χ2 = 17.4, df = 1, p < 0.001]. DISCUSSION & CONCLUSION: As expected, these preliminary results suggest that the time awake and circadian phase each contribute to performance reductions as measured by the PVT. Future analyses will explore the contributions of other factors to the model.

sleep deprivation↗
NASA NTRS1975

A study of attitude control concepts for precision-pointing non-rigid spacecraft

Attitude control concepts for use onboard structurally nonrigid spacecraft that must be pointed with great precision are examined. The task of determining the eigenproperties of a system of linear time-invariant equations (in terms of hybrid coordinates) representing the attitude motion of a flexible spacecraft is discussed. Literal characteristics are developed for the associated eigenvalues and eigenvectors of the system. A method is presented for determining the poles and zeros of the transfer function describing the attitude dynamics of a flexible spacecraft characterized by hybrid coordinate equations. Alterations are made to linear regulator and observer theory to accommodate modeling errors. The results show that a model error vector, which evolves from an error system, can be added to a reduced system model, estimated by an observer, and used by the control law to render the system less sensitive to uncertain magnitudes and phase relations of truncated modes and external disturbance effects. A hybrid coordinate formulation using the provided assumed mode shapes, rather than incorporating the usual finite element approach is provided.

Likins, P. W.↗
NASA NTRS1997

Piezoelectric Power Requirements for Active Vibration Control

This paper presents a method for predicting the power consumption of piezoelectric actuators utilized for active vibration control. Analytical developments and experimental tests show that the maximum power required to control a structure using surface-bonded piezoelectric actuators is independent of the dynamics between the piezoelectric actuator and the host structure. The results demonstrate that for a perfectly-controlled system, the power consumption is a function of the quantity and type of piezoelectric actuators and the voltage and frequency of the control law output signal. Furthermore, as control effectiveness decreases, the power consumption of the piezoelectric actuators decreases. In addition, experimental results revealed a non-linear behavior in the material properties of piezoelectric actuators. The material non- linearity displayed a significant increase in capacitance with an increase in excitation voltage. Tests show that if the non-linearity of the capacitance was accounted for, a conservative estimate of the power can easily be determined.

Brennan, Matthew C.↗
NASA NTRS1979

Surprisal analysis of rotational-translational energy transfer - Non-linear versus linear rotors

Surprisal versus energy gap analyses of state-to-state cross sections are presented for a number of linear rigid rotors excited by collisions with atoms for H2-H, H2-He, HCl-He, HCl-Ar, CO-He, CS-H2 (j=0) OCS-H2 (j=0) and HN2(+)-He, where (j=0) indicates that the hydrogen molecule was constrained to remain in its lowest level. Different systems exhibit wide variations in the slope of the surprisal plot and in certain cases, enough to indicate that the energy gap may not be the static dynamical constraint. Similar analyses are presented for nonlinear rotors excited by atoms for H2CO-He and H2O-He. For these, the data show a great deal of scatter, indicating that the reduced energy gap is probably not the appropriate independent variable.

Green, S.↗
NASA NTRS2003

Cranial diameter pulsations measured by non-invasive ultrasound decrease with tilt

INTRODUCTION: Intracranial pressure (ICP) may play a significant role in physiological responses to microgravity by contributing to the nausea associated with microgravity exposure. However, effects of altered gravity on ICP in astronauts have not been investigated, primarily due to the invasiveness of currently available techniques. We have developed an ultrasonic device that monitors changes in cranial diameter pulsation non-invasively so that we can evaluate ICP dynamics in astronauts during spaceflight. This study was designed to demonstrate the feasibility of our ultrasound technique under the physiological condition in which ICP dynamics are changed due to altered gravitational force. METHODS: Six healthy volunteers were placed at 60 degrees head-up, 30 degrees headup, supine, and 15 degrees head-down positions for 3 min at each angle. We measured arterial blood pressure (ABP) with a finger pressure cuff, and cranial diameter pulsation with a pulsed phase lock loop device (PPLL). RESULTS: Analysis of covariance demonstrated that amplitudes of cranial diameter pulsations were significantly altered with the angle of tilt (p < 0.001). The 95% confidence interval for linear regression coefficients of the cranial diameter pulsation amplitudes with tilt angle was 0.862 to 0.968. However, ABP amplitudes did not show this relationship. DISCUSSION: Our noninvasive ultrasonic technique reveals that the amplitude of cranial diameter pulsation decreases as a function of tilt angle, suggesting that ICP pulsation follows the same relationship. It is demonstrated that the PPLL device has a sufficient sensitivity to detect changes non-invasively in ICP pulsation caused by altered gravity.

NASA Program Biomedical Research and Countermeasur↗
NASA NTRS2012

Particle Filtering for Model-Based Anomaly Detection in Sensor Networks

A novel technique has been developed for anomaly detection of rocket engine test stand (RETS) data. The objective was to develop a system that postprocesses a csv file containing the sensor readings and activities (time-series) from a rocket engine test, and detects any anomalies that might have occurred during the test. The output consists of the names of the sensors that show anomalous behavior, and the start and end time of each anomaly. In order to reduce the involvement of domain experts significantly, several data-driven approaches have been proposed where models are automatically acquired from the data, thus bypassing the cost and effort of building system models. Many supervised learning methods can efficiently learn operational and fault models, given large amounts of both nominal and fault data. However, for domains such as RETS data, the amount of anomalous data that is actually available is relatively small, making most supervised learning methods rather ineffective, and in general met with limited success in anomaly detection. The fundamental problem with existing approaches is that they assume that the data are iid, i.e., independent and identically distributed, which is violated in typical RETS data. None of these techniques naturally exploit the temporal information inherent in time series data from the sensor networks. There are correlations among the sensor readings, not only at the same time, but also across time. However, these approaches have not explicitly identified and exploited such correlations. Given these limitations of model-free methods, there has been renewed interest in model-based methods, specifically graphical methods that explicitly reason temporally. The Gaussian Mixture Model (GMM) in a Linear Dynamic System approach assumes that the multi-dimensional test data is a mixture of multi-variate Gaussians, and fits a given number of Gaussian clusters with the help of the wellknown Expectation Maximization (EM) algorithm. The parameters thus learned are used for calculating the joint distribution of the observations. However, this GMM assumption is essentially an approximation and signals the potential viability of non-parametric density estimators. This is the key idea underlying the new approach.

Solano, Wanda↗