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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Forward modeling approach to nuclear reaction cross sections: Applications in neutron inelastic scattering

The development of nuclear reaction models for the production of evaluated nuclear data has traditionally been performed by comparing measured cross sections with predictions from reaction model codes whose physical input parameters are adjusted to obtain the best agreement between measured and modeled results. To more directly probe reaction model inputs, this work introduces a forward modeling approach to experimental reaction cross-section determination, where the most important physical input parameters to reaction model calculations are obtained via 𝜒 2 minimization between measured and calculated observables. This was demonstrated using data collected by the Gamma Energy Neutron Energy Spectrometer for Inelastic Scattering (GENESIS) at the 88-inch cyclotron at Lawrence Berkeley National Laboratory, a detection array consisting of organic liquid scintillators and high-purity germanium (HPGe) detectors. Using a broad-spectrum neutron beam and a 99.98%-enriched 56 Fe target, GENESIS was used to perform a simultaneous measurement of 56 Fe 𝛾-ray production cross sections and secondary neutron energy and angle distributions. The results of the forward modeling approach to the determination of energy-differential 𝛾-ray production cross sections for the yrast 4 + → 2 + and 6 + → 4 + transitions, as well as eight other off-yrast transitions, were compared against those obtained using conventional techniques, and the results are in good agreement. In addition to discrete 𝛾-ray yield total scattered neutron energy-angular distributions as a function of incident neutron energy were also obtained using forward modeling and found to agree with evaluated data, with the exception of elastic scattering at small angles. The fitted reaction model parameters obtained through forward modeling were also used to calculate the cross section for the unobserved (𝑛, 2⁢𝑛) reaction; excellent agreement with the current evaluation was obtained, providing a validation of the predictive capabilities of the forward model approach. This work bridges the gap between nuclear data experiment and evaluation by providing a new means for extracting inelastic neutron-scattering cross sections and neutron-induced 𝛾-ray production data while directly probing reaction model physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2021

Model-based evaluation of school- and non-school-related measures to control the COVID-19 pandemic

The role of school-based contacts in the epidemiology of SARS-CoV-2 is incompletely understood. We use an age-structured transmission model fitted to age-specific seroprevalence and hospital admission data to assess the effects of school-based measures at different time points during the COVID-19 pandemic in the Netherlands. Our analyses suggest that the impact of measures reducing school-based contacts depends on the remaining opportunities to reduce non-school-based contacts. If opportunities to reduce the effective reproduction number (R e ) with non-school-based measures are exhausted or undesired and R e is still close to 1, the additional benefit of school-based measures may be considerable, particularly among older school children. As two examples, we demonstrate that keeping schools closed after the summer holidays in 2020, in the absence of other measures, would not have prevented the second pandemic wave in autumn 2020 but closing schools in November 2020 could have reduced R e below 1, with unchanged non-school-based contacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Adiabatic quantum linear regression

Abstract A major challenge in machine learning is the computational expense of training these models. Model training can be viewed as a form of optimization used to fit a machine learning model to a set of data, which can take up significant amount of time on classical computers. Adiabatic quantum computers have been shown to excel at solving optimization problems, and therefore, we believe, present a promising alternative to improve machine learning training times. In this paper, we present an adiabatic quantum computing approach for training a linear regression model. In order to do this, we formulate the regression problem as a quadratic unconstrained binary optimization (QUBO) problem. We analyze our quantum approach theoretically, test it on the D-Wave adiabatic quantum computer and compare its performance to a classical approach that uses the Scikit-learn library in Python. Our analysis shows that the quantum approach attains up to $${2.8 \times }$$ 2.8 × speedup over the classical approach on larger datasets, and performs at par with the classical approach on the regression error metric. The quantum approach used the D-Wave 2000Q adiabatic quantum computer, whereas the classical approach used a desktop workstation with an 8-core Intel i9 processor. As such, the results obtained in this work must be interpreted within the context of the specific hardware and software implementations of these machines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Fully generalized, turbulent trace impurity transport with Gkeyll and Flan in the DIII-D far-SOL

The Monte Carlo trace impurity turbulent transport code Flan is introduced for the first time. Flan follows impurities in a turbulent background plasma from Gkeyll using the Lorentz force to resolve the full particle gyro orbit. Collisions are handled using the Nanbu collision algorithm (Nanbu 1997 Phys. Rev. E 55 4642–52), and ionization/recombination is handled via ADAS coupling. The far-SOL of a generic DIII-D L-mode is simulated with and without the collision model to show how collisions affect radial tungsten transport. Anomalous diffusion coefficient (D r ) and pinch velocity (v p ) profiles are extracted from fits to the results. With collisions, D r and v p are between 0–1.0 m 2 s −1 and −100–100 m s −1 , respectively. Without collisions, D r and v p are between 0–0.3 m 2 s −1 and −50–50 m s −1 , respectively. Exponential fits to the radial W density profiles and experimental data from W deposition along a collector probe are in good agreement, demonstrating Flan as a useful interpretive modeling tool. Additional simulations show that impurity transport away from the wall increases with atomic number, though it is not clear why. Flan has the potential to better interpret existing data and improve reactor scale predictions of core contamination because the underlying physics model is very general and does not rely on arbitrary user-defined transport coefficients.

DIII-D↗
DOE OSTI2022

The MaNGA firefly Value-Added-Catalogue: resolved stellar populations of 10,010 nearby galaxies

We present the Mapping Nearby Galaxies at Apache Point Observatory survey (MaNGA) FIREFLY Value Added Catalogue (VAC) – a catalogue of ~3.7 million spatially resolved stellar population properties across 10 010 nearby galaxies from the final data release of the MaNGA survey. The full spectral fitting code FIREFLY is employed to derive parameters such as stellar ages, metallicities, stellar and remnant masses, star formation histories, star formation rates, and dust attenuation. In addition to Voronoi-binned measurements, our VAC also provides global properties, such as central values and radial gradients. Two variants of the VAC are available: presenting the results from fits using the M11-MILES and the novel MaStar stellar population models. MaStar allows to constrain the fit over the whole MaNGA wavelength range, extends the age-metallicity parameter space, and uses empirical spectra from the same instrument as MaNGA. The fits employing MaStar models find, on average, slightly younger ages, higher mass-weighted metallicities, and smaller colour excesses. These differences are reduced when matching the wavelength range and converging template grids. We further report that FIREFLY stellar masses are systematically lower by ~0.3 dex than masses from the MaNGA PCA and Pipe3D VACs, but match masses from the NSA best with only ~0.1-dex difference. Finally, we show that FIREFLY stellar ages correlate with spectral index age indicators H δ A and D n (4000), though with a clear additional metallicity dependence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Measurement of Reactor Antineutrino Oscillation at SNO+

SNO+ Collaboration reports its second spectral analysis of reactor antineutrino oscillation using 286 ton-yr of new data. The measured energies of reactor antineutrino candidates were fitted to obtain the second-most precise determination of the neutrino mass-squared difference Δ⁢𝑚$^{2}_{21}$=⁢(7.96$^{+0.48}_{ −0.42}$) × 10 −5 eV2. Constraining Δ⁢𝑚$^{2}_{21}$ and sin 2 ⁡𝜃 12 with measurements from long-baseline reactor antineutrino and solar neutrino experiments yields Δ⁢𝑚$^{2}_{21}$=(7.58$^{+0.18}_{−0.17}$) × 10 −5 eV 2 and sin 2 ⁡𝜃 12 = 0.308 ± 0.013. This fit also yields a first measurement of the flux of geoneutrinos in the Western Hemisphere, with 73$^{+47}_{−43}$ TNU at SNO+.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Anchor-based optimization of energy density functionals

A new anchor-based optimization method of defining energy density functionals (EDFs) is proposed. In this approach, the optimization of the parameters of EDFs is carried out for a selected set of spherical anchor nuclei, the physical observables of which are modified by the correction function, which takes into account the global performance of EDFs. It is shown that the use of this approach leads to a substantial improvement in the global description of binding energies for several classes of covariant EDFs. Here, the computational cost of defining a new functional within this approach is drastically lower as compared with the one for the optimization which includes the global experimental data on spherical, transitional, and deformed nuclei in the fitting protocol.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

SMEFT probes in future precision DIS experiments

We analyze the potential of future high-energy deep-inelastic scattering (DIS) experiments to probe new physics within the framework of the Standard Model effective field theory (SMEFT). We perform a detailed study of SMEFT probes at a future Large Hadron-electron Collider (LHeC) and a Future Circular lepton-hadron Collider (FCC-eh) machine, and extend previous simulations of the potential of an electron-ion collider (EIC) to include 𝑍-boson vertex corrections. Precision 𝑍-pole constraints on vertex corrections suffer from numerous degeneracies in the Wilson-coefficient parameter space. We find that both the LHeC and the FCC-eh can help remove these degeneracies present in the existing global fits of precision 𝑍-pole observables and LHC data. The FCC-eh and LHeC will in many cases improve upon the existing precision electroweak bounds on the SMEFT parameter space. This highlights the important role of precision DIS measurements for new physics studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Decays of N ∗ and ∆ ∗ resonances into Nρ , ∆ π , and N $f$ 0 (500)

The decays of N ∗ and ∆ ∗ resonances into Nρ , ∆ π , and N $f$ 0 (500), and final states are investigated in a coupled-channel analysis of data from pion- and photoinduced reactions. Improvements in the fit quality are observed upon the inclusion of additional resonance contributions. Branching ratios for the intermediate isobars ∆(1232) π , Nρ , and N $f$ 0 (500) are reported.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Hydrolytic Aging in Rubber-Like Materials: A Micro-Mechanical Approach to Modeling

Abstract The effect of hydrolytic aging on mechanical responses of Rubber likes materials, in particular, Mullins effect and the permanent set has been modeled. Hydrolytic aging is considered as the result of the competition between two phenomena (1) chain scission and (2) cross-link scission/reformation. Both phenomena were modeled and thus, the strain energy of the polymer matrix is written with respect to three independent mechanisms; i) the shrinking original matrix which has not been attacked by water, ii) conversion of the first network to a new network due to the reduction of the crosslinks, and iii) energy loss from network degradation due to water attacks to ester groups. The model is validated with respect to a set of experimental data. Besides accuracy, the simplicity and few numbers of fitting parameters make the model a good choice for further implementations.

Bahrololoumi, Amir↗
DOE OSTI2021

Understanding Hepatitis B Virus Dynamics and the Antiviral Effect of Interferon Alpha Treatment in Humanized Chimeric Mice

Whereas the mode of action of lamivudine (LAM) against hepatitis B virus (HBV) is well established, the inhibition mechanism(s) of interferon alpha (IFN-α) is less completely defined. To advance our understanding, we, for this study, mathematically modeled HBV kinetics during 14-day pegylated IFN-α-2a (pegIFN), LAM, or pegIFN-plus-LAM (pegIFN+LAM) treatment of 39 chronically HBV-infected humanized uPA/SCID chimeric mice. Serum HBV DNA and intracellular HBV DNA were measured frequently. We developed a multicompartmental mathematical model and simultaneously fit it to the serum and intracellular HBV DNA data. Unexpectedly, even in the absence of an adaptive immune response, a biphasic decline in serum HBV DNA and intracellular HBV DNA was observed in response to all treatments. Kinetic analysis and modeling indicate that the first phase represents inhibition of intracellular HBV DNA synthesis and secretion, which was similar under all treatments with an overall mean efficacy of 98%. In contrast, there were distinct differences in HBV decline during the second phase, which was accounted for in the model by a time-dependent inhibition of intracellular HBV DNA synthesis, with the steepest decline observed during pegIFN+LAM treatment (1.28/day) and the slowest (0.1/day) during pegIFN monotherapy. Reminiscent of observations in patients treated with pegIFN and/or LAM, a biphasic HBV decline was observed in treated humanized mice in the absence of an adaptive immune response. Interestingly, combination treatment did not increase the initial inhibition of HBV production but rather enhanced second-phase decline, providing insight into the dynamics of HBV treatment response and the mode of action of IFN-α against HBV.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Mass-Energy Compensation Effect of 3$\alpha$ Hamiltonian

The 3α phenomenological model describes the structure of the carbon-12 nucleus as a cluster of three alpha particles. This model includes a pairwise α–α interaction and a three-body force. To fit the three-body potential, the 12 C data are used, while ensuring that the pair potential reproduces the α–α scattering data. Alternatively, the mass-energy compensation (MEC) effect can be used to simulate the effect of the three-body potential by adjusting the mass of the α particle within the effective-mass approach. We demonstrate the MEC effect for the 3α ground state by numerically solving the differential Faddeev equation, in which the α–α interaction is described by the Ali-Bodmer potential. The effective masses of α particles are evaluated for the ground and excited 0 + and bound 2 + states. Here, we demonstrate a coupling between the ground and first excited 0 + states, indicated by an anti-crossing of these energy levels in the energy–mass coordinates. A correspondence between the effective mass and a three-body potential is demonstrated. We discuss the results of the 0$^{+}_{2}$ calculations for various models of the α–α interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Structure functions from renormalization group improved small x evolution

Abstract We perform the fit to the structure function $$F_2$$ F 2 data from HERA including terms due to the resummation at small x . The equation for the unintegrated gluon density is solved, previously established within the renormalization group improved small x framework. We find very good description of the structure function $$F_2$$ F 2 and its charm component $$F_2^c$$ F 2 c . The resulting unintegrated gluon density is found to be consistent with the calculations based on similar approaches available in the literature, with only slightly higher intercept value.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

listenr

SAND2024-13902O listenr is an R package containing code that allows users to fit echo state networks (a machine learning model) on data to obtain predictions. The code also allows users to compute spatio-temporal feature importance on the echo state network, as described in “Characterizing Climate Pathways Using Feature Importance on Echo State Networks.” The package also provides a function for computing principal components.

Ries, Daniel↗
DOE OSTI2025

jeanloz_method

A simple Jupyter notebook program for fitting isothermal compression curves to pressure-versus-volume data.

Brennan, Matthew↗
DOE OSTI2020

On the Solution of ℓ 0 -Constrained Sparse Inverse Covariance Estimation Problems

The sparse inverse covariance matrix is used to model conditional dependencies between variables in a graphical model to fit a multivariate Gaussian distribution. Estimating the matrix from data are well known to be computationally expensive for large-scale problems. Sparsity is employed to handle noise in the data and to promote interpretability of a learning model. Although the use of a convex ℓ 1 regularizer to encourage sparsity is common practice, the combinatorial ℓ 0 penalty often has more favorable statistical properties. In this paper, we directly constrain sparsity by specifying a maximally allowable number of nonzeros, in other words, by imposing an ℓ 0 constraint. Here, we introduce an efficient approximate Newton algorithm using warm starts for solving the nonconvex ℓ 0 -constrained inverse covariance learning problem. Numerical experiments on standard data sets show that the performance of the proposed algorithm is competitive with state-of-the-art methods.

$\ell_0$-Constrained↗
DOE OSTI2022

Community Based Data of Uranium Adsorption onto Quartz

This data upload includes a compilation of experiments for quantifying uranium adsorption onto quartz. The provided .csv file has been compiled from the literature in a findable, accessible, interoperable, reusable (FAIR) data format. This was accomplished using the Lawrence Livermore National Laboratory Surface Complexation Database Converter (SCDC) code written in the R programming language (free licensing available at https://ipo.llnl.gov/technologies/software/llnl-surface-complexation-database-converter-scdc). This constitutes all current data compiled on uranium-quartz interactions in the L-SCIE (LLNL Surface Complexation/Ion Exchange) database (as of 02/28/2022). This data was used to develop surface complexation models that fit the global community dataset (https://doi.org/10.1021/acs.est.1c07109). The FAIR-formatted dataset also enables the implementation of alternative machine-learning approaches that can be explored in the future.

54 ENVIRONMENTAL SCIENCES↗