SEARCH · Engineering Papers

Results for “heating”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Comparative Analysis of Austenitized Steel Microstructures Due to Induction and Conventional Heating in High Magnetic Field

While magnetic field-assisted processing shows promise in improving alloys’ mechanical properties, insufficient descriptions of heating methods hinder development of tailored processing pathways. This study investigates how induction and conventional heating affect microstructures of near-eutectoid Fe–C steel with and without 9 T magnetic field. Conventionally heated samples in field showed more austenite grain growth compared to inductively heated counterparts. Annealing twin frequency increased under field, with twin density increases not fully explained by larger grain sizes.

Hurley, Megan [University of Florida, Gainesville]↗

DOE OSTI2022

What Can Plastic Flow Fields Tell Us About Heat Sources in Deformation Processing?

We characterize primary (shape-change) and secondary (friction) deformation, and associated temperature fields, in metal cutting and forming processes, using in situ imaging and simulation. The experimental configurations enable access to the deformation zones and die contact-interfaces, for measuring deformation, temperature, and frictional drag. Infra-red thermography reveals that the plastic strain-rate field is an excellent proxy for the deformation-induced heat sources. Both spatially confined and diffuse strain-rate fields occur, depending on the initial workpiece deformation state. When the strain rate is confined, as in pre-hardened material, the temperature modeling is much simplified, as the heat source is also now localized. However complex, microstructure-driven deformation modes, like sinuous flow in annealed metals, result in spatially diffuse strain-rate and body heat sources, more challenging to analyze. Furthermore, our unified measurements should be of value for accurately estimating the fraction of plastic dissipation that is converted into heat in large-strain deformation processes.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Experimental and analytical study of the hydrodynamic and single and two-phase convective heat transfer performance of flexible PDMS microchannels with micropillar arrays

Various copper and silicon based thermal management systems are used in the cooling of electronics. However, the rigid nature of these materials along with their high thermal and electrical conductivity pose a difficulty in developing direct contact embedded flexible cooling systems that can offer robust cooling performance. The low density, thermal stability, chemical inertness, and electrical insulation of Polydimethylsiloxane (PDMS) make it an ideal material to develop lightweight direct contact thermal management systems for electronics. Its ease of fabrication with tunable flexibility provides the opportunity to go beyond traditional electronics and develop advanced active and passive thermal management systems for a wide–range of applications in foldable and wearable electronics, liquid cooling garments, microgravity, and electric motors. In this study, a flexible PDMS based microchannel with micropillar arrays, which enhance the thermal performance of the device through capillary-assisted flow, has been developed. The hydrodynamic and convection heat transfer performance of three PDMS wick pillar geometries, ranging from a porosity of 0.8–0.91, are investigated and compared under single-phase and two-phase conditions. Dielectric coolant FC-3283 is employed and permeability measurements are made for mass fluxes ranging from 53 kg/m 2 s to 369 kg/m 2 s. Given its conformability, the device demonstrates a deviation from Darcy’s Law, within the laminar regime, with an increasing permeability with mass flux at the rate of ~0.5–0.8 Darcy/(kg/m 2 s). A semi-analytical model has been developed and reported to quantify the conformability of the device. The heat transfer performance is experimentally evaluated using the same dielectric fluid for mass flux ranging from 105 kg/m 2 s to 420 kg/m 2 s with heat fluxes ranging from 1.5 W/cm 2 to 16 W/cm 2 . Heat transfer coefficients of up to 7000 W/m 2 K are observed, which are comparable to copper and silicon microchannels. The effect of porosity on the single phase thermal performance has been evaluated against the pumping power to provide a basis for thermal management system design. Finally, high-speed imaging is performed to study the two-phase flow characteristics to provide insight into the vapor formation and removal.

42 ENGINEERING↗

DOE OSTI2024

Implementation and validation of optimal start control strategy for air conditioners and heat pumps

Commercial buildings are responsible for approximately 20 % of the total energy consumption and greenhouse gas emissions in the United States. Over 85 % of these buildings lack building automation systems, and many are small (<50,000 square feet), underserved, and use rooftop units (RTUs) for heating, ventilation, and air-conditioning needs. Because these buildings lack proper energy management systems, several operational deficiencies lead to excess energy consumption. Studies have shown that managing the RTUs’ heating and cooling set points, schedules, setbacks, and optimal start can result in a 20 % to 25 % reduction in electricity consumption in small commercial buildings. These buildings typically use fixed schedules to start the RTUs 60 to 120 min before occupancy begins, which results in excess energy consumption. This paper presents research that demonstrates and evaluates the performance of four optimal start methods, which utilize data-based modeling as a key element in facilitating adaptive control in response to time-varying inputs while requiring minimal sensor inputs. The evaluation found energy savings in two commercial buildings equipped with RTUs by periodically alternating four different optimal start models during the cooling and heating season. The resulting energy savings are positive for all models and range from 2 to 5 kWh/day/unit. The units on the east side of the building showed higher savings, while interior units showed greater variability in savings due to the differences in capacities and room sizes. Savings were considerably greater during the heating season compared to the cooling season. The performance of all four models on Mondays was poor; models suggested a shorter optimal start time, which resulted in relatively larger errors. Finally, the future work will look at using a different model for the days after weekends and holidays.

