SEARCH · Engineering Papers

Results for “probability generating function”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Development of a Particle Flow Control Mechanism

The next generation of Concentrated Solar Power (CSP) systems are utilizing solid particles as the heat transfer medium. The state-of-the-art particle flow control mechanism currently utilized is the slide-gate mechanism which functions by linearly actuating a gate across a tear-drop shaped opening. The sensitivity, which is a measure of the change in mass flow rate per unit movement of the slide-gate, is 0.354 g/s-mm[1]. The slide-gate’s main shortcomings include particle lodging and particle leakage between the plates, and difficulty actuating under the pressure from the particles. The goal of a new device is to minimize or eliminate these flaws. The new proposed device functions similar to a chuck mechanism in a drill where a rotation of an outer ring produces a change in outlet area through moving “jaws” diagonally. This mechanism aims to solve the issues of the slide-gate by having the jaws on the outside of the particle flow volume to reduce the probability of particle lodging, and the jaw tips are designed to overlap as they actuate to reduce particle leakage. Additionally, the vertical orientation of the device as well as the inclined movement of the jaws allows for easier actuation of the device under pressure. Experiments involving 3D printed prototypes have been carried out, producing a max mass flowrate of about 50 g/s and a sensitivity curve with the equation 𝑦 (𝑔/𝑠) = 0.0851𝑥^2(𝑔/𝑠−𝑚𝑚2) +3.6993𝑥(𝑔/𝑠−𝑚𝑚). These experiments validated the device’s capability to address the shortcomings of the slide-gate. Currently, an aluminum model of the device is being tested at room temperature, and a stainless-steel model is being developed to be tested at 700oC, which is the expected operating temperature of the CSP system. Furthermore, simulations are on-going to predict the mass flowrate of different device configurations at room and high temperatures.

14 SOLAR ENERGY↗

DOE OSTI2024

The reproduction number and its probability distribution for stochastic viral dynamics

We consider stochastic models of individual infected cells. The reproduction number, R, is understood as a random variable representing the number of new cells infected by one initial infected cell in an otherwise susceptible (target cell) population. Variability in R results partly from heterogeneity in the viral burst size (the number of viral progeny generated from an infected cell during its lifetime), which depends on the distribution of cellular lifetimes and on the mechanism of virion release. We analyse viral dynamics models with an eclipse phase: the period of time after a cell is infected but before it is capable of releasing virions. The duration of the eclipse, or the subsequent infectious, phase is non-exponential, but composed of stages. We derive the probability distribution of the reproduction number for these viral dynamics models, and show it is a negative binomial distribution in the case of constant viral release from infectious cells, and under the assumption of an excess of target cells. In a deterministic model, the ultimate in-host establishment or extinction of the viral infection depends entirely on whether the mean reproduction number is greater than, or less than, one, respectively. Here, the probability of extinction is determined by the probability distribution of R, not simply its mean value. In particular, we show that in some cases the probability of infection is not an increasing function of the mean reproduction number.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

All-electron molecular tunnel ionization based on the weak-field asymptotic theory in the integral representation

Tunnel ionization (TI) underlies many important ultrafast processes, such as high-harmonic generation andstrong-field ionization. Among the existing theories for TI, many-electron weak-field asymptotic theory (ME-WFAT) is by design capable of accurately treating many-electron effects in TI. An earlier version of ME-WFATrelied on an accurate representation of the asymptotic tail of the orbitals, which hindered its implementation inGaussian-basis-set-based quantum chemistry programs. In this work, we reformulate ME-WFAT in the integralrepresentation, which makes the quality of the asymptotic tail much less critical, hence greatly facilitating itsimplementation in standard quantum chemistry packages. The integral reformulation introduced here is thereforemuch more robust when applied to molecules with arbitrary geometry. Here, we present several case studies, amongwhich is the CO molecule where some earlier theories disagree with experiments. Here we find that ME-WFATproduces the largest ionization probability when the field points from C to O, as experiments suggest. Anattractive feature of ME-WFAT is that it can be used with various types of multielectron methods whether ofdensity functional or multiconfiguration types, this inturn facilitates tunnel ionization calculation in systems exhibiting a strong multireference character.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Ecovoltaic solar energy development can promote grassland bird communities

