Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Computational Fluid Dynamics Modeling of Solar Thermal Dry Reforming of Methane in a Parabolic Trough

Computational fluid dynamics simulations of solar-thermal dry reforming of methane using a parabolic trough configuration were performed. Parametric simulations of different combinations of gas flow rate, receiver tube emissivity, and geometric concentration ratio were conducted to determine configurations that could achieve the required catalyst temperatures of at least 700 °C to achieve high conversion of CH4 and CO2 to H2 and CO. Results showed that the concentration ratio of the parabolic trough collector had to be increased from ~70 to ~120 and the receiver-tube emissivity had to be reduced to ~0.2 to achieve bulk average catalyst temperatures of greater than 700 °C. Lower gas flow rates also reduced enthalpic heat losses and increased catalyst temperatures.

Ho, Clifford↗
DOE OSTI2021

High thermoelectric figure of merit of porous Si nanowires from 300 to 700 K.

Thermoelectrics operating at high temperature can cost-effectively convert waste heat and compete with other zero-carbon technologies. Among different high-temperature thermoelectrics materials, silicon nanowires possess the combined attributes of cost effectiveness and mature manufacturing infrastructures. Despite significant breakthroughs in silicon nanowires based thermoelectrics for waste heat conversion, the figure of merit (ZT) or operating temperature has remained low. Here, we report the synthesis of large-area, wafer-scale arrays of porous silicon nanowires with ultra-thin Si crystallite size of ~4 nm. Concurrent measurements of thermal conductivity (κ), electrical conductivity (σ), and Seebeck coefficient (S) on the same nanowire show a ZT of 0.71 at 700 K, which is more than ~18 times higher than bulk Si. This ZT value is more than two times higher than any nanostructured Si-based thermoelectrics reported in the literature at 700 K. Experimental data and theoretical modeling demonstrate that this work has the potential to achieve a ZT of ~1 at 1000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Solid-State Transformation of an Additive Manufactured Inconel 625 Alloy at 700 °C

Inconel 625, a nickel-based superalloy, has drawn much attention in the emerging field of additive manufacturing (AM) because of its excellent weldability and resistance to hot cracking. The extreme processing condition of AM often introduces enormous residual stress (hundreds of MPa to GPa) in the as-fabricated parts, which requires stress-relief heat treatment to remove or reduce the internal stresses. Typical residual stress heat treatment for AM Inconel 625, conducted at 800 °C or 870 °C, introduces a substantial precipitation of the δ phase, a deleterious intermetallic phase. In this work, we used synchrotron-based in situ scattering and diffraction methods and ex situ electron microscopy to investigate the solid-state transformation of an AM Inconel 625 at 700 °C. Our results show that while the δ phase still precipitates from the matrix at this temperature, its precipitation rate and size at a given time are both smaller when compared with their counterparts during typical heat treatment temperatures of 800 °C and 870 °C. A comparison with thermodynamic modeling predictions elucidates these experimental findings. Our work provides the rigorous microstructural kinetics data required to explore the feasibility of a promising lower-temperature stress-relief heat treatment for AM Inconel 625. The combined methodology is readily extendable to investigate the solid-state transformation of other AM alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetoelastic coupling and spin contributions to entropy and thermal transport in biferroic yttrium orthochromite

Direct engineering of material properties through exploitation of spin, phonon, and charge-coupled degrees of freedom is an active area of development in materials science. However, the relative contribution of the competing orders to controlling the desired behavior is challenging to decipher. In particular, the independent role of phonons, magnons, and electrons, quasiparticle coupling, and relative contributions to the phase transition free energy largely remain unexplored, especially for magnetic phase transitions. Here, we study the lattice and magnetic dynamics of biferroic yttrium orthochromite using Raman, infrared, and inelastic neutron spectroscopy techniques, supporting our experimental results with first-principles lattice dynamics and spin-wave simulations across the antiferromagnetic transition at T N ~ 138 K. Spectroscopy data and simulations together with the heat capacity (C p ) measurements, allow us to quantify individual entropic contributions from phonons (0.01 ± 0.01k B atom -1 ), dilational (0.03 ± 0.01k B atom -1 ), and magnons (0.11 ± 0.01k B atom -1 ) across TN. High-resolution phonon measurements conducted in a magnetic field show that anomalous T-dependence of phonon energies across TN originates from magnetoelastic coupling. Additionally, phonon scattering is primarily governed by the phonon–phonon coupling, with little contribution from magnon–phonon coupling, short-range spin correlations, or magnetostriction effects; a conclusion further supported by our thermal conductivity measurements conducted up to 14 T, and phenomenological modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Meteorological and Soil Data from Ecohydrology Sensor Towers at Pump House and Snodgrass Mountain in East River Watershed, Colorado, 2019-2025