42 ENGINEERING↗

DOE OSTI2026

Virtual refrigerant charge sensor for variable-speed heat pumps based on feature selection

The refrigerant charge level in heat pump systems significantly impacts their energy efficiency. Virtual refrigerant charge (VRC) sensing technology has been comprehensively investigated and well-established due to its lower cost compared to physical sensors. However, the previous VRC research often relied on expert judgment and physical reasoning for their variable selection, which can potentially select redundant (or highly correlated) or insignificant features, and it is also primarily focused on single-speed systems. To address these challenges, this study proposes a VRC algorithm for variable-speed heat pumps that selects features through a rigorous feature selection method in combination with physical insights. We also propose a piecewise linear model structure segmented by subcooling temperature to accurately predict charge levels, particularly when subcooling temperatures are substantially low. The proposed algorithm was evaluated using experimental data of a residential R410A heat pump, and the performance was compared with two baseline VRC algorithms. The results are: (1) The proposed algorithm outperforms for the case with subcooling temperature less than 1 °C. (2) The proposed algorithm achieves a tested mean absolute percentage error (MAPE) of 4.23%, and improves the overall accuracy for cooling conditions by approximately 60%, compared with the two baseline algorithms. (3) The proposed algorithm uses two fewer features and improves the accuracy for undercharge cooling conditions by 68.0%, compared with baseline algorithm 2. These improvements enhance prediction accuracy and prevent overfitting, providing a more reliable refrigerant charge level prediction and helping improve the heat pump energy efficiency.

Liang, Chenjiyu↗

DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Bimodal particle distributions with increased thermal conductivity for solid particles as heat transfer media and storage materials

Solid particles are being considered in several high temperature thermal energy storage systems and as heat transfer media in a variety of advanced power generation systems, particularly in concentrated solar power plants. The downside of such an approach is the low overall heat transfer coefficients caused by the inherently low thermal conductivity values of the low-cost solid media when coupled to heat exchanger for the power cycle working fluid. Choosing the right particle size distribution, emittance, and material of the solid media can all make a substantial difference in packed bed thermal conductivity. Current research though exclusively focuses on continuous unimodal distributions of particles. Here, in this work, we propose the use of a binary particle system with a bimodal size distribution to significantly increase packed bed thermal conductivity by reducing packed bed porosity. This is the first study related to ceramic solid particle heat transfer that has considered the thermal conductivity of non-unimodal size distributions at room and elevated temperatures. The following study found that for the binary particle system using Carbo particles CP 16/30 – CP 70/140 where the large particle volume fraction was 50% there was an 17–47% increase in packed bed thermal conductivity when compared with a nearly unimodal particle size distribution of CP 16/30 between 50 and 300 °C. Two different porosity and effective thermal conductivity models were studied, with one providing better prediction of porosity but both effective thermal conductivity models providing less predictive capacity. Importantly this approach can have a substantial impact of thermal performance, with little to no impact on the particle cost.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nature-inspired lotus-shaped fins combined with hybrid nanoparticles and metal foam for high-performance latent heat thermal energy storage

Latent heat thermal energy storage (LHTES) systems play a critical role in renewable energy integration by providing high energy density and nearly isothermal operation during phase transitions. However, their performance is often limited by slow melting/charging rates, which motivates the search for enhanced heat transfer designs. This study investigates the melting behavior of RT-82 phase change material (PCM) using novel lotus-shaped fins combined with copper metal foam and conductive graphene nanoparticles and carbon nanotubes. A two-dimensional enthalpy-porosity model in ANSYS Fluent was developed to simulate the charging/melting process, capturing non-thermal equilibrium between the foam and PCM/nano-PCM. In this study, effects of fin geometry, nanoparticle concentration, and foam porosity on melting dynamics and cost-performance trade-offs were investigated. Results showed that natural convection accelerated melting by ~12% compared to conduction-only scenarios. Optimized lotus-shaped fins with higher fin density (T3F4 and T3F10) achieved up to 63% faster melting relative to sparse configurations. Graphene nanoparticles improved thermal conductivity, with a 6% volume fraction, by reducing melting time by ~6.9%, while their combination with 75% porosity foam achieved a maximum reduction in the melting time of ~51% compared to pure PCM. Cost-performance analysis identified T3F4 as the most balanced design, offering rapid thermal response without excessive material costs, while moderate-density designs like T3S6 provided economical alternatives with acceptable performance. These results highlight the performance enhancement that can be achieved by integrating bio-inspired fins, nanoparticles, and foams, into compact and efficient LHTES for solar heating, building thermal management, and industrial waste-heat recovery applications.

