Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 505 records · Page 28

Strangeness-correlations on the pseudocritical line in ( 2 + 1 )-flavor QCD

We present some lattice QCD results on first ( χ 1 i ) and second ( χ 2 i ) cumulants of and correlations ( χ 11 i j ) among net baryon-number ( B ), strangeness ( S ) and electric charge ( Q ) along the pseudocritical line [ T p c ( μ B ) ] in the temperature ( T )–baryon chemical potential ( μ B ) phase diagram of ( 2 + 1 )-flavor QCD. We point out that violations of sum rules among second order cumulants, which hold in the isospin symmetric limit of vanishing electric charge chemical potential, are small along the T p c ( μ B ) for the entire range of μ B covered in the RHIC beam energy scan. For the strangeness neutral matter produced in heavy-ion collisions this leads to a close relation between χ 11 B S and χ 11 Q S . We compare lattice QCD results for χ 11 B S / χ 2 S along the T p c ( μ B ) line with preliminary experimental measurements of χ 11 B S / χ 2 S for collision energies 7.7 GeV ≤ s N N ≤ 62.4 GeV . While we find good agreements for s N N ≥ 39 GeV , differences are sizeable at smaller values of s N N . Moreover, we compare lattice QCD results for the ratio of the strangeness ( μ S ) to baryon ( μ B ) chemical potentials, which define a strangeness neutral system with fixed electric charge to baryon number density, with experimental results obtained by the STAR collaboration for μ S / μ B using strange baryon yields on the freeze-out line. Finally, we determine the baryon chemical potential at the freeze-out ( μ B f ) by comparing χ 1 B / χ 2 B along the T p c ( μ B ) with the experimentally measured net-proton cumulants χ 1 p / χ 2 p . We find that { μ B f , T p c ( μ B f ) } are consistent with the freeze-out parameters of the statistical-model fits to experimentally measured hadron yields for s N N ≥ 11.5 GeV . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SPT clusters with DES and HST weak lensing. II. Cosmological constraints from the abundance of massive halos

We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel’dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range 0.25–1.78 over a total sky area of 5200 deg 2 . We use DES Year 3 weak-lensing data for 688 clusters with redshifts 𝑧 < 0.95 and HST weak-lensing data for 39 clusters with 0.6 < 𝑧 < 1.7. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable-mass relation without having to make strong assumptions about, e.g., the hydrodynamical state of the clusters. For a flat Λ⁢ CDM cosmology, and marginalizing over the sum of massive neutrinos, we measure Ω m = 0.286 ± 0.032, 𝜎 8 = 0.817 ± 0.026, and the parameter combination 𝜎 8 ⁢(Ω m /0.3) 0.25 = 0.805 ± 0.016. Our measurement of 𝑆 8 ≡ 𝜎 8 ⁢$\sqrt{Ω_{m}/0.3}$ = 0.795 ± 0.029 and the constraint from Planck CMB anisotropies (2018 TT, TE, EE+lowE) differ by 1.1⁢𝜎. In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses ∑𝑚 𝜈 < 0.18 eV. When additionally allowing the dark energy equation of state parameter 𝑤 to vary, we obtain 𝑤 = −1.45 ± 0.31 from our cluster-based analysis. In combination with Planck data, we measure 𝑤 =−1.3⁢4$^{+0.22}_{−0.15}$, or a 2.2⁢𝜎 difference with a cosmological constant. We use the cluster abundance to measure 𝜎8 in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the Λ⁢ CDM model fit to Planck primary CMB data.

79 ASTRONOMY AND ASTROPHYSICS↗

Light and Strange Vector Resonances from Lattice QCD at Physical Quark Masses

