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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Appraising constrained second-order power corrections in HQET with Λ b →Λ c Ι⁢ν

We derive the Λ b →Λ c form factors for the Standard Model and beyond at second order in heavy quark effective theory, applying the recently proposed residual chiral expansion to reduce the set of unknown subsubleading hadronic functions to a single, highly constrained function, that is fully determined by hadron mass parameters at zero recoil. We fit a form factor parametrization based on these results to all available lattice QCD predictions and experimental data. We find that the constrained and predictive structure of the form factors under the residual chiral expansion is in excellent agreement with lattice QCD predictions and experimental data, as well as prior heavy-quark-effective-theory-based fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Introducing a New Spitzer Master BLISS Map to Remove the Instrument Systematic Phase-curve-parameter Degeneracy, as Demonstrated by a Reanalysis of the 4.5 μm WASP-43b Phase Curve

While Spitzer Infrared Array Camera (IRAC) systematics are generally well understood, each data set can provide its own challenges that continue to teach us about the underlying functional form of these systematics. Multiple groups have analyzed the phase curves of WASP-43b with varying detrending techniques, each obtaining different results. In this work, we take another look at WASP-43b while exploring the degenerate relation between Bilinearly Interpolated Subpixel Sensitivity (BLISS) mapping, point response function (PRF)–FWHM detrending, and phase curve parameters. We find that there is a strong correlation between the detrending parameters in the two models, and best-fit phase curve amplitudes vary strongly when the data are temporally binned. To remove this degeneracy, we present a new Gaussian centroided intrapixel sensitivity map (hereafter fixed sensitivity map), generated using 3,712,830 exposures spanning 5 yr, for a variety of aperture sizes at 4.5 μm. We find evidence for time variability in the sensitivity at 3.6 μm and do not generate a visit-independent map for this channel. With the fixed 4.5 μm intrapixel sensitivity map, the best fits for WASP-43b no longer vary strongly with bin size and PRF–FWHM detrending is no longer required to remove correlated noise. For data sets that do not fall completely within the sweet spot, temporal binning should not be used in the analysis of Spitzer phase curves. We confirm nightside emission for WASP-43b with a disk-integrated nightside temperature of 806 ± 48 K at 4.5 μm. The 4.5 μm maps are available at github.com/kevin218/POET.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

A Denoising Autoencoder for Improved Kikuchi Pattern Quality and Indexing in Electron Backscatter Diffraction

The rapid collection and indexing of electron diffraction patterns as produced via electron backscatter diffraction (EBSD) has enabled crystallographic orientation and structural determination, as well as additional property-determining strain and dislocation density information with increasing speed, resolution, and efficiency. Pattern indexing quality is reliant on the noise of the collected electron diffraction patterns, which is often convoluted by sample preparation and data collection parameters. EBSD acquisition is sensitive to many factors and thus can result in low confidence index (CI), poor image quality (IQ), and improper minimization of fit, which can result in noisy datasets and misrepresent the microstructure. In an attempt to enable both higher speed EBSD data collection and enable greater orientation fit accuracy with noisy datasets, an image denoising autoencoder was implemented to improve pattern quality. Here, we show that EBSD data processed through the autoencoder results in a higher CI, IQ, and a more accurate degree of fit. In addition, using denoised datasets in HR-EBSD cross correlative strain analysis can result in reduced phantom strain from erroneous calculations due to the increased indexing accuracy and improved correspondence between collected and simulated patterns.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Joint Factorization of QCD and QED Radiation in Lepton-Hadron Scattering

