Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal effects”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Site specific porosity-thermal performance correlations in neutron irradiated U-10Zr fuel

In this work, we present a site-specific three-dimensional analysis of porosity evolution in neutron-irradiated U–10Zr metallic fuel using high-resolution synchrotron X-ray tomography. Focused ion beam cubic lift-outs from four radial positions—fuel center, middle, edge, and fuel-cladding interaction (FCCI) zones—were imaged, reconstructed, and segmented to quantify pore volume, density, morphology, and pore connectivity. Porosity increased modestly from 5.5% at the center to 10.5% at the edge, yet pore number density increased by over two orders of magnitude in the fuel portion near the FCCI interface (from 4.8 x 10 4 to 6.0 x 10 6 pores/mm 3 ). Morphological classification revealed a progression from small spherical pores at the center to equiaxed and tortuous networks at the periphery, with enhanced orientation along the radial direction. Connectivity and permeability analysis reveal that the FCCI region maintains dense, highly interconnected pores despite a moderate volume fraction, enabling rapid fission gas transport and lanthanide migration. Effective thermal conductivity models incorporating sodium logging confirm that these site-specific pore features critically influenced heat transport during reactor irradiation. These findings demonstrate that pore topology—not just porosity fraction—controls thermal performance and thermal pathways in U–10Zr fuels, impacting fuel thermal performance in a reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Quantum Monte Carlo determination of the principal Hugoniot of deuterium

Here, we present coupled electron-ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock-wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the potential energy surface by wave-function optimization within variational Monte Carlo and projection quantum Monte Carlo methods. Compared to a previous study, our calculations also include low pressure-temperature (P,T) conditions resulting in close agreement with experimental data, while our revised results at higher (P,T) conditions still predict a more compressible Hugoniot than experimentally observed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Uncertainty in Thermal Modeling of Spent Nuclear Fuel Casks

Uncertainty is a key metric in computational modeling that must be evaluated for results to have wide ranging applicability. A well characterized uncertainty range is ideal with clear error bars on results that can be presented to stakeholders. In the field of spent fuel cask modeling, this ideal has been historically difficult to achieve in practice because of the computationally intensive nature of the models used and the difficulty assigning reasonable uncertainties to quantities in as-built systems. The work in this report has been conducted to evaluate the overall state of uncertainty and sensitivity in spent fuel cask models and develop methodologies for evaluating these uncertainties. These methodologies must be practical for engineering applications. They should not require excessive computational resources or calendar time to achieve results. In engineering, the model must be on a scale such that it can be changed and adapted throughout a project as new information is discovered and project goals evolve. This report covers three major modeling task areas that provide an overview of the types of sensitivity and uncertainty present in a spent fuel storage and transportation system. Section 3 discusses sensitivity and uncertainty analysis in the effective thermal conductivity model for the fuel region and applies these results to a single assembly model. Section 4 shows sensitivity analysis of a full cask model in the TN-32B and Section 5 demonstrates the overall uncertainty workflow using Coolant Boiling in Rod Arrays – Spent Fuel Storage and STAR-CCM+ developed from the sensitivity work in the preceding sections.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Numerical investigation of the effect of ultrasound on paper drying

The paper drying process is very energy inefficient. More than two-thirds of the total energy used in a paper machine is for drying paper. Novel drying technologies, such as ultrasound (US) drying, can be assessed numerically for developing next-generation drying technologies for the paper industry.This work numerically illustrates the impact on drying process energy efficiency of US transducers installed on a two-tiered dryer section of a paper machine. Piezoelectric transducers generate ultrasound waves, and liquid water mist can be ejected from the porous media. The drying rate of handsheet paper in the presence of direct-contact USis measured experimentally, and the resultant correlation is included in the theoretical model.The drying section of a paper machine is simulated by a theoretical drying model. In the model, three scenarios are considered. In the first scenario, the US modules are positioned in the dryer pockets, while in the second scenario, they are placed upstream of the drying section right after the press section. The third case is the combination ofthe first and second scenarios. The average moisture content and temperature during drying, enhancement of total mass flux leaving the paper by the US mechanism, total energy consumption, and thermal effect of heated US transducers are analyzed for all cases. Results show that the application of the US can decrease the total number of dryer drums for drying paper. This numerical study is based on the US correlation obtained with the US transducer in direct contact with the paper sample. Thus, future work should include US correlation based on a non-contact US transducer.

