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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nuclear Global Talent Development Program (Final Report)

Lawrence Livermore National Laboratory (LLNL) have developed sorption database named as L-SCIE (LLNL Surface Complexation/Ion Exchange). This framework (L-SCIE) includes not only a community database of sorption data but also a fitting workflow to produce optimized surface complexation reaction constants. Through this fellowship program we have updated the L-SCIE database for selenium (Se) sorption and potentiometric titration data of iron oxides. In addition, the relevant surface complexation models (SCMs) have been updated and details in model construction are given in the following section.

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

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗
DOE OSTI2021

Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: II. Fully vibrationally-resolved electronic excitation of the isotopologues of H 2 (X 1 $Σ^{+}_{g}$)

Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of the isotopologues of molecular hydrogen (D 2 , T 2 , HD, HT, and DT) initially in the ground electronic state. We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to calculate cross sections from threshold to 500 eV for excitation of all bound vibrational levels and dissociative excitation of the B 1 $Σ^{+}_{g}$, C 1 Π u , EF 1 $Σ^{+}_{g}$, B' 1 $Σ^{+}_{g}$, GK 1 $Σ^{+}_{g}$, I 1 Π g , J 1 Δ g , D 1 Π u , H 1 $Σ^{+}_{g}$, b 1 $Σ^{+}_{g}$, c 3 Π u , a 3 $Σ^{+}_{g}$, e 3 $Σ^{+}_{u}$, d 3 Π u , h 3 $Σ^{+}_{g}$, g 3 $Σ^{+}_{g}$, i 3 Π g , and j 3 Δ g electronic states from all bound vibrational levels of the ground electronic (X 1 $Σ^{+}_{g}$) state. Including the previously-published MCCC e-H2 cross sections the data set contains cross sections for over 60,000 electronic and vibrational transitions. The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files. The data can also be downloaded from the MCCC database at mccc-db.org.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Assessment of a Detailed Biomass Pyrolysis Kinetic Scheme in Multiscale Simulations of a Single-Particle Pyrolyzer and a Pilot-Scale Entrained Flow Pyrolyzer

A detailed biomass pyrolysis kinetic scheme was assessed in the multiscale simulations of a single-particle pyrolyzer with slow pyrolysis and a pilot-scale entrained flow pyrolyzer with fast pyrolysis. The detailed kinetic scheme of biomass pyrolysis developed by the CRECK group consists of 32 reactions and 58 species. A multiscale simulation model was developed, where the CRECK kinetics was employed to simulate biomass pyrolysis reactions, a one-dimensional particle model was utilized to simulate the intraparticle transport phenomena, and the particle-in-cell (PIC) model was employed to simulate the hydrodynamics. The multiscale model was first applied to simulate a single-particle pyrolysis experiment. The simulation with nonisothermal particles matched the experimental data better than the simulation with isothermal particles. Then the multiscale model was applied to simulate the pilot-scale entrained flow pyrolyzer. In this case, the simulation with isothermal particles matched the experimental data better than the simulation with nonisothermal particles. The reason for this difference might be that the kinetics itself already partially included the intraparticle transport effect as it was fitted using both TGA data (slow pyrolysis of small size biomass) and fluidized bed data (fast pyrolysis of relatively large size biomass). This study provides some insights into biomass pyrolysis kinetics development and pyrolyzer multiscale simulation for a future study.

09 BIOMASS FUELS↗
DOE OSTI2024

Constraining models for the origin of ultra-high-energy cosmic rays with a novel combined analysis of arrival directions, spectrum, and composition data measured at the Pierre Auger Observatory

The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. The model takes into account a rigidity-dependent magnetic field blurring and an energy-dependent evolution of the catalog contribution shaped by interactions during propagation. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around 20° for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of 4.5σ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the γ-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Studying baryon acoustic oscillations using photometric redshifts from the DESI Legacy Imaging survey DR9

