Next Article in Journal
Fully Superconducting Josephson Bolometers for Gigahertz Astronomy
Next Article in Special Issue
Reuse of Decommissioned Hydrocarbon Wells in Italian Oilfields by Means of a Closed-Loop Geothermal System
Previous Article in Journal
2-D Cross-Plot Model Analysis Using Integrated Geophysical Methods for Landslides Assessment
Previous Article in Special Issue
Well-Doublets: A First-Order Assessment of Geothermal SedHeat Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation on the Influence of Surface Flow Direction on the Heat Transfer Characteristics in a Granite Single Fracture

1
College of Construction Engineering, Jilin University, Changchun 130012, China
2
Key Lab of Groundwater Resource and Environment, Ministry of Education, Jilin University, Changchun 130026, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(2), 751; https://doi.org/10.3390/app11020751
Submission received: 2 December 2020 / Revised: 7 January 2021 / Accepted: 11 January 2021 / Published: 14 January 2021
(This article belongs to the Special Issue New Trends in Enhanced, Hybrid and Integrated Geothermal Systems)

Abstract

:
As fluid passes through the fracture of an enhanced geothermal system, the flow direction exhibits distinct angular relationships with the geometric profile of the rough fracture. This will inevitably affect the heat transfer characteristics in the fracture. Therefore, we established a hydro-thermal coupling model to study the influence of the fluid flow direction on the heat transfer characteristics of granite single fractures and the accuracy of the numerical model was verified by experiments. Results demonstrate a strong correlation between the distribution of the local heat transfer coefficient and the fracture morphology. A change in the flow direction is likely to alter the transfer coefficient value and does not affect the distribution characteristics along the flow path. Increasing injection flow rate has an enhanced effect. Although the heat transfer capacity in the fractured increases with the flow rate, a sharp decline in the heat extraction rate and the total heat transfer coefficient is also observed. Furthermore, the model with the smooth fracture surface in the flow direction exhibits a higher heat transfer capacity compared to that of the fracture model with varying roughness. This is attributed to the presence of fluid deflection and dominant channels.

1. Introduction

Geothermal energy is a clean and environmentally friendly renewable energy source with a wide distribution range, large reserves, and long duration [1,2,3]. The high temperature rock mass buried 3–10 km underground is characterized by low permeability and low porosity [4]. In order to extract thermal energy from such high temperature rock mass, several countries (led by the United States) proposed the use of artificial hydraulic fracturing, which led to the development of reservoir-transforming hydro-shearing techniques to establish fracture network channels and improve the heat transfer capacity [5]. This geothermal project is often referred to as an enhanced geothermal system (EGS). The fracture surface formed by hydraulic fracturing is typically rough, which directly affects the heat transfer performance of the working medium (water or CO2) [6].
The key to extracting geothermal resources by EGS is the flow and heat exchange process of working fluid in fractures of high temperature rock mass [7]. However, when the fluid begins to flow from the injection well into the fracture, the initial direction of flow often has an angle equal to the direction of the rough fracture surface, which inevitably affects the heat exchange effect. The accurate determination of the convection heat transfer coefficient between the fluid and fracture is crucial to the optimization of the reservoir transformation and productivity predictions [8]. In addition, a clear heat transfer law within the fractures is of great significance to the establishment of EGS heat recovery models [9].
Numerous studies have attempted to quantify the roughness of fracture surfaces. For example, Patton [10] proposed the concept of the wave angle and established the relationship between wave angle and fracture morphology. Mandelbrot [11] proposed the concept of fractal geometry and developed a fractal dimension framework to describe the geometrical characteristics of fracture rough surfaces. Barton [12] used the joint roughness coefficient (JRC) to characterize the section geometry of fracture surfaces to determine 10 fracture types with an approximate length of 10 cm. Xie [13] demonstrated the fractal dimension to be the measured values of the joint roughness coefficients.
Experimental and numerical simulations have demonstrated the fluid flow and heat transfer processes in fractures to be affected by many factors, including the aperture, roughness, type of fluid, injection flow rate, initial temperature of rock, and rock thermophysical parameters [14]. Early scholars used the plate model to describe the heat transfer characteristics in fracture. However, Tsang and Brown [15,16] pointed out that the parallel plate theory may cause flow estimation errors of 1–2 orders of magnitude. In an experimental study on the seepage heat transfer in a single fracture of granite, Zhao and Tso [17] revealed the single fracture with a certain roughness to have a better heat transfer performance than the smooth-plate fracture. Similarly, the two-dimensional coupled heat transfer model of a single fracture established by He [18] also demonstrated a strong correlation between the local convection heat transfer coefficient and the fracture profile. Neuville [19] established a self-affine rough fracture coupling model that assumes a reduced heat transfer efficiency in a rough fracture with constant temperature due to the channel effect. Research carried out by Huang [20] also showed that a larger relative roughness increases the flow friction and significantly reduces the heat exchange. Luo [21] investigated the influence of initial rock temperature on the heat transfer efficiency of fractures, demonstrating a positive relationship between the heat transfer rate and rock temperature. The two-dimensional single-fracture heat transfer model established by He [22] shows that the effect of fracture surface roughness on heat transfer intensity decreases with an increasing injection flow rate. Confining pressure has a significant impact on the fracture aperture. Shu [23] employed a novel experimental device to simulate the evolution of the hydraulic and heat transfer characteristics of fractures under a confining pressure. Results show that both the heat transfer rate and total heat transfer coefficient decrease with an increasing confining pressure. Zhang [24] established a three-dimensional numerical model based on three-dimensional laser scanning to evaluate the effects of rock temperature, water flow velocity, roughness, and fracture aperture on the heat transfer coefficient. Simulation tests reveal the dominant influence of the flow rate on the rock roughness, followed by the aperture size. Andrade et al. [25] simulated the flow and heat transfer between two-dimensional parallel rough surfaces, demonstrating the minimal effect of wall roughness on the heat transfer for the characteristic parameter Pe = u(d/2)2/L < 200. Zhang [26] performed convection heat transfer tests of carbon dioxide in a single fracture to establish a numerical model and proposed the rough fracture channel effect. Experimental and numerical simulations can increase our understanding of the geothermal field model. Fox et al. [27] established a discrete fracture network and investigated the relationship between hydraulic aperture and fluid runoff, concluding that variations in the fracture aperture can cause flow channeling and reduce the heat transfer area. Furthermore, several studies [28,29] have employed Monte-Carlo [30,31] stochastic simulations to build fracture network models with different geometric parameters, and in particular, the EGS site model to predict the heat recovery performance. Bruel [32] integrated experimental data into a single model used to generate an estimation model of the main hydraulic parameters. Chen [33] established a geothermal reservoir model with a rough fracture surface based on small-scale research results, revealing the constant heat transfer coefficient (HTC) recommended in previous studies to underestimate the final outlet fluid temperature in the case of rough fractures. These studies demonstrate that research on the influence of fracture morphology on the fluid flow and heat transfer process at the experimental scale can be applied to the actual site model [28,34].
In summary, previous studies have shown that the geometry of the fracture surface has a profound influence on the fluid flow and heat transfer process in the fracture [35,36]. According to the different influencing factors, many scholars [23,37,38] have designed experiments and carried out corresponding numerical simulations to determine how these influence factors affect the heat transfer process in fractures. The results of such studies are very important for the development of EGS stimulation technology.
The aforementioned research generally focuses on the influence of fracture morphology on fluid heat transfer, while studies on the effect of the angle between the direction of the fluid flow and the morphology of the fracture wall on the heat transfer performance are limited. In addition, because of the geological history, the fracture aperture of EGS reservoirs range from the micron to centimeter scale [39]. This makes it difficult to obtain the temperature distribution of rock fractures and thermal mediums under laboratory conditions. Hence, in the current study, we established a single fracture heat transfer model with a random geometry profile that was subsequently verified by experiments. Four cases with fracture profiles and varying angles between flow direction were set up to simulate and explore the heat transfer performance of distilled water through fractures. The influence of the fluid flow direction on the fracture heat transfer characteristics was then discussed.

