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 235 records · Page 13

Determination of sounding rocket position and attitude from radar and magnetometer data

The techniques are described that were used to determine the trajectories and orientations of three sounding rockets instrumented to study the aurora. The radar plot board data were fitted to a near earth expansion of the force of gravity to determine the trajectory. Only onboard magnetometer data were used to determine the attitude of the payload with respect to the earth's magnetic field. Computer programs in the FORTRAN language are available which generate the trajectory and attitude parameters.

Pongratz, M.↗

Local area analysis of high-degree solar oscillations: New ring fitting procedures

The local area analysis of five-minute solar oscillations using ring diagrams to determine subphotospheric velocity flows is a tool for convection zone dynamics. In relation to the problem of the large computational task of fitting the rings, a faster method is presented that carries out the ring fitting using data obtained with a high l helioseismometer. Noise sources are eliminated, and a perturbation approach is used to fit the azimuthally averaged spectrum. The parameters determined in this way are held constant while the ring diagram is fitted. The results obtained are presented and discussed.

Haber, Deborah A.↗

Application of a Bayesian Framework for Plasticity Model Selection

Interpretable Machine Learning (IML) has performed well when tasked with deriving constitutive material models. However, IML has been shown to prefer models that overfit noise in data, which tends to lead to bloat and a decrease in interpretability. Due to these issues, the ability of IML to reliably derive models that fit the data and are both interpretable and generalizable is limited. A method developed recently has shown promise to improve upon traditional IML by using a Bayesian fitness definition for the evolution of free-form models with non-deterministic parameters. This framework was developed for genetic-programming-based symbolic regression(GPSR) and involves model parameter estimation using Sequential Monte Carlo sampling (SMC).The method has demonstrated a reduction in bloat when dealing with noisy data in comparison to conventional GPSR. The results of this framework applied to stress-strain data for copper show models that more effectively predict the experimental data better than was previously shown with GPSR.

plasticity↗

A model of Saturn's magnetic field based on all available data

A model of Saturn's magnetic field, called Saturn Pioneer Voyager (SPV), is developed on the basis of an analysis of three sets of data: those from the Pioneer 11, the Voyager 1, and the Voyager 2 magnetometers. It is shown that the SPV model fits the data observed between 1.3 and 8.0 Saturn radii from the planet's center with a 1.13 percent weighted rms average of the percent differences between the observed and modeled fields, which is substantially better than the fits yielded by any of the previous models.

Davis, Leverett, Jr.↗

Properties of Martian Hematite at Meridiani Planum by Simultaneous Fitting of Mars Mossbauer Spectra

Mossbauer spectrometers [1] on the two Mars Exploration Rovers (MERs) have been making measurements of surface rocks and soils since January 2004, recording spectra in 10-K-wide temperature bins ranging from 180 K to 290 K. Initial analyses focused on modeling individual spectra directly as acquired or, to increase statistical quality, as sums of single-rock or soil spectra over temperature or as sums over similar rock or soil type [2, 3]. Recently, we have begun to apply simultaneous fitting procedures [4] to Mars Mossbauer data [5-7]. During simultaneous fitting (simfitting), many spectra are modeled similarly and fit together to a single convergence criterion. A satisfactory simfit with parameter values consistent among all spectra is more likely than many single-spectrum fits of the same data because fitting parameters are shared among multiple spectra in the simfit. Consequently, the number of variable parameters, as well as the correlations among them, is greatly reduced. Here we focus on applications of simfitting to interpret the hematite signature in Moessbauer spectra acquired at Meridiani Planum, results of which were reported in [7]. The Spectra. We simfit two sets of spectra with large hematite content [7]: 1) 60 rock outcrop spectra from Eagle Crater; and 2) 46 spectra of spherule-rich lag deposits (Table 1). Spectra of 10 different targets acquired at several distinct temperatures are included in each simfit set. In the table, each Sol (martian day) represents a different target, NS is the number of spectra for a given sol, and NT is the number of spectra for a given temperature. The spectra are indexed to facilitate definition of parameter relations and constraints. An example spectrum is shown in Figure 1, together with a typical fitting model. Results. We have shown that simultaneous fitting is effective in analyzing a large set of related MER Mossbauer spectra. By using appropriate constraints, we derive target-specific quantities and the temperature dependence of certain parameters. By examining different fitting models, we demonstrate an improved fit for martian hematite modeled with two sextets rather than as a single sextet, and show that outcrop and spherule hematite are distinct. For outcrop, the weaker sextet indicates a Morin transition typical of well-crystallized and chemically pure hematite, while most of the outcrop hematite remains in a weakly ferromagnetic state at all temperatures. For spherule spectra, both sextets are consistent with weakly ferromagnetic hematite with no Morin transition. For both hematites, there is evidence for a range of particle sizes.

