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 307 records · Page 17

GPS measurements of strain accumulation across the Imperial Valley, California: 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 +/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 +/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Global Positioning System measurements of strain accumulation across the Imperial Valley, California - 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 =/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 =/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Computer codes for the evaluation of thermodynamic properties, transport properties, and equilibrium constants of an 11-species air model

The computer codes developed provide data to 30000 K for the thermodynamic and transport properties of individual species and reaction rates for the prominent reactions occurring in an 11-species nonequilibrium air model. These properties and the reaction-rate data are computed through the use of curve-fit relations which are functions of temperature (and number density for the equilibrium constant). The curve fits were made using the most accurate data believed available. A detailed review and discussion of the sources and accuracy of the curve-fitted data used herein are given in NASA RP 1232.

Thompson, Richard A.↗

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

Infrared Communications for Small Spacecraft: From a Wireless Bus to Cluster Concepts

Nanosatellites operating singly or in clusters are anticipated for future space science missions. To implement this new communications paradigm, we are approaching cluster communications by first developing an infrared (IR) intra-craft wireless bus capability, following initially the MIL-STD-1553B protocol. Benefits of an IR wireless bus are low mass, size, power, and cost, simplicity of implementation, ease of use, minimum EMI, and efficient and reliable data transfer. Our goals are to maximize the reliable link margin in order to afford greater flexibility in receiver placement, which will ease technology insertion. We have developed a concept demonstration using a high-speed visible-band silicon PIN photodiode and a high-efficiency visible LED operating at a data rate up to 4 Mb/sec. In designing an internal IR wireless bus, we have characterized various candidate materials, emitters, and geometries, assuming a single reflection. Thus, we have measured the bidirectional reflectance distribution function (BRDF) for five different materials characteristic of typical spacecraft structures, which range from nearly Lambertian to highly specular. We have fit our data to empirical BRDF functions and modeled the detected irradiance anywhere in the plane of incidence for a divergent (LED) emitter. We have also determined the angular limits on the link geometry to remain within the required bit error rate by determining the received signal-to-noise ratio (SNR) for minimum values of irradiance received at the detector.

Webb, Suzanne C.↗

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

Time-average fluxes of heavy ions at the geostationary orbit

Measurements of average proton, helium, carbon, and oxygen fluxes at 6.6 earth radii are reported. The data represent averages obtained on ATS 6 between June 15 and Oct. 3, 1974. The energy range covered was 0.36-1.1 MeV for protons and approximately 1-4 MeV for the heavier ions. The results indicate that above about 1 MeV the heavier ion fluxes dominate over the proton flux on the energy/ion scale. Using two different spectral dependencies to fit the data, the carbon to oxygen concentration ratio for energies above 1 MeV was found to be 0.43 for the power law spectrum and 0.44 for the exponential spectrum. Thus in either case the abundance ratio is consistent with the solar origin of the particles. Similarly, the helium to oxygen concentration ratio is found to be 0.17 for the power law spectrum and 0.28 for the exponential spectrum. This is inconsistent with the solar wind ratio, which is about 10. The results quoted above are based on a very small portion of the distribution function above 1 MeV.

Konradi, A.↗

Application of HCMM data to regional geologic analysis for mineral and energy resource evaluation

Experimentation with several potentially promising techniques led to the selection of a fairly simple scheme for registration of data from the HCMM using an affine transformation. A method based solely on remote sensing was developed to estimate those meteorological effects which are required for thermal inertia mapping. It assumes that the atmospheric fluxes are spatially invariant and that the solar, sky and sensible heat fluxes can be approximated by a simple mathematical form. Coefficients are determined from least squares method by fitting observational data to the thermal model.

Watson, K.↗

ISO Key Project: Exploring the full range of QUASAR/AGN properties

