搜索

x

留言板

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

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

纳米通道粗糙内壁对流体流动行为的影响

梅涛 陈占秀 杨历 王坤 苗瑞灿

引用本文:
Citation:

纳米通道粗糙内壁对流体流动行为的影响

梅涛, 陈占秀, 杨历, 王坤, 苗瑞灿

Effect of rough inner wall of nanochannel on fluid flow behavior

Mei Tao, Chen Zhan-Xiu, Yang Li, Wang Kun, Miao Rui-Can
PDF
HTML
导出引用
  • 纳米流动系统具有高效、经济等优势, 在众多领域具有广泛的应用前景. 因该类系统具有极高的表面积体积比, 致使界面滑移效应对流动具有显著影响. 本文采用分子动力学方法以两无限大平行非对称壁面组成的Poiseuille流动为对象, 分析了壁面粗糙度与润湿性变化对通道内流体流动的影响. 对于不同结构类型的壁面, 需要通过水动力位置来确定固液界面位置, 准确计算固液界面位置有助于更好地分析界面滑移效应. 研究结果表明, 上下壁面不对称会引起通道内流场参数分布的不对称, 壁面粗糙度及润湿性的变化会影响近壁面附近流体原子的流动特性, 由于壁面凹槽的存在, 粗糙壁面附近的数密度分布低于光滑壁面一侧. 壁面粗糙度及润湿性的变化会影响固液界面位置, 肋高变化及壁面润湿性对通道中速度分布影响较大, 界面滑移速度及滑移长度随肋高和润湿性的增大而减小; 肋间距变化对通道内流体流动影响较小, 界面滑移速度和滑移长度基本保持恒定.
    Flow system on a nano scale, as an effective and economic system, has been widely employed. While on a macroscopic scale, for the non-slip boundary, the velocity of the fluid at the surface is assumed to be consistent with the surface. This approximation may become invalid on a smaller length scale pertinent to the operation of microfluid devices. The interface slip effect has a significant influence on the flow because of its higher ratio of surface to volume. In this paper, the Poiseuille flow, which is composed of two infinite parallel asymmetric walls, is studied by the molecular dynamics method. The influence of wall roughness and surface wettability of channel on fluid flow in the channel are analyzed. The results show that the asymmetric upper and lower wall can lead to an asymmetric distribution of flow parameters. The change of wall roughness and wettability would affect the flow characteristics of fluid atoms near the wall. Due to the influence of wall grooves, the number density distribution near the rough wall is lower than that on the smooth wall side. As the rib height and wall wettability increase, the number density of fluid atoms in the groove increases gradually, and the change of the rib spacing does not substantially affect the number density distribution of fluid atoms near the rough wall. For different structure types of walls, the real solid-liquid boundary positions are determined by simulating the velocity field distribution in the channel under both Couette flow and Poiseuille flow, which can help us to better analyze the interface slip effect. The variation of wall roughness and wettability can affect the position of the solid-liquid interface. The change of rib height and wettability can greatly influence the velocity distribution in channel, and the position of the solid-liquid boundary as well. Conversely, the rib spacing has a less effect on the boundary position. The difference in boundary position can affect the interface slip effect. We can find the slip velocity and the slip length on one side of the rough wall to be smaller than those on the smooth wall side, and as the rib height and wall wettability increase, the slip velocity and the slip length significantly decrease near the rough wall side. The effect of rib spacing on fluid flow is trivial, and the interface slip velocity and length are relatively stable.
      通信作者: 陈占秀, zhanxiu_chen@hebut.edu.cn
    • 基金项目: 国家重点基础研究发展计划(批准号: 2018YFB0605101-1)资助的课题.
      Corresponding author: Chen Zhan-Xiu, zhanxiu_chen@hebut.edu.cn
    • Funds: Project supported by the National Basic Research Program of China (Grant No. 2018YFB0605101-1).
    [1]

    黄桥高, 潘光, 宋保维 2014 物理学报 63 054701Google Scholar

    Huang Q G, Pan G, Song B W 2014 Acta Phys. Sin. 63 054701Google Scholar

    [2]

    葛宋, 陈民 2013 工程热物理学报 34 1527

    Ge S, Cheng M 2013 J. Eng. Therm. 34 1527

    [3]

    顾骁坤, 陈民 2010 工程热物理学报 31 1724

    Gu X K, Cheng M 2010 J. Eng. Therm. 31 1724

    [4]

    Sun M, Ebner C 1992 Phys. Rev. Lett. 69 3491Google Scholar

    [5]

    Voronov R S, Papavassiliou D V, Lee L L 2007 Chem. Phys. Lett. 441 273Google Scholar

    [6]

    胡海豹, 鲍路瑶, 黄苏和 2013 力学学报 45 507Google Scholar

    Hu H B, Bao L Y, Huang S H 2013 Chin. J. Theor. Appl. Mech. 45 507Google Scholar

    [7]

    Nagayama G, Cheng P 2004 Int. J. Heat Mass Transfer 47 501Google Scholar

    [8]

    Barisik M, Beskok A 2011 Microfluid. Nanofluid. 11 269Google Scholar

    [9]

    Shi Z, Barisik M, Beskok A 2012 Int. J. Therm. Sci. 59 29Google Scholar

    [10]

    张冉, 谢文佳, 常青 2018 物理学报 67 084701Google Scholar

    Zhang R, Xie W J, Chang Q 2018 Acta Phys. Sin. 67 084701Google Scholar

    [11]

    Cieplak M 2001 Phys. Rev. Lett. 86 803Google Scholar

    [12]

    闫寒, 张文明, 胡开明, 刘岩, 孟光 2013 物理学报 62 174701Google Scholar

    Yan H, Zhang W M, Hu K M, Liu Y, Meng G 2013 Acta Phys. Sin. 62 174701Google Scholar

    [13]

    张程宾, 许兆林, 陈永平 2014 物理学报 63 214706Google Scholar

    Zhang C B, Xu Z L, Chen Y P 2014 Acta Phys. Sin. 63 214706Google Scholar

    [14]

    Rahmatipour H, Azimian A R, Atlaschian O 2017 Phys. A 465 159Google Scholar

    [15]

    Toghraie D, Mokhtari M, Afrand M 2016 Physica E 84 152Google Scholar

    [16]

    Zhang Z Q, Yuan L S, Liu Z, Cheng G G, Ye H F, Ding J N 2018 Comput. Mater. Sci. 145 184Google Scholar

    [17]

    Cao B Y, Sun J, Chen M, Guo Z Y 2009 Int. J. Mol. Sci. 10 4638Google Scholar

    [18]

    Cao B Y, Chen M, Guo Z Y 2006 Phys. Rev. E 74 066311Google Scholar

    [19]

    Cao B Y, Chen M, Guo Z Y 2006 Int. J. Eng. Sci. 44 927Google Scholar

    [20]

    Xie J F, Cao B Y 2017 Mol. Simul. 43 65

    [21]

    Tretyakov N, Müller M 2013 Soft Mater. 9 3613Google Scholar

    [22]

    Plimpton S 1995 J. Comput. Phys. 117 1Google Scholar

    [23]

    曹炳阳, 陈民, 过增元 2006 物理学报 55 5305Google Scholar

    Cao B Y, Chen M, Guo Z Y 2006 Acta Phys. Sin. 55 5305Google Scholar

    [24]

    Schmatko T, Hervet H, Léger L 2006 Langmuir 22 6843Google Scholar

    [25]

    Thompson P A, Robbins M O 1990 Phys. Rev. A 41 6830Google Scholar

  • 图 1  (a) 模拟系统图; (b) 纳米结构示意图

    Fig. 1.  (a) Simulation system; (b) schematic of nanostructure.

    图 2  模型结构示意图 (a) Couette流动; (b) Poiseuille流动

    Fig. 2.  Schematic of nanostructure: (a) Couette flow; (b) Poiseuille flow.

    图 3  不同势能系数c下流体沿z方向的密度分布 (a) c = 2.0; (b) c = 1.0; (c) c = 0.6; (d) c = 0.2

    Fig. 3.  Density profiles in the z-direction with different energy coefficient c: (a) c = 2.0; (b) c = 1.0; (c) c = 0.6; (d) c = 0.2.

    图 4  肋高h对壁面附近流体数密度分布的影响 (a)粗糙壁面; (b)光滑壁面

    Fig. 4.  Effect of rib height h on the distribution of fluid number density near wall surface: (a) Rough wall surface; (b) smooth wall surface.

    图 5  肋间距a对壁面附近流体数密度分布的影响 (a)粗糙壁面; (b)光滑壁面

    Fig. 5.  Effect of rib spacing a on the distribution of fluid number density near wall surface: (a) Rough wall surface; (b) smooth wall surface.

    图 6  不同肋高h下流体沿y方向的速度分布 (a) Couette流动; (b) Poiseuille流动

    Fig. 6.  Velocity profiles in the y-direction with different rib height h: (a) Couette flow; (b) Poiseuille flow.

    图 7  不同肋高h下滑移长度标准差分布

    Fig. 7.  Standard deviation distribution of slip length with different rib height h.

    图 8  (a) 肋高h对滑移长度的影响; (b) 肋高h对滑移速度的影响

    Fig. 8.  (a) Effect of rib height h on the slip length; (b) effect of rib height h on the slip velocity.

    图 9  不同肋间距a下流体沿y方向的速度分布 (a) Couette流动; (b) Poiseuille流动

    Fig. 9.  Velocity profiles in the y-direction with different rib spacing a: (a) Couette flow; (b) Poiseuille flow.

    图 10  不同肋间距a下滑移长度标准差分布

    Fig. 10.  Standard deviation distribution of slip length with different rib spacing a.

    图 11  (a) 肋间距a对滑移长度的影响; (b) 肋间距a对滑移速度的影响

    Fig. 11.  (a) Effect of rib spacing a on the slip length; (b) effect of rib spacing a on the slip velocity.

    图 12  势能系数c对壁面附近流体数密度分布的影响 (a)粗糙壁面; (b)光滑壁面

    Fig. 12.  Effect of energy coefficient c on the distribution of fluid number density near wall surface: (a) Rough wall surface; (b) smooth wall surface.

    图 13  不同势能系数c下流体沿y方向的速度分布 (a) Couette流动; (b) Poiseuille流动

    Fig. 13.  Velocity profiles in the y-direction with different energy coefficient c: (a) Couette flow; (b) Poiseuille flow.

    图 14  不同势能系数c下滑移长度标准差分布

    Fig. 14.  Standard deviation distribution of slip length with different energy coefficient c.

    图 15  (a) 势能系数c对滑移长度的影响; (b) 势能系数c对滑移速度的影响

    Fig. 15.  (a) Effect of energy coefficient c on the slip length; (b) effect of energy coefficient c on the slip velocity.

    表 1  不同模拟工况下对应的粗糙度与接触角

    Table 1.  Corresponding roughness and contact angle under different simulation conditions.

    模拟工况rθ/(°) (c = 1.0)θ/(°) (c = 0.75)θ/(°) (c = 0.5)θ/(°) (c = 0.25)
    h = 0006090120
    h = 0.45σ, a = 3.6σ1.143055.2190124.79
    h = 0.9σ, a = 3.6σ1.28605090130
    h = 1.35σ, a = 3.6σ1.429044.4490135.56
    h = 1.8σ, a = 3.6σ1.571038.2290141.78
    h = 1.8σ, a = 2.7σ1.667033.6490146.36
    h = 1.8σ, a = 4.5σ1.500041.3890138.62
    h = 1.8σ, a = 5.4σ1.444043.8090136.20
    下载: 导出CSV
  • [1]

    黄桥高, 潘光, 宋保维 2014 物理学报 63 054701Google Scholar

    Huang Q G, Pan G, Song B W 2014 Acta Phys. Sin. 63 054701Google Scholar

    [2]

    葛宋, 陈民 2013 工程热物理学报 34 1527

    Ge S, Cheng M 2013 J. Eng. Therm. 34 1527

    [3]

    顾骁坤, 陈民 2010 工程热物理学报 31 1724

    Gu X K, Cheng M 2010 J. Eng. Therm. 31 1724

    [4]

    Sun M, Ebner C 1992 Phys. Rev. Lett. 69 3491Google Scholar

    [5]

    Voronov R S, Papavassiliou D V, Lee L L 2007 Chem. Phys. Lett. 441 273Google Scholar

    [6]

    胡海豹, 鲍路瑶, 黄苏和 2013 力学学报 45 507Google Scholar

    Hu H B, Bao L Y, Huang S H 2013 Chin. J. Theor. Appl. Mech. 45 507Google Scholar

    [7]

    Nagayama G, Cheng P 2004 Int. J. Heat Mass Transfer 47 501Google Scholar

    [8]

    Barisik M, Beskok A 2011 Microfluid. Nanofluid. 11 269Google Scholar

    [9]

    Shi Z, Barisik M, Beskok A 2012 Int. J. Therm. Sci. 59 29Google Scholar

    [10]

    张冉, 谢文佳, 常青 2018 物理学报 67 084701Google Scholar

    Zhang R, Xie W J, Chang Q 2018 Acta Phys. Sin. 67 084701Google Scholar

    [11]

    Cieplak M 2001 Phys. Rev. Lett. 86 803Google Scholar

    [12]

    闫寒, 张文明, 胡开明, 刘岩, 孟光 2013 物理学报 62 174701Google Scholar

    Yan H, Zhang W M, Hu K M, Liu Y, Meng G 2013 Acta Phys. Sin. 62 174701Google Scholar

    [13]

    张程宾, 许兆林, 陈永平 2014 物理学报 63 214706Google Scholar

    Zhang C B, Xu Z L, Chen Y P 2014 Acta Phys. Sin. 63 214706Google Scholar

    [14]

    Rahmatipour H, Azimian A R, Atlaschian O 2017 Phys. A 465 159Google Scholar

    [15]

    Toghraie D, Mokhtari M, Afrand M 2016 Physica E 84 152Google Scholar

    [16]

    Zhang Z Q, Yuan L S, Liu Z, Cheng G G, Ye H F, Ding J N 2018 Comput. Mater. Sci. 145 184Google Scholar

    [17]

    Cao B Y, Sun J, Chen M, Guo Z Y 2009 Int. J. Mol. Sci. 10 4638Google Scholar

    [18]

    Cao B Y, Chen M, Guo Z Y 2006 Phys. Rev. E 74 066311Google Scholar

    [19]

    Cao B Y, Chen M, Guo Z Y 2006 Int. J. Eng. Sci. 44 927Google Scholar

    [20]

    Xie J F, Cao B Y 2017 Mol. Simul. 43 65

    [21]

    Tretyakov N, Müller M 2013 Soft Mater. 9 3613Google Scholar

    [22]

    Plimpton S 1995 J. Comput. Phys. 117 1Google Scholar

    [23]

    曹炳阳, 陈民, 过增元 2006 物理学报 55 5305Google Scholar

    Cao B Y, Chen M, Guo Z Y 2006 Acta Phys. Sin. 55 5305Google Scholar

    [24]

    Schmatko T, Hervet H, Léger L 2006 Langmuir 22 6843Google Scholar

    [25]

    Thompson P A, Robbins M O 1990 Phys. Rev. A 41 6830Google Scholar

  • [1] 王孟妮, 刘珍, 顾浩, 张福建, 张忠强, 丁建宁. 多孔石墨烯时间维度反渗透滤盐机理研究. 物理学报, 2022, 71(13): 138201. doi: 10.7498/aps.71.20212283
    [2] 梅涛, 陈占秀, 杨历, 朱洪漫, 苗瑞灿. 非对称纳米通道内界面热阻的分子动力学研究. 物理学报, 2020, 69(22): 224701. doi: 10.7498/aps.69.20200491
    [3] 陈仙, 张静, 唐昭焕. 纳米尺度下Si/Ge界面应力释放机制的分子动力学研究. 物理学报, 2019, 68(2): 026801. doi: 10.7498/aps.68.20181530
    [4] 李瑞涛, 唐刚, 夏辉, 寻之朋, 李嘉翔, 朱磊. 二维随机蜂巢网格熔断动力学过程和熔断面标度性质的数值模拟. 物理学报, 2019, 68(5): 050301. doi: 10.7498/aps.68.20181774
    [5] 李锐, 刘腾, 陈翔, 陈思聪, 符义红, 刘琳. 界面结构对Cu/Ni多层膜纳米压痕特性影响的分子动力学模拟. 物理学报, 2018, 67(19): 190202. doi: 10.7498/aps.67.20180958
    [6] 叶学民, 张湘珊, 李明兰, 李春曦. 液滴在不同润湿性表面上蒸发时的动力学特性. 物理学报, 2018, 67(11): 114702. doi: 10.7498/aps.67.20180159
    [7] 张永建, 叶芳霞, 戴君, 何斌锋, 臧渡洋. 纳米粗糙度对胶体液滴蒸发图案的影响机制. 物理学报, 2017, 66(6): 066101. doi: 10.7498/aps.66.066101
    [8] 王成龙, 王庆宇, 张跃, 李忠宇, 洪兵, 苏折, 董良. SiC/C界面辐照性能的分子动力学研究. 物理学报, 2014, 63(15): 153402. doi: 10.7498/aps.63.153402
    [9] 张程宾, 许兆林, 陈永平. 粗糙纳通道内流体流动与传热的分子动力学模拟研究. 物理学报, 2014, 63(21): 214706. doi: 10.7498/aps.63.214706
    [10] 唐翠明, 赵锋, 陈晓旭, 陈华君, 程新路. Al与α-Fe2O3纳米界面铝热反应的从头计算分子动力学研究. 物理学报, 2013, 62(24): 247101. doi: 10.7498/aps.62.247101
    [11] 姚祎, 周哲玮, 胡国辉. 有结构壁面上液滴运动特征的耗散粒子动力学模拟. 物理学报, 2013, 62(13): 134701. doi: 10.7498/aps.62.134701
    [12] 周化光, 林鑫, 王猛, 黄卫东. Cu固液界面能的分子动力学计算. 物理学报, 2013, 62(5): 056803. doi: 10.7498/aps.62.056803
    [13] 邱丰, 王猛, 周化光, 郑璇, 林鑫, 黄卫东. Pb液滴在Ni基底润湿铺展行为的分子动力学模拟. 物理学报, 2013, 62(12): 120203. doi: 10.7498/aps.62.120203
    [14] 王陶, 李俊杰, 王锦程. 界面润湿性及固相体积分数对颗粒粗化动力学影响的相场法研究. 物理学报, 2013, 62(10): 106402. doi: 10.7498/aps.62.106402
    [15] 杨平, 吴勇胜, 许海锋, 许鲜欣, 张立强, 李培. TiO2/ZnO纳米薄膜界面热导率的分子动力学模拟. 物理学报, 2011, 60(6): 066601. doi: 10.7498/aps.60.066601
    [16] 石自媛, 胡国辉, 周哲玮. 润湿性梯度驱动液滴运动的格子Boltzmann模拟. 物理学报, 2010, 59(4): 2595-2600. doi: 10.7498/aps.59.2595
    [17] 张程宾, 陈永平, 施明恒, 付盼盼, 吴嘉峰. 表面粗糙度的分形特征及其对微通道内层流流动的影响. 物理学报, 2009, 58(10): 7050-7056. doi: 10.7498/aps.58.7050
    [18] 王晓冬, 董 鹏, 陈胜利, 仪桂云. 亚微米聚苯乙烯微球在气-液界面组装的机理研究. 物理学报, 2007, 56(5): 3017-3021. doi: 10.7498/aps.56.3017
    [19] 王晓冬, 董 鹏, 陈胜利, 仪桂云. 亚微米聚苯乙烯微球在气-液界面组装的机理研究. 物理学报, 2007, 56(3): 1831-1836. doi: 10.7498/aps.56.1831
    [20] 戴永兵, 沈荷生, 张志明, 何贤昶, 胡晓君, 孙方宏, 莘海维. 金刚石/硅(001)异质界面的分子动力学模拟研究. 物理学报, 2001, 50(2): 244-250. doi: 10.7498/aps.50.244
计量
  • 文章访问数:  7939
  • PDF下载量:  90
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-02
  • 修回日期:  2019-01-30
  • 上网日期:  2019-05-01
  • 刊出日期:  2019-05-05

/

返回文章
返回