We present the first ab initio calculation at physical quark masses of scattering amplitudes describing the lightest pseudoscalar mesons interacting via the strong force in the vector channel. Using lattice quantum chromodynamics, we postdict the defining parameters for two short-lived resonances, the ρ(770) and K*(892) which manifest as complex energy poles in ππ and Kπ scattering amplitudes, respectively. The calculation proceeds by first computing the finite-volume energy spectrum of the two-hadron systems and then determining the amplitudes from the energies using the Lüscher formalism. The error budget includes a data-driven systematic error, obtained by scanning possible fit ranges and fit models to extract the spectrum from Euclidean correlators, as well as the scattering amplitudes from the latter. The final results, obtained by analytically continuing multiple parametrizations into the complex energy plane, are M ρ = 796(5)(50)MeV, Γ ρ = 192(10)(31) MeV, M K* = 893(2)(54) MeV, and Γ K = 51(2)(11) MeV, where the subscript indicates the resonance and M and Γ stand for the mass and width, respectively, and where the first bracket indicates the statistical and the second bracket the systematic uncertainty.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ultraviolet to far-infrared dielectric function of n -doped cadmium oxide thin films

Spectroscopic ellipsometry and Fourier transform infrared spectroscopy were applied to extract the ultraviolet to far-infrared (150 – 33333 cm -1 ) complex dielectric functions of high-quality, sputtered indium-doped cadmium oxide (In:CdO) thin crystalline films on MgO substrates possessing carrier densities (N d ) ranging from 1.1 x 10 19 cm -1 to 4.1 x 10 20 cm -1 . A multiple oscillator fit model was used to identify and analyze the three major contributors to the dielectric function and their dependence on doping density: interband transitions in the visible, free-carrier excitations (Drude response) in the near- to far-infrared, and IR-active optic phonons in the far-infrared. More specifically, values pertinent to the complex dielectric function such as the optical band gap (E g ), are shown here to be dependent upon carrier density, increasing from approximately 2.5–3 eV, while the high-frequency permittivity (ε ∞ ) decreases from 5.6 to 5.1 with increasing carrier density. The plasma frequency (ω p ) scales as $\sqrt{N_d}$, resulting in ω p values occurring within the mid- to near-IR, and the effective mass (m *) was also observed to exhibit doping density-dependent changes, reaching a minimum of 0.11 m o in unintentionally doped films (1.1 x 10 19 cm -1 ). Good quantitative agreement with prior work on polycrystalline, higher-doped CdO films is also demonstrated, illustrating the generality of the results. The analysis presented here in this paper will aid in predictive calculations for CdO-based next-generation nanophotonic and optoelectronic devices, while also providing an underlying physical description of the key properties dictating the dielectric response in this atypical semiconductor system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Construction of the damped Ly⁢𝛼 absorber catalog for DESI DR2 Ly⁢𝛼 BAO

We present the Damped Ly⁢𝛼 Toolkit for automated detection and characterization of damped Ly⁢𝛼 absorbers (DLAs) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance (Δ⁢𝜒$^{2}_{𝑟}$ >0.03), the technique achieves an estimated 80% purity and 79% completeness when evaluated on simulated spectra with S/N>2 that are free of broad absorption lines (BALs). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21 719 were found in S/N>2 spectra with predicted log 10 ⁡(𝑁 𝙷𝙸 )>20.3 and detection significance Δ⁢𝜒$^{2}_{𝑟}$ >0.03. We compare the Damped Ly⁢𝛼 Toolkit to two alternative DLA finders based on a convolutional neural network and Gaussian process models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly⁢𝛼 forest baryon acoustic oscillation measurement. The combined catalog contains 41 152 candidate DLAs with log 10 ⁡(𝑁 𝙷𝙸 )>20.3 from quasar spectra with S/N>2. We estimate this sample to be approximately 85% pure and 79% complete when BAL quasars are excluded.

79 ASTRONOMY AND ASTROPHYSICS↗

Storms Are an Important Driver of Change in Tropical Forests

Tropical forest dynamics and composition have changed over recent decades, but the proximate drivers of these changes remain unclear. Investigations into these trends have focused on increasing drought stress, CO 2 , temperature, and fires, whereas convective storms are generally overlooked. We argue that existing literature provides clear support for the importance of storms as drivers of forest change. We reanalyze the largest plot-based study of tropical forest carbon dynamics to show that lightning frequency—an indicator of storm activity—strongly predicts forest carbon storage and residence time, and its inclusion improves model fit and weakens evidence for the effects of high temperatures. Convective storm activity has increased 5%–25% per decade over the past half century. Extrapolating from historic trends, we estimate that storms likely contribute ca. 50% of the reported increases in biomass mortality across Amazonia, with all realistic combinations of assumptions indicating a possible range of 12%–118%. Spatial variation in storm activity shows weak relationships with drought, demonstrating that forests can experience high drought stress, high storm activity, or both. Accordingly, we hypothesise that convective storms are among the most important drivers of tropical forest change, and as such, they require significant research investment to avoid misguiding science, policy, and management.

