SEARCH · Engineering Papers

Results for “process 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 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A hybrid method to evaluate the life cycle climate performance of heat pumps

Heat pumps have a significant impact on the climate, and numerous studies exist discussing the environmental impact of heat pumps quantitatively. Life Cycle Climate Performance (LCCP) is a widely accepted metric and can evaluate the carbon footprint of heat pumps from cradle to grave. LCCP calculation typically requires inputs like annual energy consumption, the material used, refrigerant type, and charge level, which are usually the intermediate or final results from heat pump simulation tools. Thus, a platform linking the heat pump simulation software, local weather conditions, and LCCP evaluation tool can simplify the heat pump design process, especially for complicated systems. This paper presents an LCCP evaluation tool accessing the outputs from the DOE Heat Pump Design Model and the building loads and temperature bin data from AHRI standard 210/240. As for case studies, this paper investigates the LCCP values and life cycle electricity costs (LCEC) of a given cooling-based heat pump with R410A in five typical US cities. In conclusion, this cooling-based heat pump has the highest LCCP value at Chicago and lowest LCCP value at Seattle. In contrast, the heat pump has the highest LCEC at Seattle and the lowest LCEC at Miami. As a future work, a more comprehensive carbon emission and cost analysis of heat pump systems will be conducted through all major cities around the world.

Shen, Bo↗

DOE OSTI2023

Design and Prototyping of a Novel Toroidal Magnet System for MOLLER Experiment at Jefferson Lab

The Thomas Jefferson National Accelerator Facility (JLab) has designed a unique spectrometer system to measure the weak interaction between electrons. The experiment— Measurement of Lepton-Lepton Electroweak Reaction (MOLLER)—requires leveraging the recent 12 GeV electron beam upgrade and will run in JLab for three years. Focusing the signal for the MOLLER experiment requires five water-cooled toroidal magnets, each with unique geometry and with 7-fold symmetry. The five magnets operate in a vacuum and provide the magnetic field required to separate the incident beam electrons scattered from the target electrons (Møller scattering) and protons (elastic e-p scattering) in a liquid hydrogen target. The conceptual design was developed by the MOLLER Collaboration and was given to JLab in the form of amp turns and physical location, with additional physics requirements. This article presents prototyping of the coils and magnet support system and discusses the lessons learned during the process along with the plans for full magnet testing and installation. The JLab Magnet Group along with the MOLLER Collaboration developed the specification document that includes keep out zones to design the set of magnets. JLab contracted the design of the first toroid magnet (TM0) of the magnet system to Massachusetts Institute of Technology. The other four toroid magnets (TM1 through TM4) have been designed by JLab and are in the process of fabrication and assembly. Prototype coils of TM1-TM4 were fabricated by Everson-Tesla Incorporated, PA (USA). Finally, this article presents the unique challenges of the design, alignment, high current density, operating range, high radiation dose, and vacuum environment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗

DOE OSTI2023

Do Ionic Liquids Slow Down in Stages?

High impact recent articles have reported on the existence of a liquid–liquid (L–L) phase transition as a function of both pressure and temperature in ionic liquids (ILs) containing the popular trihexyltetradecylphosphonium cation (P 666,14 + ), sometimes referred to as the “universal liquifier”. The work presented here reports on the structural-dynamic pathway from liquid to glass of the most well-studied IL comprising the P 666,14 + cation. We present experimental and computational evidence that, on cooling, the path from the room-temperature liquid to the glass state is one of separate structural-dynamic changes. The first stage involves the slowdown of the charge network, while the apolar subcomponent is fully mobile. A second, separate stage entails the slowdown of the apolar domain. Whereas it is possible that these processes may be related to the liquid–liquid and glass transitions, more research is needed to establish this conclusively.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Influence of open ocean biogeochemistry on aerosol and clouds: Recent findings and perspectives

Aerosols and clouds are key components of the marine atmosphere, impacting the Earth’s radiative budget with a net cooling effect over the industrial era that counterbalances greenhouse gas warming, yet with an uncertain amplitude. Here we report recent advances in our understanding of how open ocean aerosol sources are modulated by ocean biogeochemistry and how they, in turn, shape cloud coverage and properties. We organize these findings in successive steps from ocean biogeochemical processes to particle formation by nucleation and sea spray emissions, further particle growth by condensation of gases, the potential to act as cloud condensation nuclei or ice nucleating particles, and finally, their effects on cloud formation, optical properties, and life cycle. We discuss how these processes may be impacted in a warming climate and the potential for ocean biogeochemistry—climate feedbacks through aerosols and clouds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Pressure safety approach for PIP-II cryogenic distribution system and cryomodules

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vesse l relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Precision measurement of the magnetic field in run 1 of the g-2 experiment

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vessel relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

AGR-5/6/7 Final Release-to-Birth Ratio Data Analysis