Ecologically informed photovoltaic (PV) developments that co-prioritize PV electricity generation with ecosystem function (‘ecovoltaics’) have emerged as a promising land sharing strategy to minimize ecological conflicts associated with PV solar energy development. While habitat-focused ecovoltaic designs can conceptually benefit biodiversity by offsetting or enhancing impacts of PV development, foundational field research is needed to examine how wildlife respond to these novel ecosystems. We conducted passive acoustic monitoring (PAM) in 2023 and 2024 at 13 solar facilities and paired control sites to investigate avian community responses to ecovoltaic facilities in the Midwestern United States. Compared to control sites (row crop agricultural fields), we found that ecovoltaic sites supported more grassland bird species throughout a 17-week monitoring period between May and September. Grassland bird communities on ecovoltaic sites were also more stable than on agricultural controls, as measured by the Jaccard dissimilarity index. We also used PAM-based weekly species occurrences in an occupancy-modelling framework to investigate the influence of PV development and other landscape variables on grassland bird occupancy. 10 out of 13 modelled grassland bird species had greater predicted occupancy probabilities (ψ) on PV sites than control sites. Synthesis and applications. Our findings suggest that properly sited and developed ecovoltaic solar facilities in human altered landscapes can improve habitat for birds and other wildlife, but further research is needed to understand which species may benefit most from these novel ecosystems.

14 SOLAR ENERGY↗

DOE OSTI2024

Photometric redshifts probability density estimation from recurrent neural networks in the DECam local volume exploration survey data release 2

Photometric wide-field surveys are imaging the sky in unprecedented detail. These surveys face a significant challenge in efficiently estimating galactic photometric redshifts while accurately quantifying associated uncertainties. In this work, we address this challenge by exploring the estimation of Probability Density Functions (PDFs) for the photometric redshifts of galaxies across a vast area of 17,000 square degrees, encompassing objects with a median 5 σ point-source depth of g = 24.3, r = 23 . 9 , i = 23.5, and z = 22.8 mag. Our approach uses deep learning, specifically integrating a Recurrent Neural Network architecture with a Mixture Density Network, to leverage magnitudes and colors as input features for constructing photometric redshift PDFs across the whole DECam Local Volume Exploration (DELVE) survey sky footprint. Subsequently, we rigorously evaluate the reliability and robustness of our estimation methodology, gauging its performance against other well-established machine learning methods to ensure the quality of our redshift estimations. Our best results constrain photometric redshifts with the bias of − 0 . 0013 , a scatter of 0.0293, and an outlier fraction of 5.1%. These point estimates are accompanied by well-calibrated PDFs evaluated using diagnostic tools such as Probability Integral Transform and Odds distribution. We also address the problem of the accessibility of PDFs in terms of disk space storage and the time demand required to generate their corresponding parameters.We present a novel Autoencoder model that reduces the size of PDF parameter arrays to one-sixth of their original length, significantly decreasing the time required for PDF generation to one-eighth of the time needed when generating PDFs directly from the magnitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

FY24 Task 5: Leachate Disposition

