Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed inference”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Inferred inflow forecast horizons guiding reservoir release decisions across the United States

Abstract. Medium- to long-range forecasts often guide reservoir release decisions to support water management objectives, including mitigating flood and drought risks. While there is a burgeoning field of science targeted at improving forecast products and associated decision support models, data describing how and when forecasts are applied in practice remain undeveloped. This lack of knowledge may prevent hydrological modelers from developing accurate reservoir release schemes for large-scale, distributed hydrology models that are increasingly used to assess the vulnerabilities of large regions to hydrological stress. We address this issue by estimating seasonally varying, regulated inflow forecast horizons used in the operations of more than 300 dams throughout the conterminous United States (CONUS). For each dam, we take actual forward observed inflows (perfect foresight) as a proxy for forecasted flows available to the operator and then identify for each week of the year the forward horizon that best explains the release decisions taken. Resulting “horizon curves” specify for each dam the inferred inflow forecast horizon as a function of the week of the water year. These curves are analyzed for strength of evidence for contribution of medium- to long-range forecasts in decision making. We use random forest classification to estimate that approximately 80 % of large dams and reservoirs in the US (1553±50 out of 1927 dams with at least 10 Mm3 storage capacity) adopt medium- to long-range inflow forecasts to inform release decisions during at least part of the water year. Long-range forecast horizons (more than 6 weeks ahead) are detected in the operations of reservoirs located in high-elevation regions of the western US, where snowpack information likely guides the release. A simulation exercise conducted on four key western US reservoirs indicates that forecast-informed models of reservoir operations may outperform models that neglect the horizon curve – including during flood and drought conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Single-Molecule Tracking Measurements Reveal the Detailed Mechanisms of Molecular Diffusion in Solvent Mixtures under Nanoconfinement

Understanding mass transport mechanisms in nanopores is important for developing advanced materials for chemical separations, chemical sensing, and energy storage. This paper reports a novel imaging platform that is employed for the first time to investigate the detailed diffusion dynamics of single rhodamine B (RhB) dye molecules confined within solution-filled cylindrical anodic aluminum oxide (AAO) nanopores. The imaging platform relies on illumination of horizontally-oriented AAO nanopores in a highly inclined and laminated optical (HILO) light sheet microscopy geometry. The method was used to investigate the translational and orientational dynamics of single rhodamine B (RhB) molecules within horizontally-oriented 5- and 10-nm diameter AAO nanopores filled with water–ethanol mixtures. The established platform enabled the observation of one-dimensional motion along the pore axis involving occasional short- or long-term immobilization at the single-molecule level. Analysis of cumulative squared-displacement distributions revealed fast (5 – 30 µm²/s), intermediate (1 – 5 µm²/s), and slow (< 1 µm²/s) diffusion components. From the effects of mixture composition and pore size on the contributions of these three components, we inferred that the fast, intermediate, and slow components could be assigned to desorption-mediated hopping, crawling, and wiggling motions, respectively. The platform based on the horizontally-oriented AAO nanopores also permitted single-molecule emission polarization measurements that revealed the negligible steric confinement of individual diffusing RhB molecules within the AAO nanopores. The imaging platform based on AAO membranes and HILO microscopy provided a unique means to investigate how solvation-mediated surface interactions and nanoconfinement govern molecular transport in nanoporous environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Flat-spectrum Radio Quasars and BL Lacs Dominate the Anisotropy of the Unresolved Gamma-Ray Background

We analyze the angular power spectrum (APS) of the unresolved gamma-ray background (UGRB) emission and combine it with the measured properties of the resolved gamma-ray sources of the Fermi-LAT 4FGL catalog. Our goals are to dissect the composition of the gamma-ray sky and to establish the relevance of different classes of source populations of active galactic nuclei in determining the observed size of the UGRB anisotropy, especially at low energies. We find that, under physical assumptions for the spectral energy distribution, i.e., by using the 4FGL catalog data as a prior, two populations are required to fit the APS data, namely flat-spectrum radio quasars at low energies and BL Lacs at higher energies. The inferred luminosity functions agree well with the extrapolation of the flat-spectrum radio quasar and BL Lac ones obtained from the 4FLG catalog. We use these luminosity functions to calculate the UGRB intensity from blazars, finding a contribution of 20% at 1 GeV and 30% above 10 GeV. Finally, bounds on an additional gamma-ray emission due to annihilating dark matter are also derived.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc − 3 yr − 1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc − 3 yr − 1 , assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc − 3 yr − 1 at a fiducial redshift ( z = 0.2 ). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to ( 1 + z ) κ with κ = 2. 9 − 1.8 + 1.7 for z ≲ 1 . Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2 − 0.2 + 0.1 to 2.0 − 0.3 + 0.3 M ⊙ . We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3 − 0.5 + 0.3 and 27.9 − 1.8 + 1.9 M ⊙ . While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60 M ⊙ , which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χ i ≈ 0.25 . While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Beyond Point Estimates: Benchmarking Uncertainty Quantification Methods on the AION-1 Astronomical Foundation Model

