Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat map”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Improving the creation of SiV centers in diamond via sub-μs pulsed annealing treatment

Silicon-vacancy (SiV) centers in diamond are emerging as promising quantum emitters in applications such as quantum communication and quantum information processing. Here, we demonstrate a sub-μs pulsed annealing treatment that dramatically increases the photoluminescence of SiV centers in diamond. Using a silane-functionalized adamantane precursor and a laser- heated diamond anvil cell, the temperature and energy conditions required to form SiV centers in diamond were mapped out via an optical thermometry system with an accuracy of ±50 K and a 1 μs temporal resolution. Annealing scheme studies reveal that pulsed annealing can obviously minimize the migration of SiV centers out of the diamond lattice, and a 2.5-fold increase in the number of emitting centers was achieved using a series of 200-ns pulses at a 50 kHz repetition rate via acousto-optic modulation. Our study provides a novel pulsed annealing treatment approach to improve the efficiency of the creation of SiV centers in diamond.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Axon-like active signal transmission

Any electrical signal propagating in a metallic conductor loses amplitude due to the natural resistance of the metal. Compensating for such losses presently requires repeatedly breaking the conductor and interposing amplifiers that consume and regenerate the signal. This century-old primitive severely constrains the design and performance of modern interconnect-dense chips. Here we present a fundamentally different primitive based on semi-stable edge of chaos (EOC), a long-theorized but experimentally elusive regime that underlies active (self-amplifying) transmission in biological axons. By electrically accessing the spin crossover in LaCoO 3 , we isolate semi-stable EOC, characterized by small-signal negative resistance and amplification of perturbations. In a metallic line atop a medium biased at EOC, a signal input at one end exits the other end amplified, without passing through a separate amplifying component. While superficially resembling superconductivity, active transmission offers controllably amplified time-varying small-signal propagation at normal temperature and pressure, but requires an electrically energized EOC medium. Operando thermal mapping reveals the mechanism of amplification—bias energy of the EOC medium, instead of fully dissipating as heat, is partly used to amplify signals in the metallic line, thereby enabling spatially continuous active transmission, which could transform the design and performance of complex electronic chips.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Low-Temperature Geothermal Play Fairway Analysis for the Denver Basin

This project is part of a national initiative to showcase the benefits of incorporating low-temperature geothermal resource assessment into the deployment of geothermal heating, combined heat and power (CHP), and geothermal direct-use (GDU) technologies. The initiative was established to accelerate the country's decarbonization efforts by identifying potential for low-temperature geothermal resource utilization (<150 degrees C, e.g., CHP and GDU) in selected sedimentary basins with numerous population centers. The play fairway analysis (PFA) methodologies in this study were adapted from previous PFA investigations of sedimentary basin geothermal play types (SBGPTs) that evaluated the potential for low-temperature resources (<150 degrees C). Workflows, relevant datasets, a new Python library, and common and composite geological criteria maps are utilized to develop low-temperature geothermal resource favorability maps for the Denver Basin, a sedimentary basin spanning Colorado, Nebraska, and Wyoming. The replication of these methodologies in other SBGPTs can evaluate potential for low-temperature resources. To facilitate future assessment of low-temperature geothermal resources in SBGPTs, this project provides PFA workflows, data, tools, and favorability maps that will ultimately support the utilization of low-temperature geothermal resources in sedimentary basins.

combined heat and power↗
DOE OSTI2023

High-Temperature Air Stability of Electrodeposited Copper Cobalt Oxide and Copper Manganese Oxide Absorber Coatings for Concentrating Solar Power

The solar absorptance of a thermal receiver surface significantly affects the photothermal conversion efficiency of concentrating solar power (CSP) plants. The stability of solar receivers and coatings applied to improve solar absorptance at high temperatures is critical for efficient and long-term operation. This study analyzes the high-temperature air stability of multiscale fractal textured electrodeposited copper manganese oxide and copper cobalt oxide solar receiver coatings over a range of temperatures. The study also investigates the effect of high temperature on phase stability, chemical and morphological evolution, and the optical properties of the two solar absorber materials. A computational model is developed to determine the absorptance of the textured absorbers and is shown to match the experimental measurements closely over a range of heat treatment temperatures from 450 °C to 850 °C and durations up to 100 h. Degradation curves and coating stability maps are developed from the study. Copper manganese oxide and copper cobalt oxide coatings exhibit good thermal stability up to ~815 °C and 630 °C, respectively, with minimal changes in their optical and structural properties with thermal efficiency greater than 90%. Furthermore, the study offers valuable insights into the high-temperature annealing effects on the optical and structural properties of multiscale copper alloyed coatings fabricated by electrodeposition and presents a means of designing the optimal temperature range for maximizing absorber efficiency in practical applications.

14 SOLAR ENERGY↗
DOE OSTI2021

Report on Mowry Shale Thermal Maturity Mapping in the Powder River Basin, Wyoming

