内角钝度对微重力下液体推进剂毛细流动特性的影响

杨恩博, 金宇鹏, 杨 光, 黄永华, 王天祥, 雷 刚, 吴静怡

(1.上海交通大学 机械与动力工程学院,上海 200240;2.航天低温推进剂技术国家重点实验室,北京 100028)

随着载人航天技术的发展以及空间探索的日趋长期化,空间推进剂的在轨管理成为一个重要课题,从而带动微重力流体科学的进一步发展.在微重力环境下,重力的影响可以忽略,此时表/界面张力成为主导流体行为的最主要因素.板式表面张力贮箱正是根据表面张力驱动下的内角自流现象而设计的推进剂空间管理装置.

作为微重力流体力学下的一个重要模型,内角流动是研究在表面张力主导下,液体沿固体二面角爬升的理论.有关内角流动的研究可以追溯到20世纪60年代,Concus等[1]提出微重力条件下内角流动液体前缘稳定性的临界条件,即Concus-Finn条件; Weislogel 等[2]对内角流动的Navier-Stokes方程进行简化,将三维问题简化为一维问题,利用滑移假设进行求解,提出流阻的理论近似解,并推广到复杂几何形状的计算[3-4]; Wang等[5]研究了微重力条件下不同初始液体体积对内角毛细流动的影响; 李京浩等[6]针对不对称内角情形,给出扇形内角情形下的计算公式; 沈逸等[7]利用磁补偿原理在地面实现微重力环境,并分析了重力水平、内角材质等因素对液面位置的影响.

尽管有关内角流动的理论与实验分析已经比较成熟,但现有理论模型大多仅适用于理想情况下的完全尖锐内角.而实际加工制造得到的内角通常不可避免地存在一定钝度; 此外,内角尖端采用圆角过渡也能够增强装置的机械承载能力.理想内角的流动模型无法直接推广到含钝度的毛细流动.有关内角钝度对液体内角流动的影响,Ransohoff等[8-9]对内角流动过程的流阻进行了数值计算,并研究了有关内角圆率对内角自流现象影响的一系列问题; Chen等[10]对内角存在一定钝度的情况进行分析,并研究了极限情况下薄层流动的问题; Zhou等[11]得出薄层流动的计算方法,并将其推广到一般情形; 魏月兴等[12]利用Ransohoff提出的方法对存在钝度的情形进行了理论分析.然而,由于微重力环境实验复杂且成本高,钝度对内角毛细流动的影响机理仍未完全解析并缺少实验验证.

建立含钝度的微重力内角毛细驱动流动模型,并首次开展基于磁补偿的微重力模拟实验以验证理论模型的准确性,定量获得钝度对毛细流动特性的影响特性.同时,将理论模型推广到以液氢和液氧为代表的低温推进剂的空间应用,为面向未来深空探测任务的推进剂管理装置的设计提供重要基础数据.

1 钝度内角中的毛细力驱动流动模型

图1 钝度内角中的毛细力驱动流动示意图Fig.1 Schematics of capillary flowing in rounded corners

内角毛细流动模型中,假定流体的流动方向为x轴方向,内角开口为2α,流体与壁面的接触角为θ,固体夹角处形成的圆弧曲率半径为r0,流体的毛细流动距离为xf,液体润湿内角的边长为D.液面在x处沿y-z平面的曲率半径记为r(x,t),其中x=0处的曲率半径记为R,t为时间,弯曲液面的圆心角为2δ,且在该截面上有δ=π/2-θ-α.实验表明,在内角流体爬升过程中,R始终为定值[2],从而定义τ=r0/R,表征决定钝度大小的相对曲率半径.微重力条件下,针对图1所示的流动过程,液体主要受到表面张力、外界压力以及流动阻力的共同作用.在气液交界面处,由Young-Laplace方程可得气液交界面的压强差为

(1)

