Research Article | Open Access

Haijun Liang, Xiaofeng Guo, Tao Gao, Lingbao Wang, Xianbiao Bu, "Performance Analysis of Single-Well Enhanced Geothermal System for Building Heating", *The Scientific World Journal*, vol. 2020, Article ID 7927052, 8 pages, 2020. https://doi.org/10.1155/2020/7927052

# Performance Analysis of Single-Well Enhanced Geothermal System for Building Heating

**Academic Editor:**Christophe Waterlot

#### Abstract

Deep borehole heat exchanger (DBHE) technology does not depend on the existence of hot water reservoir and can be used in various regions. However, the heat extraction from DBHE can hardly be improved due to poor thermal conductivity of rocks. Here, a single-well enhanced geothermal system (SWEGS) is proposed, which has a larger heat-exchange area of artificial reservoir created by fracturing hydrothermal technology. We find that, due to heat convection between rocks and fluid, the extracted thermal output for SWEGS is 4772.73 kW, which is 10.64 times of that of DBHE. By changing the injection water temperature, volume flow rate, and artificial reservoir volume, it is easy to adjust the extracted thermal output to meet the requirement of building thermal loads varying with outdoor air temperature. Understanding these will enable us to better apply SWEGS technology and solve the fog and haze problem easily and efficiently.

#### 1. Introduction

In recent years, one of the greatest challenges in China is the fog and haze problem. About 83% of space heating areas in north China use coal, and annual coal consumptions are 400 million ton-coal-equivalent where decentralized coals in rural areas account for half of total consumptions [1]. A coal-based heating system contributes significantly to air pollution in winter. There is thus an urgent need for clean energy building heating technology at the present stage. Geothermal energy has showed significant potential as renewable energy resource because of its low environmental impact, low greenhouse gas emissions, and technical feasibility [2, 3].

The hydrothermal system for geothermal heating relies heavily on the existence of rich hot water reservoir [3]. Moreover, the recharging problem also limits its application and promotion. Without the resource constraint and recharging problem, the technology of deep borehole heat exchanger (DBHE), as shown in Figure 1, is widely accepted in China, and it can be applied almost everywhere. Some researchers have focused on the utilization of geothermal energy acquired from DBHE for building heating [4–16]. Bu et al. carried out experimental and simulation studies of DBHE for building heating, and the research results are shown in Figure 2 [17, 18]. In Figure 2, DBHE is configured in a geothermal well with a depth of 2605 m, diameter 0.1778 m, bottomhole temperature 87.42°C, and average geothermal gradient 27.8°C/km. The volume flow rate is about 30 m^{3}/h throughout the experimental process. It is obvious from Figure 2(a) that *Q* always decreases with time. The average extracted thermal output in the first heating season is 448.49 kW. As shown in Figure 2(b), the rock temperature decreases due to being cooled by the injection water, particularly the rocks close to well wall having the maximum temperature drop. The reason for this is that the heat loss from rocks to fluid cannot rapidly be compensated by heat conduction of rocks due to its poor thermal conductivity, resulting in a lower thermal output. It can thus be concluded that one of the effective means of improving the heat extraction from DBHE is to enhance the heat transfer of rocks, especially in near well bore zone. According to the principle of heat transfer, the heat transfer coefficient of the convection is far greater than that of the heat conduction. To conduct the convective heat exchange between fluid and rocks, the enhanced geothermal system (EGS) with artificial fractures created by the hydraulic fracturing was proposed [19, 20]. The conventional EGS has two or multiple vertical or directional wells, which causes the high drilling cost and thus limits its extensive applications [19, 20]. Moreover, the creation of the artificial fractures has a high risk, especially in the hard rock [21, 22]. In addition, the various physical and chemical processes take place in the geothermal reservoir during the process of fluid injection and extraction [23, 24].

**(a)**

**(b)**

**(c)**

