Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “highly deviated”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A new generation of uranium coincidence fast neutron collars for assay of LWR fresh fuel assemblies

The active uranium neutron coincidence collar provides a means of non-destructively assaying the fissile linear density of Light Water Reactor fresh fuel assemblies containing low enriched uranium. These neutron collars can operate in two modes: a thermal and a fast mode. In fast mode, a neutron collar has an added cadmium (Cd) liner in the sample cavity of the detector to reduce the impact of the burnable poison (thermal neutron absorber) on the detector signal (doubles). The main advantage for operating in fast mode is a detected signal that is less dependent of the burnable neutron poison content and thus less dependent on facility operator declarations. The drawback is that operating in fast mode requires a longer measurement time (~hour vs tens of minutes for thermal mode) to achieve the statistically needed precision in the measurements. The trend in the modern reactor fuel assemblies is moving to higher burnup by using higher initial enrichment and, consequently, a higher number of burnable poison rods to compensate the initial neutron reactivity. The increase of the burnable poison loading has motivated the development of a new generation of high efficiency fast neutron collars to allow practical measurements in-field by nuclear inspectors. This paper describes the development and performance evaluation of a new generation of neutron collars, for both boiling water reactor (BWR) and pressurized water reactor (PWR) fuels, jointly developed between the US Department of Energy, through Los Alamos National Laboratory, and the Euratom Safeguards Directorate of the European Commission. In this work, we present here calibrations with reference fuel assemblies at Los Alamos National Laboratory as well as the results of in-field measurement campaigns in fuel fabrication plants with modern commercial fuel assemblies. The experimental results show that a typical PWR verification can be made in a total time of 30 min with an uncertainty in the measured mass of 2% at one standard deviation (1σ). A BWR verification can be made in 47 min with an uncertainty in the measured mass of 1.9% at 1σ, or a total time of 20 min with 1σ uncertainty in the measured mass of 2.5%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Magnetic field-temperature phase diagram of the spin-$\frac{1}{2}$ triangular lattice antiferromagnet KYbSe 2

A quantum spin liquid (QSL) is a state of matter characterized by fractionalized quasiparticle excitations, quantum entanglement, and a lack of long-range magnetic order. However, QSLs have evaded definitive experimental observation. Several Yb 3+ -based triangular lattice antiferromagnets with effective 𝑆 = $\frac{1}{2}$ have been suggested to stabilize the QSL state as the ground state. Here, in this work, we build a comprehensive magnetic temperature phase diagram of a high-quality single crystalline KYbSe 2 via heat capacity and magnetocaloric effect down to 30 mK with magnetic field applied along the 𝑎 axis. At zero magnetic field, we observe the magnetic long-range order at 𝑇 N =0.29 K entering 120 degrees ordered state in heat capacity, consistent with neutron scattering studies. Analysis of the low-temperature (𝑇) specific heat (𝐶) at zero magnetic field indicates linear 𝑇 dependence of 𝐶/𝑇 and a broad hump of 𝐶/𝑇 in the proximate QSL region above 𝑇 N . By applying magnetic field, we observe the up-up-down phase with 1/3 magnetization plateau and oblique phases, in addition to two new phases. These observations strongly indicate that while KYbSe 2 closely exhibits characteristics resembling an ideal triangular lattice, deviations may exist, such as the effect of the next-nearest-neighbor exchange interaction, calling for careful consideration for spin Hamiltonian modeling. Further investigations into tuning parameters, such as chemical pressure, could potentially induce an intriguing QSL phase in the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Breakup dynamics in a pressure-swirl injector for urea-water solution applications: A computational study