AGR-5/6/7 is the last of a series of Advanced Gas Reactor (AGR) experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) in support of development and qualification of tristructural isotropic (TRISO) low-enriched fuel for use in the high-temperature gas cooled reactor (HTGR). AGR configuration and irradiation conditions are based on prismatic HTGR technology that is distinguished primarily through the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO-coated particle fuel. The AGR tests provide valuable irradiation-performance data to support fuel process development, qualify fuel for normal operation and accident conditions, and support development and validation of fuel performance and fission-product (FP) transport models and codes. Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. Release-to-birth ratios (R/B) for fission-gas isotopes released from each capsule are calculated from release rates, measured by germanium detectors in the Fission Product Monitoring System (FPMS) installed downstream from each capsule, and birth rates calculated using numerical models of FG generation. The R/Bs are a critical measure of the ability of the fuel kernel, the particle coating layers, and the compact matrix to retain fission-gas atoms, preventing their release into the sweep-gas flow, and the impact of initially defective particles and/or particle-coating failures that occur during irradiation. For fission-gas isotopes, particle failure is defined as failure of all coating layers, allowing gaseous fission atoms to escape from a particle. During the first five cycles (162B ? 165A), R/Bs were stable in the 10-8?10-6 range, and no in-pile particle failures were observed, based on the gross gamma counts. The maximum R/B value of around 2 ? 10-6 for Kr-85m resulted from the presence of as-fabricated exposed kernels (based on the high exposed kernel fraction), the dispersed uranium, and high fuel particle temperatures in Capsule 1. Comparison of capsule-measured R/Bs from these early cycles to predictions using the previously developed AGR R/B model demonstrated FG release from the AGR-5/6/7 TRISO fuel was comparable to that of previous experiments. In addition, the Kr-85m R/B per-exposed-kernel values are comparable to R/B values obtained in AGR-3/4 irradiation experiment and four irradiation experiments performed during 1980s: (1) HRB-17/18, (2) COMEDIE-BD1, (3) HFR B1, and (4) HRB-21. In contrast, all measured R/B values are lower than predictions by the commonly used Richards and German models, which are intentionally conservative. A large number of in-pile particle failures occurred in Capsule 1 by the end of Cycle 166A. During the final four cycles (166A ? 168A), apparent damage to the Capsule 1 gas line appeared to cause FG leakage from that capsule into the other four capsules, resulting in an increase in fission gas (FG) detected in the effluent for all capsules. Isolation of the Capsule 1 gas line during the last three cycles also prevented measurement of its FG release. Thus, R/Bs in all capsules after Cycle 166A are highly uncertain because of undefined amount of leakage from Capsule 1, especially for long-lived isotopes. A few hundred in-pile particle failures were estimated for Capsule 1 before the end of Cycle 166A, but the total number of failures is unknown due to the lack of FG release data in the later cycles. Based primarily on evidence from the gross gamma counts during Cycle 168A, approximately 15 particles failed in Capsule 3 and four particles failed in Capsule 2. In-pile failures in Capsule 3 were anticipated because this capsule was designed to operate beyond the HTGR normal operating temperature range. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on the absence of the typical spikes in gross gamma counts and low failure estimates using the AGR model, developed in INL/EXT-14-32970, for R/B of the short-lived isotopes (Kr-89 and Xe-137) with minimal leakage from Capsule 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Pressure safety approach for PIP-II cryogenic distribution system and cryomodules

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vessel relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Computational Modeling, Simulation, and Potential Applications of Optical Stochastic Cooling