To circumvent these problems, a single-well enhanced geothermal system (SWEGS), as illustrated in Figures 3(a) and 3(b), is proposed in this study. There are many hydrothermal systems in urban areas at the depth of less than 3000 m, which usually have small water output (usually less than 30 m^{3}/h or even less). If the shallow depth hydrothermal system can be transformed into SWEGS for building heating, this not only has a good use of natural structural fissures in hydrothermal system but also reduces the risk and investment costs of EGS, and thus broadening the application range of EGS and effectively solving the fog and haze problem caused by winter heating using fossil energy. The artificial reservoir for SWEGS is created by two steps (1) drilling one main well and multilateral wells, as shown in Figure 3(a); (2) implementing hydraulic fracturing in the multilateral wells and transporting proppants to produce artificial reservoir. The function of the multilateral wells is to disperse fluid uniformly in the artificial reservoir and to make easy to fracture rocks. The common problems in both SWEGS and conventional EGS are, respectively, (1) the loss of circulating fluid; (2) the detection of formation stress distribution, fracture orientations, and tracer-swept volume; and (3) the risk of seismicity induced by hydraulic fracturing and so on. Compared to conventional EGS implemented in the hard rock, the creation of fractures is easy in hydrothermal system due to having many natural fractures. The achievements from SWEGS are expected to apply to conventional EGS and thus reducing its risk and costs. Although the spatial-temporal variation of fracture aperture during injection/production significantly influences the heat extraction from the reservoir, this is not the focus of this study [25]. The main aim of this study is to evaluate the performance of SWEGS so that more people know its performance and apply it for buildings heating to solve the fog and haze problem caused by using fossil energy.

**(a)**

**(b)**

#### 2. Methods

##### 2.1. Physical Model

The physical model of main well and multilateral wells and SWEGS is shown in Figures 3(a) and 3(b). In Figure 3, the main well has a depth of 2605 m, diameter 0.1778 m, bottomhole temperature 87.42°C, and average geothermal gradient 27.8°C/km. The density, specific heat, and thermal conductivity for rocks are, respectively, 2800 kg/m^{3}, 0.92 kJ/kg/K, and 3.49 W/m/K. The rock physical property, well depth, and geothermal gradient for SWEGS are the same as those of DBHE. The diameter and length (distance from the main well) of multilateral wells are 0.08 m and 50 m, respectively. The artificial reservoir located at a depth of 2100–2600 m has a thickness of 500 m and a radius of 50 m. In general, there are more than four multilateral wells in the engineering projects.

Facing the engineering application, the following questions related to SWEGS technology need to be answered: (1) is the extracted thermal output of SWEGS stable and sustainable? (2) Can the extracted thermal output be easily adjusted as the variation of the building heating loads occurs? (3) What is the appropriate distance (well spacing) between two wells and appropriate thickness for artificial reservoir? These questions are hard to be addressed by the experimental studies due to the high cost. Alternatively, as a conventional method, numerical simulation can be used to analyze SWEGS performance in order to answer the above questions.

##### 2.2. Model Equations

Assuming that the rocks enclosing the artificial reservoir are impermeable, that is, there is no fluid loss in the SWEGS. In addition, the influence of fluid-rock reaction is neglected during the process of injection and extraction fluid. The artificial reservoir is treated as equivalent to a porous medium with uniform porosity *ε* and permeability *K* [26–28]. Local thermal nonequilibrium between the porous rock matrix and fluid is considered, and thus, two energy equations are employed. The governing equations describing the conservation of continuity, momentum, and energy are formulated as the following [29, 30].

Fluid continuity equation is

Fluid momentum equation is

Energy equation in the porous rock matrix is

Energy equation for heat transport in the fluid iswhere , , , and are the fluid velocity vector, fluid pressure, rock matrix temperature, and fluid temperature, respectively; , , and with subscript “” and “”, respectively, indicate the density, specific heat capacity, and heat conductivity of fluid and rock matrix, respectively. The symbol is the viscosity of fluid. As annotated by the superscript, , the heat conductivity of fluid and rock matrix in the artificial reservoir is implemented with the effective forms, i.e., and , where denotes the artificial reservoir porosity [26]. The symbol denotes the volumetric heat exchange coefficient between the porous rock matrix and fluid.

##### 2.3. Initial and Boundary Conditions

The artificial reservoir is treated as equivalent to a porous medium, and its porosity and permeability is, respectively, 0.1 and 10^{−12} m^{2}. The injection water temperature and volume flow rate are, respectively, 5°C and 50 m^{3}/h. The heating time per year is 140 days, and the rest time is used for heat recovery.

It is assumed that the temperature of rocks is constant and not influenced by the fluid when *r* is larger than 200 m.