式中:σ为液体的表面张力系数;r1和r2分别为y-z平面和x-y平面的曲率半径.在流动假设中,认为流动的长度远远大于截面尺度,此时r2趋于无穷大,因此只考虑r1对流动的影响.

当液体从一端进入内角时,在表面张力作用下沿x方向的曲率半径逐渐减小,从而在液体内部形成压强梯度.在流动的任意位置x,气液交界面的压强差[12]可以表示为

(2)

根据Weislogel等[2]对内角流动模型的简化,连续性方程有如下微分形式:

(3)

式中:ρ为液体密度;v为流动速度;液体在某x处截面的截面积设为S.

假设液体的密度为常数,可得

(4)

式中:q为液体的体积流量.

在图1中任意位置x处取y-z平面的横截面,可以得到该处曲率半径r(x,t)和S的关系为

此外，我国疆土东西跨度大，南北迥异，地质条件复杂多变，且不同的地质构造单元中地壳物质组成差异较大，导致不同的城市地下管线铺设方法差异较大，铺设深度和管线材质选用方面，都使得在探测过程中应该根据管线材质及用途来选择(表1)。因此，在不同的城市探测地下管线时，应结合当地地下管线的材质的探测技术，才能取得较好探测结果。

(5)

(6)

式中:μ为流体的动力黏度.

体积流量可以写成如下的微分形式:

(7)

将式 (7) 代入式(4)可得在二维流动假设条件下关于S的控制方程:

(8)

求解式(8)可得[13]:

(9)

(10)

并由此得到流量的表达式:

(11)

其中,根据薄层流动理论可得βv的取值[11, 15]如下:

(12)

φ1=cos2(α+θ)+sin(α+θ)cos(α+θ)tanα

将式(12)代入式(10),即可得到钝度内角中爬升距离以及流量的动态表达式:

(13)

(14)

当钝度τ=0时,可得到理想尖锐内角下毛细驱动距离以及流量的动态表达式:

(15)

(16)

2 基于磁补偿原理的微重力地面模拟实验

在前期已建立的磁流体微重力模拟实验台[7]中开展实验,其原理是利用磁流体作为实验工质,利用梯度磁场产生的磁场力抵消重力从而获得等效微重力环境.实验中利用亥姆霍兹-麦克斯韦线圈作为梯度磁场发生装置,可实现优于±2.5%的磁场梯度纵向非均匀度,具有成本低、重复性好、重力水平可控等显著优点.

图2为实验装置示意图.利用高精度三轴定位装置对样品位置进行调节定位,从而将样品腔置于线圈的中心位置,并保证整个样品处于磁流体的重力被磁场力完全补偿的区域内.利用高速照相机捕捉磁流体在内角整个流动过程中液面位置的动态变化,并利用光源进行辅助照明.

图2 基于磁补偿原理的微重力地面模拟实验系统Fig.2 Microgravity simulation experimental system with magnetic fluid

内角样品材料为光学透明的石英玻璃,肋板直径d= 10.0 mm,D=7.80 mm,2α=45°, 并设置不同的内角钝度值.实验所采用磁流体的物理性质如下:σ=36.2×10-3N/m,μ= 29.0×10-3Pa·s,ρ=1 163.33 kg/m3,其在石英玻璃上的接触角测量值为38° (见图3).

图3 实验样品及表面接触角Fig.3 Experimental sample and contact angle

初始状态下,内角样品在磁流体中的浸没高度为3 mm.待液面稳定后,开启线圈的电流开关.利用高速照相机记录磁流体的液面高度变化,待磁流体达到最大爬升距离且液面稳定后关闭线圈电流开关,完成一次实验记录.所得图像的位置分辨率为0.2 mm/像素,照相机的采样帧率为870帧/s.因此,实验数据中的液面位置测量误差为 ±0.2 mm,时间误差为1.1 ms.实验过程中,内角样品中液面的升高会导致样品腔中液位下降,根据体积守恒可计算得其对毛细爬升距离测量结果的影响小于4%.每个工况的实验至少独立测试3次,以保证数据结果的可重复性.

