粗糙纳通道内流体流动与传热的分子动力学模拟研究

为研究粗糙表面对纳尺度流体流动和传热及其流固界面速度滑移与温度阶跃的影响,本文建立了粗糙纳通道内流体流动和传热耦合过程的分子动力学模型,模拟研究了粗糙通道内流体的微观结构、速度和温度分布、速度滑移和温度阶跃并与光滑通道进行了比较,并分析了固液相互作用强度和壁面刚度对界面处速度滑移和温度阶跃的影响规律. 研究结果表明,在外力作用下,纳通道主流区域的速度分布呈抛物线分布,由于流体流动导致的黏性耗散使得纳通道内的温度分布呈四次方分布. 并且,在固体壁面处存在速度滑移与温度阶跃. 表面粗糙度的存在使得流体剪切流动产生了额外的黏性耗散,使得粗糙纳通道内的流体速度水平小于光滑通道,温度水平高于光滑通道,并且粗糙表面的速度滑移与温度阶跃均小于光滑通道. 另外,固液相互作用强度的增大和壁面刚度的减小均可导致界面处速度滑移和温度阶跃程度降低. 关键词： 速度滑移 温度阶跃 流固界面 粗糙度     

纳米通道滑移流动的分子动力学模拟研究

本文采用非平衡分子动力学方法对平板纳米通道滑移流动进行了非平衡分子动力学模拟,获得了不同壁面势能和不同温度时流体的速度分布及密度分布。研究结果表明滑移速度在很大程度上决定于流体温度和壁面吸引力作用强度的大小。由于不同壁面吸引力时流体的密度分布受温度的影响规律不同,使得不同壁面吸引力时流体的滑移速度受温度影响规律也不一致。而且,流体结构受壁面流速的影响要受到壁面势能的制约。     

纳米通道内液态微流动密度分布特性数值模拟研究

微通道内流动因表面积/体积比值极大, 造成许多微尺度效应, 进而使微通道内出现完全不同于宏观流动的流体密度分布特性. 本文以纳米通道内液态Poiseuille流为对象, 采用非平衡分子动力学模拟方法研究了流体原子间相互作用强度ε_LL, 流体原子间平衡距离σ_LL以及壁面原子与流体原子间平衡距离σ_LS对通道内流体密度分布的影响规律. 数值模拟中, 统计系综取微正则系综, 势能函数选用LJ/126模型, 壁面设为Rigid-atom壁面, 温度校正使用速度定标法, 牛顿运动方程的求解则采用Verlet算法. 模拟结果表明, 随ε_LL的减弱, 近壁面区密度分布的振荡幅度则逐渐增大; 而σ_LL 则同时影响流体原子的存在形态和密度分布, 较大的σ_LL 会造成流体原子在整个通道内呈现面心立方结构的类似固体排列, 较小的σ_LL会使得流体原子呈现不断变化的 "团簇" 结构; 随σ_LS的变大, 近壁面区流体密度振荡幅度增大, 且流体密度分布起点离壁面越远. 另外, 本文还从近壁面区流体原子的 "俘获-逃逸" 行为角度, 初步解释了原子间相互作用强度对密度分布的影响规律. 关键词： 纳米通道 微流动 密度分布 分子动力学     

双尺度疏油结构纳米通道滑移效应的研究

针对双尺度结构表面疏油特性的优异性,采用分子动力学的方法建立油液流体正十六烷烃分子模型,研究双尺度结构壁面润湿性影响下的纳米通道内流体的流动特性,通过对通道壁面亲疏油性下的双尺度结构的构建,与光滑壁面和单尺度壁面进行比较来探究双尺度纳米通道表面结构影响下油液流体在纳米通道内密度分布、速度分布、速度滑移和滑移长度的影响.模拟结果表明:对于亲油通道壁面,双尺度结构壁面亲油性明显加强,主流区域流体密度、流体速度和速度滑移都减小,甚至出现负滑移;而对于疏油通道壁面,双尺度分层结构能加强壁面的疏油性,通道内壁面形成稳定的气层使流体主流区域的密度增大,并且通道内流体的速度、速度滑移和滑移长度明显大于光滑和单尺度结构壁面.因此,纳米通道内双尺度结构能改变通道壁面的润湿性,并且能够加强流体在纳米疏油通道内的滑移减阻效应,为纳米通道内油液运输以及润滑薄膜减阻提供了设计基础.     

