Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VI”

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 613 records · Page 34

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗

Nuclear Software Validation for 7% Enriched UO 2 Fuel Lattices

This analysis performs validation for the Virtual Environment for Reaction Application’s (VERA’s) neutronics codes MPACT and Shift against measured data from the Seven Percent Critical Experiment (7uPCX) at Sandia National Laboratory. The benchmarking supports the future application of VERA for analysis of high-assay low-enriched uranium (HALEU) fuel in commercial pressurized water reactors (PWRs). VERA demonstrates very good agreement with measured data, and the deterministic code MPACT also agrees very well with higher-fidelity stochastic methods KENO-VI and Shift in both total neutron flux and fission rate distribution comparisons. Despite some significant differences between the 7uPCX and existing commercial PWRs that challenge the neutron cross section data and transport solver methods used by MPACT, all of the results are excellent. The agreement increases the confidence in the applicability and accuracy of VERA for reactor analysis of HALEU fuels in commercial reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Application of SCALE to Molten Salt Fueled Reactor Physics in Support of Severe Accident Analyses

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a molten salt-fueled reactor (MSR). SCALE capabilities for the modeling of MSR physics were demonstrated based on the Molten Salt Reactor Experiment (MSRE). Of primary interest were the determination of the system’s nuclide inventory, as well as the inventories in the various regions of the loop, considering that the fuel is continuously pumped through the system. This report contains discussions on the following: 1. Determination of the system-average fuel salt inventory considering fission gas removal in the off-gas system and noble metal removal through plating out at the heat exchanger using recent enhancements in SCALE’s depletion sequence TRITON, 2. Assessment of the nuclide spatial distribution throughout the loop using SCALE’s depletion solver ORIGEN, 3. Calculation of the core’s power profile, flux profile, temperature reactivity coefficients, and xenon reactivity using full-core calculations with SCALE’s Monte Carlo code KENO-VI. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the system, as well as the inventory and decay heat of individual regions in the loop, a zone-wise power profile, temperature feedback coefficients, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SCALE 6.2.4 Validation: Nuclear Criticality Safety

The computational bias of criticality safety computer codes must be established through the validation of the codes to critical experiments. A large collection of suitable experiments has been vetted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP) and made available in the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook). More than 600 cases from this handbook have been prepared and reviewed within the Verified, Archived Library of Inputs and Data (VALID), which is maintained by the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory. The performance of the KENO V.a and KENO-VI Monte Carlo codes within the SCALE 6.2.4 code system is assessed using the VALID models of benchmark experiments. A range of nuclear cross section libraries based on Evaluated Nuclear Data File (ENDF)/B-VII.1 in both multigroup (MG) and continuous energy (CE) formats is considered. The critical experiments available to validate the KENO V.a code cover 15 broad categories of systems. These systems use a range of fissile materials, including a range of uranium enrichments, various plutonium isotopic vectors, and some mixed uranium/plutonium oxides. The physical forms of the fissile material also vary and are represented as metal, solutions, or arrays of rods or plates in a water moderator. The neutron energy spectra of the systems also vary and cover fast, intermediate, mixed, and thermal spectra. Over 550 of the total cases use the KENO V.a code for the four nuclear data libraries considered in this report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modular HF Isolated MV String Inverters Enable a New Paradigm for Large PV Farms