The PIA (PHOT Interactive Analysis) software was upgraded as new releases were made available by VILSPA. We have continued to analyze our data but, given the large number of still outstanding problems with the calibration and analysis (listed below), we remain unable to move forward on our scientific program. We have concentrated on observations with long (256 sec) exposure times to avoid the most extreme detector responsivity drift problems which occur with a change in observed flux level, ie. as one begins to observe a new target. There remain a significant number of problems with analyzing these data including: (1) the default calibration source (FCS) observations early in the mission were too short and affected by strong detector responsivity drifts; (2) the calibration of the FCS sources is not yet well-understood, particularly for chopped observations (which includes most of ours); (3) the detector responsivity drift is not well-understood and models are only now becoming available for fitting chopped data; (4) charged particle hits on the detector cause transient responsivity drifts which need to be corrected; (5) the "flat-field" calibration of the long-wavelength (array) detectors: C1OO, C200 leaves significant residual structure and so needs to be improved;(6) the vignetting correction, which affects detected flux levels in the array detectors, is not yet available; (7) the intra-filter calibrations are not yet available; and (8) the background above 60 microns has a significant gradient which results in spurious positive and negative "detections" in chopped observations. ISO Observation planning, conferences and talks, ground based observing and other grant related activities are also briefly discussed.

Wilkes, B.↗

Application of statistical distribution theory to launch-on-time for space construction logistic support

The ability to launch-on-time and to send payloads into space has progressed dramatically since the days of the earliest missile and space programs. Causes for delay during launch, i.e., unplanned 'holds', are attributable to several sources: weather, range activities, vehicle conditions, human performance, etc. Recent developments in space program, particularly the need for highly reliable logistic support of space construction and the subsequent planned operation of space stations, large unmanned space structures, lunar and Mars bases, and the necessity of providing 'guaranteed' commercial launches have placed increased emphasis on understanding and mastering every aspect of launch vehicle operations. The Center of Space Construction has acquired historical launch vehicle data and is applying these data to the analysis of space launch vehicle logistic support of space construction. This analysis will include development of a better understanding of launch-on-time capability and simulation of required support systems for vehicle assembly and launch which are necessary to support national space program construction schedules. In this paper, the author presents actual launch data on unscheduled 'hold' distributions of various launch vehicles. The data have been supplied by industrial associate companies of the Center for Space Construction. The paper seeks to determine suitable probability models which describe these historical data and that can be used for several purposes such as: inputs to broader simulations of launch vehicle logistic space construction support processes and the determination of which launch operations sources cause the majority of the unscheduled 'holds', and hence to suggest changes which might improve launch-on-time. In particular, the paper investigates the ability of a compound distribution probability model to fit actual data, versus alternative models, and recommends the most productive avenues for future statistical work.

Morgenthaler, George W.↗

An Investigation Into Possible Systematic Effects on Neutron Star Radius Estimates using NICER-like Synthetic Data

Neutron star cores contain the densest matter in the observable universe. The state of this matter is of interest in numerous fields, but laboratory experiments cannot explore this matter. Although the composition of the matter would be of great interest, macroscopic observables such as the neutron star mass-radius relation depend primarily on the equation of state (EOS). As a result, precise and reliable radius measurements would be valuable in constraining the EOS. However, most attempts at radius measurements are susceptible to systematic errors, meaning the inferred radius can be significantly biased even though the fit to the data appears to be statistically good. Previous studies suggested that radii inferred using X-ray data provided by NASA’s NICER mission may be more immune from such systematic errors. This is in part because, compared with previous measurements that obtained averaged spectra and fluxes, NICER observes millisecond pulsars by timing the photons so precisely, it is possible to obtain the spectrum as a function of rotational phase and see variations such as heated regions on the star rotate into and out of view hundreds of times per second. This extra information seems promising to break degeneracies and mitigate systematic errors, but a more in-depth study is necessary. We report the first steps of that study, in which we generate NICER-like synthetic data and determine the quality of fit and bias in radius obtained when we fit the data using a model different from the model used to generate the data.

Isiah Holt↗

A Census of X-ray Gas in NGC 1068: Results From 450 ks of CHANDRA High Energy Transmisson Grating Observations