2. Numerical Simulations

2.1. Mathematical and Physical Models

A three-dimensional single fracture model was established to study the influence of flow direction on the heat transfer characteristics in a granite single fracture. The geometry profile of the fracture was constructed using a random method (Figure 1a), while the rough fracture wall was built by stretching in the vertical direction (Figure 1b,c).
The surface flow of the fracture surface was assumed to have four different directions, with the 0° direction taken as the z-axis direction (Figure 1b). The angle between the surface flow direction and the positive direction of the z-axis is denoted as α. The surface flow direction had a gradient of 30° and was rotated counterclockwise around the z-axis by 30°, 60°, and 90° to obtain the fracture wall. The rough single-fracture numerical model shown in Figure 2 was obtained by embedding these fracture walls into a cylindrical model with a size of 100 × 50 mm (length × diameter). The vertical displacement of the fracture in the model was determined as 0.3 mm by translation method. Table 1 lists the geometric parameters of the four fracture models. The surface area and inlet area deviations of the four fractures were 1.87% and 0.33%, respectively. Therefore, at the scale of this study, the effects on the simulation results due to differences in surface area and inlet area are negligible.
For natural fracture, the Reynolds number is much smaller than the critical Reynolds number (Re < 1800 [35]) due to its small aperture and low flow rate of underground fluid. Thus, the fluid flow in a granite fracture remains laminar. Convective heat transfer in three-dimensional rock fractures was modeled using FLUENT 6.3 [40] by employing a finite volume method (FVM) to solve the Navier–Stokes equations.
The governing equations for the conservation of momentum, mass, and energy in the fracture are described as follows:
Continuity equation:
ρ w t + · ( ρ w V ) = 0
Momentum equation:
( ρ w V ) t + ( ρ w V V ) = P + · ( μ [ ( V + V T ) 2 3 · V I ] )
Energy equation:
t ( ρ w E ) + · ( V ( ρ w E + p ) ) = · ( λ w T + μ [ ( V + V T ) 2 3 · V I ] · V )
where ρ w (kg/m3) and ρ r (kg/m3) are the densities of water and rock, respectively; V (m/s) is the velocity vector; p (Pa) is the pressure; I is the unit tensor; E is the total energy; λ w (W/(m·°C)) is the fluid thermal conductivity; μ (Pa s) is the fluid viscosity.
In EGS frameworks, water is typically used as a heat exchange medium. However, the physical properties of water are more susceptible to temperature than pressure. Hence, the relationship between the physical properties of water and temperature was described by the following empirical formula [41]:
Dynamic viscosity of water:
μ w = 1.3799 0.0212 T + 1.3604 × 10 4 T 2 4.6454 × 10 7 T 3 + 8.9043 × 10 10 T 4 9.0791 × 10 13 T 5 + 3.8457 × 10 16 T 6
Specific heat capacity of water:
C p , w = 12010 80.4 × T + 0.3 × T 2 5.4 × 10 4 T 3 + 3.6 × 10 7 T 4
Thermal conductivity of water:
λ w = 0.8691 + 8.9 × 10 3 T 1.5837 × 10 5 T 2 + 7.9754 × 10 9 T 3
Water density:
ρ w = 838.4661 + 1.4005 × T 3 × 10 3 T 2 + 3.7182 × 10 7 T 3
where T (°C) is the water temperature and which ranges from 10 to 100 °C.
The fracture inlet and outlet employ mass flow rate and outflow boundary conditions, respectively. The outer surface of the rock is considered to be an adiabatic boundary, while the inner wall of the fracture is considered as an impervious and no-slip boundary.

2.2. Model Meshing and Mesh Independence Verification

In order to focus on fluid flow and heat transfer in fractures, we employed structured hexahedral grids in fracture spaces (Figure 3a) and an increased grid density around fractures (Figure 3b,c), while an unstructured tetrahedral grid was used for the computational domain outside the fracture. Furthermore, we selected α = 0° as the verification sample to verify that the simulation results are not affected by the number of meshes. Table 2 lists the four grid numbers for α = 0° and Figure 3 presents the grid characteristics of the fourth grid number.
Under the same simulation conditions, four models with different numbers of grids were constructed and the average temperature of the fluid at the outlet was monitored (Figure 4). Mesh case 1 is distorted due to the presence of grids with a high value and aspect ratio. When the number of meshes exceeds 1.3 million, the outlet temperature of different mesh cases is generally consistent with time. In order to accurately reflect the geometry of the fracture surface and save computer resources, we selected mesh case 3 for further numerical calculations. Table 3 details the grid numbers for the four models.