分子动力学模拟纳米通道内混合气体流动的温度效应

温度对纳米通道内流体的流动有显著的作用。运用分子动力学方法,模拟了不同温度下气体混合物在纳米通道内的Poiseuille流动。结果表明:气体混合物化学成分和物理结构都是非均匀的,固壁附近亲水粒子密度随着温度的升高而降低,疏水粒子随着温度的升高逐渐能够到达固壁附近。纳米通道内混合气体在温度较低时有明显的分层现象,而随着温度的升高,密度分布趋于一致。同时在固体壁面从温度较低时的无表观滑移到表观滑移速度随着温度的升高而逐渐增大,而在通道中心混合气体的流动速度随着温度的升高而降低。     

纳米通道内液体流动的滑移现象

采用分子动力学模拟方法研究了液态氩在铂纳米通道内的流动，通过改变流体和壁面之间的势能作用获得了流体和通道表面之间浸润性质不同时的滑移现象. 研究发现：液体分子在亲水性通道表面附近呈类固体性质，数密度和有序性较大，而在疏水性表面附近的平均数密度降低，形成一个低密度层；液体流动在固体表面的速度滑移随着液体与表面势能作用的增强而减小，当液体和表面的浸润性不同时可以发生滑移、表观无滑移和负滑移现象；液体在固体表面的表观滑移是液体在固体表面的速度滑移、粘附和流体内部滑移的综合作用的结果. 关键词： 纳米尺度流动 速度滑移 浸润性 分子动力学模拟     

非对称纳米通道内流体流动与传热的分子动力学

采用分子动力学方法研究了流体在非对称浸润性粗糙纳米通道内的流动与传热过程,分析了两侧壁面浸润性不对称对流体速度滑移和温度阶跃的影响,以及非对称浸润性组合对流体内部热量传递的影响.研究结果表明,纳米通道主流区域的流体速度在外力作用下呈抛物线分布,但是纳米通道上下壁面浸润性不对称导致速度分布不呈中心对称,同时通道壁面的纳米结构也会限制流体的流动.流体在流动过程中产生黏性耗散,使流体温度升高.增强冷壁面的疏水性对近热壁面区域的流体速度几乎没有影响,滑移速度和滑移长度基本不变,始终为锁定边界,但是会导致近冷壁面区域的流体速度逐渐增大,对应的滑移速度和滑移长度随之增大.此时,近冷壁面区域的流体温度逐渐超过近热壁面区域的流体温度,流体出现反转温度分布,流体内部热流逆向传递.随着两侧壁面浸润性不对称程度增加,流体反转温度分布更加明显.     

微通道内电解质溶液离子分布动电学效应模拟

动电学效应对微通道内流体流动特性影响很大,其对通道内粒子分布特性的影响使得通道近壁面流体流动特性极不稳定。本文采用分子动力学方法模拟了二维矩形微通道内NaCl稀电解质溶液的流动特性,考虑存在于不同粒子间的Lennard-Jones势能、静电力、以及带电离子与水分子间的相互作用,得到了粒子在通道内的分布特征。结果显示在动电学效应主要作用于通道壁面附近,而主流区域影响极小。Na~+离子在无量纲通道高度达到0.08和0.91时其浓度达到最大值,沿远离壁面其浓度逐渐降低,与壁面电性相反的Cl~-离子则在无量纲通道高度达到0.15和0.84附近浓度最高。其结果与基于连续介质解理论的Boltzamnn统计分布一致。水分子的浓度在壁面附近也较通道中心处高。     

纳米通道滑移流动的分子动力学研究

本文采用非平衡分子动力学方法对纳米通道内的Couette流动进行了研究,计算获得了不同壁面和流体间势能作用强度下通道内流体的速度分布和密度分布,反映出可能存在的速度滑移、表观无滑移和负滑移现象,并探讨了速度滑移程度和壁面作用强度之间的关系。     

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

