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 721 records · Page 40

A comparison of mapped and measured total ionospheric electron content using global positioning system and beacon satellite observations

Total ionospheric electron contents (TEC) were measured by global positioning system (GPS) dual-frequency receivers developed by the Jet Propulsion Laboratory. The measurements included P-code (precise ranging code) and carrier phase data for six GPS satellites during multiple five-hour observing sessions. A set of these GPS TEC measurements were mapped from the GPS lines of sight to the line of sight of a Faraday beacon satellite by statistically fitting the TEC data to a simple model of the ionosphere. The mapped GPS TEC values were compared with the Faraday rotation measurements. Because GPS transmitter offsets are different for each satellite and because some GPS receiver offsets were uncalibrated, the sums of the satellite and receiver offsets were estimated simultaneously with the TEC in a least squares procedure. The accuracy of this estimation procedure is evaluated indicating that the error of the GPS-determined line of sight TEC can be at or below 1 x 10 to the 16th el/sq cm. Consequently, the current level of accuracy is comparable to the Faraday rotation technique; however, GPS provides superior sky coverage.

Lanyi, Gabor E.↗

Stellar Occultation Probe of Triton's Atmosphere

The goals of this research were (i) to better characterize Triton's atmospheric structure by probing a region not well investigated by Voyager and (ii) to begin acquiring baseline data for an investigation of the time evolution of the atmosphere which will set limits on the thermal conductivity of the surface and the total mass of N2 in the atmosphere. Our approach was to use observations (with the Kuiper Airborne Observatory) of a stellar occultation by Triton that was predicted to occur on 1993 July 10. As described in the attached reprint, we achieved these objectives through observation of this occultation and a subsequent one with the KAO in 1995. We found new results about Triton's atmospheric structure from the analysis of the two occultations observed with the KAO and ground-based data. These stellar occultation observations made both in the visible and infrared, have good spatial coverage of Triton including the first Triton central-flash observations, and are the first data to probe the 20-100 km altitude level on Triton. The small-planet light curve model of Elliot and Young (AJ 103, 991-1015) 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 homogeneous 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 +/- 1 K and the atmospheric pressure at 1400 km radius (approximately 50 km altitude) is 1.4 +/- 0.1 microbar. Both of these are not consistent with a model based on Voyager UVS and RSS observations in 1989 (Strobel et al, Icarus 120, 266-289). 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.

Elliot, James L.↗

The Radio luminosity Function of Radio-Loud Quasars from the 7C Redshift Survey

We present a complete sample of 24 radio-loud quasars (RLQs) from the new 7C Redshift Survey. Every quasar with a low-frequency (151 MHz) radio flux-density S(sub 151) > 0.5 Jy in two regions of the sky covering 0.013 sr is included; 23 of these have sufficient extended flux to meet the selection criteria, 18 of these have steep radio spectra (hereafter denoted as SSQs). The key advantage of this sample over most samples of RLQs is the lack of an optical magnitude limit. By combining the 7C and 3CRR samples, we have investigated the properties of RLQs as a function of redshift z and radio luminosity L(sub 151). We derive the radio luminosity function (RLF) of RLQs and find that the data are well fitted by a single power-law with slope alpha(sub 1) = 1.9 +/- 0.1 (for H(sub 0) = 50 km/s.Mpc, OMEGA(sub M) = 1, OMEGA(sub DELTA) = 0). We find that there must be a break in the RLQ RLF at log(sub 10)(L(sub 151)/W Hz.sr) approximately < or = 27, in order for the models to be consistent with the 7C and 6C source counts. The z-dependence of the RLF follows a one-tailed gaussian which peaks at z = 1.7 +/- 0.2. We find no evidence for a decline in the co-moving space density of RLQs at higher redshifts. A positive correlation between the radio and optical luminosities of SSQs is observed, confirming a result of Serjeant. We are able to rule out this correlation being due to selection effects or biases in our combined sample. The radio-optical correlation and best-fit model RLF enable us to estimate the distribution of optical magnitudes of quasars in samples selected at low radio frequencies, We con- clude that for samples with S(sub 151) approximately < or = 1 Jy one must use optical data significantly deeper than the POSS-I limit (R approximately equal 20), in order to avoid severe incompleteness.

Willott, Chris J.↗

Model of visual contrast gain control and pattern masking

