纳米通道粗糙内壁对流体流动行为的影响
|
摘要
纳米流动系统具有高效、经济等优势,在众多领域具有广泛的应用前景.因该类系统具有极高的表面积体积比,致使界面滑移效应对流动具有显著影响.本文采用分子动力学方法以两无限大平行非对称壁面组成的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.
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.
引文
[1] Huang Q G, Pan G, Song B W 2014 Acta Phys. Sin. 63054701(in Chinese)[黄桥高,潘光,宋保维2014物理学报63054701]
[2] Ge S, Cheng M 2013 J. Eng. Therm. 34 1527(in Chinese)[葛宋,陈民2013工程热物理学报34 1527]
[3] Gu X K, Cheng M 2010 J. Eng. Therm. 31 1724(in Chinese)[顾骁坤,陈民2010工程热物理学报31 1724]
[4] Sun M, Ebner C 1992 Phys. Rev. Lett. 69 3491
[5] Voronov R S, Papavassiliou D V, Lee L L 2007 Chem. Phys.Lett. 441 273
[6] Hu H B, Bao L Y, Huang S H 2013 Chin. J. Theor. Appl.Mech. 45 507(in Chinese)[胡海豹,鲍路瑶,黄苏和2013力学学报45 507]
[7] Nagayama G, Cheng P 2004 Int. J. Heat Mass Transfer 47 501
[8] Barisik M, Beskok A 2011 Microfluid. Nanofluid. 11 269
[9] Shi Z, Barisik M, Beskok A 2012 Int. J. Therm. Sci. 59 29
[10] Zhang R, Xie W J, Chang Q 2018 Acta Phys. Sin. 67 084701(in Chinese)[张冉,谢文佳,常青2018物理学报67 084701]
[11] Cieplak M 2001 Phys. Rev. Lett. 86 803
[12] Yan H, Zhang W M, Hu K M, Liu Y, Meng G 2013 Acta Phys.Sin.62 174701(in Chinese)[闫寒,张文明,胡开明,刘岩,孟光2013物理学报62 174701]
[13] Zhang C B, Xu Z L, Chen Y P 2014 Acta Phys. Sin. 63214706(in Chinese)[张程宾,许兆林,陈永平2014物理学报63 214706]
[14] Rahmatipour H, Azimian A R, Atlaschian O 2017 Phys. A465 159
[15] Toghraie D, Mokhtari M, Afrand M 2016 Physica E 84 152
[16] Zhang Z Q, Yuan L S, Liu Z, Cheng G G, Ye H F, Ding J N2018 Comput. Mater. Sci. 145 184
[17] Cao B Y, Sun J, Chen M, Guo Z Y 2009 Int. J. Mol. Sci. 104638
[18] Cao B Y,Chen M,Guo Z Y 2006 Phys.Rev. E 74 066311
[19] Cao B Y, Chen M, Guo Z Y 2006 Int. J. Eng. Sci. 44 927
[20] Xie J F, Cao B Y 2017 Mol. Simul. 43 65
[21] Tretyakov N, Muller M 2013 Soft Mater. 9 3613
[22] Plimpton S 1995 J. Comput. Phys. 117 1
[23] Cao B Y, Chen M, Guo Z Y 2006 Acta Phys. Sin. 55 5305(in Chinese)[曹炳阳,陈民,过增元2006物理学报55 5305]
[24] Schmatko T, Hervet H, Leger L 2006 Langmuir 22 6843
[25] Thompson P A, Robbins M O 1990 Phys. Rev. A 41 6830
[2] Ge S, Cheng M 2013 J. Eng. Therm. 34 1527(in Chinese)[葛宋,陈民2013工程热物理学报34 1527]
[3] Gu X K, Cheng M 2010 J. Eng. Therm. 31 1724(in Chinese)[顾骁坤,陈民2010工程热物理学报31 1724]
[4] Sun M, Ebner C 1992 Phys. Rev. Lett. 69 3491
[5] Voronov R S, Papavassiliou D V, Lee L L 2007 Chem. Phys.Lett. 441 273
[6] Hu H B, Bao L Y, Huang S H 2013 Chin. J. Theor. Appl.Mech. 45 507(in Chinese)[胡海豹,鲍路瑶,黄苏和2013力学学报45 507]
[7] Nagayama G, Cheng P 2004 Int. J. Heat Mass Transfer 47 501
[8] Barisik M, Beskok A 2011 Microfluid. Nanofluid. 11 269
[9] Shi Z, Barisik M, Beskok A 2012 Int. J. Therm. Sci. 59 29
[10] Zhang R, Xie W J, Chang Q 2018 Acta Phys. Sin. 67 084701(in Chinese)[张冉,谢文佳,常青2018物理学报67 084701]
[11] Cieplak M 2001 Phys. Rev. Lett. 86 803
[12] Yan H, Zhang W M, Hu K M, Liu Y, Meng G 2013 Acta Phys.Sin.62 174701(in Chinese)[闫寒,张文明,胡开明,刘岩,孟光2013物理学报62 174701]
[13] Zhang C B, Xu Z L, Chen Y P 2014 Acta Phys. Sin. 63214706(in Chinese)[张程宾,许兆林,陈永平2014物理学报63 214706]
[14] Rahmatipour H, Azimian A R, Atlaschian O 2017 Phys. A465 159
[15] Toghraie D, Mokhtari M, Afrand M 2016 Physica E 84 152
[16] Zhang Z Q, Yuan L S, Liu Z, Cheng G G, Ye H F, Ding J N2018 Comput. Mater. Sci. 145 184
[17] Cao B Y, Sun J, Chen M, Guo Z Y 2009 Int. J. Mol. Sci. 104638
[18] Cao B Y,Chen M,Guo Z Y 2006 Phys.Rev. E 74 066311
[19] Cao B Y, Chen M, Guo Z Y 2006 Int. J. Eng. Sci. 44 927
[20] Xie J F, Cao B Y 2017 Mol. Simul. 43 65
[21] Tretyakov N, Muller M 2013 Soft Mater. 9 3613
[22] Plimpton S 1995 J. Comput. Phys. 117 1
[23] Cao B Y, Chen M, Guo Z Y 2006 Acta Phys. Sin. 55 5305(in Chinese)[曹炳阳,陈民,过增元2006物理学报55 5305]
[24] Schmatko T, Hervet H, Leger L 2006 Langmuir 22 6843
[25] Thompson P A, Robbins M O 1990 Phys. Rev. A 41 6830