Agresti, D. G.↗

Effective Tree Scattering at L-Band

For routine microwave Soil Moisture (SM) retrieval through vegetation, the tau-omega [1] model [zero-order Radiative Transfer (RT) solution] is attractive due to its simplicity and eases of inversion and implementation. It is the model used in baseline retrieval algorithms for several planned microwave space missions, such as ESA's Soil Moisture Ocean Salinity (SMOS) mission (launched November 2009) and NASA's Soil Moisture Active Passive (SMAP) mission (to be launched 2014/2015) [2 and 3]. These approaches are adapted for vegetated landscapes with effective vegetation parameters tau and omega by fitting experimental data or simulation outputs of a multiple scattering model [4-7]. The model has been validated over grasslands, agricultural crops, and generally light to moderate vegetation. As the density of vegetation increases, sensitivity to the underlying SM begins to degrade significantly and errors in the retrieved SM increase accordingly. The zero-order model also loses its validity when dense vegetation (i.e. forest, mature corn, etc.) includes scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. The tau-omega model (when applied over moderately to densely vegetated landscapes) will need modification (in terms of form or effective parameterization) to enable accurate characterization of vegetation parameters with respect to specific tree types, anisotropic canopy structure, presence of leaves and/or understory. More scattering terms (at least up to first-order at L-band) should be included in the RT solutions for forest canopies [8]. Although not really suitable to forests, a zero-order tau-omega model might be applied to such vegetation canopies with large scatterers, but that equivalent or effective parameters would have to be used [4]. This requires that the effective values (vegetation opacity and single scattering albedo) need to be evaluated (compared) with theoretical definitions of these parameters. In a recent study [9], effective vegetation opacity of coniferous trees was compared with two independent estimates of the same parameter. First, a zero-order RT model was fitted to multiangular microwave emissivity data in a least-square sense to provide effective vegetation optical depth as done in spaceborne retrieval algorithms. Second, a ratio between radar backscatter measurements with a corner reflector under trees and in an open area was calculated to obtain measured tree propagation characteristics. Finally, the theoretical propagation constant was determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). Results indicated that the effective attenuation values are smaller than but of similar magnitude to both the theoretical and measured values. This study will complement the previous work [9] and will focus on characterization of effective scattering albedo by assuming that effective vegetation opacity is same as theoretical opacity. The resultant effective albedo will not be the albedo of single forest canopy element anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy including multiple scattering as described.

Kurum, Mehmet↗

Kinetics of the reaction of nitric oxide with hydrogen

The reaction of nitric oxide with hydrogen has been studied in the temperature range 2400-4500 K using a shock-tube technique. Mixtures of NO and H2 diluted in argon or krypton were heated by incident shock waves, and the infrared emission from the fundamental vibration-rotation band of NO at 5.3 microns was used to monitor the time-varying NO concentration. The decomposition of nitric oxide behind the shock was found to be modeled well by a fifteen-reaction system. A principal result of the study was the determination of the rate constant k1 for the reaction H + NO yields N + OH, which may be the rate-limiting step for NO removal in some combustion systems. Experimental values of k1 were obtained for each test through comparisons of measured and numerically predicted NO profiles. The data are fit closely by the expression k1 = 1.34 times 10 to the fourteenth power exp(-49 200/RT) cu cm/mole-sec. These data appear to be the first available for this rate constant.

Flower, W. L.↗

A simulation of air pollution model parameter estimation using data from a ground-based LIDAR remote sensor

One way to obtain estimates of the unknown parameters in a pollution dispersion model is to compare the model predictions with remotely sensed air quality data. A ground-based LIDAR sensor provides relative pollution concentration measurements as a function of space and time. The measured sensor data are compared with the dispersion model output through a numerical estimation procedure to yield parameter estimates which best fit the data. This overall process is tested in a computer simulation to study the effects of various measurement strategies. Such a simulation is useful prior to a field measurement exercise to maximize the information content in the collected data. Parametric studies of simulated data matched to a Gaussian plume dispersion model indicate the trade offs available between estimation accuracy and data acquisition strategy.

Kibler, J. F.↗

Present-day plate motions

A data set comprising 110 spreading rates, 78 transform fault azimuths and 142 earthquake slip vectors was inverted to yield a new instantaneous plate motion model, designated RM2. The mean averaging interval for the relative motion data was reduced to less than 3 My. A detailed comparison of RM2 with angular velocity vectors which best fit the data along individual plate boundaries indicates that RM2 performs close to optimally in most regions, with several notable exceptions. On the other hand, a previous estimate (RM1) failed to satisfy an extensive set of new data collected in the South Atlantic Ocean. It is shown that RM1 incorrectly predicts the plate kinematics in the South Atlantic because the presently available data are inconsistent with the plate geometry assumed in deriving RM1. It is demonstrated that this inconsistency can be remedied by postulating the existence of internal deformation with the Indian plate, although alternate explanations are possible.

