SEARCH · Search NASA

Results for “module cooling”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

The Subchalcogenides Ir 2 In 8 Q (Q = S, Se, Te): Dirac Semimetal Candidates with Re-entrant Structural Modulation

Subchalcogenides are uncommon compounds where the metal atoms are in unusually low formal oxidation states. They bridge the gap between intermetallics and semiconductors and can have unexpected structures and properties because of the exotic nature of their chemical bonding as they contain both metal–metal and metal–main group (e.g., halide, chalcogenide) interactions. Finding new members of this class of materials presents synthetic challenges as attempts to make them often result in phase separation into binary compounds. We overcome this difficulty by utilizing indium as a metal flux to synthesize large (millimeter scale) single crystals of novel subchalcogenide materials. Herein, we report two new compounds Ir 2 In 8 Q (Q = Se, Te) and compare their structural and electrical properties to the previously reported Ir 2 In 8 S analogue. Ir 2 In 8 Se and Ir 2 In 8 Te crystallize in the $P_{4_2}$/mnm space group and are isostructural to Ir2In8S, but also have commensurately modulated (with q vectors q = 1/6a* + 1/6b* and q = 1/10a* + 1/10b* for Ir 2 In 8 Se and Ir 2 In 8 Te, respectively) low-temperature phase transitions, where the chalcogenide anions in the channels experience a distortion in the form of In–Q bond alternation along the ab plane. Both compounds display re-entrant structural behavior, where the supercells appear on cooling but revert to the original subcell below 100 K, suggesting competing structural and electronic interactions dictate the overall structure. Notably, these materials are topological semimetal candidates with symmetry-protected Dirac crossings near the Fermi level and exhibit high electron mobilities (~1500 cm 2 V –1 s –1 at 1.8 K) and moderate carrier concentrations (~10 20 cm –3 ) from charge transport measurements. This work highlights metal flux as a synthetic route to high quality single crystals of novel intermetallic subchalcogenides with Dirac semimetal behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

An oceanic pathway for Madden–Julian Oscillation influence on Maritime Continent Tropical Cyclones

While the Madden–Julian Oscillation (MJO) has been shown to affect tropical cyclones (TCs) worldwide through its modulation of large-scale circulation in the atmosphere, little or no role for the ocean has been identified to date in this influence of MJO on TCs. Using observations and numerical model simulations, we demonstrate that MJO events substantially impact TCs over the Maritime Continent (MC) region through an oceanic pathway. While propagating across the MC region, MJO events cause significant sea surface cooling with an area-averaged value of about 0.35 ± 0.12 °C. Hence, TCs over the MC region immediately following the passage of MJO events encounter considerably cooler sea surface temperatures. Consequently, the enthalpy fluxes under the storms are reduced and the intensification rates decrease by more than 50% on average. These results highlight an important role played by the ocean in facilitating MJO-induced sub-seasonal variability in TC activity over the MC region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗

DOE OSTI2025

Thermal Modeling and Limitations for Power Electronics Embedded in Medium-Voltage Cables

As next-generation energy technologies gain traction and power demand increases, the existing electrical infrastructure faces significant stress, prompting innovative solutions to enhance the grid's capacity and lifespan. This work explores the possibility of embedding medium-voltage (MV) power electronics directly inline with the cable, and the resulting thermal challenges. Since the majority of power distribution cables installed in the U.S. are passively cooled, the work focuses primarily on passive cooling, with an emphasis on the limitations of axial heat spreading within the cable. To date, literature on axial spreading of high incident heat loads on cables and cable environments is limited, typically reporting cases with <10 W of incident heat load. This work will explore the considerations, limits, and tradeoffs of cable-embedded heat loads significantly larger than the cable losses. Both external and internal effects are modeled analytically in nondimensional terms via a Biot number analysis, allowing fundamental limits and tradeoffs to be derived. The work culminates in the design and experimental validation of a cable-embedded thermal system capable of passively dissipating 300 W of heat from a coaxial SiC mosfet switch module over a length of 20 cm, thus validating the possibility of MV cable-embedded power electronics from a thermal standpoint.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Q-dependent collective relaxation dynamics of glass-forming liquid Ca0.4K0.6(NO3)1.4 investigated by wide-angle neutron spin-echo

