Search

Heat transfer in jointed rock masses exhibits characteristics of discontinuous, heterogeneous, and anisotropic due to the effects of rock discontinuities. In this work, thermal contact resistance (TCR) evolution of granite joints under normal stress was investigated. A novel method was proposed for the measurement of TCR for rough rock joints. A numerical calculation program was established for the determination of contact area based on the initial morphology of the joint surfaces and the monitored joint closure during the experiment. An estimation model of TCR for granite joints based on the equivalent joint aperture was established. It is found that TCR of granite joints shows decreasing trend with normal stress. The decreasing rate is relatively large at the initial loading stage and gradually reduces with normal stress. The values of TCR eventually stabilize. TCR of granite joint has no obvious correlation with joint roughness coefficient. It is more significantly affected by the contact area. Plastic deformation and rupture of the micro-convex body on the joint surfaces occurred due to the application of normal stress, which improve the heat transfer efficiency and lead to the decrease in TCR. In general, predication results of the TCR estimation model show good consistency with the experimental data. The model is more suitable for the cases with small joint apertures. This study can provide knowledge to the better understanding of heat transfer in rock masses containing different kinds of discontinuities. Highlights: A novel method was proposed for the measurement of thermal contact resistance for rough rock joints. A numerical calculation program was established for the determination of contact area for rough rock joints. Evolution of thermal contact resistance with normal stress and joint contact area was obtained and analyzed. An estimation model of thermal contact resistance for granite joints based on the equivalent joint aperture was established. [ABSTRACT FROM AUTHOR]

In this work, effective thermal conductivity (λEff) of fractured rock masses was numerically investigated. A two-dimensional Discrete Fracture Network (DFN) model of the fractured rock masses was established based on the statistic results of natural fracture development in a potential area for high level radioactive waste disposal in China. Steady state heat transfer processes in the fractured granite rock masses were numerically simulated using finite element method (FEM). The calculated λEff values of the fractured granite rock masses in dry and saturated conditions are 1.99 W/(m K) and 2.31 W/(m K), respectively. Compared with the thermal conductivity of intact granite [λIntact, 2.5 W/(m K)], the drop rates are 20.4% and 7.6%, respectively. Sensitivity analysis was conducted on the main model parameters including fracture density (FDensity), trace length (FLength), thermal contact resistance (FTCR), and λIntact. The results indicate the relation between λEff and three fracture parameters (FDensity, Flength and FTCR) can be fitted using power law or negative exponent functions with good consistency. When fracture network parameters remain unchanged, λEff is in linear positive correlation to λIntact. The slop of the fitted line is determined by the fracture network parameters. Due to the fact that distribution of generated fractures in different directions are quite uniform, λEff did not show significant difference in different directions. On the basis of the above-mentioned results, an estimation model was proposed for the determination of λEff of fractured rock masses using P21 (total length of fracture traces per unit area), FTCR, and λIntact. The proposed estimation model shows good consistency to the calculated results of FEM model. [ABSTRACT FROM AUTHOR]

The natural fracture system exerts a significant influence on the physical and mechanical properties of rock masses. The heat transfer in fractured rock masses is much more complicated than that of intact rock due to the thermal contact resistance induced by natural fractures with varying size, aperture, roughness, and orientation. This work describes a method for determining the representative element volume (REV) of heat transfer and the 3D thermal conductivity tensor of fractured rock masses. The method is based on the element-wise averaging of temperature and heat flux, which can be obtained from a finite element solution to a steady state heat transfer problem with specific boundary conditions. A case study was conducted based on the geological data of a potential area for high level radioactive waste disposal in China. The results indicate that the size of heat transfer REV of the rock masses in the study area is approximately 18 m × 18 m × 18 m. 3D thermal conductivity tensor, the magnitudes and directions of the three principal thermal conductivities, and anisotropy of the fractured rock masses were obtained. Influences of main fracture parameters including density, disc diameter, and thermal contact resistance on the heat transfer characteristics were investigated. The research results in this work can provide a better understanding of the discontinuous heat transfer processes in fractured rock masses and their homogenization methods. [ABSTRACT FROM AUTHOR]

Convection heat transfer characteristics in a single rock fracture was experimentally investigated in this work. Distilled water was used as working medium for heat transfer in manmade smooth and rough fractures of sandstone and granite samples. The influence of fracture surface roughness on heat transfer intensity was analyzed. During seepage heat transfer process, evolution of samples external surface temperatures was monitored using six resistance temperature detectors to study the local heat transfer characteristics along the flow direction. Subsequently, heat transfer correlations were proposed on the basis of the experimental results. The results indicate that the roughness of rock fracture surface improves the overall heat transfer intensity to a certain extent, whereas lithology has little influence on the heat transfer intensity. For smooth fractures, the local heat transfer intensity decreases with the increase of the distance from the inlet. For rough fractures, evolution of sample external temperatures along the flow direction exhibits strong nonlinear characteristics with multiple peaks. Nu/Pr 1/3 is close to a power function of Re under the experimental conditions in this work. With the increase of experimental temperature, the relationship gradually tends to be linear. [ABSTRACT FROM AUTHOR]