We present models for the X-ray spectrum of the Seyfert 2 galaxy NGC 1068. These are fitted to data obtained using the High Energy Transmission Grating on Chandra. The data show line and radiative recombination continuum emission from a broad range of ions and elements. The models explore the importance of excitation processes for these lines including photoionization followed by recombination, radiative excitation by absorption of continuum radiation, and inner shell fluorescence. The models show that the relative importance of these processes depends on the conditions in the emitting gas and that no single emitting component can fit the entire spectrum. In particular, the relative importance of radiative excitation and photoionization/recombination differs according to the element and ion stage emitting the line. This in turn implies a diversity of values for the ionization parameter of the various components of gas responsible for the emission, ranging from log(E ) = 1 to 3. Using this, we obtain an estimate for the total amount of gas responsible for the observed emission. The mass flux through the region included in the HETG extraction region is approximately 0.3M/yr, assuming ordered flow at the speed characterizing the line widths. This can be compared with what is known about this object from other techniques.

galaxies: Seyfert↗

Analysis of resistance data on a good ceramic sample of YBa2Cu3O(7-x)

Resistance data between 91.4 and 200 K on a ceramic sample of YBa2Cu3O(7-x) are analyzed using a two-band model for normal-state conductance plus the Lawrence-Doniach theory for effects of superconducting fluctuations. The data is fitted to a rms accuracy of 0.16 percent of the resistance in the middle of the temperature range. The steepest slope of the resistance of the sample occurs at 91.13 K, and the analysis shows that the critical region extends to less than 0.3 K above Tc. Values of parameters found indicate that defect or impurity scattering within grains is unusually small in this sample. Some difficulties in understanding Hall data using either simple two-band models or more sophisticated one-band models are discussed briefly.

Eagles, David M.↗

Computational Bayesian Methods Applied to Complex Problems in Bio and Astro Statistics

In this dissertation we apply computational Bayesian methods to three distinct problems. In the first chapter, we address the issue of unrealistic covariance matrices used to estimate collision probabilities. We model covariance matrices with a Bayesian Normal-Inverse-Wishart model, which we fit with Gibbs sampling. In the second chapter, we are interested in determining the sample sizes necessary to achieve a particular interval width and establish non-inferiority in the analysis of prevalences using two fallible tests. To this end, we use a third order asymptotic approximation. In the third chapter, we wish to synthesize evidence across multiple domains in measurements taken longitudinally across time, featuring a substantial amount of structurally missing data, and fit the model with Hamiltonian Monte Carlo in a simulation to analyze how estimates of a parameter of interest change across sample sizes.

Elrod, Chris↗

WMAP Detection of Early Reionization

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systematic errors. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales the WMAP data show the temperature-polarization correlation expected from adiabatic density perturbations. The data for $\ell > 20$ agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. The existence of correlations on super-horizon scales provides significant support for inflationary cosmologies. We detect excess power on large angular scales compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift $11 < z < 30$ at 95\% confidence. A model-independent fit to reionization optical depth yields results consistent with the best-fit $\Lambda$CDM model, with best fit value $\tau = 0.17 \pm 0.04$ at 68\% confidence, including systematic and foreground uncertainties. This trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history. I will discuss the WMAP data and its implications for reionization in the early universe.

Kogut, A.↗

Soft X-ray observation of supernova remnant IC 443

The low-energy X-ray spectrum (0.15-3.0 keV) and X-ray centroid position of the supernova remnant IC 443 have been determined using data obtained by the HEAO-1 A-2 experiment. The spectrum exhibits a low-energy excess over and above that predicted by an extrapolation of the higher-energy spectrum determined by Parkes et al (1976). The excess can be attributed to line emission, predominately from high-ionization states of iron from a plasma at a temperature of 7.7 + or - 0.6 or -4.3 x 10 to the 6th K. Other models which fit the data are discussed, in particular, a lower-temperature component required to explain the Fe X optical emission.

Galas, C. M. F.↗

The nature of cometary dust as determined from infrared observations