The “Modular HF Isolated MV String Inverters Enable a New Paradigm for Large PV Farms” project focuses on exploring alternative power converter and system-level plant configurations to achieve the lowest cost and highest energy output for a given solar plus storage (e.g., PV + battery) plant, including the use of medium voltage (MV) collection while taking into account detailed models of all elements. To realize this objective, four approaches were utilized (i) employ a novel Medium Voltage String Inverter (MVSI) topology (soft switching solid state transformer – S4T) to convert 1000 Vdc to 4.16 kVac; (ii) plant collection using standard, low-cost overhead MV distribution network; (iii) enable energy storage integration without additional converter cost to achieve dispatchability of the PV resource; and (iv) provide advanced functionality (autonomous operation, track ISO signals for dynamic balancing and ancillary services, and PV farm operation as a virtual grid resource). Subsequently and in alignment with the previously mentioned approaches, the project was structured in five efforts (i) S4T MVSI simulation and design; (ii) system analysis and storage optimization; (iii) financial analysis; (iv) power converter prototype build and test; and (v) regulatory and commercial impact study. The outcomes provided by each effort can be summarized as follows (i) Project explored the use of MV AC distribution architecture for hybrid PV+storage utility-scale PV farms; (ii) Detailed loss and LCOE analysis for AC and DC side BESS architecture, including multiple converter topologies, as well as for proposed MVSI/MDCT systems; (iii) MVSI was built and holds promise but needs lower-cost high-voltage Si-C devices, which does not seem possible in the near term; (iv) MDCT provides a simpler modular building block – validated through HIL and farm level modeling, simulation and experimental validation; (v) 300 kVA MDCT prototype built and tested, technology is being commercialized; and (vi) Regulatory model of utility building PV plants, where PV panels are treated as DC generation (IPP), seems viable and can allow improved grid integration.

14 SOLAR ENERGY↗

NDCS Group Meeting Work overview: NCSP Task ORNL IP&D-5 [Slides]

The presentation discusses the primary goal of the NCSP Task IP&D-5 which is to use available data from HPRR operation to create a benchmark report for inclusion in the ICSBEP, as a Criticality Accident Alarm System (CAAS) shielding benchmark. The presentation elaborates on the HPRR operation and how it relates to the KENO VI model, sulfur pellets measurements, and calculation methodology between KENO and MAVRIC FW-CADIS. The presentation concludes by stating that some work is still needed to evaluate uncertainties and to explain the discrepancy computation/experiment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Novel Chalcopyrites for Advanced Photoelectrochemical Water Splitting

With the support of DoE’s EERE office, our team has established a unique tool-chest of capabilities, including theoretical modeling (Lawrence Livermore National Laboratory: LLNL), state-of-the-art synthesis (Hawaii Natural Energy Institute: HNEI, Stanford, and the National Renewable Energy Laboratory: NREL) and advanced materials and interfaces characterization (University of Nevada, Las Vegas: UNLV, and Lawrence Berkeley National Laboratory: LBNL), to accelerate the development of high efficiency and durable chalcopyrite materials for advanced photoelectrochemical (PEC) water splitting. Using this synergistic approach, we have successfully created new wide bandgap chalcopyrite photocathodes generating over 10 mA/cm 2 , developed innovative strategies to protect them from corrosion, and engineered novel integration methods to circumvent thin film materials mechanical, chemical and thermal incompatibility. In Task 1 “Modeling and synthesis of chalcopyrite photocathodes”, we expanded our library of wide bandgap chalcopyrites for PEC water splitting. With support from LLNL’s “Computational Materials Diagnostics and Optimization of PEC Devices”, LBNL’s “photophysical” and NREL’s “I-III-VI Compound Semiconductors for Water-Splitting” nodes, we investigated two new chalcopyrite candidates for PEC water splitting: Cu(In,Al)Se 2 and Cu(In,B)Se 2 . We also further developed ordered vacancy compounds, such as CuGa 3 Se 5 , with unprecedented durability during PEC waters splitting in acidic solutions. In Task 2 “Interfaces engineering for enhanced efficiency and durability”, we addressed both the non-ideal band-edge positions of chalcopyrites with respect to water redox potentials, as well as their chemical instability under PEC water splitting, with a buried-junctions approach. With help from NREL’s “High-Throughput Experimental Thin Film Combinatorial Capabilities” and “Corrosion Analysis of Materials” nodes, we engineered environmentally friendly n-type buffers, including Mn x Zn 1-x O, to adjust the chalcopyrite band-edge positions and achieved photovoltages as high as 925 mV. Also, we integrated non-precious catalytic-protecting layers, such as WO 3 , to enhance the water splitting long-term stability of chalcopyrite absorbers. Finally, in Task 3 “Hybrid photoelectrode device integration”, we proposed an innovative method to bond wide bandgap photocathodes onto narrow bandgap PV drivers at room temperature using conductive polymers. Our semi-monolithic approach addressed fundamental processing incompatibility issues, as both the photocathode and the PV driver are processed separately. Proof-of-concept whole-chalcopyrite tandems were obtained by consecutive exfoliation and transfer of fully integrated 1.85 eV CuGa 3 Se 5 and 1.13 eV CuInGaSe 2 stacks from their Mo/SLG substrates onto a new single FTO host substrate.

