Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed model predictive control”

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 811 records · Page 45

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO 2 ) and methane (CH 4 ) release from HCP soils are needed to support modeling and evaluation of warming-induced CO 2 and CH 4 fluxes from tundra soils. This study quantifies the distribution of microbial CO 2 and CH 4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at −2, +4, and +8°C were used to quantify temporal trends in CO 2 and CH 4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2–9 days), higher CO 2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO 2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH 4 -producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO 2 and CH 4 production with elevated water content, but HCP trough mineral soils produced more CH 4 than LCP mineral soils. HCP mineral soils also released substantial CO 2 but did not show a strong trend in CO 2 and CH 4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

59 BASIC BIOLOGICAL SCIENCES↗

Bioinformatics resource manager v2.3: an integrated software environment for systems biology with microRNA and cross-species analysis tools

Background: MicroRNAs (miRNAs) are noncoding RNAs that direct post-transcriptional regulation of protein coding genes. Recent studies have shown miRNAs are important for controlling many biological processes, including nervous system development, and are highly conserved across species. Given their importance, computational tools are necessary for analysis, interpretation and integration of high-throughput (HTP) miRNA data in an increasing number of model species. The Bioinformatics Resource Manager (BRM) v2.3 is a software environment for data management, mining, integration and functional annotation of HTP biological data. In this study, we report recent updates to BRM for miRNA data analysis and cross-species comparisons across datasets. Results: BRM v2.3 has the capability to query predicted miRNA targets from multiple databases, retrieve potential regulatory miRNAs for known genes, integrate experimentally derived miRNA and mRNA datasets, perform ortholog mapping across species, and retrieve annotation and cross-reference identifiers for an expanded number of species. Here we use BRM to show that developmental exposure of zebrafish to 30 uM nicotine from 6–48 hours post fertilization (hpf) results in behavioral hyperactivity in larval zebrafish and alteration of putative miRNA gene targets in whole embryos at developmental stages that encompass early neurogenesis. We show typical workflows for using BRM to integrate experimental zebrafish miRNA and mRNA microarray datasets with example retrievals for zebrafish, including pathway annotation and mapping to human ortholog. Functional analysis of differentially regulated (p<0.05) gene targets in BRM indicates that nicotine exposure disrupts genes involved in neurogenesis, possibly through misregulation of nicotine-sensitive miRNAs. Conclusions: BRM provides the ability to mine complex data for identification of candidate miRNAs or pathways that drive phenotypic outcome and, therefore, is a useful hypothesis generation tool for systems biology. The miRNA workflow in BRM allows for efficient processing of multiple miRNA and mRNA datasets in a single software environment with the added capability to interact with public data sources and visual analytic tools for HTP data analysis at a systems level. BRM is developed using Java™ and other open-source technologies for free distribution (http://www.sysbio.org/dataresources/brm.stm).

59 BASIC BIOLOGICAL SCIENCES↗

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Influence of African Easterly Waves on Atlantic Tropical Cyclone Tracks and Landfall in Large Ensembles

African easterly waves (AEWs) are an important precursor or “seed” for Atlantic tropical cyclones (TCs), with 60%–80% of major hurricanes observed to originate from AEWs. However, climate model simulations indicate that AEWs are not necessary to maintain annual Atlantic TC frequency. Furthermore, small ensembles suggest that AEWs may impact the spatial distribution and landfall of Atlantic TCs. Here, we investigated the influence of AEWs on the spatial distribution of Atlantic TC tracks and landfall using 50-member ensembles of TC-permitting regional model simulations for five hurricane seasons characterized by different levels of TC activity. The control simulations are seasonal hindcasts in which AEWs were prescribed through the eastern lateral boundary condition using reanalysis. In the experiments, we suppressed AEWs by applying a 2–10 day filter to the eastern lateral boundary condition. In response to AEW suppression, we discovered statistically significant increases in Atlantic TC frequency (10%–26%) and landfall (16%–44%), a westward shift in TC genesis location and tracks with increased landfall over the Caribbean Islands, southwestern Gulf Coast, and southeastern US coast, and increases in mid-tropospheric relative humidity in the main development region. In addition, we evaluated TC genesis mechanisms in the absence of AEWs and found evidence that the intertropical convergence zone (ITCZ) intensified and extended northward, resulting in ITCZ wave-breaking that shed vortices which served as TC seeds. By uncovering the connections between TC seed types and the likelihood of TC landfall, this research can provide Atlantic coastal and island communities with useful information to prepare for TC impacts.

African easterly waves↗

Synchronous Machine Governor Upgrade

Conventional generation sources play a critical role in the stability and reliability of the electrical grid, particularly as we transition towards more renewable energy sources. To understand and accurately emulate their behavior for optimizing grid operations and ensuring seamless integration with renewable technologies, it is essential to better emulate the grid- and plant-level impacts of conventional generation sources, such as natural gas (NG) driven heat recovery steam generators (HRSGs) and combustion turbines (CTs). Therefore, a governor model is developed in a programmable logic controller (PLC) to investigate the performance of the conventional generator under various dynamic operating conditions and to identify the impact on grid stability in a controlled environment. The governor model aims to enable the hardware-in-the-loop (HIL) based emulation of these conventional generation sources using the existing 2 MVA synchronous machine/generator that is driven by a flexible 2.5 MW variable speed drive. This setup will allow us to replicate the dynamic characteristics and response behaviors of NG-driven HRSGs and CTs. The controls for the emulated conventional plants follow the industry standard and are adjustable, ensuring they accurately reflect the operational capabilities and limitations of real-world systems. These controls include load-following capabilities, ramp rates, startup and shutdown sequences, and emissions characteristics. By incorporating these adjustable controls, we aim to capture the nuanced impacts of conventional generation, such as their ability to provide ancillary services like frequency regulation, voltage support, and spinning reserve. In this report, we simulate two types of dynamic operations: grid-connected and islanding. For each dynamic operation, representative starting sequences are tested, including turbine purge, ignition, speed ramping up, generator excitation and synchronizing, and breaker close. The HIL based tests provides insights for field deployment, specifically the high-fidelity governor model provides results to predict the potential stability and reliability risk and suggest possible integration measures (e.g., generation and load balancing, tuning of governor control parameters). Ultimately, this enhanced emulation capability will be integrated into our Advanced Research on Integrated Energy Systems (ARIES), enabling us to conduct comprehensive studies on the interactions between conventional and renewable energy sources. By better understanding these interactions, we can develop strategies to optimize the overall performance and reliability of the grid. This will support the deployment of advanced grid management techniques, such as demand response, grid-forming inverters, and energy storage systems. The main contributions are summarized as follows: (1) This report introduces a PLC-based governor model for gas turbines. This model accurately simulates the dynamic behavior of conventional generation sources under various operational scenarios; (2) The model is integrated with an HIL testbed that includes a 2.5 MW variable speed drive and a 2 MVA synchronous machine. This setup enables realistic, real-time emulation of conventional power plants, particularly NG driven HRSGs and CTs; (3) The developed model is adaptable to various gas turbine configurations and allows for precise control over parameters such as MW ramp rates. This flexibility makes it a valuable tool for future research and industry collaboration; and (4) By incorporating the model into the National Renewable Energy Laboratory's Advanced Research on Integrated Energy Systems, the report lays the groundwork for future studies on interactions between conventional and renewable energy sources, enhancing the ability to develop advanced grid management strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Peak Reduction Using Mode Adjustment of Heat Pump Water Heaters in a Residential Neighborhood

Building electrification is putting pressure on distribution grid worldwide. Peak reduction is an important concern that can help reduce the growing stress and allow to defer investments in new capacity. Water heaters represent a convenient way of reducing peak because they are less dependent on weather, and their storage volume allows for asynchronous water heating and hot water use. Previous empirical studies investigated the ability of water heaters to reduce peak through the adjustment of the temperature setpoint. However, not all equipment vendors offer this option. This study aims at understanding the feasibility of peak reduction with an alternative configuration available in the market - by adjusting the device mode rather than temperature setpoint. The peak reduction methodology is tested in an occupied 46-townhome neighborhood located in Atlanta, GA. We find that peak shifting is possible with the adjustable mode approach, with the change in the peak load by 30-60%.

demand response↗

Roll-to-Roll Advanced Materials Manufacturing DOE Laboratory Collaboration - FY2018 Final Report

High-value R2R processing is used to support a wide range of products in applications which span many industrial business sectors. The overall R2R process technology can be considered as “mature” as the process methodology and has been in use for decades. This continuous processing technique traditionally involves deposition of material(s) onto moving webs or carriers or other continuous R2R, belt-fed, or conveyor-based processes that enable successive steps to build a final version which serves to support the deposited materials. Current process technologies, which typify “roll to roll”, include tape casting, silk-screen printing, reel-to-reel vacuum deposition/coating and R2R lithography. Products supported by R2R manufacturing include micro-electronics, electro-chromic window films, PVs, fuel cells for energy conversion, battery electrodes for energy storage, and barrier materials. Due to innovation in materials and process equipment, high-quality yet very low-cost multi-layer technologies can be manufactured on a very cost-competitive basis. To move energy-related products from high-cost niche applications to the commercial sector, a means must be available to enable manufacture of these products in a cost-competitive manner that is affordable by the general consumer. Fortunately, products such as fuel cells, thin- and mid-film PVs, batteries, electrochromic and piezoelectric films, water separation membranes, and other energy saving technologies readily lend themselves to manufacture using R2R approaches.Within the DOE EERE AMO, it was recognized that establishment of a program supported at the DOE National Laboratories, along with their immense design of materials and equipment modelling capabilities enabled with use of high-performance computing, could take advantage of available R2R infrastructure to manufacture new technologies. In FY 2016, a R2R Consortium was established and provided with initial “seed” funding to take an approach that was envisioned to be supportable by advanced manufacturing R2R processes. This collaborative approach was designed to foster identification and development of materials and processes related to R2R for clean-energy product development. Using computational and experimental capabilities by acknowledged subject matter experts within the supported National Laboratory system, this collaborative project would leverage the capabilities and expertise at each of four laboratories to further the development of an enabling high-volume cost-competitive platform technology. The collaboration team that is comprised of ORNL, ANL, NREL and LBNL, coordinating with EPB and other selected industry partners, was formed in April 2016 to initially address enhancing battery electrode performance and R2R manufacturing challenges. The research efforts were to predict and measure changes and results in electrode morphology and performance based on process condition changes; to evaluate mixed, active, particle size deposition and drying for novel electrode materials; to model various process condition changes and the resulting morphology and electrode performance; and to develop and validate NDE techniques for in-line measurement of battery electrode material properties. These efforts carried through FY 2017 and completed at the end of FY 2018.The approach was to look at compositions of materials with different particle sizes to make electrode samples using a R2R manufacturing process. The shape, size, and morphology of the materials, the chemistry of the formulation, the nature of slurries, their coating rate, the rate of drying all play a role in determining the final coating architecture, quality, and performance. A commercial cathode material was selected to make a series of cathodes and anodes by single pass, dual pass and slot die methods. Analysis of all the compiled results for this battery electrode development were that the best performing cathodes included: dual-pass electrode with large particles near the foil, mixed small and large particles, and small particles only. Whereas, the best performing anodes were with a mix of small and large graphite particles. An additional core project was added in FY 2017 to conduct studies of fuel cell materials that can be produced using R2R processes. The goal of this project is to explore, understand and optimize material and process parameters to support increased throughput, increased quality, and reduced cost for high volume production of gas-diffusion electrodes (GDEs) for PEMFCs. Project work in FY 2018 were to use a R2R process to fabricate electrodes without ionomer overlayer that can produce the equivalent mass activity as spray-coated electrodes with an ionomer overlayer. Oxygen-limiting current measurements were utilized to optimize oxygen mass transport. Alcohol-rich solvents were investigated, and results were that the alcohol and water ratio can be tuned to control ionomer distribution. Water-rich solvents produce a more dispersed ink which results in better high-current density performance.

36 MATERIALS SCIENCE↗

2015 Advanced Site Investigation and Monitoring Report Riverton, Wyoming, Processing Site September 2016

The U.S. Department of Energy conducted initial groundwater characterization of the Riverton, Wyoming, Processing Site in the 1990s. The characterization culminated in a Site Observational Work Plan in 1998 that recommended a natural flushing compliance strategy. Results of verification monitoring indicated that natural flushing was generally progressing as expected until June 2010, when significant increases in contaminant concentrations were measured in several monitoring wells downgradient of the site after the area flooded. In response to the unexpected results following the flood, an enhanced characterization of the surficial aquifer was conducted in 2012, which included installation of 103 boreholes along nine transects with a Geoprobe, collection of 103 water samples and 65 soil samples, laboratory tests on the soil samples, and additional groundwater modeling. This advanced site investigation report summarizes additional investigation in 2015 through the use of backhoe trenching, sonic drilling, multilevel monitoring wells, direct-push drilling, and temporary well points to collect soil and groundwater samples. Additional surface water measurements were made included the installation of a stilling well and the measurement of stream elevation along the Wind River to approximate upgradient groundwater heads. Groundwater sampling included the addition of geochemical constituents and isotopes that have not been sampled in the past to better understand post-flood conditions and the possibility of additional or ongoing contaminant sources. This sampling was performed to (1) better define the contaminant plumes, (2) verify the occurrence of persistent secondary contaminant sources, (3) better understand the reason for the contaminant spikes after a 2010 flood, and (4) assess contaminant plume stagnation near the Little Wind River. This report provides data analyses and interpretations for the 2015 site investigation that addresses these issues and provides recommendations for future efforts. Observations from trenches and sonic drilling indicate the general lithology of the shallow, unconsolidated sediments consists of a silt zone at the surface that ranges from 2.5 to 4.8 feet below ground surface, underlain by sand and gravel, underlain by the top of the weathered bedrock (Wind River Formation). Soil data from trenches and sonic drilling indicate (1) elevated concentrations of several constituents in the silt zone, likely due to the formation of evaporites, (2) uranium is the only measured element that appears to be concentrated in the silt over the groundwater contaminant plume, (3) in the former tailings impoundment area, there may be a thin unsaturated zone with elevated uranium in the native material just below the fill, (4) in the former tailings impoundment area, slightly higher uranium concentrations occur in the underlying saturated sand and gravel, and (5) several bedrock samples have a unique geochemical signature, generally related to a higher silt content. Assessment of groundwater flow included measuring river elevations along the Wind River and installing the temporary well points adjacent to the Little Wind River that provided additional data points to refine contours for water table elevations. These data confirm past interpretations of groundwater flow to the southeast across the site toward the Little Wind River. Hydraulic head elevations between paired surficial and semiconfined aquifer wells indicate variable vertical gradients across the site with the potential for upward and downward flow. Additional direct-push drilling and groundwater sampling confirmed the contaminant plume configuration, but it also revealed a low-sulfate-concentration zone at the edge of the former tailings impoundment. Temporary well points provided better definition of plume concentrations at the bank of the Little Wind River, and data from these wells indicate plume discharge to the river. Additional sampling in an area southwest of the plume that had elevated uranium groundwater concentrations in the past did not have any uranium concentrations above the U.S. Environmental Protection Agency maximum concentration limit for uranium. Results from multilevel monitoring wells indicate some geochemical differences with depth, but overall concentrations are similar to those in nearby conventional monitoring wells in the long-term monitoring program. Geochemistry data from these multilevel monitoring wells confirm the general increase in contaminant concentrations toward the river and toward the plume centerline for chloride, sulfate, and uranium but highlight geochemical controls on calcium. Iron data indicate slightly reducing conditions, especially near water table and bedrock surfaces, with more oxidizing conditions in the middle of the sand and gravel. Uranium activity ratios (234U/238U) confirm the uranium plume in the surficial aquifer as being mill related, and the area to the southwest outside the plume as natural, non-mill related. In the semiconfined aquifer, evidence of aquifer connection and impacts from the mill is inconclusive. Values of δD and δ 18 O suggest water is derived from different sources and tritium data confirm that the semiconfined aquifer water is generally older than the surficial aquifer water. However, these data do indicate some groundwater communication from the surficial aquifer into the semiconfined aquifer, which resulted in δD and δ 18 O and tritium values in the semiconfined aquifer that are more similar to those in the surficial aquifer. Values of δ 34 S sulfate in the semiconfined aquifer combined with sulfate concentrations indicate the potential for some mill- related sulfate in the semiconfined aquifer, albeit limited to an area near and beneath the former tailings impoundment. Uranium and molybdenum concentrations in the semiconfined aquifer are below groundwater standards in all wells. However, the elevated molybdenum concentrations at one semiconfined aquifer well underneath the former tailings impoundment suggests a mill- derived source for the uranium and molybdenum in that well. It is possible that aquifer cross- communication occurred when the tailings impoundment was active and created a higher head. Current cross-communication appears unlikely given the large differences in tritium values and an upward hydraulic head at this location. Nine domestic wells are located within the institutional control boundary (eight in the confined aquifer and one in the semiconfined aquifer). Uranium and molybdenum concentrations in all samples collected from these wells were one or two orders of magnitude below the groundwater standards. Surface water flow in the Little Wind River in September 2015 was low compared to historical averages for that time of year. As a result, the uranium concentration measured in the Little Wind River was at a historical maximum at the sampling location upstream of the site. However, the impact of uranium discharge from the groundwater plume into the Little Wind River was not measureable. Elevated sulfate concentrations were observed in an outfall ditch related to an active sulfuric acid plant. Uranium concentrations in the oxbow lake remain at concentrations above the groundwater standard. Plume contaminant concentrations had returned to levels found prior to the 2010 flood by the end of 2015. However, these concentrations still exceed model predictions for natural flushing, and the current data indicate that natural flushing to achieve remediation goals within the 100-year time period is not likely, especially with the high potential for additional floods in the update to the conceptual site model (CSM), soil data indicate additional contaminant sources, specifically uranium, in evaporites within the silt layer over the uranium plume and in naturally reduced zones (NRZs). Additional zones of slightly elevated uranium concentration are in the native sediments just above the water table but below the fill layer in the former tailings impoundment area. This area also has slightly elevated uranium in the sand and gravel below the water table. Mass balance calculations indicate that small amounts of dissolution in any of these zones with increased uranium in the solid phase can produce groundwater uranium concentrations above the groundwater standard and could account for the post-flood uranium spike. The additional uranium near the former tailings impoundment provides a mechanism for a continuing source for the uranium plume that was not considered in earlier natural flushing models. In addition, uranium released from the silt layer or the NRZs seasonally and during flooding may add uranium to the groundwater plume near the Little Wind River. These mechanisms provide a possible explanation for plume persistence, along with spikes in concentrations during floods, that creates the current plume configuration. Additional updates to the CSM include (1) chloride flushes more rapidly than uranium beneath the former mill site, (2) chloride in the silt layer provides a scenario in which chloride cannot be used as a conservative tracer (especially in areas prone to flooding), (3) uranium concentrations with depth can be variable (especially below NRZs), and (4) calcite and gypsum solubility limits appear to provide important geochemical controls on groundwater geochemistry. The conclusion of this study provides several recommendations for additional work to refine the CSM and continue assessment of the natural flushing compliance strategy. Recommendations for additional work include targeted soil and groundwater sampling to assess geochemical conditions, distribution of contaminants, and groundwater/surface water interaction; additional column tests to provide data for geochemical modeling; and development of an updated groundwater flow model, which will be used in conjunction with a geochemical model to assess the viability of the natural flushing compliance strategy.

2015 Monitoring Riverton, WY↗

Recovery Act: Advanced Load Identification and Management for Buildings

In response to the U.S. Department of Energy (DoE)’s goal of achieving market ready, net-zero energy residential and commercial buildings by 2020 and 2025, Eaton partnered with the Department of Energy’s National Renewable Energy Laboratory (NREL) and Georgia Institute of Technology to develop an intelligent load identification and management technology enabled by a novel “smart power strip” to provide critical intelligence and information to improve the capability and functionality of building load analysis and building power management systems. Buildings account for 41% of the energy consumption in the United States, significantly more than either transportation or industrial. Within the building sector, plug loads account for a significant portion of energy consumption. Plug load consumes 15-20% of building energy on average. As building managers implement aggressive energy conservation measures, the proportion of plug load energy can increase to as much as 50% of building energy leaving plug loads as the largest remaining single source of energy consumption. This project focused on addressing plug-in load control and management to further improve building energy efficiency accomplished through effective load identification. The execution of the project falls into the following three major aspects; An intelligent load modeling, identification and prediction technology was developed to automatically determine the type, energy consumption, power quality, operation status and performance status of plug-in loads, using electric waveforms at a power outlet level. This project demonstrated the effectiveness of the developed technology through a large set of plug-in loads measurements and testing; A novel “Smart Power Strip (SPS) / Receptacle” prototype was developed to act as a vehicle to demonstrate the feasibility of load identification technology as a low-cost, embedded solution; and Market environment for plug-in load control and management solutions, in particular, advanced power strips (APSs) was studied. The project evaluated the market potential for Smart Power Strips (SPSs) with load identification and the likely impact of a load identification feature on APS adoption and effectiveness. The project also identified other success factors required for widespread APS adoption and market acceptance. Even though the developed technology is applicable for both residential and commercial buildings, this project is focused on effective plug-in load control and management for commercial buildings, accomplished through effective load identification. The project has completed Smart Receptacle (SR) prototype development with integration of Load ID, Control/Management, WiFi communication, and Web Service. Twenty SR units were built, tested, and demonstrated in the Eaton lab; eight SR units were tested in the National Renewable Energy Lab (NREL) for one-month of field testing. Load ID algorithm testing for extended load sets was conducted within the Eaton facility and at local university campuses. This report is to summarize the major achievements, activities, and outcomes under the execution of the project.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗

Endomembrane Cation Transporters and Membrane Trafficking

Multicellular, as well as unicellular, organisms have evolved mechanisms to regulate ion and pH homeostasis in response to developmental cues and to a changing environment. The working hypothesis is that the balance of fluxes mediated by diverse transporters at the plasma membrane and in subcellular organelles determines ionic cellular distribution, which is critical for maintenance of membrane potential, pH control, osmolality, transport of nutrients, and protein activity. An emerging theme in plant cell biology is that cells respond and adapt to diverse cues through changes of the dynamic endomembrane system. Yet we know very little about the transporters that might influence the operation of the secretory system in plants. Here we focus on transporters that influence alkali cation and pH homeostasis, mainly in the endomembrane/ secretory system. The endomembrane system of eukaryote cells serves several major functions: i) sort cargo (e.g. enzymes, transporters or receptors) to specific destinations, ii) modulate the protein and lipid composition of membrane domains through remodeling, and iii) determine and alter the properties of the cell wall through synthesis and remodeling. We had uncovered a novel family of predicted cation/H + exchangers (CHX) and K + efflux antiporters (KEA) that are prevalent in higher plants, but rare in metazoans. We combined phylogenetic and transcriptomic analyses with molecular genetic, cell biological and biochemical studies, and have published the first reports on functions of plant CHXs and KEAs. CHX studied to date act at the endomembrane system where their actions are distinct from the better-studied NHX (Na/K-H + exchangers). Arabidopsis thaliana CHX20 in guard cells modulate stomatal opening, and thus is significant for vegetative survival. Other CHXs ensure reproductive success on dry land, as they participate in organizing pollen walls, targeting of pollen tubes to the ovule or promoting fertilization. Based on localization and mutant analyses, we conclude that CHXs modulate the ion balance, pH or both in micro-regions of endoplasmic reticulum, endosomes and prevacuolar compartment (PVC), and so influence membrane trafficking and signaling resulting in proper osmoregulation in guard cells and seed formation. We also demonstrated for the first time that AtKEA2 associates with chloroplasts, especially at the two poles of developing plastids. These results show that AtKEA1 and AtKEA2 transporters in specific microdomains of the inner envelope link local osmotic, ionic, and pH homeostasis to plastid division and thylakoid membrane formation. The first 3-D structure model of AtCHX was generated, and architecture-directed mutagenesis identified critical residues of the transport core giving insights to the transport mode of this family. Thus we have revealed for the first time crucial roles of an unknown K + /H + transport family on plant growth (KEA), gas exchange, pollen cell wall, and different phases of reproduction (CHXs). The dynamic endomembrane of plant cells is integral to cytokinesis, cell expansion, defense, and cell wall formation, thus these studies are directly relevant to the mission of the Department of Energy and to a better understanding of determinants for enhancing plant biomass and plant tolerance to abiotic stress.