3 结果与讨论

3.1 内角流动实验及其与计算结果的对比

图4为在内角钝度τ=0.26恒定的条件下,磁流体在表面张力驱动下沿内角爬升的位置动态变化情况.其中t=0时刻为在磁场开启之前,液体在地面常重力时的状态.当开启磁场后,液面位置随时间逐渐升高.通过改变样品内角钝度的大小, 即可通过实验定量获得钝度对动态内角毛细流动特性的影响.

图4 不同时刻下磁流体的内角爬升图像(τ=0.26)Fig.4 Movement of magnetic fluid along interior corner at different moments (τ=0.26)

为进一步验证基于磁补偿原理的微重力模拟实验的准确性,首先开展存在微小钝度τ= 0.13时的内角流动实验,并与理想内角流动模型式(15)和钝度内角流动模型式(13)进行对比.由图5可见,实验过程中磁流体的毛细驱动距离xf与时间的1/2次方(t1/2)近似呈正比例关系,该结论与文献[3,12,16]中结果一致.此外, 图5中实验结果、式(13)和式(15)之间的平均相对误差在15%以内,误差线为3次独立重复实验的标准差,下同.

图5 τ = 0.13时磁流体的毛细爬升曲线Fig.5 Movement of magnetic fluid at τ =0.13

图6为在不同钝度条件下液体爬升距离动态变化的实验与计算结果.由图可见,当钝度一定时,xf与t1/2始终近似符合线性关系,且钝度值越小线性度越高.在本文所研究的钝度范围内,液体在内角中的爬升速度随着钝度值的增大而减小,该规律与Zhou等[11]的分析结果一致.对比实验结果与计算模型可知,两者的偏差随着钝度的增加而增大,其原因可能在于实验加工的圆角仍不可避免地存在粗糙度的影响,因此会导致流速降低.然而,多种不同钝度下实验与计算的平均相对偏差仍在20%以内.因此考虑到实验误差,可以认为式(13)能够定量反映在不同钝度条件下液体毛细爬升距离的动态变化特性.

图6 不同钝度下的内角爬升曲线Fig.6 Fluid motion for different roundedness of interior corners

3.2 钝度对内角毛细流动流量的影响

从前面的分析得出,在理想内角以及含钝度内角的毛细流动中,xf与t1/2始终保持线性关系.为了定量描述液体的毛细流动速度,定义毛细爬升系数为

k=xf/t1/2

(17)

图7为在不同钝度条件下毛细爬升系数的实验值与计算值.从图中可以看出,实验与计算结果均表明毛细爬升速度随钝度的增加而单调减小.当内角钝度从0.13增加到0.33时,实验测得的毛细爬升系数从22.68 mm/s1/2降低至12.69 mm/s1/2,下降约44%,而该条件下通过计算模型获得的毛细爬升系数从18.85 mm/s1/2降低至14.68 mm/s1/2.

图7 不同钝度条件下的毛细爬升系数Fig.7 Values of k of different corner roundednesses

图8为在不同钝度条件下内角毛细流动流量的动态变化计算结果.在钝度一定时,液体的流量在初始阶段随时间推移迅速减小,后又逐渐趋于稳定.其原因在于随着液体前端位置不断向前发展,液体流动的截面积逐渐减小,且爬升速度在初始阶段迅速下降,两方面因素共同加剧了流量随时间的下降.此外,在同一时刻,钝度越大,液体流量越小.

图8 不同钝度条件下内角毛细流动流量的动态变化Fig.8 Variation of capillary flow rate at different corner roundednesses

3.3 低温推进剂的内角流动特性