Materials Science↗
DOE OSTI2023

A high-performance lung-inspired ceramic 3D-printed heat exchanger for high-temperature energy-efficient systems

The thermal-to-electric conversion efficiency of next-generation supercritical carbon dioxide and concentrated solar power plants is significantly augmented by operating at temperatures exceeding 1000 °C if durable, compact, high-temperature heat exchangers (HXs) are developed. Whereas metal 3D-printed HXs fail to operate under extreme temperature conditions, ceramic 3D-printed HXs are deemed promising candidates for high-temperature, highly oxidizing environments. Ceramic 3D-printed HXs, however, exhibit a high gas permeability issue through hot and cold separating walls due to an excessive non-volatile photopolymer content utilized. Here, in this work, a novel leak-free lung-inspired ceramic 3D-printed HX employing a highly complex heat transfer topology is introduced for high-temperature energy-efficient systems. The high permeability issue of ceramic 3D-printed HXs is fully eliminated through a uniform zinc-based coating. It was determined that the effective thermal conductivity of the sintered ceramic heat exchangers is affected by the coated layer. The high-temperature thermal and hydraulic characteristics of the lung-inspired ceramic 3D-printed HX are comprehensively examined and compared against a millichannel 3D-printed HX module. Experimental results indicated that the lung-inspired 3D-printed ceramic HX significantly outperforms its millichannel counterpart. The lung-inspired HX shows a volume-based power density of 8.2 MW/m3 at a hot-side inlet air temperature of 700°C, a 71 % improvement compared with that of the millichannel HX. More importantly, the higher thermal duty of the lung-inspired 3D-printed HX is realized at a lower normalized pressure drop penalty of 510 Pa/W, 22 % lower than that of the millichannel ceramic HX. This study accelerates the advancement of high-performance ceramic 3D-printed HXs with complex topologies for high-temperature energy-efficient systems and/or corrosive environments.

42 ENGINEERING↗
DOE OSTI2025

Nanoscale and Element-Specific Lattice Temperature Measurements Using Core-Loss Electron Energy-Loss Spectroscopy

Measuring nanoscale local temperatures, particularly in vertically integrated and multicomponent systems, remains challenging. Spectroscopic techniques like X-ray absorption and core-loss electron energy-loss spectroscopy (EELS) are sensitive to lattice temperature, but understanding thermal effects is nontrivial. This work explores the potential for nanoscale and element-specific core-loss thermometry by comparing the Si L2,3 edge’s temperature- dependent redshift against plasmon energy expansion thermometry (PEET) in a scanning TEM. Using density functional theory (DFT), time-dependent DFT, and the Bethe−Salpeter equation, we ab initio model both the Si L 2,3 and plasmon redshift. We find that the core-loss redshift occurs due to bandgap reduction from electron−phonon renormalization. Our results indicate that despite lower core-loss signal intensity compared to plasmon features, core-loss thermometry has key advantages and can be more accurate through standard spectral denoising. Specifically, we show that the Varshni equation easily interprets the core-loss redshift for semiconductors, which avoids plasmon spectral convolution for PEET in complex junctions and interfaces. We also find that core-loss thermometry is more accurate than PEET at modeling thermal lattice expansion in semiconductors, unless the specimen’s temperature-dependent dielectric properties are fully characterized. Furthermore, core-loss thermometry has the potential to measure nanoscale heating in multicomponent materials and stacked interfaces with elemental specificity at length scales smaller than the plasmon’s wave function.

Bethe–Salpeter equation↗
DOE OSTI2020

Foam Degradation Model during Thermal Hypothetical Accident Conditions

