SEARCH · Engineering Papers

Results for “Distribution parameters”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Electric fields and stationary drift flows in the island divertor SOL of Wendelstein 7-X

Abstract In the island divertor scrape-off layer (SOL) of the Wendelstein 7-X stellarator, 2D-resolved Langmuir probe measurements uncover a complex spatial structure of electron temperature, density, and floating potential. Depending on details of magnetic island size and position, significant poloidal gradients and non-monotonic radial profiles are observed. The resulting electric fields imply the presence of sheared stationary drift flows with poloidal velocities of a few km/s, with flow channels as narrow as one cm. Such flows are observed with a gas puff imaging diagnostic in qualitative agreement with the probe results and imply a significant transport channel in the SOL. The 2D distribution of plasma parameters is intricately associated to the unique magnetic structure of the island divertor SOL and highlights the necessity for more accurate magnetic field reconstructions.

Killer, Carsten (ORCID:0000000177473066)↗

DOE OSTI2021

Landau Damping with Electron Lenses in Space-Charge Dominated Beams

Further progress of fundamental physics requires accelerated beams of high intensity. The intensities, however, are limited by many factors, including coherent beam instabilities. Usual methods to control the instabilities, such as octupole magnets, beam feedback dampers and employment of chromatic effects, may be ineffective or insufficient. In this study, electron lenses were proposed as a means to provide stabilizing spread in the beam betatron frequencies. It was shown that electron lenses are uniquely effective for Landau damping of transverse beam instabilities in high energy particle accelerators, and that their employment does not compromise incoherent (single particle) stability, dynamic aperture and the beam lifetime. Here we consider effectiveness of the Landau damping with electron lenses when the space charge tune shift cannot be neglected. We demonstrate that the desired stability can be assured with proper choice of the electron beam parameters and current distributions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Three-body capture, ejection, and the demographics of bound objects in binary systems

ABSTRACT We study the capture of light objects of arbitrary velocity by binary systems. Extending results for the capture of comets in the Solar system, we develop a simple geometric characterization of the capture cross-section, leading directly to the distribution of orbital parameters of captured objects. We use the same framework to study the lifetimes of these bound orbits prior to ejection, and find that a simplified version of the Öpik–Arnold approach readily yields a closed-form estimate for the ejection rate that agrees well with numerical experiments. Without any detailed-balance assumptions, our results make manifest the characteristics of close encounters leading to capture and ejection. As an application of our results, we demonstrate the estimation of the equilibrium population of captured dark matter particles in a binary system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗

DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Observation of time-dependent C P violation and measurement of the branching fraction of B 0 → J / ψ π 0 decays

We present a measurement of the branching fraction and time-dependent charge-parity ( C P ) decay-rate asymmetries in B 0 → J / ψ π 0 decays. The data sample was collected with the Belle II detector at the SuperKEKB asymmetric e + e − collider in 2019–2022 and contains ( 387 ± 6 ) × 10 6 B B ¯ meson pairs from ϒ ( 4 S ) decays. We reconstruct 392 ± 24 signal decays and fit the C P parameters from the distribution of the proper-decay-time difference of the two B mesons. We measure the branching fraction to be ( B 0 → J / ψ π 0 ) = ( 2.00 ± 0.12 ± 0.09 ) × 10 − 5 and the direct and mixing-induced C P asymmetries to be C C P = 0.13 ± 0.12 ± 0.03 and S C P = − 0.88 ± 0.17 ± 0.03 , respectively, where the first uncertainties are statistical and the second are systematic. We observe mixing-induced C P violation with a significance of 5.0 standard deviations for the first time in this mode. Published by the American Physical Society 2025

Adachi, I. (ORCID:0000000322870173)↗

DOE OSTI2025

Data Assimilation for Robust UQ Within Agent-Based Simulation on HPC Systems

Agent-based simulation provides a powerful tool for in silico system modeling. However, these simulations do not provide built-in methods for uncertainty quantification (UQ). Within these types of models a typical approach to UQ is to run multiple realizations of the model then compute aggregate statistics. This approach is limited due to the compute time required for a solution. When faced with an emerging biothreat, public health decisions need to be made quickly and solutions for integrating near real-time data with analytic tools are needed. We propose an integrated Bayesian UQ framework for agent-based models based on sequential Monte Carlo sampling. Given streaming or static data about the evolution of an emerging pathogen this Bayesian framework provides a distribution over the parameters governing the spread of a disease through a population. These estimates of the spread of a disease may be provided to public health agencies seeking to abate the spread. By coupling agent-based simulations with Bayesian modeling in a data assimilation, our proposed framework provides a powerful tool for modeling dynamical systems in silico. We propose a method which reduces model error and provides a range of realistic possible outcomes. Moreover, our method addresses two primary limitations of ABMs: the lack of UQ and an inability to assimilate data. Our proposed framework combines the flexibility of an agent-based model with UQ provided by the Bayesian paradigm in a workflow which scales well to HPC systems. We provide algorithmic details and results on a simulated outbreak with both static and streaming data.

Spannaus, Adam [ORNL] (ORCID:0000000225213657)↗

DOE OSTI2025

time-resolved spectroscopy fit (trspecfit) v0.01