This data package includes hourly meteorological and soil sensor data at eight ecohydrology monitoring sites in East River Watershed, Colorado as part of the Watershed Function Scientific Focus Area (WFSFA) research led by Lawrence Berkeley National Lab (LBNL). Four field sites were located on the hillslope of East River (ER) near Pump House (PH) at Mount Crested Butte (ER-PHS1 to 4), and the other four are in the Snodgrass Mountain (SG) area (SG-EHS5 to 8). In terms of vegetation cover, three sites are in montane grasslands (ER-PHS1, ER-PHS2, and SG-EHS5), three are below evergreen conifer canopy (ER-PHS3, SG-EHS6, and SG-EHS7), and two are below deciduous aspen canopy (ER-PHS4 and SG-EHS8). The monitoring period began in October 2019 at the East River sites, in October 2020 at SG-EHS5 and SG-EHS6, and in October 2021 at SG-EHS7 and SG-EHS8. In September 2024, all four East River sites were fully retired. The four Snodgrass Mountain sites remain active. Each site is equipped with a comprehensive suite of meteorological sensors on a tripod and soil sensors that measure weather, energy fluxes, and soil variables. This data package includes measurements from ten different types of sensors and up to thirteen individual sensors per site, including (1) a weather station (measurement height ranges from 2.8~3.8 meters (m) above ground), (2) a quantum sensor for photosynthetic active radiation (PAR) (2.4~3.3m), (3) a net radiometer (1.7~2.1m), (4) an infrared radiometer (1.6~2.2m), (5) a sonic distance sensor (1.5~1.9m), (6) a soil carbon dioxide (CO2) flux chamber (0m), (7) a soil heat flux plate (-0.05m below ground), (8) a soil oxygen sensor (-0.3m), (9) a soil water potential sensor (-0.3m), and (10) soil water content sensors at 3~4 depths (-1.15 ~ -0.1m). A total of twenty-three variables is reported in this data package, including (1) atmospheric variables: air temperature (TA), atmospheric pressure (PA), vapor pressure (VP), and vapor pressure deficit (VPD), (2) precipitation variables: rain precipitation (P) and snow depth (D_SNOW), (3) energy fluxes variables: four-component net radiation (NETRAD) (shortwave/longwave incoming/outgoing radiation, SW_IN, SW_OUT, LW_IN, LW_OUT), photosynthetic photon flux density (PPFD), and soil heat flux (G), (4) soil variables: soil water content (SWC), soil water potential (SWP), soil temperature (TS), soil bulk electrical conductivity (COND_SOIL), and soil gaseous oxygen concentration (O2_SOIL), (5) wind variables: two-dimensional wind speed (WS), gust speed (WS_MAX), and wind direction (WD), and (6) surface variables: surface infrared temperature (T_CANOPY) and soil CO2 flux (CO2_SOIL). Please see the Methods section for data processing and QA/QC steps taken to generate the hourly datasets. The following files are included in this data package (notes on version: v{x}-{y}, where x is the metadata version, and y is the data version, when applicable): (1) “metadata_site_v{x}-{y}.csv” - a site metadata file that summarizes location information of all sites, including site ID, description, coordinates, timeframe, elevation, and vegetation cover, (2) “metadata_instrument_v{x}-{y}.csv” - an instrument metadata file that summarizes sensor information of all sites, including sensor manufacturer and model, measurement height, and sampling and averaging interval of all variables, (3) "data_{SITE_ID}_v{x}-{y}.csv" - eight data files that contain hourly data of each site indicated by {SITE_ID} in the filename, (4) “/figure/data_{SITE_ID}_v{x}-{y}.png" - eight figures that help visualize data of each site indicated by {SITE_ID} in the filename, (5) “/photo/*” - photos of each site indicated by {SITE_ID} in the filename, and (6) four file level metadata (flmd.csv) and data dictionary (*_dd.csv) files that summarize file, header, column, and variable information of all files. Notes: (1) Measurement height: Each variable name is followed by conventional positional qualifiers “H_V_R”, where H indicates the relative horizontal positions of that specific variable, V the vertical positions, and R the replicates. In this data package, only the vertical qualifier V varies, and V increases from the highest vertical position (V=1) to the lowest. Variables with the same qualifier are not necessarily measured by the same sensor, and the same variable with the same qualifier across different sites are not necessarily measured at the same height. Please refer to “metadata_instrument.csv” for the sensor information and measurement heights, and whether a variable is measured below the canopy. (2) Variable availability: Snow depth is not available at ER-PHS3 and SG-EHS7. SWC, soil temperature, and soil bulk EC at the deepest depth (<-1m) are not available at SG-EHS6 and SG-EHS7. The missing value code for numeric variables is -9999, except for SWP. For SWP, the missing value code is +9999, because SWP values are negative. (3) Sampling frequency: Please refer to “metadata_instrument.csv” for the increase of sampling frequency of some variables from 30-min to 1-min at ER-PHS1 to 4 in July 2020. (4) Sensors: While the methods of each sensor are not detailed, all sensors are commercially available, and their methods can be found in their manuals. Please refer to “metadata_instrument.csv” for the sensor manufacturer and model information. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Creep Deformation and Damage Mechanisms in an Advanced High-Temperature Additively Manufactured Nickel-Base Superalloy

