水中颗粒孔洞流的最大休止倾角和流量公式*

谢文韬 李若如 彭政† 蒋亦民

1) (中南大学物理与电子学院, 长沙 410083)2) (中南大学粉末冶金研究院, 长沙 410083)(2020年2月13日收到; 2020年3月9日收到修改稿)

实验测量了完全浸泡在水中的玻璃珠颗粒样品在重力驱动下通过不同倾角和孔径的圆形孔洞的流量.发现与空气中的情形类似, 不同孔径时的流量均与倾角余弦呈良好的线性关系; 线性外推得到的零流量角,即流量休止临界角, 随颗粒粒径与孔洞直径之比d/D 的减小而线性增加; 在无穷大孔径极限下, 此临界角在实验误差范围内与样品的休止角一致. 此外,所有测量都可用Beverloo 公式 Q =C0ρg1/2(D −kd)5/2 很好地拟合;其中参数 C 0 和k 仅与倾角余弦有关, 分别呈线性和平方反比关系. 与文献报道的空气孔洞流测量结果对比, 发现差别主要来自浮力和流体拖曳力对参数C0 的影响. 这些结果表明用倾斜孔洞流测量颗粒材料休止角的方法和Beverloo 公式具有一定的普适性. 无论颗粒间隙中填充的是水还是空气, 孔洞流的行为在定性上是一样的.

1 引 言

人们早就知道, 颗粒孔洞流(granular orifice flow, GOF)的流量非常稳定, 用沙漏计时比用水漏计时更加精确[1,2]. 最有名的三维颗粒孔洞流量的经验公式由Beverloo 在1961年给出[1,3], 对于底部开口尺寸为D的料仓, 只要满足颗粒高度H> 2D和料仓直径F> 2D, 那么孔洞的流量就几乎只是孔洞尺寸D的函数, 而与料仓中剩余的颗粒高度和施加于顶部的荷载压力无关.Beverloo 流量表达式为

其中,ρ为颗粒材料的体密度,g为重力加速度,d是颗粒的平均粒径,C0和k是取决于颗粒材料的常数. 该经验公式有一个基于“自由下落拱”(free fall arch, FFA)假说给出的解释[1,4−7], 这个假说认为在出口附近存在一个由出口尺寸D决定的“拱”(arch), 承受住了上方的所有压力, 从而保证了在“拱”的下方, 颗粒是自由下落(free fall)的, 即认为料仓内的应力分布在“拱”的地方出现了较大的跃变甚至不连续的间断. 对于直径为D的孔洞, 颗粒从“拱”上掉落至出口的典型速度而孔洞的面积有三维孔洞流的流量至于(1)式中的kd则来自于颗粒尺寸d对孔大小D的修正. 这个图像听上去有一定道理, 所以孔洞附近的颗粒速度场和应力场的分布开始受到关注[7,8], 但对FFA 在GOF 中是否真实存在仍然有争议[7,9]. 最近的计算机模拟工作显示, FFA 即使存在也不会是传统意义上的应力不连续的拱[7].

为了从更多侧面理解GOF 的流量为什么这么稳定, 一些新的实验方法和几何装置开始出现, 例如有研究人员实验和模拟研究了外加压力对GOF 流量的影响[10−12], 以及孔洞平面的倾斜角θ对GOF 流量的影响[13−17]. 文献[13]发现, 对于大孔洞GOF 来说, 流量与 c osθ呈线性关系, 并认为可以保留Beverloo 公式(1)对倾斜孔洞情况的适用性, 只是将其中的系数C0和k理解为下面的倾角余 弦函数:

其中a1,a2,a3,a4是常数. 文献[16]基于此进一步发现, 将流量对 c osθ的线性关系外推得到的零流量休止角θc与粒径-孔径比d/D也呈很好的线性关系, 并且在趋向零(即D趋向无穷)的极限下,给出了颗粒材料的休止角θr. 这显示动态的颗粒流动和静态的颗粒休止角之间存在联系. 由于GOF 的流量很稳定, 上述倾斜GOF 的线性外推提供了一个精确测量颗粒材料休止角θr的新方法.文献[17]在更宽广的孔径范围验证了这一结论, 并证实了GOF 测得的颗粒休止角与颗粒自由堆积角以及由库伦屈服测得的颗粒材料内摩擦角三者的确具有一致性.