Current practice in Safety Analysis Report for Packaging (SARP) thermal analysis for hypothetical accident conditions (HAC) of packages is to not explicitly model the consequence to the foam impact limiting material, but instead to make conservative estimates of those consequences and begin calculations from that point. This approach is typically used for the foam regression distance, which is the thickness of the degraded foam layer after the 30-minute fire event. Estimates for regression distance can come directly from package burn test results or from correlations based on one-dimensional testing by the foam manufacturer. In either case, this is a subjective process, and although it can be used to produce a conservative result, it is a significant approximation in the analysis. This paper describes testing and refinement of a simplified foam regression model. Comparisons are made against available data from laboratory experiments. Model refinements were tested until satisfactory agreement was achieved. The final modeling approach was then compared with post-test inspection results from HAC burn test data. A new model with an effective thermal conductivity based on foam to char density variation was shown successful in representing laboratory-scale foam thermal degradation experiments and full-scale burn tests of a current package. This model offers a technically defensible and predictive approach for modeling thermal degradation of rigid polyurethane insulation in transport packages.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of N 2 /O 2 plasma interaction with Pt-catalyst: effect of metastable adsorbates on product hysteresis

The coupling of catalysts and atmospheric-pressure plasma has the potential to improve the efficiency of certain catalytic reactions. Understanding the changes that the catalyst surface undergoes during exposure to plasma is key to improving plasma–catalytic performance. In this work, long term exposure of Pt–Al 2 O 3 powder catalyst to an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet was investigated. Products produced by the interaction were analyzed downstream with Fourier-transform infrared spectroscopy while surface species were analyzed operandi with diffuse reflectance infrared Fourier transform spectroscopy. During exposure, the catalyst temperature was ramped cyclically between 100 °C and 350 °C to understand how substrate temperature affects the plasma–catalyst interaction. Long-lasting changes were revealed to take place on the catalyst surface during plasma exposure. At low temperatures, Pt–O and Pt–NO accumulate on the surface which react at elevated temperatures to form NO 2 . NO 2 initially appears to spill on to the Al 2 O 3 support as nitrites and nitrates instead of desorbing. Stable surface conditions are only achieved after prolonged plasma exposure, when nitrate sites on the Al 2 O 3 support are filled. By changing the catalyst temperature at various rates, the impact of total plasma species flux to the surface was analyzed. It was found that decreasing the heating rate increased the hysteresis in the pattern of NO 2 formation during thermal cycling. The variation with temperature demonstrates that plasma exposure results in a buildup of surface NO x and oxygen species which react or desorb at high temperatures. The observed changes are discussed from the generic viewpoint that a non-equilibrium plasma interacting with a catalyst at low temperature introduces metastable steady-state surface conditions. Upon heating above a threshold temperature, the introduced surface modifications can change either due to thermal effects, or, for a plasma environment, by additional interaction with the incident plasma species flux. The surface/material changes take place in a highly predictable fashion and after sufficient time above the threshold temperature reach a steady-state condition that is different from the transient behavior that is observed during initial heating. During cooling the plasma-surface interaction exhibits a different behavior than during heating, and this results in hysteresis of diverse observables. The metastability/hysteresis description appears quite generic and analogous to hysteresis behavior seen for different systems. Furthermore, it is expected to be useful for understanding the consequences of plasma–catalyst surface interactions for various systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Effect of Heat Treatment on Microstructure and Mechanical Property of 316L Stainless Steel Produced by Laser Powder Bed Fusion