54 ENVIRONMENTAL SCIENCES↗

Theory and Simulation of Ultrafast Multidimensional Nonlinear X-ray Spectroscopy of Molecules (Final Report)

Emerging X-ray free electron laser (XFEL) beam sources offer new types of probes of matter with unprecedented spatial and temporal resolutions. These experimental advances must be met by robust theoretical and computational tools that provide predictive modeling capacity of the underlining electronic and structural dynamics. The latter will be essential for the design of sophisticated multi-pulse experiments and for their interpretation. The proposed research effort will focus on developing cutting- edge simulation tools for nonlinear multidimensional X-ray/optical spectroscopies and aims to address key questions in Priority Research Opportunities 1 (Probing and controlling electron motion within a molecule) and 3 (Capturing rare events and intermediate states in the transformation of matter) as mentioned in the BES Roundtable Report “Opportunities for Basic Research at the Frontiers of XFEL Ultrafast Science”. XFEL multidimensional nonlinear techniques, which combine sequences of X-ray and possibly optical pulses, provide a unique experimental toolbox for probing the dynamics of core and valence electronic excitations, as well as material structure. Predictive modeling of these dynamical processes requires the combination of analytical theory for nonlinear interactions of light and matter, robust quantum-chemical methodologies for the accurate description of electronic structure of various materials, and multiscale ab initio electron and nuclear dynamics techniques operating beyond Born-Oppenheimer approximation. These challenges will be addressed with three research thrusts (i) Develop and implement theoretical apparatus for modeling a broad range of multidimensional spectroscopic techniques enabled by present and upcoming XFEL facilities. This thrust also includes the incorporation of a computational module in the DOE supported open-source NWChem computational chemistry package as well as the development of other open-source codes ready for dissemination across a broad user base; (ii) Propose and design new multi-pulse experiments that make use of the capabilities of the incoming LCLS-II facility; (iii) Perform selected applications to specific molecular systems that can be carried out at LCLS-II and demonstrate how these X-ray sources may be used to study nonadiabatic dynamics through conical interactions, electronic correlations in multi-core excitons, and charge transfer/energy transfer processes. The proposed research will be carried out by a multi-disciplinary four-institution research team which combines academia and national laboratories and spans the broad and necessary expertise in theoretical spectroscopy, nonlinear optics, quantum chemistry, molecular non-adiabatic dynamics and code development. The work will be performed in a highly interactive team environment with junior researchers shared between institutions thus cementing cross- disciplinary interactions. The developed simulation tools will be immediately deployed for XFEL facility users, both experimentalists and theorists, via freely distributed codes and databases. Altogether, this project will facilitate establishing XFEL-based multidimensional spectroscopies as a novel diagnostic tool for monitoring electronic and structural dynamics in molecular materials.