文献[13−17]考虑的都是空气中的颗粒样品,近来水中的GOF 开始受到关注. 文献[18,19]测量了没有倾斜(θ=0 )情况的水中GOF 流量, 发现流量会在料仓接近流完时加大(即流得更快), 并认为这是间隙流体的影响. 同时还认为对于水中的GOF, Beverloo 公式(1)不再正确, 式中流量Q随D的变化关系不再是5/2次方, 而是变为D2关系.这是由于水的粘性较大, 与空气中自由下落的颗粒不同, 水中颗粒很容易达到一个与D无关的饱和速度vt, 从而使得GOF 的流量Q=ρvts ∝D2. 为了考察利用GOF 测量休止角θr的方法对间隙流体为水的情况是否依然有效, 本文通过将文献[16]的装置完全浸入水中, 测量了完全饱和水情况下的GOF 流量随D和θ的变化关系. 结果显示, 水与空气的行为在定性上是一样的, 水中的倾斜颗粒孔洞流仍然能在D趋向无穷时逼近颗粒体的休止角.特别是与文献[18,19]不同, 测量认为Beverloo 公式(1)及其5/2 幂率仍然适用.

近年来, 对微重力环境下颗粒材料的力学行为的研究越来越受到关注[20−22]. 因水的浮力可以等效为抵消部分甚至全部重力, 所以水中实验常被用来模拟费用昂贵的太空实验. 有工作报道, 与休止角关系密切的颗粒材料内摩擦角在微重力极限下有增大的趋势[21,22]. 本文的水中GOF 休止角测量方法对此有一定的参考价值. 这里没有看到微重力实验报道的休止角增大现象, 应该是使用的颗粒材料体密度与水相差较大, 对重力的抵消不足所致.不排除利用与水密度接近的颗粒, 可观察到休止角增大现象. 但注意随着流量的减小, 测量装置的精度需要大幅度提升.

2 实验装置与实验结果

如图1(a)所示, 实验的主要结构是一个部分浸泡在水中的铝合金方形料仓, 料仓横截面为100 mm × 100 mm, 高700 mm, 壁厚3 mm. 颗粒样品采用直径d=0.9±0.1 mm 的球形玻璃珠,在实验过程中颗粒样品始终位于水面以下, 处于饱和水状态. 悬挂着料仓的刚性悬臂由两个调节过水平的电子秤支撑(精度为0.5 g, 量程为30 kg, 采样频率为10 Hz). 电子秤固定在刚性支架上, 示数之和就是料仓的总质量M(t) (包含了不随时间变化的仓壁自重和随时间变化的颗粒质量). 实验中随着颗粒的流出, 电子秤示数M(t) 随时间减少, 减少率即为GOF 的流量Q.

当颗粒从料仓中流出时, 料仓内部减少的颗粒体积需要及时由水填补, 以维持料仓内外水面高度不变和样品的水饱和状态, 从而防止内外的液压差对颗粒孔洞流量的影响. 为此料仓的四个侧壁均镶嵌了密排冲孔(孔径为0.8 mm)的不锈钢板(图1(b))用以透水, 以保证实验过程中料仓内外液面保持高度一致(也保证了料仓壁在实验过程中受到水的浮力不变).

与文献[16]类似, 实验测量了不同倾角θ下,7个孔径D=4, 6, 8, 10, 14, 18, 20 mm 的M(t) 结果. 倾角θ定义为圆孔所在平面与外水平面的夹角,如图1(c)和图1(d)所示. 倾角大小由精度为0.1 度的倾角仪(bevel box)测出. 利用底部开孔的装置(图1(c))可测量θ∈[0◦,45◦] 的颗粒流量. 对于θ>45◦的情形, 需要采用图1(d)所示的料仓侧壁开孔装置. 在距离底部约100 mm 的高度开一个100 mm × 100 mm 的方孔, 用以嵌入圆孔D所在的硬铝板. 为尽量减小板厚度对流量的影响, 孔洞外侧加工成了如图1(e)所示的楔形.

