微尺度电渗混合强化的数值分析研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093498 微尺度电渗混合强化的数值分析 唐桂华，王斐斐, 叶培兴，陶文铨 （西安交通大学，能源与动力工程学院，动力工程多相流国家重点实验室，西安，710049） Tel:029Email:ghtang@mail.xjtu.edu.cn 摘要：应用统一的格子-Boltzmann 模型对微尺度平板均匀通道、壁面粗糙元通道、壁面非均匀非对称 电势分布以及壁面非均匀对称电势分布四类通道中电渗驱动流的外加电场分布、溶液电势分布、流场 以及混合物浓度在同一框架内进行了数值模拟。数值分析表明：与平板均匀通道相比，三类强化通道 都能在一定程度上强化电渗混合，但是由于流动阻力增加，混合时间的延长对混合强化的作用不可忽 略。所研究的三类强化通道中，壁面布置粗糙元通道能够在增加混合时间较小的前提下，获得较高的 混合强度。NPUC 关键词：微尺度；电渗；混合；格子-Boltzmann 方法 0 前言 [1] 由于没有运动部件，外加电场的电渗驱动(EOF)在微流动系统中得到了广泛应用 。 但是，由于EOF 流动系统通常局限于很小的Reynolds 数，惯性影响弱，因此实际应用 中不同流体组分的有效混合变得非常重要。其混合主要由扩散主导，需要较长的混合通 道或者较小的流量以保证足够的混合时间以获得出口的均匀混合。为了强化混合，可采 用主动式强化方法，但是主动式强化方法除了增加损耗外，由于存在运动组件，实际应 用中结果并不理想。因此被动式强化方法广泛应用，其主要措施为在壁面布置粗糙元或 采用弯曲通道，以及壁面布置非均匀电势。以前文献中主要讨论了各种强化方法对混合 效率的影响，而忽略了采取强化措施后，由于流动阻力增加，导致混合时间延长对混合 效果的影响。本文采取格子-Boltzmann 方法（LBM ）统一框架程序，对微尺度平板均匀 通道、壁面粗糙元通道、壁面非均匀非对称电势分布以及壁面非均匀对称电势分布四类 通道中电渗驱动流的外加电场分布、溶液电势分布、流场以及混合物浓度进行统一数值 模拟，重点讨论了混合时间对强化混合效果的影响规律。 1 统一LBM 模型 [2] 采用BGK 模型后的统一格式的离散LBM 方程形式如下 ： δ f i (r + ciδt ,t + δt ) f i (r,t ) =− t [f i (r,t ) − f ieq (r,t )]+ A (1) τ 上式中f 为分布函数，τ 是松弛时间，δt 为时间步长，ci 是粒子速度。对D2Q9 方形格 子，九个速度分量定义为：c0 0 , ci (cos[(i − 1)π / 2],sin[(i − 1)π / 2]) c ：i=1,2,3,4 以及 ci (cos[(i =− 5)π 2 + π 4],sin[(i − 5)π 2 +π 4]) 2c ：i=5,6,7,8 并且满足c δx δt ， δx 为 国家自然科学基金资助项目(No, 新世纪优秀人才支持计划(NCET-07-0676)资助 格子间距。A 是表征外力项的替代因子。(1)式中，f ieq (i 0,1,...8) 表示平衡态分布函数， 对于D2Q9 方形格子，可以表示为统一格式f eq ω B ，权因子ω = 4/9, ω = 1/9 ：=1, 2,