Abstract The relaxation behavior of glass formers exhibits spatial heterogeneity and dramatically changes upon cooling towards the glass transition. However, the underlying mechanisms of the dynamics at different microscopic length scales are not fully understood. Employing the recently developed wide-angle neutron spin-echo spectroscopy technique, we measured the Q -dependent coherent intermediate scattering function of a prototypical ionic glass former Ca 0.4 K 0.6 (NO 3 ) 1.4 , in the highly viscous liquid state. In contrast to the structure modulated dynamics for Q < 2.4 Å −1 , i.e., at and below the structure factor main peak, for Q > 2.4 Å −1 , beyond the first minimum above the structure factor main peak, the stretching exponent exhibits no temperature dependence and concomitantly the relaxation time shows smaller deviations from Arrhenius behavior. This finding indicates a change in the dominant relaxation mechanisms around a characteristic length of 2π/(2.4 Å −1 ) ≈ 2.6 Å, below which the relaxation process exhibits a temperature independent distribution and more Arrhenius-like behavior.

42 ENGINEERING↗

DOE OSTI2023

Phonon softening and slowing-down of charge density wave fluctuations in BaNi 2 As 2

BaNi 2 As 2 is a nonmagnetic analog of the iron pnictide superconductors, and exhibits an incommensurate charge density wave (IC-CDW) and a sizable elastoresistance. Here, in this Letter, phonons in BaNi 2 As 2 associated with the IC-CDW and uniform in-plane lattice distortions are investigated using high-resolution inelastic x-ray scattering. The in-plane transverse acoustic phonons reveal no softening at temperatures where the elastoresistance increases strongly, indicating the latter to be electronically driven. Systematic phonon measurements suggest the IC-CDW occurs in two stages upon cooling: Underdamped phonons first soften to zero energy well above the IC-CDW ordering temperature, then the resulting quasielastic IC-CDW fluctuations gradually slow down and coalesce into the static IC-CDW order. A possible origin for our observations is the IC-CDW in BaNi 2 As 2 being uniaxial, which provides an additional Ising degree of freedom favorable for disordered IC-CDW modulations, and accounts for the elastoresistance through a weak coupling to the lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Longitudinal space charge kick in Coherent Electron Cooling

In recent EIC Coherent electron Cooling (CeC) design, we adopted 591 MHz SRF LINAC to accelerate electron beam to certain energy(150 MeV, 54.1 MeV and 22.3 MeV). To reduce the energy spread caused by the RF curvature, a short bunch is preferred and results in high peak current. The longitudinal space charge kick will cause the energy spread increasing and lengthen the bunch length when through the R56 elements. This slippage of modulated microbunches may misalign with hadrons providing imprinting at kicker section. In this note, we will estimate this effect and discuss possible solution.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Mobility of a spatially modulated electron liquid on the helium surface

Here, we present the results of large-scale numerical simulations of the mobility of a two-dimensional electron liquid on the helium surface in the presence of a one-dimensional periodic potential. Even where the potential is much weaker than the electron-electron interaction, it can strongly change the mobility. The effect depends on the interrelation between the potential period and the mean interelectron distance. It is most pronounced where the period is close to the period of the Wigner crystal that would form if the liquid were cooled to a lower temperature. The results suggest, in particular, that the correlation length in the electron liquid can be found by measuring the mobility in a weak periodic potential. The simulations are based on the microscopic model of the electron scattering by the excitations in helium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Next Gen High Efficiency Boosted Engine Development

