双圆柱坐标下垂直埋管土壤热湿耦合迁移理论模型与数值模拟 theoretical model and simulation on heat and moisture transfer around u-vertical underground heat exchanger in double cylinder coordinates.pdfVIP

9期 李新国等：双圆柱坐标下垂直埋管土壤热湿耦合迁移理论模型与数值模拟 文章编号：02s．0096(20憾)09．1067-05 ． 双圆柱坐标下垂直埋管土壤热湿耦合迁移 理论模型与数值模拟 李新国1，李 伟1，李丽梅2 (1．天津大学机械工程学院热能工程系，天津300072；2．北京龙源冷却技术有限公司，北京100176) 摘要：以土壤热湿迁移理论为基础，建立双圆柱坐标系下u型垂直埋地换热器周围土壤内三维非稳态热湿耦合 迁移模型，并考虑土壤内热湿迁移的相互影响，编制相应的计算软件。进行不同土壤物性、不同运行工况、不同初 始含湿量对U型垂直埋地换热器周围土壤内热湿迁移影响的研究。将模拟值与实验结果进行比较，表明该模型可 用于实际应用的模拟分析。 关键词：U型垂直埋地换热器；土壤内热湿迁移；双圆柱坐标系 中图分类号：TKSl2 文献标识码：A 匀一致； 0前言 4)流体与固体骨架之间瞬间达到局部热平衡； 埋地换热器与周围土壤问通过热交换向土壤排 5)u型埋管与周围土壤接触紧密，相互之间不 放或提取热量，埋地换热器作用下土壤的热湿迁移 存在接触热阻； 有其特殊性【l刮。土壤内温度梯度是土壤水气运动 6)忽略粘性耗散和毛细滞后； 的主要推动力，温度直接影响土壤中水气的密度和 7)忽略土壤基质的压力势和溶质势，且不考虑 水气压，温度越高，水气压越高，土水势降低，土壤中 土壤内水分的宏观流动。 水分子间的作用力一氢键作用减弱，水分子的热运 1．2物理模型 动加剧，因此必然导致水从温度高处向温度低处移 实验研究表明，单U型埋地换热器周围土壤温 动。在温度梯度的作用下，土壤水以水汽状态由暖 度在距埋地换热器管中心半径10m以外可视为恒 端流向冷端，水汽运移到冷端凝结，使得冷端土壤含 定【4]。据此本文设定计算区域为以埋地换热器井孔 水量增高，土壤中水分发生迁移，土壤水分的迁移又 中心为圆心，10m为半径，深为90m的圆柱体，如图l 会使得土壤热物性发生变化，进一步影响了土壤的 所示。此处以外土壤温度为初始温度，不受埋地换 温度分布。因此，u型垂直埋管周围土壤内热湿迁 热器取、排热的影响。对于计算区域，选埋管中心为 移过程是相互耦合的三维非稳态传热传质过程。 z轴建立双圆柱坐标系，且取向下为正。柱坐标系中 1埋管周围土壤内热湿迁移物理数学 任一点由3个坐标给出：垂直坐标Z，半径r和角度 模型 节，如图2所示。且在z方向考虑土壤分层现象。 1．3网格划分 1．1模型假定 在双圆柱坐标系下对计算区域内进出水管周围 鉴于埋地换热器周围土壤内热湿迁移过程的复 土壤划分网格，将坐标系的竖轴与进出水埋管的中 杂性，本文作如下假设： 心线重合，在z方向按照lm间距划分网格；水平面 1)土壤固体骨架为刚性，不发生变形，不与液体 内的网格划分，考虑到在埋管周围温度变化较为剧 及气体发生化学反应； 烈，温度梯度较大，计算要求精度较高，网格划分较 2)水在土壤中为纯水，不考虑冻结现象，气体为 密，距离埋管较远的区域，土壤温度受埋管传热的影 理想气体； 响较小，网格划分较疏，桩孔内r方向采用等间距 3)将地下土壤分为若干个层，每一层内物性均 0．01m来划分网格，桩孔外r方向上的网格采用变长 收穆日期：20074)1．26 基金项目：国家高新技术研究发展计划(863)项目(2007AA057A32)；天津市科技计划项目(08ZCKFGX03600)