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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
DOE OSTI2022

Thermoelastic properties of zircon: Implications for geothermobarometry

Abstract A thermal-pressure equation of state has been determined for zircon (ZrSiO4) that characterizes its thermoelastic behavior at metamorphic conditions. New pressure-volume (P-V) data from a “Mud Tank” zircon have been collected from 1 bar to 8.47(1) GPa using X-ray diffraction, and elastic moduli were measured from room temperature up to 1172 K by resonance ultrasound spectroscopy. These data were fitted simultaneously with temperature-volume (T-V) data from the literature in EosFit7c using a new scaling technique. The parameters of a third-order Birch-Murnaghan EoS with a Mie-Grüneisen-Debye model for thermal pressure have compressional EoS parameters K0T = 224.5(1.2) GPa, K0T' = 4.90(31) with a fixed initial molar volume V0 = 39.26 cm3/mol and thermal parameters γ0 = 0.868(15), q = 2.37(80), and ΘD = 848(38) K. EoS parameters that describe the variation of unit-cell parameters with pressure and temperature were determined using an isothermal-type EoS. This new EoS confirms that zircons are stiffer than garnets and exhibit a much lower thermal expansion. This results in steep isomekes between zircon and garnets, which makes zircon trapped as inclusions in garnets at metamorphic conditions a good piezothermometer.

Geochemistry & Geophysics↗
DOE OSTI2021

Critical and Subcritical Californium Source Driven Noise Analysis Experiments with Fresh PWR Fuels Pins

These experiments were performed in 1983 between August 1 and September 16 and used 35 operational days of critical facility time at the Babcock & Wilcox critical facility in Lynchburg, Va. These measurements were to assess the capability of the Cf source driven noise analysis method to measure the subcriticality of light water reactor fresh fuel configurations. This report documents the experiments of 1983 that were not reported at that time. These measurement with 2.459 wt. % 235 U enriched uranium oxide fresh fuel pins started at the delayed critical configuration of the fully assembled pressurized water reactor (PWR) core configuration with 1511 parts per million (ppm) natural boron (19.8 wt. % 10 B) in the water. This was a cylindrical configuration of 4961 fuel pins. After the initial measurements at delayed criticality, subcritical measurements by the Cf source driven noise analysis method were performed in which the number of detectors and Cf source locations were varied. Additional boron was added to the water moderator and reflector to reduce the neutron multiplication factor, with resulting boron concentrations of 1561, 1613, 1765, 1880, 2104, 2384, 2975, 3606, and 4303 ppm. After the completion of these measurements the increased, the boron concentration was returned to that at delayed critical (1511 ppm) and the reactivity was decreased by reducing the number of fuel pins in steps to 3713, 2553, 1281, 749, and 333 fuel pins and then to that (289) of a s ingle square 17 x 17 pin array, corresponding to a single PW fuel element. For all measurement that started at delayed criticality and then reduced the neutron multiplication factor to various subcritical states, inverse kinetics rod drop measurements at the higher reactivities were also performed but not reported here. At all subcritical states the break frequency noise analysis data was fitted to obtain the prompt neutron decay constant. Some data presented in this report are from notes not in the logbook. The purpose of this report is to document the experimental information for the measurements performed so that a later date someone could perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Program (INCSBEP) or Nuclear Energy Agency (NEA) benchmark. The data from these measurements is available from the Laboratory Records Department of Oak Ridge National Laboratory and the logbook is available from INCSBEP at Idaho National Laboratory.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Exploring the Whole Set of Accurate Sparse Interpretable Models

In data science applications, there are often many models that fit the data well. This phenomenon was called the Rashomon Effect by Leo Breiman. The set of good models is called the Rashomon Set, and the goal of this project is to locate, store, and study the Rashomon sets for classes of interpretable models, including decision trees and generalized additive models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Perspective on Traditional and Data Driven Electrochemical Modeling and Analysis