This work represents an advanced engineering research project partially funded by the U.S. Department of Energy (DOE). Ford Motor Company, FEV North America, and Oak Ridge National Laboratory collaborated to develop a next generation boosted spark ignited engine concept. The project goals, specified by the DOE, were 23% improved fuel economy and 15% reduced weight relative to a 2015 or newer light-duty vehicle. The fuel economy goal was achieved by designing an engine incorporating high geometric compression ratio, high dilution tolerance, low pumping work, and low friction. The increased tendency for knock with high compression ratio was addressed using early intake valve closing (EIVC), cooled exhaust gas recirculation (EGR), an active pre-chamber ignition system, and careful management of the fresh charge temperature. Engine weight reduction measures were implemented throughout the engine system making use of composite materials, advanced manufacturing techniques, and architectural choices. This report outlines the analytical, design, fabrication, and test work conducted for the duration of the project. The combustion system stability, EGR tolerance, and knock resistance were validated on a single cylinder engine. An inline six-cylinder engine was then designed targeting application in the Ford F150. Multi-cylinder engines were produced and tested achieving the target vehicle fuel economy improvement of 23% assessed using measured engine fuel consumption combined with a vehicle drive cycle simulation. Actions were identified and designs were demonstrated to achieve the 15% weight reduction target. This project included items covering a range of technology readiness levels. Some of the technologies explored are production ready, while others were investigated to understand the limitations for what can be achieved in a stoichiometric, gasoline-fueled, spark-ignited internal combustion engine.

42 ENGINEERING↗

DOE OSTI2021

High efficiency compact microwave sources using electrostatic potential depression

In this article, we report theoretical and design studies on the application of an electrostatic potential depression (EPD) on two types of traditional microwave sources, the klystron and the inductive output tube (IOT), both operating in the UHF band. Additionally, the EPD is implemented by a section of the metallic beam pipe on which a direct current (dc) negative high voltage is applied. Inside the EPD section, the modulated electron beam develops into bunches over a very short longitudinal distance, indicating a possible design approach for a compact device. With the application of an EPD, the peak value of the first harmonic current in a microwave tube can be further increased compared to that in a conventional design, suggesting higher circuit efficiency. As the electron beam leaves the EPD section, the reacceleration of the beam provides a longitudinal cooling effect. With the addition of an EPD, circuit power extraction efficiency above 80% is achieved in both the modified design of the IOT and the novel conceptual circuit design of the klystron.

42 ENGINEERING↗

DOE OSTI2023

Magnetic Energy Dissipation and γ-Ray Emission in Energetic Pulsars

Abstract Some of the most energetic pulsars exhibit rotation-modulated γ -ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power ( γ -ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower γ -ray efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗

DOE OSTI2025

Atmospheric River Frequency-Category Characteristics Shape U.S. West Coast Runoff

Abstractrunoff response to atmospheric rivers (ARs) over the U.S. West Coast. We focused on runoff time series variations impacted by AR characteristics (e.g., category and frequency) and land preconditions during Northern Hemisphere cool seasons in the period of 1940–2023. Results show that high-category ARs significantly increase local runoff with higher hourly precipitation rates leading to a greater incremental rate and peak runoff. Extreme runoff increases greatly with the AR category with an increase rate up to 12.5 times stronger than non-extreme runoff. Besides the AR category, land preconditions such as soil moisture and snowpack also play crucial roles in modulating runoff response. We found that runoff induced by weak-category ARs is more sensitive to land preconditions than high-category ARs, with high peak runoff occurring when soil is nearly saturated. Additionally, more than 50% of high-peak-runoff events in snow-covered grid cells are associated with rain-on-snow events particularly for the events associated with weaker ARs. Regression analysis reveals that AR precipitation and land preconditions jointly influence runoff, emphasizing the importance of including soil moisture and snowpack levels in AR impact assessments. The study also highlights the intensified runoff response to back-to-back ARs with short intervals, which may become more frequent with climate warming, posing increased flood risks via facilitating wet soil conditions. Our findings have significant implications for AR risk predictions and the development of prediction models for AR-induced runoff.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