The advanced non-light water reactor designs (Gen IV reactors), including molten salt/ very high temperature/ sodium-cooled and lead-cooled fast reactors, typically operate at higher temperatures and more extreme radiation conditions than light water reactors. An intrinsic part of the deployment and progress of Gen IV reactor designs is selecting the most suitable structural material for a specific application. Additive manufacturing (AM), a fairly new process of making physical, three-dimensional objects from a computer design file, is going to completely change the way of design, build and certify nuclear systems. It offers a range of opportunities to produce complex geometries from existing materials, offers new routes for processing of previously difficult to process materials, allows for design of new high-performance materials, and finally facilitates hybridization of dissimilar materials. This emerging technology has successfully produced cars, wind turbine blade molds and even live cells. It could also open up big opportunities for the nuclear industry to quickly deploy technologies at a fraction of the cost. So far, AM techniques have been preliminarily applied in the field of nuclear reactors, including the classical parts such as the pressure vessel of a small reactor with 508-III steel, the bottom nozzle of a fuel assembly with 304L steel, the fuel cladding with zirconium alloy and the integrated impeller of a pump and the multi-channel valve body with 316L steel [6,7]. The AM applications for operating nuclear reactors started in auxiliary plant components and have slowly migrated to metallic reactors and core components, but many of these are not safety critical components. Although many parts used for nuclear reactors have been fabricated by AM techniques, practical applications in engineering are still a long way off due to the uncertainty factors focused on the processing, material properties, analysis methods and application standards, which feeds the safety and life-cycle of the nuclear reactor. Due to rapid, repeated heating and cooling during production, a high dislocation density was present in the AM material. This microstructure feature is unstable at elevated temperature while high temperature is one of the typical operation environments for nuclear reactors. Thus, it is important to understand the thermal effect on the microstructure of AM material. The objectives of this study are to investigate the effect of heat treatment on the microstructure and mechanical properties of 316L stainless steel produced by laser powder bed fusion additive manufacturing, and to determine an appropriate heat treatment practice that will be applied to the lightweight AM lattice-structured material with the same chemistry. The heat treatment study consisted of annealing the samples at a temperature range of 800 to 1200 oC with a 50 oC increment for different times (1-24 hours), followed by vacuum or air cooling. Microstructural characterization was carried out by Scanning Electron Microscope (SEM). Grain size and crystallographic orientation were investigated by Electron Backscatter Diffraction (EBSD). Vickers hardness tests with a 0.5 kg load were employed to determine the hardness of samples after different heat treatments. After heat treatment, the random crystallographic orientation was preserved, and the volume fraction of high-angle grain boundaries (grain boundary misorientation =15 oC) remained the same. The dislocation density decreased with annealing temperature due to recovery. The fine subgrain structures in the as-printed specimen were quite stable up to 1200 oC. Minimal recrystallization was observed up to 1200 oC. Recrystallization initiated only after 8.5 hours at 1200 oC. The SEM images did not show obvious dependence of microstructure on cooling rate. The hardness of the specimens decreased with increasing annealing temperature as a result of the decrease in dislocation density. It is interesting to note that the AM material showed very similar hardness to the wrought material when annealing at similar temperature, although the microstructures are very different. Annealing at 1050 oC for 1 hour followed by air cooling was selected as the heat treatment procedure for the lattice designed lightweight AM 316L material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct Observation of Room–Temperature Magnetic Skyrmion Motion Driven by Ultra–Low Current Density in Van Der Waals Ferromagnets

The recent discovery of room-temperature ferromagnetism in 2D van der Waals (vdW) materials, such as Fe 3 GaTe 2 (FGaT), has garnered significant interest in offering a robust platform for 2D spintronic applications. Various fundamental operations essential for the realization of 2D spintronics devices are experimentally confirmed using these materials at room temperature, such as current-induced magnetization switching or tunneling magnetoresistance. Nevertheless, the potential applications of magnetic skyrmions in FGaT systems at room temperature remain unexplored. In this work, the current-induced generation of magnetic skyrmions in FGaT flakes employing high-resolution magnetic transmission soft X-ray microscopy is introduced, supported by a feasible mechanism based on thermal effects. Furthermore, direct observation of the current-induced magnetic skyrmion motion at room temperature in FGaT flakes is presented with ultra-low threshold current density. This work highlights the potential of FGaT as a foundation for room-temperature-operating 2D skyrmion device applications.

2D ferromagnets↗
DOE OSTI2021

Dark matter detection, Standard Model parameters and Intermediate Scale Supersymmetry

The vanishing of the Higgs quartic coupling at a high energy scale may be explained by Intermediate Scale Supersymmetry, where supersymmetry breaks at (10 9 -10 12 ) GeV. The possible range of supersymmetry breaking scales can be narrowed down by precise measurements of the top quark mass and the strong coupling constant. On the other hand, nuclear recoil experiments can probe Higgsino or sneutrino dark matter up to a mass of 1012 GeV. We derive the correlation between the dark matter mass and precision measurements of standard model parameters, including supersymmetric threshold corrections. The dark matter mass is bounded from above as a function of the top quark mass and the strong coupling constant. The top quark mass and the strong coupling constant are bounded from above and below respectively for a given dark matter mass. We also discuss how the observed dark matter abundance can be explained by freeze-out or freeze-in during a matter-dominated era after inflation, with the inflaton condensate being dissipated by thermal effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Lepto-axiogenesis in minimal $\mathrm{SUSY}$ $\mathrm{KSVZ}$ model

