大洋勘探船岩芯储存库内温度场的数值模拟

郭月姣，吴毅昊，马双，陈旭，孟博，雷淑雅，石大川，冯国增*

1 江苏科技大学 能源与动力学院，江苏 镇江 212003

2 上海外高桥造船有限公司，上海 200137

3 中国船舶及海洋工程设计研究院，上海 200011

4 重庆大学 土木工程学院，重庆 400045

0 引 言

岩芯是研究地下地质和矿产情况的重要材料。岩芯开采出来后，经清洗、干燥、密封等过程，放置在大洋勘探船的岩芯储存库中。岩芯热容较大，在入库或出库时温度变化小，具有种类多、岩性变化大、灵敏度高等特点，因此储存条件非常严格[1]。当岩芯种类和研究目的不同时，对库温的要求也不同，通常松软的沉积物岩芯样品、冰心样品等需进行冷冻保存，而固体矿产岩芯一般进行冷藏储存。例如，“东方红3”号勘探船岩芯储存库分为样品冷冻库和样品冷藏库，库温分别为-20～-18 ℃和-4～10 ℃；“科学”号考察船的岩芯储存库库温分别为4 ℃和-20 ℃；“雪龙2”号科考船岩芯储存库分为低温样品冷藏室和样品冷冻室，库温分别为2～10 ℃和-20～-18 ℃。在岩芯储存库设计时，若库内冷风机放置位置不当，存在温度死区，造成温度分布不均匀或温度过高过低，都不利于岩芯储存，故需研究岩芯库冷风机合理的布置位置以长时间、高质量保存岩芯。

岩芯热容较大，对库温影响大，岩芯储存库的设计必须满足岩芯存储目标温度的均匀性，但是大洋勘探船岩芯储存库的运输机制十分复杂，很难以经验法来确定最优的冷库设计方案。随着计算机运算能力的提高和计算流体动力学（CFD）的发展，数值模拟作为方案择优的一种辅助手段，被广泛应用于船舶设计领域，例如舱室通风[2]、热量变化过程分析[3]、船型优化[4-5]及设计[6]等方面。同样，在空调系统设计中，数值模拟方法也得到了普遍应用。权崇仁等[7]和敖庆章等[8]通过模拟计算舰船典型舱室的温度场分布，验证了空调系统布置的合理性；亓海青等[9]采用数值模拟方法对舱室内部空气环境的舒适性进行了分析；倪崇本等[10]运用CFD技术对液化天然气船（LNG船）机舱内通风气流进行了数值模拟，消除了通风死角产生的安全隐患；Wang等[11]运用CFD技术模拟了制冷柜中的空气三维温度变化，有效改善了制冷系统性能；Söylemez等[12]和Rusly等[13]通过CFD模拟提供的更完整流场信息，确定了制冷系统中设备构件的最佳设置位置。

但是，上述研究均未涉及大洋勘探船岩芯储存库中的数值模拟应用。本文拟针对大洋勘探船陆上模拟岩芯储存库，采用CFD技术研究库内温度分布，从而改变冷风机设置位置；通过分析3种不同冷风机设置位置的壁面温度分布，以确定有利于岩芯储存的冷风机布置方案。

1 岩芯储存库结构及冷风机布置方式

大洋勘探船陆上模拟岩芯储存库高、低温冷库外观如图1所示，长5 m，宽13.6 m，高2.5 m，总容量为170 m3，其中低温库库容量70 m3，长5 m，宽5.6 m，高2.5 m；高温库库容量100 m3，长5 m，宽8 m，高2.5 m。冷库隔热材料为聚乙烯泡沫等。高温库工况：负荷6.345 kW，送风量2.71 m3/s。低温库工况：负荷6.915 kW，送风量2.67 m3/s。高、低温库冷风机分别采用横向、纵向和斜角3种布置方式，如图2所示。

图 1 岩芯储存库外观Fig. 1 Exterior view of refrigeratory for core storage

图 2 冷风机布置方式Fig. 2 Layout of air coolers

2 冷库数值模型

冷库流体流动需符合连续性方程、动量方程和能量方程。由于库内流体流动状态基本为湍流，故选用标准k-ε模型解决湍流模型的不封闭问题，且考虑重力影响。连续性方程、动量方程（以x方向为例）、能量方程、k方程和ε方程分别为：

3 模拟计算与实验条件

模型边界条件为：1)冷风机出风口流速设为冷库进风口速度；2)冷风机回风口采用压力出口；3)冷库壁面边界条件选用第2类边界条件，冷库库板保温绝热，热流密度为0，将壁面的传热量折算成内热源，高温库内热源为63.45 W/m3，低温库内热源为98.78 W/m3。室内设计参数及送风参数如表1所示。在模拟计算时假设高、低温冷库完全独立，2个库作为单独的库区来进行计算。以纵向吊顶布置为例，整个计算区域网格数约1.79×106，网格整体划分情况如图3(a)所示；对进风口网格进行加密处理，结果如图3(b)所示。

