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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The CXSFIT spectral fitting code: Past, present and future

Magnetically confined plasma experiments generate a wealth of spectroscopic data. The first step toward extracting physical parameters is to fit a spectral model to the often complex spectra. The CXSFIT (Charge eXchange Spectroscopy FITting) spectral fitting code was originally developed for fitting charge exchange spectra on JET from the late 1980s onward and has been further developed over decades to keep up with the needs of the users. The primary use is to efficiently fit a large number of spectra with many constrained Gaussian spectral lines of which the physical parameters can be coupled in a user-friendly manner. More recent additions to the code include time-dependent couplings between parameters, flexible background subtraction, and a non-linear coupling scheme between fit parameters. The latter was a pre-requisite for implementing Zeeman and motional Stark effect multiplets in the library of spectral features. The ability to save and replay “fit recipes,” even when multiple iterations are required, has ensured the traceability of the results and is one of the keys to the longevity and success of the code. The code is also in use on other tokamaks (AUG, ST-40) and to fit data from other spectroscopic diagnostics on JET. In this paper, we document the current capabilities and philosophy behind the structure of the code, including some of the algorithms used to calculate spectral features numerically efficiently. We also provide an outline of how CXSFIT could be transferred into a framework that would be able to meet the spectral fitting requirements of future devices, such as ITER.

Delabie, Ephrem G.↗
DOE OSTI2023

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The Dark Energy Survey Supernova Program: an updated measurement of the Hubble constant using the inverse distance ladder

We measure the current expansion rate of the Universe, Hubble’s constant $H_0$, by calibrating the absolute magnitudes of supernovae to distances measured by baryon acoustic oscillations (BAO). This ‘inverse distance ladder’ technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning $0.01\lt z\lt 1.13$ anchored to the recent baryon acoustic oscillation measurements from Dark Energy Spectroscopic Instrument (DESI) spanning $0.30 \lt z_{\mathrm{eff}}\lt 2.33$. To trace cosmology to $z=0$, we use the third-, fourth-, and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find $H_0=67.19^{+0.66}_{-0.64}\mathrm{~km} \mathrm{~s}^{-1} \mathrm{~Mpc}^{-1}$, in agreement with the value found by Planck without the need to assume Flat-$\Lambda$CDM. However, the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Investigation of secondary γ -ray angular distributions using the N 15 ( p , α 1 γ ) C * 12 reaction

The observation of secondary γ-rays provides an alternative method of measuring cross sections that populate excited final states in nuclear reactions. The angular distributions of these γ-rays also provide information on the underlying reaction mechanism. Despite a large amount of data of this type in the literature, publicly available R-matrix codes do not have the ability to calculate these types of angular distributions. In this report, the mathematical formalism derived in Brune and deBoer is implemented in the R-matrix code AZURE2 and calculations are compared with previous data from the literature for the 15 N( p, α1γ ) 12 C* reaction. In addition, new measurements, made at the University of Notre Dame Nuclear Science Laboratory using the Hybrid Array of Gamma Ray Detectors (HAGRiD), are reported that span an energy range from E p = 0.88 to 4.0 MeV. Excellent agreement between data and the phenomenological fit is obtained up to the limit of the previous fit at E p = 2.0 MeV and the R-matrix fit is extended from E x ≈ 13.5 up to E x ≈ 15.3 MeV, where 15 N+p and 12 C+α reactions are fit simultaneously for the first time. An excellent reproduction of the 15 N( p, α1γ ) 12 C* and 12 C(α, α) 12 C data is achieved, but inconsistencies and difficulty in fitting other data is encountered and discussed.

6 ≤ A ≤ 19↗
DOE OSTI2022

Isotherm Model for Moisture-Controlled CO 2 Sorption

Moisture-controlled sorption of CO 2 , the basis for moisture-swing CO 2 capture from air, is a novel phenomenon observed in strong-base anion exchange materials. Prior research has shown that Langmuir isotherms provide an approximate fit to moisture-controlled CO 2 sorption isotherm data. However, this fit still lacks a governing equation derived from an analytic model. As such, in this paper we derive an analytic form for an isotherm equation from a bottom-up approach, starting with a fundamental theory for an alkali liquid. In the range of interest relevant to CO 2 capture from air, an isotherm equation for an alkali liquid reduces to a simple analytic form with a single parameter, $K_{\text{eq}}$. In the limit $K_{\text{eq}} \gg 1$, a 2nd order approximation simplifies to a Langmuir isotherm that, however, deviates from experimental data. The isotherm theory for an alkali liquid has been generalized to a strong-base anion exchange material. In a strong-base anion exchange material, water concentration inside a sorbent, [H 2 O], is not large enough to be regarded as constant, which allows us to extend $K_{\text{eq}}$ to $K_{\text{eq(AEM)eff}} = K_{\text{eq(AEM)}} \times$ [H 2 O]$^{–n}$ according to the law of mass action. The final isotherm formula has been validated by experimental data from the literature. For a moisture-controlled CO 2 sorbent, $K_{\text{eq(AEM)eff}}$ varies significantly with moisture content of the sorbent. Depending on moisture level, the observed $K_{\text{eq(AEM)eff}}$ in a specific sorbent ranges from a few times to a few thousand times the value of $K_{\text{eq}}$ of a 2 mol L –1 alkali liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Automatic Multiple Experiment Simulation and Fitting (Ames-Fit)