74 ATOMIC AND MOLECULAR PHYSICS↗

Predictive performance of multi-model ensemble forecasts of COVID-19 across European nations

Background: Short-term forecasts of infectious disease burden can contribute to situational awareness and aid capacity planning. Based on best practice in other fields and recent insights in infectious disease epidemiology, one can maximise the predictive performance of such forecasts if multiple models are combined into an ensemble. Here, we report on the performance of ensembles in predicting COVID-19 cases and deaths across Europe between 08 March 2021 and 07 March 2022. Methods: We used open-source tools to develop a public European COVID-19 Forecast Hub. We invited groups globally to contribute weekly forecasts for COVID-19 cases and deaths reported by a standardised source for 32 countries over the next 1–4 weeks. Teams submitted forecasts from March 2021 using standardised quantiles of the predictive distribution. Each week we created an ensemble forecast, where each predictive quantile was calculated as the equally-weighted average (initially the mean and then from 26th July the median) of all individual models’ predictive quantiles. We measured the performance of each model using the relative Weighted Interval Score (WIS), comparing models’ forecast accuracy relative to all other models. We retrospectively explored alternative methods for ensemble forecasts, including weighted averages based on models’ past predictive performance. Results: Over 52 weeks, we collected forecasts from 48 unique models. We evaluated 29 models’ forecast scores in comparison to the ensemble model. We found a weekly ensemble had a consistently strong performance across countries over time. Across all horizons and locations, the ensemble performed better on relative WIS than 83% of participating models’ forecasts of incident cases (with a total N=886 predictions from 23 unique models), and 91% of participating models’ forecasts of deaths (N=763 predictions from 20 models). Across a 1–4 week time horizon, ensemble performance declined with longer forecast periods when forecasting cases, but remained stable over 4 weeks for incident death forecasts. In every forecast across 32 countries, the ensemble outperformed most contributing models when forecasting either cases or deaths, frequently outperforming all of its individual component models. Among several choices of ensemble methods we found that the most influential and best choice was to use a median average of models instead of using the mean, regardless of methods of weighting component forecast models. Conclusions: Our results support the use of combining forecasts from individual models into an ensemble in order to improve predictive performance across epidemiological targets and populations during infectious disease epidemics. Our findings further suggest that median ensemble methods yield better predictive performance more than ones based on means. Our findings also highlight that forecast consumers should place more weight on incident death forecasts than incident case forecasts at forecast horizons greater than 2 weeks. Funding: AA, BH, BL, LWa, MMa, PP, SV funded by National Institutes of Health (NIH) Grant 1R01GM109718, NSF BIG DATA Grant IIS-1633028, NSF Grant No.: OAC-1916805, NSF Expeditions in Computing Grant CCF-1918656, CCF-1917819, NSF RAPID CNS-2028004, NSF RAPID OAC-2027541, US Centers for Disease Control and Prevention 75D30119C05935, a grant from Google, University of Virginia Strategic Investment Fund award number SIF160, Defense Threat Reduction Agency (DTRA) under Contract No. HDTRA1-19-D-0007, and respectively Virginia Dept of Health Grant VDH-21-501-0141, VDH-21-501-0143, VDH-21-501-0147, VDH-21-501-0145, VDH-21-501-0146, VDH-21-501-0142, VDH-21-501-0148. AF, AMa, GL funded by SMIGE - Modelli statistici inferenziali per governare l'epidemia, FISR 2020-Covid-19 I Fase, FISR2020IP-00156, Codice Progetto: PRJ-0695. AM, BK, FD, FR, JK, JN, JZ, KN, MG, MR, MS, RB funded by Ministry of Science and Higher Education of Poland with grant 28/WFSN/2021 to the University of Warsaw. BRe, CPe, JLAz funded by Ministerio de Sanidad/ISCIII. BT, PG funded by PERISCOPE European H2020 project, contract number 101016233. CP, DL, EA, MC, SA funded by European Commission - Directorate-General for Communications Networks, Content and Technology through the contract LC-01485746, and Ministerio de Ciencia, Innovacion y Universidades and FEDER, with the project PGC2018-095456-B-I00. DE., MGu funded by Spanish Ministry of Health / REACT-UE (FEDER). DO, GF, IMi, LC funded by Laboratory Directed Research and Development program of Los Alamos National Laboratory (LANL) under project number 20200700ER. DS, ELR, GG, NGR, NW, YW funded by National Institutes of General Medical Sciences (R35GM119582; the content is solely the responsibility of the authors and does not necessarily represent the official views of NIGMS or the National Institutes of Health). FB, FP funded by InPresa, Lombardy Region, Italy. HG, KS funded by European Centre for Disease Prevention and Control. IV funded by Agencia de Qualitat i Avaluacio Sanitaries de Catalunya (AQuAS) through contract 2021-021OE. JDe, SMo, VP funded by Netzwerk Universitatsmedizin (NUM) project egePan (01KX2021). JPB, SH, TH funded by Federal Ministry of Education and Research (BMBF; grant 05M18SIA). KH, MSc, YKh funded by Project SaxoCOV, funded by the German Free State of Saxony. Presentation of data, model results and simulations also funded by the NFDI4Health Task Force COVID-19 ( https://www.nfdi4health.de/task-force-covid-19-2 ) within the framework of a DFG-project (LO-342/17-1). LP, VE funded by Mathematical and Statistical modelling project (MUNI/A/1615/2020), Online platform for real-time monitoring, analysis and management of epidemic situations (MUNI/11/02202001/2020); VE also supported by RECETOX research infrastructure (Ministry of Education, Youth and Sports of the Czech Republic: LM2018121), the CETOCOEN EXCELLENCE (CZ.02.1.01/0.0/0.0/17-043/0009632), RECETOX RI project (CZ.02.1.01/0.0/0.0/16-013/0001761). NIB funded by Health Protection Research Unit (grant code NIHR200908). SAb, SF funded by Wellcome Trust (210758/Z/18/Z).

