摘要:为了探究电沉积法制备金属泡沫微观孔胞结构对渗透特性的内在作用机理,针对不同结构与材料的金属泡沫完成了空气渗透流动实验。采用基于正五边形宽骨架结构的正十二面体的金属泡沫孔胞模型,建立了金属泡沫半经验渗透模型;分析了金属泡沫内流动阻力的变化规律,获得了包括孔隙率和孔密度等不同微观结构参数下金属泡沫内流体渗透特性。研究结果表明:小流速范围内,金属泡沫内流体流动为Darcy流态,受黏性作用影响;大流速范围内,金属泡沫内流体流动为Forhheimer流态,受惯性作用影响。电沉积金属泡沫内阻力因子的常数项为0.099 5,较粉末烧结金属泡沫内阻力因子的常数项明显增大。提出了基于孔隙率与孔密度的金属泡沫水力直径修正计算模型,与正五边形宽骨架模型预测结果相比误差在±5%以内;获得了单位长度压差的Darcy-Forhheimer型预测公式,预测值与实验值的最大误差在±15%以内,并提出了以孔隙率和水力直径为变量的金属泡沫渗透率K和惯性系数F的预测公式。
关键词:金属泡沫;渗透率;惯性系数;孔胞结构
金属泡沫凭借独特的微观结构,具备以下特点:①孔隙率高,ε大于85%;②密度低;③比表面积大,约为2 000~10 000 m2/m3;④导热系数高;⑤流道复杂;⑥孔胞结构可控;⑦机械加工性能优良。近年来金属泡沫被广泛应用于强化传热、油气水分离、催化反应、消减噪声等工程领域,特别在空气深度冷却技术中,应用金属泡沫材料的换热器可以满足飞行器轻质、紧凑、高效的工作要求[1]。
Xu等分别采用热平衡模型与非热平衡模型研究了金属泡沫内流动与传热规律,发现对于完全填充金属泡沫的流动区域,热平衡模型适用于金属泡沫孔隙率或流固介质导热系数较大的工况[2]。Chen等将金属泡沫层包覆在圆管管束外,研究了金属泡沫作为一种新型翅片形式对流动与换热特性的影响规律[3]。Yang等探索了矩形通道内金属泡沫翅片填充相变材料的储热性能,发现金属泡沫翅片可以提升储能温度、缩短相变材料的相变时长[4];管束间金属泡沫翅片填充相变材料可以缩短相变时长,达到无金属泡沫填充管束的88.548%[5]。在微通道强化传热领域,Xu等研究了多孔介质内渗透流动与微通道内边界滑移的耦合作用[6];对于微通道内部分填充金属泡沫,考虑多孔介质与微通道纯流体区域的复杂工况,Xu等获得了考虑双滑移效应的流动与传热解析解,分析了流动惯性作用、速度与热边界层滑移效应、流固非热平衡效应对热传递过程的影响机理[7-8]。
基于对矩形通道和圆管内填充金属泡沫的流动阻力实验结果分析,大量学者对不同加工制备方法得到的金属泡沫的阻力因子进行了研究[8-10],不同制备工艺下金属泡沫的渗透阻力因子随雷诺数变化总体符合反比例函数曲线,但流动阻力相差较大;相同制备工艺下的金属泡沫流动阻力也有所差别。
由于金属泡沫内部的微观几何结构比较复杂,金属骨架通过在空间中多段相连构成孔胞单元。因此,为了研究泡沫金属内部的热传导和流动现象,以上金属泡沫的应用研究中,从宏观角度将金属泡沫简化为一种大尺度下各相均匀的多孔介质。多孔介质的有效导热系数、渗透率等宏观参数通过对金属泡沫微观孔隙尺度的孔胞结构进行建模分析获得。目前金属泡沫孔胞结构的理想模型有很多种,例如平行圆柱结构[9]、正六面体结构[10-11]、正五边形组成的正十二面体金属泡沫孔胞结构[12]等规则几何体孔胞结构。简化的孔胞模型适用范围和精度受孔隙率和孔密度的影响较大,而Calmidi基于正十二面体金属泡沫孔胞结构模型获得了以金属泡沫结构参数为变量的金属泡沫渗透率和惯性系数[13],可适用于较大孔隙率范围的金属泡沫,是目前认可度比较高的孔胞结构简化模型。
金属泡沫的微观结构参数与物理特性是金属泡沫工程应用研究的基础,因此本文通过实验与数值建模分析对比金属泡沫的渗透阻力特性,同时基于微观孔胞结构模型,得到以金属泡沫孔胞结构参数为变量的金属泡沫内渗透特性参数,为金属泡沫在强化传热领域的延伸应用提供一定的参考。
1 实验系统与实验样件
1.1 实验系统
本文的实验工作是在西安交通大学动力工程多相流国家重点实验室空气深度预冷实验平台上进行的。图1为实验平台系统示意图。
1—双螺杆式空气压缩机;2—稳压罐;3—干燥器;4—过滤器;
5—质量流量计;6—空气加热器;7—PID控制器;8—调节阀;
9—压力传感器;10—热电偶;11—压差传感器;12—控制阀;
13—水箱;14—离心水泵。
图1 实验平台系统示意图
实验以空气作为工作介质。空气由螺杆式空气压缩机提供,压缩机内的稳压罐稳定系统气体压力,空气由高准质量流量计F010计量,电加热器可以控制空气的温度。通道入口布置罗斯蒙特微压力传感器和Honeywell微差压传感器,分别测量入口压力和进出口压差;通道入口和出口分别布置了3×3的T型Omega热电偶测量温度,实验段尺寸为140 mm×100 mm×18 mm,测试金属泡沫尺寸为40 mm×100 mm×18 mm,图2给出了金属泡沫示意图。实验段外壳由耐高温、低导热率PC板制成。
图2 金属泡沫示意图
空气质量流量mair的范围为6~130 kg/h;空气温度T的范围为20 ℃~110 ℃。
1.2 实验样件
对于金属泡沫微观结构,形态参数主要包括了孔径、丝径、孔隙率和孔密度。孔隙率采用排水称重法测量,其他参数采用电子显微镜测量。由于金属泡沫制备方法不同,相同形态参数下的物理特性有所不同。本文对7组不同孔胞结构参数的金属泡沫样件的流动特性进行了实验研究。实验中所使用的金属泡沫材料的物理特性见表1。
表1 金属泡沫材料物理特性
2 实验结果及讨论
本文对电沉积法制备的金属泡沫微观结构参数进行了测量,深入分析了其流体渗透特性,并建立了更符合金属泡沫孔胞结构的数学模型,得到了相关特性参数。
2.1 金属泡沫阻力特性分析
空气在金属泡沫中的流动阻力主要表现在多孔介质的黏性作用与流动的惯性作用。Darcy渗透模型为
(1)
式中:u为金属泡沫样件入口空气流速;Δp为金属泡沫样件的流动阻力;L为金属泡沫样件沿流动方向的长度;μ为空气的动力黏度;K为金属泡沫渗透率。单位长度压差Δp/L与空气流速为线性关系,与金属泡沫渗透率成反比。这种仅考虑多孔介质中流体黏性作用的模型明显与金属泡沫中的流动阻力特性不相符合。因此,将表征高速流动时的速度二次项的Forchheimer项引入Darcy方程中,即多孔介质渗透流动模型中经典的Darcy-Forchheimer方程,表达式为
(2)
式中:F为金属泡沫的无量纲惯性系数;ρ为空气密度。考虑了金属泡沫中的惯性作用,流速越快,金属泡沫中的惯性作用越显著。图3给出了单位长度下压差随空气流速的变化,从图3中可以看出,随着流速的增大,单位长度压差也明显增大。
图3 单位长度压差随空气流速的变化
同时,孔隙率与孔密度作为金属泡沫最重要的两个微观结构形态参数,对单位长度压差的影响也十分明显。对于金属泡沫2、3、4、7,孔密度均为2.54 mm-1,单位长度压差随着孔隙率的增大而减小。这主要是由于较大孔隙率的金属泡沫具有较大的孔径和较小的流体固体接触面积,因此耗损的泵功也相对较小。对于金属泡沫1、3、5、6,孔密度从1.27 mm-1逐渐增大到5.08 mm-1,单位长度压差相应急剧增大,这主要是由于孔密度的增大导致相同长度内更多的泡沫骨架和更复杂的微观结构,会引起空气流动更大的扰动和摩擦阻力。
图4给出了在较宽空气流速范围内的金属泡沫单位长度压差的二次项特性。根据Darcy方程式(1)和Darcy-Forchheimer方程式(2),不同流型下的渗透率K和惯性系数F可以由拟合曲线得到。
(a)完整流速范围
(b)小流速范围
图4 不同流速范围内的单位长度压差梯度
但是Darcy流态仅在流体流速较小的工况下存在,为了识别Darcy流态与Post-Darcy流态之间的转换区域,对式(1)(2)的左右两侧除以空气流速u,可以得到单位长度压差梯度
(3)
(4)
对于Darcy流态,式(3)的右侧为常数,速度的二次项影响可以忽略;Darcy流态的流体流动动能较小,流体附着在金属泡沫的骨架上,因此流动压降仅受到黏性力的影响,而惯性力和尾迹区域的影响几乎可以忽略。Darcy流态下的流体流动属于稳定的层流。从图4b中可以看出,当空气在小流速区域内流动时,单位长度压差梯度几乎不变。对于孔密度P为5.08 mm-1的金属泡沫,由于其流动阻力最大,在较小流速时,其单位长度压差梯度出现了下降的变化趋势,这一区域的流动属于Pre-Darcy流态。目前对于Pre-Darcy流动的流动机理和边界划分还处于初步探索的阶段。
对于Post-Darcy流态,式(4)是关于空气流速u的线性关系,速度的二次项对压差的影响十分显著。图4为不同流速范围内的单位长度压差梯度,从图4中可以看出,在空气流速大于1 m/s时,单位长度压差梯度随空气流速线性增长。因此,空气流速0.13 m/s是本文使用金属泡沫的Darcy流态向Post-Darcy流态转换的边界。
Post-Darcy流态根据流动特性的不同区分为Forchheimer流态和湍动流态。对于Forchheimer流态,流体流动的能量耗散主要包括黏性力和惯性力的综合影响。流体在金属骨架外形成了边界层。流动中运动能量的衰减主要是由于孔胞结构的收缩,即流动区域的减小和流动方向的改变,即环绕金属泡沫骨架的流动。但是Forchheimer流态下的流体流动仍然属于稳定的层流。湍动流态下的流体流动是非稳态和高度无序的,主要受到了非均匀的孔胞大小分布、接触面的粗糙程度、金属泡沫内的流体聚合与分离的综合影响。Rode等认为金属泡沫内的湍动流态属于一种表象湍流,而非充分发展的湍流[14]。Dukhan通过对单位长度压差梯度的分析,认为对于相同的金属泡沫,在不同流态下的渗透率K和惯性系数F是不同的[15]。因此应当根据不同流速工况,分别讨论金属泡沫的流动特性参数的影响。
2.2 金属泡沫内阻力因子
根据Darcy方程式(1)和Darcy-Forchheimer方程式(2),金属泡沫内的流动阻力可以无量纲化为阻力因子f,用于表征不同结构金属泡沫的流动特性。以为特征长度的阻力因子fK可以表达为
(5)
式中:ReK是以金属泡沫渗透率为特征长度的无量纲流动参数;为特征长度,用于表征金属泡沫的微观孔结构;C为常数项。
(6)
图5 阻力因子fK与ReK的关系
图5给出了阻力因子fK与ReK的关系图。从图5中可以看出,在较小空气ReK时,阻力因子fK与ReK成线性关系;随着ReK的增大,空气流动进入Forchheimer流态,空气流速的二次项起到显著作用。对于整体空气流速范围,阻力因子fK与ReK成指数变化关系,通过对实验数据的回归统计,得到本文实验所用的电化学沉积法制备金属泡沫的阻力因子fK经验关联式为
(7)
很多学者对不同加工制备方法得到的金属泡沫的阻力因子进行了研究分析,对于式(5)中的C得到了不同的结果。其中Kim对电化学沉积法制备铝泡沫研究得到的C值为0.105[16],Beavers和Sparrow得到的C值为0.074[17],Wang对粉末烧结制备不锈钢泡沫得到的C值为0.024[18]。本文常数项为0.099 5,与Kim[16]的结果十分接近,主要是由于两种金属泡沫采用相同的制备工艺,而Wang采用的粉末烧结制备法得到的金属泡沫在孔胞结构与骨架连接面都有较大的差异,导致其流动阻力特性截然不同[18]。
2.3 金属泡沫渗透特性参数分析
本文在Darcy-Forchheimer方程的基础上,基于正十二面体的金属泡沫孔胞模型,用于预测金属泡沫的渗透率K和惯性系数F。
根据Ergun对堆积床类多孔介质中阻力系数的定义[19]
(8)
式中:dh为水力直径。阻力因子f可以表达为
(9)
Re可以定义为以水力直径dh为特征长度的无量纲流动参数
(10)
根据Dietrich的定义[20]可知,水利直径dh的表达式为
(11)
式中asf为金属泡沫的比表面积。根据式(8)~(10),提出了单位长度压差的Ergun型关联式
(12)
从式(11)可以看出,单位长度压差的Ergun型并联式转换为了Darcy-Forchheimer型关联式。为了得到金属泡沫的比表面积,Huu建立了正五边形的正十二面体孔胞结构[21],根据杨肖虎等的研究发现,从电子显微镜下的图像可以得到金属泡沫骨架节点处的直径是骨架直径的1.25倍[22]。因此,本文采用正五边形宽骨架结构的正十二面体孔胞结构,如图6所示,其中df为金属泡沫圆柱骨架的直径,dp为正五边形内切圆的直径,ds为正五边形的边长。
图6 基于正五边形宽骨架结构的正十二面体模型[20]
根据正五边形宽骨架结构的几何特性,可以得到其比表面积的表达式为
(13)
式中:黄金比例φ=(51/2+1)/2;dc为金属泡沫的孔胞直径;c为骨架直径df与骨架长度ds的比值,可以由下式得到
(14)
孔胞直径dc与孔密度P的几何关系为
(15)
Huu在对金属泡沫微观结构的图像分析中发现,在孔隙率小于0.9时,金属泡沫骨架节点处为圆形,而在孔隙率大于0.9时,金属泡沫骨架节点处为三角形[21]。因此,随着孔隙率的增大,骨架直径与骨架长度的比值c和比表面积是不连续的。本文引入了金属泡沫形状参数G,它的定义为
(16)
综合上述式(12)~(16)可以得到金属泡沫孔胞单元的水力直径
(17)
针对孔隙率对金属泡沫骨架节点形状的影响,本文比较了基于正五边形宽骨架模型和Dietrich模型[20]的dh和修正水力直径dh。如图7所示,修正水力直径dh与基于正五边形宽骨架结构模型的水力直径dh的误差在±5%以内,可以更简洁的根据金属泡沫的孔隙率与孔密度两个微观结构参数确定水力直径。
图7 基于正五边形宽骨架结构模型和Dietrich模型[20]的dh和修正水力直径dh的对比结果
根据式(11)和式(17),对实验数据进行线性回归拟合,可以得到式(12)中常数A、B分别为326和1.43,即以金属泡沫孔胞单元的水力直径dh为特征长度的单位长度压差预测模型为
(18)
图8为单位长度压差的模型预测误差,从图8中可以看出,所有金属泡沫实验点的最大误差在±15%以内。基于陶瓷泡沫的Dietrich预测模型[20]的预测值与实验值相对误差也在±40%以上。本文基于正五边形宽骨架结构的正十二面体孔胞单元模型,可以更为精确地预测金属泡沫流动渗透特性。
图8 单位长度压差的模型预测误差
由式(12)可以看出,单位长度压差的预测公式为Darcy-Forchheimer型速度的二次式。由Darcy-Forchheimer方程式(2)和式(18)可以得到金属泡沫的渗透率K和惯性系数F
(19)
(20)
3 实验测量的不确定度分析
本文实验中所有数据均为间接被测量,压力、压差和流量传感器输出信号为4~20 mA电流,通过250 Ω的标准电阻电流信号转换为1~5 V的电压信号;温度信号输出为毫伏电压信号。间接被测量的不确定度受到多个直接测量的物理量误差影响,假设共有n个,且各分量彼此独立,则间接被测量的不确定度可以表示为
(21)
式中:U为间接测量物理量;δU为间接测量物理量的不确定度,δxi为各直接测量值xi的不确定度。
本文实验中,空气流量测量精度等级为0.5级,量程范围是0~100 kg/h,最大不确定度为2.41%。空气压力测量精度为0.1级,量程范围为0~6 220 Pa,根据最小测量值确定的最大不确定度为1.56%。空气压差测量精度为0.037 5级,量程范围为0~6 kPa,根据最小测量值确定的最大不确定度为0.57%。
4 结 论
本文通过实验研究的方法得到了矩形通道内金属泡沫的流动阻力特性,给出了孔隙率与孔密度对金属泡沫内流动阻力的影响规律,拟合了电沉积法制金属泡沫的阻力因子关联式,对比不同加工工艺的金属泡沫流动阻力特性,发现随着制备工艺的成熟与标准化,电沉积法制金属泡沫的阻力特性较为稳定,但不同材料的金属泡沫仍需要进行阻力因子校正。
本文分析了金属泡沫内渗透流动的流态分布与转换,由于金属泡沫高渗透率的特性,使得渗透流态主要分布在Darcy流态与Forchheimer流态,并确定了两种流态的转换边界。
本文基于正五边形宽骨架结构的正十二面体孔胞单元模型,得到了金属泡沫孔胞单元的水力直径和比表面积,根据孔胞单元模型与实验数据,获得了单位长度压差的预测公式,并将预测值与实验值进行对比,相对误差在±15%以内。与之前学者们提出的不同预测模型对比发现,本文提出的预测模型更适用于金属泡沫类型的多孔介质。在孔胞模型的基础上,得到了以孔隙率和水力直径为变量的金属泡沫渗透率K和惯性系数F。