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

Glass Property-Composition Models Update for use in Direct Feed High-Level Waste Flowsheet Development: EWG2.6

A set of preliminary glass property models and constraints were developed and augmented by models from literature for use in design of direct-feed high-level waste (DFHLW) glasses for flowsheet evaluation, testing, and design of the Waste Treatment and Immobilization Plant (WTP) high-level waste (HLW) Facility. These models and constraints are meant to be used as a place-holder while glass property-composition data gaps are filled and final plant operating models are developed. This report describes the motivation and intended use of the models, the compilation of data, model fitting and selection, methods to apply the models and constraints in glass design and offers example calculations demonstrating their intended use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Hierarchical Fractional Advection-Dispersion Equation (FADE) to Quantify Anomalous Transport in River Corridor over a Broad Spectrum of Scales: Theory and Applications

Fractional calculus-based differential equations were found by previous studies to be promising tools in simulating local-scale anomalous diffusion for pollutants transport in natural geological media (geomedia), but efficient models are still needed for simulating anomalous transport over a broad spectrum of scales. This study proposed a hierarchical framework of fractional advection-dispersion equations (FADEs) for modeling pollutants moving in the river corridor at a full spectrum of scales. Applications showed that the fixed-index FADE could model bed sediment and manganese transport in streams at the geomorphologic unit scale, whereas the variable-index FADE well fitted bedload snapshots at the reach scale with spatially varying indices. Further analyses revealed that the selection of the FADEs depended on the scale, type of the geomedium (i.e., riverbed, aquifer, or soil), and the type of available observation dataset (i.e., the tracer snapshot or breakthrough curve (BTC)). When the pollutant BTC was used, a single-index FADE with scale-dependent parameters could fit the data by upscaling anomalous transport without mapping the sub-grid, intermediate multi-index anomalous diffusion. Pollutant transport in geomedia, therefore, may exhibit complex anomalous scaling in space (and/or time), and the identification of the FADE’s index for the reach-scale anomalous transport, which links the geomorphologic unit and watershed scales, is the core for reliable applications of fractional calculus in hydrology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Condensed Phase Guerbet Reactions of Ethanol/Isoamyl Alcohol Mixtures

The self-condensation and cross-condensation reactions of ethanol and isoamyl alcohol are examined to better understand the potential routes to value-added byproducts from fuel ethanol production. Reactions have been carried out in both batch autoclave and continuous condensed-phase reactors using a lanthanum-promoted, alumina-supported nickel catalyst at near-critical condensed phase conditions. Analysis of multiple candidate kinetic models led to a Langmuir–Hinshelwood rate expression that is first-order in alcohol with water as the strongly adsorbed species. This model provides the best fit of data from both batch and continuous reactor experiments. Activation energies for primary condensation reactions increase as carbon chain lengths increase. Selectivities to higher alcohols of 94% and 87% for ethanol and isoamyl alcohol, respectively, were observed at different operating conditions.

Nezam, Iman (ORCID:0000000277598069)↗
DOE OSTI2021

Observations of the Bright Star in the Globular Cluster 47 Tucanae (NGC 104)

The Bright Star in the globular cluster 47 Tucanae (NGC 104) is a post-asymptotic giant branch (post-AGB) star of spectral type B8 III. The ultraviolet spectra of late-B stars exhibit myriad absorption features, many due to species unobservable from the ground. The Bright Star thus represents a unique window into the chemistry of 47 Tuc. We have analyzed observations obtained with the Far Ultraviolet Spectroscopic Explorer, the Cosmic Origins Spectrograph aboard the Hubble Space Telescope, and the Magellan Inamori Kyocera Echelle Spectrograph on the Magellan Telescope. By fitting these data with synthetic spectra, we determine various stellar parameters (T {sub eff} = 10,850 ± 250 K, logg=2.20±0.13) and the photospheric abundances of 26 elements, including Ne, P, Cl, Ga, Pd, In, Sn, Hg, and Pb, which have not previously been published for this cluster. Abundances of intermediate-mass elements (Mg through Ga) generally scale with Fe, while the heaviest elements (Pd through Pb) have roughly solar abundances. Its low C/O ratio indicates that the star did not undergo third dredge-up and suggests that its heavy elements were made by a previous generation of stars. If so, this pattern should be present throughout the cluster, not just in this star. Stellar-evolution models suggest that the Bright Star is powered by a He-burning shell, having left the AGB during or immediately after a thermal pulse. Its mass (0.54 ± 0.16M {sub ⊙}) implies that single stars in 47 Tuc lose 0.1–0.2 M {sub ⊙} on the AGB, only slightly less than they lose on the red giant branch.

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