AMES-Fit is a program used to automatically fit the multi-field solid-state NMR spectra of half-integer quadrupolar nuclei. Due to the high dimensional space, gradient algorithms have failed to address the fitting of such data, which is at present done manually. AMES-Fit diverges from these approaches by using an adaptive step size random search algorithm to fit the NMR spectra to consistently find the global best fit parameters.

Perras, Frederic↗
DOE OSTI2022

Flat-spectrum Radio Quasars and BL Lacs Dominate the Anisotropy of the Unresolved Gamma-Ray Background

We analyze the angular power spectrum (APS) of the unresolved gamma-ray background (UGRB) emission and combine it with the measured properties of the resolved gamma-ray sources of the Fermi-LAT 4FGL catalog. Our goals are to dissect the composition of the gamma-ray sky and to establish the relevance of different classes of source populations of active galactic nuclei in determining the observed size of the UGRB anisotropy, especially at low energies. We find that, under physical assumptions for the spectral energy distribution, i.e., by using the 4FGL catalog data as a prior, two populations are required to fit the APS data, namely flat-spectrum radio quasars at low energies and BL Lacs at higher energies. The inferred luminosity functions agree well with the extrapolation of the flat-spectrum radio quasar and BL Lac ones obtained from the 4FLG catalog. We use these luminosity functions to calculate the UGRB intensity from blazars, finding a contribution of 20% at 1 GeV and 30% above 10 GeV. Finally, bounds on an additional gamma-ray emission due to annihilating dark matter are also derived.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Characterization of the atomic-level structure of γ-alumina and (111) Pt/γ-alumina interfaces

The atomic-level structure of platinum/γ-alumina interfaces is characterized in a model system of dense γ-alumina embedded with faceted Pt NPs produced by implantation of platinum ions into sapphire followed by thermal annealing in air at 800 °C. Aberration-corrected scanning transmission electron microscopy (STEM) was used to collect atomic-resolution images, which are compared to STEM image simulations of two experimentally-based bulk models of γ-alumina by Smrčok et al. and Zhou and Snyder. A density functional theory (DFT) based model of (111) interfaces with different chemical terminations (O, Al 1 , Al 2 ) of the γ-alumina developed by Oware Sarfo et al. is also compared to experimental STEM data from Pt/γ-alumina interfaces. Further, the Smrčok γ-alumina model provides a better fit than the Zhou structure to the bulk of the γ-alumina. The oxygen-terminated Oware Sarfo model best fits the experimental data and is a very good model close to the interface. However, the fit of the interface model to the experimental data is poorer beyond the third atomic layer in the γ-alumina. This is attributed to compromises required in the design of the model to limit the cell size and computational time for DFT calculations. Understanding the accuracy and limits of the structural models of γ-alumina and Pt/γ-alumina interfaces is important to further the understanding of the structure/property relationships in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Outbreak analysis from case count and survey data

This code enables fitting a model of an infectious disease outbreak, customized for when available data is a combination of case counts and symptom survey responses. This is just a small bit of code being used only for a publication; we need to release it to fulfill a journal's requirement that our work be reproducible.

Goldberg, Emma↗
DOE OSTI2020

Measurement of 3-Flavour Neutrino Oscillation Parameters in the NOvA Experiment

