1. 引言
路面积雪结冰,给交通出行带来巨大的困扰,甚至导致严重的交通事故[1]。循环流体换热系统于路面结构中水平布置换热管,通过换热管中的循环热流向路面传递热量融雪化冰,可有效保障交通顺畅和防止交通事故[1] [2]。循环流体换热系统通过地源热泵提供热源;地源热泵在建筑供暖和制冷等领域都有广泛的应用,技术相对比较成熟;而循环流体换热系统缺少大量的工程应用,技术体系建设尚未完成[3] [4]。
模型试验和现场试验是研究循环流体换热系统的最直接方式,如:Tan Y Q等[5]铺设了一段循环流体加热的机场融雪路面,进行融雪效果的研究;王华军[3]建立了小型循环流体融雪试验模型,对融雪效果随时间的变化规律展开分析;林艳艳等[6]开展了路面循环流体融雪化冰现场实验研究,就布管形式和布管参数对融雪效果的影响进行了初步的分析;陈鑫等[4]通过现场试验,就换热管的位置对桥面循环流体融雪化冰效果的影响进行了较为深入的研究。循环流体换热系统的模型试验和现场试验需要搭建复杂的试验模型,耗费大量的时间、人力、物力和技术成本,限制了大量对比试验的开展和试验数据的全方位全过程的采集。
有限元数值分析方法针对求解域的传热理论模型(偏微分方程组),将计算域离散为有限数量的单元,将连续域的偏微分方程组求解问题转化为离散域的线性方程组求解问题,可以实现循环流体换热系统全方位、全过程的计算和采集[7]-[12],如:Chen X等[13]基于现场试验搭建了路面融雪化冰的二维瞬态有限元数值分析模型,精确刻画了路面融雪化冰的传热过程;彭建国等[8]建立了机场跑道的融雪化冰三维有限元分析模型,完成了传热过程的模拟,同时分析了温度场随时间的变化。总体而言,应用循环流体换热系统进行融雪化冰过程的数值模拟分析,相对比较成熟[14] [15]。
可见,现有研究应用模型试验、现场试验和数值模拟计算等方法,对循环流体换热系统进行较为深入的研究。但是,一方面,相对于路面,桥面和机场跑道等更加容易积雪结冰,且导致的交通事故更加严重,所以现有研究中的循环流体换热系统主要布置于桥面结构或机场跑道中,布置于公路路面结构中的较少;另外一方面,为了配合桥面板和机场跑道结构中钢筋网的布置,现有研究中换热管的布管形式主要为曲流式布管和蛇式布管,存在着温度场分布不均匀、融雪效果不佳等问题。
鉴于此,本研究提出一种新型的、与路面结构匹配的双螺旋布管形式,以提升循环流体换热系统融雪化冰的均匀性、融化效果和融化效率;以现场试验为依托,构建路面循环流体换热系统的精细化有限元数值分析模型,验证双螺旋布管方式的有效性和优越性,掌握其结构参数的技术特征,进一步论证其工程适用性。
2. 循环流体换热系统数值分析模型的构建
2.1. 循环流体换热系统的工作原理
在路面结构的面层、基层或垫层中埋设循环流体换热管;通过地源热泵将远高于环境温度的流体泵入换热管;流体在换热管中循环流动向周围介质传递热量,热量传至路表进行融雪化冰;失去热量的流体从换热管的出口流出,如图1所示。
Figure 1. Circulating fluid heat exchange system operating principle
图1. 循环流体换热系统工作原理
路表面持续从循环流体中吸收热量,同时由于融雪化冰作用持续向大气中释放热量(路表热通量);当循环流体提供的热量大于融雪化冰所需的路表热通量时,换热系统能够实现融雪化冰;反之,则无法实现融雪化冰。
循环流体在流动的过程中逐渐损失热量,温度逐步降低;距离换热管的入口越近,流体温度及其周围介质温度越高;反之,温度越低,如图1所示。所以,在循环流体换热管及其周围介质中形成了高温段(区)和低温段(区)。换热管的布置,要克服温度场的不均匀分布对融雪化冰效果的影响。为了适应钢筋网的布置,工程中常用的循环流体换热系统的布置形式为曲流式或蛇式,如图1所示。曲流式和蛇式循环流体换热系统中,循环流体的流入路径(高温段)与流出路径(低温段)分离,导致周围介质的高温区与低温区也发生分离,温度场不均匀,影响融雪化冰的效果。
因此,本研究提出一种新型的双螺旋布管形式。双螺旋布管条件下,循环流体由布管区域的右下角呈螺旋状由外向内流动,到达布管区域中心后再呈螺旋状由内向外流动,最终从右下角流出,如图1所示。双螺旋布管条件下,高温段和低温段交错布置,克服了曲流式和蛇式布管条件下的高/低温段分离、局部温度集中的问题;换热管在两个水平方向的布置形式相同,且都是对称布置,呈正交对称的特征,克服了曲流式和蛇式布管在不同方向的正交异性导致的温度场分布不均匀。
2.2. 有限元模型的构建
文献[4]中于某公路段进行了融雪化冰现场试验,采用的是曲流式布管形式。应用有限元数值分析软件COMSOL进行循环流体换热系统数值分析模型的构建,如图2所示。试验路段路面结构的长度和宽度为15 m,路面高0.81 m,由上到下依次为18 cm的沥青面层、43 cm的素混凝土基层和20 cm的砂砾垫层,如图2(a)所示。曲流式循环流体换热系统布置于路表面以下15 m处的沥青面层中,布管间距为15 cm;换热管为外径25 mm、内径21 cm的聚乙烯升温管(PERT),如图2(b)所示。
(a)
(b)
(c)
Figure 2. Finite element model of circulating fluid heat exchange system. (a) Road layer structure and boundary conditions; (b) Heat exchange tubes; (c) Mesh generation
图2. 循环流体换热系统的有限元模型。(a) 路面结构和边界条件;(b) 换热管;(c) 网格划分
根据现场数据资料[4],路面结构、换热管和流体的初始温度取最低环境温度−5℃。进行模型验证时,换热管入口温度按照图3所示的试验入口温度进行设置,入口流量取为1.1 m3/h [4];后续分析中,流体入口温度设置为40℃,入口流量是设置为1.1 m3/h。为了保证融雪化冰的效果,会对路面结构的两侧和底部进行隔热处理,所以这几个表面设置为绝缘边界,如图2(a)所示。道路表面设置为热通量边界条件。现场试验中没有进行热通量的测试,根据ASHRAE (美国采暖制冷与空调工程师协会)出版的《2019 ASHRAE Handbook》中推荐的Chapman模型[8]来计算热通量,计算得到路表面热通量q0值为108.68 W/m2。
Figure 3. In-site monitored inlet/outlet temperature and simulated outlet temperature
图3. 现场监测进/出口温度与模拟出口温度
管道尺寸相对于路面结构是很小的;进行管道传热模型构建时,忽略管道中流体流动的状态,将管道简化为一条空间曲线,降低网格剖分难度。因此,对管道采用极细化网格进行剖分,对周围路面结构采用常规尺寸进行自由四面体剖分,管道与周围路面结构的网格尺寸采用自由过渡的方式进行衔接;将路面结构分为3层,对顶面进行网格加密处理,以增加路面温度场提取的准确性;对底部进行网格稀疏处理,以减小计算量;按照上述方法进行模型的网格划分,共划分出164,433个单元,如图2(c)所示。图2(c)所示的模型是曲流式循环流体换热系统的数值分析模型;进行蛇式或双螺旋循环流体换热系统模型构建时,只需要对图2(b)的换热管布置形式进行替换即可。
2.3. 传热控制方程
循环流体换热系统的热传递路径主要包括了换热管流体流动过程中的管道传热、管道流体到周围介质的耦合传热、周围介质中的固体传热和路表融雪化冰的传热,本质上是一个多物理场耦合的求解问题。如前所述,本研究通过Chapman模型[8]计算融雪化冰所需的路表热通量,将计算所得的热通量作为路表边界条件来表征路表的融雪化冰传热过程。因此,循环流体换热系统的传热控制方程只剩下换热管流体的管道传热、路面结构的固体传热和二者之间的耦合作用。
换热管内流体的管道传热理论模型如公式(1)所示。
(1)
式中:
–流体密度,kg/m3;
–管道内径表示的截面积,m2;
–流体比热容,kg/m3;
为时间;
表示温度,℃;
–流速单位向量;
表示流速方向的变量变化率;
–流体导热系数,W/(m∙k);
–Darcy摩擦因子,取默认值;
–管道长度,m;
–流体流速,m/s;
–管壁内外部热交换,W/m;Q表示一般热源,W/m。
路面结构的固体传热理论模型如公式(2)所示。
(2)
式中:
–材料的密度,kg/m3;
–固体比热容,J/(kg·K);
–温度对时间的偏导数,表征温度场随时间变化情况;
–固体内的平均流速,m/s,本模型中为0;
–变量的梯度,表征变量的空间变化率;
、
、
–分别为
、
、
方向单位向量;
–热流密度,W/m2,其计算如公式(3)所示;
为热源,W/m3。
(3)
式中:
–固体导热系数,W/(m∙K)。
公式(1)中的
表示管道内壁处传递的热量,其计算如公式(4)所示。公式(4)中,
表示管道内壁的温度,与公式(1)中的
一致;
表示管道外壁的温度,与公式(2)中的
一致。将公式(1)计算得到的管道内壁处流体的温度T代入公式(4)中的
,将公式(2)计算得到的管道外壁处的温度T代入公式(4)中的
,即可计算得到
;再将计算得到的
代入公式(1),即完成了固体传热物理场与管道传热物理场的耦合。
(4)
式中:
–管道外壁温度,℃;
–管道内壁的温度,℃;
–等效热阻,由管道结构的热物性参数计算得到。
公式(1~4)中的参数包括路面各层结构、管道结构和管道流体的密度、比热容和导热系数等,其设置如表1所示。
Table 1. Material parameters for numerical models
表1. 数值模型的材料参数
参数 材料 |
导热系数[W/(m∙K)] |
比热容[J/(kg∙˚C)] |
密度(kg/m3) |
面层 |
1.4 |
900 |
2400 |
基层 |
2.2 |
850 |
2500 |
垫层 |
1.55 |
1100 |
2240 |
管道结构 |
0.44 |
2300 |
950 |
管道流体 |
0.6 |
4182 |
998.2 |
对于构建的循环流体换热系统有限元分析模型,计算得到的换热管出口温度与现场试验的出口是吻合的,如图3所示。模型计算结果与试验结果的误差保持在5%以内,满足工程应用要求。
3. 布管形式的对比分析
3.1. 温度场分析
系统运行5 d后,双螺旋、曲流式和蛇式等三种布管布置形式条件下的路表和埋管层的温度场如图4所示。
双螺旋布管条件下,换热管交错布置,所以周围介质中的高低温分布也是交错呈现的;埋管层的高/低温区由外向内交错分布,温度逐步降低;没有出现明显的温度局部集中和相对分离的现象。曲流式布管条件下,换热管的流出段始终处于流入段的内侧,流入段从左向右流动,流出段从右向左流动,所以埋管层的左侧温度明显高于右侧,上侧温度高于下侧;流入段区域和流出段区域都存在着严重的温度局部集中现象,温度场突变严重;高低温区域相对分离。蛇式布管条件下,换热管的流入段始终处于换热区域的下侧,流出段始终处于换热区域的上侧,所以埋管层下侧的温度显著高于上侧的温度,下侧高温集中,上侧低温集中;高低温区域相对分离严重。双螺旋布管条件下的高低温最大差异值(40.4℃)显著小于曲流式(41.97℃)和蛇式(43.17℃),进一步说明双螺旋循环流体换热系统的温度更加均匀。
双螺旋布管条件下,换热管是正交对称布置的,所以温度场也呈现正交对称的特征,埋管层任意两个垂直方向上的温度场是相同的。曲流式布置条件下,换热管的布置是正交异性的,任何1个水平方向的温度场不对称,且2个水平方向的温度场的分布特征也不同,所以埋管层的温度场极为扭曲,高低温区域相对分离严重。蛇式布管条件下,换热管的布置是正交异性的,左右侧呈现一定的对称性,上下侧完全不对称,埋管层的高低温区域相对分离严重。
路表层中,由于埋管层与路表中间介质的“差异放大”效应,温度场的不均匀性在路表表现得更加明显。双螺旋布管区域外侧介质的温度高于内侧介质的温度,温度场由外向内均匀减小,但是减小的幅度不大;双螺旋布管的温度场整体比较均匀。曲流式布管区域左侧的温度显著高于右侧、上侧温度显著高于下侧,温度场更加扭曲,且不易通过工程措施进行控制。蛇式布管区域的温度由下至上逐步降低,呈现明显的两极化差异;相对于双螺旋和曲流式布管,蛇式布管形式的温度场更加难以控制。双螺旋布管条件下的最高温(9.53℃)低于其曲流式(9.94℃)和蛇式(11.8℃),最低温(−3.46℃)高于曲流式(−4.2℃)和蛇式(−4.6℃),说明路表层双螺旋循环流体换热系统的温度差异更小,温度场更加均匀。
(a) 双螺旋
(b) 曲流式
(c) 蛇式
Figure 4. Temperature fields under different tube arrangement conditions/˚C
图4. 不同布管形式条件下的温度场/℃
3.2. 融雪化冰效果分析
通过换热功率(Ps)、路表平均温度(Ts)、路表平均温度达到0℃的运行时间(Et)和融化比(R)等指标对融雪化冰效果进行评价。换热功率Ps表征换热管单位时间内传递热量的能耗,由管内流体的质量流速、流体的比热容和进/出口温度等指标计算而得,如公式(5)所示。Ts表示路表各处温度的平均值,反映换热作用下路表的整体升温效果。R表示路表温度大于0的面积与整个面积的比值,直观表示换热管融雪化冰的效果。Et表征换热管融雪化冰的换热效率。
(5)
式中:
为换热功率,W;m为管内流体的质量流速,kg/s;
为流体的比热容,J/(kg∙℃);Tin、Tout为入口、出口温度,℃。
图5为不同布管形式条件下,Ps、Ts和R等指标随时间的变化关系曲线。对于三种布管方式,系统开始运行的瞬间,Ps是很大的;随着时间的增长,Ps逐渐减小,且最终趋于稳定;在系统运行的初始阶段,Ts会急剧增大,并且很快到达0度;随着运行时间的增加,Ts的增速会放缓,且最终趋于稳定;R与Ts的变化规律相似。
(a) (b)
(c)
Figure 5. Heat transfer process under different tube arrangement conditions, (a) Ps; (b) Ts; (c) R
图5. 不同布管形式条件下的换热过程,(a) Ps;(b) Ts;(c) R
由图5(a)可知,在传热初期,入口温度远高于初始值,3种布管形式条件下的Ps均为81.5 kW;随着系统的运行,出口温度迅速升高,埋管进出口温差变小,Ps逐渐降低;在不同的时间节点,蛇式布管条件下的Ps最大,双螺旋布管条件下的Ps最小。同样的入口温度条件下,出口温度越小,Ps越大,表明系统的能耗越大:相对于曲流式和蛇式布管,双螺旋布管条件下的能耗分别降低3.21%和12.35%。
由图5(b)可知,不同的时间节点上,蛇式布管的Ts最大,双螺旋布管的Ts最小,但是差异不是很大;蛇式布管、曲流式布管和双螺旋布管等条件下的Et分别为8.3 h、9.1 h、10 h。
由图5(c)可知,在系统运行的初期,曲流式布管条件下的R最大,蛇式布管条件下的R最小;在系统运行的中后期,双螺旋布管条件下的R最大。系统运行1 d后,曲流式、双螺旋和蛇式布管条件下,R分别为62.92%、59.05%和54.71%;系统运行5 d结束后,曲流式、双螺旋和蛇式布管等条件下,R分别达到90.32%、92.28%和75.08%;相对于曲流式布管和蛇式布管,双螺旋布管条件下的融雪面积比分别提升2.16%和22.91%。可见,短期内曲流式布管条件下融化速度最快,长期时间内双螺旋布管条件下的融化效果更好。
综上所述,3种布管形式各有优劣,蛇式布管条件下路表平均升温快,曲流式布管条件下短期内融雪效果好,但是双螺旋布管换热系统的优点更为明显:整体温度场更加均匀,可提高换热效果的同时减少不均匀温度应力导致的结构破坏;以较低的换热能耗带来更好的融雪效果,效率更高;长期融化效果最好,短期融化效果也较好。可见,双螺旋布管明显优于工程上常用的曲流式布管和蛇式布管。
4. 双螺旋循环流体换热系统结构参数分析
4.1. 布管深度对换热效果的影响
图6为循环流体换热系统在运行稳定时,不同布管深度(10、15、20、25、30 cm)条件下的各项评价指标。由图可知,随着布管深度的增加,Ps减小、Ts减小、R减小、Et增大。随着布管深度的增大,Ps和Ts随着布管深度减小的曲线形式类似于双折线:布管深度 ≤ 20 cm时,下降的梯度保持不变;布管深度 ≥ 20 cm时,下降的梯度减小。R随布管深度减小的曲线呈现明显的两段特征:布管深度 < 20 cm时,R保持不变;布管深度 ≥ 20 cm时,R急剧下降,下降的梯度逐步减小。Et随布管深度增大的曲线类似线性曲线,上升的趋势保持不变。可见,由于双螺旋循环流体换热系统呈平面形式布置,布管深度越大,循环流体换热系统与路表的传热路径越长,路表热通量对换热管的影响越小,所以循环流体换热系统的换热能耗减小,路表融雪状态变差、融雪效率降低。
Figure 6. Tube depth arrangement effect on heat transfer
图6. 布管深度对换热效果的影响
4.2. 布管间距对换热效果的影响
图7为循环流体换热系统在运行稳定时,不同布管间距(10、15、20、25、30 cm)条件下的各项评价指标。由图可知,随着布管间距的增大,Ps、Ts和R先增大后减小,变化曲线类似于反向的抛物线,抛物线的拐点位于布管间距20 cm处;Et先减小后增大,变化曲线类似于抛物线,抛物线的拐点也位于布管间距20 cm处。Ps和Ts的反向抛物线型变化曲线的左右两侧相对比较对称,即:布管间距 < 20 cm和>20 cm时,Ps和Ts随着布管间距的增大而增大或减小的幅度相近。R的反向抛物线型变化曲线的左右两侧不对称:R随布管间距增大的梯度远大其减小的梯度。Et的抛物线变化曲线的左右两侧不对称:左侧的下降梯度大于右侧的上升梯度。
Figure 7. Tube spacing effect on heat transfer
图7. 布管间距对换热效果的影响
双螺旋循环流体换热系统的进路与出路交错布置;布管间距越大,换热管支路之间的相互作用距离越大,进路与进路之间的“热干扰”效应越小,同时进路与出路之间的“热叠加”效应也越小——“热干扰”效应导致换热效果降低,“热叠加”效应换热效果提升。在布管间距 ≤ 临界值(20 cm)时,支管之间存在着较大“热干扰”效应;因此,随着布管间距的增大,支管之间的“热干扰”效应降低,循环流体换热系统的能耗增大、融雪化冰状态变好、融雪效率增高。在布管间距 ≥ 临界值(20 cm)时,支管之间的“热干扰”作用弱化,支管之间的“热叠加”效应占主导地位;因此,随着布管间距的增大,支管间的“热叠加”效应给换热效果的强化作用下降,进而导致系统的换热效果下降,也导致系统的能耗降低。从抛物线的形态可以看出,“热干扰”效应对系统换热效果的影响,大于“热叠加”效应对系统换热效果的影响。
4.3. 换热管尺寸对换热效果的影响
图8为循环流体换热系统在运行稳定时,不同换热管尺寸(换热管壁厚为4 cm,外径分别为20、25、32、40、50 cm)条件下的各项评价指标。由图可知,随着换热管外径的增大,Ps有一定程度的减小,但是减小的幅度不大;Ts呈线性减小;R呈非线性减小,其变化曲线呈双曲线形式,拐点位于外径 = 32 cm处;Et呈线性增大。
Figure 8. Heat exchange tube size effect on heat transfer
图8. 换热管尺寸对换热效果的影响
双螺旋循环流体换热系统的进口支路与出口支路交错布置,管内流体与周围介质的有效换热时间较长。换热管尺寸的变化,主要影响换热管内循环流体的流速:在循环流体流量(1.1 m3/h)不变的情况下,管径越大,循环流体流动的速度越小。在换热时间较长的情况下,流速的变化对换热效果的影响较大;管道流速降低,埋管循环流体向周围介质传递的热量减少,管道流体出口的温度增加,所以Ps减小;流速越小,流体与周围介质的换热效率更低,所以Et随之增大;同样的系统运行时间内,换热效率越低,则换热效果越差,所以Ts和R也随之降低。由图8可知,换热管外径的增加,造成Ts和Et等指标的变化较小,但是导致R发生了约35%的降低——如果系统运行时间继续增加,其他3个指标的变化幅度也会增大。
5. 结论
本研究提出一种新型的布管形式——双螺旋式,以克服循环流体换热系统融雪均匀性差、稳定性差等问题。基于典型工程案例,构建双螺旋循环流体换热系统的精细化三维有限元数值分析模型,分析其有效性及其主要结构参数的技术特征,得到以下的主要结论:
(1) 相对于常用的曲流式和蛇式,双螺旋布管条件下循环流体换热系统的流动路径更合理,温度场更均匀、更对称,以较低的换热能耗带来更好的融雪效果和更高的融雪效率。
(2) 布管深度越小,路表热通量对换热管的影响越大,双螺旋循环流体换热系统的换热能耗越大、换热效果越好、换热效率越高。双螺旋布管形式存在着“热干扰”弱化效应和“热叠加”强化效益:布管间距小于临界值的时候,换热效果随着布管间距的增大而提升;布管间距大于临界值时,换热效果随布管间距的增大而降低。换热管流量恒定的条件下,换热管断面积的增大,导致换热效果的大幅降低。
基金项目
山区道路工程与防灾减灾技术国家地方联合工程实验室开放基金(SQDL-2023-02);重庆市教委科学技术研究项目(KJQN202301530);重庆市自然科学基金面上项目(CSTB2023NSCQ-MSX0828)。
NOTES
*通讯作者。