疏水纳米通道内微流动数值计算方法研究 - scholarmate.pdfVIP

疏水纳米通道内微流动数值计算方法研究 胡海豹，黄苏和，宋东，宋保维 (西北工业大学 航海学院，陕西 西安 710072) 摘 要：疏水表面水下减阻研究是当前减阻研究领域的热门课题之一，该类表面一般分布着微纳米尺度的 微结构，其附近流动属微流动问题。以疏水纳米通道内微流动为例，探索了将现有微流动相关理论和方法 应用于疏水性表面减阻研究的可行性，并结合疏水表面界面特性，给出了疏水纳米通道内微流动的分子动 力学数值计算方法（MDS ）。其中，在数学模型中，统计系综为正则系综（NVT ），势能函数采用Lennard-Jones 模型，壁面设定为刚性壁面，壁面疏水性则通过在势能函数中引入修正因子来表征；而在数值计算中，温 度校正使用速度标定法，牛顿运动方程求解则采用 Verlet 算法。实例计算结果表明，疏水微通道内流体密 度呈现内高、外低的振荡型对称分布特点，且随着壁面疏水性的增强，振荡程度减弱；同时疏水壁面处存 在滑移现象，且滑移量随壁面疏水性增强而增大。 关 键 字：减阻，分子动力学，微流动，势能函数，疏水表面，滑移 中图分类号：O357.1 文献标识码：A 自然界中荷叶、水黾等动植物体表面 为此，论文以疏水性纳米通道流动为 均呈现良好的疏水性，低表面能化学组成和 例，开展了疏水微结构表面微流动数值计算 表面微观结构是决定固体表面疏水性的两 方法研究。 个主要因素。近十几年来，疏水表面在流动 1 疏水纳米通道内微流动的可行求解 减阻方面的应用潜力引起了人们广泛的关 方法 注[1] [2] 。Ou J.等 在层流状态下测量了疏水表 微流动学是一门涉及流体力学、热理 面微通道的压降和流量，结果表明其无量纲 物、表面物理、材料学等学科的新兴交叉学 压降比超过了40% 。Choi C. H.等[3]通过对疏 科。在一般的微流动系统中，其通道的水力 水表面微通道内水流动时滑移效果的研究， 直径非常小，表面积/体积的比值极大，出 发现在相同压力下疏水表面微通道内的流 现了许多微尺度效应，如稀薄效应、表面优 率大于常规微通道内的流率，并把流率增大 势效应等[6] 。这些因素使得用于宏观流体运 的原因归结为疏水表面壁面处存在滑移速 动分析的N-S方程不再适用。与常规的微流 度。Lauga E.等[4]基于N-S方程，数值计算了 动问题相比，疏水性纳米通道内微流动具有 圆管内的滑移流动，但由于其引入了充分发 以下特点： 展层流假设等简化模型，致使给出的滑移量 （1）疏水性纳米通道内微流动的介质 明显偏大。陈大融等[5]利用铝基底、聚四氟 为液体，相同体积内其分子数要远大于稀薄 乙烯制备了几种低表面能和表面微观结构 气体，这就要求所采用模拟方法须适用于液 相结合的疏水涂层，并应用到水下航行体的 体介质的微流动。 减阻实验，在小型水洞中获得了最大约 25 （2 ）疏水性纳米通道的壁面具有较低 ％的减阻效果。 表面能，其对液体介质分子的作用力相对较 从目前的研究情况来看，国内外学者 小，模拟方法中须有改变固体壁面内分子对 已充分证明了疏水表面的减阻作用，并得到 液体分子作用力的途径。 了典型试验或数值模拟条件下，特定疏水性 （3 ）疏水性表面的微形貌复杂，且水 微纳米结构表面的基本减阻规律，但相关研 为极性分子，模拟时就需要考虑非常复杂的 究仍存在诸多不足： 分子之间的作