The factorization theorem plays an important role in the analysis of high energy quantum chromodynamic (QCD) processes, separating the nonperturbative hadronic interaction into the universal parton distribution functions (PDFs) and fragmentation functions (FFs) and the process-dependent interactions into short distance perturbative calculations, with any interference power suppressed. With a virtual photon exchange, lepton-hadron deep inelastic scattering (DIS) provides an electromagnetic hard probe for the partonic structure of colliding hadrons and has played an important role in the development of QCD factorization. However, the collision induced QED radiation can change the momentum of the exchanged but unobserved virtual photon, making the photon-hadron frame, where the factorization formalism for DIS and semi-inclusive DIS (SIDIS) was derived, ill defined. A new analogous factorization approach has been introduced to separate the leading power process-independent QED radiative contributions to the single photon exchange by introducing lepton distribution functions (LDFs) and lepton fragmentation functions (LFFs), while process-dependent effects are perturbatively calculated with large logarithms removed [J. High Energ. Phys. 2021, 157 (2021)]. These LDFs and LFFs are considered global, as they appear in many different interactions, such as e+e-, DIS and SIDIS, so data from experiments can be used to fit and describe these functions across a wide range of lepton scattering. In this work, I will apply this new hybrid factorization approach to lepton-hadron DIS and SIDIS. For DIS, I derive the NLO short distance perturbative contribution to the cross section and demonstrate the effects the QED radiation has on the cross section using this approach using the CTEQ parameterization for the QCD functions. As part of the SIDIS analysis, I study the cross-section in two different kinematic regions: (1) the scattered lepton and observed hadron are not near back-to-back, and (2) they are close to back-to-back, where collinear QCD factorization works for (1) and TMD QCD factorization for (2) while collinear QED factorization works for both. As part of this work, I show the effects on the SIDIS cross section using fixed order calculations for the unpolarized structure function by first showing the effect of the radiative corrections on the main kinematic variables, especially how the internal transverse momentum is significantly correlated to the external angular dependence, and then the unpolarized structure function (or cross section) with matching between the descriptions for low and high transverse momentum. This work will impact the calculations for predictions for data from COMPASS and various Jefferson Lab experiments.