纳米流动系统具有高效、经济等优势,在众多领域具有广泛的应用前景.因该类系统具有极高的表面积体积比,致使界面滑移效应对流动具有显著影响.本文采用分子动力学方法以两无限大平行非对称壁面组成的Poiseuille流动为对象,分析了壁面粗糙度与润湿性变化对通道内流体流动的影响.对于不同结构类型的壁面,需要通过水动力位置来确定固液界面位置,准确计算固液界面位置有助于更好地分析界面滑移效应.研究结果表明,上下壁面不对称会引起通道内流场参数分布的不对称,壁面粗糙度及润湿性的变化会影响近壁面附近流体原子的流动特性,由于壁面凹槽的存在,粗糙壁面附近的数密度分布低于光滑壁面一侧.壁面粗糙度及润湿性的变化会影响固液界面位置,肋高变化及壁面润湿性对通道中速度分布影响较大,界面滑移速度及滑移长度随肋高和润湿性的增大而减小;肋间距变化对通道内流体流动影响较小,界面滑移速度和滑移长度基本保持恒定.     

Depolarised Rayleigh light scattering in argon layer confined between graphite plains: MD simulation

The interaction-induced light scattering many-body correlation functions and their spectra in a thin argon layer located between two parallel graphite walls have been investigated by molecular dynamics simulation method. The calculations have been performed for three different distances between graphite plates. Our simulations show the increased intensity of the interaction-induced light scattering spectra at low frequencies for argon atoms in confined space, in comparison to the bulk (unconfined) sample. Moreover, we show a substantial dependence of the interaction-induced light scattering correlation functions of argon on the distances between graphite walls, that is, on the density of argon layer. The mean square displacement and related diffusion coefficient of argon atom in the confined space has been also investigated. Moreover, the structural feature of the thin layer has been studied by calculating the argon density profile, perpendicular to the graphite walls. An interesting observation is the development of a fluid phase in the innermost region of the confined argon layer.     

A novel fluid–wall heat transfer model for molecular dynamics simulations

The ‘fluid–wall thermal equilibrium model’, to numerically simulate heating/cooling of fluid atoms by wall atoms, is used to compare molecular dynamics simulation results to the analytical solution of 1-D heat equation. Liquid argon atoms are placed between two platinum walls and simultaneous heating and cooling is simulated at the walls. Temperature gradient in liquid argon is evaluated and the results are found to match well with the analytical solution showing the physical soundness of the proposed model. Additional simulations are done where liquid argon atoms are heated by both the walls for two different channel heights and it is shown that in such cases, heat transfer occurs at a faster rate than predicted by heat equation with decreasing channel heights.     

Pressure-driven molecular dynamics simulations of water transport through a hydrophilic nanochannel

ABSTRACT

Transport of fluids inside porous materials is relevant to many fields of application. Non-equilibrium molecular dynamics simulation is a powerful technique to explore fluid transport through porous media at the molecular scale. In this work, we compared two commonly used methods for studying pressure-driven transport. The first method was based on the application of an external force field on each fluid particle. The second method made use of two movable walls, acting as pistons, so as to generate transport. These two methods were used to study water transport inside a cylindrical hydrophilic silica nanopore. Several pressure differences were considered from 20 bar to 1000 bar. The results were compared to the theoretical Poiseuille fluid flow. No significant difference was found between the two methods. However, a substantial water flow enhancement was observed compared with the theoretical flow. Both the structural and dynamical properties of water remained unaffected by the applied pressure difference.     

单壁碳纳米管内流体流动的手性效应

纳米尺度范围内流道分子结构形成的粗糙度将会影响其中的流体流动。本文采用分子动力学方法,以氩为工质,模拟了不同手性的单壁碳纳米管（SWCNT）内流体的流动。结果表明：由于范德瓦耳斯力的作用,流体分子与壁而间有一定的距离。可以看出在壁面附近有两个明显的分层,这是纳米流动中持有的流体密度分层现象。由于不同手性碳纳米管的分子排...     

高温壁面润湿性对气层稳定性及其壁面滑移性能的分子动力学研究