Context. The Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Survey DR9 (DR9 hereafter), with its extensive dataset of galaxy locations and photometric redshifts, presents an opportunity to study baryon acoustic oscillations (BAOs) in the region covered by the ongoing spectroscopic survey with DESI. Aims. We aim to investigate differences between different parts of the DR9 footprint. Furthermore, we want to measure the BAO scale for luminous red galaxies within them. Our selected redshift range of 0.6–0.8 corresponds to the bin in which a tension between DESI Y1 and eBOSS was found. Methods. We calculated the anisotropic two-point correlation function in a modified binning scheme to detect the BAOs in DR9 data. We then used template fits based on simulations to measure the BAO scale in the imaging data. Results. Our analysis reveals the expected correlation function shape in most of the footprint areas, showing a BAO scale consistent with Planck’s observations. Aside from identified mask-related data issues in the southern region of the South Galactic Cap, we find a notable variance between the different footprints. Conclusions. We find that this variance is consistent with the difference between the DESI Y1 and eBOSS data, and it supports the argument that that tension is caused by sample variance. Additionally, we also uncovered systematic biases not previously accounted for in photometric BAO studies. We emphasize the necessity of adjusting for the systematic shift in the BAO scale associated with typical photometric redshift uncertainties to ensure accurate measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

UNESE Data Analysis-Disko Elm Gas Transport Characteristics (NA-22 Quarterly Report)

Recently published papers have demonstrated the impact of sorption mechanisms on gas transport in the subsurface following an underground nuclear explosion. To include sorption processes in the multi-physics models of noble gas production and transport, researchers have conducted experiments of gas sorption in geologic materials and calibrated a linear (Henry’s law) sorption model for all gases in various geologic materials. As a part of Task #4, we requested the experimental data from Sandia Laboratories and derived a dual-mode (Henry’s absorption and Langmuir adsorption) sorption model. The derived model better describes the nonlinear sorption mechanism with much higher fitness with the experimental data. When gas concentrations are relatively low (in the conditions of field experiments or underground nuclear explosions) compared to those used sorption experiments, the nonlinear sorption model can be approximated as a linear sorption model. The equivalent Henry’s constant, which equals the product of Langmuir capacity and affinity, is systematically higher than the literature values. Further analyses and comparison of sorption models are planned in FY2020 Q4 and the sorption model will be coupled in models of nonisothermal multiphase transport for studying noble gas detectability.

58 GEOSCIENCES↗
DOE OSTI2020

Hyperon III: $$K^{-}p - \pi \varSigma $$ coupled-channel dynamics in the $$\varLambda (1405)$$ mass region

Abstract Low-energy data on the three charge states in $$\gamma p \rightarrow K^+(\varSigma \pi )$$ γ p → K + ( Σ π ) from CLAS at JLab, on $$K^-p\rightarrow \pi ^0\pi ^0\varLambda $$ K - p → π 0 π 0 Λ and $$\pi ^0\pi ^0\varSigma $$ π 0 π 0 Σ from the Crystal Ball at BNL, bubble chamber data on $$K^-p\rightarrow \pi ^-\pi ^+\pi ^{\pm }\varSigma ^{\mp }$$ K - p → π - π + π ± Σ ∓ , low-energy total cross sections on $$K^-$$ K - induced reactions, and data on the $$K^-p$$ K - p atom are fitted with the BnGa partial-wave-analysis program. We find that the data can be fitted well with just one isoscalar spin-1/2 negative-parity pole, the $$\varLambda (1405)$$ Λ ( 1405 ) , and background contributions. In a fit with one isocsalar state, the $$\varLambda (1405)$$ Λ ( 1405 ) structure can be determined as a dominantly SU(3) singlet state. A fit with two isoscalar singlet states, with imposed properties of the low-mass state, is, however, also not incompatible with data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

gaia: An R package to estimate crop yield responses to temperature and precipitation

