Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydraulic conductivity”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Electrical monitoring of saline tracers to reveal subsurface flow pathways in a flat ditch-drained field

Time-lapse electrical resistivity imaging (ERI) combined with salt tracers can enhance our understanding of subsurface flow and transport in ditch-drained agroecosystems. In a nearly level field adjacent to a drainage ditch, we installed 192 electrodes in a 72-m2 plot along with five hydrometric stations to enable time-lapse ERI monitoring during storms. Immediately prior to a wintertime rainstorm, we applied 5 kg of KBr to a shallow trench located roughly 11 m upgradient from the ditch. Antecedent conditions prior to the storm were wet, assuring a hydrologically connected state that favored tracer dissolution and movement in the subsurface. We monitored the evolution of the tracer plume with ERI over a seven-day period. Four-dimensional (4D) delineation of precipitation and tracer transport using time-lapse ERI showed that the tracer plume moved preferentially downward through a discontinuity in the argillic horizon within the first 57 hrs. Over the following 79 hrs, the tracer plume exhibited steady, lateral movement towards the drainage ditch within a more permeable sandy horizon situated just below the argillic. A spatial moment analysis of logarithmic changes in conductivity in the sand layer revealed seepage velocities ranging from 0.0027 to 0.0482 m d-1; these values were positively correlated with hydraulic gradients, suggesting laminar flow processes typical of sandy aquifers. Increasing groundwater velocities with hydraulic gradients suggest that solute transfers from groundwater to ditch are likely to be faster in storms, especially during wet periods. Even so, predominantly slow lateral transport within the sandy horizon, even during storms, indicates that critical source areas of pollutant loss are likely within a few meters of the ditch.

Robinson, Judith L.↗
DOE OSTI2023

Northern Hemisphere Snow Drought in Earth System Model Simulations and ERA5‐Land Data in 1980–2014

Abstract Low snow levels over the past few decades and predictions of a low‐to‐no snow future have spurred research into snow droughts, which pose a threat to water security and management. Systematic data‐model comparisons of snow drought have been lacking, hindering our understanding of the drivers of snow drought in the past. To address this gap, we analyzed snow drought events using standardized snow water equivalent index derived from monthly results of four numerical experiments using the E3SM Land Model (ELM) and ERA5‐Land data during the period of 1980–2014. Additionally, we compared snow drought duration calculated from models with those from the ERA5‐Land data during selected El Niño‐Southern Oscillation (ENSO) years. The numerical experiments were conducted with ELM driven by two prescribed atmospheric forcings, and with the coupled land‐atmosphere configuration of E3SM with and without plant hydraulics scheme feedback. Analysis reveals that 20%–30% of snow droughts occur due to factors other than above‐normal temperature and low snowfall, such as low soil moisture, warm soil temperature, and low relative humidity, etc., especially in high latitudes (50° North). Furthermore, our study highlights the exacerbating effect of ENSO events on snow drought conditions in various regions, despite some discrepancies between model and ERA5‐Land results. We also identified limitations of the coupled land‐atmosphere models in our current configuration in capturing the spatial patterns of snow droughts. This study underscores the challenge of predicting and mitigating snow drought and the need for a comprehensive understanding of the factors contributing to snow drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

DEMO-FTES: Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (CRADA Final Report)

The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.

25 ENERGY STORAGE↗
DOE OSTI2022

Initial development of a generic fluoride salt-cooled reactor model

Fluoride high-temperature reactors (FHRs) are high-temperature, low-pressure reactor concepts that use tri-structural isotropic (TRISO) fuel and molten fluoride salt coolant. These reactors have the potential to provide both electrical power and high-temperature process heat. We used generic FHR parameters for a pebble-bed FHR to develop an initial model with fresh fuel for a generic FHR (gFHR) in MELCOR and SCALE (NEWT and KENO). In this paper, we present the development of our gFHR models, which will serve as the baseline for a sensitivity and uncertainty analysis to quantify the range of possible source terms for FHRs in severe accidents. We present MELCOR results for fuel and coolant temperatures through the core, a nodalization study for the steady-state thermal hydraulic model, and development of reactor physics models in SCALE. As this work progresses, these models will be used to calculate source terms for a loss-of-forced-flow accident and to conduct a sensitivity study on this accident to establish a range of possible source terms. SCALE will provide reactor physics parameters like isotopic inventory, decay heat generation, and temperature coefficients of reactivity. Using the uncertainty quantification tools within SCALE, we will generate distributions for those parameters and will use the uncertainty quantification code RAVEN or DAKOTA to sample those distributions in MELCOR to quantify the impact of reactor physics and thermal hydraulic uncertainties on FHR source terms. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