To understand the behavior of electrochemical systems, we need to reduce the dimensionality of the measured current-voltage-time (I-V-t) data by fitting models, thus enabling us to analyze and compare the governing physics. Traditionally, the process for this is an 'expert first' approach: defining the model and its explicit assumptions based on inductive reasoning or empirical observation, fitting small portions of the I-V-t data where assumptions are most valid or carefully designing experiments to enforce key assumptions, and then interpreting the model parameters. However, modern data-driven methods enable a new paradigm: a 'data first' approach, where the latent behaviors governing the system's measured response are identified directly using machine-learning models that optimize both model structure and parameters from the I-V-t data, guaranteeing that the learned model explains as much of the observed system response as possible. After model identification, the model can then be interrogated by an expert to connect observed behaviors with underlying physics. This talk will review several different types of electrochemical analysis (electrochemical impedance, differential voltage-capacity, electrochemical kinetics) and compare the traditional and data-driven methods for analyzing the data.

42 ENGINEERING↗
NASA NTRS1982

Universal approach to analysis of cavitation and liquid-impingement erosion data

Cavitation erosion experimental data was analyzed by using normalization and curve-fitting techniques. Data were taken from experiments on several materials tested in both a rotating disk device and a magnetostriction apparatus. Cumulative average volume loss rate and time data were normalized relative to the peak erosion rate and the time to peak erosion rate, respectively. From this process a universal approach was derived that can include data on specific materials from different test devices for liquid impingement and cavitation erosion studies.

Rao, P. V.↗
NASA NTRS1987

Orbital atmospheric physics and dynamics

There are two ways of modeling the upper atmosphere. One is the empirical model that makes use of experimental data on means and excursions from the mean and fits the data in a self-consistent manner. The other approach is to deal directly with the physical processes. This is difficult since what is happening is extremely complex. Data measured using an interferometer to give Doppler shifts of airglow lines showed 300 to 800 m/sec winds with a complex structure in the upper region of the thermosphere at high latitudes. Ionospheric electric fields, strongly influenced by interaction with the solar wind, drive the ionized component and large neutral winds result due to momentum transfer between the charged particles and the neutrals. Frictional heating results from movement of ions through the neutrals, which also influences the compositional structure. These are examples of the complex interactions involved. The NCAR General Circulation Model (tropospheric) was adapted for use at thermospheric altitudes: the Thermospheric General Circulation Model (TGCM). The model makes use partly of primitive equations and partly of empirical data for some quantities such as electron density, magnetic field, and ion drift.

Roble, Raymond↗
NASA NTRS1994

Predicting Fatigue Lives Of Metal-Matrix/Fiber Composites

Method of prediction of fatigue lives of intermetallic-matrix/fiber composite parts at high temperatures styled after method of universal slopes. It suffices to perform relatively small numbers of fatigue tests. Data from fatigue tests correlated with tensile-test data by fitting universal-slopes equation to both sets of data. Thereafter, universal-slopes equation used to predict fatigue lives from tensile properties.

Bartolotta, Paul A.↗
DOE OSTI2023

Encoding nonlinear and unsteady aerodynamics of limit cycle oscillations using nonlinear sparse Bayesian learning

This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Streak camera based SLR receiver for two color atmospheric measurements

To realize accurate two-color differential measurements, an image digitizing system with variable spatial resolution was designed, built, and integrated to a photon-counting picosecond streak camera, yielding a temporal scan resolution better than 300 femtosecond/pixel. The streak camera is configured to operate with 3 spatial channels; two of these support green (532 nm) and uv (355 nm) while the third accommodates reference pulses (764 nm) for real-time calibration. Critical parameters affecting differential timing accuracy such as pulse width and shape, number of received photons, streak camera/imaging system nonlinearities, dynamic range, and noise characteristics were investigated to optimize the system for accurate differential delay measurements. The streak camera output image consists of three image fields, each field is 1024 pixels along the time axis and 16 pixels across the spatial axis. Each of the image fields may be independently positioned across the spatial axis. Two of the image fields are used for the two wavelengths used in the experiment; the third window measures the temporal separation of a pair of diode laser pulses which verify the streak camera sweep speed for each data frame. The sum of the 16 pixel intensities across each of the 1024 temporal positions for the three data windows is used to extract the three waveforms. The waveform data is processed using an iterative three-point running average filter (10 to 30 iterations are used) to remove high-frequency structure. The pulse pair separations are determined using the half-max and centroid type analysis. Rigorous experimental verification has demonstrated that this simplified process provides the best measurement accuracy. To calibrate the receiver system sweep, two laser pulses with precisely known temporal separation are scanned along the full length of the sweep axis. The experimental measurements are then modeled using polynomial regression to obtain a best fit to the data. Data aggregation using normal point approach has provided accurate data fitting techniques and is found to be much more convenient than using the full rate single shot data. The systematic errors from this model have been found to be less than 3 ps for normal points.