08 HYDROGEN↗

(U) Average Scattering Angle Cosine Sensitivities

The Evaluated Nuclear Data File (ENDF) provides covariances of the average scattering angle cosine $\overline{μ}$ for 117 nuclides, and the NJOY nuclear data processing code can process these covariances. To propagate them to the uncertainty of a response due to uncertainties in $\overline{μ}$ requires sensitivities of the response with respect to $\overline{μ}$. This report uses NJOY’s definition of $\overline{μ}$ and derives a relationship between the sensitivity of a response to $\overline{μ}$ and the sensitivity of the response to the anisotropic P 1 elastic scattering cross section. This relationship has been implemented in the multigroup neutron sensitivity code SENSMG, which relies on the discrete ordinates code PARTISN for the neutron transport calculations. Kodeli has also calculated sensitivities with respect to $\overline{μ}$ and associated uncertainties. However, the formal derivation of this paper shows that the definition of the sensitivity with respect to $\overline{μ}$ is ambiguous because there are multiple ways to change $\overline{μ}$ by changing scattering cross sections. Here we derive two such ways. Both are implemented in SENSMG. This report is organized as follows. Section II discusses $\overline{μ}$ and its significance. Section III presents the equations used for the sensitivity of a response with respect to $\overline{μ}$, and Section IV discusses how to compute these sensitivities with SENSMG and MCNP6.3. Section V is a fast-spectrum test problem, and Sec. VI is a thermal-spectrum problem. Section VII is a summary. The example problem inputs are listed in the appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Plutonium Oxidation State Distribution under WIPP Relevant Conditions (Rev. 2)

The oxidation state of plutonium in the Waste Isolation Pilot Plant (WIPP) environment has been a topic of interest since the initial compliance certification application. Plutonium (Pu) was initially expected to be present primarily as Pu(IV), but since the presence of Pu(III) could not be ruled out, it was also included in performance assessment (PA) calculations. The redox of the other variant actinides, uranium (U) and neptunium (Np), are also coupled to the plutonium redox. This led to a position in performance assessment calculations that is commonly referred to as “50/50”. In which half of the PA realizations assume solubility is dominated by the low oxidation state [U(IV), Np(IV), Pu(III)] and the other half assume solubility is dominated by the higher oxidation state [U(VI), Np(V), Pu(IV)]. Recently, a joint group of WIPP project participants from Los Alamos National Laboratory (LANL), Sandia National Laboratories (SNL), and the Department of Energy’s Carlsbad Field Office (DOE-CBFO) agreed on a new path forward for the oxidation state model in PA which will decouple the redox active actinides and will instead consider their redox as a function of E h . This will take the current binary approach and turn it into a system with four possible redox states.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Soliton formation and topology manipulation of coupled spins via ultrafast re-magnetization