Minster, J. B.↗

Extension of the MSIS thermosphere model into the middle and lower atmosphere

The MSIS-86 empirical model has been revised in the lower thermosphere and extended into the mesosphere and lower atmosphere to provide a single analytic model for calculating temperature and density profiles representative of the climatological average for various geophysical conditions. Tabulations from the Handbook for MAP 16 (Barnett and Corney, 1985) are the primary guide for the lower atmosphere and are supplemented by historical rocket and incoherent scatter data in the upper mesosphere and lower thermosphere. Low-order spherical harmonics and Fourier series are used to describe the major variations throughout the atmosphere including latitude, annual, semiannual, and simplified local time and longitude variations. While month to month details cannot be completely represented, lower atmosphere temperature data are fit to an overall standard deviation of 3 K and pressure to 2 percent. Comparison with rocket and other data indicates that the model represents current knowledge of the climatological average reasonably well, although there is some conflict as to details near the mesopause.

Hedin, A. E.↗

Safety Aboard the International Space Station

As with any task that NASA takes on, safety is of utmost importaqce. There are pages of safety codes and procedures that must be followed before any idea can be brought to life. Unfortunately, the International Space Station s (ISS) safety regulations and procedures are based on lg standards rather than on Og. To aide in making this space age home away from home a less hazardous environment, I worked on several projects revolving around the dangers of flammable items in microgravity. The first task I was assigned was to track flames. This involves turning eight millimeter video recordings, of tests run in the five second drop tower, into avi format on the computer. The footage is then compressed and altered so that the flame can be seen more clearly. Using another program called Spotlight, line profiles were used to collect data describing the luminescence of the flame at different points. These raw data are saved as text files and run trough a macro so that a Matlab program can analyze it. By fitting the data to a curve and determining the areas of brightest luminescence, the behavior of the flame can be recorded numerically. After entering the data into a database, researchers can come back later and easily get information on flames resulting from different gas and liquid mixtures in microgravity. I also worked on phase two of the FATE project, which deals with safety aboard the ISS. This phase involves igniting projected droplets and determining how they react with secondary materials. Such simulations represent, on a small scale, the spread of onboard fires due to the effervescence of burning primary materials. I set up existing hardware to operate these experiments and ran tests with it, photographing the results. I also made CAD drawings of the apparatus and the area available on the (SF)2 rig for it to fit into. The experiment will later be performed on the KC-135, and the results gathered will be used to reanalyze current safety standards for the ISS, including the distance of required separation for flammable materials. Additional information is included in the original extended abstract.

Mintz, Shauna M.↗

Reanalysis of X-ray emission from M87. 1: The single-phase medium

We reanalyze the density and temperature profiles of the X-ray-emitting gas around M87, assuming a spherically symmetric, single-phase model. For given assumed density and temperature profiles, we predict Einstein High Resolution Imager (HRI) and Imaging Proportional Counter (IPC) surface brightness distribution as well as Focal Point Crystal Spectrometer (FPCS) line fluxes, which we compare with the data. We find that a good fit to these data can be obtained and that, with a suitable adjustment of the abundances, the equivalent width of the iron complex at 7 keV as observed by wide-beam instruments can also be explained. In contrast to this, the Einstein Solid State Spectrometer (SSS) spectrum and optically determined mass profiles cannot be accounted for simultaneously with the previous X-ray data. We show that the disagreement of the Solid State Spectrometer with our models is due to an inconsistency of the Solid State Spectrometer data with that of the other X-ray instruments. Because of this, we find no evidence for the existence of absorption above the Galactic value, contrary to an earlier claim based on Solid State Spectrometer data. The disagreement of our models with optical mass data is due either to the assumption of a single-phase model, thus rendering the model unacceptable, or possibly to maximal systematic errors in the FPCS data. We also show that thermal conduction at the Spitzer rate is important in our models, although this conclusion is of limited importance in view of the poor fit of our models to the data.

Tsai, John C.↗

Culling of redundant data

The inversion of a large data set, with it errors, is demonstrated and the tradeoff curve, or resolving power calculations, are discussed. Consideration is given to inverse problems of the earth, represented in the idealized form of a single scaler function of a single coordinate. The difference between the functional for the real earth and the functional for some model is formulated. A procedure is given for using the linear relationships between the differences in the data and the differences in the model as side conditions. Through a process of diagonalization a model is created that fits the data. The data are culled by accepting only those relative standard errors that are less than 100 percent.