The Enhanced Oil Recovery Institute (EORI) has completed a map of thermal maturity of the Mowry Shale in the Powder River Basin (PRB). Thermal maturity, the degree to which the total organic carbon in a formation has been transformed from kerogen to producible hydrocarbons due to heat and pressure, is an important measure for not just the quality of a source rock, but also helps for delineating areas more favorable for unconventional drilling. There is a dearth of published Mowry maturity maps available for the PRB and this publication endeavors to provide a detailed map based on a large, public data set covering the Wyoming portion of the PRB. Additional data, currently held confidential by operators in the basin, would be of benefit to future iterations of this map.

58 GEOSCIENCES↗
DOE OSTI2023

Chapter 6: Abuse Response of Batteries Subjected to Mechanical Impact

Electrochemical and thermal models to simulate nominal performance and abuse response of lithium-ion cells and batteries have been reported widely in the literature. Studies on mechanical failure of cell components and how such events interact with the electrochemical and thermal response are relatively less common. This chapter outlines a framework developed under the Computer Aided Engineering for Batteries program to couple failure modes resulting from external mechanical loading to the onset and propagation of electrochemical and thermal events that follow. Starting with a scalable approach to implement failure criteria based on thermal, mechanical, and electrochemical thresholds, we highlight the practical importance of these models using case studies at the cell and module level. The chapter also highlights a few gaps in our understanding of the comprehensive response of batteries subjected to mechanical crash events, the stochastic nature of some of these failure events, and our approach to build safety maps that help improve robustness of battery design by capturing the sensitivity of some key design parameters to heat generation rates under different mitigation strategies.

abuse simulations↗
DOE OSTI2026

Extension of Clad Damage Propagation Model for Fission Gas Dispersal and Two-Phase Flow Effects in MOOSE SubChannel Module

This report presents an extension of the Clad Damage Propagation (CDAP) model implemented in the MOOSE SubChannel Module (SCM) to capture post-failure fission-gas dispersal and two-phase flow effects in sodium-cooled fast reactor assemblies. The extended model tracks discharged gas axially and radially, computes channel-averaged flow quality and void fraction using a Lockhart–Martinelli framework, evaluates two-phase frictional pressure-drop multipliers, determines inlet mass-flow degradation under fixed core pressures, and applies an intensified-void-based heat-transfer degradation to affected fuel pins. Radial plume expansion is parameterized using mineral-oil jet experiments mapped to sodium conditions via Reynolds–Weber similarity. Implementation details are documented, along with the new methods and user inputs needed to control plume mapping and two-phase behavior. Demonstration simulations for 19- and 37-pin bundles show that breach size and inlet velocity strongly influence propagation potential: small breaches (≤0.5 mm) produce limited degradation while larger breaches (~1 mm) can drive oscillatory temperature spikes and enhanced failure propagation, especially at higher velocities. These results demonstrate that the extended CDAP model provides a more complete framework for quantifying cladding damage propagation and evaluating propagation potential in transient scenarios. The approach remains computationally efficient, consistent with subchannel-level analysis, yet incorporates sufficient physics to bridge localized post-failure effects with bundle- and assembly-scale degradation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

A Comprehensive Comparison of Building Height Products at a Global Level

The built environment continues to expand due to increasing human populations with some projections indicating additional growth of 2.5 billion persons by 2050. Buildings are a key component of the built environment and information on building height and urban morphology in general can help to better understand and address a host of development and urban planning issues including energy consumption, urbanization, population dynamics, urban heat islands, environmental degradation, and others. To date, several approaches have been developed to address the need for mapping building height at scale. These approaches can be grouped into two main categories: 1) directly measure height from remote sensing products or 2) indirectly estimate height from footprint morphology features at the building-by-building level. Due to a number of these global height datasets now being available for consumption and use across domains, there is a critical need to understand the differences displayed between these products and how the use of these may influence the outcome when applied. In this analysis, we compare and evaluate 5 global building height datasets: Global Human Settlement, World Settlement Footprint 3D, Global Building Morphology dataset, three-dimensional Global Building Footprints, and an approach developed at Oak Ridge National Laboratory (ORNL). Our analysis shows that, across 5 cities at a global scale with high-quality building height reference data, the approach developed by ORNL of inferring height at the building-by-building level is more accurate and consistent. This was evidenced by the ORNL product improving upon the baseline in all 5 cities, the only product to do so. These findings will aid researchers in choosing a global building height product that best meets their needs and provides global availability, consistency, and accuracy.

42 ENGINEERING↗
DOE OSTI2025