图1 (a)实验装置示意图; (b)料仓透水侧壁照片; (c)倾角小于45 度时采用的实验装置; (d)倾角大于45 度时采用的实验装置;(e)楔形孔洞示意图Fig. 1. (a) Schematic of the setup; (b) photograph of the permeable side wall of the silo; (c) the experimental devices used when the inclination is less than 45 degrees; (d) the experimental devices used when the inclination is greater than 45 degrees; (e) schematic of the wedge-shaped orifice D.

实验中发现, 如果将颗粒直接从空气置入水中的料仓, 颗粒表面会裹挟大量空气泡, 这些气泡在流动时会逐渐释出, 对孔洞流量带来显著影响. 在孔洞较小时气泡甚至会直接阻塞颗粒流. 因此在开始实验之前, 需要提前将颗粒放入水桶中, 让其在水中反复多次流动, 使气泡充分释出直至肉眼不可见. 之后再将去气泡后的颗粒通过开口的侧壁在液面以下加入料仓, 并且在加料的过程中始终保证颗粒不再暴露到空气中. 通过去气泡处理后的颗粒流量将十分稳定和可重复.

实验开始前先将孔D塞住, 在料仓中加入去气泡处理过的颗粒直至接近水面, 然后开放孔D, 同时采集电子秤随时间变化的M(t) . 典型的M(t) 数据如图2 所示. 由图2 可以看出, 在初期GOF 有一段短暂的不稳定流动期(包括开放孔D带来水面波动的影响), 但之后M(t) 呈现出良好的直线关系, 其斜率的负值就是GOF 的流量Q(由于测量的是料仓中颗粒的减少量, 故M(t) 斜率为负). 这个流量在很长一段时间内是稳定不变的, 如图2 插图所示, 40—80 s 和80—120 s 计算出来的流量Q均为6.53 g/s. 在流动的尾声(140 s 左右),M(t) 曲线只有非常微弱(小于2%)的变陡趋势. 没有明显地观察到文献[19]中提到的在料仓快流完时流量变大的现象. 这个分歧的原因目前不详. 但快流完的情况不是关注的重点, 本文仅研究中间段的稳定流量Q随D和θ的变化关系.

图2 D=14mm,θ=90◦时典型的 M (t) 数据;左下和右上插图分别为从主图中摘出的40−80s 及80−120s的 M (t) 数据, 均呈良好的线性关系.由这两段数据计算得到的流量 Q 没有差别, 均为6.53 g/s,表明流量非常稳定Fig. 2. Typical data of M(t) at D = 14 mm, ;the lower left and upper right insets are the data of 40−80 s and 80−120 s extracted from the main graph, both of which show good linearity. Both flow rates Q calculated from these two insets are 6.53 g/s, indicating that the flow is very stable.

对于每一给定的D和θ重复三次实验, 平均值记为对应的流量Q(D,θ) , 结果见图3(a)所示. 虽然水中操作的扰动比空气中的实验略大, 但在去除气泡后水中的GOF 流量仍然非常稳定, 实验重复性好, 三次测量的偏差均小于数据点图标的大小,因此图3a 的实验数据点没有标出误差棒.

与文献[16]的空气实验结论一样, 对不同孔径D, 流量Q均随着角度θ的增大而逐渐减小, 并与cosθ呈良好的线性关系(图3(a)). 注意对较小的D(如D= 4 mm)而言, 因GOF 在倾角较小时流动就停止了, 所以测量时 c osθ的取值范围也相应变 小.

图3 (a)不同孔径D 下流量Q 随倾角余弦 c os θ 的变化,实线为直线拟合; (b)用水平( θ =0◦ )流量 Q 0 归一化的流量 Q /Q0 随倾角余弦 c os θ 的变化关系, 实线为公 式(3)的拟合结果; (c)临界流量休止角 θc 随粒径-孔径比d/D 的变化关系, 实线和 θ0 为直线拟合结果Fig. 3. (a) The variation of flow rateQ with the inclination cosine c os θ at different orifices D, where the solid line is a linear fit; (b) variation of the normalized flow rate Q/Q0 with c os θ , where Q 0 is the rate at θ =0◦ , and the solid line is the fitted result of equation (3); (c) the relationship between the critical angle of flow ceasing θc and the ratio d/D, where the solid line and θ0 are results of linear fitting.

3 流量与休止角

