Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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 703 records · Page 39

The analytical representation of viscoelastic material properties using optimization techniques

This report presents a technique to model viscoelastic material properties with a function of the form of the Prony series. Generally, the method employed to determine the function constants requires assuming values for the exponential constants of the function and then resolving the remaining constants through linear least-squares techniques. The technique presented here allows all the constants to be analytically determined through optimization techniques. This technique is employed in a computer program named PRONY and makes use of commercially available optimization tool developed by VMA Engineering, Inc. The PRONY program was utilized to compare the technique against previously determined models for solid rocket motor TP-H1148 propellant and V747-75 Viton fluoroelastomer. In both cases, the optimization technique generated functions that modeled the test data with at least an order of magnitude better correlation. This technique has demonstrated the capability to use small or large data sets and to use data sets that have uniformly or nonuniformly spaced data pairs. The reduction of experimental data to accurate mathematical models is a vital part of most scientific and engineering research. This technique of regression through optimization can be applied to other mathematical models that are difficult to fit to experimental data through traditional regression techniques.

Hill, S. A.↗

Constraints on the nucleus and dust properties from mid-infrared imaging of comet Hyakutake

We use Mie scattering theory to determine the expected thermal emission from dust grains in cometary comae and apply these results to mid-infrared images of comet Hyakutake (C/1996 B2) obtained preperihelion in 1996 March. Calculations were performed for dust grains in the size range from 0.1 to 10 micrometers for two different compositions: amorphous olivine (a silicate glass) and an organic residue mixture. The resulting emission efficiencies are complicated functions of wavelength and particle size and are significantly different for the two materials in question. The Hyakutake data set consists of three nights of high-resolution imaging (100-150 km pixel-1 at the comet) of the inner coma at 8.7, 11.7, 12.5, and 19.7 micrometers. Attempts to fit the observed colors (ratios of fluxes at different wavelengths) using a single grain composition failed. However, fits to the data were achieved for all three nights using a mixture of approximately 1 micrometer olivine grains and approximately 7 micrometers organic grains. The resulting olivine mass fraction was between 8% and 16% of the total dust mass-loss rate. We also estimate the radius of the nucleus to be r = 2.1 +/- 0.4 km.

Non-NASA Center↗

Temperature and EUV Intensity in a Coronal Prominence Cavity and Streamer

We analyze the temperature and EUV line emission of a coronal cavity and surrounding streamer in terms of a morphological forward model. We use a series of iron line ratios observed with the Hinode Extreme-ultraviolet Imaging Spectrograph (EIS) on 2007 Aug. 9 to constrain temperature as a function of altitude in a morphological forward model of the streamer and cavity. We also compare model prediction of the EIS EUV line intensities and polarized brightness (pB) data from the Mauna Loa Solar Observatory (MLSO) MK4. This work builds on earlier analysis using the same model to determine geometry of and density in the same cavity and streamer. The fit to the data with altitude dependent temperature profiles indicates that both the streamer and cavity have temperatures in the range 1.4-1.7 MK. However, the cavity exhibits substantial substructure such that the altitude dependent temperature profile is not sufficient to completely model conditions in the cavity. Coronal prominence cavities are structured by magnetism so clues to this structure are to be found in their plasma properties. These temperature substructures are likely related to structures in the cavity magnetic field. Furthermore, we find that the model overestimates the line intensities by a factor of 4-10, while overestimating pB data by no more than a factor of 1.4. One possible explanation for this is that there may be a significant amount of material at temperatures outside of the range log T(K) approximately equals 5.8 - 6.7 in both the cavity and the streamer.

Kucera, T. A.↗

Analysis of Dual Rotating Rake Data from the NASA Glenn Advanced Noise Control Fan Duct with Artificial Sources

The Rotating Rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. Initially, the mode amplitudes and phases were quantified from a single rake measurement at one axial location. To directly measure the modes propagating in both directions within a duct, a second rake was mounted to the rotating system with an offset in both the axial and the azimuthal directions. The rotating rake data analysis technique was then extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode levels at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode amplitudes for the modes propagating in both directions within the duct. Validation experiments have been conducted using artificial acoustic sources. Results are shown for the measurement of the standing waves in the duct from sound generated by one and two acoustic sources that are separated into the component modes propagating in both directions within the duct. Measured reflection coefficients from the open end of the duct are compared to analytical predictions.

Noise↗

Dual Rotating Rake Measurements of Higher-Order Duct Modes: Validation Using Experimental and Numerical Data

