微通道管液滴形成机理影响分析.docVIP

微通道管液滴形成机理影响分析 　　【摘 要】在微通道管液滴形成的数值模拟过程中通常采用固定速度入口，而且主要集中于研究液滴外形的变化。为了更加系统全面的研究液滴的形成，本文使用固定压力入口，主要对液滴内部轴线速度的变化进行细致研究，运用有限体积法求解液滴形成的N-S方程，采用VOF方法以及基于PLIC 的几何重构方法追踪液滴运动过程的自由表面， 捕捉气液界面的演化过程。为微流控设备的设计和研究提供新的数值预测方法和参考。 　　【关键词】计算流体力学方法；微通道管；液滴形成；数值模拟 　　0 引言 　　近年来，随着微流体技术的不断发展，基于液滴的微流控系统已经在生物、化学和医学上得到了广泛的应用，能够实现液滴在微小通道中的流动控制。由于该技术应用极广，理论尚未完善，至今仍是各国学者研究的重点。对液滴的研究始于十九世纪末期。早期的研究主要集中在实验方面， 张略等设计了由音圈电机驱动的非接触喷射机构，在喷墨打印中能有效减少液滴拖尾现象[1]；Furbank 等进行了含微小颗粒流体的液滴形成实验，过程复杂且实验研究的成本非常高[2]。近年来，实验方法的局限性在研究微通道几何尺寸对液滴形成的影响时十分明显，需要设计制作大量不同尺寸和形状的芯片，使得实验研究的成本过高，导致对微通道中微液滴形成机制的研究进展缓慢[3-5]。随着流体力学理论的发展，计算机仿真技术的提高，研究的重点越来越倾向于数值实验， Chan 等运用遗传算法和模糊神经网络建立滴化过程模型；Sirotkin 提出运用CSPH方法模拟两相流，并分析其可行性。通过数值模拟不仅可以节约成本、节省时间，对液滴形成过程中跟随液滴出现一些较为复杂的问题能够进一步理解。 　　1 研究方法 　　本文运用基于有限元思想的VOF模型求解N-S方程和连续性方程。连续相（continuous phase）为液相，离散相（dispersed phase）为气相，并且假设两相均为黏性不可压的，两相互不相容，只存在动量的传递，不涉及传热问题。模拟过程中采用固定压力入口，不虑重力的影响。 　　1.1 控制方程 　　由于本文的计算模型是轴对称的，因此只需要取经过对称轴的一个面进行计算即可，将三维问题转化为二维问题，研究的流体为不可压缩牛顿流体，流动处于稳态，对二维模型连续方程为： 　　（3）中v代表速度矢量，m?s-1；P为压力，Pa；σ为表面张力系数，N?m-1；μ为动力学黏性系数，kg?m-1?s-1；ρ为液体密度，kg?m-3；M（x，y，z，t）为每个控制体内各相的体积分数，是用VOF方法定义的流体体积标量函数；K为界面曲率。 　　1.2 数值计算 　　针对所研究的为轴对称问题，取经过对称轴的一个平面进行计算，模型的网格如图1所示。 　　微管的口1设置为固定压力入口边界条件。2、3设置为壁面无滑移、无穿透边界条件，其中2 的接触角设置为60°，3 的接触角设置179°。计算域出口4、5定义为远场（Far field）边界条件，参考压力点位于4、5的交点处，其压力值为0.1MPa，6为对称轴。模拟过程中连续相为甲醇，分散相为空气，整个过程为恒温300K。数值模拟采用有限体积法（Finite Volume Method）对液滴形成与破碎过程N-S方程进行离散，压力场采用体积力分数计算方法（ body-force-weighted） ，这种方法在计算面压力时假设体积力和压力间差值梯度为常数，当体积力在动量方程中可得到时，应用此方法可以获得较好的效果。本文采用轴对称模型进行计算适合应用这种方法。速度压力耦合计算采用基于交错网格的PISO算法，这种方法允许使用较大时间步长，有利于节省计算时间，并且该算法具有两个修正步骤（corrector step ） ，对于非稳态流动问题采用这种算法更为精确。动能项的离散方式采用二阶迎风格式（second order upwind）。对于自由表面的处理， 使用VOF的几何重构对交界进行处理，并引入连续表面张力模型。几何重构基于PLIC方法，这种方法能较好地跟踪自由表面。 　　2 液滴形成过程分析 　　本文采用有限差分法配合体积分数法VOF模拟液滴的形成和下落。图2（a）到图2（f）为数值模拟液滴形成过程中几个典型的序列，可以将液滴的形成分为初始阶段、颈缩阶段和断裂阶段。 　　2.1 液滴形成的初始阶段 　　2.2 颈缩阶段 　　状液滴开始呈现出类椭圆的形状，并且，在此时出现了颈缩区域，如图2（d）所示。 　　2.3 断裂阶段 　　3 结果与讨论 　　3.1 不同时刻液滴颈缩前轴线速度变化 　　在模拟的过程中，本文所采用的是固定压力入口，图（3）是在固定压力的作用下，液滴在颈缩产生之前几个时刻轴线速度的变化。从图中可以明显的看出，随