Results for “Parametric Capability” — Page 10

DOE OSTI2023

FY23 Report on Water NSTF Testing at Two-Phase Conditions: Off-normal Scenarios

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report details the experimental activities and upgrades performed during the program’s fifth year of water-based operation. The theme of this year of testing was examining off-normal scenarios. In practice, that meant examining scenarios with blockages or other design basis conditions. From a maintenance and capability standpoint that meant assessing valving options for implementing blockages, assessing drain and refill capabilities, and ensuring in all cases that operators could safely, efficiently, and repeatably bring the facility to an off-normal condition and return it to the normal operating state. In one notable case, no solution existed, nor was a solution commercially available, that met the facility needs for implementing blockages within the two-phase region of the chimney. A valve and feedthrough were custom-designed for NSTF to operate submerged in saturated water while being safely actuated from outside the tank. The fifth year of testing included eight matrix tests consisting of 162 hours of active heating, 9,735 kWh of electrical heating, with seven tests classified as Accepted per NQA-1 and one as Trending. All matrix tests were performed at two-phase flow conditions. The power parametric series from previous years was extended to include another scenario with a decay load equivalent to 1.75 MWt, full scale, further resolving the system flow oscillation response as a function of input power. A depletion test was performed to extend previous work and examine the system’s approach to stagnation and geysering that were not previously observed. This test began at an initial inventory level of 50%, prototypic power of 2.1 MWt, and employed an accelerated drain of 0.75 gpm plus boiloff until stagnation conditions were achieved and multiple geysering events were recorded. Further exploring the geysering phenomenon was a separate effects test examining static boiling in the risers and chimney. This static boiling test was performed at 48 kWe, a higher power than testing in the previous year, and included enhancements to the drain system to remove inventory in a better-controlled and more reliable manner. Two tests were performed, at prototypic conditions of 1.4 MWt and 2.4 MWt, respectively, examining the system response to throttled conditions in the singlephase region of the piping. A ball valve immediately before he riser inlet header was increasingly throttled, progressing the system through a repeatable pattern of two-phase oscillations and stable regimes before ultimately inducing stagnation and geysering. Concluding the year, two tests were performed at baseline conditions of 70% initial inventory level and 2.1 MWt of prototypic decay heat to examine both repeatability and the response to throttling within the two-phase region, at the inventory tank inlet using the custom valve mentioned above. The first such test repeated transient operating conditions and revealed significant sensitivity t

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Advanced Fuels Campaign Execution Plan

