Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Circulating Water System”

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 145 records · Page 8

Tracking Aerosol Convection Interactions Experiment (TRACER) Field Campaign Report

Convective clouds serve a critical role in the Earth’s energy and water cycles through their transport of heat, moisture, momentum, and chemical species through the troposphere driving the global circulation (e.g., Hartmann et al. 1984, Del Genio et al. 2012, Su et al. 2014). On more local scales, convective clouds impact the atmospheric heating profile through diabatic heating effects, removal of water from the atmospheric column through precipitation, and conditioning of the local environment impacting further development of clouds (e.g., Sullivan and Voigt 2021). These critical roles underscore the importance of realistic representation of convective processes across scales of models from large-eddy simulation (LES), to convection-permitting models (CPM; e.g., Kendon et al. 2020, Marinescu et al. 2021), to numerical weather prediction (NWP) models used for operational weather forecasting, to Earth system models used to predict climate sensitivity (Sanderson et al. 2011, Sherwood et al. 2014, Tomassini et al. 2014, Zhao et al. 2016, Cronin et al. 2017). A key component of improving model representation of convective clouds is better quantification and parameterization of updraft microphysics and dynamics, including their interactions with the surrounding environment and storm organization (Bony et al. 2015, Hagos and Houze 2016, Donner et al. 2016, Morrison et al. 2020). Aerosol is an important environmental factor that could affect convective clouds and precipitation since cloud droplet and ice formation processes are initiated by it. Andrae et al. (2004) hypothesized that aerosols associated with increased biomass burning particles acting as cloud condensation nuclei (CCN) result in smaller and more monodisperse cloud droplets leading to suppression of warm rain formation, ultimately leading to more cloud water being lofted above the freezing level based on observations in the Amazon region. The subsequent increase in latent heat release increases the buoyancy of rising convective parcels invigorating the deep convection. This work was followed by a description of the theoretical basis for this “cold-phase invigoration” by Rosenfeld et al. (2008), who argued that it could have a significant effect for deep convective clouds with warm cloud-bases. Several modeling studies (e.g., Khain et al. 2005, 2009, van den Heever et al. 2006, Fan et al. 2007, 2009, 2012, Lee et al. 2008, Storer et al. 2010, Lebo et al. 2012, Storer and van den Heever 2013, Chen et al. 2020, Dagan et al. 2022) have investigated these aerosol-convection interactions and the environmental factors that influence their relative importance and magnitude. More recently, several studies have indicated that “warm-phase invigoration”, the enhancement of convection through condensational heating, also appears to play a role in enhancing both shallow cumuli (Seiki and Nakajima 2014, Saleeby et al 2015) and deeper tropical convection (Lebo and Seinfeld 2011, Khain et al. 2012, Sheffield et al 2015, Fan et al. 2018, Igel and van den Heever 2021), as well as Houston thunderstorms (Fan et al. 2007, 2020). However, still other studies have provided additional evidence of systematic biases in simulated convective outflow ice size distribution properties, which are consistent with a lack of poorly understood secondary ice production within convective updrafts (e.g., Fridlind et al. 2017). To help address these critical gaps in our understanding of cloud processes, aerosol processes and aerosol-cloud interactions, the Tracking Aerosol Convection Interactions Experiment was designed building upon efforts by the Aerosol, Cloud, Precipitation and Climate (ACPC) Initiative (http://acpcintiative.org/), a joint effort of the International Geosphere-Biosphere Programme (IGBP) and the World Climate Research Program (WCRP) that focused on resolving uncertainties in the interactions between aerosol and clouds towards better understanding the role that these interactions play in the climate system. The TRACER campaign was motivated by recommendations from a number of pilot studies undertaken by ACPC (van den Heever et al. 2017, Fridlind et al. 2019, Hu et al. 2019, Fan et al. 2020, Marinescu et al. 2021, Hernandez-Deckers et al. 2022) that pointed towards the southeastern Texas region as a locale where aerosol-convection interactions could be studied owing to the copious occurrence of isolated convection during the summer months accompanied by diverse and significant sources of aerosols from both anthropogenic and natural sources. The TRACER campaign began on 01 October 2021 and extended through 30 September 2022 with an intensive operational period (IOP) during June-September 2022. Three main sites (Table 1) were managed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility.

54 ENVIRONMENTAL SCIENCES↗

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗

Shallow-to-Deep Convective Transition in the Amazon (Final Report)

The frequent and extensive deep cloud formations over the Amazon Basin provide a vast area of atmospheric heating in the Tropics, forcing large-scale global circulations on daily to secular time scales. Motivated by the importance of these systems to the global climate system and the well-documented seasonal cycle in rainfall over the Amazon Basin, this study focused on improving understanding of the processes governing these systems on seasonal time scales and their impact on the typical growth of convection from early afternoon through evening, termed the shallow-to-deep transition. This study is based on observations from the Department of Energy (DOE) ARM Mobile Facility (AMF) near Manacapuru, Brazil collected during the Observations and Modeling of the Green Ocean Amazon (GOAmazon) field campaign from 1 January 2014 to 30 November 2015, together with total column water vapor (CWV) and surface meteorological observations from a dense Global Navigation Satellite System (GNSS) network within 50 km of the AMF site provided by the CHUVA Program collected in coordination with the GOAmazon campaign. The specific objectives of this study were to address the following questions: 1) What environmental conditions support the shallow-to-deep convective transition over the Central Amazon? 2) How does the relative importance of these conditions compare during the wet and dry seasons? 3) How are these conditions modulated by basin-scale transient forcing? 4) How well are these processes captured by regional model simulations convective resolving scale? Our main hypothesis was that intra-seasonal atmospheric disturbances were a primary driver of rainfall during the wet season in the Central Amazon, while localized surface heating of the land was the primary driver of enhanced rainfall during the dry season. Overall, our findings supported our hypothesis, with the convectively active phase of atmospheric Kelvin waves strongly favoring the development of larger, more organized rain systems during the wet season, while localized isolated deep convection dominated the types of rain events seen during the dry season. These results have implications for the predictability of rainfall over the Central Amazon and for its representation in climate and weather models.

