Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “water evaluation”

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 1,063 records · Page 59

Temporal sequence transformer to advance long-term streamflow prediction

Accurate streamflow prediction is crucial for understanding climate change impacts on water resources and for effective management of extreme hydrological events. While Long Short-Term Memory (LSTM) networks have been the dominant data-driven approach for streamflow forecasting, recent advancements in transformer architectures for time series tasks have shown promise in outperforming traditional LSTM models. This study introduces a transformer-based model that integrates historical streamflow data with climatic variables to enhance streamflow prediction accuracy. We evaluated our transformer model against a benchmark LSTM across five diverse basins in the United States. Results demonstrate that the transformer architecture consistently outperforms the LSTM model across all evaluation metrics, highlighting its potential as a more effective tool for hydrological forecasting. This research contributes to the ongoing development of advanced AI techniques for improved water resource management and climate change adaptation strategies.

Singh, Ruhaan [Farragut High School]↗

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↗

Sensitivity of deep convection and cross-tropopause water transport to microphysical parameterizations in WRF

Deep convective storms can overshoot the tropopause, thus altering the composition of the stratosphere by vertically transporting tropospheric air. The transport of water vapor and ice particles into a sub-saturated environment can hydrate the stratosphere, with implications for radiative forcing and ozone chemistry. Cloud-resolved models, if employed at high spatial resolutions, are used to probe process-level questions about cross-tropopause deep convective hydration and its controls. There is considerable diversity in model representations of processes associated with water transport and transformation, and the choice of a microphysics scheme affects model skill in simulating deep convective events. This motivates our evaluation of state-of-the-art, as well as widely used standard schemes, in a high spatial- and temporal-resolution framework. Six bulk microphysics schemes were employed in a WRF-LES setup, initialized with a sounding profile representative of a tropopause-overshooting storm. Here, we used an idealized framework to isolate the effect of microphysics on the dominant processes that control the reach of deep convection and stratospheric hydration. All schemes produced the highest reaching updrafts 8-12 hours into the simulation but the strength and persistence of updrafts varied across the schemes; maximum storm heights ranged 9.1-12.6 km across the schemes. Varying microphysics produced large differences in the vertical extent and horizontal aggregation of convection, and an order of magnitude spread in above-tropopause water vapor concentrations.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Influence of Plants on Hydrologic Cycling: Quantifying and Validating the Role of Plant Processes and Stomatal Conductance (Final Report)

Plants can exert strong controls on water cycling, both locally and across the globe. The impacts are driven by how plants control the exchange of water and energy between the land and the atmosphere, which in turn are influenced by how plants photosynthesize and grow (e.g. biogeochemical cycling). Our team has shown previously that plant responses to increasing CO 2 in the atmosphere can influence the water cycle through changes in rainfall and the amount of water on land with implications for drought. Plant responses can also alter both average and extreme flow of water in rivers, with implications for freshwater availability and flood frequency. Although the impact of plant responses on water cycling has been demonstrated to exist, significant uncertainty remains in our understanding of the magnitude and form of the plant responses themselves, as well as their ultimate impact on water availability for people and ecosystems. In this project we worked to quantify the role of plant processes in regulating water cycling (e.g., precipitation, evapotranspiration, runoff, droughts, etc.). We addressed the following objectives: quantify the control of physiological vs. radiative effects on water cycling across many coupled models in idealized and realistic scenarios, as well as to quantify how our assumptions about leaf-level processes (e.g. coupling between stomatal conductance and photosynthesis), and organism-level processes (e.g. leaf area response to high CO 2 ) contributed to uncertainties and biases in water cycling. We focused on the impact of these processes on water cycling to assess the sensitivities of precipitation, evapotranspiration, and runoff to model assumptions about plant processes at the leaf, organism, and community level. As part of this project we used measurements of stable carbon isotopes to generate observationally based estimates of plant functioning during the historical period. This work is directly relevant to DOE’s primary scientific research questions in the RGMA topic area (b) on Biogeochemical processes, feedbacks, and interactions in the Earth system.

59 BASIC BIOLOGICAL SCIENCES↗

Development and Evaluation of Embrittlement Resistant Alloys for Advanced LWR Cladding

