格林函数法在竖直u型地埋管传热计算中的应用 application of green function method in calculation on heat conduction of vertical u-tubes heat exchanger.pdfVIP

第31卷第2期 太阳能学报 V01．31．No．2 2010年2月 AcrAENERGU正SOLARISSINICA F出．．20lO 格林函数法在竖直U型地埋管传热计算中的应用 张长兴1，胡松涛2，李绪泉2 (1．西安建筑科技大学环境与市政工程学院，西安710055；2．青岛理工大学环境与市政工程学院，青岛266033) 摘要：以竖直u型地埋管传热的数学模型为基础，采用格林函数法确定了钻孔外IX．域的温度分布解析式。通过 求解钻孔壁温度，与经典的“负荷叠加”G函数法进行了比较。结果表明，该方法计算速度更快，结果精确。 关键词：哆直U型地埋管；传热模型；格林函数法；计算耗时 中图分类号：TK52 文献标识码：A 型地埋管传热的准确、快速计算提供了必要的参考。 0 引 言 竖直U型埋管换热器是应用最为普遍的一种地 1竖直U型埋管传热数学模型 源热泵换热器形式。目前，对于其换热性能的分析， 1．1数学模型 一般是以钻孔为界，分成两个计算区域：孔内区域按 u型埋管在土壤中的传热较为复杂，在进行如 照稳态传热计算…，孔外区域进行非稳态传热计算。 下简化的基础上，按照能量平衡的原则，建立了柱坐 这两个区域通过钻孔的壁温耦合起来，组成了整个 标竖直U型埋管的传热模型。 u型地埋管的数学模型。孔外区域的非稳态传热计 ①忽略热湿迁移的影响，认为土壤中的传热为 算可分为数值法和分析法两种类型，其中数值方法 纯导热I’口J题；②忽略钻孑L内填料沿竖直方向上的导 主要通过有限差分法或有限元法进行数值计 热，导热只沿径向进行。 算【2《]。而分析法则以“柱热源”或“线热源”理论进 1 aTa2T l ar ，．、 行计算[4’6圳。日前，采用分析法进行动态逐时负荷 i万2万+7万 u， 下的竖直u型地埋管传热计算时，充分考虑计算时 边界条件为： 刻前所有时刻的负荷对该时刻传热计算的影响，通 m0)=t；一J|}筹|，=r=掣 过进行“负荷叠加”的方法求解№’8’9J。此时，需结合 土壤特性、钻孑L的几何特性等因素计算在不同时刻 式中，兄——土壤的远边界温度，℃；矗——土壤的导 的响应系数被称为G函数№’81或g函数法[1’9]。该 热系数，w／(m·oc)；q，——单位钻孔长度的负荷， 方法在进行地源热泵系统长时间(1a甚至20a)的运 W／m；rb——钻孑L半径，m；a——土壤的导温系数， 行模拟时，计算时间很长。因此，为了缩短计算时 m2／h；t——地埋管换热的持续时间，h；r——土壤中 间，文献E6，9～11]中采用“负荷块法”进行计算，取 某点至热源的径向距离，m，r≥“； 得了令人满意的计算结果。 1．2传热模型的Laplace变换求解 本文通过引入格林函数法u2|，以“柱热源”模型 针对当q。为定热流q。时对应的数学模型，文 为基础，提出了一种地埋管外区域传热计算的新方 法。该方法只考虑前一时刻负荷对当前计算时刻的 g。作用下对应的温度分布，可用式(2)表示： 影响，与G函数法相比，缩短了计算时间，为竖直U 收稿日期：2008-06-12 基金项目：青岛If『科技发展计划项目(06-2-2-22-jch) 通讯作者：张长兴(19r77一)，男，博士研究生，主要从