Cammarota, Justin [Univ. of Kentucky, Lexington, K↗

DOE OSTI2025

Reconstruction of 2D line-integrated electron density using angular filter refractometry and a fast marching Eikonal solver

Refraction of an optical probe beam by a plasma can be measured with angular filter refractometry (AFR), which produces an image of the beam’s 2D spatial profile that contains intensity contours corresponding to curves of constant refraction angle. Further analysis is required to reconstruct the underlying line-integrated electron density. Most prior efforts to calculate density from AFR data have been limited to 1D analysis or forward-fitting techniques. Here, in this paper, we detail the use of a fast-marching Eikonal solver to directly invert AFR data and obtain the full 2D line-integrated electron density. The analysis method is first verified with synthetic data and then applied to experimental measurements of single and colliding plasma plumes collected at the OMEGA EP Laser Facility. The calculated densities agree with 1D results and are shown to be consistent with the original AFR measurements via forward modeling. We also discuss ways to improve the precision of this technique.

McCluskey, B. [Princeton Univ., NJ (United States)↗

DOE OSTI2023

Optimal sensor placement for reconstructing wind pressure field around buildings using compressed sensing

Deciding how to optimally deploy sensors in a large, complex, and spatially extended structure is critical to ensure that the surface pressure field is accurately captured for subsequent analysis and design. In some cases, reconstruction of missing data is required in downstream tasks such as the development of digital twins. Here, this paper presents a data-driven sparse sensor selection algorithm, aiming to provide the most information contents for reconstructing aerodynamic characteristics of wind pressures over tall building structures parsimoniously. The algorithm first fits a set of basis functions to the training data, then applies a computationally efficient QR algorithm that ranks existing pressure sensors in order of importance based on the state reconstruction to this tailored basis. The findings of this study show that the proposed algorithm successfully re- constructs the aerodynamic characteristics of tall buildings from sparse measurement locations, generating stable and optimal solutions across a range of conditions. As a result, this study serves as a promising first step toward leveraging the success of data-driven and machine learning algorithms to supplement traditional genetic algorithms currently used in wind engineering.

42 ENGINEERING↗

DOE OSTI2022

Single Gaussian process method for arbitrary tokamak regimes with a statistical analysis

Abstract Gaussian process regression is a Bayesian method for inferring profiles based on input data. The technique is increasing in popularity in the fusion community due to its many advantages over traditional fitting techniques including intrinsic uncertainty quantification and robustness to over-fitting. This work investigates the use of a new method, the change-point method, for handling the varying length scales found in different tokamak regimes. The use of the Student’s t-distribution for the Bayesian likelihood probability is also investigated and shown to be advantageous in providing good fits in profiles with many outliers. To compare different methods, synthetic data generated from analytic profiles is used to create a database enabling a quantitative statistical comparison of which methods perform the best. Using a full Bayesian approach with the change-point method, Matérn kernel for the prior probability, and Student’s t-distribution for the likelihood is shown to give the best results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Geophysical and Environmental Monitoring Data, and Subsurface Flow Modelling Results for Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO

This dataset includes geoelectrical monitoring data acquired between October 2021 and November 2022, soil moisture and temperature data, groundwater data obtained from borehole SNIB covering the period from June 2021 to September 2022, and hydrological modelling results. The data were acquired to investigate how variations in bedrock type and topography, and vegetation cover control subsurface flow dynamics. To provide insights into the subsurface flow dynamics and their controls, a monitoring transect was installed at the Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO, measuring the spatio-temporal variations of soil moisture, soil and snow temperature, subsurface electrical resistivity variations, and groundwater dynamics. Field data are organized in a folder structure, with Electrical Resistivity Tomography (ERT) data being provided as one file per measurement, and data of the soil moisture and temperature sensors being provided as text files covering the entire monitoring period. The ‘Locations.csv’ file contains the location of all sensors, given in NAD83 – UTM Zone 13N. ERT monitoring data has been processed to filter data based on reciprocal errors (data with errors > 30% were removed), a linear error model was fitted to each survey, and to ensure a constant set of measurements for time-lapse inversion, filtered data were interpolated and assigned a 100% measurement error. Soil moisture and temperature data were acquired at 15 min intervals, and averaged to provide 1h data. Weather data and borehole data (groundwater depth, conductivity and temperature) were acquired at 30 min intervals, and are provided as daily measurements; all measurements are averaged, except of precipitation values, which are given as daily accumulation. The hydrological model was set up along the ERT monitoring transect, and net infiltration was used as surface boundary condition and derived from the weather data. Four different results are provided, (1) results for a parameterization using hydraulic permeability and porosity as derived from the ERT data through petrophysical relationships, and (2) three simplified model results, using 1 to 3 geological layers above the bedrock. Modelling was performed using PFLOTRAN, and for each model the PFLOTRAN input files are provided. The result files include weekly hydrological modelling results (e.g., saturation, velocities, pressures), as well as the model parameterization. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Streaming Generalized Canonical Polyadic Tensor Decompositions

In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced usability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Inclusive Search for Anomalous Single-Photon Production in MicroBooNE

The MicroBooNE experiment is an 85-ton active volume liquid argon time projection chamber (LArTPC) neutrino detector situated in the Fermilab Booster Neutrino Beam. One of the main goals for MicroBooNE is to investigate the “low energy excess” (LEE) of electromagnetic events observed by the MiniBooNE experiment by leveraging the unique capabilities of the LArTPC technology to distinguish photons from electron electromagnetic showers. In this poster, we present an inclusive search for anomalous production of single-photon events from neutrino interactions in the MicroBooNE experiment, motivated by the MiniBooNE LEE. The analysis uses the Wire-Cell reconstruction framework to select a sample of inclusive single-photon final-state interactions and leverages simultaneous measurements of sidebands to constrain signal and background predictions and reduce uncertainties. We perform a blind analysis using a dataset collected from February 2016 to July 2018 from the Booster Neutrino Beam (BNB) at Fermilab. In the full signal region, we observe agreement between the data and the prediction, with a goodness-of-fit p-value of 0.11. We then isolate a sub-sample of these events containing no visible protons, and observe 93 $/pm$ 22(stat.) $\pm$ 35(syst.) data events above prediction, corresponding to just above 2$\sigma$ local significance, concentrated at shower energies below 600 MeV.

Yandel, Erin [Los Alamos] (ORCID:0000000277123709)↗

DOE OSTI2022

Examining polymer‐protein biophysical interactions with small‐angle x‐ray scattering and quartz crystal microbalance with dissipation

Abstract Polymer‐protein hybrids can be deployed to improve protein solubility and stability in denaturing environments. While previous work used robotics and active machine learning to inform new designs, further biophysical information is required to ascertain structure–function behavior. Here, we show the value of tandem small‐angle x‐ray scattering (SAXS) and quartz crystal microbalance with dissipation (QCMD) experiments to reveal detailed polymer‐protein interactions with horseradish peroxidase (HRP) as a test case. Of particular interest was the process of polymer‐protein complex formation under thermal stress whereby SAXS monitors formation in solution while QCMD follows these dynamics at an interface. The radius of gyration ( R g ) of the protein as measured by SAXS does not change significantly in the presence of polymer under denaturing conditions, but thickness and dissipation changes were observed in QCMD data. SAXS data with and without thermal stress were utilized to create bead models of the potential complexes and denatured enzyme, and each model fit provided insight into the degree of interactions. Additionally, QCMD data demonstrated that HRP deforms by spreading upon surface adsorption at low concentration as shown by longer adsorption times and smaller frequency shifts. In contrast, thermally stressed and highly inactive HRP had faster adsorption kinetics. The combination of SAXS and QCMD serves as a framework for biophysical characterization of interactions between proteins and polymers which could be useful in designing polymer‐protein hybrids.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Dynamical Mass Estimates of the β Pictoris Planetary System through Gaussian Process Stellar Activity Modeling

Nearly 15 yr of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star β Pictoris. The δ Scuti pulsations of the star, which overwhelm planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than available in the literature to model the activity of the star. We show that a Gaussian process (GP) with a stochastically driven damped harmonic oscillator kernel can model the δ Scuti pulsations. It provides similar results to parametric models but with a simpler framework, using only three hyperparameters. It also enables us to model poorly sampled RV data that were excluded from previous analyses, hence extending the RV baseline by nearly five years. Altogether, the orbit and mass of both planets can be constrained from RV only, which was not possible with the parametric modeling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for β Pic b favor a low eccentricity of 0.029$_{−0.024}^{+0.061}$ and a relatively short period of 21.1$_{−0.8}^{+2.0}$ yr. The orbit of β Pic c is eccentric with 0.206$_{−0.063}^{+0.074}$ with a period of 3.36 ± 0.03 yr. We find model-independent masses of 11.7 ± 1.4 and 8.5 ± 0.5 M Jup for β Pic b and c, respectively, assuming coplanarity. The mass of β Pic b is consistent with the hottest start evolutionary models, at an age of 25 ± 3 Myr. A direct detection of β Pic c would provide a second calibration measurement in a coeval system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Binary classification with noisy quantum circuits and noisy quantum data

We study the effects of single-qubit noises in the quantum circuit and the corruption in the quantum training data to the performance of binary classification problem. We find that under the presence of errors, the measurement at a qubit is affected only by the noise in the corresponding qubit, and that the errors on other qubits do not affect this outcome. Furthermore, for the task where we fit a binary classifier using a quantum training data, we show that the noise in the data can work as a regularizer, implying that we can benefit from the noise in certain cases. We support our findings with simulations.

Lee, Yonghoon↗

DOE OSTI2024

Latest Three Flavor Neutrino Oscillation Results from the NOvA Experiment.

NOvA, is a two-detector, long-baseline neutrino oscillation experiment located at Fermilab, Batavia, IL, USA. It is designed primarily to constrain neutrino oscillation parameters such as the atmospheric mass squared splitting, $\Delta m^2_{32}$, the mixing angle, $\theta_{23}$, neutrino mass hierachy, and the CP-violating phase, $\delta_{CP}$, using $\nu_\mu \ (\bar{\nu}_\mu)$ disappearance and $\nu_e \ (\bar{\nu}_e)$ appearance data. NOvA receives a high purity instense beam of neutrinos and anti-neutrinos from Fermilab's Neutrinos at Main Injector (NuMI) beamline. The NOvA near detector (ND) observes un-oscillated $\nu_\mu \ (\bar{\nu}_\mu)$ and beam $\nu_e \ (\bar{\nu}_e)$ events, while the far detector (FD), which is situated 810 km away from the near detector, records un-oscillated $\nu_\mu \ (\bar{\nu}_\mu)$ and oscillated $\nu_e \ (\bar{\nu}_e)$ events. The best fit values of the oscillation parameters are extracted by performing a joint fit of the far detector predicted events to data from $\nu_\mu \to \nu_\mu$, $\bar{\nu}_\mu \to \bar{\nu}_\mu$, $\nu_\mu \to \nu_e$, and $\bar{\nu}_\mu \to \bar{\nu}_e$ oscillation channels. We present latest joint fit results from NOvA based on a neutrino-beam exposure of $27.24 \times 10^{20}$ POT and an anti-neutrino beam exposure of $12.55\times 10^{20}$ POT and improved simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Latest Three Flavor Neutrino Oscillation Results from the NOvA Experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Fundamental differences in the properties of red and blue quasars: measuring the reddening and accretion properties with X-shooter

ABSTRACT We have recently found fundamental differences in the radio properties of red quasars when compared to typical blue quasars. In this paper, we use data from the X-shooter spectrograph on the Very Large Telescope, providing spectral coverage from ∼3000–$25\,000\,$ Å, of a sample of 40 red and blue luminous quasars at 1.45 < z < 1.65 to explore the connections between the radio, emission-line, and accretion-disc properties. We fit various dust-extinction curves to the data and find that dust reddening can fully explain the observed colours for the majority of the red quasars in our sample, with moderate extinctions ranging from AV ∼ 0.06–0.7 mag. We confront our spectra with a simple thin accretion-disc model and find this can describe the continua of both the blue and red quasars, once corrected for dust extinction; we also find no significant differences in the accretion properties. We detect ionized outflows in a number of red and blue quasars, but do not find any significant evidence that they are more prevalent in the red quasar population. Overall our findings imply that the radio emission is more closely connected to circumnuclear/ISM opacity rather than accretion disc or outflow differences.

Fawcett, V. A. (ORCID:000000031251532X)↗

DOE OSTI2022

Uptake and speciation of trace metal inputs to Wetland Soils from Illinois and South Carolina and Stream Sediments from Tennessee

Metals occur in all ecosystems, although their concentrations vary depending on their natural geologic conditions and surrounding human activities. In addition to metals intrinsically present from the geology of a particular system, metals can enter environmental systems from a wide variety of natural and anthropogenic sources, such as sediment re-suspension, mining operations, industrial processes, agricultural activities, and atmospheric deposition. While metals can be toxic at high concentrations, some metals serve as essential micronutrients for biogeochemical processes. Metal transport and availability in engineered and natural water systems depend on processes of adsorption/desorption, oxidation/reduction, dissolution/precipitation, and ligand complexation. Insights into the speciation of metals and their bioavailability will also help advance understanding of the roles of metals in the biogeochemical cycling of nutrients. We conducted batch experiments under anoxic conditions on soils and sediments collected from three different natural aquatic systems to understand their response to influxes of dissolved Cu, Ni and Zn. While soils and sediments from all sites could strongly bind added trace metals, there were substantial differences in trace metal uptake trends between different sites, especially for Cu. There was no distinct correlation between trace metal uptake and the total organic matter, iron, and sulfur content present in the samples. X-ray absorption spectroscopy indicated that the speciation of the freshly added metals taken up by the solids differs substantially from the speciation of the metals originally present in unamended samples. Cu sulfides dominated speciation at low loadings (1 µmol/g), whereas complexation to thiol groups and formation of metallic Cu governed speciation at high loadings (10 µmol/g). For Ni and Zn, adsorption to mineral surfaces and organic matter governed their speciation in materials from most sites. This study suggests that the background speciation of metals in natural aquatic systems is a poor predictor of the speciation and lability of metals introduced to terrestrial aquatic systems from anthropogenic or natural processes. Our findings imply that geochemical processes controlling trace metal speciation may vary considerably with metal loading in different natural systems. Data are provided for marsh wetland soils at Argonne National Laboratory (Marsh 1 and Marsh 2), riparian wetland soils (Riparian 1 and Riparian 2) in the Tims Branch watershed at Savannah River National Laboratory, and stream bed sediments (Stream 1 and Stream 2) from East Fork Poplar Creek near Oak Ridge National Laboratory. Data package includes the results of trace metal uptake experiments conducted on the selected sites for determining their capacity to immobilize metals under different loadings. XANES spectra for Cu, Ni, ad Zn at different loadings are included in the package. The abundance of different trace metal species obtained using linear combination fitting in ATHENA are also included in the data.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Pursuit of paired dijet resonances in the Run 2 dataset with ATLAS

New particles with large masses that decay into hadronically interacting particles are predicted by many models of physics beyond the Standard Model. A search for a massive resonance that decays into pairs of dijet resonances is performed using 140 fb –1 of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. Resonances are searched for in the invariant mass of the tetrajet system, and in the average invariant mass of the pair of dijet systems. A data-driven background estimate is obtained by fitting the tetrajet and dijet invariant mass distributions with a four-parameter dijet function and a search for local excesses from resonant production of dijet pairs is performed. No significant excess of events beyond the Standard Model expectation is observed, and upper limits are set on the production cross sections of new physics scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