Abstract This research investigates the processing–structure–properties–performance relationship in a novel nickel-base superalloy, ABD ® -900AM, designed for extreme environments. Specifically tailored for additive manufacturing (AM), ABD ® -900AM maintains mechanical integrity at high temperatures and is comparable to other nickel-based superalloys with a 30–40% gamma-prime volume fraction. A comprehensive study was conducted using laser-beam powder bed fusion and electron-beam powder bed fusion methods. Factors such as heat treatment, porosity, build orientation, and hot isostatic pressing were evaluated to understand their effects on microstructure and mechanical performance. Microstructural characterization revealed significant differences in grain size and orientation across build processes and heat treatments. High-temperature mechanical testing indicated that grain size, heat treatment, and orientation significantly influence creep behavior. A super-solvus heat treatment led to recrystallization and grain growth, significantly improving creep properties compared to a near-solvus heat treatment. Various creep mechanisms were identified across different conditions, and creep rupture models were developed for each build process. Post-test microstructural analysis showed grain boundary damage, with differences in creep cavitation morphology under varying stress conditions. It was shown that MC carbides grow at the expense of gamma-prime near grain boundaries, leading to precipitate-free zones in specimens tested at higher temperatures. This study fills a significant gap in fundamental research by offering a deeper insight into the high-temperature mechanical behavior of additively manufactured nickel-base superalloys. It also explores critical research questions regarding the role of carbides and the significance of heat treatment. The insights gained enhance confidence in the industry adoption of ABD ® -900AM and similar alloys for high-temperature applications, bridging the knowledge gap and supporting the development of reliable AM processes for extreme environments.

Bridges, Alex (ORCID:000000030338759X)↗
DOE OSTI2023

Evaluation of ITER divertor shunts as a synthetic diagnostic for detachment control

Reliable diagnostics that measure the detached state of the ITER divertor plasma will be necessary to control heat flux to the divertor targets during steady state, burning plasma operation. This paper conducts an initial exploration into the feasibility of the divertor shunt diagnostic as a lightweight, robust, and real-time detachment sensor. This diagnostic is a set of shunt lead pairs that measure the voltage drop along the divertor cassette body, from which the plasma scrape-off layer (SOL) current is calculated. Using SOLPS-ITER simulations for control-relevant ITER plasma scenarios, the thermoelectric current magnitude along the SOL is shown to decrease significantly with the onset of partial detachment at the outer divertor target. Electromagnetic modelling of a simplified divertor cassette is used to develop a control-oriented inductance-resistance circuit model, from which SOL currents can be calculated from shunt pair voltage measurements. The sensitivity and frequency-response of the resulting system indicates that the diagnostic will accurately measure SOL thermoelectric currents during ITER operation. These currents will be a good measure of the detached state of the divertor plasma, making the divertor shunt diagnostic a potentially extremely valuable and physically robust sensor for real-time detachment control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Gas Generation Models for Thermal Runaway in Lithium Ion Batteries