NOvA is a long-baseline neutrino oscillation experiment consisting of two functionally identical tracking calorimeters, a Near and Far Detector, that measure neutrino interactions induced by the Fermi National Accelerator Laboratory’s NuMI beam at baselines of 1 km and 810 km, respectively. The NuMI beam can be configured to produce either a primary ν µ neutrino or $\bar{ν}$ µ anti-neutrino beam. Neutrino oscillations are observed and measured by the analysis of ν µ + $\bar{ν}$ µ disappearance and νe + $\bar{ν}$ e appearance in the beam, comparing the neutrino energy spectra in the Near and Far detectors means that neutrino oscillation parameters sin 2 θ 23 , |Δ$m^2_{32}$|, and δ CP can be constrained. This thesis presents the 2018 NOvA ν µ +$\bar{ν}$ µ disappearance, ν e +$\bar{ν}$ e appearance, and combined analyses using both neutrino and anti-neutrino data, where oscillation fits have been performed, where possible, on an event-by-event basis rather than on a bin-by-bin basis as has conventionally been used in NOvA oscillation analyses. This allows for better precision in applying both neutrino oscillation probabilities and systematic uncertainties. Furthermore, oscillation analyses for the disappearance, appearance, and combined channels are presented using an unbinned likelihood fit and compared with the equivalent binned χ 2 likelihood fit used in the standard analysis. The 14 ktonne detector equivalent beam exposures used for this thesis are 8.85 × 10 20 and 6.91 × 10 20 protons on target for neutrino and antineutrino data respectively, corresponding to 5 years of NOvA data taking. A combined ν µ + $\bar{ν}$ µ disappearance and ν e + $\bar{ν}$ e appearance fit to the Far Detector data, assuming normal mass ordering and using the event-by-event oscillation and unbinned fitting methodologies, produces oscillation parameter constraints of Δ$m^2_{32}$ = (2.50$^{+0.08}_{-0.06}$) × 10 -3 eV 2 , sin 2 θ 23 = 0.59$^{+0.02}_{-0.04}$, and δ CP = 0.72$^{+0.5}_{-0.9}$π.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Stochastic Models, Indices & Optimization Algorithms for Pricing & Hedging Reliability Risks in Modern Power Grids: Data Plan - Princeton

We collected and cleaned the synthetic grid data produced by NREL for the Texas and New York synthetic grids. We developed a high dimensional joint stochastic model for load at the zone level, and solar and wind power productions at the asset level, capturing the spatial and temporal dependencies between all the variables, and demonstrated how such a model could be fitted to historical data. We designed and implemented a simulation engine which can produce Monte Carlo scenarios for the hourly day-ahead values of load, and solar and wind power productions at the spatial and temporal resolutions of the historical data used to fit the model. Finally we developed an open-source Python package which can, from an input grid model, efficiently use forecasts and large numbers of Monte Carlo scenarios to provide unit commitment and economic dispatch for each of these scenarios. The high dimensional stochastic model and the subsequent Monte Carlo simulation engine were implemented in the package PGscen and the corresponding UC and ED optimization programs in the package Vatic.

14 SOLAR ENERGY↗
DOE OSTI2024

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗
DOE OSTI2022

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Fast and flexible analysis of direct dark matter search data with machine learning

We present the results from combining machine learning with the profile likelihood fit procedure, using data from the Large Underground Xenon (LUX) dark matter experiment. This approach demonstrates reduction in computation time by a factor of 30 when compared with the previous approach, without loss of performance on real data. We establish its flexibility to capture non-linear correlations between variables (such as smearing in light and charge signals due to position variation) by achieving equal performance using pulse areas with and without position-corrections applied. Its efficiency and scalability furthermore enables searching for dark matter using additional variables without significant computational burden. We demonstrate this by including a light signal pulse shape variable alongside more traditional inputs such as light and charge signal strengths. Furthermore, this technique can be exploited by future dark matter experiments to make use of additional information, reduce computational resources needed for signal searches and simulations, and make inclusion of physical nuisance parameters in fits tractable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Transverse momentum dependent feed-down fractions for bottomonium production

We extract transverse-momentum-dependent feed-down fractions for bottomonium production using a data-driven approach. We use data published by the ATLAS, CMS, and LHCb Collaborations for √𝑠 =7 TeV proton-proton collisions. Based on this collected data, we produce fits to the differential cross sections for the production of both 𝑆- and 𝑃-wave bottomonium states. Combining these fits with branching ratios for excited state decays from the Particle Data Group, we compute the feed-down fractions for both the ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) as a function of transverse momentum. Our results indicate a strong dependence on transverse momentum, which is consistent with prior extractions of the feed-down fractions. When evaluated at the average momentum of the states, we find that approximately 75% of ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) states are produced directly. Our results for the transverse-momentum-dependent feed-down fractions are provided in tabulated form so that they can be used by other research groups.

Astronomy & Astrophysics↗
DOE OSTI2025

Finite width effects in nonleptonic D -meson decays

Many analyses of two-body nonleptonic decays of D -mesons rely on flavor S U ( 3 ) symmetry relations and fits of experimental data of decays rates to extract the universal transition amplitudes. Such fits assume that the final state mesons are well-defined asymptotic states of QCD. We develop a technique to take into account the finite width effects of the final state mesons and study their effects on the extracted values of transition amplitudes. Published by the American Physical Society 2025

Kumar, Girish (ORCID:0000000160512495)↗
DOE OSTI2022

Substructure at High Speed. II. The Local Escape Velocity and Milky Way Mass with Gaia eDR3

Abstract Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be v esc = 445 − 8 + 25 km s −1 . A fit with a single component to the same data gives v esc = 472 − 12 + 17 km s −1 . Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of c 200 = 19 − 7 + 11 and a mass of M 200 = 4.6 − 0.8 + 1.5 × 10 11 M ⊙ , which is considerably lighter than previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