RELAP5-3D in the OECD-NEA HTGR Thermal Hydraulics Benchmark

The OECD-NEA HTGR Thermal Hydraulics benchmark uses data from the High Temperature Test Facility for a series of code-to-code and code-to-data exercises aimed at improving the state of knowledge on existing thermal hydraulics modeling and simulation tools for prismatic HTGR applications. This benchmark includes problems representing hot gas mixing in the lower plenum, the depressurized conduction cooldown accident and the pressurized conduction cooldown accident. This presentation discusses the benchmark and RELAP5-3D's role in that benchmark, including as a tool for predicting behavior in the core and as a tool for providing boundary conditions to computational fluid dynamics analysis in the lower plenum.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Inclusion of a cold hardening scheme to represent frost tolerance is essential to model realistic plant hydraulics in the Arctic–boreal zone in CLM5.0-FATES-Hydro

As temperatures decrease in autumn, vegetation of temperate and boreal ecosystems increases its tolerance to freezing. This process, known as hardening, results in a set of physiological changes at the molecular level that initiate modifications of cell membrane composition and the synthesis of anti-freeze proteins. Together with the freezing of extracellular water, anti-freeze proteins reduce plant water potentials and xylem conductivity. To represent the responses of vegetation to climate change, land surface schemes increasingly employ “hydrodynamic” models that represent the explicit fluxes of water from soil and through plants. The functioning of such schemes under frozen soil conditions, however, is poorly understood. Nonetheless, hydraulic processes are of major importance in the dynamics of these systems, which can suffer from, e.g., winter “frost drought” events. In this study, we implement a scheme that represents hardening into CLM5.0-FATES-Hydro. FATES-Hydro is a plant hydrodynamics module in FATES, a cohort model of vegetation physiology, growth, and dynamics hosted in CLM5.0. We find that, in frozen systems, it is necessary to introduce reductions in plant water loss associated with hardening to prevent winter desiccation. This work makes it possible to use CLM5.0-FATES-Hydro to model realistic impacts from frost droughts on vegetation growth and photosynthesis, leading to more reliable projections of how northern ecosystems respond to climate change.

58 GEOSCIENCES↗
DOE OSTI2022

The influence of increasing atmospheric CO 2 , temperature, and vapor pressure deficit on seawater-induced tree mortality

We report increasing seawater exposure is killing coastal trees globally, with expectations of accelerating mortality with rising sea levels. However, the impact of concomitant changes in atmospheric CO 2 concentration, temperature, and vapor pressure deficit (VPD) on seawater-induced tree mortality is uncertain. We examined the mechanisms of seawater-induced mortality under varying climate scenarios using a photosynthetic gain and hydraulic cost optimization model validated against observations in a mature stand of Sitka-spruce (Picea sitchensis) trees in the Pacific Northwest, USA, that were dying from recent seawater exposure. The simulations matched well with observations of photosynthesis, transpiration, nonstructural carbohydrates concentrations, leaf water potential, the percentage loss of xylem conductivity, and stand-level mortality rates. The simulations suggest that seawater-induced mortality could decrease by ~16.7% with increasing atmospheric CO 2 levels due to reduced risk of carbon starvation. Conversely, rising VPD could increase mortality by ~5.6% because of increasing risk of hydraulic failure. Across all scenarios, seawater-induced mortality was driven by hydraulic failure in the first two years after seawater exposure began, with carbon starvation becoming more important in subsequent years. Changing CO 2 and climate appear unlikely to have a significant impact on coastal tree mortality under rising sea levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hydraulic and Thermal Stimulation Program at Raft River Idaho, A DOE EGS