gaia is an open-source R package designed to estimate crop yield shocks in response to annual weather variations and CO 2 concentrations at the country scale for 17 major crops. This innovative tool streamlines the workflow from raw climate data processing to projections of annual shocks to crop yields at the country level, using the response surfaces from an empirical econometric model developed and documented in Waldhoff et al. (2020), which leverages historical weather, CO 2 , and crop yield data for robust empirical fitting for 17 crops. gaia uses these response surfaces with monthly temperature and precipitation projections (e.g., from the Coupled Model Intercomparison Project Phase 6 (CMIP6) (O’Neill et al., 2016) climate data bias-adjusted and statistically downscaled by the ISIMIP3BASD approach (Lange, 2019) in the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) (Warszawski et al., 2014)) to project yield shocks that can be applied to agricultural productivity changes at the country level for use in multisectoral economic models. The historical and future projections use gridded, country-and-crop specific monthly growing season precipitation and temperature data, aggregated to the national level, and weighted by cropland area derived from the global Monthly Irrigated and Rainfed Crop Areas around the year 2000 (MIRCA2000) dataset (Portmann et al., 2010). These annual, country, and crop-specific yield shocks can be aggregated to different definitions of regions, crop commodities, and time periods, as needed by specific multisectoral economic models. gaia serves as a lightweight, powerful tool that can aid exploration of crop yield responses under a broad range of future climate projections, enhancing human-Earth system analysis capabilities.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Deriving spatially explicit direct and indirect interaction networks from animal movement data

Abstract Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems [“GPS”]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous‐time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions—individuals occurring at the same location, but at different times—while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease‐relevant interaction networks for two behaviorally differentiated species, wild pigs ( Sus scrofa ) that can host African Swine Fever and mule deer ( Odocoileus hemionus ) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30‐min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM‐Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine‐scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer–resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Detecting an underground tunnel by applying joint traveltime and waveform inversion

Underground tunnel and void detection is a challenging geophysical problem, and many methods have been proposed. Seismic techniques are promising because of the large seismic velocity contrasts between the air-filled void and the surrounding sediment and concrete. We apply a joint seismic traveltime and waveform inversion method to image a buried tunnel with concrete walls and a void space inside. The joint inversion images the top of the concrete tunnel as a high-velocity anomaly and the void space as a low-velocity anomaly. The location of the velocity anomalies predicted by the method agrees with the known location of the tunnel. As a comparison, the stand-alone full waveform inversion is also applied to the data. The first-arrival traveltime tomography shows weak nonlinearity but fails to image the hidden low-velocity layer. Full waveform inversion is able to image complex near-surface structures, but fitting waveforms may not honor the traveltime fit, especially when the data contain noise. Synthetic and real data tests show that the joint inversion method retains the advantages of both traveltime inversion and full waveform inversion and overcomes their respective drawbacks at the same time. The field example shows the joint inversion provides a better reconstruction of the high-velocity feature representing the top of the concrete wall in terms of its magnitude and location.

58 GEOSCIENCES↗
DOE OSTI2022

Improved limits on the tensor-to-scalar ratio using BICEP and P l a n c k data

In this study, we present constraints on the tensor-to-scalar ratio r using a combination of BICEP/Keck 2018 (BK18) and Planck PR4 data allowing us to fit for r consistently with the six parameters of the Λ CDM model. We discuss the sensitivity of constraints on r to uncertainties in the Λ CDM parameters as defined by the Planck data. In particular, we are able to derive a constraint on the reionization optical depth τ and thus propagate its uncertainty into the posterior distribution for r . While Planck sensitivity to r is slightly lower than the current ground-based measurements, the combination of Planck with BK18 and baryon-acoustic-oscillation data yields results consistent with r = 0 and tightens the constraint to r < 0.032 at 95% confidence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neural Networks for Supercapacitor Modeling