2.3. Initial Conditions and Simulations

The initial temperature of the rock was set to 90 and 70 °C to prevent a water phase transformation. The initial flow rates were 10, 20, 30, and 40 mL/min, respectively, depending on the performance of the plunger pump. The injection temperature of water was 25 °C without considering the change of room temperature. Table 4 lists the simulation conditions.
Calculations were performed following the unsteady double precision method and the second-order implicit mode to improve calculation accuracy. The SIMPLEC (semi-implicit method for pressure-linked equations consistent) algorithm, an improved SIMPLE algorithm, was adopted to amend the fluid pressure and flow rate, while the PRESTO (PREssure STaggering Option) scheme was used for the pressure discretization. The second order discretization scheme was used for the convection terms. The residual convergence criteria for the computational equations was set to 10−5 and the time step size for the transient simulations time was equal to 1 s. The total solution time was taken as 30 min.

3. Model Verification

Although the establishment of a realistic fracture model proves to be a complicated task, a single-fracture granite model with a flat and smooth surface can validate the proposed model. Thus, experimental and numerical simulation results for smooth-plate fractures are compared.
A fracture heat transfer laboratory simulation system was independently developed by our research group for the experiments. Figure 5 depicts the experimental system, which can be divided into five key components including the seepage, confining pressure, heating, core holder, and data measurement systems. Six temperature sensor mounting holes were preset on the top of the holder to ensure direct contact between the sensor and the outer surface of the rock sample. The data measurement system monitored the rock and fluid temperatures on the outer surface of the core and at the core outlet in real time, respectively.
The rock sample with smooth-plate fractures shown in Figure 6a,b was tested to verify the correctness of the numerical model, while a numerical model (Figure 6c) was established to simulate the test results. The rocks used in this study were obtained from important geothermal resource targets in the northern Songliao Basin of China. The granite was processed into a cylindrical specimen with a diameter and height of 50 and 100 mm, respectively. The thermophysical properties of the granite were obtained from tests. The specific heat capacity measuring device (BBR series, made in Xiangtan Xiangyi Instrument Co., Ltd., Xiangtan, China) based on cooling method and the thermal conductivity scanner (TCS, made in Lippmann and Rauen GbR, Schaufling, Germany) based on optical scanning principle are used to measure the specific heat capacity and thermal conductivity of rock, respectively. Figure 7 presents the variations in specific heat capacity and the thermal conductivity of rock with temperature.
The experimental flow rates were set as 5, 10, 20, and 30 mL/min for temperatures of 60, 70, 80, and 90 °C, respectively. Both the boundary and initial conditions followed those detailed in Section 2.3. Numerical simulations were performed based on the experiments to ensure that the results were not affected by other factors. The unsteady numerical simulation of fluid temperature Toutlet at the fracture outlet for a given time was calculated as follows:
T o u t l e t = A o u t T w ( x , y , z ) μ w ( x , y , z ) ρ w ( x , y , z ) c p , w ( x , y , z ) d A A o u t μ w ( x , y , z ) ρ w ( x , y , z ) c p , w ( x , y , z ) d A
where Aout (m2) is the area of the fracture outlet.
Figure 8 clearly describes the comparison between the numerical simulation and experimental results for a smooth single fracture. Although there are some errors in the outlet temperature and the rock surface temperature obtained from the experiments and numerical simulation, the maximum errors are only 0.8% and 1.3%, respectively. This indicates that the numerical model established using FLUENT 6.3 can accurately reflect the convection heat of single-fracture granite.

4. Data Processing

4.1. Heat Extraction Rate

The heat extraction rate can measure the heat absorbed by water flowing through a rock fracture per unit time and is calculated as follows [26]:
Q = ( H o u t l e t H i n l e t ) M ,
where M (kg/s) is the flow rate; H i n l e t (kJ/kg) and H outlet   (kJ/kg) are the water enthalpy at the inlet and outlet, respectively. The enthalpy of water can be determined by numerical simulations of the inlet and outlet temperature and pressure.

4.2. Total Heat Transfer Coefficient

For fracture channels with complex surface morphology, the total heat transfer coefficient can describe the overall convective heat transfer performance within the fracture. The total heat transfer coefficient can be calculated from Newton’s law of cooling [42]:
h t o t a l = q f T f T w ,
where q f (W/m2) is the average heat flux on the inner wall of the fracture; T f (°C) is the average temperature on the inner wall of the fracture; T w (°C) is the average temperature of the water in the fracture. These parameters can be calculated from the numerical simulation results as follows:
q f = A f q ( x , y , z ) d A A f d A ,
T f = A f T ( x , y , z ) d A A f d A ,
  T w = V f T ( x , y , z ) ρ ( x , y , z ) C p ( x , y , z ) u ( x , y , z ) d V V f ρ ( x , y , z ) C p   ( x , y , z ) u ( x , y , z ) d V ,
where A f   (m2) is the surface area of the fracture; V f   (m3) is the volume of the fracture.

4.3. Local Heat Transfer Coefficient

The local heat transfer coefficient can describe the local heat transfer characteristics in fractures and is an important parameter reflecting the heat transfer characteristics in the flow direction. Here, we calculated the discrete form of the local heat transfer coefficient via Equation (10) as follows:
h x = q f x T f x T w x ,
where q f x (W/m2), T f x (°C), and T w x (°C) are the heat flux, temperature, and fluid temperature at position x of the fracture wall, respectively. Furthermore, the Reynolds number (Re) can be calculated using the following formula for seepage in single fracture of rock:
R e = 2 M d μ w ,
where M (kg/s) is the mass flow rate and d (m) is the fracture aperture.
The overall Nusselt number (Nu) can be determined as follows:
N u = h t o t a l d e λ w ,
where d e (m) is the equivalent hydraulic diameter of a channel with arbitrary shape.

5. Results and Discussion

5.1. Temperature Field Distribution