Foundation models for astronomical surveys offer powerful learned representations that can be transferred to downstream regression tasks such as galaxy property estimation. However, point predictions alone are insufficient for scientific inference; reliable uncertainty quantification (UQ) is essential. We compare seven UQ methods on galaxy property regression using frozen AION-1 foundation-model embeddings, predicting redshift, stellar mass, stellar-population age, gas-phase metallicity, and specific star-formation rate, from Legacy Survey photometry/imaging and DESI spectra, with PROVABGS-derived labels. Distribution-free conformal methods achieve marginal coverage within $\sim$1 pp of the nominal 90% across all properties, while non-conformal baselines (Deep Ensembles, MC~Dropout) fail to calibrate reliably. Among conformal approaches, Conformalized Quantile Regression (CQR) delivers the best coverage in the bin with the poorest model predictions. More importantly, only the Locally Valid and Discriminative (LVD) framework -- particularly when operating on AION-1 embeddings -- also provides finite-sample \emph{local validity}, producing intervals that adapt to each galaxy's local prediction difficulty rather than relying on marginal guarantees alone. These results establish conformal prediction, and LVD in particular, as the preferred UQ framework for uncertainty-aware inference on foundation-model embeddings in astrophysics.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗
DOE OSTI2024

Applications of emulation and Bayesian methods in heavy-ion physics

Abstract Heavy-ion collisions provide a window into the properties of many-body systems of deconfined quarks and gluons. Understanding the collective properties of quarks and gluons is possible by comparing models of heavy-ion collisions to measurements of the distribution of particles produced at the end of the collisions. These model-to-data comparisons are extremely challenging, however, because of the complexity of the models, the large amount of experimental data, and their uncertainties. Bayesian inference provides a rigorous statistical framework to constrain the properties of nuclear matter by systematically comparing models and measurements. This review covers model emulation and Bayesian methods as applied to model-to-data comparisons in heavy-ion collisions. Replacing the model outputs (observables) with Gaussian process emulators is key to the Bayesian approach currently used in the field, and both current uses of emulators and related recent developments are reviewed. The general principles of Bayesian inference are then discussed along with other Bayesian methods, followed by a systematic comparison of seven recent Bayesian analyses that studied quark-gluon plasma properties, such as the shear and bulk viscosities. The latter comparison is used to illustrate sources of differences in analyses, and what it can teach us for future studies.

Paquet, Jean-François (ORCID:0000000187368171)↗
DOE OSTI2021

The differences between mass- and light-derived structural parameters over time for MaNGA elliptical galaxies

ABSTRACT We apply stellar population synthesis analysis to obtain spatially resolved archaeological inferences for a large sample of ‘red and dead’ Elliptical galaxies (Classical Ellipticals, CLEs) from the MaNGA/SDSS-IV DR15 survey. From their 2D stellar light and mass maps, we explore the differences between the radial mass and light distributions in the rest-frame bands g, r, and i as functions of look-back time, tlb, or redshift, z. We characterize these differences through the ratios between the following mass- and light-derived global properties: sizes, concentrations, and effective surface densities. We find that the mass-to-light ratios of these properties change with tlb, more the more massive the galaxies are. The CLE galaxy archaeological progenitors are, on average, less compact, concentrated, and dense in light than in mass as z decreases. However, at later times, when also the evolution of the progenitors becomes passive at all radii, there is an upturn in these trends and the differences between mass and light in compactness/concentration decrease towards z ∼ 0. The trends in the ratios of mass-to-light sizes agree qualitatively with results from direct observations in galaxy surveys at different redshifts. We discuss the caveats and interpretations of our results, and speculate that the strong structural evolution found in some previous studies for early-type galaxies could be explained partially by photometric changes rather than by intrinsic structural changes.