相较于常规推进剂,以液氢/液氧为代表的低温推进剂具有高比冲、无毒、无污染等诸多优势,是目前以及未来很长一段时间内空间工程的首选推进剂[17].然而,低温推进剂沸点低、汽化潜热低等特殊物性,为其长期空间贮存和在轨管理技术带来巨大挑战[18].因此,对于低温推进剂,采用表面张力式流体液体管理装置是解决其空间应用难题的关键突破点.基于前文的毛细流动模型,对液氢、液氧两种流体的内角流动情况进行计算分析.液氢和液氧的物理性质参数如表1所示.为了与空间贮箱中的应用一致,采用低温流体在不锈钢表面上的接触角(近似为0°)[19-22].

表1 液氢与液氧的物性参数(0.1 MPa)Tab.1 Physical properties of liquid hydrogen and liquid oxygen (0.1 MPa)

图9为在τ= 0.1,R=6×10-3m,α=15° 的内角条件下,液氢、液氧以及磁流体的毛细爬升距离随时间的变化曲线.由图可见,液氢与液氧的流动速度均显著高于前文实验中的磁流体,且液氢的运动速度高于液氧.其原因在于尽管低温推进剂的表面张力较小,但由于其黏滞系数同样较小,导致式(13)中因子σ/μ反而增大; 且上述两种流体对不锈钢有很好的润湿效果,最终导致爬升能力的提高.此外,对于上述两种低温推进剂,液氢的表面张力以及黏度均小于液氧,但由于黏度对流动的作用更显著,所以液氢的运动速度始终高于液氧.图10为在与图9相同的工作条件下, 3种不同流体流量动态变化规律的计算结果.在同一时刻, 液氢的内角流动流量高于液氧约40% 以上, 而两者流量均高出磁流体流量一个数量级.

图9 不同流体的动态毛细爬升特性对比Fig.9 Movement of different fluids along the corner at microgravity

图10 不同流体的动态流量大小对比Fig.10 Comparison of flow rate of different fluids

图11(a)和11(b)分别为钝度对液氢和液氧内角毛细流动流量的影响.随着钝度的增大,两种推进剂的流量均减小.当钝度τ从0.1增大到0.5时,液氢和液氧流量均分别下降约 17.3%,由此可见减小内角钝度对提高推进剂流量的重要性.

图11 液氢和液氧流量随钝度的变化曲线Fig.11 Flow rates of liquid hydrogen and liquid oxygen versus roundedness of interior corners

4 结论

针对含钝度条件下内角毛细流动建立数学模型,通过引入流阻的概念,得出内角毛细流动液面前端位置和液体流量的计算表达式.搭建磁补偿微重力模拟装置,在地面上开展了基于磁流体的微重力流动实验,并将实验结果与数学模型进行比较.并在此基础上,针对低温推进剂的内角流动特性进行分析,主要结论如下:

(1) 获得微重力条件下流体沿内角爬升距离与流量的计算模型,通过与微重力模拟实验数据对比,验证了该计算模型的可行性,且平均相对误差在20%以内.理论与实验数据均表明,在内角毛细流动中,xf与t1/2始终近似保持线性关系.

(2) 在内角钝度τ∈[0.13, 0.33]范围内,流体沿内角爬升的速度、流量均随着内角钝度的增大而单调减小.实验结果表明,当内角钝度从0.13增加到0.33时,磁流体的毛细爬升系数(h/t1/2)下降约44%.在相同条件下,流量随着液面高度的增加而减小.

(3) 液氢、液氧等低温推进剂由于黏度小,其在内角毛细流动流量上高于磁流体一个数量级以上.尽管液氢的表面张力小,但其低黏度的特性在微重力流动中起主导性作用,导致其内角流动速度高于液氧和磁流体.

由于微重力流动机理的复杂性,所以计算模型的精度仍需要通过空间搭载等实验途径进一步验证.此外,低温推进剂具有低沸点等特殊物理性质,其在流动过程中的热学耦合特性还需进一步探索.