Analyze 2D time- and energy-resolved data, such as from a pump-probe spectroscopy experiment. User can select and input different peak shapes/ functions and background types to first fit a ground state/ unperturbed spectrum. This would be similar to how standard spectroscopy data is fit. Subsequently, to describe the time domain, users can choose functions that describe the temporal dynamics of one or more spectral features, such as a peak amplitude, peak position, etc. These time dynamics functions can be added or convoluted (e.g. describing an instrument response function) with each other. Functionality to integrate implicit variables leading to distributions of certain parameters/ functions is in development. Alternatively, 2D data can be analyzed one time step at a time to get an idea of the time dynamics of the system before deploying the global 2D fit described above. Typically people write custom software for this purpose. During my PhD I've seen five internal LBL and external researchers write one-off code in different languages to analyze time- and energy-resolved spectra. While this was specifically was for a laser pump - X-ray probe spectroscopy experiment, I'm trying to write a general package for the time-resolved spectroscopy community.

Mahl, Johannes [Lawrence Berkeley National Laborat↗

DOE OSTI2022

Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code

The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

A Generative Model for Synthetic Electroluminescence Images

This work will develop modular, open-source model and analysis components including crack detection workflow and parameterization for quantitative inspection of large EL large datasets. These tools will allow users to quickly and accurately assess the extent and types of cracking in their modules. Measured statistical distributions of crack parameters, together with the imposed stress and electrical properties will be used to generate models to predict future crack behavior and power loss.

Pierce, Benjamin Garrett↗

DOE OSTI2026

Probing Non-Standard Interactions in NOvA with Bayesian MCMC Methods

We present a search for non-standard neutrino interactions (NSI) using NOvA's joint $\nu_e$ appearance and $\nu_\mu$ disappearance samples in both neutrino and antineutrino beam modes. A Bayesian Markov Chain Monte Carlo approach is used to map the posterior distribution over NSI parameters $\varepsilon_{e\mu}$, $\varepsilon_{e\tau}$, and $\varepsilon_{\mu\tau}$, simultaneously with standard oscillation parameters. We investigate the effect of prior choice for the complex NSI parameters and present projected sensitivity to off-diagonal NSI.

Huang, Xiaoyan [Mississippi U.]↗

DOE OSTI2023

Signals of a New Gauge Boson from IceCube and Muon g-2

Gauging the muon lepton flavor minus the tau lepton flavor number, which is a global symmetry of the SM, introduces a new gauge boson Z' upon symmetry breaking. Interestingly, it offers an economical solution to the long-standing g_µ-2 anomaly, confirmed and strengthened by recent measurements at Fermilab. Here, we revisit the impact of such a Z' on the spectrum of high-energy astrophysical neutrinos, as measured by the IceCube experiment. This spectrum has been observed to exhibit a dip-like feature at sub-PeV energies, which could plausibly arise from the physics of the sources themselves, but could also be the consequence of high-energy neutrinos resonantly scattering with the cosmic neutrino background, mediated by a Z' with a mass on the order of m_Z’ ~ 10 MeV. In this study, we calculate the impact of such a Z' on the high-energy neutrino spectrum for a variety of model parameters and source distributions. For couplings that can resolve the g-2 anomaly, we find that this model could self-consistently produce a spectral feature that is consistent with IceCube's measurement, in particular if the neutrinos observed by IceCube predominantly originate from high-redshift sources. We also briefly discuss a possible scenario where Z' could act as a portal to a dark sector.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Incremental Fermi Large Area Telescope Fourth Source Catalog

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi Large Area Telescope (LAT) catalog of γ-ray sources. Based on the first 12 years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss the degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Signals of a New Gauge Boson from IceCube and Muon $g-2$

A $Z'$ boson associated with a broken $U(1)_{L_{\mu} - L_{\tau}}$ gauge symmetry offers an economical solution to the long-standing $g_\mu-2$ anomaly, confirmed and strengthened by recent measurements at Fermilab. Here, we revisit the impact of such a $Z'$ on the spectrum of high-energy astrophysical neutrinos, as measured by the IceCube experiment. This spectrum has been observed to exhibit a dip-like feature at $E_{\nu} \sim 0.2-1 \, {\rm PeV}$, which could plausibly arise from the physics of the sources themselves, but could also be the consequence of high-energy neutrinos resonantly scattering with the cosmic neutrino background, mediated by a $Z'$ with a mass on the order of $m_{Z'} \sim 10 \, {\rm MeV}$. In this study, we calculate the impact of such a $Z'$ on the high-energy neutrino spectrum for a variety of model parameters and source distributions. For couplings that can resolve the $g_{\mu}-2$ anomaly, we find that this model could self-consistently produce a spectral feature that is consistent with IceCube's measurement, in particular if the neutrinos observed by IceCube predominantly originate from high-redshift sources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Lens Modeling of STRIDES Strongly Lensed Quasars Using Neural Posterior Estimation

Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory’s Legacy Survey of Space and Time and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, γ lens , to be $\mathcal{M}_γ$ lens = 2.13 ± 0.06. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone toward cosmological inference with large samples of strongly lensed quasars.

79 ASTRONOMY AND ASTROPHYSICS↗