A rotating rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. After analysis of the measured data, the mode coefficient amplitudes and phases were quantified. Early studies using this system found that mode power levels computed from rotating rake measured data would agree with the far-field power levels. However, this agreement required that the sound from the noise sources within the duct propagated outward from the duct exit without reflection and previous studies suggested conditions could exist where significant reflections could occur. This paper shows that mounting a second rake to the rotating system, with an offset in both the axial and the azimuthal directions, measures the data necessary to determine the modes propagating in both directions within a duct. The rotating rake data analysis technique was extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode coefficients at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode coefficients for the modes propagating in both directions within the duct while accounting for the presence of evanescent modes. The validation of the dual-rotating-rake measurements was conducted using data from a combination of experiments and numerical calculations to compute reflection coefficients and other mode coefficient ratios. Compared to results from analytical and numerical computations, the results from dual-rotating-rake measured data followed the expected trends when frequency, mode number, and duct termination geometry were changed.

Validation↗

Analysis of Physical Properties of Dust Suspended in the Mars Atmosphere

Methods for iteratively determining the infrared optical constants for dust suspended in the Mars atmosphere are described. High quality spectra for wavenumbers from 200 to 2000 1/cm were obtained over a wide range of view angles by the Mariner 9 spacecraft, when it observed a global Martian dust storm in 1971-2. In this research, theoretical spectra of the emergent intensity from Martian dust clouds are generated using a 2-stream source-function radiative transfer code. The code computes the radiation field in a plane-parallel, vertically homogeneous, multiply scattering atmosphere. Calculated intensity spectra are compared with the actual spacecraft data to iteratively retrieve the optical properties and opacity of the dust, as well as the surface temperature of Mars at the time and location of each measurement. Many different particle size distributions a-re investigated to determine the best fit to the data. The particles are assumed spherical and the temperature profile was obtained from the CO2 band shape. Given a reasonable initial guess for the indices of refraction, the searches converge in a well-behaved fashion, producing a fit with error of less than 1.2 K (rms) to the observed brightness spectra. The particle size distribution corresponding to the best fit was a lognormal distribution with a mean particle radius, r(sub m) 0.66 pm, and variance, omega(sup 2) = 0.412 (r(sub eff) = 1.85 microns, v(sub eff) =.51), in close agreement with the size distribution found to be the best fit in the visible wavelengths in recent studies. The optical properties and the associated single scattering properties are shown to be a significant improvement over those used in existing models by demonstrating the effects of the new properties both on heating rates of the Mars atmosphere and in example spectral retrieval of surface characteristics from emission spectra.

Snook, Kelly↗

Prediction Model for the Life of Nickel-cadmium Batteries in Geosynchronous Orbit Satellites

A mathematical model is described which predicts the service life of nickel-cadmium batteries designed for geosynchronous orbit satellites. Regression analysis technique is used to analyze orbital data on second generation trickle charged batteries. The model gives average cell voltage as a function of design parameters, operating parameters and time. The voltage model has the properties of providing a good fit to the data, good predictive capability, and agreement with known battery performance characteristics. Average cell voltage can be predicted to within 0.02 volts for up to 8 years. This modeling shows that these batteries will operate reliably for 10 years. Third-generation batteries are expected to operate even longer.

Engleman, J. H.↗

Analysis of Wind Tunnel Oscillatory Data of the X-31A Aircraft

Wind tunnel oscillatory tests in pitch, roll, and yaw were performed on a 19%-scale model of the X-31A aircraft. These tests were used to study the aerodynamic characteristics of the X-31A in response to harmonic oscillations at six frequencies. In-phase and out-of-phase components of the aerodynamic coefficients were obtained over a range of angles of attack from 0 to 90 deg. To account for the effect of frequency on the data, mathematical models with unsteady terms were formulated by use of two different indicial functions. Data from a reduced set of frequencies were used to estimate model parameters, including steady-state static and dynamic stability derivatives. Both models showed good prediction capability and the ability to accurately fit the measured data. Estimated static stability derivatives compared well with those obtained from static wind tunnel tests. The roll and yaw rate derivative estimates were compared with rotary-balanced wind tunnel data and theoretical predictions. The estimates and theoretical predictions were in agreement at small angles of attack. The rotary-balance data showed, in general, acceptable agreement with the steady-state derivative estimates.

Smith, Mark S.↗

Graphical User Interface Development for Representing Air Flow Patterns