biomass carbon↗

Physics-Based Continuum Modeling for an Aqueous Rechargeable Zn/MnO 2 Battery

This study introduces a framework for modeling the aqueous Zn/MnO 2 rechargeable battery. A reaction system and a physics-based continuum model are proposed based on two reaction types, one involving insertion and the second related to dissolution and deposition of a solid reaction product. The model, fitted to empirical data, predicts voltage behavior and capacity limitations during cycling, identifying electrolytic zinc depletion as a limiting mechanism, depending on the original cell construction. The research suggests the need for further material characterization and reaction analysis, which will advance our understanding and facilitate the development of grid-scale energy storage solutions.

Electrochemistry↗

Generalized Tensor-on-Tensor Regression (GToTR)

SAND2026-23069O Generalized Tensor-on-Tensor Regression (GToTR) is a Python-based tool for conducting generalized tensor-on-tensor regression. It provides Canonical Polyadic (CP)-based generalized tensor regression models, support for generalized linear model-like families and links, alternating-optimization model fitting methods, and a standard statistics software interface. The tool supports tensor-valued responses and covariates using the open-source Python Tensor Toolbox (pyttb) software package. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Dunlavy, Daniel [Sandia National Lab. (SNL-CA), Li↗

Observation of spin‑wave altermagnetic splitting in MnF2

-Contents of the Data Repository: - Polarized neutron diffraction data acquired in the (HK0) scattering plane. - Inelastic neutron scattering (INS) data from both unpolarized and polarized measurements with an incident neutron energy of Ei =9 meV.- - Reduced multidimensional single-crystal datasets (MDE) and corresponding S(Q,ω) slices used to generate all figures presented in the manuscript. - Julia source code and supporting input files used for spin-wave calculations, model fitting, and simulation of neutron scattering intensity maps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of burning on soil properties and respiration; Wood Buffalo National Park, Canada; June 2019

In this study, we used burn simulations on boreal forest soil under high and low moisture conditions to investigate the effect of burning on soil temperature, soil properties such as pH, total C and total C, and soil respiration. Location and site characteristic data of 19 sampling locations within Wood Buffalo National Park, Canada are included in "sampling-site-characteristics.csv". File "Temp-all-burns.csv" contains temperature data collected 2 to 9 cm below the soil surface during 2 min. heat treatment (i.e. burn simulation) and the following 5-24 hrs. Following the burn simulations, total respiration from burned and unburned soil was collected via KOH base traps for up to 6 months ("raw-conductivity-readings.csv"). The metadata file ("all-sample-metadata.csv") contains measurements of soil pH, total C, and total N of O horizon and underlying mineral soil as well as model fits for respiration data.

54 ENVIRONMENTAL SCIENCES↗

Linear viscoelastic measurements of commerical therapy putties for design and intuition

This dataset complements the publication Mapping Linear Viscoelasticity for Design and Tactile Intuition (Corman & Ewoldt, 2019). It contains rigorous shear rheometry measurements for six commercial silicone putties, including creep recovery, SAOS, LAOS, and stress relaxation. Additional, unpublished data was generated as part of a graduate rheology course taught by Randy Ewoldt at University of Illinois Urbana-Champaign. Together, these datasets provide a comprehensive example of linear viscoelasticity of a commercial material. The collection is intended both as a pedagogical resource for developing physical intuition in linear viscoelasticity and as a reference dataset for data interpretation, model fitting, and comparative rheological analysis.

C. Marsh, Maxwell↗

Uncertainty Quantification in High Explosives Equations of State [Slides]

In order to fit models to data, we need a way to simulate the data given a model. Our projects involved simulators for the Manganin gauge experiment and for cylinder tests, which we use to calibrate models to experiments. We call the parameters we seek to optimize the degrees of freedom (DOF). We use a Bayesian approach with its associated UQ benefits.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Calibrating the SPECTACULAR constitutive model using legacy Sandia data for two filled epoxy systems: 828/CTBN/DEA/GMB and 828/DEA/GMB

The SPECTACULAR model is a development extension of the Simplified Potential Energy Clock (SPEC) model. Both models are nonlinear viscoelastic constitutive models used to predict a wide range of time-dependent behaviors in epoxies and other glass-forming materials. This report documents the procedures used to generate SPECTACULAR calibrations for two particulate-filled epoxy systems, 828/CTBN/DEA/GMB and 828/DEA/GMB. No previous SPECTACULAR or SPEC calibration exists for 828/CTBN/DEA/GMB, while a legacy SPEC calibration exists for 828/DEA/GMB. To generate the SPECTACULAR calibrations, a step-by-step procedure was executed to determine parameters in groups with minimal coupling between parameter groups. This procedure has often been deployed to calibrate SPEC, therefore the resulting SPECTACULAR calibration is backwards compatible with SPEC (i.e. none of the extensions specific to SPECTACULAR are used). The calibration procedure used legacy Sandia experimental data stored on the Polymer Properties Database website. The experiments used for calibration included shear master curves, isofrequency temperature sweeps under oscillatory shear, the bulk modulus at room temperature, the thermal strain during a temperature sweep, and compression through yield at multiple temperatures below the glass transition temperature. Overall, the calibrated models fit the experimental data remarkably well. However, the glassy shear modulus varies depending on the experiment used to calibrate it. For instance, the shear master curve, isofrequency temperature sweep under oscillatory shear, and the Young's modulus in glassy compression yield values for the glassy shear modulus at the reference temperature that vary by as much as 15 %. Also, for 828/CTBN/DEA/GMB, the temperature dependence of the glassy shear modulus when fit to the Young's modulus at different temperatures is approximately four times larger than when it is determined from the isofrequency temperature sweep under oscillatory shear. For 828/DEA/GMB, the temperature dependence of the shear modulus determined from the isofrequency temperature sweep under oscillatory shear accurately predicts the Young's modulus at different temperatures. When choosing values for the shear modulus, fitting the glassy compression data was prioritized. The new and legacy calibrations for 828/DEA/GMB are similar and appear to have been calibrated from the same data. However, the new calibration improves the fit to the thermal strain data. In addition to the standard calibrations, development calibrations were produced that take advantage of development features of SPECTACULAR , including an updated equilibrium Helmholtz free energy that eliminates undesirable behavior found in previous work. In addition to the previously mentioned experimental data, the development calibrations require data for the heat capacity during a stress-free temperature sweep to calibrate thermal terms.

36 MATERIALS SCIENCE↗

Finding my drumbeat: applying lessons learned from Remo Ruffini to understanding astrophysical transients [Slides]

Lessons Learned: Bandwagon science can miss key physics; Bandwagon scientists cling to their bandwagon, oftentimes denying fundamental physics; and, Better approach: Be aware of the physical limitations of a paradigm. What assumptions have been made? How can this be different in Nature?; Be cognizant of alternative models; Encourage the development of those models; and, Although validation is important, fundamental physics forms the basis for Occam’s razor, not whether a toy model fits the data well.

79 ASTRONOMY AND ASTROPHYSICS↗

An Applied Strategy for Using Empirical and Hybrid Models in Online Monitoring

The monitoring of plant equipment for failure prediction is one of the key contributors to operation and maintenance (O&M) costs for a nuclear power plant (NPP) because O&M monitoring depends on labor-intensive activities that are required to meet high equipment reliability standards. These activities rely primarily on humans for information gathering, condition diagnosis, and predictive analysis. Online monitoring aims to automate these activities by relying on sensors to replace human information gathering and machine learning to replace human analysis and decision making. To facilitate automated monitoring, a systematic strategy for anomaly detection is needed to optimally use the available sensor data, empirical models, and physics-supported models. This strategy is essential to provide credible reasoning on why and when an empirical (i.e., purely data-driven) versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. The extant methods usually adopt an ad hoc trial-and-error approach that, in addition to being time-consuming and costly, is also highly subjective; it is impacted by the background and the skill set of the personnel making the decisions. Thus, such an approach cannot guarantee an optimum outcome. This represents the motivation of the current research effort, which is focused on devising a scientifically supported strategy for the optimum selection of anomaly detection methods. This report presents a detailed assessment of the main anomaly detection techniques within the empirical or hybrid method streams. Empirical methods include pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The listed techniques within these two streams represent the vast majority of techniques performed for anomaly detection. Using the techniques as outcomes, a strategy was developed to enable a systematic decision-making process to lead to one of these techniques. The strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability and suitability for training and testing, cause-effect, entropy inference, and model fitting. Each of these decision points in the strategy is explained in detail in this report with examples, along with the scientific basis behind the decisions and outcomes in common and simplified terminology. The strategy is developed for use by any NPP staff with basic engineering or science knowledge. A user-friendly graphical state flow diagram was also developed as a visual presentation of the strategy. The strategy was tested and demonstrated through two pilot projects for the application of anomaly detection at an NPP. Each pilot had two use cases: an initial case in which certain decisions were made that resulted in one or more empirical techniques and a revised use case where one or more key decisions were modified resulting in using a set of hybrid methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Yttrium solid phase equation of state with uncertainty quantification

We discuss the development of an equation of state (EOS) for solid yttrium up to roughly 80 GPa. The EOS makes use of different experimental datasets, including measurements of the isobaric density, isobaric specific heat, room temperature isotherm, and principal shock Hugoniot. The fitting procedure is carried out using Markov Chain Monte Carlo (MCMC), where we fit model parameters for both the cold curve and ion thermal models simultaneously. The results show close agreement with experimental data and provide meaningful uncertainty estimates on the model parameters. This work serves as a first step towards a multiphase EOS for yttrium, which will include higher pressure solid phases (> 80 GPa), as well as modeling of the liquid phase.

36 MATERIALS SCIENCE↗

Uncertainty-Driven Rapid Thermodynamic Assessment of Nb-Ta-Zr System and Effects of C impurities (L25GF9298S): Annual Progress Report

An integrated computational materials engineering (ICME) method is in development for rapid thermodynamic experimental investigation and high-fidelity computational modeling of refractory multi-principal element alloys (RMPEAs). These ultra-high-temperature (UHT) alloys are of interest for structural applications in extreme environments, but deficiency of reliable data, especially melting temperatures, impedes the prediction of alloys with favorable properties. The method leverages UHT capabilities and computational expertise of LLNL’s Materials Science Division and the McCormack Lab’s UHT conical nozzle levitation (CNL) system to iteratively map the Nb-Ta-Zr phase space, with focus on the liquidus surface, through targeted experiments selected by quantifying uncertainty in the thermodynamic model fitting parameters. This method will reduce the time to map uncharted RMPEA phase space and thereby accelerate discovery and development of advanced materials for applications in extreme environments.

36 MATERIALS SCIENCE↗

Variation in surface properties as error sources in dynamic temperature measurements

Most temperature measurement methods rely upon the assumption of a single uniform source as part of the analysis in inferring an observed temperature T from the surface of a material. This is true for the radiometric and spectroscopic methods we rely upon in shock physics. However, there are obvious cases where this is a poor assumption that can cause large errors in analysis, such as the measurement of emission from a surface composed of more than one material, such as a rusted surface of Fe and Fe 2 O 3 , or a multi-material surface that could be produced through alloying methods (e.g., mokume gane) or 3D printing. Here we perform an initial examination of the effects of varied surface properties upon the analysis using idealized test data to constrain when we may have better confidence in the resulting calculations. We also examine the structures of the fits to test whether, in principle, comparing the fit with ideal gray bodies would reveal such a structure and enable the researcher to identify when a two-temperature fitting model should be used instead of a single- T model or a wavelength-dependent emissivity model.

36 MATERIALS SCIENCE↗