54 ENVIRONMENTAL SCIENCES↗

Signatures of aerosol-cloud interactions in GiOcean: a coupled global reanalysis with two-moment cloud microphysics

Aerosols influence the Earth's radiative balance through direct interactions with radiation and by affecting cloud properties. Anthropogenic aerosols have led to cooling during the industrial era through aerosol–cloud interactions (ACI), including aerosol effects on cloud microphysical properties and the subsequent adjustments. However, large uncertainties remain in Earth system models (ESMs) regarding the magnitude of this cooling. In part, ESMs substantially disagree on cloud properties, thermodynamics, the hydrological cycle, and general circulation. Reanalysis provides a useful avenue for exploring the impact of ACI on clouds and radiation because its atmosphere is forced to match realistic conditions through the assimilation of observations. Here, we explore the impact of ACI on clouds in the GiOcean reanalysis – the first to incorporate aerosol-cloud interactions. We contrast variables important for ACI between GiOcean and satellite observations and develop 2-dimensional lookup tables of ACI for both using a source-sink budget perspective to attribute the changes in cloud droplet number (Nd) and liquid water path (LWP) to aerosol and meteorology. A compositing analysis using lookup tables shows that GiOcean captures key aspects of aerosol–cloud–precipitation interactions, including (1) activation of aerosol into cloud droplets, (2) effective precipitation scavenging of Nd, (3) suppression of precipitation by high Nd in regions with heavy aerosol emissions. In contrast, satellite observations do not exhibit clear patterns for processes (2) and (3). Random Forest analysis shows that interannual variability in Nd and LWP over the Northern Hemisphere ocean in GiOcean is primarily driven by precipitation, consistent with satellite observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vertical Variations of Cloud Microphysical Relationships in Marine Stratocumulus Clouds Observed During the ACE‐ENA Campaign

