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 1,153 records · Page 64

Simulation of EAS properties on the basis of high energy interaction model deduced from the accelerator data

Calculations of extensive air showers in atmosphere were performed using formulae describing p-p and p-air nucleus interactions. The formulae fitted to the accelerator data were extrapolated taking the same trend up to 10 to the 16 eV. Above that energy it was assumed that the degree of scaling violating/alpha-parameter/ is saturating or even decreasing. The latter assumption follows from earlier work where it was found that without this restriction shower maxima at the highest energies are located too high in the atmosphere. Results of calculations have been compared with experimental data. The comparison was made separately for the curves obtained from the so called equal intensity cuts and for the Cerenkov data.

Kubiak, G.↗

The response of the thermospheric nitric oxide to an auroral storm. I - Low and middle latitudes

The response of low- and mid-latitude NO to the September 19, 1984 magnetic storm is studied. In particular, the effect of variations in the temperature and composition of the lower thermosphere on NO is investigated. Much of the observed NO response is adequately explained by including these variations in a photochemical model. In particular, Joule heating at auroral latitudes can cause temperature enhancements at nonauroral latitudes, which in turn can lead directly to an increase in the NO. The NO is most directly sensitive to heating at altitude above 120 km. Below 120 km, NO increased through downward diffusion from above. Sensitivity studies were performed to investigate the effect that uncertainties in the NO chemical scheme and in the state of the neutral atmosphere can have on the comparison of the model to the data. It is shown that the best fit to the midlatitude American orbit data is for O/O2 ratio twice that in the original TGCM simulation and for an N(2D) quenching value of 5 x 10 to the -13th.

Siskind, D. E.↗

Triaxial Probe Magnetic Data Analysis

The Triaxial Magnetic Moment Analysis software uses measured magnetic field test data to compute dipole and quadrupole moment information from a hardware element. It is used to support JPL projects needing magnetic control and an understanding of the spacecraft-generated magnetic fields. Evaluation of the magnetic moment of an object consists of three steps: acquisition, conditioning, and analysis. This version of existing software was extensively rewritten for easier data acquisition, data analysis, and report presentation, including immediate feedback to the test operator during data acquisition. While prior JPL computer codes provided the same data content, this program has a better graphic display including original data overlaid with reconstructed results to show goodness of fit accuracy and better appearance of the report graphic page. Data are acquired using three magnetometers and two rotations of the device under test. A clean acquisition user interface presents required numeric data and graphic summaries, and the analysis module yields the best fit (least squares) for the magnetic dipole and/or quadrupole moment of a device. The acquisition module allows the user to record multiple data sets, selecting the best data to analyze, and is repeated three times for each of the z-axial and y-axial rotations. In this update, the y-axial rotation starting position has been changed to an option, allowing either the x- or z-axis to point towards the magnetometer. The code has been rewritten to use three simultaneous axes of magnetic data (three probes), now using two "rotations" of the device under test rather than the previous three rotations, thus reducing handling activities on the device under test. The present version of the software gathers data in one-degree increments, which permits much better accuracy of the fit ted data than the coarser data acquisition of the prior software. The data-conditioning module provides a clean data set for the analysis module. For multiple measurements at a given degree, the first measurement is used. For omitted measurements, the missing field is estimated by linear interpolation between the two nearest measurements. The analysis module was rewritten for the dual rotation, triaxial probe measurement process and now has better moment estimation accuracy, based on the finer one degree of data acquisition resolution. The magnetic moments thus computed are used as an input to summarize the total spacecraft field.

Shultz, Kimberly↗

An interstellar cloud density from Copernicus observations of CO in the spectrum of Zeta Ophiuchi

Interstellar CO absorption bands in Copernicus spectra of Zeta Oph have been studied. Absorption profiles, computed under the assumption that excitation is due to collisions with H2 molecules and interaction with the 3-K background radiation field, were fitted to the reduced data of nine bands. When a gas kinetic temperature of 56 K is assumed, the best-fit condition implies a hydrogen-nucleus density of 120 per cu cm, a CO column density of 1.2 by 10 to the 15th power per sq cm, and a radial-velocity dispersion of 0.9 km/s. The relevance of these results to existing ideas concerning the Zeta Oph interstellar clouds is discussed. It is suggested that the strongest interstellar component is not circumstellar in origin but is instead part of a supernova remnant. Simple calculations are made to establish the plausibility of the supernova-remnant identification. This suggestion is also supported by Heiles's (1976) 21-cm pictures.