Imaging Bragg Edge Analysis TooLs for Engineering Structures (iBeatles)

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) provides pulsed neutrons with energies varying from epithermal to cold. In preparation for VENUS, the neutron imaging beamline to be located at beam port 10, we have performed a series of experiments focused on wavelength-dependent radiography and computed tomography for a broad range of applications, from materials science to biological tissues.One of the time-of-flight (TOF) techniques that is of interest to the scientific community is the 2-dimensional mapping of phases and average crystalline plane orientation in samples both ex-situ and during applied stresses such as tensile loading and heating. This technique is known as Bragg edgeimaging and relies on the identification of changes of transmission values, fitting of the edge to measure its displacement, and thus identify the shift in lattice parameter due to stresses. One of the challenges of TOF imaging measurements is the amount of data and the inability to observe Bragg edge shifts in real time during an experiment. Thus, we have been focusing on creating a Python-based interface that allows fast data processing and instantaneous mapping and fitting of the Bragg edges, and their evolution through time. Python libraries and Jupyter notebooks have been implemented to facilitate decision making during an experiment. The advantage of the notebooks is the possibility to guide an experiment as they can quickly process and display Bragg edge data. These notebooks can be used independently, or can be combined in a Python Graphical User Interface (GUI) tool called iBeatles. This interface permits visualization and fitting of the Bragg edges, and ultimately back-projects the fitting results onto the radiographs to display a strain map. Assuming data collection has sufficient statistics, the strain mapping analysis can be performed on a pixel-by-pixel basis. This development is a step forward toward a better user experience at the future VENUS beamline in terms of live feedback and productivity. Analysis that used to take days of switching between different applications can now be done in minutes within the

Bilheux, JeanChristophe [Oak Ridge National Labora↗
DOE OSTI2023

Process-dependent anisotropic thermal conductivity of laser powder bed fusion AlSi10Mg: impact of microstructure and aluminum-silicon interfaces

Purpose AlSi10Mg alloy is commonly used in laser powder bed fusion due to its printability, relatively high thermal conductivity, low density and good mechanical properties. However, the thermal conductivity of as-built materials as a function of processing (energy density, laser power, laser scanning speed, support structure) and build orientation, are not well explored in the literature. This study aims to elucidate the relationship between processing, microstructure, and thermal conductivity. Design/methodology/approach The thermal conductivity of laser powder bed fusion (L-PBF) AlSi10Mg samples are investigated by the flash diffusivity and frequency domain thermoreflectance (FDTR) techniques. Thermal conductivities are linked to the microstructure of L-PBF AlSi10Mg, which changes with processing conditions. The through-plane exceeded the in-plane thermal conductivity for all energy densities. A co-located thermal conductivity map by frequency domain thermoreflectance (FDTR) and crystallographic grain orientation map by electron backscattered diffraction (EBSD) was used to investigate the effect of microstructure on thermal conductivity. Findings The highest through-plane thermal conductivity (136 ± 2 W/m-K) was achieved at 59 J/mm 3 and exceeded the values reported previously. The in-plane thermal conductivity peaked at 117 ± 2 W/m-K at 50 J/mm 3 . The trend of thermal conductivity reducing with energy density at similar porosity was primarily due to the reduced grain size producing more Al-Si interfaces that pose thermal resistance. At these interfaces, thermal energy must convert from electrons in the aluminum to phonons in the silicon. The co-located thermal conductivity and crystallographic grain orientation maps confirmed that larger colonies of columnar grains have higher thermal conductivity compared to smaller columnar grains. Practical implications The thermal properties of AlSi10Mg are crucial to heat transfer applications including additively manufactured heatsinks, cold plates, vapor chambers, heat pipes, enclosures and heat exchangers. Additionally, thermal-based nondestructive testing methods require these properties for applications such as defect detection and simulation of L-PBF processes. Industrial standards for L-PBF processes and components can use the data for thermal applications. Originality/value To the best of the authors’ knowledge, this paper is the first to make coupled thermal conductivity maps that were matched to microstructure for L-PBF AlSi10Mg aluminum alloy. This was achieved by a unique in-house thermal conductivity mapping setup and relating the data to local SEM EBSD maps. This provides the first conclusive proof that larger grain sizes can achieve higher thermal conductivity for this processing method and material system. This study also shows that control of the solidification can result in higher thermal conductivity. It was also the first to find that the build substrate (with or without support) has a large effect on thermal conductivity.

Engineering↗
DOE OSTI2026

Disentangling multiple gas kinematic drivers in the Perseus galaxy cluster

Galaxy clusters, the Universe’s largest halo structures, are filled with an X-ray-emitting gas with a temperature between 10 million and 100 million degrees. Their evolution is shaped by energetic processes such as feedback from supermassive black holes (SMBHs) and mergers with other cosmic structures. The imprints of these processes on gas kinematics remain largely unknown, restricting our understanding of energy conversion within clusters4. High-resolution spectral mapping with the X-Ray Imaging and Spectroscopy Mission (XRISM) observatory offers a way forward. Here we present XRISM kinematic measurements of the Perseus cluster, radially covering the extent of its cool core. We find direct evidence for at least two dominant drivers of gas motions operating on distinct physical scales: a small-scale driver in the inner approximately 60 kpc, probably associated with the SMBH feedback; and a large-scale driver in the outer core, powered by mergers. This finding suggests that, during the active phase, SMBH feedback drives gas motions, which, if fully dissipated into heat, could have a substantial role in offsetting radiative cooling losses in the Perseus core. Furthermore, our study underscores the necessity of kinematic mapping observations of extended sources to robustly characterize the velocity fields and their role in the evolution of massive halos. It further offers a kinematic diagnostic for SMBH feedback models.

Astronomy and AstroPhysics↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
DOE OSTI2021

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