The co-optimization of in-cylinder combustion and after-treatment technology has become a major aspect in engine design and development, with the goal of meeting the increasingly restrictive emission regulations in the transportation industry. Selective Catalytic Reduction is a robust technology to control the emission of NO x , and the injection of urea in water solution is the exhaust tailpipe is a key aspect of its operation. The proposed work uses high-fidelity Computational Fluid Dynamics to characterize the atomization dynamics of the liquid jet in relevant cross-flow conditions. The study focuses on a commercial low-pressure (9 bar) pressure-swirl injector which is characterized in its internal geometry through high-resolution X-ray micro-computational tomography. The internal two-phase flow has been modeled according to the volume-of-fluid approach in a large eddy simulation framework and validated against near-nozzle X-ray radiography measurement. Moreover, characterizing the breakup dynamics for the swirling hollow cone formation, and assessing the influence of the cross-flow in the breakup dynamics was completed. The results have been reported proposing Re-Oh maps and probability density functions of the spray kinematics. Higher cross-flow momentum generates an increase in the jet intact length and a reduction of the liquid droplet diameters. The axial momentum of the jet is affected by the cross-flow already in the near-nozzle region, determining a relevant deviation of the spray velocities. In conclusion, this work aims to inform the initialization of Eulerian-Lagrangian spray models through the assignment of droplet kinematics and static one-way coupling between volume-of-fluid results and Lagrangian spray parcels, to be used for system-size domain simulations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

A propagation-based fault detection and discrimination method and the optimization of sensor deployment

Industrial processes can be affected by faults having a serious impact on operation when not promptly detected and diagnosed. Here in this paper, a propagation-based fault detection and discrimination(PFDD) method is proposed to develop a strategy for fault diagnosis while in the design phase of a system. The PFDD method constructs the system model using the Integrated System Fault Analysis(ISFA) technique. Based on the system model, the propagation of hardware and software faults are simulated qualitatively. Given the results of the simulation, the process by which a fault propagates can be characterized using the qualitative features of system variables including the deviation of the system variables from their expected values, the variation of the system variables over time, and the order in which each variable is influenced during the propagation of the fault. The strategy by which a fault can be detected and discriminated is defined using those features. The PFDD method supports the detection and discrimination of faults in both steady states and transient states. Based on the PFDD method, the optimization of sensor deployment in a system is discussed. A brute force algorithm is developed to examine the system’s capability at diagnosing faults and the cost of sensor deployment for all possible configurations of sensors. The optimal sensor deployment strategy can be derived accordingly. However, the brute force method is only applicable to small-scale systems due to its high computational cost. A genetic algorithm is used to optimize sensor deployment in large-scale systems. The PFDD and sensor deployment optimization methods are applied to the Experimental Breeder Reactor II (EBR-II) for verification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗
DOE OSTI2023

The Spatial Clustering of ROSAT All-Sky Survey Active Galactic Nuclei. V. The Evolution of Broad-line AGN Clustering Properties in the Last 6 Gyr

This is the fifth paper in a series of investigations of the clustering properties of luminous, broad-emission-line active galactic nuclei (AGNs) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this work we measure the cross-correlation function (CCF) between RASS/SDSS Data Release 14 AGNs with the SDSS CMASS galaxy sample at 0.44 < z < 0.64. We apply halo occupation distribution (HOD) modeling to the CCF along with the autocorrelation function of the CMASS galaxies. We find that X-ray-selected and optically selected AGNs at 0.44 < z < 0.64 reside in statistically identical halos with a typical dark matter halo (DMH) mass of $M$$^{typ,AGN}_{DMH}$ ~ 10 12.7 h –1 M ⊙ . The acceptable HOD parameter space for these two broad-line AGN samples have only statistically marginal differences caused by small deviations of the CCFs in the one-halo-dominated regime on small scales. In contrast to optically selected AGNs, the X-ray AGN sample may contain a larger population of satellites at MDMH ~ 1013h–1M⊙. We compare our measurements in this work with our earlier studies at lower independent redshift ranges, spanning a lookback time of 6 Gyr. The comparison over this wider redshift range of 0.07 < z < 0.64 reveals (i) no significant difference between the typical DMH masses of X-ray-selected and optically selected AGNs, (ii) weak positive clustering dependencies of $M$$^{typ,AGN}_{DMH}$ with L X and M BH , (iii) no significant dependence of $M$$^{typ,AGN}_{DMH}$ on Eddington ratio, and (iv) the same DMH masses host more-massive accreting black holes at high redshift than at low redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Development & Experimental Validation of a Generalized Resistance-Capacitance Model for Numerical Simulation of Phase-Change Material Embedded Heat Exchangers