Table 5 lists the outlet water temperature determined from the 32 simulations. The outlet water temperature is maximized at α = 90°, followed by models α = 0°, α = 30°, and α = 60°.
Figure 9 depicts the distribution characteristics of the temperature fields on the fracture surfaces of four different model cases at several injection flow rates. With the injection of low temperature water, the temperature inside the fracture begins to decrease and the cold front pushes towards the outlet. This phenomenon is consistent with previous studies [37,43,44]. The temperature at the center of the fracture typically exceeds that on both sides of the fracture. This is attributed to the greater thermal replenishment at the center of the fracture by the cylinder model compared to the two sides. Increasing the flow rate enhances the rate of the temperature reduction in the fracture. The temperature on the fractured surface of model α = 0° is distributed symmetrically along the z-axis due to the symmetry of the model. At α = 30° and α = 60°, the fracture surface temperature in the positive direction of the x-axis exceeds that in the negative direction. This indicates that the angle between the fracture morphology and flow direction has an important influence on the convection heat transfer process.
Figure 10 presents the local temperature distribution along the z-axis at x = 0.02 m and x = −0.02 m. At α = 30° and α = 60°, the local wall temperature at x = 0.02 m is lower than that at x = −0.02 m. This is related to the migration of the channel formed by the fracture profile towards the x-axis, which makes the fracture geometry profile no longer perpendicular to the flow direction. The fluid preferentially flows through the fractures in an easily accessible path where the flow direction is deflected. In addition, Figure 10 shows that the temperature curve at low flow rates (10 mL/min) reaches its peak value, while the temperature curve at high flow rates (30 mL/min) does not. This indicates that the cryogenic fluid injected at high flow rates penetrates the fractures more quickly and results in a thermal breakthrough in a short time, which causes rapid wall temperature reduction and no timely thermal replenishment. Fluids with low flow rates are more likely to utilize matrix heat transfer to compensate for the decrease in wall temperature caused by convection heat transfer.
Figure 11 illustrates the streamline distribution characteristics of the models with α = 30° and α = 60°. The flow direction is observed to slightly deflect in the fracture towards the x-axis, which results in the temperature difference shown in Figure 9 and Figure 10. This deflection is attributed to the presence of a gently inclined microfractured surface in front of the flow field, preventing the flow from passing. Figure 12 simplifies fracture surface morphology to account for the deflection of the fluid flow direction. A fractured surface consists of a limited number of microfractured surfaces. When the direction of the fluid injection is not parallel to the microfracture dip, the fluid will slightly deflect towards a blunt angle between the flow direction and the microfracture surface strike.

5.2. Heat Extraction Rate

The average Reynolds and Nusselt numbers in the fracture are calculated based on Equations (15) and (16), respectively (Figure 13). Based on the simulation conditions, the Reynolds number calculated using Equation (16) ranges from 14.5 to 53.3, which is much lower than the critical Reynolds number leading to the change of laminar flow state [35]. The results of Reynolds number calculation also show that the assumption of using laminar flow model in Section 2.1 is reasonable. The difference of thermophysical properties of water results in a greater Reynolds number at initial temperature T0 = 90 °C compared to that with initial temperature T0 = 70 °C. The maximum Nusselt number at α = 90° is 15.8 (Figure 13b).
Considering the small scale of the model in this study, we used the cumulative heat extraction rate to represent the heat absorbed by water from the model. The cumulative heat extraction rate can be calculated as follows:
Q s u m = ( H o u t l e t , t H i n l e t , t   ) M .
Figure 14 presents the cumulative heat extraction rate Qsum (W) under different operating conditions. The cumulative heat extraction rate is maximized at α = 90°, followed by α = 0°, α = 30°, and α = 60°. Within the local range, increasing the roughness of the fracture does not necessarily enhance its heat transfer capacity. This is inconsistent with the results of other studies because other researchers use two-dimensional models [22,45]. Moreover, although the cumulative heat extraction increases with the injection flow rate under the same model [24,43], the increasing rate (ΔQsum, see Figure 14a) will gradually decrease. Thus, blindly increasing the injection flow rate during EGS operation may be detrimental to heat recovery. At the same flow rate level, an increase of the initial temperature from 70 to 90 °C will enhance the cumulative heat extraction rate by approximately 45%. Figure 15 demonstrates the variation of heat extraction rate with time, indicating that a larger initial flow rate at the initial stage of the heat transfer will improve the heat extraction capacity. The heat loss on the fracture wall caused by high convection heat exchange caused by high flow rate is much greater than that supplied by the matrix by heat conduction, which causes a sharp decrease in heat extraction rate. The higher the flow rate, the faster the descent rate.

5.3. Total Heat Transfer Coefficient

Figure 16 shows the variation of the total convective heat transfer coefficient with time. The total heat transfer coefficient curve can be divided into increasing and decreasing stages. During the increasing stage, the low temperature water injection causes severe convection heat exchange on the fracture wall, which can prevent cold water from penetrating the whole fracture. Therefore, the total heat transfer coefficient increases rapidly with the injection of fluid until it reaches its peak value. At the decreasing stage, the low temperature water enters the fractures continuously and the heating effect of the rock matrix begins to weaken with the decrease of temperature. In addition, the fluid and fracture wall temperatures further decrease, and the low temperature water begins to penetrate the fractures. Although increasing the initial flow rate and temperature can amplify the total heat transfer coefficient, increasing the flow rate will shorten the time required by the total heat transfer coefficient to reach its peak value. This is obviously not conducive to the sustainability of heat recovery. Figure 16 shows that the total heat transfer coefficient in descending order is α = 90°, α = 0°, α = 30°, and α = 60°.

5.4. Local Heat Transfer Coefficient