We have implemented a model of contrast gain and control in human vision that incorporates a number of key features, including a contrast sensitivity function, multiple oriented bandpass channels, accelerating nonlinearities, and a devisive inhibitory gain control pool. The parameters of this model have been optimized through a fit to the recent data that describe masking of a Gabor function by cosine and Gabor masks [J. M. Foley, "Human luminance pattern mechanisms: masking experiments require a new model," J. Opt. Soc. Am. A 11, 1710 (1994)]. The model achieves a good fit to the data. We also demonstrate how the concept of recruitment may accommodate a variant of this model in which excitatory and inhibitory paths have a common accelerating nonlinearity, but which include multiple channels tuned to different levels of contrast.

NASA Discipline Space Human Factors↗

Dynamic sea surface topography, gravity and improved orbit accuracies from the direct evaluation of SEASAT altimeter data

A method for the simultaneous solution of dynamic ocean topography, gravity and orbits using satellite altimeter data is described. A GEM-T1 based gravitational model called PGS-3337 that incorporates Seasat altimetry, surface gravimetry and satellite tracking data has been determined complete to degree and order 50. The altimeter data is utilized as a dynamic observation of the satellite's height above the sea surface with a degree 10 model of dynamic topography being recovered simultaneously with the orbit parameters, gravity and tidal terms in this model. PGS-3337 has a geoid uncertainty of 60 cm root-mean-square (RMS) globally, with the uncertainty over the altimeter tracked ocean being in the 25 cm range. Doppler determined orbits for Seasat, show large improvements, with the sub-30 cm radial accuracies being achieved. When altimeter data is used in orbit determination, radial orbital accuracies of 20 cm are achieved. The RMS of fit to the altimeter data directly gives 30 cm fits for Seasat when using PGS-3337 and its geoid and dynamic topography model. This performance level is two to three times better than that achieved with earlier Goddard earth models (GEM) using the dynamic topography from long-term oceanographic averages. The recovered dynamic topography reveals the global long wavelength circulation of the oceans with a resolution of 1500 km. The power in the dynamic topography recovery is now found to be closer to that of oceanographic studies than for previous satellite solutions. This is attributed primarily to the improved modeling of the geoid which has occurred. Study of the altimeter residuals reveals regions where tidal models are poor and sea state effects are major limitations.

Marsh, J. G.↗

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗

Modelling atmospheric scatterers using spacecraft observations

Voyager images of Triton indicate considerable spatial variability in the concentration of at least two different scattering components in the atmosphere. Data from high phase angle limb scans were fit to Mie scattering models to derive mean particle sizes, number densities, and vertical extent for both types of scattering material at ten different locations between 15 deg S and 70 deg S. These fits reveal a thin haze at latitudes equatorward of 25-30 deg S. The imaging data can be fit reasonably well by both conservatively scattering and absorbing hazes with particle sizes near 0.18 micron and optical depths of order 0.001-0.01. Rayleigh scattering haze fits the imaging data somewhat less well, and can be totally ruled out by combining the imaging and UVS measurements. At high southern latitudes, Triton displays clouds below an altitude of approximately 8 km, as well as the haze at higher altitudes. The clouds have particle sizes which may range from 0.7-2.0 microns, or may be near 0.25 micron. The atmospheric optical depth poleward of 30 deg S must be generally greater than 0.1, but need not be more than 0.3. Horizontal inhomogeneities are quite noticeable, especially at longitudes east of (i.e., higher than) 180 deg.

Rages, Kathy A.↗

Visual Data Analysis for Satellites

The Visual Data Analysis Package is a collection of programs and scripts that facilitate visual analysis of data available from NASA and NOAA satellites, as well as dropsonde, buoy, and conventional in-situ observations. The package features utilities for data extraction, data quality control, statistical analysis, and data visualization. The Hierarchical Data Format (HDF) satellite data extraction routines from NASA's Jet Propulsion Laboratory were customized for specific spatial coverage and file input/output. Statistical analysis includes the calculation of the relative error, the absolute error, and the root mean square error. Other capabilities include curve fitting through the data points to fill in missing data points between satellite passes or where clouds obscure satellite data. For data visualization, the software provides customizable Generic Mapping Tool (GMT) scripts to generate difference maps, scatter plots, line plots, vector plots, histograms, timeseries, and color fill images.

Lau, Yee↗

Remote and Ground Truth Spectral Measurement Comparisons