The primary challenge preventing the deployment of reactor designs that leverage wrought FeCrAl as an advanced light-water reactor (LWR) cladding is irradiation hardening and embrittlement. Wrought FeCrAl alloys experience a loss of post-irradiation ductility and fracture toughness under low-temperature neutron irradiation (< 350°C)resulting from the combined effects of dislocation loop formation and the irradiation-enhanced precipitation of Cr-rich alpha-prime precipitates throughout the microstructure. Over the past decade, significant improvements in our understanding regarding the effect of Cr and Al content have been enabled through neutron irradiations, but even optimized wrought alloys such as C26M remain vulnerable to brittle failure during storage, transportation, and handling following irradiation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development and Evaluation of Embrittlement Resistant Alloys for Advanced LWR Cladding

The primary challenge preventing the deployment of reactor designs that leverage wrought FeCrAl as an advanced light-water reactor (LWR) cladding is irradiation hardening and embrittlement. Wrought FeCrAl alloys experience a loss of post-irradiation ductility and fracture toughness under low-temperature neutron irradiation (< 350°C) resulting from the combined effects of dislocation loop formation and the irradiation-enhanced precipitation of Cr-rich alpha-prime precipitates throughout the microstructure. Over the past decade, significant improvements in our understanding regarding the effect of Cr and Al content have been enabled through neutron irradiations, but even optimized wrought alloys such as C26M remain vulnerable to brittle failure during storage, transportation, and handling following irradiation. This report summarizes multiple irradiation campaigns initiated over the past 7 years and provides recent insights into the effect of (1) minor alloying composition, (2) alloy processing strategy, (3) crystal structure, and (4) the use of advanced oxide dispersion strengthened (ODS) alloys. The first irradiation campaign irradiated body-centered cubic FeCrAl alloy C26M with and without elements such as Mo and Y. It also included face-centered cubic alumina-forming austenitic (AFA) alloys that have a steeper neutronic penalty but that have far superior high-temperature strength than wrought FeCrAl. Finally, this campaign compared conventional wrought C26M with the same material produced using modern powder metallurgical hotisostatic pressing (PM-HIP). The second irradiation campaign, initiated in 2018, irradiated wrought and FeCrAl-ODS materials up to 50 dpa to assess their irradiation resistance to high-doses. The results of this work indicate that wrought C26M suffers significant degradation following neutron irradiation at LWR-relevant temperatures regardless of Mo and Y content. The PM-HIP variant appears to be more resilient to irradiation-induced ductility loss in comparison with wrought variants but still loses almost all ductility by 8 dpa. Notably, the AFAs retained at least 5% total elongation after irradiation at the same dose and temperature condition, although additional deterioration may be expected at higher fluence levels. FeCrAl-ODS materials show the greatest promise as a transformative longer-term accident-tolerant fuel (LT-ATF) cladding material. The extruded 106ZY10C alloy retained 10% total elongation after 16 dpa irradiation and retained 8% total elongation after 50 dpa. Additional challenges with FeCrAl-ODS alloys also remain, including optimization of end cap joining methodologies, enhancement of fracture toughness, scaling production to prove economic viability, and optimizing post-pilger heat treatments to maximize ductility and irradiation resistance. This series of irradiations demonstrates the efficacy of the LWR testbed available within the United States for the rapid irradiation and down-selection of LT-ATF candidates. As the country rapidly accelerates its timeline for the deployment of advanced reactor concepts, the effective utilization of the High Flux Isotope Reactor for separate-effects style irradiations should continue to be prioritized to answer the final questions pertaining to LT-ATF candidates necessary for the deployment of advanced boiling water reactors and small modular reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Enhanced Pool Boiling of Low-Pressure Refrigerants on Round Tubes- An Experimental Evaluation

Pool boiling heat transfer is one of the most effective heat transfer methos and crucial for several energy conversion processes such as electronic cooling, steam generation, water purification and thermal management. Enhancing boiling heat transfer around the tubes enables higher heat flux, which can improve the efficiency and lead to compact heat exchangers. Delaying the onset of critical heat flux (CHF) allows for higher CHFs and protects equipment during anomalous transient conditions. Growing interest in process intensification and delayed CHF has motivated the current study. An experimental approach has been deployed to investigate a potential technique to enhance the pool boiling performance around small diameter tubes by open cell metal foams. The findings for nucleate pool boiling of Novec HFE 7000 on the externally enhanced and smooth tubes are reported. The enhanced tube comprised of compressed aluminum foam brazed around. All the samples were horizontally oriented, and the experiments were performed under ambient pressure conditions. The heat transfer measurements for the smooth tubes were compared with several existing correlations and established the baseline. The preliminary results showed that the aluminum foam effectively enhanced the heat transfer coefficient and delayed the onset of critical heat flux. In addition, the pool boiling process around the round tube sample was visualized from both side and front views to comprehend the heat transfer mechanisms underlying substantial enhancement observed during experiments.