We study the lepto-axiogenesis scenario in the minimal supersymmetric KSVZ axion model. Only one Peccei-Quinn (PQ) field and vector-like fields are introduced besides the MSSM with the type-I see-saw mechanism. The PQ field is stabilized by the radiative correction induced by the Yukawa couplings with the vector-like fields introduced in the KSVZ model. We develop a way to follow the dynamics of the PQ field, in particular we found a semi-analytical solution which describes the rotational motion under the log-arithmic potential with including the thermalization effect via the gluon scattering which preserves the PQ symmetry. Based on the solution, we studied the baryon asymmetry, the effective number of neutrino, and the dark matter density composed of the axion and the neutralino. We found that the baryon asymmetry is successfully explained when the mass of PQ field is $\mathcal{O}$(10 6 GeV) ($\mathcal{O}$(10 5 GeV)) with the power of the PQ breaking term being 10 (8).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Correlated signals of first-order phase transitions and primordial black hole evaporation

Fermi balls produced in a cosmological first-order phase transition may collapse to primordial black holes (PBHs) if the fermion dark matter particles that comprise them interact via a sufficiently strong Yukawa force. We show that phase transitions described by a quartic thermal effective potential with vacuum energy, 0.1 ≲ B 1/4 /MeV ≲ 10 3 , generate PBHs of mass, 10 –20 ≲ M PBH /M ⊙ ≲ 10 –16 , and gravitational waves from the phase transition (at THEIA/μAres) can be correlated with an isotropic extragalactic X-ray/γ-ray background from PBH evaporation (at AMEGO-X/e-ASTROGAM).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Correlated gravitational wave and microlensing signals of macroscopic dark matter

Fermion dark matter particles can aggregate to form extended dark matter structures via a first-order phase transition in which the particles get trapped in the false vacuum. We study Fermi balls created in a phase transition induced by a generic quartic thermal effective potential. We show that for Fermi balls of mass, 3 x 10 -12 M ⊙ ≲ M FB ≲ 10 -5 M ⊙ , correlated observations of gravitational waves produced during the phase transition (at SKA/THEIA/μAres), and gravitational microlensing caused by Fermi balls (at Subaru-HSC), can be made.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗
DOE OSTI2023

Thermo-Hydro-Mechanical Modeling of Brine Migration in a Heated Borehole Test in Bedded Salt

Abstract This research paper focuses on the thermo-hydro-mechanical (THM) modeling of brine migration in a heated borehole test conducted as part of the ongoing Brine Availability Test in Salt (BATS) at the Waste Isolation Pilot Plant (WIPP) in New Mexico. It is a component of the international collaboration project DECOVALEX-2023 (DEvelopment of COupled models and their VALidation against EXperiments), which aims to understand the THM processes governing brine flow in heated rock salt repositories through collaborative analysis by multiple research teams. Using the TOUGH–FLAC simulator, THM simulations were performed and compared with data from the BATS phase 1a. This experiment involved two identical horizontal-borehole arrays, one heated and one serving as a control, both equipped with sensor arrays. Analysis of measurements revealed water flow rate surges during heater power transitions, with the highest jump observed during cooling. Acoustic emission activity exhibited distinct patterns in response to heater power changes, suggesting that damage to rock salt is particularly pronounced during the cooling phase. The THM simulations successfully captured these phenomena, highlighting the significance of thermal effects, brine migration, and mechanical behavior in predicting brine availability in heated and damaged rock salt. Our modeling also revealed the critical interplay between heating and cooling-induced damage and its influence on flow properties, particularly affecting brine inflow estimation. Notably, we found that cooling-induced brine inflow spikes result from increased permeability due to tensile dilatancy. These findings have important implications for the development of robust containment strategies and enhance our understanding of the complex processes involved in repository performance.

58 GEOSCIENCES↗