Latent heat thermal energy storage (LHTES) using phase change material (PCM) has attracted increased attention as a viable solution for overcoming the mismatch between energy supply and demand for renewable energy-based systems. PCM-embedded heat exchangers (PCM-HX) have the potential to significantly improve thermal performance due to high storage capacity and low temperature variation during the phase change process. Most models for simulating LHTES heat transfer use Computational Fluid Dynamics (CFD) simulations, which have high computational costs resulting from considering the complex and time-dependent physics relevant to PCM-HXs. In this paper, a Generalized Resistance Capacitance-based Model (GRCM) was developed to predict the thermal performance of arbitrary PCM-HXs in a computationally efficient manner without compromising modeling accuracy. The GRCM is exercised for three case studies: (i) verification for a single-slabbed finned PCM-HX, (ii) verification and validation for a copper foam/paraffin composite PCM-HX, and (iii) validation for a straight tube annular finned PCM-HX. The copper foam PCM-HX uses an electric heater at the top of HX, while the other two configurations utilize water as heat transfer fluid. For the single-slabbed finned PCM-HX melting case, the mean deviation in average PCM temperature predicted by the GRCM compared to the CFD model was between 0.56 – 0.73 K, with maximum temperature deviation of 2.68 K. For the HTF outlet temperature, the validation results showed that GRCM prediction matches very well with experimental data, with mean temperature deviation of 0.24 K during melting case, while for solidification case was 0.34 K. These results showcase the GRCM’s capability for accurately reproducing the thermal characteristics of PCM-HXs with considerably lower computational effort.

42 ENGINEERING↗
DOE OSTI2026

Grid and Wind Transient Ride Through Assessment of Type-5 Synchronous Wind Turbine

The integration of the inverter-based renewable energy resources is increasing and replacing conventional synchronous generators in the evolving bulk power system (BPS), impacting the system inertia along the way. Also, the limited over-current capability of power electronic converters continues to deteriorate the grid strength of the evolving power systems. Type-5 wind turbines have the potential to mitigate these problems because this configuration interfaces with the grid through a synchronous generator; hence, its operation and consequent grid impacts might mimic the conventional synchronous generator. The drivetrain of this type of turbines, however, is different from the more common Type-3 and Type-4 turbines. A hydraulic torque converter and gearbox with torque limiting feature are integral parts of a Type-5 wind turbine drivetrain, to isolate the wind speed variations at the wind turbine rotor from the grid connected generator side. This supports the synchronous operation of the generator. However, high fidelity models of Type-5 wind turbine drivetrain including these core components are not openly and widely available for grid integration and transient stability studies. Moreover, Type-5 turbines are not commonly deployed and hence its power generation and frequency stability characteristics are not widely publicized. This hinders appropriate assessment of Type-5 wind power plant's contribution to bulk grid resilience. This work provides transient stability assessment of a standalone Type-5 wind turbine. Normal operation and grid transient studies are performed to see the step load change and wind speed disturbance impacts on the Type-5 drivetrain and generation dispatch. Preliminary results show that Type-5 wind turbines can support load changes of between 0.6 p.u to 1.15 p.u. Beyond this point, the synchronization of the generator starts to deviate from its set point.

17 - WIND ENERGY↗
DOE OSTI2026

Pixelated plastic scintillator array manufacturing using fast-, photo-curable resin