Directly feeding sludge solids to the high-level waste (HLW) Waste Treatment Plant represents an alternative flowsheet seeking to initiate sludge processing as soon as possible. Key processing functions previously captured during baseline pretreatment operations include leaching and washing prior to solids concentration. These operations should be considered in the potential direct feed flowsheets to maximize waste feed loading, minimize HLW volume, and mitigate corrosion challenges associated with vitrification of high phosphate and fluoride concentrations. Additionally, single-shell tank (SST) retrievals and waste transfers to double-shell tanks (DSTs) in a direct feed flowsheet would likely also benefit from some level of leaching, washing, and solids concentration in order to reduce DST space and mission duration. These operations could occur in a new facility or potentially in available DSTs. If washing and leaching are utilized, an effective disposition pathway for the wash water and leachate solutions are needed. Three target species that benefit significantly from leaching and washing are phosphate, fluoride and aluminum. Phosphate (PO 4 3- ) and fluoride (F - ) can contribute substantially to the amount of carrier fluid needed for dissolution, and the resulting volume of liquid generated. Disposition of this retrieval solution should be evaluated in order to prevent crystallization of these anions throughout system processing. Since there is a high probability that any retrieval solutions will be at or near their PO 4 3- and F - solubility limits, evaporation or blending with a high Na supernate (>3.5 M) is not recommended for the wash water streams without a method to remove precipitants prior to solution disposal. Additionally, aluminum present in the southeast quadrant of the Hanford site represents roughly 60% of the waste solids in the initial processing tanks. These aluminum solids are in the form of gibbsite (Al(OH) 3 ) and can pose significant challenges for processing due to the fast-settling times and high solids loading associated with these materials. Easily remediated by caustic addition to the solids, these wash solutions could be processed through crystalline silicotitanate (CST) ion exchange columns to prepare the supernate solutions for disposition. The current target for feed conditions to the Low Activity Waste (LAW) melter are waste streams that contain nominally 5-6 M Na. Fractions within the tanks contain upwards of 0.2 M phosphate and fluoride in solution at 3.5 M Na. Concentrating these solutions above 5 M Na would result in an exceedance of the solubility limits, and potential for uncontrolled precipitation of the phosphate and fluoride crystal material. The resulting crystalline salt material is typically sodium fluoride phosphate, also referred to as natrophosphate (Na 7 FPO 4 ·19H 2 O). Salt phases are of importance in tank waste due to their chemical reactivity, which can result in precipitation, dissolution, or transformation, impacting any downstream processes (Bolling et al. 2020, Russell, Snow, and Peterson 2010). Salt generation and precipitation could pose challenges by causing system plugging and melter corrosion if left in the supernate stream, or limit sodium molarity of the supernate that would be accepted without incident in waste operations. To understand the impact of this salt generation, the crystallization of natrophosphate in multiple simulant feed matrices was studied to understand the implications of various tank waste supernate chemistries. Three matrices were examined: high PO 4 3- /low F - , low PO 4 3- /high F - , and an average matrix. Subsequent testing was performed with the average matrix with the inclusion of CsNO 3 , and a final run with the average matrix including CsNO 3 and a 137 Cs spike for tracer purposes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Reduced Order Modeling conditioned on monitored features for response and error bounds estimation in engineered systems

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to provide numerical error bounds associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Conditional VAEs↗

DOE OSTI2023

The DESI One-Percent Survey: Constructing Galaxy–Halo Connections for ELGs and LRGs Using Auto and Cross Correlations

In the current Dark Energy Spectroscopic Instrument (DESI) survey, emission line galaxies (ELGs) and luminous red galaxies (LRGs) are essential for mapping the dark matter distribution at z ~ 1. We measure the auto and cross correlation functions of ELGs and LRGs at 0.8 < z ≤ 1.0 from the DESI One-Percent survey. Following Gao et al., we construct the galaxy–halo connections for ELGs and LRGs simultaneously. With the stellar–halo mass relation for the whole galaxy population (i.e., normal galaxies), LRGs can be selected directly by stellar mass, while ELGs can also be selected randomly based on the observed number density of each stellar mass, once the probability P sat of a satellite galaxy becoming an ELG is determined. We demonstrate that the observed small scale clustering prefers a halo mass-dependent P sat model rather than a constant. With this model, we can well reproduce the auto correlations of LRGs and the cross correlations between LRGs and ELGs at r p > 0.1 Mpc h –1 . We can also reproduce the auto correlations of ELGs at r p > 0.3 Mpc h –1 (s > 1 Mpc h –1 ) in real (redshift) space. Although our model has only seven parameters, we show that it can be extended to higher redshifts and reproduces the observed auto correlations of ELGs in the whole range of 0.8 < z ≤ 1.6, which enables us to generate a lightcone ELG mock for DESI. With the above model, we further derive halo occupation distributions for ELGs, which can be used to produce ELG mocks in coarse simulations without resolving subhalos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Why is it still too warm or cold in my house? Examining the relationships between energy efficient capital and household energy insecurity