Figure 17 depicts the local heat transfer coefficients of the rock fracture surfaces. The local convection heat transfer coefficient is observed to gradually decrease along the outlet direction, which is consistent with the previous study [22,42]. This is attributed to the gradual decrease in the temperature difference between the fluid and fracture wall along the outlet direction. The fracture model shows that the fracture roughness of α = 30° and α = 60° are different in the three axes, while that of α = 0° and α = 90° are different in the two axes. The flow direction on the fractured surface is slightly deflection, and consequently the local convection heat transfer coefficient in the negative direction of the x-axis exceeds that of the positive direction (Figure 11 and Figure 12). In addition, the distribution of the local heat transfer coefficient is strongly correlated with the rough fracture profile [45]. Figure 18 shows the distribution of the local heat transfer coefficients on the fracture profile at x = 0 m. Figure 18 shows that the model of fracture roughness from large to small at x = 0 m is α = 0°, α = 30°, α = 60°, and α = 90°. The fluctuation of local convection heat transfer coefficient increases with the increase of fracture roughness. At α = 90°, the local heat transfer coefficients exhibit a logarithmic distribution in the smooth section along the z-axis. Similar to the fracture profile, the distribution curve of the local heat transfer coefficient is characterized by peaks and troughs. The local heat transfer coefficient curves at α = 0°, α = 30°, and α = 60° reveal that the local heat transfer coefficient at the peaks of the fracture profile is consistently greater than that at the troughs. Moreover, increasing the injection flow rate not only improves the local heat transfer coefficient, but also has an enhancing effect that makes the difference between the peaks and troughs of the local heat transfer coefficient more evident. The enhanced effect can be observed via the increased slope of the local heat transfer coefficient curve. Figure 19 shows the enhancement effect of increasing the flow rate on the convective heat transfer intensity.
Figure 20 presents the distribution of the local heat transfer coefficient, wall temperature, and average flow rate on the fracture profile at z = 0.05 m. The curves of the average flow rate and local heat exchange coefficient become increasingly complex as α increases from 0° to 90°. Since the fracture model is constructed by translating along the y-axis, the distribution of fracture aperture in each direction of the model is not uniform. Therefore, a position with a higher absolute gradient results in a smaller fracture aperture, which results in a low flow rate as the water flows through the path of the small aperture. In contrast, fractures with larger aperture will form a dominant channel where the flow rate is greater. Compared with the temperature, the difference of flow rates at fractures with different apertures is the main factor causing uneven distribution of local heat transfer coefficient. In addition, a strong positive correlation is observed between the local heat exchange coefficient and the flow rate distribution, which indicates that the flow rate is the main factor affecting the heat transfer intensity.

5.5. Limitations of the Present Study

The rough fracture surface used in this study is actually a pseudo-3D fracture surface that ignores the anisotropy of roughness to a certain extent. Therefore, models for this study are not common in actual natural environments. Moreover, the scale of the fracture model may weaken the impact of the proposed flow direction deflection on the heat transfer effect. The heat transfer characteristics of inhomogeneous three-dimensional rough fractures need to be further studied. However, the focus of this work is to study how different angles between flow direction and fracture profile affect the heat transfer process in fractures. In fact, a certain size fracture surface is composed of a limited number of microfractures with occurrence, which have different angles with the flow direction of fluid. From this point of view, this work can actually provide basic knowledge for EGS research.

6. Conclusions

In this study, a hydro-thermal coupling model was developed to investigate the influence of fluid flow direction on the heat transfer characteristics of fractures. The accuracy of the proposed numerical model was verified by laminar convection heat transfer tests of smooth-plate fractured rock samples. A detailed analysis of the temperature field distribution, total heat transfer coefficient, and local heat transfer coefficient of the model were provided based on the numerical simulations. The main conclusions are as follows:
(1) At α = 30° and α = 60°, the temperature of the fracture surface in the negative direction of the x-axis is higher than that in the positive direction. This is because the flow direction of the fluid is not parallel to the orientation of the microfracture surface that prevents the fluid from flowing forward, resulting in fluid deflection.
(2) When the initial temperature is raised from 70 to 90 °C, the cumulative heat extraction rate Qsum increases the yield by 45%. Although the cumulative heat extraction increases with the injection flow rate, the rate of increase gradually decreases. An increase in the injection flow rate from 10 to 40 mL/min enhances the peak heat extraction rate. However, the increase in flow rate speeds up the rate of reduction of the heat extraction rate.
(3) Variations in the total heat transfer coefficient with time can be divided into increasing and decreasing stages. The total heat transfer coefficients decrease with α. This indicates that the presence of rough fractured surfaces in the forward direction of the fluid does not necessarily increase the convection heat transfer intensity between the fluid and the fractured surfaces. Increasing the flow rate and temperature can enhance the total heat transfer coefficient. However, a higher flow rate will shorten the time required by the total heat transfer coefficient to reach its peak value, while also increasing the rate of reduction of the total heat transfer coefficient in the decreasing stage.
(4) The local heat transfer coefficient gradually decreases along the direction of the water flow, which is consistent with the conclusions drawn by Bai [46]. An increasing flow rate not only increases the local heat transfer coefficient, but also has an enhancing effect that highlights the influence of the fracture profile on the local heat transfer coefficient distribution.
(5) There is a consistently strong correlation between the local heat transfer coefficient and the fracture profile, independent of the axis. In particular, although a change in the angle between the flow direction and the fracture profile can affect the local heat transfer coefficient, it does not affect the distribution characteristics of the local heat transfer coefficient.

Author Contributions