Pixelated plastic scintillator arrays can serve as high efficiency and high resolution neutron imaging detectors. Manufacturing these arrays is intensive in both time and labor. This article presents a fabrication method based on additive manufacturing for two-dimensional plastic organic scintillator arrays using a custom-built automated assembly machine and a custom photocurable resin that has significant non-aromatic acrylate oligomer content. The process involves two main stages: fully autonomous production of one-dimensional layered arrays, followed by semi-autonomous cutting and stacking to form two-dimensional pixel arrays. One-dimensional arrays were manufactured at a rate of around 4 layers per hour with minimal defects and tight dimensional tolerances, while two-dimensional arrays up to 7 x 7 pixels and 70 mm in length were completed in approximately 3.5 hours. Final arrays exhibited dimensional deviations of less than 0.5 mm. Two-dimensional arrays read out by a multi-anode photomultiplier tube demonstrated per-pixel position resolution and pulse-shape discrimination, enabling gamma–neutron interaction separation in mixed radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Investigating the role of strangeness in baryon–antibaryon annihilation at the LHC

Annihilation dynamics plays a fundamental role in the baryon–antibaryon interaction ($B–\bar{B}$) at low-energy and its strength and range are crucial in the assessment of possible baryonic bound states. Experimental data on annihilation cross sections are available for the $p–\bar{p}$ system but not in the low relative momentum region. Data regarding the $B–\bar{B}$ interaction with strange degrees of freedom are extremely scarce, hence the modeling of the annihilation contributions is mainly based on nucleon–antinucleon ($N–\bar{N}$) results, when available. In this letter we present a measurement of the $p–\bar{p}$, $p–\bar{Λ}$⊕$p–{Λ}$ and $Λ–\bar{Λ}$ interaction using correlation functions in the relative momentum space in high-multiplicity triggered pp collisions at $\sqrt{s}$=13 TeV recorded by ALICE at the LHC. In the $p–\bar{p}$ system the couplings to the mesonic channels in different partial waves are extracted by adopting a coupled-channel approach with recent χEFT potentials. The inclusion of these inelastic channels provides good agreement with the data, showing a significant presence of the annihilation term down to zero momentum. Predictions obtained using the Lednický–Lyuboshits formula and scattering parameters obtained from heavy-ion collisions, hence mainly sensitive to elastic processes, are compared with the experimental $p–\bar{Λ}$⊕$p–{Λ}$ and $Λ–\bar{Λ}$ correlations. The model describes the $Λ–\bar{Λ}$ data and underestimates the $p–\bar{Λ}$⊕$p–{Λ}$ data in the region of momenta below 200 MeV/c. The observed deviation indicates a different contribution of annihilation channels to the two systems containing strange hadrons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Investigation of Nanoparticle Degradation in Hydrogen Fuel Cell Systems through Automated Electron Microscopy

Proton exchange membrane fuel cells (PEMFC) are promising devices for the deployment of hydrogen-powered heavy-duty vehicles, providing a higher efficiency for similar driving range and fueling time than the existing ones. However, PEMFCs still encounter durability challenges mainly due to catalyst degradation in the cathode. Mitigating these performance losses requires a better understanding of the degradation mechanisms under heavy-duty accelerated stress tests (ASTs) [1]. Scanning transmission electron microscopy (STEM) combined with energy dispersive X-ray spectroscopy (EDS) are key tools for the analysis of Pt and PtCo nanoparticle size, spatial distribution and composition [2]. Electron tomography is also used to determine the rate and type of degradation of catalyst nanoparticles as a function of their position on the carbon support. In this work, automated data acquisition software, paired with a custom Python code, have been used to study the effect of different accelerated stress tests (ASTs) on nanoparticle coarsening [2]. Figure 1 shows high-angle annular dark-field (HAADF)-STEM images and EDS maps comparing the cathodes of membrane electrode assemblies (MEAs) following an electrocatalyst AST performed under H2/N2 with that of the heavy-duty AST performed under H2/air. We will discuss how AST conditions affect considerably the spatial distribution of the nanoparticles across the electrode between the membrane and microporous layer. Although the median particle size increased more in the MEA aged under the heavy-duty AST, as determined using a high-throughput image analysis, the quantitative EDS measurements demonstrate that the electrocatalyst AST resulted in more Pt and Co dissolution from the cathode, which is another important indicator of electrocatalyst degradation. We will further present the impact of the relative humidity (% RH) on the degradation mechanisms demonstrated using the same approach. Electron tomography has been used to distinguish the Pt nanoparticles residing on the carbon support surface (exterior) from those within the pore structure (interior) in order to determine the relative stability of interior and exterior nanoparticles. As shown in Figure 2, we will compare the Pt catalyst particle size at the beginning of test (BOT) and end of test (EOT), and discuss the importance of automating the electron tomography workflow, i.e. acquisition, reconstruction, and visualization, to increase sampling and determine the standard deviation of these measurement. The outlook for utilizing low-dose cryo-tomography for limiting damage to the catalyst, support, and especially proton-conducting ionomer will also be discussed [3].