The heat exchange between the interface of artificial reservoir and surrounding rocks is given as follows:where is the convective heat transfer coefficient, which is calculated according to Dittus–Boelter formula [31].where *d*_{e} is the hydraulic diameter, *d*_{e} = 100 m.

##### 2.4. Numerical Method

Equations (1)–(4) together with the initial and boundary conditions are solved using MATLAB software. The algorithm of SIMPLE (Semi-Implicit Method for Pressure Linked Equation) is used to address the pressure-velocity coupling [32, 33].

#### 3. Results and Discussion

In Figure 4, the annual mean extracted thermal output for SWEGS ranges from 4431.38 to 4772.73 kW during ten heating seasons, with percentage reduction about 0.715% per year. The average *Q* for ten years is 4591.43 kW, which can provide heat for 131183.71 m^{2} buildings with a specific heat load of 35 W/m^{2} and can reduce carbon dioxide emissions by 13884.48 tons per heating season. The heating price is 4.32 USD/m^{2} each heating season at Qingdao in China, and thus, the total income from building heating is 566713.65 USD/year. In the first heating season, the annual mean extracted thermal output for SWEGS is 4772.73 kW, which is 10.64 times of that of DBHE. The extracted thermal output from SWEGS is much greater than that of DBHE, which is mainly caused by two reasons. First, the convective heat transfer is formed between fluid and rocks in the artificial reservoir, whose coefficient of heat transfer is much higher than that of heat conduction in DBHE. Second, large heat transfer area is created in the artificial reservoir by implementing hydraulic fracturing. The extracted heat is from two parts: one is heat convection between rocks and fluid in artificial reservoir, and the other is heat conduction from surrounding rocks to artificial reservoir. In Figure 2(c), the annual mean *Q* decreases drastically during the first few years, and then, the decreasing rate slows down gradually. The average *Q* for ten years in DBHE is 424.45 kW. The annual mean *Q* for DBHE is, respectively, 448.49 and 413.63 kW in the first and tenth heating seasons, which has a percentage change about 7.77%.

Figure 5(a) shows that there is a lower temperature gradient for the rocks in near well bore zone for SWEGS due to having a higher coefficient of convective heat transfer between fluid and rocks. Compared with undisturbed (initial) rocks temperature, the maximum rocks temperature drop for SWEGS is 8.63°C at a depth of 2500 m for ten years’ operation, while it is 40.16°C for DBHE.

**(a)**

**(b)**

Figure 5(b) implies that more heat at the upper artificial reservoir is extracted by circulating water, as the injection water temperature is only 5°C. The temperature drop of artificial reservoir at a depth of 2400 m is big, while it becomes very small at a depth of 2500 m, which indicates that for ten years’ operation, the cold front surface has not yet reached the depth of 2500 m, that is, to say, the service life for SWEGS is more than 10 years in terms of the reservoir thickness of 500 m, a multilateral well length of 50 m, the volume flow rate of 50 m^{3}/h, and an injection water temperature of 5°C.

From Figure 5(b), the rocks temperature keeps unchanged (undisturbed) for *r* > 100 m at the end of tenth heating season, which indicates that the thermal influence radius (impact scope) is about 100 m for ten heating seasons. This leads us to conclude that in terms of the above parameters of the artificial reservoir shape and size and operation condition, the well spacing in the practical projects is recommended to be no less than 200 m with 10 years’ operation period in order to avoid the mutual interference.

From view point of heat transfer principle, three main methods can be used to enhance heat transfer between rocks and fluid. The first method is to improve heat transfer coefficient, which is reflected by improving the volume flow rate in Figure 6; the second method is to increase the heat transfer temperature difference, which is reflected by changing *T*_{in} in Figure 6; the third method is to increase heat exchange area, which is reflected by increasing the reservoir volume in Figure 7.

In Figures 6 and 7, the symbols V and *L*, respectively, represent the volume flow rate and the length of multilateral well. In Figure 6, the abscissa is the heating time at the first heating season (3360 h or 140 days). As evident in Figure 6, a higher volume flow rate will result in an increase in *Q*. A higher V will cause a higher convective heat transfer coefficient between rocks and fluid, thus leading to a greater *Q*. Generally speaking, the pump power has almost a cubic relationship with flow velocity. Consequently, an increase in V will require a rapidly increasing level of pump power. In the actual project, the pump power and *Q* should, therefore, be comprehensively considered so as to decide a reasonable V.