A Swing of the Pendulum: The Chemodynamics of the Local Stellar Halo Indicate Contributions from Several Radial Merger Events

We find that the chemical abundances and dynamics of APOGEE and GALAH stars in the local stellar halo are inconsistent with a scenario in which the inner halo is primarily composed of debris from a single massive, ancient merger event, as has been proposed to explain the Gaia-Enceladus/Gaia Sausage (GSE) structure. The data contain trends of chemical composition with energy that are opposite to expectations for a single massive, ancient merger event, and multiple chemical evolution paths with distinct dynamics are present. We use a Bayesian Gaussian mixture model regression algorithm to characterize the local stellar halo, and find that the data are fit best by a model with four components. We interpret these components as the Virgo Radial Merger (VRM), Cronus, Nereus, and Thamnos; however, Nereus and Thamnos likely represent more than one accretion event because the chemical abundance distributions of their member stars contain many peaks. Although the Cronus and Thamnos components have different dynamics, their chemical abundances suggest they may be related. We show that the distinct low- and high-α halo populations from Nissen & Schuster are explained by VRM and Cronus stars, as well as some in situ stars. Because the local stellar halo contains multiple substructures, different popular methods of selecting GSE stars will actually select different mixtures of these substructures, which may change the apparent chemodynamic properties of the selected stars. We also find that the Splash stars in the Solar region are shifted to higher v $\phi$ and slightly lower [Fe/H] than previously reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Autodifferentiable Spectrum Model for High-dispersion Characterization of Exoplanets and Brown Dwarfs

We present an autodifferentiable spectral modeling of exoplanets and brown dwarfs. This model enables a fully Bayesian inference of the high-dispersion data to fit the ab initio line-by-line spectral computation to the observed spectrum by combining it with the Hamiltonian Monte Carlo in recent probabilistic programming languages. An open-source code, ExoJAX, developed in this study, was written in Python using the GPU/TPU compatible package for automatic differentiation and accelerated linear algebra, JAX. We validated the model by comparing it with existing opacity calculators and a radiative transfer code and found reasonable agreements for the output. As a demonstration, we analyzed the high-dispersion spectrum of a nearby brown dwarf, Luhman 16 A, and found that a model including water, carbon monoxide, and H 2 /He collision-induced absorption was well fitted to the observed spectrum (R = 10 5 and 2.28-2.30 μm). As a result, we found that T0 = 1295$^{-32}_{+35}$ K at 1 bar and C/O = 0.62 ± 0.03, which is slightly higher than the solar value. This work demonstrates the potential of a full Bayesian analysis of brown dwarfs and exoplanets as observed by high-dispersion spectrographs and also directly imaged exoplanets as observed by high-dispersion coronagraphy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reverse-mode differentiation in arbitrary tensor network format: with application to supervised learning.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations of tensor networks that may have arbitrary and unconventional network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve efficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sterile Neutrino Search with the NOvA Detectors

NOvA is a two-detector experiment designed to observe neutrino oscillations from νµ toνe flavor states. Additionally, the location of the NOvA Near Detector is well situated to additionally observe short base line oscillations to sterile neutrino states. We expand upon previous NOvA sterile neutrino searches by using a covariance matrix fitting technique. This analysis is the first with NOvA to directly include neutrino interactions in the Near Detector while fitting the data. The NOvA Near Detector’s greater statistics to better constrain our Far Detector uncertainties in addition to increasing the range of Δ m$^2_{41}$ to which we are sensitive.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Determination of Neutrino Oscillation Parameters through the Feldman-Cousins Method by the NOvA Experiment