A Department of Energy Enhanced Geothermal System (EGS) stimulation program has injected over 254 million gallons of water into the well RRG-9 ST1 since the summer of 2013. Three major stimulations have been conducted during the program increasing injection flow rates from less than 20 gpm to 550 gpm. Geologic, water chemistry, microseismic activity, and borehole imaging data have been used to develop a conceptual model describing possible flow paths of this injected water. This model contains two major fracture zones one of which intersects the RRG-9 ST1 wellbore. Modified Hall and injectivity index plots constructed using injection flow rates, surface temperatures, and wellhead pressures show steady improvement in the injectivity of the well. Here, the injectivity index has risen from 0.15 gpm/psi to 2.0 gpm/ psi. A pressure falloff test conducted on April 28, 2015 indicates a reservoir permeability of 1,220 md and -5.38 skin factor. The well stimulation program was simulated numerically using an Idaho National Laboratory reservoir simulation code, FALCON. These simulations show a significant increase in the permeability of connecting fracture pathways after each stimulation event.

Enhanced geothermal system↗
DOE OSTI2025

Experimental Comparison of Elastomeric Materials for Hydraulic Seal Durability under Reciprocating Conditions

Wave Energy Converters (WECs) rely on hydraulic Power Take-Off (PTO) systems to transform the kinetic energy of ocean waves into electricity, where hydraulic seals play a critical role in preventing fluid leakage and ensuring watertight operation under harsh marine conditions. However, seal degradation due to wear, fatigue, and corrosion remains a key challenge limiting system reliability and increasing maintenance costs. In this study, the tribological performance of elastomeric materials for hydraulic sealing applications that are based on commercial polyurethane and custom-compounded Ethylene Propylene Diene Monomer (EPDM) is investigated. Using a CSM tribometer, reciprocating wear tests were conducted in accordance with ASTM G133 under dry sliding conditions, with normal loads ranging from 3 to 10 N, and sliding speeds of 10 to 30 mm/s. Coefficient of friction (COF), specific wear rate, and worn surface morphology are evaluated to understand how load, velocity, and material composition affect seal durability. Improved wear stability and high Shore A hardness (93–94) were demonstrated by the commercial thermoplastic polyurethane (Duralast 4203 and 4758). EPDM-based compounds showed different filler-dependent behaviors, formulations filled with carbon black (E3) had the lowest wear rate (1.45 × 10?4 mm³ N?¹ m?¹) and increased load-bearing capacity, while samples modified with silica and plasticizer (E4, E5) had better conformability and damping but less hardness. Compounded EPDMs showed wider deformation zones because of their lower modulus and filler-matrix interfacial fluctuations, while commercial grades exhibited shallower wear scars (= 380 µm), according to optical profilometry. The Archard model-based lifetime estimation showed that, under 10 N and 30 mm/s, operational endurance ranged from 3.1 to 6.2 years. Because of its low wear rate and steady interfacial behavior, E3 had the longest anticipated life. For medium-load hydraulic sealing applications, these results demonstrate the possibility of customized EPDM formulations as affordable substitutes for commercial polyurethanes. They also offer a quantitative foundation for material selection and lifespan prediction.

Durability↗
DOE OSTI2021

Extreme Miller cycle with high intake boost for improved efficiency and emissions in heavy-duty diesel engines

This study experimentally investigates the impact of extreme Miller cycle strategies paired with high intake manifold pressures on the combustion process, emissions, and thermal efficiency of heavy-duty diesel engines. Well-controlled experiments isolating the effect of Miller cycle strategies on the combustion process were conducted at constant engine speed and load (1160 rpm, 1.76 MPa net IMEP) on a single cylinder research engine equipped with a fully-flexible hydraulic valve train system. Late intake valve closing (LIVC) timing strategies were compared to a conventional intake valve profile under either constant cylinder composition, constant engine-out NO x emission, or constant overall turbocharger efficiency (η TC ) to investigate the operating constraints that favor Miller cycle operation over the baseline strategy. Furthermore, utilizing high boost with conventional intake valve closing timing resulted in improved fuel consumption at the expense of sharp increases in peak cylinder pressures, engine-out NO x emissions, and reduced exhaust temperatures. Miller cycle without EGR at constant λ demonstrated LIVC strategies effectively reduce engine-out NO x emissions by up to 35%. However, Miller cycle associated with very aggressive LIVC timings led to fuel consumption penalties due to increased pumping work and exhaust enthalpy. LIVC strategies allowed for increased charge dilution at the baseline NO x constraint of 3.2 g/kWh, resulting in significant fuel consumption benefits over the baseline case without compromising exhaust temperatures or peak cylinder pressures. As Miller cycle implementation was shown to affect the boundary conditions dictating η TC , the LIVC and conventional IVC cases were studied at an equivalent η TC point representative of high boost operation. With high boost, LIVC yielded reduced NO x emissions, reduced peak cylinder pressures, and elevated exhaust temperatures compared to the conventional IVC case without compromising fuel consumption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Non-Fouling, Low Cost Electrolytic Coagulation & Disinfection for Treating Flowback and Produced Water for Reuse