The major goal of the project was to explore the properties of magnetic order driven far out of equilibrium by optical excitations. These include fundamental questions related to the interplay of magnetic, structural and electronic degrees of freedom in materials that are optically excited. These problems break down in short term quenching of the magnetization as a results of the energy of the optical pulse studying how the energy of the optical pulse and the angular momentum of the magnetization flow between the different degrees of freedom. This is followed by the longer time re-emergence of the magnetism as the system cools. Within this goal we explored the formation of solitons, both with and without non-trivial topology, via rapid re-magnetization processes after optical-driven ultrafast demagnetization. Guided by theory, we predict that turbulence and modulational instabilities will drive the formation of solitons, including dispersive shock waves, magnon droplets, and skyrmions. In the case of topological defect formation, i.e. skyrmion generation, related processes have long been predicted based on general principles of symmetry-breaking phase transitions; the density of topological defects in a long-range-ordered phase can be controlled by varying the quench rate through the second-order phase transition, i.e. the Kibble-Zurek (KZ) mechanism. We, more broadly, had the goal to study a range of other emergent magnetic behaviors after optical excitation. While great attention has been given to ultrafast demagnetization, less is known of the subsequent spin dynamics and coupling to the lattice. What is increasingly appreciated is that ultrafast demagnetization leads to spin currents that can carry angular momentum from the rapidly demagnetized sample. These spin currents appear in many ways but fundamentally controls demagnetization, drives interactions between regions of the material, and control the spin resulting structure. They can be probed many ways including THz emission which, in turn, gives insight into the demagnetization processes. Simultaneously if there is magneto-elastic coupling ultrafast demagnetization can drive structural excitations that are not expected from thermal energy added to the lattice. We explored a range of complex phenomena that arise from photo-excitation of magnetic systems with the coupling of electronic, magnetic and structural degrees of freedom. We have made a number of fundamental discoveries that are reflected in our publications list with additional work still being prepared for publication. Highlights of this work include (i) ultra-efficient, nonlinear THz surface acoustics (ii) spin-current-mediated rapid magnon localization and coalescence, (iii) spin-wave soliton formation in ferromagnetic FePt nanoparticles, (iv) dynamic phonon coupling in elemental antiferromagnetic Cr, (v) THz emission from Co/Pt bilayers and FeRh/Pt bilayers, (vi) theoretical investigation of spin hydrodynamics, solitons and shock waves, and (vii) ultrafast perturbation of magnetic domains by optical pumping.

36 MATERIALS SCIENCE↗

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗

Development of Advanced Solid Sorbents for Direct Air Capture (Final Report)

RTI International is partnering with Creare, and Mohammed VI Polytechnic University (UM6P) successfully met all major technical objectives as outlined below: Identified one MOF adsorbent for DAC; Identified one amine-P-dendrimer adsorbent for DAC; Performed Computational fluid dynamics simulations of the MOF and amine-P-dendrimer adsorbents and validate them with experimental data; Selected one adsorbent for DAC; Demonstrated the scale-up of selected candidate and perform cost review; Performed a preliminary process design; Performed technology maturation plan and environmental health and safety risk assessment. After a thorough comparison between both adsorbents (MOFs and P-dendrimers) and based on their technical merits, the amine-based P-dendrimers sorbent was selected as the best performing adsorbent for DAC. In fact, the amine-based P-dendrimers have demonstrated superior performance over the MOF sorbents under relevant DAC conditions. The selected sorbent has been used in subsequent projects (DE-FE0032099 and DE-AR0001412) and has demonstrated similar superior performance with no loss of activity after long-term sorption-desorption testing (e.g., 100 cycles). This sorbent has demonstrated a very high and stable sorption performance, excellent kinetics, and fast regeneration at 80 °C, specifically at 75% relative humidity (RH). The incorporation of RTI’s high-performance, high-durability amine-based P-dendrimers sorbent and Creare’s hybrid additive manufacturing (H-AM) technology that will produce high performance, compact heat and mass exchange structures at low cost using methods, will lead to a contactor that is optimized for wind-driven operation. This contactor technology provides a technically viable pathway for reducing the cost of DAC to <$100/tonne of CO 2 .

36 MATERIALS SCIENCE↗

Surface Complexation/Ion Exchange Hybrid Model for Radionuclide Sorption to Clay Minerals (M4SF-23LL010301062)