The Advanced Fuels Campaign (AFC) Execution Plan outlines the strategy, mission, scope, near-term and long-term goals, structure, and organization associated with nuclear fuels and materials research, development, and demonstration activities within the Department of Energy’s (DOE) Nuclear Fuel Cycle and Supply Chain (NFCSC) program. NFCSC has been given responsibility to identify and mature advanced fuel technologies for the DOE using a science-based approach, focused on developing a fundamental understanding of nuclear fuels and materials to drive development of integrated nuclear fuel and materials technology. This science-based approach combines theory, experiments, and multiscale modeling and simulation to achieve a predictive understanding of relevant behaviors ranging from fuel fabrication processes (and their resulting fuel microstructures) through fuel/cladding performance under irradiation (in contrast to more empirical, observation-based approaches frequently used in fuel performance modeling and fuel qualification). The traditional scope of AFC includes the evaluation and development of multiple fuel forms to support two fuel cycle options: once-through and full recycle. The word “fuel” is used generically to include conventional fuels, transmutation targets, and any associated cladding or duct materials. The once-through fuel cycle addresses advanced light water reactor fuels with enhanced performance, extended burnup, and reduced waste generation. In fiscal year (FY) 2012, AFC’s scope expanded to include research, development, and demonstration (RD&D) for light water reactor (LWR) fuels with enhanced accident tolerance. Fuel fabrication activities include the development of innovative methods to enhance process efficiencies, reduce waste, and improve control over as-fabricated fuel microstructural properties to achieve desired in-reactor performance. Using modern modeling and simulation approaches, the objective is to predict fresh fuel properties given the feedstock characteristics and fabrication process parameters. The performance-related activities include small-scale, in-reactor, and out-of-reactor phenomenological testing (distinct from, but synergistic with, integral prototypic testing) and extensive, quantitative characterization (focusing on characterization of fuel and cladding materials at the scale of microstructure) both before and after testing. Larger-scale, prototypic experiments are conducted in concert with phenomenological testing to drive a Fuel Development and Qualification program, incorporating a fundamental understanding of fuel behavior performance characteristics. Then, using the tools developed under the productive science-based approach, fuels will be optimized to meet specific performance requirements, thereby minimizing the need to repeatedly perform large-scale, integral experiments over a wide parametric range as a means of experimental exploration. Two significant initiatives are underway within AFC. First, a gap analysis completed in early FY 2019 identified critical irradiation testing needs that are lacking within the national light water reactor (LWR) fuels testbed since the shutdown of the Halden Reactor in 2018. The identified gaps are for instrumented, prototypic testing of LWR fuels, especially under boiling water reactor conditions, ramp conditions, and conditions leading to fuel failure; these needs exist for supporting current LWR fuels and their possible extension to higher burnups, but are especially urgent relative to near-term development and qualification of accident-tolerant fuels. Recommendations that resulted from the Halden Gap Analysis focused on enhancements at Advanced Test Reactor (ATR) and Transient Reactor Test Facility (TREAT) to fill gaps in testing capabilities relative to these needs. Second, a concerted effort to develop and demonstrate a systematic approach to accelerating the development, testing, and qualification of new fuel systems has been initiated. This is highlighted by a test strategy that combines the considerable advances in multiscale, mechanistic fuel modeling of recent years with a MiniFuel separate effects test program in the High Flux Isotope Reactor (HFIR) and a Fission Accelerated Steady-state Testing (FAST) semi-integral accelerated test program in ATR. This approach is being tested/demonstrated using the metallic fuel system, but if successful it is expected to be extensible to multiple fuel types and diverse applications. This document includes an overview of the NFCSC program, a definition of science-based development of nuclear fuels, near-term goals for Advanced LWR fuels (ALFs), and longer-term goals for Advanced Reactor Fuels (ARFs) RD&D. This includes the activities that will be conducted to achieve success toward the grand challenge, as well as the goals and milestones to be achieved over the next few decades of research and development. Long-term goals are based on the DOE Office of Nuclear Energ

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Advanced Fuels Campaign Execution Plan