In the Turbine Branch, scientists carry out experimental and computational work to advance the efficiency and diminish the noise production of jet engine turbines. One way to do this is by decreasing the heat that the turbine blades receive. Most of the experimental work is carried out by taking a single turbine blade and analyzing the air flow patterns around it, because this data indicates the sections of the turbine blade that are getting too hot. Since the cost of doing turbine blade air flow experiments is very high, researchers try to do computational work that fits the experimental data. The goal of computational fluid dynamics is for scientists to find a numerical way to predict the complex flow patterns around different turbine blades without physically having to perform tests or costly experiments. When visualizing flow patterns, scientists need a way to represent the flow conditions around a turbine blade. A researcher will assign specific zones that surround the turbine blade. In a two-dimensional view, the zones are usually quadrilaterals. The next step is to assign boundary conditions which define how the flow enters or exits one side of a zone. way of setting up computational zones and grids, visualizing flow patterns, and storing all the flow conditions in a file on the computer for future computation. Such a program is necessary because the only method for creating flow pattern graphs is by hand, which is tedious and time-consuming. By using a computer program to create the zones and grids, the graph would be faster to make and easier to edit. Basically, the user would run a program that is an editable graph. The user could click and drag with the mouse to form various zones and grids, then edit the locations of these grids, add flow and boundary conditions, and finally save the graph for future use and analysis. My goal this summer is to create a graphical user interface (GUI) that incorporates all of these elements. I am writing the program in Java, a language that is portable among platforms, because it can run on different operating systems such as Windows and Unix without having to be rewritten. I had no prior experience of programming in Java at the start of my internship; I am continuously learning as I create the program. I have written the part of the program that enables a user to draw several zones, edit them, and store their locations. The next phase of my project is to allow the user to click on the side of a zone and create a boundary condition for it. A previous intern wrote a program that allows the user to input boundary conditions. I can integrate the two programs to create a larger, more usable program. After that, I will develop a way for the user to save the graph for future reference. Another eventual goal is to make the GUI capable of creating three-dimensional zones as well. Researchers such as my mentor, Dr. David Ashpis, need a quick, user-friendly

Chaudhary, Nilika↗

SEE Rate Estimation: Model Complexity and Data Requirements

Statistical Methods outlined in [Ladbury, TNS20071 can be generalized for Monte Carlo Rate Calculation Methods Two Monte Carlo Approaches: a) Rate based on vendor-supplied (or reverse-engineered) model SEE testing and statistical analysis performed to validate model; b) Rate calculated based on model fit to SEE data Statistical analysis very similar to case for CREME96. Information Theory allows simultaneous consideration of multiple models with different complexities: a) Model with lowest AIC usually has greatest predictive power; b) Model averaging using AIC weights may give better performance if several models have similar good performance; and c) Rates can be bounded for a given confidence level over multiple models, as well as over the parameter space of a model.

Ladbury, Ray↗

OSSE observations of 3C 273

We report results of multiple observations of the quasar 3C 273 with the Oriented Scintillation Spectrometer Experiment (OSSE) instrument on the Compton Gamma Ray Observatory. These observations span the period from 1991 June through 1993 January and represent the most sensitive observations to date in low-energy gamma rays. The source was detected at historically weak 100 keV fluxes compared with previous measurements. Variability by factors of approximately 3 on timescales of approximately equal 2 months was observed in the energy band 50-150 keV. The data are well described by a single power law with a proton number index Gamma = 1.7 +/- 0.1. No significant change of Gamma was observed during changes in intensity. Thermal models do not provide acceptable fits to the data. When the OSSE data are combined with contemporaneous measurements by COMPTEL and EGRET, the spectrum is seen to break at an energy of 1.0(+0.9, -0.4) MeV to a softer power law with Delta Gamma = 0.7(+0.06, -0.11), forming a power law with Gamma = 2.4 between approximately 1 MeV and several GeV.

Johnson, W. N.↗

Preprocessing of multispectral data and simulation of ERTS data channels to make computer terrain maps of a Yellowstone National Park test site

The possibility of improving the accuracy of terrain classification by preprocessing spectral data was investigated. Terrain maps were made using the following techniques: 1) preprocessing by scan angle function transformation, using the computer-selected best set of three channels; and 2) preprocessing by ratio transformation, using the specified ERTS data channels, simulated by fitting the spectral response of each of the 12 data channels to the ERTS channels by a set of weighting coefficients. By using a simple technique during printout, the maps were produced in color. The normalized scan angle function transformation resulted in the most accurate classification. The best ratio transformation for the Yellowstone Park data was the ratio of each channel to the sum of all channels. A supervised training program involving maximum likelihood decision for selecting the best spectrometer channels and similar techniques for digitizing the data of the analog magnetic tapes were used. Cloud shadows were recognized in addition to eight classes of terrain. Preprocessing of data resulted in more accurate maps, required fewer training areas (hence less preparation and computer time), and enabled much of the area formerly classified as shadow to be reclassified according to actual terrain type.