Here, this paper examines the relationships between energy efficient (EE) capital technology and household energy insecurity in the United States. The theoretical model of these relationships employs household production theory to capture the demand for and production of household energy services, and a stochastic production frontier approach to describe how having access to and the usage intensity of EE capital technology could help alleviate inefficiency in the production of household energy services. A working hypothesis formulated from the theoretical model posits that having EE capital technology in the home will reduce the level of household energy insecurity experienced. The extent of energy insecurity experienced is inferred from an energy insecurity index value assigned to each household, generated via the application of a dichotomous Rasch model to questions contained in the 2015 Residential Energy Consumption Survey. Noting the potential simultaneous relationship that exists between a household having access to and the usage intensity of EE capital technology and the experience of being energy insecure, an instrumental variables approach was employed to estimate a series of ordered logit models. Results suggest access to EE capital technology in the form of Energy Star® appliances, Energy Star® windows, or a SMART thermostat does not reduce the probability of experiencing a greater level of energy insecurity. Nor does the usage intensity of EE capital. Thus, policy instruments designed to alleviate household energy insecurity may need to go beyond simply helping households obtain EE capital technology.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Multiprobe cosmology from the abundance of SPT clusters and DES galaxy clustering and weak lensing

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the Universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3 × 2 ⁢pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining Λ cold dark matter (Λ ⁢CDM) parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure Ω m = 0.300 ± 0.017 and 𝜎 8 = 0.797 ± 0.026. Compared to constraints from Planck primary cosmic microwave background (CMB) anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 (1.2⁢𝜎) for the two-parameter difference. We further obtain 𝑆 8 ≡𝜎 8 ⁢(Ω m /0.3) 0.5 = 0.796 ± 0.013 which is lower than the Planck measurement at the 1.6⁢𝜎 level. The combined SPT cluster, DES 3 ×2 ⁢pt, and Planck datasets mildly prefer a nonzero positive neutrino mass, with a 95% upper limit ∑ 𝑚 𝜈 < 0.25 eV on the sum of neutrino masses. Assuming a 𝑤⁢CDM model, we constrain the dark energy equation of state parameter 𝑤 = −1.1⁢5$^{+0.23}_{−0.17}$ and when combining with Planck primary CMB anisotropies, we recover 𝑤 = −1.2⁢0$^{+0.15}_{−0.09}$, a 1.7⁢𝜎 difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology and our analysis paves the way for upcoming joint analyses of next-generation datasets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Seasonal Cost-Benefit Analysis of Automated Distribution Feeder Upgrades with Advanced Mitigation Technologies