The Advanced Fuels Campaign (AFC) Execution Plan outlines the strategy, mission, scope, near-term and long-term goals, structure, and organization associated with nuclear fuels and materials research, development, and demonstration activities within the Department of Energy’s (DOE) Nuclear Fuel Cycle and Supply Chain (NFCSC) program. NFCSC has been given responsibility to identify and mature advanced fuel technologies for the DOE using a science-based approach, focused on developing a fundamental understanding of nuclear fuels and materials to drive development of integrated nuclear fuel and materials technology. This science-based approach combines theory, experiments, and multiscale modeling and simulation to achieve a predictive understanding of relevant behaviors ranging from fuel fabrication processes (and their resulting fuel microstructures) through fuel/cladding performance under irradiation (in contrast to more empirical, observation-based approaches frequently used in fuel performance modeling and fuel qualification). The traditional scope of AFC includes the evaluation and development of multiple fuel forms to support two fuel cycle options: once-through and full recycle. The word “fuel” is used generically to include conventional fuels, transmutation targets, and any associated cladding or duct materials. The once-through fuel cycle addresses advanced light water reactor fuels with enhanced performance, extended burnup, and reduced waste generation. In fiscal year (FY) 2012, AFC’s scope expanded to include research, development, and demonstration (RD&D) for light water reactor (LWR) fuels with enhanced accident tolerance. Fuel fabrication activities include the development of innovative methods to enhance process efficiencies, reduce waste, and improve control over as-fabricated fuel microstructural properties to achieve desired in-reactor performance. Using modern modeling and simulation approaches, the objective is to predict fresh fuel properties given the feedstock characteristics and fabrication process parameters. The performance-related activities include small-scale, in-reactor, and out-of-reactor phenomenological testing (distinct from, but synergistic with, integral prototypic testing) and extensive, quantitative characterization (focusing on characterization of fuel and cladding materials at the scale of microstructure) both before and after testing. Larger-scale, prototypic experiments are conducted in concert with phenomenological testing to drive a Fuel Development and Qualification program, incorporating a fundamental understanding of fuel behavior performance characteristics. Then, using the tools developed under the productive science-based approach, fuels will be optimized to meet specific performance requirements, thereby minimizing the need to repeatedly perform large-scale, integral experiments over a wide parametric range as a means of experimental exploration. Two significant initiatives are underway within AFC. First, a gap analysis completed in early FY 2019 identified critical irradiation testing needs that are lacking within the national light water reactor (LWR) fuels testbed since the shutdown of the Halden Reactor in 2018. The identified gaps are for instrumented, prototypic testing of LWR fuels, especially under boiling water reactor conditions, ramp conditions, and conditions leading to fuel failure; these needs exist for supporting current LWR fuels and their possible extension to higher burnups, but are especially urgent relative to near-term development and qualification of accident-tolerant fuels. Recommendations that resulted from the Halden Gap Analysis focused on enhancements at Advanced Test Reactor (ATR) and Transient Reactor Test Facility (TREAT) to fill gaps in testing capabilities relative to these needs. Second, a concerted effort to develop and demonstrate a systematic approach to accelerating the development, testing, and qualification of new fuel systems has been initiated. This is highlighted by a test strategy that combines the considerable advances in multiscale, mechanistic fuel modeling of recent years with a MiniFuel separate effects test program in the High Flux Isotope Reactor (HFIR) and a Fission Accelerated Steady-state Testing (FAST) semi-integral accelerated test program in ATR. This approach is being tested/demonstrated using the metallic fuel system, but if successful it is expected to be extensible to multiple fuel types and diverse applications. This document includes an overview of the NFCSC program, a definition of science-based development of nuclear fuels, near-term goals for Advanced LWR fuels (ALFs), and longer-term goals for Advanced Reactor Fuels (ARFs) RD&D. This includes the activities that will be conducted to achieve success toward the grand challenge, as well as the goals and milestones to be achieved over the next few decades of research and development. Long-term goals are based on the DOE Office of Nuclear Energ

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

A Python Library for Radiance Matrix-based Simulation Control and EnergyPlus Integration

Radiance matrix-based methods enable efficient parametric simulations, allowing users to vary sky conditions, fenestration systems, and other model parameters at a minimal cost to computation. However, the steep learning curve and complex workflow hinder the widespread adoption of matrix-based methods. The frads Python library with a series of command-line tools was developed to automate the entire matrix-based simulation process, lowering entry barriers and reducing human error. Co-simulation between EnergyPlus and Radiance was also enabled using the Python library from EnergyPlus. Key Innovations • Command-line based automation of Radiance matrix-based simulation methods • Python library facilitates broader adoption of Radiance matrix-based simulation methods • Radiance EnergyPlus run-time integration enabling the modeling of advanced control systems Practical Implications The frads library, with associated command-line tools, provides practitioners with the capability to easily adopt and use Radiance matrix-based simulation methods for various daylighting, solar control, and energy-related evaluations. Frads' current form is designed for 1) users familiar with a command-line interface and 2) software developers to integrate the matrix-based methods into existing software packages.

Wang, Taoning↗

DOE OSTI2022

Accurate and Efficient Parametric Model-Order Reduction for Turbulent Thermal Transport

This project produced new algorithms and software for low-cost exploration of turbulent thermal-fluids behavior under parametric variation by using reduced-order models (ROMs). The ROMs numerically solve the governing equations for fluid motion by using a small set of basis functions (typically, N=20-200 modes) to represent the solution. The base modes are computed as optimal combinations of solutions from expensive high-fidelity "anchor-point" solutions involving millions of unknowns, which are typically generated by solving the full Navier-Stokes equations on a supercomputer. The ROM solution is itself a combination of the base modes, where the basis coefficients are determined by evolving an NxN system of nonlinear equations. The overarching idea is to use the inexpensive ROM to predict solutions under conditions where the parameters differ from the anchor-point conditions. Several ingredients are required to make ROMs useful for thermal hydraulics analysis. These include: a stable and accurate ROM that is capable of reproducing the large-scale dynamics of turbulent flow, error indicators than can guide the choice of anchor points, and low-cost mechanisms for evaluating nonlinear terms in the reduced equations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