Abstract This study examines the vertical variations of cloud microphysical relationships and their implications to cloud microphysical processes in marine stratocumulus clouds using in‐situ aircraft observations during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE‐ENA) field campaign. A new diagram with a coordinate system based on cloud droplet liquid water content (L c ) and phase relaxation time scale is proposed to investigate mixing mechanisms. This new diagram analysis shows that the inhomogeneous mixing trait is dominant near the cloud top, but homogeneous mixing trait is stronger at lower altitudes. The relevant scale parameters (i.e., transition length scale and transition scale number) also indicate a high likelihood of inhomogeneous mixing. The relationship between L c and standard deviation of droplet radius (σ R ) clearly shows the vertical transition: the correlation between L c and σ R is positive at lower cloud altitudes, but it becomes negative as altitude increases. Such a vertical transition is consistent with the vertical circulation mixing, modulating the cloud microphysical relationships to suggest homogeneous mixing at a significant depth from the cloud top.

54 ENVIRONMENTAL SCIENCES↗

A Macroalgal Cultivation Modeling System (MACMODS): Evaluating the Role of Physical-Biological Coupling on Nutrients and Farm Yield

Offshore aquaculture has the potential to expand the macroalgal industry. However, moving into deeper waters requires suspended structures that will present novel farm-environment interactions. Here, we present a computational modeling framework, the Macroalgal Cultivation Modeling System (MACMODS), to explore within-farm modifications to light, seawater flow, and nutrient fields across time and space scales relevant to macroalgae. A regional ocean model informs the site-specific setting, the Santa Barbara Channel in the Southern California Bight. A fine-scale hydrodynamic model predicts modified flows and turbulent mixing within the farm. A spatially resolved macroalgal growth model, parameterized for giant kelp, Macrocystis pyrifera , predicts kelp biomass. Key findings from model integration are that regional ocean conditions set overall farm performance, while fine-scale within-farm circulation and nutrient delivery are important to resolve variation in within-farm macroalgal performance. Therefore, we conclude that models resolving within-farm dynamics can provide benefit to farmers with insight on how farm design and regional ocean conditions interact to influence overall yield. Here, the presence of repeating longlines aligned with the mean current generate flow diversions around the farm as well as attached Langmuir circulations and increased turbulence intensity. These flow-induced phenomena lead to less biomass in the interior portion of the farm relative to the edges. We also find that there is an effluent “footprint” that extends as much as 20 km beyond the farm. In this regard, MACMODS can be used to not only evaluate farm design and cultivation practices that maximize yield but also explore interactions between the farm and ecosystem in order to minimize impacts.

54 ENVIRONMENTAL SCIENCES↗

Calculating the Effects of Solids Input and Removal as a Temperature Control in the Advanced Scale Up Reactor Experiment (ASURE) Facility at NETL Using Aspen

The Advanced Scale Up Reactor Experiment (ASURE) facility at NETL is being designed to be a fuel-flexible multi-purpose reactor that can be used for pyrolysis/gasification or evaluation of other high pressure “circulating fluidized bed” (CFB) chemical processes. The initial system design calculations for pyrolysis/gasification are presented in this work showing the expected performance of the ASURE facility when used as a biomass conversion reactor. Several other areas of application include gasification of any carbonaceous fuel including biomass, coal, plastics, and other waste materials. The reactor can therefore be used to produce SYNGAS of various compositions and hydrogen as well as other high value chemicals resulting from a typical tuned gasification process. This paper discusses an ASPEN model of the facility, focusing on the riser of the CFB reactor and the solids recirculation loop. The ASPEN model divides the riser into two sections. A bottom section which receives ash, char and sand which have been recirculated from a return loop. In this section an inert fluidization gas, (N2 or CO2), is introduced which acts as the primary mover of the solids through the system. The bottom section is equipped with a restricted air feed so that the recirculated char can be partially oxidized. This oxidation process along with the inventory of recirculating sand are used to effectively control the temperature in the following two chemical conversion sections of the reactor which are the pyrolysis zone followed by a tar cracking zone. Fresh fuel is added to the pyrolysis zone and undergoes drying and devolatilization. The products ash, char, volatile matter, and water vapor exit the pyrolysis zone and enter the reaction block for tar cracking. Steam and CO2 gasification reactions will be incorporated into the tar cracking zone, however at the design operating temperature, conversion from these reactions is expected to be essentially zero. The unit when completed in 2026 will test mixtures of biomass, plastics, and waste coal. This presentation discusses the basic ASPEN engineering design model for this project and provides preliminary sensitivity studies to determine how the various parts of the reactor will perform.