Varghese, Thomas K.↗
DOE OSTI2024

Light Detecting and Ranging (LIDAR) for in-situ heliostat optical error assessment

A system and method for optical assessment of a heliostat includes obtaining a point cloud data representing an image of the heliostat; isolating the data; filtering and fitting the filtered heliostat data to a bounding box; translating the heliostat data to a plane to aid in segmentation; segmenting a plurality of facets of the heliostat fitting each of the segmented facets to a respective plane; generating normal vectors characterizing each of the plurality of facets; and calculating a canting angle associated with each respective facet of the plurality of facets. A heliostat with mirrored facets and a scanner are provided. The scanner captures point cloud data representing the heliostat, which is segmented for each facet. Normal vectors characterize the facets and a canting angle is calculated for the respective facet.

Small, Daniel E.↗
DOE OSTI2022

Surface analysis insight note: Differentiation methods applicable to noisy data for determination of sp2‐ versus sp3‐hybridization of carbon allotropes and AES signal strengths

The derivatives of the spectra are commonly used for quantification in Auger Electron Spectroscopy (AES) spectra, while the derivative of the KLL C Auger line has proven to be valuable in obtaining a measure of the relative proportions of sp 2 ‐ and sp 3 ‐hybridization using the D‐parameter in both AES and X‐ray Photoelectron Spectroscopy (XPS). Differentiation of X‐ray Photoelectron Spectroscopy (XPS) and Auger Electron Spectroscopy (AES) spectra by numerical means is presented and illustrated for polymeric, such as PEEK and Nylon, as well as for graphitic materials including highly ordered pyrolytic graphite and graphene oxide. The most commonly available Savitzky–Golay method is explained mathematically and developed through the case of constructing a 5‐point quadratic polynomial convolution kernel suitable for differentiating spectra of adequate signal to noise. The concept of differentiation of spectra where signal to noise is less than adequate is also developed. Two alternative strategies to Savitzky–Golay differentiation are presented, which fit curves to data that allow derivatives to be obtained where Savitzky–Golay would otherwise fail. These alternative methods involve constructing a parametric curve that fits data over the entire energy interval of interest. Derivatives of spectra are then obtained by differentiating these parametric curves directly. A comparison of results for different materials for which specific sp 2 ‐ vs sp 3 ‐hybridized carbon proportions are of interest is used to emphasize the importance of characterizing methods used to differentiate spectra and understanding the characteristics of instrumentation used to measure spectra. The case for using Principal Component Analysis noise reduction with C KLL spectra is made for spectra collected from a heterogeneous graphene oxide sample.

Fairley, Neal↗
DOE OSTI2023

Measurement of neutrino oscillations using neutrino and antineutrino beams in the NOvA experiment

