1. 引言
物质系统内各点温度的集合称为温度场,反映温度在时间和空间上的分布,是时间和空间的函数 [1] 。关于温度场的研究对我们深入了解与解决传热相关问题至关重要且有较强的理论与实际意义。由于数值模拟对实际的物理、化学、工程等问题能进行预见性、合理性、可能性的模拟处理,能为解决实际问题提供启示性信息 [2] ,再加之使用传统实际实验方法耗时费力且材料成本较高等原因,越来越多地温度场问题使用数值模拟的方法进行解决。Ansys软件具有强大的温度场模拟仿真功能,是温度场计算的实用工具,如运用Ansys软件进行实际加热过程的模拟得到可靠结果 [3] [4] ;运用Ansys软件可对混凝土浇注过程的温度场和温度应力进行仿真计算 [5] 等等。本课题启发于《热学》中的热力学第零定律 [6] ,《热学》中的热平衡是指在没有外界条件影响下热力学系统的温度处处相等,有特定的温度场反映系统内部分子无规则热运动情况。系统在传热与加热过程中,非平衡态下不同区域的温度场也反映系统内部分子无规则热运动剧烈程度不同的机制,系统在平衡态和非平衡态时有不同特点的温度场,分析在传热与加热过程中不同温度物体(或物体的不同温度部位)从非平衡态到达平衡态的微观过程很有意义。本文简化模拟环境,未采用流固耦合模拟,主要是结合Siemens UG与Ansys软件完成对稳定热源 [7] 在一维热传导下加热不同质量的圆柱、棱柱形纯水的温度场进行数值模拟,分析温度场的分布特征和变化规律。由于《热学》教学中传热过程较为抽象,难于形象讲解,而本文对温度场的模拟结果使抽象的传热过程具体、直观明了的呈现,可使学生深入理解传热微过程,对热学教学极有帮助,同时通过模拟分析,还可为加热器形状的择优选择提供参考信息。
2. 数值模拟内容及方法
通过稳定热源在一维热传导下加热质量为1 kg、2 kg、3 kg、4 kg、5 kg几何形状为圆柱体、长方体、五棱柱形纯水,对其温度场和热通量进行数值模拟。
2.1. 纯水三维几何模型建立
采用Siemens UG软件建立三维模型,建模时保持底面积不变,利用控制变量法研究对纯水加热的温度场和热通量随几何体形状和质量变化的特点。圆柱体、长方体、五棱柱模型尺寸设定见表1。

Table 1. The relevant size settings of pure water of cylinder, cuboid, and pentagonal prism
表1. 圆柱形、长方体、五棱柱纯水的相关尺寸设定
2.2. 求解方法及边界条件设定
本文利用Ansys Icem对几何体进行网格划分,使用拓扑结构逼近几何模型,再在拓扑结构上生成六面体网格,并对局部网格细化,优化计算结果。本文纯水温度场模拟采用Ansys稳态热分析进行,系统初始温度为20˚C,密度9974 kg/m3,比热容4200 J/(kg∙k),导热系数0.604 w/(m∙k),热量垂直于底面流入,热量传导为一维,大小为20 w/m2。
3. 计算结果与分析
前提分析:一、选择质量为1 kg、2 kg、3 kg、4 kg、5 kg的纯水原因便于对不同质量的纯水模拟结果进行比较;二、选择圆柱、长方体、五棱柱纯水原因:1) 当纯水几何体侧壁设为绝热时,选择圆柱和棱柱形纯水进行加热,仅考虑沿轴向的一维传热情形,使模拟环境简化。2) 加热圆柱、长方体、五棱柱纯水便于对圆柱和棱柱纯水的加热温度场、轴向热通量的特点进行比较和不同棱柱纯水的棱边散热效应的比较。限于篇幅,本文只列出质量1 kg、2 kg时圆柱体、长方体、五棱柱纯水温度场模拟图,其余质量和形状的温度场、轴向热通量模拟图及数据略。
3.1. 温度场分析
3.1.1. 对一维稳定热源加热不同质量和几何形状纯水的温度场模拟图
如图1和图2所示。

Figure 1. The chart of the temperature field of 1 kg pure water of cylinder, cuboid, and pentagonal prism
图1. 1 kg圆柱、长方体、五棱柱纯水温度场图

Figure 2. The chart of the temperature field of 2 kg pure water of cylinder, cuboid, and pentagonal prism
图2. 2 kg圆柱、长方体、五棱柱纯水温度场图
3.1.2. 温度场模拟图数据
对稳定热源在一维热传导条件下加热质量为1 kg、2 kg、3 kg、4 kg、5 kg的圆柱、长方体、五棱柱形纯水的温度场模拟图进行数据整理,见表2~4。

Table 2. The simulated data of the temperature field of 1 kg pure water of cylinder, cuboid, and pentagonal prism
表2. 1 kg圆柱形、长方体、五棱柱纯水的温度场模拟数据

