Search NASA⌕ Search

SEARCH · Search NASA

Results for “activity coefficient model”

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.

286 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2023

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗
DOE OSTI2022

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management facility (EMF) stack monitor location was qualified using the LV-S1 scale model stack as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale models and its sampling locations were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the EMF stack was performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on the LV-S1 scale model test DV values, the corresponding stack flow rates for the EMF stack are as listed in Table S1. Table S1. Effluent Management Facility Stack Qualified Flow Range. Stack Parameter EM-1 Minimum Qualified Stack Flow (scfm) 781 Maximum Qualified Stack Flow (scfm) 53,432 The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were less than 5°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, which were all within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes each combination of ACV fan with the DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

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

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI1995

Nuclear Safety [Vol. 36, No. 1, January-June 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 1 The Chornobyl Accident Revisited, Part II: The State of the Nuclear Fuel Located Within the Chornobyl Sarcophagus, A A. Borovoi and A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 33 Nuclear Power Safety in Central and Eastern Europe, R. Wilson; 46 Safety of Nuclear Power Reactors in the Former Eastern European Countries, S. Chakraborty; 53 Technical Note: On the Definition of Common-Cause Failures, H. Paula; ACCIDENT ANALYSIS: 58 Modeling and Analysis of Core-Debris Recriticality During Hypothetical Severe Accidents in the Advanced Neutron Source Reactor, S.-H. Kim, V. Georgevich, D. B. Simpson, C. O. Slater, and R. P. Taleyarkhan; 68 Ignitability of Hydrogen/Oxygen/Diluent Mixtures in the Presence of Hot Surfaces, R. K. Kumar and G. W. Koroll; 94 Coupled RELAP5 and CONTAIN Accident Analysis Using PVM, K. A. Smith, A. J. Baratta, and G. E. Robinson; CONTROL AND INSTRUMENTATION: 109 Application of Fuzzy Logic in Nuclear Reactor Control Part I: An Assessment of State-of-the-Art, A. S. Heger, N. K. Alang-Rashid, and M. Jamshidi; DESIGN FEATURES: 122 Twenty-Third DOE/NRC Nuclear Air-Cleaning and Treatment Conference, R. R. Bellamy, J. J. Hayes, and M. W. First; ENVIRONMENTAL EFFECTS: 135 Atmospheric Dispersion and the Radiological Consequences of Normal Airborne Effluents from a Nuclear Power Plant, D. Fang, C. Z. Sun, and L. Yang; 142 Calculation of Distribution Coefficients for Radionuclides in Soils and Sediments, I. Puigdomenech and U. Bergstrom: OPERATING EXPERIENCES: 155 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 158 Operating Experience Feedback Report—Reliability of Safety-Related Steam Turbine-Driven Standby Pumps Used in U.S. Commercial Nuclear Power Plants, J. R. Boardman; 166 Turbine Building Hazards, H. L Ornstein; RECENT DEVELOPMENTS: 169 Reports, Standards, and Safety Guides, D. S. Queener; 175 Proposed Rule Changes as of Dec. 31,1994; ANNOUNCEMENTS: 32 Harvard School of Public Health In-Place Filter Testing Workshop; 134 International Conference on Advances in the Operational Safety of Nuclear Power Plants; 193 30th Tennessee Industries Week; 193 DOE Technical Standards Program 1995 Workshop; 194 Multiphase Flow Experiments and Instrumentation; 180 The Authors; 185 Indexes to Nuclear Safety, Volumes 34 and 35.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Multiphase Reacting Flow Simulations and Optimization of Commercial-Scale Aerobic Bioreactors

Aerobic bioreactors are widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries, wherein microbial bioreactions aid in the conversion of sugars to macromolecules. These reactors are also being actively investigated as a cost-effective pathway for the production of low-value commodities such as biofuels and animal feed. Reactor scale-up is one of the important challenges when designing these systems. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of commercial-scale (~ 500 m3) aerobic bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for oxygen uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM, to perform these simulations, which has been validated against small-scale reactors in our previous work. This work examines the performance of three different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor oxygen mass-transfer coefficient, gas hold-up, and oxygen distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented.

09 BIOMASS FUELS↗
DOE OSTI2022

Hanford Waste Treatment Plant Low Activity Waste Facility Stack Effluent Monitoring - Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant low activity waste (LAW) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAW stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on these DV values, the corresponding stack flow rates for each of the LAW stacks are 815–55,758 scfm for LV-S1, 980–112,078 scfm for LV-S2, 264–22,901 scfm for LV-S3, and 981–79,832 scfm for LV-C2. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack tests at the LAW facility. Flow angle results were primarily less than 10°, except for one LV-S2 Fan A result, which was 13.2°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 1.5 COV and 9.2% COV, which were all within the range of the target % COV values from the scale model tests. Based on these stack verification test results, the four LAW filtered exhaust stack sampling locations meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes single-fan operating conditions for LV-S1 and LV-S2, dual-fan operations for LV-S3 at both the continuous air monitor and record sampler locations, and both the single-fan as well as the dual-fan operations for LV-C2. Further changes to the system configuration or operating conditions that are outside the bounds described in this report may require additional tests or analyses to determine compliance with the standard.

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

Energy density functional and sensitivity of energies of giant resonances to bulk nuclear matter properties

The development of a modern and more realistic nuclear energy density functional (EDF) for accurate predictions of properties of nuclei is the subject of enhanced activity, since it is very important for the study of properties of nuclear matter (NM), giant resonances and, in particular, of properties of rare nuclei with unusual neutron-to-proton ratios. Here, we provide a short review of the current status of the nuclear EDF and the theoretical results obtained for properties of nuclei and nuclear matter. We will first describe a method for determining the parameters of the EDF, associated with the Skyrme type effective interaction, by carrying out a Hartree-Fock based fit to extensive set of data of ground state properties and constraints. Next we will describe the fully self-consistent Hartree-Fock (HF) based random-phase-approximation (RPA) theory for calculating the strength functions S(E) and centroid energies E CEN of giant resonances and the folding model (FM) distorted wave Вorn approximation (DWBA) to calculate the excitation cross section of giant resonances by α scattering. Then we will provide results for: (i) the Skyrme parameters of the KDE0v1 EDF; (ii) consequences of violation of self-consistently in HF-based RPA; (iii) FM-DWBA calculation of excitation cross section; (iv) values of the E CEN of isoscalar and isovector giant resonances of multipolarities L=0–3 for a wide range of spherical nuclei, using 33 EDFs associated with standard form of the Skyrme type interactions, commonly employed in the literature; and (v) the sensitivities E CEN of the giant resonances to bulk properties of NM. We also determine constraints on NM properties, such as the incompressibility coefficient and effective mass, by comparing with experimental data on E CEN of giant resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

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

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

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↗
DOE OSTI2023

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management Facility (EMF) stack monitor location was qualified using a combination of scale model stacks to mitigate the risk of identifying that the sampling location does not meet the qualification criteria on the full-scale stack. The LV-S1 scale model stack was used as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling locations were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. The LV-S1, LB-S1, and LV-S2 scale model stack tests have met the criteria of the ANSI/HPS N13.1-1999 standard to demonstrate the stack sampling locations are well mixed. Verification tests of the EMF stack were performed at normal operating conditions. The minimum 1/6 DV value and the maximum 6 DV value from the scale model testing determine the range of stack flow rates for which the full-scale stack may be operated while remaining in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rate along with the average and maximum qualified stack flow rate based on the LV-S1 scale model DV values. The operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were <5°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, all of which were within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all planned fan operating configurations. This includes each combination of ACV fans with DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

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

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Power balance and divertor asymmetries in the Super-X divertors of MAST-U using SOLPS-ITER

Spherical tokamaks (STs) present unique challenges and opportunities in the area of particle and power exhaust, intensified due to their more compact sizes. Substantial efforts are underway in STs to determine the limits in dissipative operational regimes and advanced divertor solutions, including at MAST-U which provides access to the Super-X divertor configuration. Power balance, and upper/lower divertor asymmetries have been studied using SOLPS-ITER simulations of the MAST-U Super-X divertor. A set of simulations with experimentally inferred transport coefficients with E x B and diamagnetic drifts activated, consisting of density and power scans, and high field side vs low field side gas puff locations, have been used for code experimentation to uncover trends beyond the current experimental parameter space. The upper biased asymmetry (U:L > 1) of the ratio of the peaks of the plasma energy flux densities at the outer targets increases with heating power and decreases with gas puff strength, going from symmetric to up to a factor of 15. The upper target electron temperature has been found to be a good ordering quantity for the magnitude of this asymmetry for all heating powers, gas puff strength, and gas puff locations. The lower divertor biased asymmetry (U:L < 1) of the radiation patterns processed through SOLPS-based bolometry synthetic diagnostics is in qualitative agreement with resistive bolometry experimental results, and it is in quantitative agreement with the trend of the total volume radiation within the divertors of SOLPS. However, radiation measurements alone are not sufficient to infer the magnitude of the asymmetry of the peaks of the power loads at the targets.

MAST-U↗
DOE OSTI2020

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

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

Kinetic and temperature sensitivity properties of soil exoenzymes through the soil profile down to one-meter depth at a temperate coniferous forest (Blodgett, CA)

This dataset contains data on kinetic and temperature sensitivity parameters of the exoenzymes β-glucosidase (BG), leucine/leucyl aminopeptidase (LAP) and acid phosphatase (AP) across a soil profile down to 90 cm depth at Blodgett forest, CA, USA. These data were generated to determine if kinetic and thermal properties of microbial exoenzymes involved in organic matter decomposition varied with soil depth, following variation in soil properties and microbial communities, as part of the study: Alves et al. (2021). Kinetic Properties of Microbial Exoenzymes Vary with Soil Depth but Have Similar Temperature Sensitivities Through the Soil Profile. Frontiers in Microbiology 12:3618. https://doi.org/10.3389/fmicb.2021.735282. This research was performed within the framework of the TES Belowground Biogeochemistry SFA project, in particular association with a long-term field warming experiment of the whole soil profile at Blodgett forest. Samples for this work were collected from locations representative of the field experimental plots. Potential enzyme activity rates were measured using laboratory fluorometric assays with soils collected at 0-10, 10-20, 30-40, 50-60, 60-70 and 80-90 cm deep in biological triplicates (i.e., three soil cores collected at different representative locations). Assays with each soil were conducted over a gradient of eight substrate concentrations per enzyme, and incubated at 4, 10, 16, 25, 35 or 50°C. Enzyme Michaelis-Menten kinetics were modeled over the eight substrate concentrations at each temperature, and the temperature sensitivity of the kinetic parameters was modeled over the six temperatures using linear Arrhenius/Q10 and non-linear Macromolecular Rate Theory (MMRT) models. The dataset includes the fully processed enzyme activity rate data used to model Michaelis-Menten kinetics, calculated kinetic and temperature sensitivity parameters, and basic soil and microbial biomass chemistry for each sample. All data is provided for each individual biological replicate, and kinetic and temperature sensitivity parameters are also provided as means of the biological replicates. The dataset also includes all raw measurement data and code used to parse, combine and perform the analyses described by Alves et al. (2021). For file descriptions, see the file-level metadata files: “enzymes_dataProcessed_flmd.xlsx” (processed data); “enzymes_dataRaw_flmd.xlsx” (compressed raw data and metadata); and “enzymes_code_flmd.xlsx” (compressed code). The experimental design, list of parameters measured, soil and microbial biomass chemistry data, and means of biological replicates for kinetic and temperature sensitivity parameters are also provided as human-readable tables in file “enzymeTraits_design_results_tables.xlsx”.

54 ENVIRONMENTAL SCIENCES↗