60 APPLIED LIFE SCIENCES↗

Laboratory and Lysimeter Experimentation and Transport Modeling of Neptunium and Strontium in Savannah River Site Sediments

The Savannah River Site (SRS) conducts performance assessment (PA) calculations to determine the appropriate amount of low-level radiological waste that can be safely disposed on site. Parameters are included in these calculations that account for the interaction between the immobile solid phase and the mobile aqueous phase. These parameters are either the distribution coefficient (K$_{d}$ value) or the apparent solubility value (K$_{sp}$). These parameters are readily found in the literature and are used throughout the DOE complex. One shortcoming of K$_{d}$ values is that they are only applicable to a given set of solid and aqueous phase conditions. Therefore, a given radionuclide may have several K$_{d}$ values as it moves between formations and comes into contact with different solids and different aqueous phases. It is expected that the K$_{d}$ construct will be appropriate to use for a majority of the PA and for a majority of the radionuclides. However, semi-mechanistic models would be more representative in isolated cases where the chemistry is especially transitory or the radionuclide chemistry is especially complex, bringing to bear multiple species of varying sorption tendencies to the sediment. Semi-mechanistic models explicitly accommodate the dependency of K$_{d}$ values, or other sorption parameters, on contaminant concentration, competing ion concentrations, pH-dependent surface charge on the adsorbent, and solute species distribution. Incorporating semi-mechanistic concepts into geochemical models is desirable to make the models more robust and technically defensible. Furthermore, these alternative models could be used to augment or validate a K$_{d}$-based DOE Order 435.1 Performance Assessment. The objectives of this study were to: 1) develop a quantitative thermodynamically-based model for neptunium sorption to SRS sediments, and 2) determine a sorption constant from an SRS 11-year lysimeter study. The modeling studies were conducted with existing data sets. The first data set used laboratory generated Np sorption data as a function of concentration (three orders of magnitude) and as a function of pH (four orders of magnitude of proton concentration). In this modeling exercise, a very simple solution was identified by assuming that all sorption occurred only to the iron oxides in the sediment and that all the added NpO$_{4}$$^{-}$ remained in the oxidized state and was not reduced to the Np(IV) state (as occurs rapidly with Pu(V)). With rather limited input data, very good agreement between experimental and modeling results was observed. This modeling approach would be easy to add to the PA with little additional data requirements. This model would be useful in a system where pH is expected to change greatly, such as directly beneath a grout or concrete structure. The second model discussed in the report was to derive strontium K$_{d}$ values from data collected in an 11-year-old field transport study. In this controlled lysimeter study, a sensitivity analysis was conducted of hydrological and chemical processes that influence contaminant transport, including diffusion coefficients, seepage velocity, and K$_{d}$ value. The best overall K$_{d}$ derived from the model fit to the data was 32 L kg$^{-1}$, which was the same value that was previously measured in traditional laboratory batch sorption studies. This was an unexpected result given the differences in experimental conditions between the batch test and the lysimeter flow through test, in particular the differences between strontium adsorption and desorption processes occurring in the latter test and not in the former. There were some trends in the lysimeter strontium data that were not predicted by the K$_{d}$ model, which suggest that other geochemical processes are likely also controlling strontium transport. Strontium release and cation exchange are being evaluated. These results suggest that future modeling efforts (e.g., PAs) could be improved by employing a more robust semi-empirical modeling approach to transient or complex conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Solar Resource Calibration, Measurement, and Dissemination: Final Report FY 2016-FY 2018