Table 3. The simulated data of the temperature field of 2 kg pure water of cylinder, cuboid, and pentagonal prism
表3. 2 kg圆柱形、长方体、五棱柱纯水的温度场模拟数据

Table 4. The relevant processing data of the temperature field of different qualities pure water of cylinder, cuboid, and pentagonal prism
表4. 不同质量圆柱形、长方体、五棱柱纯水的温度场相关处理数据
3.1.3. 温度场分析
图1和图2中,一种色带表征一个温度值。同质量的圆柱、长方体、五棱柱形纯水的温度场,相邻色带温度的差值为固定值即1 kg时为0.368 K,2 kg时为0.736 K (其余质量的情况本文未列出),传热方向长度见表2、表3即1 kg时为0.01111 m,2 kg时为0.02222 m,即得出温度梯度33.1233/(K/m)相同,说明在纯水初、终温度设定的加热过程中,不同质量、不同形状的纯水体从加热区到整体温度达到一致的温度变化阶段相同,这是因为不同质量、不同形状的纯水体从加热区到整体温度达到一致时内在机制是一样的;从图1和图2还看出,圆柱形纯水的温度分布场图为交叠浑浊的状态,温度梯度条带呈现渐变情形,而长方体、五棱柱等纯水的温度分布场图则界限清晰、分段明确,温度条带呈现跃变情况,这表示同质量纯水的初、终温度设定加热过程中,从加热区到圆柱纯水整体温度达到一致的速度较慢,而从加热区到棱柱纯水整体温度达到一致速度较快,出现此情况是由于盛水容器形状的不同,液体层间热传递时温度变化有梯度差别,圆柱形纯水的液体分子“相互搅拌”(传热)较快,整体温度达到一致的时间较短,而棱柱形纯水液体分子“相互搅拌”(传热)较慢,整体温度达到一致的时间较长的原因。
3.2. 轴向热通量分析
3.2.1. 对不同质量和几何形状纯水加热的轴向热通量模拟图
如图3和图4所示。

Figure 3. The chart of the axial heat flux of 1 kg pure water of cylinder, cuboid, and pentagonal prism
图3. 1 kg圆柱、长方体、五棱柱纯水的轴向热通量图

Figure 4. The chart of the axial heat flux of 2 kg pure water of cylinder, cuboid, and pentagonal prism
图4. 2 kg圆柱、长方体、五棱柱纯水的轴向热通量图
3.2.2. 轴向热通量分析
从图3、图4看出,用稳定热源在一维导热条件下对不同质量圆柱、长方体和五棱柱纯水分别加热时,它们的热通量总体趋势均沿轴向且较稳定,但加热同体积(本文水的密度相同,也可以认为是同质量)底面为四边形、五边形等棱柱形纯水体时有明显的棱边散热显示,随着棱边增多,侧面散热量增大,而加热同体积底面为圆形时的圆柱纯水体时,侧面散热量较小。这是由于加热同体积的棱柱形纯水体时,随着棱边增多,棱柱侧面积增大(即周边散热面积增大),热量集中性较差,轴向热通量较差,周边热量散失相应较多,而加热底面为圆形时的圆柱纯水体时,与同体积的棱柱纯水体相比,侧面积最小,即周边散热面积最小,热量集中程度较强,轴向热通量较集中,周边热量散失较少的原因。
4. 结论
本文使用Siemens UG和Ansys软件对稳定热源一维传热条件下加热不同质量、不同形状纯水温度场、热通量进行模拟并作分析研究。研究结论:1) 同质量、不同几何形状纯水的温度分布和温度梯度不同,柱形纯水的温度分布场图为交叠浑浊状,温度条带呈现渐变情形,而棱柱形纯水的温度分布场图则界限清晰、分段明确,温度条带是跃变情况;2) 每种形状纯水的温度分布场和温度梯度不随质量而变;3) 同质量的圆柱和棱柱形纯水的温度分布场和温度梯度情况反映出加热圆柱纯水和棱柱形纯水时从加热区到纯水整体温度达到一致的速度不同,前者较快,后者较慢,这也是人们平时选择类圆柱容器对水加热的主要原因;4) 初、终温度设定时加热不同质量的各形状纯水至整体温度达到一致所经温度变化阶段均相同。5) 总热通量确定时,轴向热通量总体趋势相同,但在棱柱纯水的棱边处有散热效应,棱边越多此效应越强,圆柱纯水最弱;另外,加热棱柱纯水时的棱边散热效应随纯水质量增加而相应增强。通过模拟分析,使抽象的传热过程具体、直观、形象的呈现,可让读者深入理解传热微过程,对热学教学极有帮助,同时还为加热器形状的择优选择提供参考信息。
参考文献
NOTES
*通讯作者。