This software implements a one-dimensional model of a supercapacitor and calibrates it to experimental data. A system of coupled partial differential equation models the salt concentration and the ionic and electronic potentials. The term in the equation modeling the capacitance contains a neural network, whose parameters are optimized to fit the provided experimental data. The software is written in python using the open-source library JAX to calculate the gradients needed for the optimization.

Salazar de Troya, Miquel↗
DOE OSTI2024

Addendum to: Combined analysis of neutrino decoherence at reactor experiments

We update our analyses to constrain neutrino decoherence induced by wavepacket separation with RENO and Daya Bay data, now including the final data sets of the two experiments. We find that while the individual bounds from Daya Bay and RENO data improve relative to our original estimates, the combined fits are still dominated by KamLAND data and are only minimally improved.

de Gouvêa, André↗
DOE OSTI2024

In-situ FTIR Spectral Deconvolution [SWR-24-60]

In-situ FTIR Spectral Deconvolution contains scripts for the fitting of gaussians to convoluted FTIR spectra. It was specifically developed to handle a time-series of in situ ATR-FTIR spectra, split the spectra into its individual gaussians, and track those peaks as they shift in intensity and energy. Along with the actual fitting script, this repository contains several support scripts for data analysis and recall of prior fit results.

Tremolet de Villers, Bertrand↗
DOE OSTI2021

Using Bayesian Methodology to Estimate Liquefied Natural Gas Leak Frequencies

This analysis provides estimates on the leak frequencies of nine components found in liquefied natural gas (LNG) facilities. Data was taken from a variety of sources, with 25 different data sets included in the analysis. A hierarchical Bayesian model was used that assumes that the log leak frequency follows a normal distribution and the logarithm of the mean of this normal distribution is a linear function of the logarithm of the fractional leak area. This type of model uses uninformed prior distributions that are updated with applicable data. Separate models are fit for each component listed. Five order-of-magnitude fractional leak areas are considered, based on the flow area of the component. Three types of supporting analyses were performed: sensitivity of the model to the data set used, sensitivity of the leak frequency estimates to differences in the model structure or prior distributions, and sufficiency of sample sized used for convergence. Recommended leak frequency distributions for all component types and leak sizes are given. These leak frequency predictions can be used for quantitative risk assessments in the future.

03 NATURAL GAS↗
DOE OSTI2025

The track-length extension fitting algorithm for energy measurement of interacting particles in liquid argon TPCs and its performance with ProtoDUNE-SP data

This paper introduces a novel track-length extension fitting algorithm for measuring the kinetic energies of inelastically interacting particles in liquid argon time projection chambers (LArTPCs). The algorithm finds the most probable offset in track length for a track-like object by comparing the measured ionization density as a function of position with a theoretical prediction of the energy loss as a function of the energy, including models of electron recombination and detector response. The algorithm can be used to measure the energies of particles that interact before they stop, such as charged pions that are absorbed by argon nuclei. The algorithm's energy measurement resolutions and fractional biases are presented as functions of particle kinetic energy and number of track hits using samples of stopping secondary charged pions in data collected by the ProtoDUNE-SP detector, and also in a detailed simulation. Additional studies describe the impact of the dE/dx model on energy measurement performance. The method described in this paper to characterize the energy measurement performance can be repeated in any LArTPC experiment using stopping secondary charged pions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Uncertainty Quantification in GADRAS Inverse Modeling

The Gamma Detector Response and Analysis Software (GADRAS) package includes an inverse modeling tool that is helpful in identifying characteristics of unknown radioactive materials. Traditionally, uncertainties in this analysis were derived solely from measurement data quality and the fit of synthetic spectra. This paper aims to rigorously quantify additional sources of uncertainty, focusing on uncertainties arising from measurements being analyzed, Detector Response Function (DRF) characterization, and DRF extrapolation. Applying these findings to the BeRPBall benchmark data set, we demonstrated the impact of these uncertainties on plutonium and polyethylene estimates. The results underscore the importance of incorporating diverse uncertainty sources to enhance the accuracy and reliability of GADRAS’s inverse modeling capabilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