ASSURE↗

Calculating the Effects of Solids Input and Removal as a Temperature Control in the Advanced Scale Up Reactor Experiment (ASURE) Facility at NETL Using Aspen

The Advanced Scale Up Reactor Experiment (ASURE) facility at NETL is being designed to be a fuel-flexible multi-purpose reactor that can be used for pyrolysis/gasification or evaluation of other high pressure “circulating fluidized bed” (CFB) chemical processes. The initial system design calculations for pyrolysis/gasification are presented in this work showing the expected performance of the ASURE facility when used as a biomass conversion reactor. Several other areas of application include gasification of any carbonaceous fuel including biomass, coal, plastics, and other waste materials. The reactor can therefore be used to produce SYNGAS of various compositions and hydrogen as well as other high value chemicals resulting from a typical tuned gasification process. This paper discusses an ASPEN model of the facility, focusing on the riser of the CFB reactor and the solids recirculation loop. The ASPEN model divides the riser into two sections. A bottom section which receives ash, char and sand which have been recirculated from a return loop. In this section an inert fluidization gas, (N2 or CO2), is introduced which acts as the primary mover of the solids through the system. The bottom section is equipped with a restricted air feed so that the recirculated char can be partially oxidized. This oxidation process along with the inventory of recirculating sand are used to effectively control the temperature in the following two chemical conversion sections of the reactor which are the pyrolysis zone followed by a tar cracking zone. Fresh fuel is added to the pyrolysis zone and undergoes drying and devolatilization. The products ash, char, volatile matter, and water vapor exit the pyrolysis zone and enter the reaction block for tar cracking. Steam and CO2 gasification reactions will be incorporated into the tar cracking zone, however at the design operating temperature, conversion from these reactions is expected to be essentially zero. The unit when completed in 2026 will test mixtures of biomass, plastics, and waste coal. This paper presents the basic ASPEN engineering design model for this project and provides preliminary sensitivity studies to determine how the various parts of the reactor will perform.

ASSURE↗

Diversion Path Analysis: A Proposed Methodology to Develop an MC&A Approach for Liquid-Fueled Molten Salt Reactors

Nuclear material control and accounting (MC&A) is a critical element of both the US Nuclear Regulatory Commission (NRC) and US Department of Energy (DOE)’s domestic safeguards and security requirements. NRC licensees are required, under Title 10 of the Code of Federal Regulations (10 CFR) Part 74 to establish and maintain an MC&A program that captures and records the quantities and locations of special nuclear material (SNM) at the facility. Along with physical protection, MC&A is a key element of domestic nuclear material safeguards that enables the NRC to ensure that SNM is controlled and accounted for. SNM, per 10 CFR Part 74, refers to plutonium, 233 U, and uranium enriched in the isotope 233 U or 235 U, but does not include source material. Periodic physical inventories, coupled with material balance evaluations, are effective and demonstrated tools to account for and detect theft or diversion of SNM in facilities containing SNM in bulk material form (i.e., not in discrete, countable items). Historically in the United States, these types of facilities have included fuel fabrication, conversion, and enrichment facilities. In comparison, reactors have relied on item counting of assemblies and control of SNM while in containment (e.g., a sealed reactor pressure vessel) because, to date, reactor fuel has been in item form. In liquid-fueled molten salt reactors (MSRs), unlike traditional light water reactors (LWRs) or bulk facilities, bulk SNM quantities can change significantly during operation as a result of depletion and transmutation. This introduces challenges to the use of traditional periodic physical inventories and material balance evaluations to detect theft or diversion of SNM in reactors that use SNM in bulk material form. Liquid-fueled (i.e., salt-fueled) MSR facilities are MSRs that use SNM within a salt eutectic as the fuel. The SNM is in a bulk material form any time it is outside of fresh or spent fuel storage containers. Some examples of when SNM will be in bulk form in the facility are during addition of fuel to the reactor system, while fuel is circulating in operation, and while fuel is in a drain tank. Periodic physical inventories and material balance evaluations can likely be effectively applied to many portions of an MSR facility, including all areas where depletion and transmutation are not significantly changing the quantities of SNM within the control area. Within an MSR facility, this would include fresh fuel receipt and loading, waste streams that may contain SNM, irradiated fuel storage outside of the reactor core, and any irradiated fuel processing that may happen after SNM has been removed from the reactor. All of these process steps could rely on measurements of SNM quantities compared with documented inventories. Any discrepancies from predicted (i.e., book) inventories and measured inventories could be quantified as inventory differences, consistent with traditional MC&A guidance from the NRC (e.g., in NUREG-1065 Revision 2, NUREG-2159 Revision 1, and RG 5.29 Revision 2). Within the reactor system, additions and removals to the book inventory include depletion of the SNM (e.g., fission of 235 U), which complicates the use of physical inventories. SNM control, however, can also likely be effectively applied to detect theft of SNM throughout a liquid-fueled MSR facility. To complement these approaches, prior technical reports have identified that a diversion path analysis may be a useful, risk-informed, and performance-based tool to determine suitable elements of an MC&A approach for the reactor system within a liquid-fueled MSR facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Origin of water masses in Floridan Aquifer System revealed by 81 Kr

