欢迎光临苏州爱尔玛特环保节能科技有限公司!联系电话:13776178696
logo

行业资讯

Product categories

联系我们

  • 苏州爱尔玛特环保节能科技有限公司
  • 手机:0512-53371508
  • 电话:0512-53371508
  • 联系人:13776178696
  • 邮箱:aiermat01@126.com
  • 地址:太仓市城厢镇常胜北路33号10幢

3.8.2数值模拟和分析

日期:2022-01-07 09:10:13点击:394

3.8.2数值模拟和分析
1.物理模型和控制方程
数值模拟的目的是确定热管周围土壤的 温度场。模拟的范围及相关条件如图3. 8.3 所示。
在图3. 8. 3中,考虑一个在热管周围充满了冻土的圆柱体。轴对称的圆柱坐标为r, 热管轴线为Z坐标,并釆取下列假定条件:
① 在冻土内只有热传导过程存在,不考虑土壤中的相变和传质过程;
② 被模拟的土壤圆柱体的顶部和底部为绝热面;
③ 只考虑单支热管的存在;
④ 圆柱体的外径足够大,其外边界可考虑为等温表面,为原冻土温度 
;
⑤ 圆柱体的内边界为热管蒸发段外径为
的外表面,并假想在此界面上存在某种第 三类边界条件,而这一条件应由热管传热的 一系列热阻推导岀来。
考虑到上述条件,热管周围土壤圆柱体内的导热过程可归结为圆柱坐标下的非稳态 一维导热过程,其控制方程为

 
如图3. 8. 3所示,从土壤至空气的传热过程由7项过程组成:①由土壤至蒸发段外表 面之间的导热;②由蒸发段外表面至内表面之间的导热;③管内介质的蒸发传热;④从蒸 发段向凝结段的蒸汽的流动过程和热量的携带(该项热阻很小,可以忽略);⑤凝结段内 介质的凝结换热;⑥从凝结段管壁内表面向外表面的导热;⑦凝结段外表面和周围空气之 间的对流换热。每一项传热过程对应一项热阻,也就是说,热量从土壤传给大气要经过串 联的7项热阻。考虑到在7项热阻中,唯有第一项热阻,即土壤至蒸发段外壁的热阻为最 大,而且考虑到热管本身的热容量相对很小,因而可以假定:在任何瞬间,在蒸发段外表面 从土壤中吸收的热量将等于在凝结段外表面散失在空气中的热量。
上述各传热过程中的第2 ~7项,可依次表示为
式中,分别表示为蒸发段外壁温度,蒸发段内壁温度,蒸 发段蒸汽温度,凝结段蒸汽温度,凝结段内壁温度,凝结段外壁温度及管外大气温度, 分别为管内蒸发换热系数、凝结换热系数和管壁导热系数;必分别为蒸发段 长度和凝结段长度;。为热流(传热量)w
将式(2) ~ (7)各式相加,即可解出热流
可表示为
式中,是以蒸发段外表面为基准的相当换热系数,由式9)确定,它可以看作所模拟的 圆柱体土壤的一个边界条件,它代表在土壤和热管蒸发器界面之间的换热强度。
根据上述解释,与式(9)和控制方程1)相结合的初始条件和边界条件为
此处,H为所模拟的圆柱体的高度;是永 久冻土的初始温度。(此处取
2.气温关联式
大气温度的变化,直接影响热管的工作状 况,从而也直接影响着土壤温度场。为了得到 在式(8)及式10)中的空气(大气)温度
,需 要根据应用现场的气象资料,绘制1年中365 天的按天平均的气温变化曲线。由青藏铁路经 过的风火山地区在1995年的气象数据,绘制了 一年中的气温变化曲线,如图3. 8.4所示。
 