Smedes, H. W.↗

LORAN-C data reduction at the US Naval Observatory

As part of its mission and in cooperation with the U.S. Coast Guard, the U.S. Naval Observatory (USNO) monitors and reports the timing of the LORAN-C chains. The procedures for monitoring and processing the reported values have evolved with advances in monitoring equipment, computer interfaces and PCs. This paper discusses the current standardized procedures used by USNO to sort the raw data according to Group Repetition Interval (GRI) rate, to fit and smooth the data points, and, for chains remotely monitored, to tie the values to the USNO Master Clock. The results of these procedures are the LORAN time of transmission values, as references to UTC(USNO) (Universal Coordinated Time) for all LORAN chains. This information is available to users via USNO publications and the USNO Automated Data Service (ADS).

Harold Chadsey↗

The Parameterization of Top-Hat Particle Sensors with Microchannel-Plate-Based Detection Systems and its Application to the Fast Plasma Investigation on NASA's Magnetospheric MultiScale Mission

The most common instrument for low energy plasmas consists of a top-hat electrostatic analyzer geometry coupled with a microchannel-plate (MCP)-based detection system. While the electrostatic optics for such sensors are readily simulated and parameterized during the laboratory calibration process, the detection system is often less well characterized. Furthermore, due to finite resources, for large sensor suites such as the Fast Plasma Investigation (FPI) on NASA's Magnetospheric Multiscale (MMS) mission, calibration data are increasingly sparse. Measurements must be interpolated and extrapolated to understand instrument behavior for untestable operating modes and yet sensor inter-calibration is critical to mission success. To characterize instruments from a minimal set of parameters we have developed the first comprehensive mathematical description of both sensor electrostatic optics and particle detection systems. We include effects of MCP efficiency, gain, scattering, capacitive crosstalk, and charge cloud spreading at the detector output. Our parameterization enables the interpolation and extrapolation of instrument response to all relevant particle energies, detector high voltage settings, and polar angles from a small set of calibration data. We apply this model to the 32 sensor heads in the Dual Electron Sensor (DES) and 32 sensor heads in the Dual Ion Sensor (DIS) instruments on the 4 MMS observatories and use least squares fitting of calibration data to extract all key instrument parameters. Parameters that will evolve in flight, namely MCP gain, will be determined daily through application of this model to specifically tailored in-flight calibration activities, providing a robust characterization of sensor suite performance throughout mission lifetime. Beyond FPI, our model provides a valuable framework for the simulation and evaluation of future detection system designs and can be used to maximize instrument understanding with minimal calibration resources.

Instrument-model↗

The Dielectric Constant at P-Band for Salinity from 0 to 150 pss

Measurements have been made at P-band (0.707 GHz) to construct a model for the dielectric constant of sea water and extend the model for the dielectric constant to high salinity (50-150 pss). The measurements are part of research to develop a model for the dielectric constant suitable for future wide-bandwidth remote sensing of salinity and for application to water bodies such as the Great Salt Lake with salinity significantly above that found in the open ocean. Measurements have been made at temperatures from 2 to 30 °C and salinity from 0 to 138 pss. The data have been fit to a Debye model for the dielectric constant with a single resonance as has been employed at L-band (1.413 GHz) where remote sensing of salinity is currently done. Comparison with contemporary models developed from data at L-band indicates that the L-band model and new P-band model do well at both frequencies for salinity less than 50 pss but at higher values of salinity the L-band models diverge from the data. The data has also been used to test at high salinity the mathematical relationship between salinity and conductivity which is the basis for the practical salinity scale (pss).

Dielectric Constant↗

The Atmospheric Structure of Triton and Pluto

The goal of this research was to better determine the atmospheric structures of Triton and Pluto through further analysis of three occultation data sets obtained with the Kuiper Airborne Observatory (KAO.) As the research progressed, we concentrated our efforts on the Triton data, as this appeared to be the most fruitful. Three papers have been prepared as a result of this research. The first paper presents new results about Triton's atmospheric structure from the analysis of all ground-based stellar occultation data recorded to date, including one single-chord occultation recorded on 1993 July 10 and nine occultation lightcurves from the double-star event on 1995 August 14. These stellar occultation observations made both in the visible and in the infrared have good spatial coverage of Triton, including the first Triton central-flash observations, and are the first data to probe the altitude level 20-100 km on Triton. The small-planet lightcurve model of J. L. Elliot and L. A. Young was generalized to include stellar flux refracted by the far limb, and then fitted to the data. Values of the pressure, derived from separate immersion and emersion chords, show no significant trends with latitude, indicating that Triton's atmosphere is spherically symmetric at approximately 50 km altitude to within the error of the measurements; however, asymmetry observed in the central flash indicates the atmosphere is not homogenous at the lowest levels probed (approximately 20 km altitude). From the average of the 1995 occultation data, the equivalent isothermal temperature of the atmosphere is 47 plus or minus 1 K and the atmospheric pressure at 1400 km radius (approximately 50 km altitude) is 1.4 plus or minus 0.1 microbar. Both of these are not consistent with a model based on Voyager UVS and RSS observations in 1989. The atmospheric temperature from the occultation is 5 K colder than that predicted by the model and the observed pressure is a factor of 1.8 greater than the model. In our opinion, the disagreement in temperature and pressure is probably due to modeling problems at the microbar level, since measurements at this level have not previously been made. Alternatively, the difference could be due to seasonal change in Triton's atmospheric structure. The second paper reports observations of a recent stellar occultation by Triton which, when combined with earlier results, show that Triton has undergone a period of global warming since 1989. The most conservative estimates of the rate of temperature and surface-pressure increase during this period imply that the atmosphere is doubling in bulk every 10 years -- significantly faster than predicted by published frost model for Triton. Our results suggests that permanent polar caps on Triton play a dominant role in regulating seasonal atmospheric changes. Similar processes should also be active on Pluto. A third paper 'Global Warming on Triton' will appear in a the January 1999 issue of Sky and Telescope.

Elliot, James L.↗

A Probabilistic Mass Estimation Algorithm for a Novel 7- Channel Capacitive Sample Verification Sensor

A document describes an algorithm created to estimate the mass placed on a sample verification sensor (SVS) designed for lunar or planetary robotic sample return missions. A novel SVS measures the capacitance between a rigid bottom plate and an elastic top membrane in seven locations. As additional sample material (soil and/or small rocks) is placed on the top membrane, the deformation of the membrane increases the capacitance. The mass estimation algorithm addresses both the calibration of each SVS channel, and also addresses how to combine the capacitances read from each of the seven channels into a single mass estimate. The probabilistic approach combines the channels according to the variance observed during the training phase, and provides not only the mass estimate, but also a value for the certainty of the estimate. SVS capacitance data is collected for known masses under a wide variety of possible loading scenarios, though in all cases, the distribution of sample within the canister is expected to be approximately uniform. A capacitance-vs-mass curve is fitted to this data, and is subsequently used to determine the mass estimate for the single channel s capacitance reading during the measurement phase. This results in seven different mass estimates, one for each SVS channel. Moreover, the variance of the calibration data is used to place a Gaussian probability distribution function (pdf) around this mass estimate. To blend these seven estimates, the seven pdfs are combined into a single Gaussian distribution function, providing the final mean and variance of the estimate. This blending technique essentially takes the final estimate as an average of the estimates of the seven channels, weighted by the inverse of the channel s variance.

Wolf, Michael↗

Interaction of the equatorial midnight pressure bulge and thermospheric zonal winds

The zonal component of the thermospheric neutral wind in the equatorial region is found to flow eastward all night, having a minimum near midnight. The data, obtained by Spencer et al. on the Dynamics Explorer-2 (DE-2) satellite, and reported more recently by Wharton et al., indicate that the minimum feature is significant. The minimum is illustrated by a 4th order Fourier series fit to the data points reported by Wharton et al. A table gives the amplitudes and phases of the four harmonic components used with the estimated error for each amplitude. It is asked whether this behavior is due to the local passage of the equatorial midnight pressure bulge associated with the midnight temperature anomaly. If so, it is to be expected that other data characterizing the midnight temperature anomaly should show consistency with this observation, especially in view of the large changes observed in the zonal velocity in the midnight sector. Consistency in momentum conservation may be checked using available data from previous independent experiments. Such data, taken under similar conditions of solar activity, was substituted into the momentum equation and found to be consistent using a priori estimates of the effect of viscosity.

Herrero, F. A.↗