An absolute recoil-carbon polarimeter for polarized light-ions at the BNL Booster

We describe a recoil-carbon polarimeter for the BNL Booster synchrotron capable of measuring the transverse polarization of both the polarized proton beam and the polarized 3 He ion beam throughout the acceleration cycle. The instrument addresses a critical gap in the BNL polarized beam program: no independent polarization measurement currently exists in the injector chain between the 200 MeV Linac polarimeter and the AGS. The proposed new polarimeter will provide continuous, absolute polarization measurements at the Booster stage, supplying additional anchor points for the polarization transmission through the accelerator chain. The polarimeter 9 is based on elastic pC and 3 He C scattering from a thin internal carbon target, with detection of the recoil 12 C nucleus in a fixed six-station silicon detector ring. A single detector geometry provides continuous kinematic coverage from injection to extraction for both beam species without mechanical adjustment. For pC scattering, the available data and the Bonin parametrization provide a well-established basis for estimating the analyzing power and figure of merit over the Booster energy range; for 3 He C the sole experimental anchor is a measurement at 443 MeV, and the Booster polarimeter itself is identified as the instrument to map the analyzing power across the remainder of the acceleration ramp by a ramp-and-return calibration that anchors the absolute 3 He polarization scale to a few percent. Provision for future deuteron polarimetry is incorporated in the chamber design without modification to the existing geometry. With appropriate detector upgrades, the same recoil-carbon technique can be extended to deuteron beams, for which analyzing power data are already available, and, with future analyzing-power measurements, also to 6 Li and 7 Li beams, making the polarimeter a versatile instrument for the full range of light polarized ion species anticipated at BNL.

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗

DOE OSTI2020

Vitality models found useful in modeling tag-failure times in acoustic-tag survival studies

Abstract Acoustic telemetry studies often rely on the assumption that premature tag failure does not affect the validity of inferences. However, in some cases this assumption is possibly or likely invalid and it is necessary to apply a correction to estimation procedures. The question of which approaches and specific models are best suited to modeling acoustic tag failures has received little research attention. In this short communication, we present a meta-analysis of 42 acoustic tag-life studies, originally used to correct survival studies involving outmigrating juvenile salmonids in the Columbia/Snake river basin. We compare the performance of nine alternative parametric models including common failure–time/survival models and the vitality models of Li and Anderson Theor Popul Biol 76:118–131, (2009) and Demogr Res 28:341–372, (2013). The tag-life studies used acoustic tags from three different tag manufacturers, had expected lifetimes between 12 and 61 days, and had dry weights ranging from 0.22 to 1.65 g. In 57% of the cases, the vitality models of Li and Anderson Theor Popul Biol 76:118–131, (2009) and Demogr Res 28:341–372, (2013) fit the tag-failure times best. The vitality models were also the second-best choices in 17% of the cases. Together, the vitality models, log-logistic, (19%), and gamma models (14%) accounted for 90% of the models selected. Unlike more traditional failure–time models (e.g., Weibull, Gompertz, gamma, and log-logistic), the vitality models are capable of characterizing both the early onset of tag failure due to manufacturing errors and the anticipated battery life. We provide further guidance on appropriate sample sizes (50–100 tags) and procedures to be considered when applying precise tag-life corrections in release–recapture survival studies.

Skalski, John R. (ORCID:0000000270702505)↗

DOE OSTI2023

Improved microwave SQUID multiplexer readout using a kinetic-inductance traveling-wave parametric amplifier