图3.8.4 一年中气温随时间的变化
3. 8. 4中的数据点代替每10天的平均气温,回归的曲线是一种三角函数,即一年 中的气温按三角函数的规律变化。其关联表达式为
式中,。代表从105日算起的天数。按照式11),最低的大气温度为,发生在 115日左右;最高气温为发生在720 -25日,当大气温度低于冻土温度 (-2℃)时,热管开始将土壤中的热量传给大气的传热过程。对于热管埋设后的第一个 年度,热管传热期间的长短可由温度为-2℃的一条水平线与气温曲线相交所形成的下部 区间来确定。一般将这一区间称为“寒季”,而将气温曲线高于冻土温度的区间,称为“暖 季”,由式11)可以算出,对于给定的这一气象条件,热管工作的寒季,从105日开始 至第二年510日结束,共215天,而暖季从510日至105日,共计145天,(每年 按360天计算)。应当指出,由该气温变化曲线所确定的寒季和暖季的时间段划分是理 想情况,事实上,随着热管周围土壤温度的不断降低,在寒季的末期,热管周围的土壤温度 已低于原冻土温度-2 ℃的水平,因此,在气温达到-2℃以前,即当气温低于-2℃时,热 管已停止了工作。所以实际的热管工作期(即寒季的长短)将少于由图确定的215天,可 能在200天左右,即在4月下旬(而非510日)热管工作的寒季就结束了。所以,在下 面的数值计算中,时间范围只计算到200天为止(即425日)。
由此可以推断,由于第一年土壤储冷降温的结果,在第二年寒季开始时的大气温度将 不再是
-2℃而是当时热管周围的实际土壤温度(小于-2℃),因而第二年的寒季比上 一年的寒季要缩短,依此类推。本文的数值模拟,将只考虑热管埋设后第一个寒季来临以 后的情况。
3.模拟方法和模拟结果
用控制容积法将控制方程离散化,并假定在控制容积中物性是均匀的。在式10) 中,相当换热系数
的计算至关重要,由式(9)可以看出,是多个变量的函数,其中,空 气侧的对流换热系数如,是一个控制因素。本身受多种因素的影响,应按翅片管管外换 热的有关公式计算,并选取由气象资料提供的当地风速。作为计算依据。在数值 计算中,选取的时间跨度是整个寒季,即气温低于-2℃的所有天数(200天)。
下面给出几个有代表性的模拟结果:
1)在整个寒季,不同半径处的土壤温度变化曲线计算结果如图3. 8. 5所示。
3.8.5中,选取土壤不同的半径坐标,从r = o. 076 mr=10. 088 mo横坐标为热 管运行的天数,概括了整个寒季200天,由图中的曲线可以看出:
① 在距离热管中心线1 m的范围内,随着寒季的到来,土壤各点温度随时间明显地下 降,在离热管很近的土层(r= 0. 076 m),土壤温度接近气温的变化,在100天左右,温度可 降至-15℃以下。
② 对于离热管较远的土层,例如r = 2 ~3 m处,在寒季开始大约30天后,才开始降温 的变化,在整个寒季,降温幅度很小,约1 ~2℃,最大降温发生在寒季即将结束的时段。 对于r>5 m的区域,热管的影响就很小了,开始发生变化的时间大约在2个月以后,发生 在寒季的末期,最大的温差在1℃左右。

 
(2)储冷的两个阶段。

计算表明,在整个200天的寒季中,土壤温度的变化规律是不一样的,可明显地分为 两个时间段,如图3.8.6中的a)、(b)所示。
3. 8.6 土壤温度的变化过程
3. 8.6(a)对应的时间段是当气温从-2℃逐渐降低到寒季的最低气温的时候所对 应的时间段,大约为100天(从105日至115日),在此时间段,各处土壤温度是逐 渐下降的,且壁面处土壤温度曲线的斜率随时间延长越来越大,说明土壤传热量即储冷量 逐渐增加,故这一时间段可称为“储冷段”。而对于图3.8.6(b)所示的第二时间段,对应 气温逐渐增加的情况,从寒季的最低气温升至寒季接近终了时的温度,时间段为从第100 天到第200天(对应从115日至425日)。在这一时间段,温度变化有不同的特点: 随着气温的回升,靠近热管的土壤温度逐渐升高,而远处的温度仍逐渐降低,温度曲线呈 现逐渐展平的趋势,故可将这一时间段称为“冷量扩散段”。
4.热管储冷量的计算
热管周围土壤在寒季终了时的温度分布代表了该时的总体致冷效果,热管周围5 m 半径内土体在一个寒季(200天)最终储冷量可表示为
其中,如叫分别为土壤的比热、不同半径处的温度和不同半径处的体积质量,取值如 下:(假定土质为亚黏土,含水量为17%)仏选取各控制容积上的 温度值;此处计算的土层高度为5 m, 土壤密度由上述计算结

果可知,热管在200天内的平均功率为82.6 W。单支热管在一个寒季中的总传热量(总 储冷量)相当于68 kg标准煤的发热量,这是一个很可观的效果。
应当指出,在第200天时的最后一条温度曲线是相当水平的一条曲线,这将是第一个 寒季结束时土壤温度的最终结果和分布形态。可以设想,在即将来临的暖季,土壤温度是 不会停止不变的,它还会继续展平,直到第二个寒季到来时,在原来储冷的效果基础上再 开始新的储冷过程。

上一篇:3.8.1应用背景

下一篇:3. 9.1分离式热管换热器的原理和优点

Copyright © 2020-2023 苏州爱尔玛特环保节能科技有限公司 All rights reserved 联系人:黄经理 地址:太仓市城厢镇常胜北路33号10幢 备案号:苏ICP备20027618号 技术支持: 核心词专业优化
  • 返回顶部
  • 13776178696
  • 0512-53371508
  • 微信公众号