Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time History”

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 397 records · Page 22

Synthesis of Optimal Nonlinear Feedback Laws for Dynamic Systems Using Neural Networks

Open-loop solutions of dynamical optimization problems can be numerically computed usingexisting software packages. The computed time histories of the state and control variables, formultiple sets of end conditions can then be used to train a neural network to 'recognize' the optimal,nonlinear feedback relation between the states and controls of the system. The 'learned' network canthen be used to output an approximate optimal control given a full set (or a partial set) of measuredsystem states. With simple neural networks, we have successfully demonstrated the efficacy of theproposed approach using a minimum-time orbit injection problem. The usefulness and limitations ofthis novel approach on real-life optimal guidance and control problems, with many state and control variables as well as path inequality constraints, remain to be seen.

Orbit injection problem↗

Time-Accurate Unsteady Pressure Loads Simulated for the Space Launch System at Wind Tunnel Conditions

A transonic flow field about a Space Launch System (SLS) configuration was simulated with the Fully Unstructured Three-Dimensional (FUN3D) computational fluid dynamics (CFD) code at wind tunnel conditions. Unsteady, time-accurate computations were performed using second-order Delayed Detached Eddy Simulation (DDES) for up to 1.5 physical seconds. The surface pressure time history was collected at 619 locations, 169 of which matched locations on a 2.5 percent wind tunnel model that was tested in the 11 ft. x 11 ft. test section of the NASA Ames Research Center's Unitary Plan Wind Tunnel. Comparisons between computation and experiment showed that the peak surface pressure RMS level occurs behind the forward attach hardware, and good agreement for frequency and power was obtained in this region. Computational domain, grid resolution, and time step sensitivity studies were performed. These included an investigation of pseudo-time sub-iteration convergence. Using these sensitivity studies and experimental data comparisons, a set of best practices to date have been established for FUN3D simulations for SLS launch vehicle analysis. To the author's knowledge, this is the first time DDES has been used in a systematic approach and establish simulation time needed, to analyze unsteady pressure loads on a space launch vehicle such as the NASA SLS.

Alter, Stephen J.↗

Calculating the High-Latitude Ionospheric Electrodynamics Using A Machine Learning-Based Field-Aligned Current Model

