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 415 records · Page 23

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

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

Tsai, John C.↗

Culling of redundant data

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

Gilbert, F.↗

Gravity - Mare Humorum.

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

Sjogren, W. L.↗

Measurement of Neutron Multiplicity in Charged Current Neutrino Interactions on Oxygen

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located 100~meters downstream in the Booster Neutrino Beam (BNB) at the Fermi National Accelerator Laboratory (Fermilab). Its primary goals are to (1) measure the neutron yield from $\nu_\mu$ interactions as a function of momentum transfer $Q^2$ so that neutrino-nucleus interaction models can be better constrained, and (2) demonstrate the power of novel, fast-timing detectors with the first deployment of Large Area Picosecond PhotoDetectors (LAPPDs). Current knowledge of neutrino-nucleus interactions fall short in modeling the topologies of such interactions, leading to inaccurate reconstruction of event kinematics such as particle energy, direction, and vertex. As a consequence, the accuracy of cross section measurements is impacted, which is necessary for precise physics measurements. Neutrons are an indication of inelasticity and affect the determination of the energy of the parent neutrino. Quantifying the neutron yield is a step towards reducing the associated uncertainties, and thus improving our understanding of these complex interactions and benefiting the next generation of long-baseline neutrino experiments. ANNIE will make use of LAPPDs to measure neutron multiplicity of CC-0$\pi$ $\nu_\mu$ interactions on oxygen, making it the first experiment to deploy an array of these photodetectors. Because the LAPPD is a novel photodetection technology, much customization is required to integrate it into existing electronics. The first half of this thesis covers the significant progress made towards the first deployment of the LAPPD system. From its test stand at Fermilab, the LAPPD system was systematically tested and put together until deployment readiness was achieved. I present my contributions to the design, fabrication, and testing of the waterproof housing and cables, and the commissioning of the LVHV board that powers the LAPPD and its readout electronics. These efforts brought the LAPPD system significantly closer to deployment, and eventually first data. The second half of this thesis presents the vertex and energy reconstruction algorithms developed to analysis the beam data with PMT-only information. While much progress has been made towards the deployment of LAPPDs, with several in the detector tank, efforts to integrate the LAPPD datastream are in progress. Thus, I developed a ring edge detection technique using PMT data to fit the muon vertex and determine its energy. The analysis in this thesis finds average neutron yields of $\Bar{n}_{data} = 0.452 \pm 0.039 (\text{stat}) \pm 0.27 (\text{sys})$ for a selection of muon neutrino candidates in the fiducial volume of ANNIE and corrected with an averaged neutron detection efficiency. An equivalent analysis for simulated beam data results in an average neutron yield of $\Bar{n}_{MC} = 0.582 \pm 0.018 (\text{stat}) \pm 0.25 (\text{sys})$. Future work includes application of efficiency corrections at a positional level, quantification of all systematic uncertainties, and neutron multiplicity measurements with other event topologies.

43 PARTICLE ACCELERATORS↗

CrysFieldExplorer : rapid optimization of the crystal field Hamiltonian

A new approach to the fast optimization of crystal electric field (CEF) parameters to fit experimental data is presented. This approach is implemented in a lightweight Python-based program, CrysFieldExplorer. The main novelty of the method is the development of a unique loss function, referred to as the spectrum characteristic loss (L Spectrum ), which is based on the characteristic polynomial of the Hamiltonian matrix. Particle swarm optimization and a covariance matrix adaptation evolution strategy are used to find the minimum of the total loss function. It is demonstrated that CrysFieldExplorer can perform direct fitting of CEF parameters to any experimental data such as a neutron spectrum, susceptibility or magnetization measurements etc. CrysFieldExplorer can handle a large number of non-zero CEF parameters and reveal multiple local and global minimum solutions. Crystal field theory, the loss function, and the implementation and limitations of the program are discussed within the context of two examples.

42 ENGINEERING↗

Predicting phase behavior in high entropy and chemically complex alloys

The interest in high entropy alloys and other metallic compounds with four or more elements at near-equiatomic ratios has drawn attention to the ability to rapidly predict phase behavior of these complex materials, particularly where existing thermodynamic data are lacking. In this paper, we discuss aspects of this from the point of view of predicting without utilizing (or fitting) experimental data. Of particular interest are heuristic approaches that provide prediction of single-phase compositions, more rigorous approaches that tackle the thermodynamics from a more fundamental point of view, and simulation approaches that provide further insight into the behaviors. This paper covers cases of all three of these, in order to examine the strengths and weaknesses of each approach, and to indicate directions where these may be utilized and improved upon. Of particular interest is moving beyond “which composition may form a solid solution,” to recognizing the importance of underlying thermodynamic realities that affect the temperature- and composition-dependent transformations of these materials.

36 MATERIALS SCIENCE↗

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

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

Wagener, R.↗

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

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

Caldwell, E. C.↗

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

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

Nisbet, J. S.↗

Assessing observational constraints on dark energy

Observational constraints on time-varying dark energy (e.g., quintessence) are commonly presented on a w 0 –w a plot that assumes the equation of state of dark energy strictly satisfies w(z) = w 0 + w a z/(1 + z) as a function of the redshift z. Recent observations favor a sector of the w 0 –w a plane in which w 0 > –1 and w 0 + w a < –1, suggesting that the equation of state underwent a transition from violating the null energy condition (NEC) at large z to obeying it at small z. In this paper, we demonstrate that this impression is misleading by showing that simple quintessence models satisfying the NEC for all z predict an observational preference for the same sector. We also find that quintessence models that best fit observational data can predict a value for the dark energy equation of state at present that is significantly different from the best-fit value of w 0 obtained assuming the parameterization above. In addition, the analysis reveals an approximate degeneracy of the w 0 –w a parameterization that explains the eccentricity and orientation of the likelihood contours presented in recent observational studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Computer Modeling of Non-Isothermal Crystallization

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

Kelton, K. F.↗

The Future is Hera! Analyzing Astronomical Over the Internet

Hera is the data processing facility provided by the High Energy Astrophysics Science Archive Research Center (HEASARC) at the NASA Goddard Space Flight Center for analyzing astronomical data. Hera provides all the pre-installed software packages, local disk space, and computing resources need to do general processing of FITS format data files residing on the users local computer, and to do research using the publicly available data from the High ENergy Astrophysics Division. Qualified students, educators and researchers may freely use the Hera services over the internet of research and educational purposes.

Valencic, L. A.↗

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↗