Groundwater in coastal aquifers serves as an essential resource in densely populated areas throughout the world. However, hydrological connection to the ocean leaves it at risk of salinization via changes in climate and hydrological cycles. Therefore, an accurate hydrological characterization of coastal aquifers is of the utmost importance. Although localized salinization of shallow aquifers has commonly been studied, geochemical constraints on the time scale of regional freshwater and seawater in coastal aquifers are limited especially in carbonate aquifers. Such information is not only crucial for water resource management, but also for understanding solute flux over the land-ocean boundary and for interpreting the geochemical signatures that each water mass carries. Toward this goal, we report tracer-based, subsurface residence times of groundwater in the Floridan Aquifer System in southern Florida, which is one of the most productive aquifers on Earth. The first application of Kr 81-an emerging tracer of old groundwater-to this aquifer has identified freshwater recharge during the last glacial period as anticipated from previous studies using other geochemical tools including radiocarbon, stable isotopes and noble gas concentrations. Moreover, freshwater of Holocene age was also detected at 100 km-distance from the recharge area, suggesting the possibility of an active flow system in the upgradient region. A contribution of fossil seawater that predates the last glacial maximum was locally identified, suggesting slow seawater circulation and limited but significant solute flux from the ocean into the aquifer for dolomitization. This also implies that the deep saline coastal groundwater can potentially serve as a paleo-seawater archive.

58 GEOSCIENCES↗

SAM Code Enhancements for Modeling of Liquid Metal-Cooled Fast Reactor Concepts

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes two major efforts in addressing the modeling gaps in SAM for liquid-metal-cooled fast reactor (LMFR) applications, i.e. thermal mixing and stratification phenomena in large pools and corrosion-oxidation of components in flowing lead. A new one-dimensional model for thermal mixing and stratification effects in large pools and enclosures is developed and implemented. Thermal mixing and stratification occur when fluid enters a pool at a temperature different than the bulk fluid itself, a scenario often encountered during transients in pool-type systems. These phenomena are critical for the safety of reactors, impacting phenomena like natural circulation, which is essential for passive cooling. The improved model in SAM addresses limitations of state-of-the-art approaches by combining one-dimensional (1D) channels, representing the coolant jet flow, with lumped-parameter zero-dimensional (0D) pools, representing the rest of coolant in the tank. Energy exchange between the 1D jet and the 0D pools is based on heat transfer correlations calibrated against 3D simulations. It is verified that this model can handle various flow configurations, including hot jets in colder pools, cold jets in hotter pools, and the presence of features like ceilings, free surfaces, and obstacles. Additionally, validation against experimental data demonstrates the ability of the model to capture mixing and stratification effects in a wide range of conditions. The flexibility and improved accuracy of the new model make it a valuable tool for reactor safety analysis, allowing for the simulation of different geometries encountered in advanced reactors. A system-level corrosion modeling capability is developed and implemented in SAM to support Lead Fast Reactor (LFR) development. Although the initial focus of this capability will be on LFR application, this can later be expanded to include other liquid metals such as Lead-Bismuth Eutectic (LBE) and PbLi. This report summarizes the common corrosion mitigation strategies and outlines the progress on implementing and validating a corrosion-oxidation model in SAM. Verification and validation of the corrosion-oxidation portion of the model was performed using analytical solution and measured data from samples tested in the non-isothermal pumped lead loop at IPPE Obninsk. The iron transport and corrosion/precipitation portion of the model was assessed using an analytical model and measured corrosion depths from a natural convection lead loop experiment performed at CEA. It is demonstrated that the model implemented in SAM performed well in these assessments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The DOE E3SM v1.2 Cryosphere Configuration: Description and Simulated Antarctic Ice‐Shelf Basal Melting