We report on the use of a kinetic-inductance traveling-wave parametric amplifier (KITWPA) as the first amplifier in the readout chain of a microwave superconducting quantum interference device (SQUID) multiplexer (μmux). This μmux is designed to multiplex signals from arrays of low temperature detectors, such as superconducting transition-edge sensor microcalorimeters. When modulated with a periodic flux-ramp to linearize the SQUID response, the flux noise improves, on average, from 1.6 μ $\Phi$ 0 /$\sqrt{Hz}$ with the KITWPA off to 0.77 μ $\Phi$ 0 /$\sqrt{Hz}$ with the KITWPA on. When statically biasing the μmux to the maximally flux-sensitive point, the flux noise drops from 0.45 μ $\Phi$ 0 /$\sqrt{Hz}$ to 0.2 μ $\Phi$ 0 /$\sqrt{Hz}$⁠. We validate this new readout scheme by coupling a transition-edge sensor microcalorimeter to the μmux and detecting background radiation. In conclusion, the combination of μmux and KITWPA provides a variety of new capabilities including improved detector sensitivity and more efficient bandwidth utilization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Coherent terahertz radiation with 2.8-octave tunability through chip-scale photomixed microresonator optical parametric oscillation

High-spectral-purity frequency-agile room-temperature sources in the terahertz spectrum are foundational elements for imaging, sensing, metrology, and communications. Here we present a chip-scale optical parametric oscillator based on an integrated nonlinear microresonator that provides broadly tunable single-frequency and multi-frequency oscillators in the terahertz regime. Through optical-to-terahertz down-conversion using a plasmonic nanoantenna array, coherent terahertz radiation spanning 2.8-octaves is achieved from 330 GHz to 2.3 THz, with ≈20 GHz cavity-mode-limited frequency tuning step and ≈10 MHz intracavity-mode continuous frequency tuning range at each step. By controlling the microresonator intracavity power and pump-resonance detuning, tunable multi-frequency terahertz oscillators are also realized. Furthermore, by stabilizing the microresonator pump power and wavelength, sub-100 Hz linewidth of the terahertz radiation with 10 –15 residual frequency instability is demonstrated. The room-temperature generation of both single-frequency, frequency-agile terahertz radiation and multi-frequency terahertz oscillators in the chip-scale platform offers unique capabilities in metrology, sensing, imaging and communications.

42 ENGINEERING↗

DOE OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part II: Continuum-Mechanics Models

The design of efficient material-handling systems for milled lignocellulosic biomass is challenging due to their complex particle morphologies and frictional interactions. Computational modeling, including the discrete element method (DEM) and continuum-based finite-element/volume methods, may offer scientific insight and predictive capabilities for the flow of milled biomass in hoppers and feeders. Herein, this article (Part II) presents a review of current state-of-the-art continuum models for the flow of milled biomass, whereas DEM models are reviewed in a companion article (Part I). Advances of numerical methods to solve the global governing equations are discussed first, followed by a comprehensive review of constitutive models for granular materials, including Drucker–Prager, hypoplastic, Cambridge-type, inertial-rheology, and nonlocal granular fluidity models. Specifically, we provide in-depth discussion on the suitability of those models for milled lignocellulosic biomass materials in terms of nonlinear elasticity, dependence of flow strength on pressure, density and shear rate, and compaction (dilation) associated with hardening (softening). Furthermore, our study shows that, despite the recent advances in continuum granular flow modeling, the most suitable constitutive models still need further development to account for material parametrization, multiflow regimes, and multiscale behavior before they can be reliably used to optimize the design and operation of biomass handling systems.

09 BIOMASS FUELS↗

DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗

DOE OSTI2022

The revised FLORIDyn model: implementation of heterogeneous flow and the Gaussian wake