按照与文献[16]类似的分析和处理方式, 用θ=0◦的水平GOF流量Q0对Q进行归一化, 得到图3(b)所示的Q/Q0与 c osθ的线性关系. 这时所有的直线均过坐标点(1, 1), 只有斜率一个参量, 所以有

其中,θc表示不同D下流量Q=0时的倾角, 是孔径D的函数. 与文献[16]的结论类似,θc与d/D的函数关系仍然是线性的(图3(c)). 用下面公式:

拟合, 得到D→∞时的零流量倾角θ0=156◦±4◦和的补角θr=180◦−θ0=24◦±4◦对应无限大颗粒自由平面极限下的休止角[16,17]. 用同样颗粒材料在水中堆积锥形沙堆, 测得水中沙堆休止角约在 23◦至 25◦之间, 与(4)式外推得到的结论相符. 这些结果表明, 无论在空气中还是水中,都可以用倾斜GOF 的方法测量颗粒材料休止角.注意由于GOF 在水中流动时受间隙流体的影响比空气中要大, 所以测量值的涨落要略大于空气中的实验.

在文献[16]中用同样颗粒材料测得的空气中的休止角为 27◦±3◦. 虽然略高于本文的水中结果,但在误差范围内还是吻合的. 水中与空气中的休止角是否存在微弱差别, 还有待仔细测量. 定性上看,空气中的休止角多少会受湿度导致的颗粒间湿粘力(cohesive force)影响而变大. 这个现象在文献[18]中也被观察和提到. 另外如引言所述, 由于水中的浮力远大于空气中的, 水中颗粒样品的应力明显小于空气中的应力, 图1 装置可用于研究航天领域关心的微重力下颗粒休止角变大的现象. 本文未观测到休止角变大现象, 原因很可能是采用的玻璃珠和水的密度比(约为5:2)不够接近1 导致的. 今后值得考虑用更接近水密度的颗粒材料来研究和逼近这个微重力现象.

4 Beverloo 公式

图3 的实验数据可以用公式(1)和(2)很好地拟合. 如图4(a)所示, 对所有不同的倾角θ, 流量Q随D的变化关系都能很好地满足5/2次幂, 实验结果并不支持文献[18,19]提到的Q与D2成比例的结论. 这还可以在图4(a)的插图中, 从Q2/5和D的直线关系中得到验证.

由于幂次2/5 和1/2 的差别不大, 如果尝试拟合Q1/2与D的直线关系, 结果也显得不错(这可能也是文献[18,19]认为其合理的原因). 尽管两种幂次公式对本文数据的拟合相关系数(Adj. R-Square)均大于0.99(相关系数越接近1 代表拟合程度越好), 2/5次幂的公式在拟合度上仍然要略高一点.因此更倾向支持Beverloo 公式(1)对于各种倾角的水中GOF 均成立, 其中参数C0和k随倾角余弦cosθ的变化关系满足公式(2)(图4(b)和图4(c)).

图4 (a)用Beverloo 公式(1)和(2)拟合图3 数据的结果,插图为不同倾角时 Q 2/5 随D 的变化关系, 实线为线性拟合; (b), (c)Beverloo 参数 C 0 和 k 随 cos θ 的变化关系, 及其用公式(2)的拟合情况. 图(c)中的插图是 k −2 随cos θ 的变化关系Fig. 4. (a) Results of fitting the data in Figure 3 using the Beverloo formula (1) and (2), the inset is the change of Q2/5 with D at different inclination, and the solid line is a linear fit; (b) and (c) variations of the parameters C 0 and k with c os θ , and solid lines are results of fits by using equation (2). The inset in (c) is the change of k −2 with cos θ.

虽然公式(2)对于水中和空气中的GOF 都有效, 但在同样条件下水中的GOF 流量明显小于空气中的流量(图5(c)—(f)). 这个流量差别, 主要由参数C0来描述(图5(a)). 参数k随 cosθ的变化情况对水和空气不敏感(图5(b)). 另外从图5(c)—(f)的插图也可以看到,k对流量的影响也基本与水还是空气无关, 只在较大倾角时才出现一点差别. 在Beverloo 公式中,k与孔径变小导致的流量阻塞停止(clogging)有关. 同样倾角下, 水中GOF 与空气中GOF 有着几乎同样的k, 显示这个阻塞机制与间隙流体的种类关系不大.