Conceptualization, writing—original draft preparation, X.G.; supervision and review, Y.Z.; formal analysis, data curation, Z.H.; software, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Nos. 41772238), by the New Energy Program of Jilin Province (Nos. SXGJSF2017-5), and the Interdisciplinary Fund of Jilin University (Nos. 101832020DJX072).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the first author or the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fridleifsson, I.B. Geothermal energy for the benefit of the people. Renew. Sustain. Energy Rev. 2001, 5, 299–312. [Google Scholar] [CrossRef] [Green Version]
  2. Mahbaz, S.B.; Dehghani-Sanij, A.R.; Dusseault, M.B.; Nathwani, J.S. Enhanced and integrated geothermal systems for sustainable development of Canada’s northern communities. Sustain. Energy Technol. Assess. 2020, 37, 100565. [Google Scholar] [CrossRef]
  3. Kinney, C.; Dehghani-Sanij, A.; Mahbaz, S.; Dusseault, M.B.; Nathwani, J.S.; Fraser, R.A. Geothermal Energy for Sustainable Food Production in Canada’s Remote Northern Communities. Energies 2019, 12, 4058. [Google Scholar] [CrossRef] [Green Version]
  4. Pettitt, R.A.; Becker, N.M. Mining Earth’s Heat: Development of Hot Dry Rock Geothermal Reservoirs. Geotherm. Energy 1983, 11. [Google Scholar]
  5. Olasolo, P.; Juárez, M.C.; Morales, M.P.; D´Amico, S.; Liarte, I.A. Enhanced geothermal systems (EGS): A review. Renew. Sustain. Energy Rev. 2016, 56, 133–144. [Google Scholar] [CrossRef]
  6. Kazemi, A.R.; Mahbaz, S.B.; Dehghani-Sanij, A.R.; Dusseault, M.B.; Fraser, R. Performance Evaluation of an Enhanced Geothermal System in the Western Canada Sedimentary Basin. Renew. Sustain. Energy Rev. 2019, 113, 109278. [Google Scholar] [CrossRef]
  7. Kolditz, O. Modeling flow and heat-transfer in fractured rocks - conceptual-model of a 3-d deterministic fracture network. Geothermics 1995, 24, 451–470. [Google Scholar] [CrossRef]
  8. Jiang, F.; Chen, J.; Huang, W.; Luo, L. A three-dimensional transient model for EGS subsurface thermo-hydraulic process. Energy 2014, 72, 300–310. [Google Scholar] [CrossRef]
  9. Shaik, A.R.; Rahman, S.S.; Tran, N.H.; Tran, T. Numerical simulation of Fluid-Rock coupling heat transfer in naturally fractured geothermal system. Appl. Therm. Eng. 2011, 31, 1600–1606. [Google Scholar] [CrossRef]
  10. Patton, F. Multiple modes of shear failure in rock. Proceeding Congress Int. Soc. 1966, 1, 509–513. [Google Scholar]
  11. Mandelbrot, B.B.; Wheeler, J.A. The Fractal Geometry of Nature. Am. J. Phys. 1983, 51, 286–287. [Google Scholar] [CrossRef]
  12. Barton, N. Review of a new shear-strength criterion for rock joints. Eng. Geol. 1973, 7, 287–332. [Google Scholar] [CrossRef]
  13. Xie, H.; Pariseau, W.G. Fractal estimation of joint roughness coefficients. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1995, 32, 266A. [Google Scholar] [CrossRef] [Green Version]
  14. Ogino, F.; Yamamura, M.; Fukuda, T. Heat transfer from hot dry rock to waterflowingthrough a circular fracture. Geothermics 1999, 28, 21–44. [Google Scholar] [CrossRef]
  15. Tsang, Y.W. The Effect of Tortuosity on Fluid Flow Through a Single Fracture. Water Resour. Res. 1984, 20, 1209–1215. [Google Scholar] [CrossRef]
  16. Brown, S. Fluid-flow through rock joints - the effect of surface-roughness. J. Geophys. Res. Solid Earth Planets 1987, 92, 1337–1347. [Google Scholar] [CrossRef]
  17. Zhao, J.; Tso, C.P. Heat-transfer by water-flow in rock fractures and the application to hot dry rock geothermal systems. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1993, 30, 633–641. [Google Scholar] [CrossRef]
  18. He, Y.; Bai, B.; Li, X. Comparative Investigation on the Heat Transfer Characteristics of Gaseous CO2 and Gaseous Water Flowing Through a Single Granite Fracture. Int. J. Thermophys. 2017, 38, 170. [Google Scholar] [CrossRef]
  19. Neuville, A.; Toussaint, R.; Schmittbuhl, J. Hydrothermal coupling in a self-affine rough fracture. Phys. Rev. E 2010, 82, 036317. [Google Scholar] [CrossRef] [Green Version]
  20. Huang, X.; Zhu, J.; Li, J.; Bai, B.; Zhang, G. Fluid friction and heat transfer through a single rough fracture in granitic rock under confining pressure. Int. Commun. Heat Mass Transf. 2016, 75, 78–85. [Google Scholar] [CrossRef]
  21. Luo, J.; Zhu, Y.; Guo, Q.; Tan, L.; Zhuang, Y.; Liu, M.; Zhang, C.; Xiang, W.; Rohn, J. Experimental investigation of the hydraulic and heat-transfer properties of artificially fractured granite. Sci. Rep. 2017, 7, 39882. [Google Scholar] [CrossRef] [Green Version]
  22. He, R.; Rong, G.; Tan, J.; Cheng, L. Numerical investigation of fracture morphology effect on heat transfer characteristics of water flow through a single fracture. Geothermics 2019, 82, 51–62. [Google Scholar] [CrossRef]
  23. Shu, B.; Zhu, R.; Elsworth, D.; Dick, J.; Liu, S.; Tan, J.; Zhang, S. Effect of temperature and confining pressure on the evolution of hydraulic and heat transfer properties of geothermal fracture in granite. Appl. Energy 2020, 272, 115290. [Google Scholar] [CrossRef]
  24. Zhang, X.; Wang, Z.; Sun, Y.; Zhu, C.; Xiong, F.; Tang, Q. Numerical Simulation on Heat Transfer Characteristics of Water Flowing through the Fracture of High-Temperature Rock. Geofluids 2020, 2020, 1–14. [Google Scholar] [CrossRef]
  25. Andrade, J.S., Jr.; Henrique, E.A.A.; Almeida, M.P.; Costa, M.H.A.S. Heat transport through rough channels. Physica. Sect. A 2004, 339, 296–310. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, L.; Jiang, P.; Wang, Z.; Xu, R. Convective heat transfer of supercritical CO2 in a rock fracture for enhanced geothermal systems. Appl. Therm. Eng. 2017, 115, 923–936. [Google Scholar] [CrossRef]
  27. Fox, D.B.; Koch, D.L.; Tester, J.W. The effect of spatial aperture variations on the thermal performance of discretely fractured geothermal reservoirs. Geotherm. Energy 2015, 3, 21. [Google Scholar] [CrossRef] [Green Version]
  28. Luo, S.; Zhao, Z.; Peng, H.; Pu, H. The role of fracture surface roughness in macroscopic fluid flow and heat transfer in fractured rocks. Int. J. Rock Mech. Min. Sci. 2016, 87, 29–38. [Google Scholar] [CrossRef]
  29. Li, Z.-W.; Feng, X.-T.; Zhang, Y.-J.; Xu, T.-F. Feasibility study of developing a geothermal heating system in naturally fractured formations: Reservoir hydraulic properties determination and heat production forecast. Geothermics 2018, 73, 1–15. [Google Scholar] [CrossRef]
  30. Elmouttie, M.K.; Poropat, G.V. A Method to Estimate In Situ Block Size Distribution. Rock Mech. Rock Eng. 2011, 45, 401–407. [Google Scholar] [CrossRef]
  31. Leung, C.T.O.; Zimmerman, R.W. Estimating the Hydraulic Conductivity of Two-Dimensional Fracture Networks Using Network Geometric Properties. Transp. Porous Media 2012, 93, 777–797. [Google Scholar] [CrossRef]
  32. Bruel, D. Heat extraction modeling from forced fluid-flow through stimulated fractured rock masses - application to the rosemanowes hot dry rock reservoir. Geothermics 1995, 24, 361–374. [Google Scholar] [CrossRef]
  33. Chen, Y.; Ma, G.; Wang, H. Heat extraction mechanism in a geothermal reservoir with rough-walled fracture networks. Int. J. Heat Mass Transf. 2018, 126, 1083–1093. [Google Scholar] [CrossRef]
  34. Li, J.; Sun, Z.; Zhang, Y.; Jiang, C.; Cherubini, C.; Scheuermann, A.; Torres, S.A.G.; Li, L. Investigations of heat extraction for water and CO2 flow based on the rough-walled discrete fracture network. Energy 2019, 189. [Google Scholar] [CrossRef]
  35. Li, Z.-W.; Feng, X.-T.; Zhang, Y.-J.; Zhang, C.; Xu, T.-F.; Wang, Y.-S. Experimental research on the convection heat transfer characteristics of distilled water in manmade smooth and rough rock fractures. Energy 2017, 133, 206–218. [Google Scholar] [CrossRef]
  36. Huang, Y.; Zhang, Y.; Yu, Z.; Ma, Y.; Zhang, C. Experimental investigation of seepage and heat transfer in rough fractures for enhanced geothermal systems. Renew. Energy 2019, 135, 846–855. [Google Scholar] [CrossRef]
  37. He, Y.; Bai, B.; Cui, Y.; Lei, H.; Li, X. 3D Numerical Modeling of Water–Rock Coupling Heat Transfer Within a Single Fracture. Int. J. Thermophys. 2020, 41, 1–22. [Google Scholar] [CrossRef]
  38. Shu, B.; Zhu, R.; Tan, J.; Zhang, S.; Liang, M. Evolution of permeability in a single granite fracture at high temperature. Fuel 2019, 242, 12–22. [Google Scholar] [CrossRef]
  39. Bai, B.; He, Y.; Li, X. Numerical study on the heat transfer characteristics between supercritical carbon dioxide and granite fracture wall. Geothermics 2018, 75, 40–47. [Google Scholar] [CrossRef]
  40. Fluent Inc. FLUENT 6.3 User’s Guide; Fluent Inc.: Lebanon, NH, USA, 2007. [Google Scholar]
  41. Qu, Z.-Q.; Zhang, W.; Guo, T.-K. Influence of different fracture morphology on heat mining performance of enhanced geothermal systems based on COMSOL. Int. J. Hydrogen Energy 2017, 42, 18263–18278. [Google Scholar] [CrossRef]
  42. Ma, Y.; Zhang, Y.; Yu, Z.; Huang, Y.; Zhang, C. Heat transfer by water flowing through rough fractures and distribution of local heat transfer coefficient along the flow direction. Int. J. Heat Mass Transf. 2018, 119, 139–147. [Google Scholar] [CrossRef]
  43. Chen, Y.; Zhao, Z. Heat transfer in a 3D rough rock fracture with heterogeneous apertures. Int. J. Rock Mech. Min. Sci. 2020, 134. [Google Scholar] [CrossRef]
  44. Tan, J.; Rong, G.; He, R.; Yang, J.; Peng, J. Numerical investigation of heat transfer effect on flow behavior in a single fracture. Arab. J. Geosci. 2020, 13, 1–16. [Google Scholar] [CrossRef]
  45. He, Y.; Bai, B.; Hu, S.; Li, X. Effects of surface roughness on the heat transfer characteristics of water flow through a single granite fracture. Comput. Geotech. 2016, 80, 312–321. [Google Scholar] [CrossRef]
  46. Bai, B.; He, Y.; Hu, S.; Li, X. An Analytical Method for Determining the Convection Heat Transfer Coefficient Between Flowing Fluid and Rock Fracture Walls. Rock Mech. Rock Eng. 2017, 50, 1787–1799. [Google Scholar] [CrossRef]