The NOvA experiment presents new measurements of the neutrino oscillation parameters obtained through a fit to data from the one megawatt NuMI neutrino beam in the NOvA detectors. The analysis uses muon-neutrino disappearance and electron-neutrino appearance in both neutrino and antineutrino beam polarities. With the addition of $\sim$ 100%\) more neutrino-mode beam exposure over the previously reported results, this analysis employs the unified approach of Feldman and Cousins to determine the confidence level intervals for the oscillation parameters $\theta_{23}, \delta_{CP}$, and $\Delta m_{32}^2$ for both neutrino mass orderings.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Carrier Dynamics and Photoluminescence in Antimony-Doped CdTe and CdSeTe

Arsenic (As)-doped graded Cd(Se,Te) photovoltaics have achieved record efficiency of 22.4%, but further increases in efficiency have been limited by doping activation of only ~1% and radiative voltage losses of ~100 meV. Antimony (Sb) can also act as a p-type dopant of Cd(Se,Te) but has received relatively little attention. We measured carrier dynamics and photoluminescence of Sb-doped CdTe and CdSeTe films grown by vapor-transport deposition, with dopant activation up to ~20% Time resolved terahertz (TRTS) and photoluminescence (TRPL) data were fit by simulating the semiconductor equations, revealing bulk lifetimes of >20 ns in CdTe:Sb and >65 ns in CdSeTe:Sb and surface recombination velocities <1000 cm/s at alumina interfaces. We also qualitatively compare transients and PL spectra of doped and undoped films, with no apparent variation associated with hole concentrations up to 1015 cm-3. Preliminary variable temperature steady state photoluminescence (VT-SSPL) indicates dominant band-gap emission in CdSeTe:Sb and no evidence of potential fluctuations with Sb-doping up to this hole concentration.

14 SOLAR ENERGY↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2021

Abridged spectral matrix inversion: parametric fitting of X-ray fluorescence spectra following integrative data reduction

Recent improvements in both X-ray detectors and readout speeds have led to a substantial increase in the volume of X-ray fluorescence data being produced at synchrotron facilities. This in turn results in increased challenges associated with processing and fitting such data, both temporally and computationally. Herein an abridging approach is described that both reduces and partially integrates X-ray fluorescence (XRF) data sets to obtain a fivefold total improvement in processing time with negligible decrease in quality of fitting. The approach is demonstrated using linear least-squares matrix inversion on XRF data with strongly overlapping fluorescent peaks. This approach is applicable to any type of linear algebra based fitting algorithm to fit spectra containing overlapping signals wherein the spectra also contain unimportant (non-characteristic) regions which add little (or no) weight to fitted values, e.g. energy regions in XRF spectra that contain little or no peak information.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Challenges of standard halo models in constraining galaxy properties from cosmic infrared background anisotropies

The halo model, combined with halo occupation distribution (HOD) prescriptions, is widely used to interpret cosmic infrared background (CIB) anisotropies and extract physical information about star-forming galaxies and their connection to large-scale structures. Recent CIB-specific implementations of the halo model have adopted more physical parameterizations. However, the extent to which these models can reliably recover meaningful physical parameters remains uncertain. We assessed whether the current parameterization of CIB halo models is sufficient to recover astrophysical quantities, such as star formation efficiency, η(M h , z), and halo mass at which the peak of star formation efficiency occurs, M max , when fit to mock data. We also assessed whether discrepancies arise from assumptions about galaxy emission (the HOD ingredients) or from more fundamental components in the halo model, such as bias and matter clustering. We fit the M21 CIB HOD model, implemented within the halo model framework, to mock CIB power spectra and star formation rate density (SFRD) data generated from the SIDES-Uchuu simulation, and compared the best-fit parameters to the known simulation inputs. We then repeated the analysis using a simplified version of the simulation (SSU), explicitly designed to match the HOD assumptions. A detailed comparison of model and simulation outputs was carried out to trace the origin of observed discrepancies. While the M21 HOD model provides a good fit to the mock data, it failed to recover the intrinsic parameters accurately, particularly the halo mass at which star formation efficiency peaks. This mismatch persists even when fitting data generated with the same model assumptions. We find strong agreement (within 5%) in the emission-related components (SFRD, emissivity), but observe a scale- and redshift-dependent offset exceeding 20% in the two-halo term of the CIB power spectrum. This likely arises from limitations in the treatment of halo bias and matter clustering within the linear approximation. Additionally, incorporating scatter in the SFR–halo mass relation and the spectral energy distribution (SED) templates significantly affects the shot noise (∼50%), but has only a modest impact (less than 10%) on the clustered component. These results suggest that recovering physical parameters from CIB clustering requires improvements to the cosmological ingredients of the halo model framework, such as adopting scale-dependent halo bias and nonlinear matter power spectra in addition to careful modeling of emission physics.