Yang, Cheng-Min↗

Data and Scripts associated with: “Laboratory evaluation of open source and commercial electrical conductivity sensor precision and accuracy”

This data package is associated with the publication “Laboratory evaluation of open source and commercial electrical conductivity sensor precision and accuracy: How do they compare?” submitted to PLOS ONE and accepted for publication (Fulton, S.G. et al. 2023; doi not yet available).Variation in electrical conductivity (EC) of water is important to reveal environmental disturbance and natural dynamics, including factors such as anthropogenic salinization. This data package supports a study addressing the need for a robust performance assessment of open source (OS) EC sensors. We evaluated the accuracy (mean error, %) and precision (sample standard deviation) of OS EC sensors in the laboratory via comparison to EC calibration standards using three different OS and OS/commercial-hybrid (OS/C) EC sensors and data logger configurations and two commercial (C) EC sensors and data logger configurations. We also evaluated the effect of cable length (7.5 m and 30 m) and sensor calibration on OS sensor accuracy and precision. This data package presents the results from the different testing laboratory configurations. It also includes an R script for statistical analysis of the data and an Arduino IDE file used to calibrate and collect data with the OS Atlas EC sensor. File types are .csv, .pdf, .R, and .ino.

54 ENVIRONMENTAL SCIENCES↗

Geothermal Reservoir Simulation Analysis in Support of Electricity Co-Production Feasibility Study at the Blackburn Oil Field, Nevada: Preprint

Geothermal electricity co-production is a viable option for oil reservoirs producing large water cuts with elevated wellhead-observed temperatures. Repurposing existing oil wells significantly reduces initial investment costs historically associated with geothermal resource utilization. The National Renewable Energy Laboratory (NREL), partnering with Gradient Geothermal, Inc. (formerly known as Transitional Energy) and Grant Canyon Oil and Gas, has been tasked to evaluate the feasibility of geothermal electricity co-production at the Blackburn Oil Field with Organic Rankine Cycle (ORC) generators. The Devonian steady-state reservoir has historically been producing high water cuts of 240 degrees F (115.6 degrees C) observed at the wellhead without documented pressure drawdown or thermal breakthrough. An estimated initial reservoir temperature of approx. 260degreesF (126.7 degrees C) has been observed in the field and history-matched in a wellbore production analysis and reservoir simulation. Our objective was to develop a conceptual geological model of the subsurface, simulate a natural-state reservoir, model production scenarios, and complete a technical feasibility analysis to accomplish this task. Through extensive modeling and the use of available proprietary and public data, it was possible simulate three scenarios that indicated minimal thermal decline over the duration of a simulated ten-year production and re-injection scheme.

Blackburn Nevada↗

Improving the particle dry deposition scheme in the CMAQ photochemical modeling system

Dry deposition of atmospheric aerosols in large-scale models is a critical, but highly uncertain, sink process with a strong dependence on particle size, meteorological conditions, and land surface properties. This study investigates the particle dry deposition scheme implemented in the standard Community Multiscale Air Quality (CMAQ) model v5.2.1, characterizes its underlying parameterized components with comparison to a similar scheme in a contemporary regional-scale model, and proposes two updated schemes that are then evaluated with available ambient particle deposition velocity (V d ) measurements. Both updated schemes reduce the surprisingly strong dependence of deposition velocity on the aerosol mode width, with one scheme further introducing a dependence on vegetation coverage that is broadly consistent with variability in observations between vegetated and non-vegetated surfaces. Compared to the base scheme, the updated scheme with vegetation dependence increases V d for submicron particles and decreases it for larger particles by an average of 37% and –66%, respectively. This scheme performs statistically better than the base scheme, reducing fractional biases by 56%–97% for vegetated land-use types and has roughly equivalent performance over water. The base and updated schemes are tested with three annual CMAQ (v5.2.1) simulations for the year 2011; predicted ambient aerosol concentrations are evaluated with routine monitoring network observations and predicted dry deposition fluxes are evaluated with data from the Clean Air Status and Trends Network (CASTNET). The updated scheme with vegetation dependence reduces negative fractional biases for PM 10 by 41% and positive fractional biases for PM 2.5 organic carbon by 15%. This scheme has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3 and beyond) to replace the scheme used in previous versions of CMAQ (v4.5 through v5.2.1).