A numerical model is implemented to study cell venting, internal pressure, and gas-phase dynamics behavior of Li-ion cells undergoing thermal runaway. A k-ɛ Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer’s equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. A series of computational fluid dynamics (CFD) simulations are conducted on a single cell at various states-of-charge (100%, 50%, 25%) to demonstrate detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is shown that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

Kim, Jinyong↗
DOE OSTI2023

Misr: Miniature Ignition Screening Rapid Compression Machine for Kinetic Measurements of Novel Fuels (Final Report)

High throughput ignition screening would be useful for evaluating fuels before scaled production. One method for doing this involves generating a map of ignition delay times over different conditions for each fuel blend. Constant condition ignition delay time maps can be used to estimate ignition for engines in prescribed thermodynamic trajectories. A Miniaturized Ignition Screening Rapid Compression Machine (MISR) has been developed and built at the University of Illinois Chicago (UIC) intending to generate those constant condition ignition delay time maps for low temperature chemistry relevant conditions. Some design features include a bellows-actuated driver section, which is lighter compared to other RCM designs, an adjustable brake section, and full optical access. The setup was designed to achieve reaction conditions (T > 600K and P < 100 bar with repetition rates of about 8 seconds per experiment). These conditions were meant to be relevant to low-to-intermediate temperature auto-ignition and for generation of data for developing relevant chemical kinetics models or ignition maps. Because of the compression characteristics and the heat loss characteristics of the MISR, the apparatus currently tests higher temperatures over shorter test times than what would be desirable, however there are some pathways to improved compression characteristics as detailed in this report. A detailed description of the MISR design and its operating mechanism is discussed. The team faced several challenges and setbacks in the construction and commissioning of the MISR, and those are highlighted, along with the modifications that were performed to get the MISR operational. To demonstrate the performance characteristics of the MISR, limited ignition delay experiments were conducted for ethanol fuel at different equivalence ratios (Φ), and those results were usable, but again test a higher temperature and shorter ignition delay time than initially intended.

02 PETROLEUM↗
DOE OSTI2022

Datasets of a Multizone Office Building under Different HVAC System Operation Scenarios

Abstract This study provides an open-source dataset of the measured weather data, building indoor data, and system data under the different test settings. The test building is the two-story Flexible Research Platform building at the US Department of Energy’s Oak Ridge National Laboratory, in Oak Ridge, Tennessee. Four heating tests and three cooling tests were conducted. The 1-min interval of weather, building indoor data, and system data from each test setting are provided. Actual weather data were collected from a weather station installed on the roof. This paper describes information on the test building and installed sensors, data collection method, and data validation. The provided dataset can be employed to understand HVAC system conditions and building indoor conditions under different HVAC system operations and the performance of building envelope without HVAC system operation using free-floating test data. Additionally, it can be used for empirical validation of the building energy modelling engine.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Transient magnetic field diffusion considerations relevant to magnetically assisted indirect drive inertial confinement fusion

We report application of a magnetic field to an indirect drive inertial confinement fusion target requires diffusion of the field through the high-Z and electrically conducting Hohlraum. The onset of the external field generates eddy currents in the Hohlraum wall that result in (1) a reduction of the peak field at the capsule, (2) heating of the Hohlraum wall through Ohmic dissipation, and (3) wall movement due to the inward force from the eddy current interacting with the field. Heating of the wall causes an increase in blackbody radiation which can preheat the capsule and frozen deuterium–tritium fuel, while wall motion leads to potential misalignment of the lasers at the Hohlraum wall. Limiting these detrimental effects sets requirements on the tolerable magnitude of each effect. We present a nonlinear model for B-field diffusion through an infinitely long thin-walled cylinder with a temperature dependent resistivity, to show that a 15 μm thick wall of pure gold fails to meet these requirements. A new Hohlraum material made from an alloy of Au and Ta has a measured resistivity of ≥60 times that of Au and is shown with the nonlinear model to meet the requirements for magnetization. We compare the nonlinear model to simulations of the actual Hohlraum target using a finite element code which includes temperature-dependent Hohlraum resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessment of RELAP5-3D Code with Molten Salt Heat Transfer Experiments

