Search NASA⌕ Search

SEARCH · Search NASA

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 First–Principles–Based Sub–Lattice Formalism for Predicting Off–Stoichiometry in Materials for Solar Thermochemical Applications: The Example of Ceria

Theoretical models that reliably can predict off-stoichiometry in materials via accurate descriptions of underlying thermodynamics are crucial for energy applications. For example, transition-metal and rare-earth oxides that can tolerate a large number of oxygen vacancies, such as CeO 2 and doped CeO 2 , can split water and carbon dioxide via a two-step, oxide-based solar thermochemical (STC) cycle. The search for new STC materials with a performance superior to that of state-of-the-art CeO 2 can benefit from predictions accurately describing the thermodynamics of oxygen vacancies. The sub-lattice formalism, a common tool used to fit experimental data and build temperature-composition phase diagrams, can be useful in this context. Here, sub-lattice models are derived solely from zero-temperature quantum mechanics calculations to estimate fairly accurate temperature- and oxygen-partial-pressure-dependent off-stoichiometries in CeO 2 and Zr-doped CeO 2 . Physical motivations for deriving some of the “excess” sub-lattice model parameters directly from quantum mechanical calculations, instead of fitting to minimize deviations from experimental and/or theoretical data, are identified. As a result, important limitations and approximations of the approach used are specified and extensions to multi-cation oxides are also suggested to help identify novel candidates for water and carbon dioxide splitting and related applications.

08 HYDROGEN↗

Flight inspection data and crack initiation times

Lockheed C-130 service-flight inspection data is reduced by means of fracture mechanics and statistical analysis into a form from which the probability of time to crack initiation and the distribution of initial flaw sizes for various locations on the aircraft can be determined. Crack sizes are normalized by extrapolating a growing crack backwards from its first recorded size to .03 inches (crack initiation size) yielding the time to crack initiation. In a similar fashion, the crack is grown backwards to a time equal to zero to yield an initial flaw size. In order to perform the discussed computations, crack growth constants were computed for each individual location and were analyzed statistically. Statistical distributions were fitted to times to crack initiation and to initial flaw sizes in order to describe their expected behavior. Weibull and Johnson distributions have been found to fit the data reasonably well but for a better fit convolution integrals of Poisson load distributions and Weibull strength distributions are needed.

Johnson, W. S.↗

Runaway electron plateau current profile reconstruction from synchrotron imaging and Ar-II line polarization angle measurements in DIII-D