54 ENVIRONMENTAL SCIENCES↗

DBD plasma-assisted ethanol steam reforming for green H 2 production: Process optimization through response surface methodology (RSM)

Herein this work investigates ethanol steam reforming (ESR) to produce hydrogen (H 2 ) in a dielectric barrier discharge (DBD) plasma reactor. A five-level, three-factor experiment design was performed using a response surface methodology (RSM) to evaluate the combined effects of the three process parameters, including discharge power, total flow rate, and ethanol-to-water (EtOH/H 2 O) molar ratio on the plasma-assisted ESR reaction. Quadratic regression models were employed in RSM to fit the experimental results and present the correlation between process parameters and targeted responses (EtOH conversion, H 2 yield, H 2 selectivity, and specific energy requirement (SER) for H 2 production). The results suggested that the EtOH/H 2 O molar ratio is considered to have the most significant effect on the EtOH conversion and H 2 , H 2 selectivity, while the total flow rate is the most significant parameter determining SER for H 2 production. Process optimization demonstrated the optimal process conditions, including a discharge power of 55.9 W, a total flow rate of 26.7 ml/min, and an EtOH/H 2 O molar ratio equal to 0.34. A validation test was performed and confirmed the feasibility of the optimization process.

08 HYDROGEN↗

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗

Atomically Dispersed Metals on Well-Defined Supports including Zeolites and Metal–Organic Frameworks: Structure, Bonding, Reactivity, and Catalysis

When metals in supported catalysts are atomically dispersed, they are usually cationic and bonded chemically to supports. Investigations of noble metals in this class are growing rapidly, leading to discoveries of catalysts with new properties. Characterization of these materials is challenging because the metal atoms reside on surfaces that are typically nonuniform in composition and structure. We posit that understanding of structures and catalytic properties of these materials is emerging most strongly from investigations of structurally uniform catalysts (metal atoms dispersed on crystalline supports) which can be characterized incisively with atomic-resolution electron microscopy, X-ray absorption spectroscopy, and infrared spectroscopy, bolstered by density functional theory. We assess the literature of such catalysts supported on zeotype materials, metal–organic frameworks, and covalent organic frameworks. Finally, assessing characterization, reactivity, and catalytic performance of catalysts for oxidation, hydrogenation, the water–gas shift reaction, and others, we consider metal–support interactions and ligand effects for various metal–support combinations, evaluating the degree of structural uniformity of exemplary catalysts and summarizing structure–reactivity and structure–catalytic property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Forum on Materials and Interfaces for Energy Storage and Conversion

This ACS Applied Materials & Interfaces forum on energy storage and conversion contains the topics of research discoveries to advance energy technology development. The topics are focused on the efficient and sustainable use of energy in the fields of fuel cells, CO 2 conversion, water splitting, batteries, and electrochemical capacitor. This forum highlights recent achievements in these emerging technologies through the design and evaluation of novel materials and devices and theoretical study. Additionally, this forum also showcases the technology at all stages of development from “embryonic” conception to industrial adaptation.

25 ENERGY STORAGE↗

Impacts of abiotic and biotic factors on tundra productivity near Utqiaġvik, Alaska

Abstract Earlier snowmelt, warmer temperatures and herbivory are among the factors that influence high-latitude tundra productivity near the town of Utqiaġvik in northern Alaska. However, our understanding of the potential interactions between these factors is limited. MODIS observations provide cover fractions of vegetation, snow, standing water, and soil, and fractional absorption of photosynthetically active radiation by canopy chlorophyll (fAPAR chl ) per pixel. Here, we evaluated a recent time-period (2001–2014) that the tundra experienced large interannual variability in vegetation productivity metrics (i.e. fAPAR chl and APAR chl ), which was explainable by both abiotic and biotic factors. We found earlier snowmelt to increase soil and vegetation cover, and productivity in June, while warmer temperatures significantly increased monthly productivity. However, abiotic factors failed to explain stark decreases in productivity during August of 2008, which coincided with a severe lemming outbreak. MODIS observations found this tundra ecosystem to completely recover two years later, resulting in elevated productivity. This study highlights the potential roles of both climate and herbivory in modulating the interannual variability of remotely retrieved plant productivity metrics in Arctic coastal tundra ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Flood Estimation under Snowmelt and Rain-on-Snow Processes in Alaska: A Military Installation Perspective