FORMOSAT III are a set of six research satellites from Taiwan that were launched in April 2006. The satellites are in 800 km, 71 degree inclination orbits and separated by 24 degrees in ascending node. Ground truth spectral measurements were taken of outer surface materials on FORMOSAT III. From those measurements, a computer model was built to predict the spectral reflectance, which included phase angle and orientation of the spacecraft relative to the observer. However, materials exposed to the space environment have exhibited spectral changes including a darkening and a reddening of the spectra. This reddening was seen as an increase in slope of the reflectance as the wavelength increases. Therefore, the model of pristine materials was augmented to include the space weathering effects. Remote data were collected on two of the six FORMOSAT satellites using the 1.6 meter telescope at AMOS (Air Force Maui Optical and Supercomputing) site with the Spica spectrometer. Due to the separation in ascending node, observations were made on whichever one of the six satellites was visible on that specific night. Three nights of data were collected using the red (6000 9500 angstroms) filter and two nights of data were collected using the blue (3200 -6600 angstroms) filter. A comparison of the data showed a good match to the pristine model for the blue filter region. The absorption feature near 5500 angstroms due to the copper colored Kapton multi-layer insulation (MLI) was very apparent in the remote samples and a good fit to the data was seen in both satellites observed. The features in the red filter regime agreed with the pristine model up through 7000 angstroms where the reddening begins and the slope of the remote sample increases. A comparison of the two satellites showed similar features in the red and blue filter regions, i.e. the satellites were aging at the same rate. A comparison of the pristine model to the first month of remote measurements showed the amount by which the satellite had reddened. This information was used in the space weathered model to predict the second month s reflectance measurement. The results showed that the reddening did not increase from the first month of observations to the second. A third month of data will be necessary to determine if the reddening occurs instantly and then stabilizes with time or if the effect continues to increase with time.

Abercromby, Kira Jorgensen↗

Magnetization of the Lunar Crust

Magnetic fields measured by the satellite Lunar Prospector show large scale features resulting from remanently magnetized crust. Vector data synthesized at satellite altitude from a spherical harmonic model of the lunar crustal field, and the radial component of the magnetometer data, have been used to produce spatially continuous global magnetization models for the lunar crust. The magnetization is expressed in terms of localized basis functions, with a magnetization solution selected having the smallest root-mean square magnetization for a given fit to the data, controlled by a damping parameter. Suites of magnetization models for layers with thicknesses between 10 and 50 km are able to reproduce much of the input data, with global misfits of less than 0.5 nT (within the uncertainties of the data), and some surface field estimates. The magnetization distributions show robust magnitudes for a range of model thicknesses and damping parameters, however the magnetization direction is unconstrained. These global models suggest that magnetized sources of the lunar crust can be represented by a 30 km thick magnetized layer. Average magnetization values in magnetized regions are 30-40 mA/m, similar to the measured magnetizations of the Apollo samples and significantly weaker than crustal magnetizations for Mars and the Earth. These are the first global magnetization models for the Moon, providing lower bounds on the magnitude of lunar crustal magnetization in the absence of multiple sample returns, and can be used to predict the crustal contribution to the lunar magnetic field at a particular location.

Lunar↗

The relative line strength and intensity of the N II 2143 doublet

The doublet emission from N II at 2139.7 A and 2143.6 A was observed by a 1.4-m scanning spectrometer with 3.1 A resolution in the daytime, high-altitude thermosphere during moderate levels of solar activity. The spectrometer viewed the earth's limb 5 deg below the local horizontal to give a nominal tangent height of 152 km. Both sub band heads of the nitric oxide gamma band system were resolved in the data at the resolution used. The emission features from N II are clearly evident on the short wavelength shoulder of the (1, 0) band. Synthetic profiles of the (1, 0) gamma band and the (0, 3) delta band of nitric oxide were fitted to the data using a chi-square analysis. These contributions were removed from the data leaving a residual emission, considered to be the N II doublet. A chi-square minimization of the data relative to a synthetic intensity profile was done. The minimum was for a line strength ratio between the 2139 A and 2143 A lines of 0.58 + or - 0.08. The mean solar EUV flux deduced from the intensity of the N II emission in this experiment is lower than other reported observations, consistent with a lower solar activity level.

Bucsela, Eric J.↗

Updating High-Order Aeroservoelastic Models from Low-Order System Identification Results

Estimating aircraft models from test data requires several simplifying assumptions that introduce biases into the parameter estimates. In this paper, these biases are defined and a method for estimating the biases is discussed. Having an estimate of the bias allows the parameters estimated from test to be integrated into a high-order model. The Integrated Adaptive Wing Technology Maturation (IAWTM) wind tunnel model is discussed and the bias is demonstrated for one of the testing configurations. The methodology was able to estimate these biases and apply corrections to high-order aeroelastic models to improve the fit to test data. The consideration of the biases allows more meaningful comparisons and avoids the erroneous differences between pretest predictions and the fitted model.