Smith, A. M.↗

Tests of General Relativity With the Binary Black Hole Signals from the LIGO-Virgo Catalog GWTC-1

The detection of gravitational waves by Advanced LIGO and Advanced Virgo provides an opportunity to test general relativity in a regime that is inaccessible to traditional astronomical observations and laboratory tests. We present four tests of the consistency of the data with binary black hole gravitational waveforms predicted by general relativity. One test subtracts the best-fit waveform from the data and checks the consistency of the residual with detector noise. The second test checks the consistency of the low- and high-frequency parts of the observed signals. The third test checks that phenomenological deviations introduced in the waveform model (including in the post-Newtonian coefficients) are consistent with 0. The fourth test constrains modifications to the propagation of gravitational waves due to a modified dispersion relation, including that from a massive graviton. We present results both for individual events and also results obtained by combining together particularly strong events from the first and second observing runs of Advanced LIGO and Advanced Virgo, as collected in the catalog GWTC-1. We do not find any inconsistency of the data with the predictions of general relativity and improve our previously presented combined constraints by factors of 1.1 to 2.5. In particular, we bound the mass of the graviton to be 𝑚𝑔≤4.7×10 −23 eV/𝑐 2 (90% credible level), an improvement of a factor of 1.6 over our previously presented results. Additionally, we check that the four gravitational-wave events published for the first time in GWTC-1 do not lead to stronger constraints on alternative polarizations than those published previously.

B P Abbott↗

Curve fits of predicted inviscid stagnation-point radiative heating rates, cooling factors, and shock standoff distances for hyperbolic earth entry

Curve-fit formulas are presented for the stagnation-point radiative heating rate, cooling factor, and shock standoff distance for inviscid flow over blunt bodies at conditions corresponding to high-speed earth entry. The data which were curve fitted were calculated by using a technique which utilizes a one-strip integral method and a detailed nongray radiation model to generate a radiatively coupled flow-field solution for air in chemical and local thermodynamic equilibrium. The range of free-stream parameters considered were altitudes from about 55 to 70 km and velocities from about 11 to 16 km.sec. Spherical bodies with nose radii from 30 to 450 cm and elliptical bodies with major-to-minor axis ratios of 2, 4, and 6 were treated. Powerlaw formulas are proposed and a least-squares logarithmic fit is used to evaluate the constants. It is shown that the data can be described in this manner with an average deviation of about 3 percent (or less) and a maximum deviation of about 10 percent (or less). The curve-fit formulas provide an effective and economic means for making preliminary design studies for situations involving high-speed earth entry.

Suttles, J. T.↗

An equation of state for oxygen and nitrogen 2

Preliminary equations of state for oxygen and nitrogen are presented which are the result of least squares fitting of the available P-density-T data, heat capacities at constant volume, and data to establish the criteria for phase equilibrium, simultaneously. The equation for nitrogen is applicable to 10,000 atmospheres at temperatures from 70 K to 1300 K, and the equation for oxygen to 350 atmospheres for temperatures from 56 K to 323 K. The development of the functional form of the equation of state, data weighting methods, and the techniques employed for least squares fitting of related data are discussed. Comparisons of the equation of state developed with selected values of experimental P-density-T data, and heat capacity values are included for both fluids.

Stewart, R. B.↗

Implementation of thermal residual stresses in the analysis of fiber bridged matrix crack growth in titanium matrix composites

In this research, thermal residual stresses were incorporated in an analysis of fiber-bridged matrix cracks in unidirectional and cross-ply titanium matrix composites (TMC) containing center holes or center notches. Two TMC were investigated, namely, SCS-6/Timelal-21S laminates. Experimentally, matrix crack initiation and growth were monitored during tension-tension fatigue tests conducted at room temperature and at an elevated temperature of 200 C. Analytically, thermal residual stresses were included in a fiber bridging (FB) model. The local R-ratio and stress-intensity factor in the matrix due to thermal and mechanical loadings were calculated and used to evaluate the matrix crack growth behavior in the two materials studied. The frictional shear stress term, tau, assumed in this model was used as a curve-fitting parameter to matrix crack growth data. The scatter band in the values of tau used to fit the matrix crack growth data was significantly reduced when thermal residual stresses were included in the fiber bridging analysis. For a given material system, lay-up and temperature, a single value of tau was sufficient to analyze the crack growth data. It was revealed in this study that thermal residual stresses are an important factor overlooked in the original FB models.