This progress report (Level 4 Milestone Number M4SF-23LL010301062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite International Collaborations Activity Number SF-23LL01030106. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. The focus of this milestone is the establishment of international collaborations for surface complexation modeling and the associated impacts of unlocking larger, community-based datasets. More specifically, we are developing a database framework for Spent Fuel and Waste and Science Technology (SFWST) that is aligned with the Helmholtz Zentrum Dresden Rossendorf (HZDR) sorption database development group in support of the database needs of the SFWST program. In our FY22 effort, we described a detailed analysis of U(VI) sorption to quartz through both traditional surface complexation modeling and through a hybrid ML framework. In FY23, effort was placed on publication of these results and expansion of the LLNL surface complexation and ion exchange database (L-SCIE) in order to assess mineral-based radionuclide retardation under a wider variety of geochemical conditions (e.g., ionic strength, varying electrolyte compositions). Efforts were initiated to expand L-SCIE to include radionuclide surface complexation and ion exchange to clays that are relevant to subsurface geochemical processes occurring at nuclear waste repositories. In particular, a large source of sorption data for clays resides at the Paul Scherrer Institute (PSI) (work primarily by Bradbury and Baeyens) and we initiated discussions on how to retrieve those data and apply FAIR principles to those datasets. In addition to L-SCIE development, two hybrid models that incorporate AI/ML were investigated and compared to discern the most promising approaches for accurate and precise estimations of radionuclide retardation. Key considerations for future model development include (1) the ability to reduce computational burden on determining retardation coefficients for PA and (2) the ability to quantify and predict radionuclide-mineral partitioning at a more efficient, rapid pace due to automated workflows. Upon the careful consideration of the most effective modeling approaches, we are identifying ways to implement these approaches into PA. Ultimately, the data science-based workflows will provide a major incentive for other institutions to adopt a FAIR-formatted, interoperable database. LLNL will play a key role in disseminating sorption data and acting as good data stewards by updating the database in a consistent format and assessing the quality of the newly assimilated data in an organized fashion. To this end, all data and workflows are open access and made available on the LLNL Seaborg research website (https://seaborg.llnl.gov/resources/geochemical-databases-modeling-codes).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Radionuclide Interaction with Hydrothermally Altered Repository Materials (M4SF-23LL010302052)

This progress report (Level 4 Milestone Number M4SF-23LL010302052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Activity Number SF-23LL01030205. The research is focused on actinide and radionuclide sequestration in hydrothermally altered repository materials. In FY23, a manuscript was in preparation for publication summarizing our analysis of radionuclide sorption and coprecipitation into Fe oxide phases and evaluation of radionuclide partitioning values across a range of radionuclides relevant to performance assessment. We demonstrated our approach in detail using Se sorption and coprecipitation with iron oxide minerals. These data were presented in our FY22 annual report and will not be repeated here. We also initiated experiments to identify radionuclide interaction with hydrothermally altered crystalline repository and backfill materials. Recent research performed at Los Alamos National Laboratory (LANL) and Sandia National Laboratory (SNL) has provided key insights regarding the hydrothermal alteration behavior of bentonite backfill in the presence of repository materials (steel, concrete, etc.). We are now examining how mineral alteration affects retardation behavior of plutonium and a suite of other radionuclides. These experiments also allow us to test the predictive ability of our component additivity approach to surface complexation and ion exchange. Our guiding hypothesis is that a robust surface complexation/ion exchange model and associated database, developed using our L-SCIE approach, can effectively predict changes in radionuclide sorption behavior resulting from the hydrothermal alteration of mineralogy in a repository near field. A short update of results to date is presented below. A manuscript was also in preparation describing a self-consistent model and approach to simulating Se(IV) and Se(VI) sorption to 5 iron oxide phases based on our L-SCIE community database. This will be the first implementation of a multi-mineral and multi-oxidation state sorption model using our new L-SCIE database and workflow. A short summary of these results is presented below.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

TRACER Carbonaceous Aerosols Thrust – University of California, Davis Field Campaign Report

The broader U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign aims to increase our understanding of convective cloud life cycles and aerosol-convection interactions. Our TRACER Carbonaceous Aerosols Thrust-University of California, Davis (TRACER-CAT-UCDavis) study complemented these broader aims by characterizing and quantifying the optical properties and composition of carbonaceous aerosols during part of the TRACER intensive sampling period (July 1- July 31, 2022) at the first ARM Mobile Facility (AMF1) main site (M1) in La Porte, Texas. Our measurements complemented the suite of instrumentation already provided by the AMF1, expanding the capabilities through deployment of unique, state-of-the-science instrumentation. The instrumentation included: (i) two cavity-attenuated phase shift spectroscopy single-scatter albedo (CAPS-SSA) instruments operating at 530 nm and 630 nm, and that were modified to characterize particle light absorption, extinction, and scattering at elevated humidities; (ii) the UC Davis dual-wavelength cavity ringdown-photoacoustic spectrometer (CRD-PAS), which characterizes dry particle extinction and absorption at 405 nm and 532 nm; (iii) a soot particle aerosol mass spectrometer (SP-AMS) that operated in “laser only” mode that characterized the size-dependent compositions of black carbon (BC)-containing particles; (v) a thermal denuder, to remove coatings on particles; and (vi) a scanning electrical mobility sizer (SEMS), to characterize particle mobility diameters from 10-1300 nm. Our measurements occurred alongside complementary observations made by Los Alamos National Laboratory (LANL) during the TRACER-CAT-LANL study, including a humidified CAPS-SSA instrument operating at 450 nm. Our primary scientific interest is in understanding the relationship between particle composition and light absorption, with a particular focus on the influence of water uptake. While it is known that coatings on BC can enhance absorption, the extent to which this occurs in the atmosphere and the specific role that water plays as a coating remain unclear. The TRACER-CAT-UCDavis measurements were made with near-complete coverage for the CRD-PAS, SP-AMS, and SEMS throughout the intensive period. The UC Davis humidified CAPS-SSA instruments had significant challenges with operation owing to the demanding conditions (large temperature fluctuations, high humidity), exacerbated by supply chain issues that delayed resolution of these challenges. However, the humidified CAPS-SSA instrument operated by LANL operated throughout the intensive period with near-complete coverage. The dry light extinction measurements from the CRD-PAS measurements and the LANL CAPS-SSA exhibited a good correlation, although the CAPS-SSA systematically measured greater extinction values than expected. While all instruments, with the exception of an aerodynamic particle sizer (APS), measured behind a common particulate matter (PM)2.5 μm cyclone, the greater extinction measured by the LANL CAPS-SSA compared to the CRD-PAS may have resulted from different losses of larger particles in the sampling lines from the cyclone to the instruments; the tubing length was shorter from the cyclone to the LANL CAPS-SSA, consistent with this idea. A summary of the TRACER-CAT-UCDavis measurements, along with some of the TRACER-CAT-LANL measurements, are shown in the figure below. Notably, there were periods when the contributions of presumed dust were substantial and even dominated the observed light extinction and absorption. Also, there was a clear shift in the behavior of submicron particles from before July 16, 2022 to after, with the prior period exhibiting regular episodes of new particle formation and the latter period exhibiting rapid variations in the concentrations of small particles.

54 ENVIRONMENTAL SCIENCES↗

Artificial Intelligence for Natural Gas Utilities: A Primer

Modern natural gas utilities face numerous challenges and competing priorities from various stakeholders. Policymakers, customers, and advocacy groups want to see gas utilities improve performance on safety, reliability, resilience, affordability, and environmental stewardship. State utility regulators — public utility commissions — are responsible for overseeing utility performance, ensuring that ratepayer funds are being spent in the public interest, and aligning utility goals with public goals. The use of new technologies is critical to enabling cost-effective performance on these attributes. Artificial intelligence (AI) is a widely used term among utilities and regulators, but the term means different things to different stakeholders, and it is often used to describe data analytics approaches that fall short of the formal definition of AI, which is: “…the ability of a machine to receive inputs and produce a behavior or reaction similar to that of an intelligent human being.” AI (and related tools, techniques, and technologies) can help utilities solve current and emerging challenges. By combining customer and system data with analytical tools and technologies, AI can augment human decision-makers by assisting in identifying problems and events before they occur, enabling resources to be more efficiently directed across utility infrastructure. The intended primary audience for this primer is state regulators, although utilities and other stakeholders might also find it useful and relevant to improve their awareness of AI. The objectives of this primer are to: (a) offer a set of broadly applicable definitions for AI and related terms, allowing regulators, utilities, and other stakeholders to speak the same language; (b) discuss how AI is currently being implemented in the gas utility sector; and (c) understand the challenges affecting AI solutions and how tools might be implemented in the future. This primer fits within NARUC’s goals of providing impartial information to improve the ability of public utility commissions to regulate in the public interest. As such, this primer does not seek to recommend AI over any other investment, nor does it endorse any particular vendor, product, or approach. It does seek to prepare state regulators to oversee AI investments by sharing information about the current landscape of commercially available tools. To these ends, the primer is organized as follows: Section I discusses the current environment in which natural gas utilities operate and how AI, when thoughtfully designed and implemented, can enable utilities to achieve performance goals; Section II offers definitions of AI and related terms within the data analytics discipline; Section III provides three current opportunities for which AI can offer solutions: replacing aging gas distribution infrastructure, preventing excavator damage to gas distribution infrastructure, and improving energy efficiency programs. This section discusses each problem statement in detail. Second, Section III includes a discussion of how costs and benefits of investments to solve each problem are measured. And third, this section offers real-world examples of utility implementation of AI solutions; Section IV discusses challenges with implementing AI, both from the perspective of utilities and regulators; Section V suggests areas in which AI could feasibly be implemented in the near future; Finally, Section VI offers concluding thoughts and areas for further research.

03 NATURAL GAS↗

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗

Building a Framework to Understand Transition Metals' Behavior in Euxinic Conditions (Final Technical Report)

This project focuses first and foremost on metal sulfide geochemistry and mineralogy as controlled by a complex matrix of environmental factors. The principal investigator’s group aim to illuminate the metal-sulfide reaction mechanisms, rates, and pathways through systematic experimentation and data collection and analyzing the relationships between the characteristics of the produced metal sulfide solid-phase/aqueous complexes and the environmental factors. This understanding is essential for obtaining a full picture of the complex cycling patterns of single or multi metal species in sulfidic environments ranging from deep-see basins, hydrothermal vents, inland seas, terrestrial water bodies, to engineered remediation systems. The specific goal of this past project was to illuminate the reaction mechanisms and kinetics of metal anions and sulfide in mixed metal cation-metal anion-sulfide systems under various aqueous conditions (which resembled a range of naturally occurring euxinic settings). For the period of this contract, we investigated the molybdenum-iron-sulfide system, with an emphasis on the conditions that caused solid phase formation. We focused on quantifying the mobility/sequestration of molybdenum under each experimental condition and identified the changes of valence states for each involved element (i.e., Mo, Fe, and S) in the precipitate. We also proposed pathways for the electron transfer that occurred in aqueous chemistry. The major analytical tools used for this study include UV-visible light spectroscopy, transmission electron microscopy, X-ray photoelectron spectroscopy, and synchrotron-based X-ray absorption spectroscopy (access gained through facility proposals to the Canadian Light Source). The biggest finding of this project was that besides pH, the iron-sulfur chemistry has a dominant control of the thiolation kinetics and subsequent reduction of Mo(VI), which are likely prerequisites of molybdenum sequestration in anoxic conditions. The results have been written up as manuscript by the end of this project (see Phillips et al.). The experimental results of this project may be critical for advancing our understanding of (1) basic chemistry involving transition metals and reduced sulfur species, (2) the validity of certain geochemical proxies, and (3) the stability and evolution of euxinic geochemical environments. It is noted that the basic results obtained through this project also have implications for Mo-S cluster-based catalyst development in inorganic chemistry and materials sciences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