In addition, the temperature difference between rock and fluid increases with the decrease of *T*_{in} and thus leading to an increase in *Q*, which shows that the lower the *T*_{in}, the higher the *Q*. Note that *T*_{in} should not be lower than zero temperature using water as working fluid; otherwise, it is easy to cause water to freeze. Therefore, the antifreeze liquid with good heat transfer performance and moderate cost can be considered as the injection fluid.

From Figure 6, *Q* can be adjusted by changing the injection water temperature and volume flow rate. In general, the outdoor air temperature has a great influence on the building heating loads. By adjusting the injection water temperature and volume flow rate, SWEGS can meet the requirement of building heating loads influenced by outdoor air temperature. Furthermore, the imbalance of the extracted thermal output among different heating seasons can also be adjusted by changing the injection water temperature and volume flow rate.

In Figure 7, the volume flow rate is 60 m^{3}/h for artificial reservoir with *L* = 150 m, and other parameters are the same as those of the reservoir with *L* = 50 m. Figure 7 indicates that the larger the reservoir volume is, the slower the attenuation degree of *Q*. However, a large reservoir volume will result in a high engineering cost. Based on the above discussion, it is thus clear that the system performance, reservoir volume, and engineering cost should be comprehensively considered so as to decide their optimum matching relation.

In Figure 2(a), the temperature of extracted water from DBHE is lower than 20°C, which does not meet the demand of the heating temperature. Therefore, heat pump with user sides supply and return water temperature of 45°C and 40°C is needed in order to provide heat for building. While for SWEGS, the extracted water temperature is higher than 80°C, which can be used efficiently by adopting the model of cascade utilization (*T*_{out} > 60°C, direct heating by radiators; 45 < *T*_{out} < 60°C, heating by the fan coil units; 30 < *T*_{out} < 45°C, floor radiation heating; and *T*_{out} < 30°C, heating by heat pump).

#### 4. Conclusions

Deep borehole heat exchanger (DBHE) technology does not depend on the existence of hot water reservoir and can be used in various regions. However, the heat extraction from DBHE can hardly be improved due to poor thermal conductivity of rocks. To conduct the convective heat exchange between fluid and rocks, the enhanced geothermal system (EGS) with artificial fractures created by the hydraulic fracturing was proposed. However, the EGS technology has a high drilling cost and a high risk for artificial fracturing. To circumvent these problems, a single-well enhanced geothermal system (SWEGS) is proposed.

There are many hydrothermal systems in urban areas at the depth of less than 3 km. In order to reduce the risk and the costs of EGS, the shallow depth hydrothermal system should be considered first. The shallow depth hydrothermal system has natural fractures, and it is easy to conduct hydraulic fracturing and stimulation to create new fractures and highly conductive zones and thus increasing the volume flow rate of SWEGS and effectively improving the heat-exchanging amount. Due to heat convection between rocks and fluid, the extracted thermal output for SWEGS is 4772.73 kW, which is 10.64 times of that of DBHE. Besides, by changing the injection water temperature, volume flow rate, and artificial reservoir volume, it is easy to adjust the extracted thermal output to meet the requirement of building thermal loads changing with outdoor air temperature. As a result, the fog and haze problem in China caused by winter heating using fossil energy can be effectively solved by applied SWEGS due to having many shallow depth hydrothermal systems in urban areas. In addition, the researchers should make great efforts to study the possibility or risk of seismicity induced by hydraulic fracturing during the process of creating artificial reservoir.

#### Data Availability

The data used to support the findings of this study are included within the supplementary information file.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

Xiaofeng Guo and Tao Gao involved in conceptualization; Haijun Liang involved in investigation; Xianbiao Bu involved in writing the original draft; Lingbao Wang involved in writing the review and editing. All authors have read and agreed to the published version of the manuscript.

#### Acknowledgments

The authors thank D. D. Zhang, H. S. Li, L. B. Wang, Y. L. Gong, and W. B. Ma for their contributions to this work. This work was supported by the National Natural Science Foundation of China (No. 41972314).

#### Supplementary Materials

The raw data used to plot Figures 2, 4, 5, 6, and 7 are provided in the supplementary file. The symbol Q represents the extracted thermal output from DBHE, *T*_{in} is the injection water temperature, *T*_{out} is the extracted water temperature, and r is the radial distance from the well tube. The symbols V and L, respectively, denote the volume flow rate and the length of multilateral well. .* (Supplementary Materials)*