从图6(a)可以看到, 水和空气的Beverloo 系数C0的比值基本上是一个与倾角无关的常数0.37.另外流量比值与0.37 的偏差会随着孔径D的减小和倾角的增加而变大图6(b). 显然这个偏差反映的是Beverloo 公式中系数k的效果. 值得指出的是,间隙流体对颗粒的作用有浮力和粘滞拖曳力. 前者的效果可通过将重力加速度乘以密度变化比(ρ −ρ流体)/ρ来描述(ρ=ρ颗粒是颗粒材料体密度).对Beverloo 公式(1)中的g做此处理, 得出水中与空气中GOF 的C0比值为0.78(忽略空气密度), 明显大于图6(b)的0.37. 这意味着水的拖曳力(即粘滞力)对GOF 流量有显著影响.

因实验数据很少, 水的粘滞系数η对C0的影响还不能给出可靠的经验公式. 受流体中匀速下落的单个颗粒的末速度反比于η的启发, 本文建议采用下 面的三参数模型

来拟合C0与η的关系. 其中C0vac为真空的C0值,η0和幂率a可理解为反映密集多颗粒效应的参数.当a趋向1,η0趋向0时, 它退化为单个颗粒的情况. 由于目前只有水和空气的两个C0数据0.37 和0.78, 还不能确定模型的参数及其适用性. 如果取幂率a=1 , 水和空气的粘滞系数分别为1.0×10−3Pa·S与1.8 × 10–5Pa·S, 得到参数为C0vac=0.796,η0= 0.868 × 10–3Pa·S.

图5 (a)和(b)是水中(实心方点)和空气中(空心圆点)GOF 的Beverloo 参数 C 0 和k 随cos θ的变化;(c)−(f)分别为θ =0◦,60◦,90◦,120◦时, 水中和空气中GOF 流量Q 随D 的变化, 插图是 Q /C0 随D 的变化. 空气中的实验数据来自文献[16]Fig. 5. (a) and (b): Beverloo parameters C 0 and k of GOF in water (solid squares) and in air (hollow circles) as a function of cos θ ;(c)−(f): the changes of GOF flow rate Q with D in water and in air when θ =0◦,60◦,90◦,120◦ , respectively, and the inset is the change of Q /C0 with D. The experimental data in air comes from ref. [16].

图6 水和文献[16]空气的(a) Beverloo 系数 C 0 和(b) GOF流量 Q 的比值Fig. 6. Ratio of (a) Beverloo coefficient C 0 and (b) GOF flow rate Q in water and in air (from Ref. [16]).

5 结 论

本文测量了不同倾角θ和孔径D下, 完全浸泡于水中的颗粒物质(玻璃珠)在重力驱动下从圆形孔洞卸载时的GOF 流量. 得到的结论如下:

1)与文献[16,17]报道的空气情形类似, 对所有孔径值, 流量Q都与倾角余弦 c osθ呈良好的线性关系. 由此线性外推给出的流量休止临界角θc随粒径-孔径比d/D的减小也呈良好的线性增加关系,并且在无穷大孔径极限d/D →0 下的线性外推值θ0与颗粒材料休止角相符. 这表明用倾斜孔洞流的方法同样可以测量材料在水中的休止角.

2)得到的水中休止角略微小于文献[16]报道的用同样方法测量的空气休止角. 这可能是因空气湿度导致的颗粒间粘力所致.

3)所有测量数据都可以用5/2 幂率的Beverloo公式很好地拟合, 并且与空气情况一样, 系数C0和k随孔洞倾角余弦cosθ的变化分别符合线性和平方反比关系.

4)通过与空气GOF 的测量结果[16]比较, 发现间隙流体种类(水或空气)对Beverloo 公式系数k的影响不大, 但对系数C0有大约0.37 的固定比例影响. 后者不能完全用浮力导致的重力加速度g减小来解释, 意味着水施加于颗粒的拖曳力对GOF 有显著影响.

5)鉴于水和空气都属于牛顿流体, 二者的GOF 行为在定性上应该保持一致, 例如Beverloo公式的有效性和给出的休止角. 故不倾向[18,19]中提出的水中GOF 的流量与D2成比例的提议.从本文的测量结果看不出有修改Beverloo 公式幂率的必要.