Accurate flood estimation in snow-dominated and high-latitude regions remains challenging due to complex interactions among rainfall, snowmelt, and rain-on-snow (ROS) processes, which are not captured in conventional precipitation-based intensity-duration-frequency (PREC-IDF) curves. This study evaluates the Next-Generation IDF (NG-IDF) framework, an extension of PREC-IDF that incorporates total water available for runoff (precipitation minus changes in snow water equivalent), in two contrasting Alaskan watersheds of Department of Defense (DoD) significance: the Little Chena River Basin (LCRB) in interior Alaska, a tributary of the Chena River that flows through Fort Wainwright and near Eielson Air Force Base, and the Upper Ship Creek Basin (USCB), which drains Joint Base Elmendorf-Richardson near Anchorage. Using long-term SNOTEL observations and event-based rainfall-runoff modeling, NG-IDF and PREC-IDF flood estimates were compared against observation-based flood frequency analyses. Results show that NG-IDF consistently reduces flood-estimation bias by 15–20% relative to PREC-IDF, particularly for snowmelt- and ROS-dominated events. In the interior LCRB, permafrost conditions can substantially amplify flood responses during snowmelt events, an effect not explicitly represented in standard design tools. These findings demonstrate that NG-IDF provides a more physically consistent and transferable framework for flood estimation in cold regions, with potential relevance to mission-critical DoD installation resilience. Projected increases in ROS frequency and permafrost degradation across Alaska further emphasize the need to integrate physics-based hydrologic models that explicitly represent snow and permafrost processes to enhance design resilience and operational readiness for military and other critical infrastructure.

Yan, Hongxiang (ORCID:000000022387403X)↗

Human Health Risk Assessment for Improper Landfill Disposal of End-of-Life CdTe PV

The present work is a continuation of the 2020 IEA PVPS Task 12 Human Health Risk Assessment Methods for PV Part 3: Module Disposal Risks. The 2020 report performed a human health risk assessment (HHRA) for disposal of a cadmium telluride (CdTe) PV module in an unlined landfill, focusing solely on risks from cadmium. This study extends the 2020 HHRA on CdTe PV, analyzing eleven constituent elements: Cd, Se, Te, Cu, Si, Cr(III), Mo, Sn, Zn, Ni, and Al. The present HHRA was performed through two methods: utilization of the U.S. Environmental Protection Agency's (USEPA) Delisting Risk Assessment Software (DRAS V.4.0) on eight exposure pathways for cancer risk and non-cancer hazards; and comparison of exposure point concentrations to federal standards for groundwater, surface water, air, and soil exposure pathways. Cancer risks and non-cancer hazards posed by elemental leaching through all evaluated exposure pathways, using both methods, were found to be several orders of magnitude below USEPA health-protective thresholds. Cadmium exhibited both the highest risks and lowest uncertainty considering data availability on chemical content, leachate, and federal screening levels.

chemicals↗

interflow: A Python package to organize, calculate, and visualize sectoral interdependency flow data

Many economic sectors rely on an uninterrupted “upstream” supply of a resource to conduct their primary functions, leaving them vulnerable to adverse effects should that resource flow be interrupted or compromised (OECD, 2017; U.S. EPA, 2010). Well-known examples of these relationships include water demand by the energy sector (e.g., thermoelectric cooling for nuclear generation) (Grubert & Sanders, 2018; Webber, 2017) and energy demand by the water sector (e.g., electricity required to treat or move water in the public water sector) (Congressional Research Service, 2017) though many others exist. Being able to calculate and document these interdependencies and evaluate where the greatest cross-sectoral intensities and flows exist can reveal opportunities to enhance the overall network. Despite the implications and potential impacts, however, these interconnections and flows have been historically complex to analyze and understand. The interflow package provides a flexible tool to organize, calculate, and visualize (using Sankey diagrams and other visualizations) sectoral interdependency flows for multiple subsectors and resources. This tool can help decision-makers, researchers, and other audiences more easily pull meaning from these interdependencies to reveal multi-faceted opportunities and risks. interflow can help investigate questions such as (1) which sectors have high cross-resource dependencies, (2) how does demand for a resource in various sectors compare across regions, and (3) where the sectoral and regional opportunities are for enhanced efficiency, security, and resiliency.

97 MATHEMATICS AND COMPUTING↗