#### References

- C. Rubio-Maya, V. M. Ambríz Díaz, E. Pastor Martínez, and J. M. Belman-Flores, “Cascade utilization of low and medium enthalpy geothermal resources: A review,”
*Renewable and Sustainable Energy Reviews*, vol. 52, pp. 689–716, 2015. View at: Publisher Site | Google Scholar - A. Keçebas, “Exergoenvironmental analysis for a geothermal district heating system: an application,”
*Energy*, vol. 94, pp. 391–400, 2016. View at: Publisher Site | Google Scholar - J. W. Lund and T. L. Boyd, “Direct utilization of geothermal energy 2015 worldwide review,” in
*Proceedings World Geothermal Congress 2015*, pp. 19–25, Melbourne, Australia, April 2015. View at: Google Scholar - Y. Noorollahi, M. Pourarshad, S. Jalilinasrabady, and H. Yousefi, “Numerical simulation of power production from abandoned oil wells in Ahwaz oil field in southern Iran,”
*Geothermics*, vol. 55, pp. 16–23, 2015. View at: Publisher Site | Google Scholar - N. M. Wight and N. S. Bennett, “Geothermal energy from abandoned oil and gas wells using water in combination with a closed wellbore,”
*Applied Thermal Engineering*, vol. 89, pp. 908–915, 2015. View at: Publisher Site | Google Scholar - R. A. Caulk and I. Tomac, “Reuse of abandoned oil and gas wells for geothermal energy production,”
*Renewable Energy*, vol. 112, pp. 388–397, 2017. View at: Publisher Site | Google Scholar - G. Cui, S. Ren, L. Zhang et al., “Geothermal exploitation from hot dry rocks via recycling heat transmission fluid in a horizontal well,”
*Energy*, vol. 128, pp. 366–377, 2017. View at: Publisher Site | Google Scholar - C. Alimonti and E. Soldo, “Study of geothermal power generation from a very deep oil well with a wellbore heat exchanger,”
*Renewable Energy*, vol. 86, pp. 292–301, 2016. View at: Publisher Site | Google Scholar - S. Gharibi, E. Mortezazadeh, S. J. Hashemi Aghcheh Bodi, and A. Vatani, “Feasibility study of geothermal heat extraction from abandoned oil wells using a U-tube heat exchanger,”
*Energy*, vol. 153, pp. 554–567, 2018. View at: Publisher Site | Google Scholar - W.-L. Cheng, T.-T. Li, Y.-L. Nian, and K. Xie, “Evaluation of working fluids for geothermal power generation from abandoned oil wells,”
*Applied Energy*, vol. 118, pp. 238–245, 2014. View at: Publisher Site | Google Scholar - Y. Shang, M. Dong, and S. Li, “Intermittent experimental study of a vertical ground source heat pump system,”
*Applied Energy*, vol. 136, no. 31, pp. 628–635, 2014. View at: Publisher Site | Google Scholar - Y.-L. Nian and W.-L. Cheng, “Insights into geothermal utilization of abandoned oil and gas wells,”
*Renewable and Sustainable Energy Reviews*, vol. 87, pp. 44–60, 2018. View at: Publisher Site | Google Scholar - X. Song, G. Wang, Y. Shi et al., “Numerical analysis of heat extraction performance of a deep coaxial borehole heat exchanger geothermal system,”
*Energy*, vol. 164, pp. 1298–1310, 2018. View at: Publisher Site | Google Scholar - L. Fang, N. Diao, Z. Shao, K. Zhu, and Z. Fang, “A computationally efficient numerical model for heat transfer simulation of deep borehole heat exchangers,”
*Energy and Buildings*, vol. 167, pp. 79–88, 2018. View at: Publisher Site | Google Scholar - D. Gordon, T. Bolisetti, D. S.-K. Ting, and S. Reitsma, “Short-term fluid temperature variations in either a coaxial or U-tube borehole heat exchanger,”
*Geothermics*, vol. 67, pp. 29–39, 2017. View at: Publisher Site | Google Scholar - D. Gordon, T. Bolisetti, D. S.-K. Ting, and S. Reitsma, “A physical and semi-analytical comparison between coaxial BHE designs considering various piping materials,”
*Energy*, vol. 141, pp. 1610–1621, 2017. View at: Publisher Site | Google Scholar - X. Bu, Y. Ran, and D. Zhang, “Experimental and simulation studies of geothermal single well for building heating,”
*Renewable Energy*, vol. 143, pp. 1902–1909, 2019. View at: Publisher Site | Google Scholar - X. Bu, K. Jiang, and H. Li, “Performance of geothermal single well for intermittent heating,”
*Energy*, vol. 186, Article ID 115858, 2019. View at: Publisher Site | Google Scholar - P. Olasolo, M. C. Juárez, M. P. Morales, S. D´Amico, and I. A. Liarte, “Enhanced geothermal systems (EGS): a review,”
*Renewable and Sustainable Energy Reviews*, vol. 56, pp. 133–144, 2016. View at: Publisher Site | Google Scholar - S.-M. Lu, “A global review of enhanced geothermal system (EGS),”
*Renewable and Sustainable Energy Reviews*, vol. 81, pp. 2902–2921, 2018. View at: Publisher Site | Google Scholar - M. W. McClure and R. N. Horne, “An investigation of stimulation mechanisms in Enhanced Geothermal Systems,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 72, pp. 242–260, 2014. View at: Publisher Site | Google Scholar - W. Zhang, Z. Qu, T. Guo, and Z. Wang, “Study of the enhanced geothermal system (EGS) heat mining from variably fractured hot dry rock under thermal stress,”
*Renewable Energy*, vol. 143, pp. 855–871, 2019. View at: Publisher Site | Google Scholar - S. N. Pandey and V. Vishal, “Sensitivity analysis of coupled processes and parameters on the performance of enhanced geothermal systems,”
*Scientific Reports*, vol. 7, p. 17057, 2017. View at: Publisher Site | Google Scholar - S. Salimzadeh and H. M. Nick, “A coupled model for reactive flow through deformable fractures in enhanced geothermal systems,”
*Geothermics*, vol. 81, pp. 88–100, 2019. View at: Publisher Site | Google Scholar - S. N. Pandey, V. Vishal, and A. Chaudhuri, “Geothermal reservoir modeling in a coupled thermo-hydro-mechanical-chemical approach: a review,”
*Earth-Science Reviews*, vol. 185, pp. 1157–1169, 2018. View at: Publisher Site | Google Scholar - C.-L. Wang, W.-L. Cheng, Y.-L. Nian, L. Yang, B.-B. Han, and M.-H. Liu, “Simulation of heat extraction from CO
_{2}-based enhanced geothermal systems considering CO_{2}sequestration,”*Energy*, vol. 142, pp. 157–167, 2018. View at: Publisher Site | Google Scholar - P. Asai, P. Panja, R. Velasco, J. McLennan, and J. Moore, “Fluid flow distribution in fractures for a doublet system in enhanced geothermal systems (EGS),”
*Geothermics*, vol. 75, pp. 171–179, 2018. View at: Publisher Site | Google Scholar - Q. Cheng, X. Wang, and A. Ghassemi, “Numerical simulation of reservoir stimulation with reference to the Newberry EGS,”
*Geothermics*, vol. 77, pp. 327–343, 2019. View at: Publisher Site | Google Scholar - W. Huang, W. Cao, and F. Jiang, “A novel single-well geothermal system for hot dry rock geothermal energy exploitation,”
*Energy*, vol. 162, pp. 630–644, 2018. View at: Publisher Site | Google Scholar - W. Cao, W. Huang, and F. Jiang, “A novel thermal-hydraulic-mechanical model for the enhanced geothermal system heat extraction,”
*International Journal of Heat and Mass Transfer*, vol. 100, pp. 661–671, 2016. View at: Publisher Site | Google Scholar - X. M. Zhang, Z. P. Ren, and F. M. Mei,
*Heat Transfer*, China Construction Industry Press, Beijing, China, 4th edition, 2001. - R. X. Li,
*Basis of Finite Volume Method*, National Defense Industry Press, Beijing, China, 2nd edition, 2008. - W. Q. Tao,
*Numerical Heat Transfer*, Xi’an Jiao Tong University Press, Xi’an, China, 2nd edition, 2001.

#### Copyright

Copyright © 2020 Haijun Liang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.