Bakuckas, John G., Jr.↗

Version 2 Ozone Monitoring Instrument SO2 Product (OMSO2 V2): New Anthropogenic SO2 Vertical Column Density Dataset

The Ozone Monitoring Instrument (OMI) has been providing global observations of SO2 pollution since 2004. Here we introduce the new anthropogenic SO2 vertical column density (VCD) dataset in the version 2 OMI SO2 product (OMSO2 V2). As with the previous version (OMSO2 V1.3), the new dataset is generated with an algorithm based on principal component analysis of OMI radiances, but features several updates. The most important among those is the use of expanded lookup tables and model a priori profiles to estimate SO2 Jacobians for individual OMI pixels, in order to better characterize pixel-to-pixel variations in SO2 sensitivity, including over snow and ice. Additionally, new data screening and spectral fitting schemes have been implemented to improve the quality of the spectral fit. As compared with the planetary boundary layer SO2 dataset in OMSO2 V1.3, the new dataset has substantially better data quality, especially over areas that are relatively clean or affected by the south Atlantic anomaly. The updated retrievals over snow/ice yield more realistic seasonal changes in SO2 at high latitudes and offer enhanced sensitivity to sources during wintertime. An error analysis has been conducted to assess uncertainties in SO2 VCDs from both the spectral fit and Jacobian calculations. The uncertainties from spectral fitting are reflected in SO2 slant column densities (SCDs) and largely depend on the signal-to-noise ratio of the measured radiances, as implied by the generally smaller SCD uncertainties over clouds or for smaller solar zenith angles. The SCD uncertainties for individual pixels are estimated to be~0.15-0.3 DU (Dobson Units) between ~40°S and ~40°N and to be~0.2-0.5 DU at higher latitudes. The uncertainties from the Jacobians are approximately ~50-100% over polluted areas, and primarily attributed to errors in SO2 a priori profiles and cloud pressures, as well as the lack of explicit treatment for aerosols. Finally, the daily mean and median SCDs over the presumably SO2-free equatorial East Pacific have increased by only~0.0035 DU and ~0.003 DU respectively over the entire 15-year OMI record; while the standard deviation of SCDs has grown by only~0.02 DU or ~10%. Such remarkable long-term stability makes the new dataset particularly suitable for detecting regional changes in SO2 pollution.

OMI, SO2, Remote Sensing↗

Definition of the Flexible Image Transport System (FITS), Version 3.0

The Flexible Image Transport System (FITS) has been used by astronomers for over 30 years as a data interchange and archiving format; FITS files are now handled by a wide range of astronomical software packages. Since the FITS format definition document (the "standard") was last printed in this journal in 2001, several new features have been developed and standardized, notably support for 64-bit integers in images and tables, variable-length arrays in tables, and new world coordinate system conventions which provide a mapping from an element in a data array to a physical coordinate on the sky or within a spectrum. The FITS Working Group of the International Astronomical Union has therefore produced this new Version 3.0 of the FITS standard, which is provided here in its entirety. In addition to describing the new features in FITS, numerous editorial changes were made to the previous version to clarify and reorganize many of the sections. Also included are some appendices which are not formally part of the standard. The FITS standard is likely to undergo further evolution, in which case the latest version may be found on the FITS Support Office Web site at http://fits.gsfc.nasa.gov/, which also provides many links to FITS-related resources.

Pence, W. D.↗

Crab Nebula observations - 0.2-10 MeV

Observations of the total emission from the Crab Nebula and also of the pulsed component were made over the 0.2 to 10-MeV range during three balloon flights in 1971 with an actively-collimated NaI scintillator. The total emission flux was positively observed over the entire interval. The observed spectrum to 1 MeV agrees with an extrapolation of the E to the -2.2 power law, which fits lower-energy data. The observations above 1 MeV are factors of 3 and 20 above this law and are better fit with a spectral index of 0.8. Confidence levels are 3 sigma or better for each half-decade band. The three observations are consistent with a constant flux level. The NP 0532 flux, detected during one flight only (August 8) between 0.2 and 0.38 MeV, agrees with the exponential power law spectrum already determined from other observations. The possibility of a rapidly rising pulsed emission fraction over the 0.1- to 1-MeV interval is excluded by this observation.