The accurate thermal-hydraulic assessment of thermal storage systems is crucial for addressing the integrity and performance of the storage system design. This is particularly important for thermal storage systems using molten salts as a heat transport and storage medium, as the chemical corrosion and erosion by the high-temperature molten salts can cause serious damage to the structural components of the storage system. To understand these effects for TerraPower’s Natrium® Demonstration Reactor design, the system thermal-hydraulic analysis code, RELAP5-3D, is utilized to analyze the thermal storage system. This study evaluates two heat transfer correlations implemented in the RELAP5-3D code—the Dittus-Boelter and Gnielinski correlations—using selected benchmark cases from two molten salt heat transfer experiments conducted at Xi’an Jiao Tong University and the German Aerospace Center. The Nusselt numbers calculated using the Gnielinski correlation agree with the experimental data within ±5% of the relative deviations at 300°C and 400°C. With the increasing Reynolds numbers, the Dittus-Boelter correlation underestimates the Nusselt numbers by up to 22% compared to the Gnielinski correlation. Comparison of the RELAP5-3D calculation results with experimental data at a high temperature of 550°C revealed that the thermo-chemical characteristics of solar salt at high temperatures surpassing the decomposition temperatures of nitrate salts reduce the accuracy of heat transfer model of the RELAP5-3D code. Specifically, it has been observed that the temperatures of the wall surface and liquid near the wall can surpass the decomposition temperatures, even though the bulk temperature of molten salts remains significantly lower than these decomposition temperatures. In summary, this study recommends the use of the Gnielinski correlation for the heat transfer analysis of molten salt, rather than the Dittus-Boelter correlation.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Advancing Mass Timber Buildings: Novel Methods Improve Thermal Assessment and Material Use

For nearly a century, thermal demand calculations for buildings have relied on simplified models developed to match the technical constraints of their era. The first standards, introduced in Germany and Austria in 1929, established climate zones and material conductivity coefficients that, with only incremental updates, still underpin many current assessments. Yet, methods such as Hot box testing, originally designed for lightweight insulation, continue to be applied for mass timber buildings, overlooking thermodynamic characteristics confer real-world advantages. Recent research at Oak Ridge National Laboratory incorporates updated methodologies, aligned with ASHRAE Standard 55 (ASHRAE, 2023) accounting for factors such as thermal inertia, inner surface temperatures, emissivity, solar gains, and dynamic outdoor conditions. These factors better reflect observed heating and cooling loads and highlight opportunities for efficient use of materials in mass timber construction. This work provides a framework for designing comfortable, resilient, and resource-efficient buildings while aligning with performance expectations in energy codes.