In this paper, a new version of the FLOw Redirection and Induction Dynamics (FLORIDyn) model is presented. The new model uses the three-dimensional parametric Gaussian FLORIS model and can provide dynamic wind farm simulations at a low computational cost under heterogeneous and changing wind conditions. Both FLORIS and FLORIDyn are parametric models which can be used to simulate wind farms, evaluate controller performance and can serve as a control-oriented model. One central element in which they differ is in their representation of flow dynamics: FLORIS neglects these and provides a computationally very cheap approximation of the mean wind farm flow. FLORIDyn defines a framework which utilizes this low computational cost of FLORIS to simulate basic wake dynamics. This is achieved by creating so-called observation points (OPs) at each time step at the rotor plane which inherit the turbine state. In this work, we develop the initial FLORIDyn framework further considering multiple aspects. The underlying FLORIS wake model is replaced by a Gaussian wake model. The distribution and characteristics of the OPs are adapted to account for the new parametric model but also to take complex flow conditions into account. To achieve this, a mathematical approach is developed to combine the parametric model and the changing, heterogeneous world conditions and link them with each OP. We also present a computationally lightweight wind field model to allow for a simulation environment in which heterogeneous flow conditions are possible. FLORIDyn is compared to Simulator for Offshore Wind Farm Applications (SOWFA) simulations in three- and nine-turbine cases under static and changing environmental conditions. The results show a good agreement with the timing of the impact of upstream state changes on downstream turbines. They also show a good agreement in terms of how wakes are displaced by wind direction changes and when the resulting velocity deficit is experienced by downstream turbines. A good fit of the mean generated power is ensured by the underlying FLORIS model. In the three-turbine case, FLORIDyn simulates 4 s simulation time in 24.49 ms computational time. The resulting new FLORIDyn model proves to be a computationally attractive and capable tool for model-based dynamic wind farm control.

17 WIND ENERGY↗

DOE OSTI2025

Five Year Comparison of Mixing Height Determinations at the Savannah River Site

Air quality dispersion modeling is performed for the Savannah River Site (SRS) to demonstrate compliance with applicable regulations. The AMS/EPA Regulatory Model (AERMOD) modeling system is an EPA recommended model for air quality applications with a data preprocessor (AERMET) to incorporate meteorological data collected on site. AERMET parameterizes or calculates meteorological variables that are not directly measured onsite. One of the parameters estimated by AERMET is the atmospheric mixing height. While the mixing height is not currently a measurement input into AERMET, SRS has the capability to measure the local mixing height. The Savannah River National Laboratory (SRNL) operates a Vaisala CL31 Lidar Ceilometer which estimates mixing height from aerosol backscatter. This study compares the parameterized mixing height from AERMET to the ceilometer estimated mixing height for the current regulatory period at SRS incorporating data from 2015-2019. Results from this study showed the average daily minimum values (morning) from AERMET were an order of magnitude lower than the commonly used Holzworth (1972) method and the ceilometer estimated mixing heights. Additionally, on average, the ceilometer exhibited a daily maximum mixing height value that occurred 1-3 hours later than the AERMET estimated maximum. This difference is likely due to the nighttime atmospheric mixing height assumptions and calculations used by AERMET. The AERMET algorithm cuts off mixing height growth at sunset while the ceilometer data show ongoing evening convection typical of the southeastern United States. These results suggest that the AERMET parametrization scheme assumptions may not be representative of a forested landscape and evening convection which could account for more mixing overnight. The results obtained in this study are significant for air dispersion modeling applications for regulatory purposes and worker safety. Mixing height can impact model estimated pollutant concentrations. A greater mixing height will provide more volume for pollutant dispersion. This report documents efforts to quantify the dependence of mixing height inputs toward a conservative estimated pollutant concentration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Robust implementation of generative modeling with parametrized quantum circuits

Although the performance of hybrid quantum-classical algorithms is highly dependent on the selection of the classical optimizer and the circuit ansätze, a robust and thorough assessment on-hardware of such features has been missing to date. From the optimizer perspective, the primary challenge lies in the solver’s stochastic nature, and their significant variance over the random initialization. Therefore, a robust comparison requires one to perform several training curves for each solver before one can reach conclusions about their typical performance. Since each of the training curves requires the execution of thousands of quantum circuits in the quantum computer, such a robust study remained a steep challenge for most hybrid platforms available today. In this work, we leverage on Rigetti’s Quantum Cloud Services (QCS™) to overcome this implementation barrier, and we study the on-hardware performance of the data-driven quantum circuit learning (DDQCL) for three different state-of-the-art classical solvers, and on two-different circuit ansätze associated to different entangling connectivity graphs for the same task. Additionally, we assess the gains in performance from varying circuit depths. To evaluate the typical performance associated with each of these settings in this benchmark study, we use at least five independent runs of DDQCL towards the generation of quantum generative models capable of capturing the patterns of the canonical Bars and Stripes dataset. In this experimental benchmarking, the gradient-free optimization algorithms show an outstanding performance compared to the gradient-based solver. In particular, one of them had better performance when handling the unavoidable noisy objective function to be minimized under experimental conditions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