Gruber, D. E.↗

Fractal analysis: A new remote sensing tool for lava flows

Many important quantitative parameters have been developed that relate to the rheology and eruption and emplacement mechanics of lavas. This research centers on developing additional, unique parameters, namely the fractal properties of lava flows, to add to this matrix of properties. There are several methods of calculating the fractal dimension of a lava flow margin. We use the 'structured walk' or 'divider' method. In this method, we measure the length of a given lava flow margin by walking rods of different lengths along the margin. Since smaller rod lengths transverse more smaller-scaled features in the flow margin, the apparent length of the flow outline will increase as the length of the measuring rod decreases. By plotting the apparent length of the flow outline as a function of the length of the measuring rod on a log-log plot, fractal behavior can be determined. A linear trend on a log-log plot indicates that the data are fractal. The fractal dimension can then be calculated from the slope of the linear least squares fit line to the data. We use this 'structured walk' method to calculate the fractal dimension of many lava flows using a wide range of rod lengths, from 1/8 to 16 meters, in field studies of the Hawaiian islands. We also use this method to calculate fractal dimensions from aerial photographs of lava flows, using lengths ranging from 20 meters to over 2 kilometers. Finally, we applied this method to orbital images of extraterrestrial lava flows on Venus, Mars, and the Moon, using rod lengths up to 60 kilometers.

Bruno, B. C.↗

Radiometer Calibrations: Saving Time by Automating the Gathering and Analysis Procedures

Mr. Abtahi custom-designs radiometers for Mr. Hook's research group. Inherently, when the radiometers report the temperature of arbitrary surfaces, the results are affected by errors in accuracy. This problem can be reduced if the errors can be accounted for in a polynomial. This is achieved by pointing the radiometer at a constant-temperature surface. We have been using a Hartford Scientific WaterBath. The measurements from the radiometer are collected at many different temperatures and compared to the measurements made by a Hartford Chubb thermometer with a four-decimal point resolution. The data is analyzed and fit to a fifth-order polynomial. This formula is then uploaded into the radiometer software, enabling accurate data gathering. Traditionally, Mr. Abtahi has done this by hand, spending several hours of his time setting the temperature, waiting for stabilization, taking measurements, and then repeating for other temperatures. My program, written in the Python language, has enabled the data gathering and analysis process to be handed off to a less-senior member of the team. Simply by entering several initial settings, the program will simultaneously control all three instruments and organize the data suitable for computer analyses, thus giving the desired fifth-order polynomial. This will save time, allow for a more complete calibration data set, and allow for base calibrations to be developed. The program is expandable to simultaneously take any type of measurement from up to nine distinct instruments.

Sadino, Jeffrey L.↗

Regression Analysis of Long-Term Profile Ozone Data Set from BUV Instruments

We have produced a profile merged ozone data set (MOD) based on the SBUV/SBUV2 series of nadir-viewing satellite backscatter instruments, covering the period from November 1978 - December 2003. In 2004, data from the Nimbus 7 SBUV and NOAA 9, ll, and 16 SBUV/2 instruments were reprocessed using the Version 8 (V8) algorithm and most recent calibrations. More recently, data from the Nimbus 4 BUT instrument, which was operational from 1970 - 1977, were also reprocessed using the V8 algorithm. As part of the V8 profile calibration, the Nimbus 7 and NOAA 9 (1993-1997 only) instrument calibrations have been adjusted to match the NOAA 11 calibration, which was established based on comparisons with SSBUV shuttle flight data. Differences between NOAA 11, Nimbus 7 and NOAA 9 profile zonal means are within plus or minus 5% at all levels when averaged over the respective periods of data overlap. NOAA 16 SBUV/2 data have insufficient overlap with NOAA 11, so its calibration is based on pre-flight information. Mean differences over 4 months of overlap are within plus or minus 7%. Given the level of agreement between the data sets, we simply average the ozone values during periods of instrument overlap to produce the MOD profile data set. Initial comparisons of coincident matches of N4 BUV and Arosa Umkehr data show mean differences of 0.5 (0.5)% at 30km; 7.5 (0.5)% at 35 km; and 11 (0.7)% at 40 km, where the number in parentheses is the standard error of the mean. In this study, we use the MOD profile data set (1978-2003) to estimate the change in profile ozone due to changing stratospheric chlorine levels. We use a standard linear regression model with proxies for the seasonal cycle, solar cycle, QBO, and ozone trend. To account for the non-linearity of stratospheric chlorine levels since the late 1990s, we use a time series of Effective Chlorine, defined as the global average of Chlorine + 50 * Bromine at 1 hPa, as the trend proxy. The Effective Chlorine data are taken from the 3-D Goddard CTM. We will show the latest trend results using this statistical model. In addition, the Nimbus 4 BUV data offer an opportunity to test the physical properties of our statistical model. From ground-based comparisons we will establish an uncertainty range for the Nimbus 4 data. We then extrapolate our statistical model fit backwards in time and compare to the Nimbus 4 data. We compare the characteristics of the residual, defined as the difference between the data and statistical regression fit, during the Nimbus 4 time period and the 1978-2003 period over which the statistical model coefficients were estimated, and present these results.