Ibarra-Medel, H. (ORCID:0000000297906313)↗
DOE OSTI2026

Modeling the Cosmological Lyman-𝛼 Forest at the Field Level

The distribution of absorption lines in the spectra of distant quasars, called the Lyman-𝛼 (Ly-𝛼) forest, is a unique probe of cosmology and the intergalactic medium at high redshifts and small scales. The statistical power of ongoing redshift surveys demands precise theoretical tools to model the Ly-𝛼 forest. We address this challenge by developing an analytic, perturbative forward model to predict the Ly-𝛼 forest at the field level for a given set of cosmological initial conditions. Our model shows a remarkable performance when compared with the Sherwood hydrodynamic simulations: it reproduces the Ly-𝛼 forest flux power spectrum, its cross-correlation with dark matter halos, and the one-point probability distribution function of both fields at the percent level down to scales of a few Mpc. Our work provides crucial tools that bridge analytic modeling on large scales with simulations on small scales, enabling field-level inference from Ly-𝛼 forest data and simulation-based priors for cosmological analyses. Furthermore, this is especially timely for realizing the full scientific potential of the Ly-𝛼 forest measurements by the dark energy spectroscopic instrument.

Cosmological parameters↗
DOE OSTI2021

Evaluation of liquid-phase thermometry in impinging jet sprays using synchrotron x-ray scattering

Liquid thermometry during primary and secondary breakup of liquid sprays is challenging due to the presence of highly dynamic, optically complex flow features. This work evaluates the use of x-ray scattering from a focused, monochromatic beam of the Advanced Photon Source at Argonne National Laboratory for the measurement of liquid temperatures within the mixing zone of an impinging jet spray. The measured scattering profiles are converted to temperature through a previously developed two-component partial least squares (PLS) regression model. Transmitive mixing during jet merging is inferred through spatial mapping of temperatures within the impingement region. In this work, the technique exhibits uncertainties of ±2K in temperature and 2% in capturing the correct scattering profile, showing its potential utility for probing liquid temperature distributions in multiphase flows.

42 ENGINEERING↗
DOE OSTI2022

Equation of State Optimization and Uncertainty Quantization: Implementation in the LANL EOS Production Code OpenSesame

We detail the approaches of particle swarm optimization and Bayesian inference through Markov chain Monte Carlo for equation of state development. This work includes formulation of the modeling for the equation of state, numeric optimization of the parametric models via particle swarm optimization, and generation of probability distributions of equations of state from Markov chain Monte Carlo.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗
DOE OSTI2022

Strong Lensing Parameter Estimation on Ground-Based Imaging Data Using Simulation-Based Inference

Current ground-based cosmological surveys, such as the Dark Energy Survey (DES), are predicted to discover thousands of galaxy-scale strong lenses, while future surveys, such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will increase that number by 1-2 orders of magnitude. The large number of strong lenses discoverable in future surveys will make strong lensing a highly competitive and complementary cosmic probe. To leverage the increased statistical power of the lenses that will be discovered through upcoming surveys, automated lens analysis techniques are necessary. We present two Simulation-Based Inference (SBI) approaches for lens parameter estimation of galaxy-galaxy lenses. We demonstrate the successful application of Neural Posterior Estimation (NPE) to automate the inference of a 12-parameter lens mass model for DES-like ground-based imaging data. We compare our NPE constraints to a Bayesian Neural Network (BNN) and find that it outperforms the BNN, producing posterior distributions that are for the most part both more accurate and more precise; in particular, several source-light model parameters are systematically biased in the BNN implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Diagnosing plasma magnetization in inertial confinement fusion implosions using secondary deuterium-tritium reactions

Diagnosing plasma magnetization in inertial confinement fusion implosions is important for understanding how magnetic fields affect implosion dynamics and to assess plasma conditions in magnetized implosion experiments. Secondary deuterium–tritium (DT) reactions provide two diagnostic signatures to infer neutron-averaged magnetization. Magnetically confining fusion tritons from deuterium–deuterium (DD) reactions in the hot spot increases their path lengths and energy loss, leading to an increase in the secondary DT reaction yield. In addition, the distribution of magnetically confined DD-triton is anisotropic, and this drives anisotropy in the secondary DT neutron spectra along different lines of sight. Implosion parameter space as well as sensitivity to the applied B-field, fuel ρR, temperature, and hot-spot shape will be examined using Monte Carlo and 2D radiation-magnetohydrodynamic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