10 沟槽微流道的牛顿流动.docVIP

第10章 沟槽微流道的牛顿流动 D.G. HASSELL* and W.B.J. ZIMMERMAN Department of Chemical and Process Engineering, University of Sheffield, Newcastle Street, Sheffield S1 3JD United Kingdom * Present address: Department of Chemical Engineering, Cambridge University, Pembroke Street, Cambridge, CB2 3RA E-mail: w.zimmerman@shef.ac.uk 对微型系统和微流道技术研究的发展趋势已经为人们对于尺度规模为“芯片上的实验室”或“板上的工厂”的化学和生物系统得到更好的控制提供了令人向往的视野。由于传统的混合（具有混沌特性的湍流）方法不能胜任系统相当高压力的需求，使得所有这些微小尺度系统的流体混合问题引起人们的兴趣。对于某些系统，仅考虑分子扩散已经足够，但对于高分子（例如DNA）和低分子扩散系数系统，它需要更精确的方法，由积极或激励的方法对起始隔离的流线进行诱导，使得该可控制对流混合成指数倍地增强。本实例分析研究由Stroock 和合作者所提出的错排人字形混合器（SHM）中的流动模式，同时确定影响这些装置中的混合的各种参数之间的关系。在比较简单平行板中流动的分析解和COMSOL Multiphysics的解之后，接着扩展到不同流动条件下的两种测试例子的三维模拟。然后通过后处理方法对所感兴趣的特性进行估算，从而对这种复杂流动的流体力学有更好的理解。 1．简介 如图1所示的错排人字形混合器（SHM），它作为一种在不产生大的压降的前提下增强微流道装置中的混合的方法被提出。 图1. Stroock 和合作者提出的微混合器例子，图中的单位为微米。薄片由两组非对称的人字形沟槽逐次排列组成，偏轴距离为流道宽度的六分之一。 它采用由两组非对称的人字形沟槽逐次循环排列组成，使流过其中的流体发生折叠，就象面包工人在揉生面团一样。这种对流式折叠使发生扩散混合的面积增大，针对这种混合的几种数值计算研究已经展开，而关于这些系统的速度分布曲线和相关的流体特性还未有研究。本章重点介绍沟槽对不同测试条件下的体相流道流动的影响。鉴于COMSOL Multiphysics 软件对于这些流动求解能力的限制，我们使用两种不同的几何体。第一种是单一槽道的流道，另外一种是对具有周期性边界条件的无限大对称沟槽进行模拟。在对更加复杂的三维区域进行分析之前，我们先从平行平板中简单流动的数值方法估算分析谈起。流体的流动由无量纲Navier-Stokes方程组来模拟，式中根据入口速度等于1来设定，被设定为流道宽度的大小。 , (1) , (2) , (3) , , . (4) 我们通过雷诺数（Re）来定义惯性力和粘性力的比值，引入雷诺数（Re）被作为一附加的常数。这种方法使COMSOL Multiphysics具有简化问题的优点，它通过减少量纲为1的变量个数，使问题更容易求解。 2．2-D牛顿流动 我们通过演示COMSOL Multiphysics对如表1所示Poiseuille流简单的2-D模拟开始。速度分布曲线如图2所示，图中还给出了流道方向上的速度分布曲线。速度分布曲线表明，流体从起始的活塞流发展为Poiseuille流需要一定的时间，出口条件也影响出口附近的速度分布。所以只有接近中心的流动才是真正的Poiseuille流。 解决该问题的另外一种方法是采用周期性边界条件，进口和出口的速度被设定为模拟无限长流道中的流动的理论值。 . (a) 速度值大小 (b) 速度分布 图2. (a) 沿2-D流道的速度场和 (b)流道中相应不同截面在x方向的速度分布. 表1. 2-D牛顿流动 Model Navigator 选择2-D空间维数 选择 COMSOL MultiPhysics|Fluid Dynamics|Incompressible|Navier-Stokes|Steady-state Add. OK 下拉菜单 选择 Draw objects. 创建高度为1和宽度为2的矩形， 左下角为坐标的原点 物理菜单： 边界条件（F7） 选择边界1条件类型：Inflow/Outflow. 输入u