Figure 1. Fracture surface and surface flow direction distribution. (a) 2D random fracture profile; (b) pseudo-three-dimensional fractured rough surface stretched towards the page direction. (c) 3D view of pseudo-three-dimensional fracture rough surface.
Figure 1. Fracture surface and surface flow direction distribution. (a) 2D random fracture profile; (b) pseudo-three-dimensional fractured rough surface stretched towards the page direction. (c) 3D view of pseudo-three-dimensional fracture rough surface.
Applsci 11 00751 g001
Figure 2. Fracture seepage heat transfer model.
Figure 2. Fracture seepage heat transfer model.
Applsci 11 00751 g002
Figure 3. Mesh features of mesh 4. (a) Fracture surface mesh characteristics; (b) overall mesh distribution; (c) mesh distribution characteristics around the fracture.
Figure 3. Mesh features of mesh 4. (a) Fracture surface mesh characteristics; (b) overall mesh distribution; (c) mesh distribution characteristics around the fracture.
Applsci 11 00751 g003
Figure 4. Variation of outlet temperature with time for the four mesh cases.
Figure 4. Variation of outlet temperature with time for the four mesh cases.
Applsci 11 00751 g004
Figure 5. The fracture heat transfer laboratory simulation system. (a) Core holder system; (b) a granite sample with a smooth-plate single fracture in the holder.
Figure 5. The fracture heat transfer laboratory simulation system. (a) Core holder system; (b) a granite sample with a smooth-plate single fracture in the holder.
Applsci 11 00751 g005
Figure 6. Smooth-plate fracture and mesh model. (a) Temperature sensor arrangement. (b) Flow direction. (c) Mesh generation characteristics of numerical model.
Figure 6. Smooth-plate fracture and mesh model. (a) Temperature sensor arrangement. (b) Flow direction. (c) Mesh generation characteristics of numerical model.
Applsci 11 00751 g006
Figure 7. Variations in thermal conductivity and specific heat capacity of rock with temperature.
Figure 7. Variations in thermal conductivity and specific heat capacity of rock with temperature.
Applsci 11 00751 g007
Figure 8. Comparison of experimental and numerical results (t = 10 min). (Left) Outlet temperature varies with injection flow rate at different initial rock temperatures. (Right) Surface temperature varies with sensor position at different injection velocities.
Figure 8. Comparison of experimental and numerical results (t = 10 min). (Left) Outlet temperature varies with injection flow rate at different initial rock temperatures. (Right) Surface temperature varies with sensor position at different injection velocities.
Applsci 11 00751 g008
Figure 9. Temperature field distribution on the fracture surface (t = 10 min).
Figure 9. Temperature field distribution on the fracture surface (t = 10 min).
Applsci 11 00751 g009
Figure 10. Local temperature distribution along the flow direction at different locations of the x-axis on the inner wall of fracture (solid line: x = −0.02 m, dashed line: x = 0.02 m, t = 2 min).
Figure 10. Local temperature distribution along the flow direction at different locations of the x-axis on the inner wall of fracture (solid line: x = −0.02 m, dashed line: x = 0.02 m, t = 2 min).
Applsci 11 00751 g010
Figure 11. Streamline distribution characteristics of models with α = 30° and α = 60° (V0 = 10 mL/min).
Figure 11. Streamline distribution characteristics of models with α = 30° and α = 60° (V0 = 10 mL/min).
Applsci 11 00751 g011
Figure 12. Schematic illustration of fluid flow direction deviation. (a) The up- and down-hill direction of the fluid migration. (b) Top view: solid and dotted lines are the ideal flow direction path and the flow deflection path, respectively.
Figure 12. Schematic illustration of fluid flow direction deviation. (a) The up- and down-hill direction of the fluid migration. (b) Top view: solid and dotted lines are the ideal flow direction path and the flow deflection path, respectively.
Applsci 11 00751 g012
Figure 13. Reynolds number in relation to (a) the flow rate and (b) Nusselt number.
Figure 13. Reynolds number in relation to (a) the flow rate and (b) Nusselt number.
Applsci 11 00751 g013
Figure 14. Histogram of cumulative heat extraction rate determined for 30 min.
Figure 14. Histogram of cumulative heat extraction rate determined for 30 min.
Applsci 11 00751 g014
Figure 15. Variation of heat extraction rate with time at several initial flow rates and initial temperatures.
Figure 15. Variation of heat extraction rate with time at several initial flow rates and initial temperatures.
Applsci 11 00751 g015
Figure 16. Variations in the total heat transfer coefficient with time for various initial conditions. Top: T0 = 90 °C. Below: V0 = 10 mL/min.
Figure 16. Variations in the total heat transfer coefficient with time for various initial conditions. Top: T0 = 90 °C. Below: V0 = 10 mL/min.
Applsci 11 00751 g016
Figure 17. Distribution of local heat transfer coefficients on rock fractures as water flows through the fracture (V0 = 10 mL/min. T0 = 90 °C. t = 10 min).
Figure 17. Distribution of local heat transfer coefficients on rock fractures as water flows through the fracture (V0 = 10 mL/min. T0 = 90 °C. t = 10 min).
Applsci 11 00751 g017
Figure 18. Distribution of local heat transfer coefficients along the outlet direction at different flow rates.
Figure 18. Distribution of local heat transfer coefficients along the outlet direction at different flow rates.
Applsci 11 00751 g018aApplsci 11 00751 g018b
Figure 19. Variation of hx with the flow rate at adjacent peak and troughs of the fracture profile.
Figure 19. Variation of hx with the flow rate at adjacent peak and troughs of the fracture profile.
Applsci 11 00751 g019
Figure 20. Local heat transfer coefficient, wall temperature, and average velocity distributed along the x-direction at z = 0.05 m (V0 = 20 mL/min, T0 = 90 °C).
Figure 20. Local heat transfer coefficient, wall temperature, and average velocity distributed along the x-direction at z = 0.05 m (V0 = 20 mL/min, T0 = 90 °C).
Applsci 11 00751 g020aApplsci 11 00751 g020b
Table 1. Geometric parameters of the fracture.
Table 1. Geometric parameters of the fracture.
Model Case (α =)30°60°90°Maximum Deviation
Surface area (mm2)5588.955620.335658.535553.841.87%
Inlet area (mm2)14.988314.939114.940514.94250.33%
Volume (mm3)1496.281495.961494.911494.250.14%
Vertical offset (mm)0.30.30.30.3\
Table 2. Four grid numbers for α = 0°.
Table 2. Four grid numbers for α = 0°.
Mesh Case1234
Mesh Number691,2571,366,5542,628,8253,120,142
Table 3. Number of simulated sample grids for the four models.
Table 3. Number of simulated sample grids for the four models.
Model Case (α =)30°60°90°
Mesh Number2,628,8252,586,7432,645,4262,653,564
Table 4. Numerical simulation conditions.
Table 4. Numerical simulation conditions.
Model Case (α =)Rock Temperature, T0 (°C)Flow Rate, V0 (mL/min)
0° 30° 60° 90°70 9010 20 30 40
Table 5. Results from 10-min numerical simulation.
Table 5. Results from 10-min numerical simulation.
T0 (°C)V0 (mL/min)Toutlet (°C)
α = 0°α = 30°α = 60°α = 90°
70°C1054.0853.9353.8754.19
2036.4736.3636.2336.51
3030.9230.8730.8030.99
4028.6728.6528.5928.73
90°C1068.4568.3868.2368.58
2042.0341.9641.8742.18
3033.7533.7033.6233.84
4030.4430.3430.2730.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, X.; Zhang, Y.; Hu, Z.; Huang, Y. Numerical Investigation on the Influence of Surface Flow Direction on the Heat Transfer Characteristics in a Granite Single Fracture. Appl. Sci. 2021, 11, 751. https://doi.org/10.3390/app11020751

AMA Style

Gao X, Zhang Y, Hu Z, Huang Y. Numerical Investigation on the Influence of Surface Flow Direction on the Heat Transfer Characteristics in a Granite Single Fracture. Applied Sciences. 2021; 11(2):751. https://doi.org/10.3390/app11020751

Chicago/Turabian Style

Gao, Xuefeng, Yanjun Zhang, Zhongjun Hu, and Yibin Huang. 2021. "Numerical Investigation on the Influence of Surface Flow Direction on the Heat Transfer Characteristics in a Granite Single Fracture" Applied Sciences 11, no. 2: 751. https://doi.org/10.3390/app11020751

APA Style

Gao, X., Zhang, Y., Hu, Z., & Huang, Y. (2021). Numerical Investigation on the Influence of Surface Flow Direction on the Heat Transfer Characteristics in a Granite Single Fracture. Applied Sciences, 11(2), 751. https://doi.org/10.3390/app11020751

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop