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 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2022

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↗

DOE OSTI2020

Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations

Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.

58 GEOSCIENCES↗

DOE OSTI2021

Hubble WFC3 Spectroscopy of the Habitable-zone Super-Earth LHS 1140 b

Atmospheric characterization of temperate, rocky planets is the holy grail of exoplanet studies. These worlds are at the limits of our capabilities with current instrumentation in transmission spectroscopy and challenge our state-of-the-art statistical techniques. Here we present the transmission spectrum of the temperate super-Earth LHS 1140b using the Hubble Space Telescope (HST). The Wide Field Camera 3 (WFC3) G141 grism data of this habitable-zone (T {sub eq} = 235 K) super-Earth (R = 1.7 R {sub ⊕}) shows tentative evidence of water. However, the signal-to-noise ratio, and thus the significance of the detection, is low and stellar contamination models can cause modulation over the spectral band probed. We attempt to correct for contamination using these models and find that, while many still lead to evidence for water, some could provide reasonable fits to the data without the need for molecular absorption although most of these cause features in the visible ground-based data which are nonphysical. Future observations with the James Webb Space Telescope would be capable of confirming, or refuting, this atmospheric detection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

A comparison of explosively driven shock wave radius versus time scaling approaches

Abstract Explosively driven shock wave radius versus time profiles are frequently used to document energy release and relative explosive performance. Recently, two universal shock wave radius versus time profiles have been presented in the literature, which demonstrate the ability to represent explosively driven shock wave profiles for all explosive sources in any fluid environment. These two universal shock wave profiles are examined here relative to each other and relative to a commonly used nonlinear shock wave profile, which is fit to experimental data for individual explosive materials. The nonlinear profile, originally developed by Dewey, is examined here, and a universal non-dimensional form of the equation is proposed. The universal shock wave profiles are all found to be relatively similar, but with slight variations in a transition region of non-dimensional radii $$0.15\lesssim R^*\lesssim 2$$ 0.15 ≲ R ∗ ≲ 2 . The variations in this region result in different estimations of energy release or blast strength between the curve fits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

DOE OSTI2026

The S-PLUS Fifth Data Release: Over 4500 Square Degrees of the Southern Sky and a Multicolor View of the Hydra and Antlia Galaxy Clusters

Abstract We present the fifth data release (DR5) of the Southern Photometric Local Universe Survey (S-PLUS), covering 4592 deg 2 across 2491 fields. Observations were conducted with the T80-South, a Brazilian robotic telescope equipped with the Javalambre 12-filter system, containing five broadband and seven narrowband filters. Data products feature FITS images and extensive catalogs containing fluxes, magnitudes, and shape parameters for over 113 million detections. In addition, several value-added catalogs are provided, offering photometric redshifts (photo- zs ), object classifications, masks, and extinction coefficients. For the first time, this release includes coverage of 110 deg 2 along the Galactic disk, facilitating new research into Galactic structure and stellar populations. The release also provides full coverage of the Hydra Supercluster and numerous other nearby clusters with improved data reduction and calibration, enhancing photo- z accuracy, which is vital for large-scale structure studies. A preliminary analysis of the Hydra and Antlia galaxy clusters up to 5 × R 200 yields an updated catalog of 1706 cluster members based on both spectroscopic and high-quality photo- zs . Our photometric data shows that both clusters have a similar proportion of galaxies with an H α excess relative to their clustercentric distance, though Hydra has a higher fraction near its center. Additionally, the spatial distribution of all objects in our sample highlights a bridge connecting both clusters. We verify that S-PLUS DR5 provides a solid foundation for future scientific investigations, ranging from solar system studies to cosmology.

Vinicius Rodrigues de Lima, Erik [Universidade de ↗

DOE OSTI2024

Polyconvex neural network models of thermoelasticity

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. Here, in this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Concordance cosmology?

ABSTRACT We propose a new intuitive metric for evaluating the tension between two experiments, and apply it to several data sets. While our metric is non-optimal, if evidence of tension is detected, this evidence is robust and easy to interpret. Assuming a flat Lambda cold dark matter (ΛCDM) cosmological model, we find that there is a modest 2.2σ tension between the Dark Energy Survey (DES) Year 1 results and the Planck measurements of the cosmic microwave background. This tension is driven by the difference between the amount of structure observed in the late-time Universe and that predicted from fitting the Planck data, and appears to be unrelated to the tension between Planck and local estimates of the Hubble rate. In particular, combining DES, baryon acoustic oscillations, big bang nucleosynthesis, and supernovae measurements recover a Hubble constant and sound horizon consistent with Planck, and in tension with local distance–ladder measurements. If the tension between these various data sets persists, it is likely that reconciling all current data will require breaking the flat ΛCDM model in at least two different ways: one involving new physics in the early Universe, and one involving new late-time Universe physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Full waveform analysis for high pressure ultrasonic measurement

In this paper we report a data analysis protocol for ultrasonic velocity measurements carried out in a synchrotron x-ray facility based multi-anvil high-pressure apparatus. Synthetic ultrasonic signals for the time period between echoes from the two ends of the sample assembly are created with a few simple parameters as the echo waveform is modeled by the first echo from the anvil-assembly boundary. Each echo is modeled as a delta function with an amplitude and arrival time and convolved with this waveform. The final waveform, fit to the data, enables analysis of signals with overlapped echoes as they are common in small cell assemblies and high pressure.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Interacting radiation after Planck and its implications for the Hubble tension

Standard cosmology predicts that prior to matter-radiation equality about 41% of the energy density was in free-streaming neutrinos. In many beyond Standard Model scenarios, however, the amount and free-streaming nature of this component is modified. For example, this occurs in models with new neutrino self-interactions or an additional dark sector with interacting light particles. Here, we consider several extensions of the standard cosmology that include a non-free-streaming radiation component as motivated by such particle physics models and use the final Planck data release to constrain them. This release contains significant improvements in the polarization likelihood which plays an important role in distinguishing free-streaming from interacting radiation species. Fixing the total amount of energy in radiation to match the expectation from standard neutrino decoupling we find that the fraction of free-streaming radiation must be ffs > 0.8 at 95% CL (combining temperature, polarization and baryon acoustic oscillation data). Allowing for arbitrary contributions of free-streaming and interacting radiation, the effective number of new non-free-streaming degrees of freedom is constrained to be Nfld < 0.6 at 95% CL. Cosmologies with additional radiation are also known to ease the discrepancy between the local measurement and CMB inference of the current expansion rate H0. We show that including a non-free-streaming radiation component allows for a larger amount of total energy density in radiation, leading to a mild improvement of the fit to cosmological data compared to previously discussed models with only a free-streaming component.

particle physics - cosmology connection↗

DOE OSTI2021

Single $\pi ^0$ production off neutrons bound in deuteron with linearly polarized photons

The quasifree γ → d → π 0 n(p) photon beam asymmetry, Σ, has been measured at photon energies, E γ , from 390 to 610 MeV, corresponding to center of mass energy from 1.271 to 1.424 GeV, for the first time. The data were collected in the A2 hall of the MAMI electron beam facility with the Crystal Ball and TAPS calorimeters covering pion center-of-mass angles from 49° to 148°. In this kinematic region, polarization observables are sensitive to contributions from the Δ(1232) and N(1440) resonances. The extracted values of Σ have been compared to predictions based on partial-wave analyses (PWAs) of the existing pion photoproduction database. Our comparison includes the SAID, MAID and Bonn–Gatchina analyses; while a revised SAID fit, including the new Σ measurements, has also been performed. In addition, isospin symmetry is examined as a way to predict π 0 n photoproduction observables, based on fits to published data in the channels π 0 p, π + n and π – p.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗

DOE OSTI2025

AQuaRef: machine learning accelerated quantum refinement of protein structures

Cryo-EM and X-ray crystallography provide crucial experimental data for obtaining atomic-detail models of biomacromolecules. Refining these models relies on library-based stereochemical data, which, in addition to being limited to known chemical entities, do not include meaningful noncovalent interactions. Quantum mechanical (QM) calculations could alleviate these issues but are too expensive for large molecules. Here we present a novel AI-enabled Quantum Refinement (AQuaRef) based on AIMNet2 machine learned interatomic potential (MLIP) mimicking QM at substantially lower computational costs. By refining 41 cryo-EM and 30 X-ray structures, we show that this approach yields atomic models with superior geometric quality compared to standard techniques, while maintaining an equal or better fit to experimental data. Notably, AQuaRef aids in determining proton positions, as illustrated in the challenging case of short hydrogen bonds in the parkinsonism-associated human protein DJ-1 and its bacterial homolog YajL.

Zubatyuk, Roman [Carnegie Mellon University, Pitts↗

DOE OSTI2024

Effects of renormalon scheme and perturbative scale choices on determinations of the strong coupling from e + e − event shapes

We study the role of renormalon cancellation schemes and perturbative scale choices in extractions of the strong coupling constant α s ( m Z ) and the leading nonperturbative shift parameter Ω 1 from resummed predictions of the e + e − event shape thrust. We calculate the thrust distribution to N L 3 L ′ resummed accuracy in soft-collinear effective theory (SCET) matched to the fixed-order O ( α s 2 ) prediction, and perform a new high-statistics computation of the O ( α s 3 ) matching in , although we do not include the latter in our final α s fits due to some observed systematics that require further investigation. We are primarily interested in testing the phenomenological impact sourced from varying amongst three renormalon cancellation schemes and two sets of perturbative scale profile choices. We then perform a global fit to available data spanning center-of-mass energies between 35–207 GeV in each scenario. Relevant subsets of our results are consistent with prior SCET-based extractions of α s ( m Z ) , but we are also led to a number of novel observations. Notably, we find that the combined effect of altering the renormalon cancellation scheme and profile parameters can lead to few-percent-level impacts on the extracted values in the α s − Ω 1 plane, indicating a potentially important systematic theory uncertainty that should be accounted for. We also observe that fits performed over windows dominated by dijet events are typically of a higher quality than those that extend into the far tails of the distributions, possibly motivating future fits focused more heavily in this region. Finally, we discuss how different estimates of the three-loop soft matching coefficient c S ˜ 3 can also lead to measurable changes in the fitted { α s , Ω 1 } values. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗

DOE OSTI2020

EXAFS investigation of the local structure in URu 2-x Fe x Si 2 : Evidence for distortions below 100 K

X-ray absorption measurements at the U L III , Ru K , and Fe K edges are reported for the hidden order (HO) material URu 2 - x Fe x Si 2 ( x = 0 , 0.05, 0.08, 0.10, 0.12, 0.15, and 0.20) as a function of x and temperature T . Furthermore, when Fe is substituted for Ru, the local structure about Fe shrinks slightly and the first neighbor Fe-Si bond length decreases by ≈ 0.05 Å . More importantly excess disorder is observed below 80–100 K (the coherence temperature T * ) in plots of the Debye-Waller factor σ 2 ( σ is the width of the pair distribution function); at low T the data deviate from the usual Einstein or correlated-Debye model plots. This excess disorder is most prominent for the Ru-Si bond, and σ 2 actually increases below 80 K. These results suggest a local orthorhombic distortion with B 1 g -like symmetry that develops below 80–100 K. A model that describes these local distortions is presented, and discussed in terms of other measurements that indicate a breaking of fourfold symmetry at low T . In addition, the square root of the difference between σ 2 ( T ) for the Ru-Si pair and a Debye fit to these data serves as an order parameter for this orthorhombic distortion, in the temperature range below 100 K. This quantity is a length related to a - b , the difference between the a and b lattice constants in the orthorhombic phase, and provides a connection between this distortion and T * . X-ray absorption near edge structure (XANES) measurements also show that there are no changes in the edge positions down to 0.1 eV for any edge as a function of x , for T in the HO regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 304L Stainless Steel Alloy

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Predictions for the Detectability of Milky Way Satellite Galaxies and Outer-Halo Star Clusters with the Vera C. Rubin Observatory

We predict the sensitivity of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) to faint, resolved Milky Way satellite galaxies and outer-halo star clusters. We characterize the expected sensitivity using simulated LSST data from the LSST Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) accessed and analyzed with the Rubin Science Platform as part of the Rubin Early Science Program. We simulate resolved stellar populations of Milky Way satellite galaxies and outer-halo star clusters over a wide range of sizes, luminosities, and heliocentric distances, which are broadly consistent with expectations for the Milky Way satellite system. We inject simulated stars into the DC2 catalog with realistic photometric uncertainties and star/galaxy separation derived from the DC2 data itself. We assess the probability that each simulated system would be detected by LSST using a conventional isochrone matched-filter technique. We find that assuming perfect star/galaxy separation enables the detection of resolved stellar systems with $M_V$ = 0 mag and $r_{1/2}$ = 10 pc with >50% efficiency out to a heliocentric distance of ~250 kpc. Similar detection efficiency is possible with a simple star/galaxy separation criterion based on measured quantities, although the false positive rate is higher due to leakage of background galaxies into the stellar sample. When assuming perfect star/galaxy classification and a model for the galaxy-halo connection fit to current data, we predict that 89 +/- 20 Milky Way satellite galaxies will be detectable with a simple matched-filter algorithm applied to the LSST wide-fast-deep data set. Different assumptions about the performance of star/galaxy classification efficiency can decrease this estimate by ~7%-25%, which emphasizes the importance of high-quality star/galaxy separation for studies of the Milky Way satellite population with LSST.

79 ASTRONOMY AND ASTROPHYSICS↗