Amichi, Lynda↗
DOE OSTI2022

Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first ~3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm -3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s -1 at 1234 days, the recently detected X-ray emission represents a ≥ 3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥ 3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥ 3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2–25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

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

AGC-4 Experiment Irradiation Monitoring Data Qualification Final Report

The Graphite Technology Development Program ran a series of six experiments to quantify the effects of irradiation on nuclear grade graphite. This report focuses on the fourth experiment, Advanced Graphite Creep 4 (AGC 4). The Advanced Reactor Development (ARD) Technology Development Office (TDO) Program for research and development activities require documentation of qualified monitoring data to design and license the first high-temperature reactor nuclear plant. Qualified data meets the requirements for use as described in the experiment planning and quality assurance documents. Failed data do not meet the requirements and provide no useable information. Trend data may not meet all requirements but still provide some useable information. Use of Trend data requires assessment of how any deficiencies affect a particular use of the data. AGC-4 began with Advanced Test Reactor (ATR) Cycle 157D on May 30, 2015. After irradiating the graphite for two cycles, the capsule was removed from the reactor after ATR Cycle 158A, which ended on January 2, 2016, due to interference with another experiment. Irradiation was resumed with Cycle 162A on October 7, 2017 after the interfering experiment was removed from the reactor. AGC-4 irradiation completed by the end of Cycle 166B on January 10, 2020. Between Cycles 162A and 166B, AGC-4 capsule was removed from the reactor core during two Powered Axial Locator Mechanism (PALM) Cycles 163A and 165A to avoid overheating. All thermocouples (TCs) have functioned throughout the AGC-4 experiment. A total of 9,256,184 out of all 11,167,213 temperature records (or 82.9%) are Qualified for use by the ARD Program and 1,911,029 missing values are Failed. Argon, helium, and total gas flow data were within expected ranges, except only eight out-of-range values occurred during outages. A total of 21,240,627 out of all 22,288,5103 gas flow rates (or 92.8%) are Qualified for use by the ARD Program and 1,644,476 are Failed records mostly due to missing values. Discharge gas line moisture values were consistently low during full ATR power, except for 162B, when moisture content was increasing to more than 200 ppmv by the end of this cycle. During outages, moisture content reached as high as 700 ppmv. 1,089,941 out of a total of 1,107,331 moisture values (or 98.4%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 17,390 missing moisture values are Failed. Graphite creep specimens were subjected to one of three loads: 393, 491, or 589 lbf. For a brief period during Cycle 157D between 12:19 on June 2, 2015, and 08:23 on June 11, 2015, the load cells were wired incorrectly, resulting in missing stack load data. Missing stack loads were estimated from measured ram pressures using regression equations developed from the existing data from Cycle 157D. Estimated stack loads during this period are considered to be an accurate representation of actual load applied to the stacks. These loads deviate slightly from the planned loads. This deviation does not prevent the data from being Qualified for use but must be taken into account when analyzing the effect of load on creep. 6,095,728 out of a total of 6,403,985 moisture values (or 95.2%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 308,257 missing load values are Failed. Stack displacement increased consistently throughout the eight cycles, with total displacement reached highest value of 2.4 in by the end of irradiation. During ATR outages, a set of pneumatic rams raised the stacks of graphite creep specimens to ensure the specimens were not stuck within the test train. This stack raising was performed seven times throughout irradiation. All stacks were raised successfully each time. 4,744,974 out of a total of 6,094,513 displacement values (or 77.9%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 1,349,539 displacement values are Failed mostly due to missing values. Analyses were conducted on correlations between TCs to look for trends and step changes that might indicate instrument degradation or failure. Correlation analysis was used to identify instances when TCs form short circuits, referred to as virtual junctions, which result in TCs reporting temperatures from some location in the capsule other than the location where they were intended to read. No evidence of virtual junctions was found. Analyses were also conducted on control charts of temperature differences of between two TCs, which are expected to behave consistently throughout the entire irradiation period. Upward or downward trend over time indicates at least one TC in the pair was drifted. Examining control charts for all 66 possible pairs out of twelve TCs installed in the AGC-4 capsules reveal no clear drift failures occurred, except unstable behavior of two TCs in Zone 3, TC-7 and TC-8, over irradiation time. In conclus

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

High-order diamond differencing schemes for the Boltzmann Fokker-Planck equation in 3D Cartesian geometries

The Boltzmann Fokker-Planck, an approximate form of the linear Boltzmann equation is commonly used to treat efficiently the transport of charged particles in matter. This paper introduces the application of high-order diamond differencing schemes (HODD), specifically the DD1 and DD2 schemes which are 4- and 6-order accurate respectively, to handle the spatial discretization of that equation in 3D Cartesian geometries. The energy deposition solutions for the coupled transport of electrons and photons presented in this work shows that HODD, compared to classical DD scheme, provides correction to the oscillations and a reduced propensity to yield negative fluxes. They are useful tools to minimize local error, notably in regions with abrupt variations of the flux solution. They also can be used to reduced execution time by decreasing the needed number of voxels to obtain a fixed accuracy. On the tested benchmarks, the DD1 scheme is 87%- 92%-91% more accurate than the classical DD scheme for total, mean per-voxels and maximum deviation of energy deposition values respectively. For comparison, a calculation with 8 times more voxels, requiring roughly 2.5 times more time to execute, is 92%-90%-77% more accurate. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Extending deterministic transport capabilities for very-high and ultra-high energy electron beams

Focused Very-High Energy Electron (VHEE, 50–300 MeV) and Ultra-High Energy Electron (UHEE, > 300 MeV) beams can accurately target both large and deeply seated human tumors with high sparing properties, while avoiding the spatial requirements and cost of proton and heavy ion facilities. Advanced testing phases are underway at the CLEAR facilities at CERN (Switzerland), NLCTA at Stanford (USA), and SPARC at INFN (Italy), aiming to accelerate the transition to clinical application. Currently, Monte Carlo (MC) transport is the sole paradigm supporting preclinical trials and imminent clinical deployment. In this paper, we propose an alternative: the first extension of the nuclear-reactor deterministic chain Njoy-Dragon for VHEE and UHEE applications. We have extended the Boltzmann-Fokker-Planck (BFP) multigroup formalism and validated it using standard radio-oncology benchmarks, complex assemblies with a wide range of atomic numbers, and comprehensive irradiation of the entire periodic table. We report that 99% of water voxels exhibit a BFP-MC deviation below 2% for electron energies under 1.5 GeV. Additionally, we demonstrate that at least 97% of voxels of bone, lung, adipose tissue, muscle, soft tissue, tumor, steel, and aluminum meet the same criterion between 50 MeV and 1.5 GeV. For water, the thorax, and the breast intra-operative benchmark, typical average BFP-MC deviations of 0.3% and 0.4% were observed at 300 MeV and 1 GeV, respectively. By irradiating the entire periodic table, we observed similar performance between lithium (Z = 3) and cerium (Z = 58). Deficiencies observed between praseodymium (Z = 59) and einsteinium (Z = 99) have been reported, analyzed, and quantified, offering critical insights for the ongoing development of the Evaluated Nuclear Data File mode in NJOY.

62 RADIOLOGY AND NUCLEAR MEDICINE↗