Stolarski, Richard S.↗

Frequency domain state-space system identification

An algorithm for identifying state-space models from frequency response data of linear systems is presented. A matrix-fraction description of the transfer function is employed to curve-fit the frequency response data, using the least-squares method. The parameters of the matrix-fraction representation are then used to construct the Markov parameters of the system. Finally, state-space models are obtained through the Eigensystem Realization Algorithm using Markov parameters. The main advantage of this approach is that the curve-fitting and the Markov parameter construction are linear problems which avoid the difficulties of nonlinear optimization of other approaches. Another advantage is that it avoids windowing distortions associated with other frequency domain methods.

Chen, Chung-Wen↗

Use of SPOT and ERS-1 SAR data to study the tectonic and climatic history of arid regions

In order to separate the effects of the different tectonic and climatic processes on the shapes of desert piedmonts, a modified conic equation was fitted to digital topographic data for individual alluvial fans in Death Valley (California, U.S.). The topographic data were obtained from a SPOT panchromatic stereo pair and from the airborne interferometric SAR (Synthetic Aperture Radar) (TOPSAR). The conic fit allows parameters for the epex position, slope, and radial curvature to be compared with unit age, uplift rate, and climatic conditions. Preliminary results indicate that slope flattens with age and radial curvature is concave up, but decreases with age. Work is continuing on correlation of fit residuals and apex position with fan unit age. This information will help in the determination of tectonic uplift rates and the climatic history of the western U.S. ERS-1 SAR images were used to study an area of western China where a large strike slip fault crosses a series of alluvial fans and stream valleys. Previous analysis of SPOT panchromatic images of the area shows that offsets fans and streams can be recognized. Measurement of the rate of motion of this fault will help in the overall model of deformation of the Asian tectonic plate in response to the collision of the Indian plate.

Farr, Tom G.↗

Photographic photometry with Iris diaphragm photometers

A general method is presented for solving problems encountered in the analysis of Iris diaphragm photometer (IDP) data. The method is used to derive the general shape of the calibration curve, allowing both a more accurate fit to the IDP data for comparison stars and extrapolation to magnitude ranges for which no comparison stars are measured. The profile of starlight incident and the characteristic curve of the plate are both assumed and then used to derive the profile of the star image. An IDP reading is then determined for each star image. A procedure for correcting the effects of a nonconstant background fog level on the plate is also demonstrated. Additional applications of the method are made in the appendix to determine the relation between the radius of a photographic star image and the star's magnitude, and to predict the IDP reading of the 'point of optimum density'.

Schaefer, B. E.↗

Limits on the radius and a possible atmosphere of Charon from its 1980 stellar occultation

Walker's (1980) stellar occultation data for Charon are presently fit by a model which encompasses the possibility of differential refraction by an atmosphere, followed by a sudden occultation behind Charon's limb. The 601.5-km Charon radius lower limit thus obtained may serve as a constraint in models of the mutual event data; while the model fits considered support a Charonian atmosphere of indeterminate composition, time resolution is insufficient for certainty and the data may be interpretable as indications of either a slight extinction near Charon or an entirely unidentified and unassociated effect.

Elliot, J. L.↗