Abstract The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS), ice‐shelf basal melting and iceberg calving, are generally poorly represented in current Earth System Models (ESMs). Here we document the cryosphere configuration of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) v1.2. This includes simulating Antarctic ice‐shelf basal melting, which has been implemented through simulating the ocean circulation within static Antarctic ice‐shelf cavities, allowing for the ability to calculate ice‐shelf basal melt rates from the associated heat and freshwater fluxes. In addition, we added the capability to prescribe forcing from iceberg melt, allowing for realistic representation of the other dominant mass loss process from the AIS. In standard resolution simulations (using a noneddying ocean) under preindustrial climate forcing, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a transition to a high basal melt regime under the Filchner‐Ronne Ice Shelf (FRIS), presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy‐induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic ice‐shelf basal melt rates across the continent that are generally within the range inferred from observations. The accurate representation of ice‐shelf basal melting within a coupled ESM is an important step toward reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea‐level rise.

54 ENVIRONMENTAL SCIENCES↗

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗

Thermal Analysis and Simulations of the Abort Kicker Magnets Systems for the Electron-Ion Collider

The Abort Kicker Magnets system used in the current Relativistic Heavy Ion Collider (RHIC) to steer the circulating beam into the dump will be subjected to higher heat loads in the Electron-Ion Collider (EIC). After analyzing the existing abort kicker magnets and running thermal simulations in ANSYS it was concluded that they may not be suitable for use as-is in the future EIC due to heat and impedance concerns in the magnets. Possible solutions include adding a round titanium-coated ceramic beam tube to help solve the heat and impedance concerns, but this will reduce the limiting aperture of the kicker magnet. Another possible solution to meet all the performance requirements for EIC would be to add titanium-coated ceramic plates with water-cooling and tapered transitions that significantly improve the impedance and lower the heat in the magnets with less reduction to the aperture.

43 PARTICLE ACCELERATORS↗

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy's Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

closed-loop geothermal systems↗

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy’s Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

Closed-loop geothermal systems, hot-dry-rock, hot-↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

National Modeling of Geothermal District Energy Systems with Ambient-Temperature Loops Using dGeo: Preprint

Geothermal district energy systems (DES) with ambient-temperature loops, also known as thermal energy networks, are one option for decarbonizing space heating and cooling loads. Geothermal fifth-generation DES include an "ambient" temperature thermal loop that connects heat pumps at each building with thermal balancing sources such as geothermal borehole fields. Heating and cooling are provided via a water-source heat pump at each end-user. This project seeks to analyze the nationwide potential for ambient-temperature loop districts by creating a new module within the Distributed Geothermal Market Demand Model (dGeo). dGeo is an agent-based modeling tool for distributed geothermal resources; it can investigate potential on a nationwide or statewide scale using geospatial data for all 50 states and thermal demands for existing buildings. This process allows for high-level estimates of technical and economic potential for ambient-temperature loop districts across the United States. Using GHEDesigner, a lookup table was created to size borehole fields for different thermal loads and ground conditions experienced across the country. A cost and financing structure, along with incentives, were applied. Cost estimates include costs for the distribution network, borehole field installation and operation, and circulation pump operation, while savings are calculated based on agent energy bills. This newly developed module can be used for assessing which areas of the country have the highest potential for agent benefits from ambient-temperature loop installation and assess the impact of different costing and pricing future scenarios. While the code is still under development and nationwide simulations are ongoing, initial results for two states are provided. Future work includes expanding the module to consider mixed residential and commercial districts and considering multiple costing scenarios.

ambient temperature loop↗