The infrared measurements of comets, the compositional information available from interplanetary dust particles (IDPs), and the recent results of flybys to Comet Halley can help in restricting the nature and composition of cometary dust models (c.f., Proceedings of the 20th ESLAB Symposium on Exploration of Halley's Comet, 1986). Researchers tried to incorporate some of these results into a coherent model to account for the observed cometary infrared emission. The presence of 10 and 3.4 micron features in Comet Halley (c.f. Bregman et al. 1987; Wickramasinghe and Allen 1986) indicated the presence of at least two components in the grain material, namely silicates and some form of amorphous carbon. These two components could reside in separate grains or may be parts of composite particles. Both these cases have been considered (see Krishna Swamy el a. 1988a, 1988b). In the absence of refractive index data for cometary analogs, the authors used the optical constants of olivine-rich lunar material 12009.48 (Perry et al. 1972) for the infrared region and that of alpha:C-H film for amorphous carbon (angus et al. 1986). For the visible region, a value of m = 1.38-0.39i was used for the silicates, and values published by Arakawa et al. (1985) were used for the amorphous carbon. These materials should give a representative behavior of the expected results. The model results were compared to observational data. The strength of the 3.4 micron and 10 micron features relative to the adjacent continuum, as well as the slope of the continuum between 2500 and 1250 cm(exp -1) (4 to 8 microns), were used as criteria for comparison. Model calculations with alpha approx. equals -3.5, and also the size distribution function inferred for Comet Halley, with a mass fraction (X) of silicate to amorphous carbon grains of about 40 to 1 can fit the data. A good match is obtained for the infrared spectra of Comets Halley and West from a 40 to 1 mixture of silicate and amorphous carbon grains with a a(exp -3.5) size distribution function. The results are consistent with compositional constraints provided by interplanetary dust particles (IPDs) and Halley flyby data. The variation of grain temperature with heliocentric distance appears to account for the major changes observed in cometary spectra.

Swamy, K. S. Krishna↗

Direct Measurement of the Spectral Structure of Cosmic-Ray Electrons+Positrons in the Tev Region With Calet on the International Space Station

Detailed measurements of the spectral structure of cosmic-ray electrons and positrons from 10.6 GeV to 7.5 TeV are presented from over 7 years of observations with the CALorimetric Electron Telescope (CALET) on the International Space Station. The instrument, consisting of a charge detector, an imaging calorimeter, and a total absorption calorimeter with a total depth of 30 radiation lengths at normal incidence and a fine shower imaging capability, is optimized to measure the all-electron spectrum well into the TeV region. Because of the excellent energy resolution (a few percent above 10 GeV) and the outstanding e/p separation (10 5 ), CALET provides optimal performance for a detailed search of structures in the energy spectrum. The analysis uses data up to the end of 2022, and the statistics of observed electron candidates has increased more than 3 times since the last publication in 2018. By adopting an updated boosted decision tree analysis, a sufficient proton rejection power up to 7.5 TeV is achieved, with a residual proton contamination less than 10%. The observed energy spectrum becomes gradually harder in the lower energy region from around 30 GeV, consistently with AMS-02, but from 300 to 600 GeV it is considerably softer than the spectra measured by DAMPE and Fermi-LAT. At high energies, the spectrum presents a sharp break around 1 TeV, with a spectral index change from −3.15 to −3.91, and a broken power law fitting the data in the energy range from 30 GeV to 4.8 TeV better than a single power law with 6.9 sigma significance, which is compatible with the DAMPE results. The break is consistent with the expected effects of radiation loss during the propagation from distant sources (except the highest energy bin). We have fitted the spectrum with a model consistent with the positron flux measured by AMS-02 below 1 TeV and interpreted the electron+positron spectrum with possible contributions from pulsars and nearby sources. Above 4.8 TeV, a possible contribution from known nearby supernova remnants, including Vela, is addressed by an event-by-event analysis providing a higher proton-rejection power than a purely statistical analysis.

O. Adriani↗