针对流体在纳米通道的小尺度效应,采用分子动力学方法模拟了传热效应以及流体流动行为,研究在壁面温度影响下,不同润湿性壁面上方气层生成状态以及流体流动时气层的稳定特性和相应的减阻性能.结果表明:当壁面为纯疏水壁面时,不能形成气层;疏水基底+亲水组合壁面形成不规则气层;纯亲水壁面和亲水基底+疏水组合壁面能形成规则气层.当流体流动时,疏水基底+亲水组合壁面气层消失,而纯亲水壁面和亲水基底+疏水组合壁面气层较为稳定.纯疏水壁面主流区域速度较大,而纯亲水壁面主流区域最低.对于壁面滑移速度,存在气层的壁面滑移速度与纯疏水表面相对接近,甚至稍优于纯属疏水表面,而疏水基底+亲水组合壁面滑移速度最小.     

Ultrasonic flow measurement and wall acoustic impedance effects

An examination of the influence of wall acoustic impedance effects on sound propagation in flowing liquids confined by cylindrical walls is presented. Special focus is given to the importance of the wall acoustic impedance value for ultrasonic flow meter performance. The mathematical model presented allows any radially-dependent axial flow profile to be examined in the linear flow acoustics regime where fluid flow speed is much smaller than the fluid sound speed everywhere in the fluid medium.     

纳米尺度流动速度计算的新方法

采用流体力学中流量与流速的计算和分子动力学相结合的方法,模拟液态氩在纳米通道内的三维Poiseuille流动和驱动方腔流动,计算流体流速.结果表明:平板形纳米通道内,该方法求得的流速与传统分子动力学方法所求流速基本吻合,可以用该方法计算不同壁面情况下的流速;对于纳米方腔通道内的流体,在不同模型下两种方法计算出的流速分布大致相同,但是其边界速度差别很大,在边界的速度计算方面新方法的精确度更高,收敛速度比传统方法快.     

Anomalous flow behavior in closed and open thin walled nanochannels

Molecular dynamics simulations have been carried out to examine water flow in symmetric and asymmetric open and closed ends nanochannels with hydrophilic surfaces. The results are counterintuitive and the opposite of what is observed in macro-systems-closed channels fill faster with fluid than do their open counterparts. In addition, hybrid closed-open asymmetric channels fill up even faster. These results can be explained on the basis of the fluid-structure interaction that arises through the different vibrational behaviors of the surface molecules that are part of the wall forming these channels. Such effects are not expected to be of significance in macro-channels, and point to an important case where macro and nanochannels exhibit contrary behavior. Since these effects results from strong interactions between the fluid molecules and solid surface, one would not expect them with hydrophobic walls, and our simulations confirm such behavior.     

Gas flow characteristics of argon inductively coupled plasma and advections of plasma species under incompressible and compressible flows

In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma properties is specified. In the incompressible flow, only the velocity variable is calculated, while in the compressible flow, both the velocity and density variables are calculated. The compressible flow is more realistic; nevertheless, a comparison of the two types of flow is convenient for people to investigate the respective role of velocity and density variables. The peripheral symmetric profile of metastable density near the chamber sidewall is broken in the incompressible flow. At the compressible flow, the electron density increases and the electron temperature decreases. Meanwhile, the metastable density peak shifts to the dielectric window from the discharge center, besides for the peripheral density profile distortion, similar to the incompressible flow.The velocity profile at incompressible flow is not altered when changing the inlet velocity, whereas clear peak shift of velocity profile from the inlet to the outlet at compressible flow is observed as increasing the gas flow rate. The shift of velocity peak is more obvious at low pressures for it is easy to compress the rarefied gas. The velocity profile variations at compressible flow show people the concrete residing processes of background molecule and plasma species in the chamber at different flow rates. Of more significance is it implied that in the usual linear method that people use to calculate the residence time, one important parameter in the gas flow dynamics, needs to be rectified. The spatial profile of pressure simulated exhibits obvious spatial gradient. This is helpful for experimentalists to understand their gas pressure measurements that are always taken at the chamber outlet. At the end, the work specification and limitations are listed.