Pickett, Robert [International Mass Timber Allianc↗
DOE OSTI2020

Comparison of Model Predictions and Performance Test Data for a Prototype Thermal Energy Storage Module

Although model predictions of thermal energy storage (TES) performance have been explored in previous investigations, relevant test data that enable experimental validation of performance models have been limited. This is particularly true for high-performance TES designs that facilitate fast input and extraction of energy. In this paper, we present a summary of experimental tests of a high-performance TES unit using lithium nitrate trihydrate phase change material as a storage medium. Performance data are presented for complete dual-mode cycles consisting of extraction (melting) followed by charging (freezing). These tests simulate the cyclic operation of a TES unit for asynchronous cooling in a variety of applications. Finally, the model analysis is found to agree reasonably well, within 10%, with the experimental data except for conditions very near the initiation of freezing, a consequence of subcooling that is required to initiate solidification.

25 ENERGY STORAGE↗
DOE OSTI2024

Simulation of Helium Flow Visualization Apparatus for Studies of Blanket Cooling in Fusion Reactors

Flow visualization is essential to understanding helium cooling performance. This article investigates helium flow visualization in an apparatus to support the design of the blanket first wall for a fusion reactor. Helium’s safety advantages make it an attractive coolant, but effective cooling and flow visualization remain challenging. Sophisticated simulations are conducted in different test sections to address three key areas in the design of a test apparatus and cooling channel enhancements: the influence of viewing glass windows on flow patterns, the impact of increased heating, and the effectiveness of cooling structures like baffles. Viewing glass windows introduce flow recirculation region bifurcation and flow asymmetry, affecting flow patterns and necessitating careful analysis in future experiments. Increased heating results in helium flow detachment from the heated surface, leading to hot spots. Cooling structures, particularly baffles, prove effective in maintaining consistent attachment to the heated surface, improving the heat transfer performance. Further, this study also examines a tunable parameter in the turbulence model, highlighting the importance of accurate model tuning for future fusion reactor cooling designs. With an imminent helium flow visualization facility, these simulations will be used to optimize cooling structures on the heated wall to improve flow attachment and heat transfer efficiency. This work serves as a first investigation of the helium flow visualization apparatus for blanket cooling enhancement in fusion reactor design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermal Characterization and Exhumation of Northwest San Juan Basin Area, NM

The San Juan Basin (SJB), located in southwestern Colorado northwestern New Mexico, containing Paleozoic through middle Cenozoic strata, formed as a partitioned basin approximately 80 Ma in response to the Laramide orogeny. The SJB is a commercially mature, petroleum-producing basin, and currently being explored for CCUS and geothermal resources. Modern heat flow within the SJB is spatially variable with higher geothermal gradients in its northern and eastern portions. The temporal history of the basin thermal history is essential for understanding capacity for carbon storage and geothermal exploration. We present thermochronometric analyses of two surface and two subsurface samples to constrain spatial-temporal thermal evolution of the northwestern San Juan Basin. Apatite (U-Th)/He thermochronometric analyses were conducted on Cliff House and Kirtland Formation outcrops and subsurface samples of the Ojo Alamo and Pictured Cliffs Formations from a DOE funded, CCUS project pilot-well. In addition, vitrinite reflectance data and 1D basin modelling constrain possible time-temperature pathways of the samples in an area where no previous studies have constrained the uplift/exhumation of the Hogback Monocline using thermochronometry. Modeling suggests that elevated subsurface temperatures developed simultaneously with regional late Oligiocene volcanism in the San Juan Volcanic field and timing of maximum burial in the SJB. Late Miocene to Pliocene cooling/exhumation through the apatite (U-Th)/He partial retention zone support stratigraphic evidence of the Colorado Plateau uplift, possibly regional epeirogenic uplift driven by mantle processes.

02 PETROLEUM↗
DOE OSTI2022

Flow in additively manufactured super-rough channels

Metal additive manufacturing has enabled geometrically complex internal cooling channels for turbine and heat exchanger applications, but the process gives rise to large-scale roughness whose size is comparable to the channel height (which is 500 $\mathrm {\mu }$ m). These super-rough channels pose previously unseen challenges for experimental measurements, data interpretation and roughness modelling. First, it is not clear if measurements at a particular streamwise and spanwise location still provide accurate representation of the mean (time- and plane-averaged) flow. Second, we do not know if the logarithmic layer survives. Third, it is unknown how well previously developed rough-wall models work for these large-scale roughnesses. To answer the above practical questions, we conduct direct numerical simulations of flow in additively manufactured super-rough channels. Three rough surfaces are considered, all of which are obtained from computed tomography scans of additively manufactured surfaces. The roughness’ trough to peak sizes are 0.1 $h$ , 0.3 $h$ and 0.8 $h$ , respectively, where $h$ is the intended half-channel height. Each rough surface is placed opposite a smooth wall and the other two rough surfaces, leading to six rough-wall channel configurations. Two Reynolds numbers are considered, namely $Re_\tau =180$ and $Re_\tau =395$ . We show first that measurements at one streamwise and spanwise location are insufficient due to strong mean flow inhomogeneity across the entire channel, second that the logarithmic law of the wall survives despite the mean flow inhomogeneity and third that the established roughness sheltering model remains accurate.

Altland, Samuel↗