25 ENERGY STORAGE↗

DOE OSTI2021

Flow condensation heat transfer performance of natural and emerging synthetic refrigerants

There is great interest in predicting flow condensation heat transfer for lower global warming potential (GWP) fluids. This paper analyzes the efficacy of common flow condensation correlations developed for particular fluids in order to identify their suitability to predict heat transfer performance of low GWP fluids. Condensation heat transfer data were extracted from the literature, including 19 papers and 1,473 data points for natural refrigerants [i.e., ammonia (R717), CO 2 (R744), propane (R290), isobutane (R600a)] and 35 papers and 5,030 data points for synthetic refrigerants [i.e., R12, R1234yf, R1234ze(E), R1234ze(Z), R22, R32, R41, R123, R125, R134a, R142b, R152a, R161, R404A, R410A, R448A, R449A, R450A, R452B, R454C, R455A, R513A] encompassing tube diameters of 0.1–11.5 mm, mass fluxes of 55–1200 kg/m2s, and saturation temperatures of -25°C–65°C. Correlations analyzed included Akers et al. (1959), Cavallini et al. (2006, 2011), Kim and Mudawar (2013), Macdonald and Garimella (2016), Shah (1979, 2009, 2013, 2016) and Traviss et al. (1973) for smooth tubes and Chamra et al. (2005) and Kedzierski and Goncalves (1999) for enhanced tubes. Since most studies did not report wall temperature, correlations which relied on wall temperature directly or indirectly were excluded from the analysis. For synthetic refrigerants, mean average error (MAE) ranged from 6%–257%, and Cavallini et al. (2011) and Kim and Mudawar (2013) were the best predictors for emerging synthetic refrigerants. Furthermore, the Kim and Mudawar (2013) correlation was found to best predict the heat transfer performance for propane and R600a data, but most correlations did not accurately predict ammonia and CO 2 flow condensation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Heat capacity and thermodynamic functions of transition metal ion (Cu 2+ , Fe 2+ , Mn 2+ ) exchanged, partially dehydrated zeolite $\mathrm{A}$ ($\mathrm{LTA}$)

Here we have measured the heat capacity from 1.8 to 300 K of partially dehydrated zeolite A (LTA), fully exchanged with Cu 2+ , Fe 2+ , and Mn 2+ ions. The samples have a broad excess heat capacity contribution centered around 4 K, which we attribute to local electric fields splitting the magnetic moments of the cations. The excess heat capacity is modelled using a sum of several Schottky anomalies. From these models, we conclude that the cations in the Cu 2+ zeolite reside in at least four distinct coordination environments, and that some of the coordination environments in all three zeolites are highly asymmetric. We also report theoretical fits of the heat capacity data, and values of the standard thermodynamic functions C P,m , Δ 0K T S m ° , Δ 0K T H m ° , and Φ m ° at smooth temperatures. The standard molar entropies at 298.15 K are 71.8 J·K -1 ·mol -1 for Cu-zeolite A (Cu 0.22 Al 0.49 Si 0.51 O 2 1.04 H 2 O), 71.1 J·K -1 ·mol -1 for Fe-zeolite A (Na 0.01 Fe 0.23 Al 0.50 Si 0.51 O 2 ∙0.77 H 2 O), and 66.0 J·K -1 ·mol -1 for Mn-zeolite A (Mn 0.26 Al 0.49 Si 0.50 O 2 ∙0.53 H 2 O).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Investigating the stability of reoriented hydrides and their reprecipitation using in-situ heating experiments

The behavior of radial hydrides in Zr-2.5Nb pressure tube samples at elevated temperature has been studied. The aim of this work is to investigate the stability of stress-reoriented hydrides under heating without an external applied load. Using in-situ heating at beamline, 1-ID at Argonne National Laboratory, USA, it was shown that two populations of hydrides existed in the sample, naturally oriented and reoriented hydrides. At a relatively low temperature of 180 °C it was observed that the reoriented hydrides completely dissolved, leaving only the naturally oriented hydrides in the sample, when heated without an applied load. On subsequent cooling of the sample, only naturally oriented hydrides were observed. These observations were confirmed using in-situ scanning electron microscopy heating experiments and post-mortem transmission electron microscopy studies. The results presented here offer a potential route for ductility recovery of pressure tube materials where radial hydrides have been shown to significantly reduce fracture toughness of components.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Microstructural evolution, defect mitigation, and precipitation behavior in AA6061 via laser powder bed fusion with high-temperature substrate heating