NOvA is a long-baseline accelerator neutrino oscillation experiment using the NuMIneutrino beam from Fermilab. Its main physics goals are to probe the 3-flavour oscillationparameters: neutrino mass hierarchy, CP-violating phase dcp and octant of .23 mixingangle by observing electron neutrino appearance and muon neutrino disappearance. Twofunctionally identical detectors are placed off-axis from the centre of the NuMI beam.The near detector at Fermilab is 100 m underground, and the far detector is locatedon the surface at Ash River, 810 km away from the beam source. The initial neutrinobeam spectra are measured using the near detector data and the oscillation parametersare extracted by fitting the observed data to the predicted neutrino spectrum in the fardetector.This thesis is centered around how to improve the sensitivity of |.m232| and .23 measurementsin the muon neutrino disappearance analysis. NOvA will take data for about12 years. The operation of the NOvA experiment for each year costs tens of millions ofdollars, thus it is valuable to maximise the sensitivity of the analysis. Three samples ofmuon neutrino events are studied in this thesis to improve the analysis sensitivities. First,higher energy muon neutrinos are investigated by extending the energy range in NOvA’scurrent standard analysis. Second, for the sample of events used in NOvA’s existing analysis,a new energy estimator which has been developed to improve the neutrino energy resolution is considered. Furthermore, in addition to binning the events as function ofenergy and hadronic energy fraction, three particle identifiers are introduced to separateneutrino events by signal purity to reduce the effects from backgrounds. Third, an additionallower purity sample of muon neutrino charged current (CC) events that look similarto neutral current events and have not been included in NOvA’s existing analyses havebeen studied.This thesis reanalyses NOvA’s data used in the 2020 analysis, corresponding to anexposure of 13.60×1020 protons on target (POT) in the neutrino beam mode recordedfrom February 6, 2014 to March 20, 2020, and 12.50×1020 protons on target in theantineutrino beam mode recorded between June 29, 2016 to February 26, 2019. Thisthesis has implemented a fit to Asimov fake data, generated where sin2 .23 = 0.59 and.m232 = 2.5 × 10-3 eV2. These sensitivity studies show that the uncertainty range of|.m232| at 1 s in the new analysis is reduced by 5.5% and the significance of maximaldisappearance rejection improves by 7.7%, compared to the standard analysis. This isequivalent to adding 11-16% more data. The best fit values of the oscillation parametersfrom fitting to the far detector (FD) data with the new analysis are found to besin2 .23 = 0.568+0.025-0.043 (sin2 .23 = 0.454+0.046-0.026) and .m232 = 2.399+0.055-0.070 × 10-3 eV2 (.m232= -2.427+0.055-0.067 × 10-3 eV2) for the normal (inverted) hierarchy. The difference in thebest fit for sin2 .23 (.m232) between the new analysis and NOvA’s 2020 analysis is around2% (1.4%). The uncertainty range at 1 s for .m232 decrease by 8% (4%) for the normalhierarchy (inverted hierarchy) compared to the standard analysis. The uncertainty rangefor sin2 .23 is close to the standard analysis. This thesis also implements the fit from combiningelectron neutrino appearance and muon neutrino disappearance. The combinedanalysis shows that the best fit values are very close to the standard analysis. However,the uncertainty range of .m232 at 1 s is reduced by 3.7% using the new analysis. The maximaldisappearance significance is not improved in the new analysis, but the new analysisslightly improves the rejection of the disfavoured octant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1979

Scale covariance and G-varying cosmology. III - The /m, z/, /theta sub m, z/, /theta sub i, z/, and /N/m/, m/ tests

The Einstein gravitational equations in atomic units are fully solved for a matter-dominated universe in the context of a recently proposed scale-covariant cosmology. The magnitude-redshift relation for elliptical galaxies is studied, the evolutionary parameter used in such a study is derived, and the relation between isophotal angular diameters and redshifts is investigated, along with the relation between metric angular diameters and redshifts, the N(m)-magnitude relation for QSOs, and the magnitude-redshift relation for QSOs. Results are presented for four gauges (i.e., relations between G and the scale function beta (t)), and no contradictions are found between the proposed theory and the observational data. It is shown that only an open universe can fit the data if certain gauges suggested by a recent analysis of the time variation of the moon's period are selected and that observations made with atomic instruments do not necessarily yield geometrical parameters unless specific assumptions are made regarding the relation between atomic and gravitational dynamics.

Canuto, V. M.↗
NASA NTRS1994

Coronal temperatures of unusually active K-dwarf binary systems

We report the results of a ROSAT pointed study of 4 BY Dra systems. Good quality pulse-height spectra are available from all four systems. Except for a required interstellar absorption component in HD 319139, the four systems have remarkably similar x-ray spectra; the two systems BD +22deg.669 and BD +23deg.635 look virtually identical in x rays. Analysis of the 4 x-ray spectra reveals that, in all cases, a single-temperature hot plasma (RS or Mewe) spectra is inadequate to fit the data, and two temperatures are required. We present examples of fitted pulse-height spectra and chi squared contours in kT(sub 1)-kT(sub 2) space.

Stern, Robert A.↗
NASA NTRS1995

The point spread function of the soft X-ray telescope aboard Yohkoh