为了获得冷库内的合理流场，在建立描述冷藏库内流动换热的数学模型时，假设如下：1）冷库内空气为不可压缩理想气体，高、低温冷库之间无空气泄漏；2）鉴于流体为低速、不可压缩流动，故忽略由流体粘性力做功引起的耗散热；3）不考虑辐射换热；4）冷库外温度保持稳定。

表 1 模型参数Table 1 Model parameters

图 3 纵向布置网格划分Fig. 3 Mesh division of longitudinal layout

实验中，在冷库内布置温度传感器，温度数据由自主开发的数据采集系统显示并记录。选取高度方向2.2 m处进行温度测量，距离冷库壁面50 mm处等距安置温度传感器（其中高温库测温点有40个，低温库测温点有25个），测温点布置情况如图4所示。

4 模拟结果分析与实验验证

利用该模型研究岩芯储存库内部冷风机的最佳摆放位置，在空库和有吊顶情况下，冷风机分别采用纵向、横向和斜角3种布置方式，对大洋勘探船陆上模拟岩芯储存库进行分析，以确定冷风机的最佳放置位置。

当冷风机采用纵向吊顶布置方式时，模拟结果如图5所示，沿着纵向远离冷风机的区域温度越来越高，高温库内冷风机对侧壁面温度最高为9 ℃，低温库内冷风机对侧壁面温度最高为-14 ℃。当冷风机采用横向吊顶布置方式时，模拟结果如图6所示，低温库内冷风机出口斜下方温度较高，为-16 ℃，大部分区域温度分布较均匀，为-20～-18 ℃，偏离程度较低；高温库内温度在3～6 ℃内波动，少数区域的温度达到7 ℃。当岩芯储存库冷风机采用斜角（轴向45°）吊顶布置方式时，部分截面模拟的温度分布云图如图7所示，由于冷风机为轴向45°放置，库内另外2个角落风量不足，从而导致角落附近区域的温度较高，其中低温库角落的温度约-12 ℃，高温库角落的温度约9 ℃，从而不利于岩芯的储存。

由于温度控制精度为±1 ℃，所以高温库温度波动范围为3～5 ℃，低温库温度波动范围为-21～-19 ℃，温度在此范围内视为合格区域。选择截面z=2.2 m处进行温度对比分析，结果如图8所示。图中，分组1为合格区域，分组2为偏离设计温度2 ℃范围内，分组3为偏离设计温度3～8 ℃范围内，分组4为偏离设计温度9 ℃及以上。冷风机采用横向吊顶布置方式时，最大偏离温度不超过8 ℃，且温度大都集中在合格区域，面积为33 m2；斜角吊顶时合格区域面积为36 m2，但冷风机采用纵向和斜角吊顶布置方式时，温度偏差9 ℃的区域面积较大，面积分别为9和7 m2，不利于岩芯储存。综上所述，最终选择将冷风机横向吊顶布置。

当冷风机横向吊顶布置时，在高、低温库稳定时测量温度，并与模拟结果对比（图9）。由图可见，各测温点的实验与模拟温度之间存在1～2℃的温差，且在大部分测温点两者误差小于10%，表明数值模拟结果较准确。而在如图9（a）所示的测温点12～15处，实验与模拟温度的误差大于10%，推测是温度传感器距离冷风机送风口较近，使其位置发生偏离导致测量结果不准确。

图 4 温度传感器布置Fig. 4 Temperature sensor layout

图 5 纵向吊顶布置时温度分布云图Fig. 5 Temperature contours of longitudinal ceiling layout

图 6 横向吊顶布置时温度分布云图Fig. 6 Temperature contours of transverse ceiling layout

图 7 斜角吊顶布置时温度分布云图Fig. 7 Temperature contours of oblique ceiling layout

图 8 z=2.2 m处温度对比结果Fig. 8 Temperature comparison results at z=2.2 m

5 结 语

本文建立了陆上模拟岩芯储存库的数值模型，并进行了数值模拟与实验分析。根据大洋勘探船陆上模拟岩芯储存库的实际温度需要，对纵向吊顶、横向吊顶和斜角吊顶3种冷风机布置方式进行了模拟。结果显示，当冷风机横向吊顶时，高温库温度大部分维持在3～6 ℃，低温库温度大部分维持在-20～-18 ℃。岩芯储存库大部分区域温度分布均匀，满足储存温度需求。根据模拟结果，冷风机最终采用横向吊顶布置方式。岩芯储存库温度场数值模拟与实验结果的对比分析表明，两者温度误差在10%以内，这说明数值模拟技术可用于指导岩芯储存库工程设计。下一步可利用数值模拟技术分析当岩芯储存库内摆放货架以及改变冷风机的布置高度时对库内温度场的影响。