Executive Summary This Final Report is composed of two major sections. The first section presents experimental results obtained in the laboratory during Budget Period 1. The second section presents results from the field test conducted during Budget Period 2. Laboratory Results This research investigated a novel electrochemical process for producing a ferric iron coagulant for use in treating flowback and produced water from hydraulic fracturing and oil production operations. The treatment system improves the effectiveness and lowers the cost of coagulation processes using Fe3+ as the coagulant. The electrolytic coagulant generation (ECG) system uses an electrochemical cell to produce acid and base from oilfield brine solutions. The acid is used to dissolve scrap iron to provide a Fe3+ coagulating agent. The base is used to neutralize the treated water. Compared to conventional electrocoagulation (EC), the main advantage of the ECG system is an order of magnitude lower cost for the source of iron. The second advantage over conventional EC is that it can deliver Fe3+ doses greater than 1 mM, since it is not limited by the amount of dissolved oxygen in the water required to oxidize ferrous to ferric iron. The capital costs for conventional EC and the ECG system are similar, but the operational costs for the ECG system are an order of magnitude lower than conventional EC. The combined costs for iron and electrical energy for treating 1 m3 of FPW with 1 mM Fe3+ is estimated to be $0.87 for conventional EC, and $0.087 for the ECG system. The estimated all-in cost for treating FPW with a 2 mM Fe3+ dose is $0.73/m3 ($0.12/bbl). Field Test Results The field test was conducted at the Paul Foster Central Tank Battery (CTB) in Lea County, New Mexico from November 10, 2022 through December 15, 2022. The feed water to the system was produced water from the Tatanka 1H formation. After approximately two weeks of testing, the initial batch of produced water had been treated and no untreated produced water was available. Thus, after this time, the feed water to the system consisted of previously treated water (i.e., recycled water). The recycled water had nearly all colloidal particles removed, and had a much lower alkalinity due to precipitation of carbonate minerals during the first pass through the system. Although the recycled water was not an ideal test solution due to its low particulate concentrations, its lower alkalinity did allow us to identify the main problem with the treatment system. The main problem with the treatment system was caused by the high alkalinity of the initial feed water (5.4 meq/L) that consumed a significant fraction of the electrochemically generated acid. This resulted in pH values exiting the iron contact tank that were too high to dissolve enough iron to effectively treat the produced water. Tests performed with recycled water with lower alkalinity did not have this issue, and dissolved iron concentrations greater than 20 mM could be achieved. One consistent observation was that effluent water from the iron contact tank was always free of particulates, even when fed with circumneutral solutions. This suggests that there is no need to dissolve high concentrations of iron if all the water to be treated is passed through the scrap iron canister. In this case, dissolved O2 and hypochlorous acid can promote sufficient iron corrosion to provide an effective coagulating agent – even in neutral pH water. This solves the problem resulting from highly alkaline produced water. Modifications to the design of the treatment system were made based on the field test results. These modifications will add minimal additional cost, and were tested in bench-scale laboratory experiments. In short-term testing, the modified treatment process was able to remove colloidal FeS particulates to levels below detection. Long-term, steady state testing will be required to determine whether the modified process is suitable for commercial treatment systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The EGS Collab Project – Stimulations at Two Depths