With the rising demand for intense particles beams, much research is being conducted in the area of particle beam cooling. One of these methods, called Stochastic Cooling (SC) (developed at CERN in the 1970's), delivered a feedback method to improve the quality and lifetime of circulating proton beams by reducing their 6D phase-space and has been widely implemented in a number of hadron machines. However, traditional stochastic cooling schemes are limited by the bandwidth of microwave frequency systems. Optical Stochastic Cooling (OSC) is a promising extension of the stochastic cooling beam cooling technique. OSC instead uses optical wavelengths which allows for improved control and increased cooling but creates its own technical challenges. This lays out work conducted toward the experimental demonstration of OSC at Fermilab's Integrable Optics Test Accelerator (IOTA) storage ring. This includes the design and characterization of parts of the optical delay system, the development and validation of a high-fidelity computational model of the OSC process, and the investigations into possible applications of the OSC mechanism to advanced beam manipulations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Remarkable enhancement in thermal conductivity of stainless-steel leveraging metal composite via laser powder bed fusion: 316L-Cu composite

This work aimed to significantly enhance the thermal conductivity of the SS matrix by fabricating a metal composite SS-Cu system utilizing a laser powder bed fusion process. The impact of process parameters on the microstructure of SS-Cu samples is investigated. Here, the processing and post-processing factors that influence the effective thermal conductivity of the metal composite SS-Cu system are discussed. Then, the effective thermal conductivity and mechanical properties of the SS-Cu composite are measured at room temperature and operating temperature of 300 °C. The optimized laser powder bed fusion parameters led to low Marangoni convection, high cooling rate, greater macro-segregation, and higher thermal conductivity values albeit at the cost of ductility in SS-Cu metal composite likely due to the presence of the unmelted powder particles. As-printed SS-60 vol % Cu and annealed SS-60 vol % Cu had 2.5X and 6.6X thermal conductivity enhancement, respectively. During annealing, phase separation, partial melting of Cu, and formation of a continuous network of Cu surrounding SS play a key role in enhancing the thermal conductivity of the SS-Cu metal composite structure.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

Individual Wave Detection and Tracking within a Rotating Detonation Engine through Computer Vision Object Detection applied to High-Speed Images

Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.

Johnson, Kristyn↗

DOE OSTI2025

Comparing XRISM Cluster Velocity Dispersions with Predictions from Cosmological Simulations: Are Feedback Models Too Ejective?

The dynamics of the intracluster medium (ICM), the hot plasma that fills galaxy clusters, are shaped by gravity-driven cluster mergers and feedback from supermassive black holes (SMBHs) in the cluster cores. XRISM measurements of ICM velocities in several clusters offer insights into these processes. We compare XRISM measurements for nine galaxy clusters (Virgo, Perseus, Centaurus, Hydra-A, PKS 0745–19, A2029, Coma, A2319, and Ophiuchus) with predictions from three state-of-the-art cosmological simulation suites, TNG-Cluster, the Three Hundred Project GADGET-X, and GIZMO-SIMBA, that employ different models of feedback. In cool cores, XRISM reveals systematically lower velocity dispersions than the simulations predict, with all 10 measurements below the median simulated values by a factor of 1.5–1.7 on average and all falling within the bottom 10% of the predicted distributions. The observed kinetic-to-total pressure ratio is also lower, with a median value of 2.2%, compared to the predicted 5.0%–6.5% for the three simulations. Outside the cool cores and in non-cool-core (NCC) clusters, simulations show better agreement with XRISM measurements, except for the outskirts of the relaxed, cool-core cluster A2029, which exhibits an exceptionally low kinetic pressure support (<1%), with none of the simulated systems in either of the three suites reaching such low levels. The NCC Coma and A2319 exhibit dispersions at the lower end but within the simulated spread. Our comparison suggests that the three numerical models may overestimate the kinetic effects of SMBH feedback in cluster cores. Additional XRISM observations of NCC clusters will clarify if there is a systematic tension in the gravity-dominated regime as well.

galaxy clusters↗

DOE OSTI2021

Individual Wave Detection and Tracking within a Rotating Detonation Engine through Computer Vision Object Detection applied to High-Speed Images

Johnson, Kristyn↗

DOE OSTI2022

Machine learning-based optimization of air-cooled heat sinks

Machine learning-based models using Artificial Neural Network (ANN) and greedy search algorithm are used to optimize air-cooled parallel plate-finned heat sinks (PPFHSs) subjected to laminar flow over an extensive range of design parameters. Here, the thermal and hydraulic performances of PPFHSs are represented by heat transfer coefficient (h) and pressure drop (ΔP), respectively. Optimization objectives for PPFHS designs can vary from industry to industry depending on their design priorities. The present study proposes a novel and generalized optimization method that defines practical optimization objectives and provides an accurate optimization process to design effective PPFHSs for a wide range of industrial applications with different design requirements. Three optimization objectives are presented in this study: (i) the largest h PΔ, (ii) the largest h within a specified maximum allowed flow rate, and (iii) the lowest weight that maximizes h for operation within the maximum allowed flow rate. While the shortcoming of the first objective is demonstrated, the other two objectives are found to be suitable for designing effective heat sinks (HSs) across different applications. Results suggest a promising trend from the third objective to develop HSs with ~ 37-68% lower weight, 80-85% reduced ΔP, and negligible penalty in h compared with optimized HSs obtained from the second objective. However, since the third objective leads to HSs with thinner fins, structural analysis should be performed to ensure reliable operation of the HSs.

42 ENGINEERING↗

DOE OSTI2026

Full-Strength Additive Manufacturing of Pure Aluminum Alloy 7075 Using Liquid Metal Jet Printing

Additive manufacturing (AM) of high-strength aluminum alloys can produce high-performance aerospace and biomedical parts with unique geometries and functions. However, many high-performance alloy systems such as AA 7075 pose challenges during AM due to thermal cycles that cause hot cracking, limiting AM to casting-friendly alloys. Stock AA 7075 alloy poses additional challenges due to the excessive growth of insoluble precipitates during casting and needs to be thermomechanically processed to obtain acceptable mechanical properties. Here, this work uses liquid metal jet printing to produce high-resolution pure AA 7075 parts at high strength (593 MPa) and ductility (10 pct), surpassing the specific strength of any type of steel. The parts are printed by depositing fine droplets on a high-temperature substrate that reduces cooling rates, preventing hot cracking. The results show that fully dense, crack-free parts can only be printed at substrate temperatures of 500 °C with a cooling rate threshold of < 2000 K/s to prevent cracking and < 500 K/s to eliminate porosity. The tensile strength, microhardness, and ductility of the parts matched or were beyond the wrought alloy specifications after heat treatment, possibly due to high solidification rates that produce finer insoluble precipitates. Our results indicate that many challenging high-performance alloy systems can now be 3D printed based on the critical thresholds for cooling rates using any AM method that can satisfy these conditions.

36 MATERIALS SCIENCE↗

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↗