搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

谢文韬 李若如 彭政 蒋亦民

引用本文:
Citation:

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

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

Maximum ceasing angle of inclination andflux formula for granular orifice flow in water

Xie Wen-Tao, Li Ruo-Ru, Peng Zheng, Jiang Yi-Min
PDF
HTML
导出引用
  • 实验测量了完全浸泡在水中的玻璃珠颗粒样品在重力驱动下通过不同倾角和孔径的圆形孔洞的流量. 发现与空气中的情形类似, 不同孔径时的流量均与倾角余弦呈良好的线性关系; 线性外推得到的零流量角, 即流量休止临界角, 随颗粒粒径与孔洞直径之比d/D的减小而线性增加; 在无穷大孔径极限下, 此临界角在实验误差范围内与样品的休止角一致. 此外,所有测量都可用Beverloo公式$Q = {C_0}\rho {g^{1/2}}{(D - kd)^{5/2}}$很好地拟合;其中参数${C_0}$k仅与倾角余弦有关, 分别呈线性和平方反比关系. 与文献报道的空气孔洞流测量结果对比, 发现差别主要来自浮力和流体拖曳力对参数C0的影响. 这些结果表明用倾斜孔洞流测量颗粒材料休止角的方法和Beverloo公式具有一定的普适性. 无论颗粒间隙中填充的是水还是空气, 孔洞流的行为在定性上是一样的.
    In previous work [Acta Phys. Sin. 60 054501 (2011)], we found that for inclined Granular Orifice Flow (GOF) in air, regardless of the orifice size, the flow rate Q had a good linear relationship with the cosine of the inclination $\cos \theta $, i.e. $\dfrac{Q}{{{Q_0}}} = 1 - \dfrac{{\cos \theta - 1}}{{\cos {\theta _{\rm c}} - 1}}$, where Q0 is flow rate at $\theta ={0^ \circ }$, and ${\theta _{\rm c}}$ is the critical angle of flow ceasing obtained by linear extrapolation. Moreover, ${\theta _{\rm c}}$ increased linearly with ratio between grain and orifice diameter d/D, and at the limit of d/D going to zero (that is, D going to infinity), the angle of repose of the sample ${\theta _{\rm r}}~( = 180^ \circ - \theta _{\rm c\infty})$ was obtained. Since the flow of GOF is very stable, we believe that the linear extrapolation of the above-mentioned inclined GOF provides a novel method for accurately measuring the angle of repose of granular materials. This method has been proved to be effective in a wider orifice size range by another work [Acta Phys. Sin. 65 084502 (2016)]; and three angles, namely the repose angle measured by GOF, the free accumulation angle of a sandpile and the internal friction angle of the granular material measured by Coulomb yielding, are confirmed to be consistent. In this work, we extend this method to underwater, measuring the mass flow rate of a granular sample (glass beads) which completely immersed in water and driven by gravity, discharged from an inclined orifice for various inclination angles and orifice diameters. It is found that similar to the case in air, regardless of the orifice size, the flow rate increase linearly with the cosine of the inclination; the critical angle of flow ceasing increases linearly with ratio between grain and orifice diameter; at the limit of infinite orifice, this critical angle is consistent with the repose angle of the underwater sample within the experimental error range. In addition, all measurements can be well fitted by using the Beverloo formula $Q = {C_0}\rho {g^{1/2}}{(D - kd)^{5/2}}$, where the parameters C0 and k are only related to the cosine of the inclination, and are linear and inversely squared, respectively. Compared with the results of GOF in air reported by previous work, it is found that the difference mainly comes from the influence of buoyancy and fluid drag forces on the parameter C0. These results show that both the method of measuring angle of repose with the inclined GOF and the Beverloo formula have certain universality. The behavior of GOF is qualitatively the same whether the interstitial fluid is water or air.
      通信作者: 彭政, zpeng@csu.edu.cn
    • 基金项目: 国家级-基于人工表面等离激元的功能集成型辐射器研究(11274390)
      Corresponding author: Peng Zheng, zpeng@csu.edu.cn
    [1]

    Nedderman R M1992 Statics and Kinematics of Granular Materials (Cambridge: Cambridge University Press) pp292–294

    [2]

    陆坤权, 刘寄星 2004 物理 33 713Google Scholar

    Lu K Q, Liu J X 2004 Physics 33 713Google Scholar

    [3]

    Beverloo W A, Lenginer H A, van de Velde J 1961 Chem. Eng. Sci. 15 260Google Scholar

    [4]

    Tian Y, Lin P, Zhang S, Wang C L, Wan J F, Yang L 2015 Adv. Powder Technol. 26 1191Google Scholar

    [5]

    Lin P, Zhang S, Qi J, Xing Y M, Yang L 2015 PhysicaA 417 29Google Scholar

    [6]

    Zhang S, Lin P, Yang G, Wan J F, Tian Y, Yang L 2019 Chin. Phys. B 28 018101Google Scholar

    [7]

    Rubio-Largo S M, Janda A, Maza D, Zuriguel I, Hidalgo R C 2015 Phys. Rev. Lett. 114 238002Google Scholar

    [8]

    Janda A, Zuriguel I, Maza D 2012 Phys. Rev. Lett. 108 248001Google Scholar

    [9]

    van Zuilichem D J, van Egmond N D, DeSwart J G 1974 Powder Technol. 10 161Google Scholar

    [10]

    Peng Z, Zheng H P, Jiang Y M 2009 arXiv: 0908.0258 v3 [cond-mat.soft]

    [11]

    Madrid M A, Darias J R, Pugnaloni L A 2018 EPL 123 14004Google Scholar

    [12]

    Ji S Y, Wang S Q, Peng Z 2019 Powder Technol. 356 702Google Scholar

    [13]

    Sheldon H G, Durian D J 2010 Granul. Matter 12 579Google Scholar

    [14]

    Thomas C C, Durian D J 2015 Phys. Rev. Lett. 114 178001

    [15]

    Thomas C C, Durian D J 2013 Phys. Rev. E 87 052201Google Scholar

    [16]

    彭政, 蒋亦民 2011 物理学报 60 054501Google Scholar

    Peng Z, Jiang Y M 2011 Acta Phys. Sin. 60 054501Google Scholar

    [17]

    张昱, 韦艳芳, 彭政, 蒋亦民, 段文山, 厚美瑛 2016 物理学报 65 084502Google Scholar

    Zhang Y, Wei YF, Peng Z, Jiang YM, Duan WS, Hou MY 2016 Acta Phys. Sin. 65 084502Google Scholar

    [18]

    Wilson T J, Pfeifer C R, Mesyngier N, Durian D J 2014 Pap. Phys. 6 060009

    [19]

    Koivisto J, Durian D J 2017 Nature Comm. 8 15551Google Scholar

    [20]

    Hu W R, Zhao J F, Long M, Zhang X W, Liu Q S, Hou M Y, Kang Q, Wang Y R, Xu S H, Kong W J, Zhang H, Wang S F, Sun Y Q, Hang H Y, Huang Y P, Cai W M, Zhao Y, Dai J W, Zheng H Q, Duan E K, Wang J F 2014 Microgravity Sci. Technol. 26 159Google Scholar

    [21]

    Cheng X H, Xiao S Z, Cao A S, Hou M Y 2019 Granul. Matter 21 104Google Scholar

    [22]

    Kobayashi T, Ochiai H, Suyama Y, Aoki S, Yasufuku N, Omine K 2009 Soils Found. 49 115Google Scholar

  • 图 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.

    图 2  $D = 14$ mm, $\theta = 90^\circ $ 时典型的$M(t)$数据; 左下和右上插图分别为从主图中摘出的40−80 s及80−120 s的$M(t)$数据, 均呈良好的线性关系. 由这两段数据计算得到的流量$Q$没有差别, 均为6.53 g/s, 表明流量非常稳定

    Fig. 2.  Typical data of M(t) at D = 14 mm, $\theta = 90^\circ $;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.

    图 3  (a)不同孔径D下流量Q随倾角余弦$\cos \theta $ 的变化, 实线为直线拟合; (b)用水平($\theta = {0^\circ }$)流量${Q_0}$ 归一化的流量$Q/{Q_0}$ 随倾角余弦$\cos \theta $ 的变化关系, 实线为公式(3)的拟合结果; (c)临界流量休止角${\theta _{\rm{c}}}$随粒径-孔径比d/D的变化关系, 实线和${\theta _0}$为直线拟合结果

    Fig. 3.  (a) The variation of flow rateQ with the inclination cosine $\cos \theta $ at different orifices D, where the solid line is a linear fit; (b) variation of the normalized flow rate $Q/{Q_0}$ with $\cos \theta $, where ${Q_0}$ is the rate at $\theta = {0^ \circ }$, and the solid line is the fitted result of equation (3); (c) the relationship between the critical angle of flow ceasing ${\theta _{\rm{c}}}$ and the ratio d/D, where the solid line and ${\theta _0}$ are results of linear fitting.

    图 4  (a)用Beverloo公式(1)和(2)拟合图3数据的结果, 插图为不同倾角时${Q^{2/5}}$D的变化关系, 实线为线性拟合; (b), (c)Beverloo参数${C_0}$$k$$\cos \theta $ 的变化关系, 及其用公式(2)的拟合情况. 图(c)中的插图是${k^{ - 2}}$$\cos \theta $的变化关系

    Fig. 4.  (a) Results of fitting the data in Figure 3 using the Beverloo formula (1) and (2), the inset is the change of ${Q^{2/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 $\cos \theta $, and solid lines are results of fits by using equation (2). The inset in (c) is the change of ${k^{ - 2}}$ with $\cos \theta $.

    图 5  (a)和(b)是水中(实心方点)和空气中(空心圆点)GOF的Beverloo参数${C_0}$k$\cos \theta $的变化; (c)−(f)分别为$\theta = {0^ \circ }, {60^ \circ }, {90^ \circ }, {120^ \circ }$时, 水中和空气中GOF流量QD的变化, 插图是$Q/{C_0}$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 \theta $; (c)−(f): the changes of GOF flow rate Q with D in water and in air when $\theta = {0^ \circ }, {60^ \circ }, {90^ \circ }, {120^ \circ }$, respectively, and the inset is the change of $Q/{C_0}$ 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])

  • [1]

    Nedderman R M1992 Statics and Kinematics of Granular Materials (Cambridge: Cambridge University Press) pp292–294

    [2]

    陆坤权, 刘寄星 2004 物理 33 713Google Scholar

    Lu K Q, Liu J X 2004 Physics 33 713Google Scholar

    [3]

    Beverloo W A, Lenginer H A, van de Velde J 1961 Chem. Eng. Sci. 15 260Google Scholar

    [4]

    Tian Y, Lin P, Zhang S, Wang C L, Wan J F, Yang L 2015 Adv. Powder Technol. 26 1191Google Scholar

    [5]

    Lin P, Zhang S, Qi J, Xing Y M, Yang L 2015 PhysicaA 417 29Google Scholar

    [6]

    Zhang S, Lin P, Yang G, Wan J F, Tian Y, Yang L 2019 Chin. Phys. B 28 018101Google Scholar

    [7]

    Rubio-Largo S M, Janda A, Maza D, Zuriguel I, Hidalgo R C 2015 Phys. Rev. Lett. 114 238002Google Scholar

    [8]

    Janda A, Zuriguel I, Maza D 2012 Phys. Rev. Lett. 108 248001Google Scholar

    [9]

    van Zuilichem D J, van Egmond N D, DeSwart J G 1974 Powder Technol. 10 161Google Scholar

    [10]

    Peng Z, Zheng H P, Jiang Y M 2009 arXiv: 0908.0258 v3 [cond-mat.soft]

    [11]

    Madrid M A, Darias J R, Pugnaloni L A 2018 EPL 123 14004Google Scholar

    [12]

    Ji S Y, Wang S Q, Peng Z 2019 Powder Technol. 356 702Google Scholar

    [13]

    Sheldon H G, Durian D J 2010 Granul. Matter 12 579Google Scholar

    [14]

    Thomas C C, Durian D J 2015 Phys. Rev. Lett. 114 178001

    [15]

    Thomas C C, Durian D J 2013 Phys. Rev. E 87 052201Google Scholar

    [16]

    彭政, 蒋亦民 2011 物理学报 60 054501Google Scholar

    Peng Z, Jiang Y M 2011 Acta Phys. Sin. 60 054501Google Scholar

    [17]

    张昱, 韦艳芳, 彭政, 蒋亦民, 段文山, 厚美瑛 2016 物理学报 65 084502Google Scholar

    Zhang Y, Wei YF, Peng Z, Jiang YM, Duan WS, Hou MY 2016 Acta Phys. Sin. 65 084502Google Scholar

    [18]

    Wilson T J, Pfeifer C R, Mesyngier N, Durian D J 2014 Pap. Phys. 6 060009

    [19]

    Koivisto J, Durian D J 2017 Nature Comm. 8 15551Google Scholar

    [20]

    Hu W R, Zhao J F, Long M, Zhang X W, Liu Q S, Hou M Y, Kang Q, Wang Y R, Xu S H, Kong W J, Zhang H, Wang S F, Sun Y Q, Hang H Y, Huang Y P, Cai W M, Zhao Y, Dai J W, Zheng H Q, Duan E K, Wang J F 2014 Microgravity Sci. Technol. 26 159Google Scholar

    [21]

    Cheng X H, Xiao S Z, Cao A S, Hou M Y 2019 Granul. Matter 21 104Google Scholar

    [22]

    Kobayashi T, Ochiai H, Suyama Y, Aoki S, Yasufuku N, Omine K 2009 Soils Found. 49 115Google Scholar

计量
  • 文章访问数:  6272
  • PDF下载量:  111
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-13
  • 修回日期:  2020-03-09
  • 刊出日期:  2020-05-20

/

返回文章
返回