Gilbert, F.↗

Gravity - Mare Humorum.

Global tracking coverage of the Apollo 15 subsatellite has provided gravity measurements from 50 km altitudes over the entire Humorum basin. An estimate of surface mass points at 2 degree intervals, which best fit the data reveals a mass distribution having a lesser mass excess for the very central area. When two different profiles were fit using a two disk model it was again found for each profile that the smaller central disk decreased the central mass by approximately 30%. The mass distribution per unit area however for the major portion of the mascon is still consistent with the other mascons (Crisium, Nectaris and Serenitatis) of approximately 800-900 kg/sq cm. The surface mass point solution seems to correlate somewhat with the dark areas on Whitaker's IR-UV map.

Sjogren, W. L.↗

Constraints on the NH3 and PH3 distributions in the Great Red Spot

Medium resolution (10 A) UV spectra were obtained for the Great Red Spot (GRS) and South Tropical Zone (STZ) of Jupiter using the low dispersion mode of the IUE spectrometers at wavelengths from 1900-2200 A. The scans were carried out to determine the coloring agent for the GRS to improve the database for developing photochemical models of the feature. The wavelengths were selected to cover the absorption features of NH3 and PH3. The resulting data were interpreted using a vertically inhomogeneous Rayleigh scattering radiative transfer model. Various NH3 concentrations were explored in an effort to fit the data, taking into account changes which would occur at different atmospheric pressure levels and due to the projected temperature fields. A forbidden NH3/forbidden H2 mixing ratio that was calculated at the 80-125 mbar pressure level in the GRS was enhanced by 3-10 percent relative to the STZ. An upper limit was obtained for the mixing ratio of PH3 in the GRS that is significantly lower than previously predicted concentrations, implying that vertical transport in the GRS is not much greater than in adjacent regions.

Wagener, R.↗

Design and Development of a Model to Simulate 0-G Treadmill Running Using the European Space Agency's Subject Loading System

Develop a model that simulates a human running in 0 G using the European Space Agency s (ESA) Subject Loading System (SLS). The model provides ground reaction forces (GRF) based on speed and pull-down forces (PDF). DESIGN The theoretical basis for the Running Model was based on a simple spring-mass model. The dynamic properties of the spring-mass model express theoretical vertical GRF (GRFv) and shear GRF in the posterior-anterior direction (GRFsh) during running gait. ADAMs VIEW software was used to build the model, which has a pelvis, thigh segment, shank segment, and a spring foot (see Figure 1).the model s movement simulates the joint kinematics of a human running at Earth gravity with the aim of generating GRF data. DEVELOPMENT & VERIFICATION ESA provided parabolic flight data of subjects running while using the SLS, for further characterization of the model s GRF. Peak GRF data were fit to a linear regression line dependent on PDF and speed. Interpolation and extrapolation of the regression equation provided a theoretical data matrix, which is used to drive the model s motion equations. Verification of the model was conducted by running the model at 4 different speeds, with each speed accounting for 3 different PDF. The model s GRF data fell within a 1-standard-deviation boundary derived from the empirical ESA data. CONCLUSION The Running Model aids in conducting various simulations (potential scenarios include a fatigued runner or a powerful runner generating high loads at a fast cadence) to determine limitations for the T2 vibration isolation system (VIS) aboard the International Space Station. This model can predict how running with the ESA SLS affects the T2 VIS and may be used for other exercise analyses in the future.

Caldwell, E. C.↗

Ultimate limits to error probabilities for ionospheric models based on solar geophysical indices and how these compare with the state of the art

An ideal model based on a given set of geophysical indices is defined as a model that provides a least squares fit to the data set as a function of the indices considered. Satellite measurements of electron content for three stations at different magnetic latitudes were used to provide such data sets which were each fitted to the geophysical indices. The magnitude of the difference between the measured value and the derived equation for the data set was used to estimate the probability of making an error greater than a given magnitude for such an ideal model. Atmospheric Explorer C data is used to examine the causes of the fluctuations and suggestions are made about how real improvements can be made in ionospheric forecasting ability. Joule heating inputs in the auroral electrojets are related to the AL and AU magnetic indices. Magnetic indices based on the time integral of the energy deposited in the electrojets are proposed for modeling processes affected by auroral zone heating.

Nisbet, J. S.↗

Computer Modeling of Non-Isothermal Crystallization

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homogeneous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, size-dependent growth rates, and surface crystallization are also included. Model predictions are compared with experimental measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the peak in the steady-state rate is demonstrated.

Kelton, K. F.↗