The increasing deployment of distributed solar photovoltaics (DPV) to meet clean energy goals can trigger adverse grid operation issues, such as voltage excursions and the violation of thermal loading constraints of the power delivery elements (e.g., lines and transformers) on the evolving electricity infrastructure. Such integration issues would require distribution upgrades with associated costs to mitigate them and to maintain reliable and resilient grid operating conditions. Traditional distribution network upgrade approaches use a specific single snapshot analysis that is overly conservative. This study considers a multi-time point analysis to capture both moderate (probable bounds) and extreme grid operating conditions using time points such as minimum load with minimum photovoltaics (PV), maximum load with maximum PV, maximum load with minimum PV, and minimum load with maximum PV. Further, this study investigates seasonal variation impacts and associated distribution upgrade costs for a spring season case (March, representing a low load and high PV scenario) and a summer case (July, representing a high load and high PV scenario). Such seasonal analysis will allow system operators to characterize upgrade requirements and associated costs across various periods. Because the spatial distribution of DPV can impact upgrade and associated costs, this study investigates three common DPV deployment scenarios - randomly deployed, close to the substation, and far from the substation - at different penetration levels. Apart from spatial distribution impacts, this project evaluates the techno-economic impacts of the nodal photovoltaic penetration factor (NPPF) for generating the various DPV deployment scenarios at increasing penetration levels. This project investigates the impact of varying nodal PV-to-load ratios using conservative and extreme NPPF values of 3 and 10, respectively. This study investigates the deployment of traditional infrastructure upgrade strategies, such as installing new voltage regulating equipment, transformers and lines replacements, and the activation of advanced inverter functionality (e.g., autonomous volt/VAR) in expanding PV hosting capacity. Existing DPV systems are assumed to operate with the legacy unity power factor, and we considered the possibility of retrofitting such systems with the activation of volt/VAR control as integration standards and regulations continue to evolve to allow such functions. The cost-benefit analysis metrics used in study include distribution upgrade costs, average cost per watt of the upgrade cost, average marginal cost per watt of the upgrade cost, and power losses.

14 SOLAR ENERGY↗

DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Overview of the fusion nuclear science facility, a credible break-in step on the path to fusion energy

The Fusion Nuclear Science Facility (FNSF) is examined here as part of a two step program from ITER to commercial power plants. This first step is considered mandatory to establish the materials and component database in the real fusion in-service environment before proceeding to larger electricity producing facilities. The FNSF can be shown to make tremendous advances beyond ITER, toward a power plant, particularly in plasma duration and fusion nuclear environment. A moderate FNSF is studied in detail, which does not generate net electricity, but does reach the power plant blanket operating temperatures. The full poloidal Dual Coolant Lead Lithium (DCLL) blanket is chosen, with alternates being the Helium Cooled Lead Lithium (HCLL) and Helium Cooled Ceramic Breeder/Pebble Bed (HCCB/PB). Several power plant relevant choices are made in order to follow the philosophy of targeted technologies. Any fusion core component must be qualified by fusion relevant neutron testing and highly integrated non-nuclear testing before it can be installed on the FNSF in order to avoid the high probability of constant failures in a plasma-vacuum system. A range of missions for the FNSF, or any fusion nuclear facility on the path toward fusion power plants, are established and characterized by several metrics. A conservative physics strategy is pursued to accommodate the transition to ultra-long plasma pulses, and parameters are chosen to represent the power plant regime to the extent possible. An operating space is identified, and from this, one point is chosen for further detailed analysis, with R = 4.8 m, a = 1.2 m, IP = 7.9 MA, BT = 7.5 T, βN Gr = 0.9, fBS = 0.52, q95 = 6.0, H98 ∼1.0, and Q = 4.0. The operating space is shown to be robust to parameter variations. A program is established for the FNSF to show how the missions for the facility are met, with a He/H, a DD and 5 DT phases. The facility requires ∼25 years to complete its DT operation, including 7.8 years of neutron production, and the remaining spent on inspections and maintenance. The DD phase is critical to establish the ultra-long plasma pulse lengths. The blanket testing strategy is examined, and shows that many sectors have penetrations for heating and current drive (H/CD), diagnostics, or Test Blanket Modules (TBMs). The hot cell is a critical facility element in order for the FNSF to perform its function of developing the in-service material and component database. The pre-FNSF R&D is laid out in terms of priority topics, with the FNSF phases driving the time-lines for R&D completion. A series of detailed technical assessments of the FNSF operating point are reported in this issue, showing the credibility of such a step, and more detailed emphasis on R&D items to pursue. These include nuclear analysis, thermo-mechanics and thermal-hydraulics, liquid metal thermal hydraulics, transient thermo-mechanics, tritium analysis, maintenance assessment, magnet specification and analysis, materials assessments, core and scrape-off layer (SOL)/divertor plasma examinations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