Jeffrey Ouellette↗

Updating High-Order Aeroservoelastic Models from Low-Order System Identification Results

Estimating aircraft models from test data requires several simplifying assumptions that introduce biases into the parameter estimates. In this paper, these biases are defined and a method for estimating the biases is discussed. Having an estimate of the bias allows the parameters estimated from test to be integrated into a high-order model. The Integrated Adaptive Wing Technology Maturation (IAWTM) wind tunnel model is discussed and the bias is demonstrated for one of the testing configurations. The methodology was able to estimate these biases and apply corrections to high-order aeroelastic models to improve the fit to test data. The consideration of the biases allows more meaningful comparisons and avoids the erroneous differences between pretest predictions and the fitted model.

Jeffrey Ouellette↗

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic (LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖, 𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, there are some instances where the fitting needs to be done programmatically. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

A method for determining crack opening load from load-displacement data

The fatigue crack opening load is determined as the tangent point on the ascending load-displacement data between the curved portion and the upper linear region. A model for the 'unzipping' behavior of the crack indicates that the curved portion of the load-deflection curve is second order. The opening load is determined by a nonlinear, least squares fit of the data to the model, which optimally locates the tangent point of the two curves. The method provided consistent results for determining opening load P(op) for 7475-T731 aluminum using data from a crack tip opening gage.

Carman, C. Davis↗

Remote and Ground Truth Spectral Measurement Comparisons of FORMOSAT III

FORMOSAT III are a set of six research satellites from Taiwan that were launched in April 2006. The satellites are in 800 km, 71 degree inclination orbits and separated by 24 degrees in ascending node. Laboratory spectral measurements were taken of outer surface materials on FORMOSAT III. From those measurements, a computer model was built to predict the spectral reflectance accounting for both solar phase angle and orientation of the spacecraft relative to the observer. However, materials exposed to the space environment have exhibited spectral changes including a darkening and a "reddening" of the spectra. This "reddening" is characterized by an increase in slope of the reflectance as the wavelength increases. Therefore, the model of pre-flight materials was augmented to include the presumed causative agent: space weathering effects. Remote data were collected on two of the six FORMOSAT satellites using the 1.6 meter telescope at the AMOS (Air Force Maui Optical and Supercomputing) site with the Spica spectrometer. Due to the separation in ascending node, observations were acquired of whichever one of the six satellites was visible on that specific night. Three nights of data were collected using the red (6000 - 9500 angstroms) filter and five nights of data were collected using the blue (3200 - 6600 angstroms) filter. A comparison of the data showed a good match to the pre-flight models for the blue filter region. The absorption feature near 5500 angstroms due to the copper colored Kapton multi-layer insulation (MLI) was very apparent in the remote samples and a good fit to the data was seen in all satellites observed. The features in the red filter regime agreed with the pre-flight model up through 7000 angstroms where the reddening begins and the slope of the remote sample increases. A comparison of the satellites showed similar features in the red and blue filter regions, i.e. the satellite surfaces were aging at the same rate. A comparison of the pre-flight model to the first month of remote measurements showed the amount by which the satellite had reddened. The second month of data observed a satellite at a higher altitude and was therefore, not compared to the first month. A third month of data was collected but of satellites at the lower altitude regime and can only be compared to the first month. One cause of the reddening that was ruled out in early papers was a possible calibration issue.

Abercromby, Kira Jorgensen↗

Recovery of atmospheric refractivity profiles from simulated satellite-to-satellite tracking data

Techniques for recovering atmospheric refractivity profiles from simulated satellite-to-satellite tracking data are documented. Examples are given using the geometric configuration of the ATS-6/NIMBUS-6 Tracking Experiment. The underlying refractivity model for the lower atmosphere has the spherically symmetric form N = exp P(s) where P(s) is a polynomial in the normalized height s. For the simulation used, the Herglotz-Wiechert technique recovered values which were 0.4% and 40% different from the input values at the surface and at a height of 33 kilometers, respectively. Using the same input data, the model fitting technique recovered refractivity values 0.05% and 1% different from the input values at the surface and at a height of 50 kilometers, respectively. It is also shown that if ionospheric and water vapor effects can be properly modelled or effectively removed from the data, pressure and temperature distributions can be obtained.

Murray, C. W., Jr.↗