The point spread function of the SXT telescope aboard Yohkoh has been measured in flight configuration in three different X-ray lines at White Sands Missile Range. We have fitted these data with an elliptical generalization of the Moffat function. Our fitting method consists of chi squared minimization in Fourier space, especially designed for matching of sharply peaked functions. We find excellent fits with a reduced chi squared of order unity or less for single exposure point spread functions over most of the CCD. Near the edges of the CCD the fits are less accurate due to vignetting. From fitting results with summation of multiple exposures we find a systematic error in the fitting function of the order of 3% near the peak of the point spread function, which is close to the photon noise for typical SXT images in orbit. We find that the full width to half maximum and fitting parameters vary significantly with CCD location. However, we also find that point spread functions measured at the same location are consistent to one another within the limit determined by photon noise. A 'best' analytical fit to the PSF as function of position on the CCD is derived for use in SXT image enhancemnent routines. As an aside result we have found that SXT can determine the location of point sources to about a quarter of a 2.54 arc sec pixel.

Martens, Petrus C.↗
NASA NTRS2013

A Stochastic Fractional Dynamics Model of Space-time Variability of Rain

Rainfall varies in space and time in a highly irregular manner and is described naturally in terms of a stochastic process. A characteristic feature of rainfall statistics is that they depend strongly on the space-time scales over which rain data are averaged. A spectral model of precipitation has been developed based on a stochastic differential equation of fractional order for the point rain rate, that allows a concise description of the second moment statistics of rain at any prescribed space-time averaging scale. The model is thus capable of providing a unified description of the statistics of both radar and rain gauge data. The underlying dynamical equation can be expressed in terms of space-time derivatives of fractional orders that are adjusted together with other model parameters to fit the data. The form of the resulting spectrum gives the model adequate flexibility to capture the subtle interplay between the spatial and temporal scales of variability of rain but strongly constrains the predicted statistical behavior as a function of the averaging length and times scales. We test the model with radar and gauge data collected contemporaneously at the NASA TRMM ground validation sites located near Melbourne, Florida and in Kwajalein Atoll, Marshall Islands in the tropical Pacific. We estimate the parameters by tuning them to the second moment statistics of radar data. The model predictions are then found to fit the second moment statistics of the gauge data reasonably well without any further adjustment.

TRMM↗
NASA NTRS1996

Development of a Computer Architecture to Support the Optical Plume Anomaly Detection (OPAD) System

The NASA OPAD spectrometer system relies heavily on extensive software which repetitively extracts spectral information from the engine plume and reports the amounts of metals which are present in the plume. The development of this software is at a sufficiently advanced stage where it can be used in actual engine tests to provide valuable data on engine operation and health. This activity will continue and, in addition, the OPAD system is planned to be used in flight aboard space vehicles. The two implementations, test-stand and in-flight, may have some differing requirements. For example, the data stored during a test-stand experiment are much more extensive than in the in-flight case. In both cases though, the majority of the requirements are similar. New data from the spectrograph is generated at a rate of once every 0.5 sec or faster. All processing must be completed within this period of time to maintain real-time performance. Every 0.5 sec, the OPAD system must report the amounts of specific metals within the engine plume, given the spectral data. At present, the software in the OPAD system performs this function by solving the inverse problem. It uses powerful physics-based computational models (the SPECTRA code), which receive amounts of metals as inputs to produce the spectral data that would have been observed, had the same metal amounts been present in the engine plume. During the experiment, for every spectrum that is observed, an initial approximation is performed using neural networks to establish an initial metal composition which approximates as accurately as possible the real one. Then, using optimization techniques, the SPECTRA code is repetitively used to produce a fit to the data, by adjusting the metal input amounts until the produced spectrum matches the observed one to within a given level of tolerance. This iterative solution to the original problem of determining the metal composition in the plume requires a relatively long period of time to execute the software in a modern single-processor workstation, and therefore real-time operation is currently not possible. A different number of iterations may be required to perform spectral data fitting per spectral sample. Yet, the OPAD system must be designed to maintain real-time performance in all cases. Although faster single-processor workstations are available for execution of the fitting and SPECTRA software, this option is unattractive due to the excessive cost associated with very fast workstations and also due to the fact that such hardware is not easily expandable to accommodate future versions of the software which may require more processing power. Initial research has already demonstrated that the OPAD software can take advantage of a parallel computer architecture to achieve the necessary speedup. Current work has improved the software by converting it into a form which is easily parallelizable. Timing experiments have been performed to establish the computational complexity and execution speed of major components of the software. This work provides the foundation of future work which will create a fully parallel version of the software executing in a shared-memory multiprocessor system.

Katsinis, Constantine↗