VADER: A Tool for Criticality Safety Validation

The purpose of criticality safety is to prevent any inadvertent criticality from occurring during the handling or storage of fissile material. Calculations are frequently used to demonstrate that a sufficient subcritical margin exists. Validation is a key aspect of the evaluation process, establishing the suitability, accuracy, and associated uncertainty of the computational method and data to be used for the intended application. The validation process is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 in order to investigate the parameters at which such a critical condition is achieved. The validation parameters that are traditionally applied to safety analysis calculations are the bias and the bias uncertainty . The bias is the deviation of the average k eff of the validation suite from unity. The bias uncertainty accounts for the statistical uncertainty in the bias based on the standard deviation, sample size, and distribution of k eff values of the validation suite. The values of bias and bias uncertainty ensure that the systems predicted to be subcritical by the computational method will indeed be subcritical. The bias and bias uncertainty are often combined to determine an upper subcritical limit (USL) or computational margin that can then be applied to safety analysis calculations. Many methods have been developed by different organizations to calculate the bias and bias uncertainty for various types of criticality analyses. Each of these methods typically requires that the validity of various underpinning statistical assumptions be confirmed to demonstrate that the method is appropriate for the analysis of a given validation suite. An example of the validation decision making flow is shown in Fig.1. As shown in Fig. 1, the analyst performing the validation fits a trend line to the data and performs a test to determine if the trend was a statistically better representation of the data than if it were treated as an uncorrelated sample. If the trend line is a better representation of the data, then the analyst uses any one of a number of trending techniques to determine the bias and bias uncertainty. If a trend is not an appropriate representation of the data, then the analyst proceeds to perform a normality assessment for the data. If the normal assumption can be shown to be acceptable, then the analyst calculates the bias and bias uncertainty with the parametric technique. If the assumption of normality cannot be justified, then the nonparametric technique is used. Once the decision flow has been followed and the appropriate technique has been selected, the bias and bias uncertainty is typically combined with an administrative margin to determine a USL below which calculated values of k eff for safety analysis models can be considered subcritical. The calculations used in each decision are often performed with spreadsheets or with small programs available at various sites performing criticality analyses. Expertise in understanding and interpreting the results must be maintained to perform these calculations. This can often be an error-prone process. Oak Ridge National Laboratory (ORNL) is currently developing the Validation and Data Evaluation Resource (VADER) to simplify and automate the criticality safety validation process and to provide a software quality assurance pedigree to the calculational methods used. This paper discusses the use of the Fulcrum user interface with VADER, the anticipated initial capabilities of VADER to perform validation analyses, and the output from the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

On the effect of mixing-driven vaporization in a homogeneous relaxation modeling framework

The homogeneous relaxation model (HRM) is one of the most widely used models to describe the liquid–gas phase transition in multiphase flows due to the occurrence of cavitation. However, in its original formulation, the HRM does not account for the presence of ambient gas species, which generally limits its applicability to the injector's internal flow where ambient gases are negligible. In this work, a mixing-driven vaporization (MDV) model was developed to extend the capability of the HRM in handling the mixing effect in the regions external to the nozzle, where vapor–liquid equilibrium for multi-species mixtures of fuel and ambient gas is considered. Herein, to assess the model performance, simulations of the Engine Combustion Network's Spray G injector were performed with the HRM and the MDV model under both flash-boiling and evaporating conditions. It was found that the MDV model led to a better match against x-ray measurements of fuel density in the near-nozzle region. In contrast to the HRM, the MDV model was able to reproduce the vaporization process in the mixing zone at the edge of the fuel jet, which aligns with the expected physics. This resulted in substantial differences in the prediction of other flow characteristics such as mixture temperature and pressure. Furthermore, this work demonstrates that evaporation timescales have a considerable effect on the MDV model's predictions, as shown by a parametric study in which a time factor was introduced to mimic the effect of different timescales due to different phase change mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