The Solar Resource Calibration, Measurement, and Dissemination project supports the U.S. Department of Energy's Solar Energy Technologies Office (SETO) initiative to make solar energy cost-competitive with other forms of electricity by improving the tools and methods to measure and model solar radiation and therefore reduce uncertainty in predicting solar output and improve the bankability of solar projects. This project has three tasks that conduct research on advancing solar resource measurements and modeling: Task 1: Applied solar Radiation Measurements. Task 2: National Solar Radiation Data Base (NSRDB). Task 3: Knowledge Sharing. The National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) activities ensure traceable solar measurements throughout the United States that are essential for project feasibility, due diligence, financing, and plant operations. The SRRL houses NREL's Baseline Measurement System (BMS), which provides a platform for improving the accuracy of solar measurements through research to improve the calibration and characterization of instruments and enable the development of new measurement technology and standards. The platform contains a variety of instruments measuring solar radiation and other meteorological parameters. Following best practices, BMS instruments are calibrated to ensure traceability to the World Radiometric Reference in compliance with ISO/IEC 17025 accreditation requirements. The BMS data sets are quality-controlled and made publicly available through the Measurement and Instrumentation Data Center website. The NSRDB is the most accessed public database of solar irradiance and associated weather parameters for use in energy modeling. It supports SETO's goals of reducing barriers to high penetrations of solar technologies by providing easy access to high-quality, foundational data essential for innovative product development and downstream modeling through the NSRDB website. Data are currently available for the contiguous United States and beyond for 20 years covering 1998-2017. The Physical Solar Model (PSM) that underpins the current NSRDB data sets is a unique physics-based model that has been developed at NREL and has opened the door to the use of next-generation satellite data sets for solar resource assessment and forecasting. Through the American Society for Testing Materials (ASTM) International and other international standards organizations, NREL continues to provide leadership in the development of standards that are relevant to the measurement and modeling of solar radiation for solar energy applications. Standards developed cover the calibration method of pyrheliometers (direct normal irradiance measurements), pyranometers (global horizontal irradiance measurements), spectroradiometers, ultraviolet radiometers, correct field deployment of these instruments, spectral distribution reference tables for photovoltaic modeling, and transfer of indoor calibrations. Further, NREL leads the International Energy Agency (IEA) Photovoltaic Power Systems Programme Task 16: Subtask 1, which focuses on the evaluation of current and emerging resource assessment methodologies. NREL and the IEA updated the Best Practices Handbook for the Collection and Use of Solar Resource Data for Solar Energy Applications, which is widely used by the solar industry.

14 SOLAR ENERGY↗