Abstract Current profile reconstructions are obtained for high current ( I p ≃ 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Density of NaAlSi 2 O 6 Melt at High Pressure and Temperature Measured by In-Situ X-ray Microtomography

In this study, the volumetric compression of jadeite (NaAlSi 2 O 6 ) melt at high pressures was determined by three-dimensional volume imaging using the synchrotron-based X-ray microtomography technique in a rotation-anvil device. Combined with the sample mass, measured using a high-precision analytical balance prior to the high-pressure experiment, the density of jadeite melt was obtained at high pressures and high temperatures up to 4.8 GPa and 1955 K. The density data were fitted to a third-order Birch-Murnaghan equation of state, resulting in a best-fit isothermal bulk modulus K T0 of $10.8$ $^{+1.9}_{–5.3}$ GPa and its pressure derivative K' T0 of $3.4$ $^{+6.6}_{–0.4}$. Comparison with data for silicate melts of various compositions from the literature shows that alkali-rich, polymerized melts are generally more compressible than alkali-poor, depolymerized ones. The high compressibility of jadeite melt at high pressures implies that polymerized sodium aluminosilicate melts, if generated by low-degree partial melting of mantle peridotite at ~250–400 km depth in the deep upper mantle, are likely denser than surrounding mantle materials, and thus gravitationally stable.

36 MATERIALS SCIENCE↗

Distribution functions of intergalactic clouds

Width spectrum data derived from high-redshift QSO Ly-alpha absorption lines are used to study the statistical properties of intergalactic clouds. Several theoretical distribution functions for size and density are derived and tested against observational data for two types of line broadening: thermal and cosmological. Good agreement was found between the theoretical distribution functions and the observational data when the clouds were thermally broadened, and when distribution functions for both cloud size and width have a nonzero width. Attempts to fit the data with a single radius or density distribution function were unsuccessful.

Sherman, R. D.↗

Geodetic evidence of the seismic potential at Parkfield, California

The average rates of line-length change and the rates of the shallow fault slip observed near Parkfield, California, were used to estimate the slip rate at depth on the San Andreas fault (SAF), using the Matsu'ura et al. (1986) procedure to model the aseismic deformation as a perturbation to steady block motion across the SAF. Five different models were tested against the data, with one model, featuring a nearly monotonic transition from creeping behavior on the northwest to locked behavior on the southeast, found to fit the data acceptably well. It is shown that stored elastic energy southeast of Parkfield could produce an earthquake of magnitude up to 8.0, that would have displacement and rupture length similar to historic earthquakes of the same size.

Sung, Li-Yu↗

Effect on Gaseous Film Cooling of Coolant Injection Through Angled Slots and Normal Holes

A study was made to determine the effect of coolant injection angularity on gaseous film-cooling effectiveness. In the correlation of experimental data an effective injection angle was defined by a vector summation of the coolant and mainstream gas flows. The cosine of this angle was used as a parameter to empirically develop a corrective term to qualify a correlating equation presented in Technical Note D-130 that was limited to tangential injection of the coolant. Data were also obtained for coolant injection through rows of holes normal to the test plate. The slot correlating equation was adapted to fit these data by the definition of an effective slot height. An additional corrective term was then determined to correlate these data.

Papell, S. Stephen↗

NuSTAR Observations of the Magnetar 1E 2259+586

We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data show pulsations above 20 keV for the first time and we report evidence that one of the pulses in the double-peaked pulse profile shifts position with energy. The pulsed fraction of the magnetar is shown to increase strongly with energy. Our spectral analysis reveals that the soft X-ray spectrum is well characterized by an absorbed double blackbody or blackbody plus power-law model in agreement with previous reports. Our new hard X-ray data, however, suggest that an additional component, such as a power law, is needed to describe the NuSTAR and Swift spectrum. We also fit the data with the recently developed coronal outflow model by Beloborodov for hard X-ray emission from magnetars. The outflow from a ring on the magnetar surface is statistically preferred over outflow from a polar cap.

X-rays: individual↗

AN AFFINE-INVARIANT SAMPLER FOR EXOPLANET FITTING AND DISCOVERY IN RADIAL VELOCITY DATA

Markov chain Monte Carlo (MCMC) proves to be powerful for Bayesian inference and in particular for exoplanet radial velocity fitting because MCMC provides more statistical information and makes better use of data than common approaches like chi-square fitting. However, the nonlinear density functions encountered in these problems can make MCMC time-consuming. In this paper, we apply an ensemble sampler respecting affine invariance to orbital parameter extraction from radial velocity data. This new sampler has only one free parameter, and does not require much tuning for good performance, which is important for automatization. The autocorrelation time of this sampler is approximately the same for all parameters and far smaller than Metropolis-Hastings, which means it requires many fewer function calls to produce the same number of independent samples. The affine-invariant sampler speeds up MCMC by hundreds of times compared with Metropolis-Hastings in the same computing situation. This novel sampler would be ideal for projects involving large data sets such as statistical investigations of planet distribution. The biggest obstacle to ensemble samplers is the existence of multiple local optima; we present a clustering technique to deal with local optima by clustering based on the likelihood of the walkers in the ensemble. We demonstrate the effectiveness of the sampler on real radial velocity data.

79 ASTRONOMY AND ASTROPHYSICS↗

Hints of an axion-like particle mixing in the GeV gamma-ray blazar data?

Axion-Like Particles (ALPs), if exist in nature, are expected to mix with photons in the presence of an external magnetic field. The energy range of photons which undergo strong mixing with ALPs depends on the ALP mass, on its coupling with photons as well as on the external magnetic field and particle density configurations. Recent observations of blazars by the Fermi Gamma-Ray Space Telescope in the 0.1–300 GeV energy range show a break in their spectra in the 1–10 GeV range. We have modeled this spectral feature for the flat-spectrum radio quasar 3C454.3 during its November 2010 outburst, assuming that a significant fraction of the gamma rays convert to ALPs in the large scale jet of this blazar. Using theoretically motivated models for the magnetic field and particle density configurations in the kiloparsec scale jet, outside the broad-line region, we find an ALP mass m{sub a} ∼ (1−3)⋅10{sup −7} eV and coupling g{sub aγ} ∼ (1−3)⋅10{sup −10} GeV{sup −1} after performing an illustrative statistical analysis of spectral data in four different epochs of emission. The precise values of m{sub a} and g{sub aγ} depend weakly on the assumed particle density configuration and are consistent with the current experimental bounds on these quantities. We apply this method and ALP parameters found from fitting 3C454.3 data to another flat-spectrum radio quasar PKS1222+216 (4C+21.35) data up to 400 GeV, as a consistency check, and found good fit. We find that the ALP-photon mixing effect on the GeV spectra may not be washed out for any reasonable estimate of the magnetic field in the intergalactic media.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-E-T-595: Output Factor Calculation for a Uniform Scanning Proton Therapy System

Purpose: To develop an analytical model that predicts dose output factors for patient treatment fields in a uniform scanning proton therapy system. Methods: A model to predict output factors for patient specific treatment fields was produced based on the methods developed by Kooy et al. (2003). The Kooy model predicts the output factor based on the ratio of the entrance dose under calibration conditions to that for a given range and modulation corrected for the inverse square effect. Field specific output factors were plotted as a function of a single parameter r, where r = (Range-Modulation)/Modulation. The model targeted user range 1 sub-span 3 through user range 2 sup-span 2 of the IBA uniform scanning proton therapy system. The data set included points measured using the 10 cm and 18 cm snout sizes to eliminate stem effects on the monitor unit chambers. The data was fit using equation 15 from Kooy et al. (2003), and the resulting model was tested against measurements that were not included in the original data set. Results: For the range and sub span investigated, 120 data points were tested against the model prediction. The model predicted the output factor within 2% for 96% of the points tested and within 2.5% for 99% of the points tested. All points were within 3% of the predicted values. Conclusion: Monitor units for patient treatment fields with proton ranges that fall within the tested interval can be predicted using a model based on the methods developed at MGH. With further evaluation, it will be possible to model all user ranges and sub-spans of the IBA system. Further testing is also needed to predict output factors using a 25 cm snout which introduces variable scanning pattern sizes and stem effects on the monitor unit chambers.

61 RADIATION PROTECTION AND DOSIMETRY↗

Thermal Contact Conductance Modeling

Elastic, plastic and elastic-plastic thermal conductance models were run for Al-6061 and SS-304, respectively. Thermal conductance models were compared to experimental data of samples being pressed together as a function of contact pressure. For SS-304, Elastic and Plastic models did not fit the experimental data well. The Elastic-Plastic model fit the data well and it was shown that surface roughness did not play a major role. For Al-6061, the Elastic model worked best and the models strongly suggested that surface roughness Al plays an important role in how well a thermal conductance model can predict conductance through two pieces of Al in contact with one another.

36 MATERIALS SCIENCE↗

Present-day plate motions

An instantaneous plate-motion model, Relative Motion 2 (RM2), is obtained by inverting a data set comprising 110 spreading rates, 78 transform fault azimuths, and 142 earthquake slip vectors. RM2 is compared with angular velocity vectors which best fit the data along individual plate boundaries and, while the model performs close to optimally in most regions, attention is directed to those regions which are not suitably described by the model. Reasons for the discrepancies between RM2 and observations for the India-Antarctica plate boundary, the Pacific-India plate boundary, and the east-west trending transform fault azimuths observed in the French-American Mid-Ocean Undersea Study area are discussed.

Minster, J. B.↗

Current QSO statistics - Implications for the intergalactic medium

The results of numerous QSO surveys have been compiled and analyzed to form a single spatially averaged QSO ionizing function that is independent of evolution mode (number density or luminosity). An intergalactic medium (IGM) that satisfies the observational constraints may now have a density with respect to closure of approximately 0.1, which is substantially less than hitherto modeled. The aggregate of QSO data also indicate that if evolution is due to number density, then, when split into luminosity classes, an exponential in look-back time fits the data better than a power law, and the evolution rates increase roughly with absolute luminosity if analyzed by the exponential.

Sherman, R. D.↗

More data on (possible) gamma ray (point) sources

The 2CG catalog of gamma ray sources was compiled before detailed knowledge was available on the fine-scale structure of the diffuse Galactic gamma-ray emission. Two independent analyses to discriminate sources which are either compact objects of due to very local and strong enhancements in the Galactic cosmic-ray distribution from those which are artifacts due to the clumpy gas distribution are about to be completed: a maximum likelihood analysis and a cross correlation analysis. Arguments are given why differences, and therefore confusion, and in resulting source lists can be expected. Detailed analysis of all COS-B gamma-ray data on Geminga (2CG195+04), reveals the existence of a drastic spectral break below 200 MeV. A power-law spectrum with index -1.88 fits the data above about 100 MeV to 3.2 GeV, however, there are also indications for a spectral break above these energies. For energies above about 100 MeV no evidence for a long-term time variability was found. The error region of Geminga was searched for a radio counterpart at wavelengths of 90, 49, 21, 6, and 2 cm using the Westerbork Synthesis Radio Telescope and the Very Large Array. So far, 16 sources were detected in this error region. In the direction of 1E0630+178, the Einstein x ray source proposed to be a Vela-like pulsar and the counterpart of Geminga, no radio source was found at 21, 49, and 90 cm with 3 sigma upper limits on the flux densities ranging from 0.5 mJy at 21 cm to 4.5 mJy at 90 cm. Detailed structures in local molecular cloud complexes are so far only resolved in gamma rays for the closet and most massive complexes, namely those in the Orion-Monoceros and the Ophiuchus regions. For both region, there is circumstantial evidence for gamma ray emission from molecular gas that was photodissociated after the passage of a SN shell.

Hermsen, W.↗

Quenching factor measurements of neon nuclei in neon gas

The NEWS-G collaboration uses spherical proportional counters (SPCs) to search for weakly interacting massive particles (WIMPs). Here, in this paper, we report the first measurements of the nuclear quenching factor in neon gas at 2 bar using an SPC deployed in a neutron beam at the TUNL facility. The energy-dependence of the nuclear quenching factor is modeled using a simple power law: $αE^{β}_{nr}$; we determine its parameters by simultaneously fitting the data collected with the detector over a range of energies. We measured the following parameters in Ne : CH 4 at 2 bar: α= 0.2801±0.0050 (fit) ±0.0045 (sys) and β=0.0867±0.020 (fit) ±0.006 (sys). Our measurements do not agree with expected values from SRIM or Lindhard theory. We demonstrated the feasibility of performing quenching factor measurements at sub-keV energies in gases using SPCs and a neutron beam.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Demonstration of Intelligent Control and Fan Improvements in Computer Room Air Handlers

This report documents a demonstration of the energy-efficiency improvement provided by a new control system for computer room air handling devices. It also analyzes measured and reported air handling device fan power associated with changing the fan type. A 135,000 square foot commercial data center was used for the demonstration. All air handling units were upgraded with improved efficiency fans, and a control system that automatically adjusts the fan speed for the air handling units was added. Power measurements were collected for a baseline and for a period with the fan speed control system active. Changing the fan type resulted in a savings of 47 percent of energy used by the air handling equipment and associated chiller plant energy needed to cool the air handlers themselves. The addition of the fan speed control resulted in an additional 37 percent savings in the same two categories. The combined savings for the two improvements for the same categories was 66 percent compared to the data center fitted with the original fans without a control system. The energy use reduction provided by the complete air handling device improvement program for the whole data center site is estimated to be 2.9 million kilowatt hours per year—an overall data center site savings of 8.0 percent. The reduced electrical energy use at the site provides a 1.9 million pound yearly reduction of carbon dioxide emissions. This demonstration showed that fan upgrades and a control system addition provide cost-effective improvements for data centers, with a payback reported to be under two years without utility incentives. In addition to the control system providing energy savings, the data collection and visual analysis capabilities provided immediate and long-term benefits. It is recommended that data center operators consider investing in fan upgrades and/or adding fan speed control for computer room air handlers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Method for the Quantification of Model Form Error Associated with Physical Systems

In the process of model validation, models are often declared valid when the differences between model predictions and experimental data sets are satisfactorily small. However, little consideration is given to the effectiveness of a model using parameters that deviate slightly from those that were fitted to data, such as a higher load level. Furthermore, few means exist to compare and choose between two or more models that reproduce data equally well. These issues can be addressed by analyzing model form error, which is the error associated with the differences between the physical phenomena captured by models and that of the real system. This report presents a new quantitative method for model form error analysis and applies it to data taken from experiments on tape joint bending vibrations. Two models for the tape joint system are compared, and suggestions for future improvements to the method are given. As the available data set is too small to draw any statistical conclusions, the focus of this paper is the development of a methodology that can be applied to general problems.

97 MATHEMATICS AND COMPUTING↗