基于LBM-MFLS耦合模型的含水层渗流-传热过程模拟  

NUMERICAL SIMULATION OF FLUID SEEPAGE AND HEAT TRANSFER IN AQUIFER WITH LBM-MFLS COUPLED MODEL

作　　者：马玖辰 吕林海[1] 杨杰 崔阿凤 魏璠 Ma Jiuchen;Lyu Linhai;Yang Jie;Cui Afeng;Wei Fan(College of Energy Safety Engineering,Tianjin Chengjian University,Tianjin 300384,China;Key Laboratory for Efficient Use of Low and Medium Grade Energy,Ministry of Education of China,Tianjin University,Tianjin 300072,China)

机构地区：[1]天津城建大学能源与安全工程学院,天津300384 [2]天津大学中低温热能高效利用教育部重点实验室,天津300072

出　　处：《太阳能学报》2023年第5期30-39,共10页Acta Energiae Solaris Sinica

基　　金：国家自然科学基金(41402228);天津市教委科研计划(2017KJ050);中国能源建设股份有限公司重大科技专项课题(CEEC2020-KJZX-10-01)。

摘　　要：基于有限长移动线热源(MFLS)传热模型,根据时空叠加原理推导出含水层非稳态过余温度解析解ΔT_(MFLS);在格子单松弛模型(LBGK)的演化方程中引入离散力源项,建立格子Boltzmann法(LBM)与ΔT_(MFLS)的耦合计算模型。通过热响应实验,验证耦合计算模型与求解方法的正确性。研究表明,在不同计算工况下含水层各区域的渗流速度均具有启动—下降—回升—稳定的4个连续阶段。随着含水层孔隙率的降低,虚拟流体粒子动能损失增大,渗流速度降幅增大,回升过程缓慢。然而随着进水流速的提高,孔隙率的变化对于含水层水动力场演化过程的影响程度减弱;含水层热量运移过程的方向性显著增强。Taking the multiple moving finite line heat sources(MFLS)operation mode as the research object,the unsteady state analytical solutions of the excess temperature in aquiferΔT_(MFLS)are obtained by applying space-time superposition principle,based upon the transient moving finite line heat source model.A calculation model coupled lattice Boltzmann method(LBM)andΔT_(MFLS)is established,according to introduce the discrete force source term into the Lattice Bhatnagar-Gross-Krook(LBGK)evolution equation of the velocity distribution function.The coupled model and the calculation method are validated by the data determined from the in-situ thermal re-sponse test.The results show that coupled LBM-ΔT_(MFLS)simulation can well reproduce the evolution process of seepage velocity in aquifer which presents four successive stages of starting,declining,rising and stabilizing.With the decreasing of the porosity of aquifer,the movement space of virtual fluid particles reduces,consequently the collision probability of inter-particles and between particles and porous media increases,which leads to the kinetic energy loss enhancing,the descend range of seepage velocity increasing,and the recovery process retarding.While the inflow velocity increases,however,the influence degree of porosity on the hydrodynamic evolution of aquifer weakens.Meanwhile the directivity of the heat transport process in the aquifer enhances significantly.

关 键 词：地热能 传热 多孔介质 数值模型 含水层 格子BOLTZMANN方法 

分 类 号：TK529[动力工程及工程热物理—热能工程]