Defects, particularly solidification cracking, remain persistent challenges in the laser powder bed fusion (LPBF) processing of AA6061 aluminum alloy. This study systematically investigates defect mitigation, microstructural evolution, and mechanical properties associated with high-temperature substrate preheating at 500 °C. Comprehensive microstructural analyses, including characterization of defects, grain structures, and precipitation behavior, were performed on samples in both as-built and T6 heat-treated states. Elevated preheating substantially reduced solidification cracking across a wide processing window, while demonstrating decreased crack sensitivity to laser parameters. Columnar cracks along the build direction were observed despite substrate preheating. Lack-of-fusion and keyhole porosity were effectively eliminated, though gas-induced microporosity persisted at higher powers. In-depth characterization of two distinct laser power and speed conditions confirmed the formation of micron-sized, non-coherent Mg 2 Si precipitates under heated substrate conditions, alongside α-AlFeCrMnSi intermetallic phases indicative of in-situ thermal effects during fabrication. Subsequent T6 heat treatment revealed the formation of fine, coherent needle-shaped β″ strengthening precipitates. Despite substantial differences in processing parameters, comparable mechanical properties were measured in the as-built samples (∼52 MPa yield strength, ∼130 MPa tensile strength), primarily due to reduced strain hardening effects and consistent precipitation characteristics. Meanwhile, the T6 heat treatment led to significant improvement in properties, enhancing yield strength by over 400%, aligning closely with the performance of conventional wrought AA6061-T6. These findings underscore that high-temperature substrate preheating offers an effective means to suppress cracks, control precipitation, and enhance mechanical performance in LPBF-processed AA6061.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Dynamic condensation model of rolling droplets for high-performance heat transfer

Condensation is a critical process during vapor-liquid phase change in relation to heat transfer. To achieve a high heat transfer coefficient, the classical model for dropwise condensation requires a low contact angle and low contact angle hysteresis, failing to align with experimental observations on a hydrophobic and slippery quasi-liquid surface (QLS). Here, we report a dynamic condensation model that incorporates high-frequency condensate removal by emphasizing the role of timescale during droplet growth and shedding. Our model agrees well with the experimental result that a surface with high contact angle and low contact angle hysteresis promotes condensation, particularly during rolling-propelled condensate removal. Particle image velocimetry reveals that rolling droplets on a hydrophobic QLS exhibit 4-fold higher shedding speeds than the sliding droplets on a hydrophilic QLS, leading to significant heat transfer enhancement. This work deepens our theoretical understanding of condensation heat transfer and provides advanced physics-informed design rationales for water and energy systems.

42 ENGINEERING↗

DOE OSTI2026

The effectiveness of D 2 pellet and gas injection in reducing intra-ELM tungsten erosion and heat flux in the DIII-D small angle slot divertor

Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.

Carbon (C)↗

DOE OSTI2024

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗

DOE OSTI2024

Enhancement of disposal efficiency for deep geological repositories based on three design factors - Decay heat optimization, increased thermal limit of the buffer and double-layer concept

This study investigates the enhancement of disposal efficiency for deep geological repositories (DGRs) based on three design factors: decay heat optimization, increased thermal limit of the buffer, and double-layer concept using coupled thermo-hydro-mechanical (THM) numerical simulations. Decay heat optimization is achieved by iteratively emplacing spent nuclear fuels having the maximum and minimum decay heat in a canister. Disposal areas can be reduced by 20 % to 40 % compared to the current reference disposal system in Korea (KRS+) in accordance with the combinations of the three design factors, alleviating challenges in site selection for the DGR. This study additionally identifies an optimal layer spacing of 500 m for the double-layer concept in the viewpoint of the buffer temperature, where thermal interaction between the upper and lower layers nearly disappears. However, determining the ultimate disposal and layer spacing requires engineering judgement, considering not only the thermal performance of the DGR but also various factors such as cost and difficulties of the construction and rock mass stability. DGRs designed with an increased thermal limit of the buffer poses a greater probability of rock mass failure around disposal tunnels and deposition holes due to elevated thermal stresses. Densely arranged heat sources for the DGRs with enhanced disposal efficiency lead to larger temperature increase even at the far-field scale, raising a possibility of thermally driven fracture shear activation with associated hydraulic, mechanical, and seismic changes.

58 GEOSCIENCES↗

DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