cosmic background radiation↗
DOE OSTI2020

ARES I: WASP-76 b, A Tale of Two HST Spectra

We analyze the transmission and emission spectra of the ultra-hot Jupiter WASP-76 b, observed with the G141 grism of the Hubble Space Telescope’s (HST) Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauREx 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-gray clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyze this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-gray cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterization with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Failure in Confined Brazilian Tests on Sandstone

Strength of rocks in the confined tension region, where the minimum principal stress is tensile, has only infrequently been measured and is not well understood. Quasi-static confined Brazilian tests under a range of confining stresses (2.76 to 27.58 MPa) where used to determine the strength of sandstone in the confined tension region. The test results indicate that the strength in the confined tension region was a strong function of the intermediate principal stress: increasing the intermediate principal stress significantly increased the strength of the sandstone. The strength data were well fit by the Mogi–Coulomb criterion, which accounts for the intermediate principal stress. Unconfined Brazilian strength data were not well fit to the Mogi–Coulomb criterion derived from the confined Brazilian test data, consistent with a transition from tensile to shear processes dominating failure with increasing confining pressure. Observations of post-failure fracture surfaces reveal more indication of shear processes with increasing confining pressure. Numerical simulations from combined finite-discrete element method are compared to the experimental results and reflect similar conditions for failure compared to the experimental tests in the confined tension region.

42 ENGINEERING↗
DOE OSTI2023

Non-parametric projections of national income distribution consistent with the Shared Socioeconomic Pathways

Understanding and projecting income distributions within countries and regions is important to understanding consumption trends and the distributional consequences of climate impacts and responses. Several global, country-level projections of income distribution are available but most project only the Gini coefficient (a summary statistic of the distribution) or utilize the Gini along with the assumption of a lognormal distribution. We test the lognormal assumption and find that it typically underestimates income in the highest deciles and over-estimates it in others. We find that a new model based on two principal components of national time series data for income distribution provides a better fit to the data for all deciles, especially for the highest and lowest. We also construct a projection model in which the first principal component is driven by the Gini coefficient and the second captures deviations from this relationship. We use the model to project income distribution by decile for all countries for the five shared socioeconomic pathways. We find that inequality is consistently higher than projections based on the Gini and the lognormal functional form, with some countries reaching ratios of the highest to lowest income deciles that are almost three times their value using the lognormal assumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

First results from search for muon-neutrino disappearance at ICARUS

The ICARUS collaboration has completed a search for muon-neutrino disappearance in the context of the 3+1 sterile neutrino model using data collected with the Booster Neutrino Beam at Fermilab during 2022-2023 (ICARUS Run 2). Events are reconstructed with two different reconstruction frameworks and we select events with 1 muon, at least 1 proton, and no pions in the final state. A new fitting framework, called PROfit, was developed for this analysis and used for the final results presented here. As a single detector oscillation search this analysis is systematics limited, but the tools shown here will serve as a building block for future SBN searches where systematics will be constrained by the addition of a near detector. I will discuss the details of the analysis with a focus on the fitting framework, mock data studies, and the final fitting procedure during the data unboxing.

Larkin, Jacob [Rochester U.]↗