We introduce a new framework called Machine Learning (ML) based Auroral Ionospheric electrodynamics Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents of Kunduri et al. (2020, https://doi.org/10.1029/2020JA027908), the FAC-derived auroral conductance model of Robinson et al. (2020, https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen and Brekke (1993, https://doi.org/10.1029/92gl02109). The ML-AIM inputs are 60-min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pedersen/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity interval on 14 May 2013 and a geomagnetic storm on 7–8 September 2017. ML-AIM produces physically accurate ionospheric potential patterns such as the two-cell convection pattern and the enhancement of electric potentials during active times. The cross polar cap potentials (ΦPC) from ML-AIM, the Weimer (2005, https://doi.org/10.1029/2004ja010884) model, and the Super Dual Auroral Radar Network (SuperDARN) data-assimilated potentials, are compared to the ones from 3204 polar crossings of the Defense Meteorological Satellite Program F17 satellite, showing better performance of ML-AIM than others. ML-AIM is unique and innovative because it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while the other traditional empirical models like Weimer (2005, https://doi.org/10.1029/2004ja010884) designed to provide a quasi-static ionospheric condition under quasi-steady solar wind/IMF conditions. Plans are underway to improve ML-AIM performance by including a fully ML network of models of aurora precipitation and ionospheric conductance, targeting its characterization of geomagnetically active times.

auroral electrodynamics↗

Flared landing approach flying qualities. Volume 1: Experiment design and analysis

An inflight research study was conducted utilizing the USAF Total Inflight Simulator (TIFS) to investigate longitudinal flying qualities for the flared landing approach phase of flight. The purpose of the experiment was to generate a consistent set of data for: (1) determining what kind of commanded response the pilot prefers in order to flare and land an airplane with precision, and (2) refining a time history criterion that took into account all the necessary variables and their characteristics that would accurately predict flying qualities. The result of the first part provides guidelines to the flight control system designer, using MIL-F-8785-(C) as a guide, that yield the dynamic behavior pilots perfer in flared landings. The results of the second part provides the flying qualities engineer with a newly derived flying qualities predictive tool which appears to be highly accurate. This time domain predictive flying qualities criterion was applied to the flight data as well as six previous flying qualities studies, and the results indicate that the criterion predicted the flying qualities level 81% of the time and the Cooper-Harper pilot rating, within + or - 1, 60% of the time.

Weingarten, Norman C.↗

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

The triple Langmuir probe is an electrostatic probe in which three probe tips collect current when inserted into a plasma. The triple probe differs from a simple single Langmuir probe in the nature of the voltage applied to the probe tips. In the single probe, a swept voltage is applied to the probe tip to acquire a waveform showing the collected current as a function of applied voltage (I-V curve). In a triple probe three probe tips are electrically coupled to each other with constant voltages applied between each of the tips. The voltages are selected such that they would represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Under the assumption of a Maxwellian plasma, one can determine the time-varying plasma temperature T(sub e)(t) and number density n(sub e)(t) from the applied voltage levels and the time-histories of the collected currents. In the present paper we examine the theory of triple probe operation, specifically focusing on the assumption of a Maxwellian plasma. Triple probe measurements have been widely employed for a number of pulsed and timevarying plasmas, including pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. While the equilibrium assumption may be justified for some applications, it is unlikely that it is fully justifiable for all pulsed and time-varying plasmas or for all times during the pulse of a plasma device. To examine a simple non-equilibrium plasma case, we return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. We compare the distribution function that was assumed a priori with the distribution function one would infer when applying standard triple probe theory to analyze the collected currents. For the assumed class of non-Maxwellian distribution functions this serves to illustrate the effect a non-Maxwellian plasma would have on results interpreted using the equilibrium triple probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗

Ion sputtering from satellite surfaces

The manifestations of the induced surface currents as observed by the magnetic ion mass spectrometer (MIMS) and retarding potential analyzer (RPA) instruments on AE-C and AE-D are described. Data from these instruments are used to infer the source of some of the different secondary ions observed, as well as their time history and emission energy. It is noted that the emitted particles have a characteristic energy of approximately 1 eV in the spacecraft frame of reference and lack the normal ram energy possessed by the ambient ions. Alkali ions Na(+) and K(+), in addition to nonionospheric NO(+) and O2(+), are observed at low altitudes. With the exception of NO(+), all these are thought to be released by impacts of neutral N2 or O2 with the surfaces. These ion currents decrease with a time constant of about six weeks for the alkali and four weeks for the O2 ions.

Hanson, W. B.↗

Prevalent properties of gamma-ray burst variability

With the aim of isolating universal characteristics of gamma-ray burst variability, we compare time histories for 20 bright bursts detected by the Burst and Transient Source Experiment (Fishman et al., 1992). Using an autocorrelation function method, we find that the durations of individual structures within a burst, as well as the burst as a whole, become shorter with increasing energy in most events. We introduce the skewness function, a measure of temporal asymmetry. We find that most bursts exhibit a net skewness, in the sense that the intensity rises more rapidly than it falls. Over short time scales, bursts exhibit no preferred asymmetry. Taken together, these properties suggest that the overall time structure of a burst is due to an explosive phenomenon in which the evolution is initially energetic and rapid and is later slower during a 'cooling' period. We cannot rule out the possibility that short time-scale variability is caused by radiation beams sweeping past the observer.

Link, Bennett↗

Applications of Power Spectral Analysis Methods to Maneuver Loads Obtained on Jet Fighter Airplanes During Service Operations

Power spectral densities of normal load factor have been obtained for two service operational training flights of a Republic F-84G airplane and three service operational training flights of a North American F-86A airplane in order to indicate the load-factor frequency content and possible uses of power spectral methods in analyzing maneuver load data. It was determined that the maneuvering load-factor time histories appeared to be described by a truncated normal distribution. The power spectral densities obtained were relatively level at frequencies below 0.03 cycle per second and varied inversely with approximately the cube of the frequency at the higher frequencies. In general, the frequency content was very low above 0.2 cycle per second. The load-factor peak distributions were estimated fairly well from the spectrum analysis. In addition, peak load data obtained during service operations of fighter-type airplanes with flight time totaling about 24,000 hours were examined and appeared to agree reasonably well with the type of equations obtained from spectrum peak-load distributions.

Mayer, John P.↗

An Overview of the Patch Integral Method (PIM), a New Heat Transfer Analysis Tool for Hypersonic Wind Tunnel Facilities at NASA Langley

NASA Langley’s hypersonic wind tunnels are heavily leveraged for planetary missions. The data collection method in these tunnels is thermography, and surface temperature measurements of the model surface are collected and reduced to produce surface heating data, as seen in Fig 1. However, during model injection, no temperature data are collected, and thus conventional, integral heat transfer methods cannot be used to solve for surface heating. A method was developed in the 1990’s to reduce this data despite the data gap, known as the step approximation method. The method assumes that the film coefficient behaves as a step function, the model is semi-infinite, and thermal properties are constant. With these simplifying assumptions, a Laplace transform can be performed to result in an equation that takes the initial temperature of the model and a temperature at some point in time to back out the film coefficient at that time. This is the method that is used in the current thermographic data reduction software, IHEAT. While computationally light-weight, the step approximation has several issues associated with it. The time-history of temperature is not accounted for, which is vital as heat transfer is an integral process. Additionally, the required semi-infinite assumption is unnecessary and might be violated during runtime. Thermal variation of material properties can have a sizeable impact on heating results and are not modeled by the method. This method also takes multiple seconds to “collapse” to a steady state, which is undesirable from both a facility and data reduction standpoint. The method is also very sensitive to the “effective time” approximation, an approximation of when heating instantaneously starts (which is a nonphysical simplification), and a small variation in this value can result in an error in heating results.

J. S. Cheatwood↗

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗

Aerodynamic control of NASP-type vehicles through vortex manipulation. Volume 3: Wing rock experiments

Free-to-roll tests were conducted in water and wind tunnels in an effort to investigate the mechanisms of wing rock on a NASP-type vehicle. The configuration tested consisted of a highly-slender forebody and a 78 deg swept delta wing. In the water tunnel test, extensive flow visualization was performed and roll angle histories were obtained. In the wind tunnel test, the roll angle, forces and moments, and limited forebody and wing surface pressures were measured during the wing rock motion. A limit cycle oscillation was observed for angles of attack between 22 deg and 30 deg. In general, the experiments confirmed that the main flow phenomena responsible for the wing-body-tail wing rock are the interactions between the forebody and the wing vortices. The variation of roll acceleration (determined from the second derivative of the roll angle time history) with roll angle clearly slowed the energy balance necessary to sustain the limit cycle oscillation. Different means of suppressing wing rock by controlling the forebody vortices using small blowing jets were also explored. Steady blowing was found to be capable of suppressing wing rock, but significant vortex asymmetrices are created, causing the model to stop at a non-zero roll angle. On the other hand, alternating pulsed blowing on the left and right sides of the fore body was demonstrated to be a potentially effective means of suppressing wing rock and eliminating large asymmetric moments at high angles of attack.

Suarez, Carlos J.↗

Kinematic arguments against single relativistic shell models for GRBs

Two main types of models have been suggested to explain the long durations and multiple peaks of Gamma Ray Bursts (GRBs). In one, there is a very quick release of energy at a central site resulting in a single relativistic shell that produces peaks in the time history through its interactions with the ambient material. In the other, the central site sporadically releases energy over hundreds of seconds forming a peak with each burst of energy. The authors show that the average envelope of emission and the presence of gaps in GRBs are inconsistent with a single relativistic shell. They estimate that the maximum fraction of a single shell that can produce gamma-rays in a GRB with multiple peaks is 10(exp (minus)3), implying that single relativistic shells require 10(exp 3) times more energy than previously thought. They conclude that either the central site of a GRB must produce (approx)10(exp 51) erg/s(exp (minus)1) for hundreds of seconds, or the relativistic shell must have structure on a scales the order of (radical)(epsilon)(Gamma)(exp (minus)1), where (Gamma) is the bulk Lorentz factor ((approximately)10(exp 2) to 10(exp 3)) and (epsilon) is the efficiency.

Fenimore, Edward E.↗

Phase Function Determination in Support of Orbital Debris Size Estimation

To recover the size of a space debris object from photometric measurements, it is necessary to determine its albedo and basic shape: if the albedo is known, the reflective area can be calculated; and if the shape is known, the shape and area taken together can be used to estimate a characteristic dimension. Albedo is typically determined by inferring the object s material type from filter photometry or spectroscopy and is not the subject of the present study. Object shape, on the other hand, can be revealed from a time-history of the object s brightness response. The most data-rich presentation is a continuous light-curve that records the object s brightness for an entire sensor pass, which could last for tens of minutes to several hours: from this one can see both short-term periodic behavior as well as brightness variations with phase angle. Light-curve interpretation, however, is more art than science and does not lend itself easily to automation; and the collection method, which requires single-object telescope dedication for long periods of time, is not well suited to debris survey conditions. So one is led to investigate how easily an object s brightness phase function, which can be constructed from the more survey-friendly point photometry, can be used to recover object shape. Such a recovery is usually attempted by comparing a phase-function curve constructed from an object s empirical brightness measurements to analytically-derived curves for basic shapes or shape combinations. There are two ways to accomplish this: a simple averaged brightness-versus phase curve assembled from the empirical data, or a more elaborate approach in which one is essentially calculating a brightness PDF for each phase angle bin (a technique explored in unpublished AFRL/RV research and in Ojakangas 2011); in each case the empirical curve is compared to analytical results for shapes of interest. The latter technique promises more discrimination power but requires more data; the former can be assembled in its essentials from fewer measurements but will be less definitive in its assignments. The goal of the present study is to evaluate both techniques under debris survey conditions to determine their relative performance and, additionally, to learn precisely how a survey should be conducted in order to maximize their performance. Because the distendedness of objects has more of an effect than their precise shape in calculating a characteristic dimension, one is interested in the techniques discrimination ability to distinguish between an elongated rectangular prism and a short rectangular prism or cube, or an elongated cylinder from a squat cylinder or sphere. Sensitivity studies using simulated data will be conducted to determine discrimination power for both techniques as a function of amount of data collected and range (and specific region) of phase angles sampled. Empirical GEODSS photometry data for distended objects (dead payloads with solar panels, rocket bodies) and compact objects (cubesats, calibration spheres, squat payloads) will also be used to test this discrimination ability. The result will be a recommended technique and data collection paradigm for debris surveys in order to maximize this type of discrimination.

Hejduk, M. D.↗

Optimal filtering in the presence of unmodeled time correlated driving disturbances.

In many realistic data filtering problems, the cross correlation of the state estimation error and the state forcing function is unknown due to the poor knowledge of the time history of the forcing function. In this paper, the conservative and minimal approximation to the cross correlation terms is presented. It requires only the knowledge of the estimation error covariance and the forcing function covariance, with the choice of an associated free parameter left to the user. If the estimation error covariance and/or the forcing function covariance are bounded from above but not known exactly, the cross correlation approximation using those upper bounds remains conservative. This cross correlation approximation leads to a conservative approximation to the estimation error covariance matrix differential equation between measurement times. The free parameter is determined as the analytic solution to an associated optimal control problem. The procedure is expanded to include discrete linear measurement incorporation.

Fraser, D. C.↗

Notes on Digital-Data Analysis

The presentation is an overview and discussion of concepts normally encountered during digital-signal processing of sampled time histories. It is not a comprehensive literature survey, but simply a practical guide for best usage. Sections of the talk cover: digital sampling and representation, the Nyquist-Shannon sampling theorem, aliasing or frequency folding, the power spectrum, windowing and spectral spreading, scalloping, and best practices. The Advanced Air Transport Technology Project supports this effort.

Aeroacoustics↗

Flared landing approach flying qualities. Volume 2: Appendices

An in-flight research study was conducted utilizing the USAF/Total In-Flight Simulator (TIFS) to investigate longitudinal flying qualities for the flared landing approach phase of flight. A consistent set of data were generated for: determining what kind of command response the pilot prefers/requires in order to flare and land an aircraft with precision, and refining a time history criterion that took into account all the necessary variables and the characteristics that would accurately predict flying qualities. Seven evaluation pilots participated representing NASA Langley, NASA Dryden, Calspan, Boeing, Lockheed, and DFVLR (Braunschweig, Germany). The results of the first part of the study provide guidelines to the flight control system designer, using MIL-F-8785-(C) as a guide, that yield the dynamic behavior pilots prefer in flared landings. The results of the second part provide the flying qualities engineer with a derived flying qualities predictive tool which appears to be highly accurate. This time-domain predictive flying qualities criterion was applied to the flight data as well as six previous flying qualities studies, and the results indicate that the criterion predicted the flying qualities level 81% of the time and the Cooper-Harper pilot rating, within + or - 1%, 60% of the time.

Weingarten, Norman C.↗

Groundspeed filtering for CTAS

Ground speed is one of the radar observables which is obtained along with position and heading from NASA Ames Center radar. Within the Center TRACON Automation System (CTAS), groundspeed is converted into airspeed using the wind speeds which CTAS obtains from the NOAA weather grid. This airspeed is then used in the trajectory synthesis logic which computes the trajectory for each individual aircraft. The time history of the typical radar groundspeed data is generally quite noisy, with high frequency variations on the order of five knots, and occasional 'outliers' which can be significantly different from the probable true speed. To try to smooth out these speeds and make the ETA estimate less erratic, filtering of the ground speed is done within CTAS. In its base form, the CTAS filter is a 'moving average' filter which averages the last ten radar values. In addition, there is separate logic to detect and correct for 'outliers', and acceleration logic which limits the groundspeed change in adjacent time samples. As will be shown, these additional modifications do cause significant changes in the actual groundspeed filter output. The conclusion is that the current ground speed filter logic is unable to track accurately the speed variations observed on many aircraft. The Kalman filter logic however, appears to be an improvement to the current algorithm used to smooth ground speed variations, while being simpler and more efficient to implement. Additional logic which can test for true 'outliers' can easily be added by looking at the difference in the a priori and post priori Kalman estimates, and not updating if the difference in these quantities is too large.

Slater, Gary L.↗

Shock Tube Studies of Soot Formation

The objective of this research is to record the time histories of the major and minor species which appear in the pyrolysis of toluene, benzene, butadiene, allene, and acetylene; to develop a set of reactions that will model the observed profiles over a wide temperature and concentration range; to identify the critical reactions that influence the pre-particle soot formation process. Toluene and benzene were chosen as two key aromatic compounds which are representative of the pyrolytic process. Butadiene, allene, and acetylene were selected to investigate the formation of aromatic compounds from non-cyclic species. The experimental apparatus used for the study consists of a shock tube coupled to a time-of-flight mass spectrometer Spectra are recorded at 30 microsecond intervals for a total observation time of 0.50 - 1.20 milliseconds. Peak heights of the species of interest in the m/e range 12-300 are measured as a function of reaction time. Calibration curves are constructed which aid the conversion of peak heights to concentrations. The mixtures range from 1 percent-6 percent fuel; the balance is neon diluent.

R D Kern↗