The EGS Collab project, supported by the US Department of Energy, is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to address challenges in implementing enhanced geothermal systems (EGS). Data and observations from the field tests are compared to simulations to understand processes and build confidence in numerical modeling of the processes. We have completed Experiment 1 (of 3), which examined hydraulic fracturing in a well-characterized underground fractured phyllite test bed at a depth of approximately 1.5 km at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. Testbed characterization included fracture mapping, borehole acoustic and optical televiewers, full waveform sonic, conductivity, resistivity, temperature, campaign p- and s-wave investigations and electrical resistance tomography. Borehole geophysical techniques including passive seismic, continuous active source seismic monitoring, electrical resistance tomography, fiber-based distributed strain, distributed temperature, and distributed acoustic monitoring, were used to carefully monitor stimulation events and flow tests. More than a dozen stimulations and nearly one year of flow tests were performed. Quality data and detailed observations were collected and analyzed during stimulation and water flow tests using ambient temperature and chilled water. We achieved adaptive control of the tests using real-time monitoring and rapid dissemination of data and near-real-time simulation. More detailed numerical simulation was performed to answer key experimental design questions, forecast fracture propagation trajectories and extents, and analyze and evaluate results. Data are freely available from the Geothermal Data Repository. Experiment 2 examines the potential for hydraulic shearing in amphibolite at a depth of about 1.25 km at SURF. This site has a different set of stress and fracture conditions than Experiment 1. The Experiment 2 testbed consists of nine subhorizontal boreholes configured in two fans of two boreholes which surround the testbed and contain grouted-in electrical resistance tomography, seismic sensors, active seismic sources and distributed fiber sensors. A “five-spot” set of test wells that extends from a custom mined alcove includes an injection well and four production/monitoring wells. The testbed was characterized geophysically and hydrologically, and three stimulations have been performed using the Step-Rate Injection Method for Fracture In-Situ Properties (SIMFIP) tool to measure strains, and a new strain quantifying tool (downhole robotic strain analysis tool -DORSA) was deployed in a monitoring hole during stimulation. Real-time data were broadcast during stimulations to allow real-time response to arising issues.

EGS Collab, Enhanced Geothermal Systems, EGS, fiel↗
DOE OSTI2023

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. 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, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗
DOE OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2021

Modeling of the High Temperature Test Facility Using RELAP-7

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is an electrically heated, helium cooled, experimental facility. The HTTF was modeled using RELAP-7 in both 1-D and 3-D for validation. Due to core symmetry, 1/6th of the HTTF core was modeled in 3-D using the coupled heat conduction and forced convection capabilities of RELAP-7. This served as a validation case for RELAP-7 and its capabilities to model advanced nuclear technologies such as high temperature gas-cooled reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

NMDQi Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST)

The Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST) is developed to support the thermal-hydraulic aspects for the design of optimal irradiation vehicles for the Nuclear Materials Discovery and Qualification initiative (NMDQi). TOAST is currently a MATLAB-based tool that utilizes a simplified steady-state analytical model where a radial thermal resistance circuit is utilized to account for conduction through up to 2 capsules, a gas gap, a thermal bond, a sample cladding, and a sample, as well as the convection from a water coolant outside the capsules. The geometry is discretized axially along the user-defined height of a basket with a user-defined number of points/nodes, essentially turning TOAST into a semi-2D heat transfer analysis utility for cylindrical irradiation vehicles with practically any configuration, solving for temperatures radially at multiple axial locations. TOAST utilizes a Graphical User Interface (GUI) in-which a user can define the geometrical layout, the materials utilized for each component, the coolant characteristics, and the required solution. The user can define the geometry by inputting the diameters, thicknesses, and heights for each component, whereas the materials are defined via the user-provided constant or variable thermal conductivities. The user can select the coolant's inlet temperature, pressure, and the pressure drop across the height of the problem (which is used to calculate the velocity of the flow). Finally, the user can choose to solve for a sample heat generation rate limit by inputting temperature limits for the sample and capsules or can choose to purely solve for the axial temperature distributions of each component by inputting a pre-selected heat generation rate. Either way, TOAST provides the user with axial temperature distributions of all the components including the coolant, and the heat generation rate limit as well as the maximum outlet coolant temperature. The user can also choose to do one of 6 sensitivity analyses in TOAST. The sensitivity analyses yield plots of the sensitivity of the heat generation rate limit, the maximum coolant temperature, and the pressure limit for the annular components due to perturbations in 1-2 unknown variables based on the selected sensitivity analysis. Benchmarks between computations in TOAST and equivalent 2D axis symmetric ABAQUS models are presented, showing that TOAST results are within less than 3% of ABAQUS results in most cases, and a maximum of 8% difference in some cases. The benchmarks also revealed that this uncertainty is tied to the selection of an appropriate Nusselt number correlation and appropriate thermal conductivities, which the user can do from the GUI. Regardless, TOAST is demonstrated as a computationally efficient, highly accessible, and accurate utility for optimization and scoping s of studies different irradiation vehicles.

42 ENGINEERING↗