Abstract：A fast and highly efficient numerical computational method has important practical significance for development and application of the technology in groundwater source heat pumps. Normally the boundary of an aquifer has an irregularshape. Compared with unstructured grids, adopting structured grids can effectively decrease the complexity in discretization of the mathematical equations and is very beneficial for developing an adaptive mesh refinement (AMR) algorithm to simulate the seepage flow field of groundwater source heat pump system. To develop a computational method based on the structured grids, a simple pretreatment method was proposed to deal with the grid cells on the aquifer boundary. In the method, the effective permeability on the interfaces between grid cells rather than those of the grid cells was involved in computations, further an algorithm for calculating the effective permeability on the six faces of a grid cell was proposed. The threedimensional computational examples show that even more coarse cells are adopted a more accurate result can be reached for a domain with irregular boundary. In the AMR algorithm the cell sizes depend on the temperature gradient, i.e. the smaller cell size the lower the gradient; as a result of this, the time consuming for simulating the underground seepage flow in a groundwater source heat pump system is reduced considerably but the accuracy remains the same.
［1］Sanner B, Karytas C, Mendrinos D, et al. Current status of ground source heat pumps and underground thermal energy storage in Europe［J］.Geothermics, 2003, 32(4/5/6): 579-588.［2］Abdeen M O. Groundsource heat pumps systems and applications［J］.Renewable and Sustainable Energy Reviews, 2008, 12(2): 344-371.［3］王秉忱. 我国利用水源、地源热泵技术开发浅层地热能呈现好形势［J］.工程建设与设计, 2008 （1）: 8-11.Wang Bingchen. The excellent situation of exploitation of shallow geothermal energy by utilizing the WSHP and GSHP technology in China［J］.Construction ＆ Design for Project, 2008 （1）: 8-11. (in Chinese)［4］周光辉, 张震.蓄能型地下水源热泵在我国推广应用的可行性分析［J］.中原工学院学报, 2005, 16(6): 26-28.Zhou Guanghui, Zhang Zhen.The possibility analysis of the storage energy ground water heat pump expanding and applying in China［J］.Journal of Zhongyuan Institute of Technology, 2005, 16(6): 26-28.(in Chinese)［5］Molz F J, Warman J C, Jones T E. Aquifer storage of heated water—Part I: Field experiment［J］. Ground Water, 1978, 16(6): 234-241.［6］Papadopulos S S, Larson S P. Aquifer storage of heated water—Part II: Numerical simulation of field results［J］. Ground Water, 1978, 16(6): 242-248.［7］张远东, 魏加华, 李宇, 等. 地下水源热泵采能的水-热耦合数值模拟［J］. 天津大学学报:自然科学版, 2006, 39(8): 907-912.Zhang Yuandong, Wei Jiahua, Li Yu, et al. Simulation of changes in Geotemperature field due to energy abstration from underground aquifers［J］.Journal of Tianjin University, 2006, 39(8): 907-912.(in Chinese)［8］Nagano K, Mochida T, Ochifuji K.Influence of natural convection on forced horizontal flow in saturated porous media for aquifer thermal energy storage［J］.Applied Thermal Engineering, 2002, 22(12): 1299-1311.［9］倪龙, 马最良. 同井回灌地下水源热泵地下水渗流理论研究［J］. 太阳能学报, 2006, 27(12): 1219-1224.Ni Long, Ma Zuiliang. Study on the seepage theory for groundwater heat pump with pumping ＆ recharge in the same well［J］. Acta Energiae Solaris Sinica, 2006, 27(12): 1219-1224.(in Chinese)［10］Paksoy H O, Andersson O, Abacic A, et al. Heating and cooling of a hospital using solar energy coupled with seasonal thermal energy storage in an aquifer［J］. Renewable Energy, 2000, 19(1/2): 117-122.［11］Berger M J, Oliger J. Adaptive mesh refinement for hyperbolic partial differential equations［J］. J Comput Phys，1984, 53(3): 484-512.［12］Hornung R D, Trangenstein J. Adaptive mesh refinement and multilevel iteration for flow in porous media［J］. J Comput Phys, 1997, 136(2): 522-545.［13］Trangenstein J. Multiscale iterative techniques and adaptive mesh refinement for flow in porous media［J］. Advances in Water Resources, 2002, 25(8/9/10/11/12): 1175-1213.［14］Wang X H, Quintard M, Darche G. Adaptive mesh refinement for onedimensional threephase flow with phase change in porous media［J］. Numerical Heat Transfer: Part B, 2006, 50(3): 231-268.［15］施安峰,王晓宏,罗海山. 三维蒸汽辅助重力驱油过程的自适应网格法［J］. 辽宁工程技术大学学报:自然科学版,2009, 28(S1): 4-7.Shi Anfeng, Wang Xiaohong, Luo Haishan. Adaptive mesh refinement for threedimensional SAGD process［J］. Journal of Liaoning Technical University:Natural Science, 2009, 28(S1): 4-7.(in Chinese)
