力学学报 2018 , 50 (2): 254-262 https://doi.org/10.6052/0459-1879-17-291

流体力学

扩展腔对方波型微混合器混合性能的影响研究

刘赵淼²⁾^,, 王文凯, 逄燕

北京工业大学机械工程与应用电子技术学院,北京 100124

INFLUENCE OF EXPANDED CAVITY ON MIXING PERFORMANCE OF SQUARE-WAVE MICRO-MIXER

Liu Zhaomiao²⁾^,, Wang Wenkai, Pang Yan

(College of Mechanical Engineering and Applied Electronics Technology, Beijing University of Technology, Beijing 100124, China)

中图分类号: O35

文献标识码: A

作者简介:

作者简介：2) 刘赵淼，教授，主要研究方向：微流体力学、计算流体力学、流固耦合分析. E-mail: lzm@bjut.edu.cn

展开

摘要

微混合器凭借节约试剂、混合强度高和易于集成等优点,在材料合成、医药制备和生化检测等领域中具有广泛的应用. 为了进一步提高混合性能,保证混合过程的安全性及生化反应结果的准确性,设计了一种带扩展腔的新型方波型微混合器. 在综合考虑混合强度和压降的前提下,通过实验研究和数值模拟分析了窄缝宽度、窄缝长度和扩展腔高度对微混合器混合性能的影响并得到了不同雷诺数 $Re$ 条件下的最优结构参数. 与方波型微混合器的混合性能进行比较,发现 $Re = 20$ 时,带扩展腔的方波型微混合器的混合强度更高,其中 $Re$ 时两者混合强度相差最多,可达12%. 在相同 $Re$ 下,带扩展腔的方波型微混合器的压降要低于方波型微混合器. 对带扩展腔的方波型微混合器进行内部流场分析,发现扩展腔结构能在流体层流状态的基础上引入涡流,使通道中流体的流动状态发生改变、对流增强,进而混合性能提高.

关键词： 微混合器 ; 混合强度 ; 压降 ; 结构优化 ; 数值模拟

Abstract

Micro-mixer has great application potentials in many fields such as material synthesis, pharmaceutical preparation and biochemical detection due to its advantages of saving reagents, higher mixing index and easy integration. In order to further improve the mixing performance, to ensure the safety of the mixing process and the accuracy of the biochemical reaction results, a new square wave micro-mixer with extended cavity was designed. Under the premise of considering the mixing index and pressure drop, the effects of the width of slit, the length of slit, and the height of extended cavity on the mixing performance of micro-mixer were analyzed by experiment and simulation. The optimal structural parameters were achieved under different Reynolds number ( $Re = 20$ ). Compared with the square-wave micro-mixer, the mixing strength of the square-wave micro-mixer with the extended cavity is higher when $Re$ . Moreover, the gap of mixing index between two square-wave micro-mixer reaches maximum, up to 12%, at $Re < 10$ . Under the same Reynolds number, the pressure drop of the square-wave micro-mixer with the extended cavity is lower than that of the square-wave micro-mixer. Meanwhile, the analysis for the internal flow field of the square-wave micro-mixer with the extended cavity was carried out. It is found that the eddy current is introduced on the basis of the laminar flow state of the fluid because of the existence of the extended cavity structure, which means the change of the flow state of the fluid in the channel and the enhancement of convection effect, thus the mixing performance is further improved.

Keywords： micro-mixer ; mixing index ; pressure drop ; structural optimization ; numerical simulation

PDF (12487KB) 元数据多维度评价相关文章收藏文章

本文引用格式导出 EndNote Ris Bibtex

刘赵淼, 王文凯, 逄燕. 扩展腔对方波型微混合器混合性能的影响研究[J]. 力学学报, 2018, 50(2): 254-262 https://doi.org/10.6052/0459-1879-17-291

Liu Zhaomiao, Wang Wenkai, Pang Yan. INFLUENCE OF EXPANDED CAVITY ON MIXING PERFORMANCE OF SQUARE-WAVE MICRO-MIXER[J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(2): 254-262 https://doi.org/10.6052/0459-1879-17-291

引言

微流控芯片是微型化、集成化的检测与分析器件,凭借其比表面积大、流体表面张力效果明显、节约试剂和安全性高等优点,广泛应用于生物分析、化工合成和临床测试等领域^{[1,2,3,4,5,6]}. 生化反应顺利进行的前提是不同试剂之间的充分混合,设计具有较高混合强度的微混合器具有重要意义^[3,7-8]. 微尺度下流体流动的典型特征是入口雷诺数 $P$ 非常低,流动状态通常为层流,分层流动造成混合困难且混合强度低,通过优化微混合器的结构,进而改变流体的流动状态成为实现微混合器高效混合的研究重点^[8,9,10,11].

被动式微混合器通过优化设计微通道结构诱导产生对流,以打破流体层流状态,使流体间的接触面积增大,缩短分子扩散的距离,进而有助于提高混合强度. 同时,其具有结构简单、运行稳定、易于集成等优势^{[3, 9, 12-15]}. 目前,对于被动式微混合器的研究多集中在T型 ^{[4, 16-19]}、Y型^[20,21]、锯齿型^{[20, 22-25]}以及方波型微混合器^{[20, 24, 26-28]}中.

T型微混合器是结构最简单的被动式微混合器,其只有在高流速、高压下才能实现在短时间内流体的充分混合,满足不了绝大多数生化反应的要求^[19]. Mengeaud等^[23]通过实验对一种Y型入口的锯齿型微混合器进行了研究,在相同雷诺数下,混合效果显著高于T型微混合器. Hossain等^[24]评估了3种不同形状(方波型、多波型和锯齿型)微通道的混合效果,发现在所研究的雷诺数范围内方波型微混合器的混合强度相对较好. Chen等^[25]通过实验和数值模拟的方法研究了6种不同通道形状(方波型、多波型、锯齿型、T型、口型和圆环型)的微混合器对混合性能的影响,发现方波型微混合器的混合强度最大. Chen等^[26]通过数值模拟分析了方波型微混合器通道比和方波单元数对其混合性能的影响,发现通道比和方波单元数的增加有助于混合强度的提高.

为了进一步提高微混合器的混合强度,保证混合过程的安全性及生化反应结果的准确性,可以在微通道中加入不同类型的扩展腔. Johnson等^[29]在T型微通道的底部布置了一系列平行矩形凹槽,发现在 $P$ 的条件下,不同组分的流体完全混合时间小于1 s. Alam等^[30]设计了一种带有矩形腔的多波型微混合器,矩形腔的引入有效地促进了流体之间快速混合. 同样地,陈方璐等^[31]在Y型微通道上设计了多个呈反向分布规则排列的扩展腔,仿真和实验结果表明：该混合器在低流速(0.01 L/s)和高流速(7 L/s)下均可实现流体之间高效快速的混合,同时,在相同入口流速下混合强度远远高于同等尺寸下的Y型微混合器.

方波型微混合器是研究较多且混合性能较好的被动式微混合器,在微通道中引入扩展腔可以有效地提高混合强度. 在实际应用中,压降 $w_{a} = w_{b} = 300$ 是衡量其性能优异的一个重要参数,当 $w_{c} = 300$ 太大时,集成较为困难且容易对微混合器造成一定程度的损坏. 目前,在优化微混合器结构时,综合考虑混合强度与压降两个因素的研究相对较少. 有鉴于此,本文以方波型微混合器为基础,在其转折处引入扩展腔,设计了一种带扩展腔的新型方波型微混合器,实验和数值研究窄缝宽度、窄缝长度和扩展腔高度对其混合性能的影响,并确定最优结构参数.

1 几何模型

本文参考并借鉴Chen等^[25]数值模拟及实验中的方波型微混合器,设计了一种带扩展腔的新型方波型微混合器,如图1所示. 微混合器的入口为相互垂直的Y型入口,不同组分的流体分别从入口1和入口2中流入. 为保证不同组分的流体等量流入微通道中,入口1和入口2的通道宽度相同. 具体通道尺寸如下：入口通道宽度 $w_{1} = 300$ m,出口通道宽度 $H = 1 350$ m,混合主通道宽度 $L_{1} = 1000$ m,混合通道高度 $T = 200$ m,入口通道长度 $w_{2}$ m,通道深度 $h$ m和微混合器总长度10.5 mm. 窄缝宽度 $s$ 、扩展腔高度 $s$ 和窄缝长度 $L_{1}$ 作为结构优化的3个参数. 为确保微混合器总长度不会变化,当 $L_{2}$ 改变时, $L_{3}$ 保持不变,混合通道长度 $s$ 和通道出口长度 $\nabla \cdot V = 0 (1)$ 随 $ρ V \cdot \nabla V = - \nabla p + μ \nabla^{2} V (2)$ 的变化而变化.

显示原图| 下载原图ZIP| 生成PPT

图1 带扩展腔的方波型微混合器示意图

Fig. 1 Schematic diagram of square-wave micro-mixer with extended cavity

2 研究方法及可行性验证

2.1 数值模拟方法

在对微混合器模型进行三维数值模拟求解中,假设流体为不可压缩牛顿流体,流体流动为稳态、层流,壁面为无滑移,忽略重力对流动的影响. 控制方程采用连续性方程、动量方程及组分浓度方程,分别表示如下

$V \cdot \nabla C = D \nabla^{2} C (3)$

$V$

$ρ$

式中, $p$ 表示速度矢量, $μ$ 为流体工质的密度, $C$ 为压力, $D$ 为动力黏度, $\nabla V$ 为组分浓度, $D = 3.23 \times 1 0^{- 10}$ 为组分扩散系数, $ρ$ 为速度梯度. 两个入口分别采用不同流体组分的去离子水作为流体工质,其扩散系数 $μ$ m²/s ^{[3, 32]},其中流体工质密度 $\cdot$ 和动力黏度 $Re$ 分别为998 kg/m³和0.000 97 kg/(m $Re = ρ Ul μ$ s),两种不同流体组分的去离子水之间不发生化学反应,忽略溶解热效应.

利用Fluent 16.0软件对微混合器内部的流体流动状态与混合性能进行三维数值仿真. 微混合器两个入口(inlet)均采用速度入口,出口(outlet)为环境大气,采用压力出口边界条件. 求解器采用SIMPLEC方法,三维双精度、非耦合隐式求解器进行求解,空间离散采用二阶迎风格式.

微尺度下流体混合的关键参数是雷诺数,雷诺数 $U$ 的计算公式如下

$l$ (4)

式中, $w_{a}$ 为入口的平均速度, $T$ 为特征长度,对于本文 $l$ 为300 m, $Re$ 为200 m的带扩展腔的方波型微混合器来说, $M = [1 - \sqrt{1 N \sum_{i}^{=} ‍ {((K_{i} - K_{m}) σ_{m})}^{2}}] \times 100 %$ 为240 m.

在微尺度范围下, $N$ 值很小,一般不超过100,所以微混合器中流体的流动状态主要为层流状态.

采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下

$K_{i}$ (5)

式中, $K_{m}$ 为微混合器出口截面上节点的个数; $σ_{m}$ 为微混合器出口截面上各个节点的组分质量分数; $M$ 为微混合器出口截面上流体完全混合后的各个节点的组分质量分数,为0.5; $0 < M < 1$ 为微混合器出口截面上初始时刻无混合时质量分数的偏差,本文为0.5. $M = 0$ 为混合强度, $M = 1$ ,即 $w_{a}$ 时,代表两种流体完全未混合; $T$ 时,代表两种流体完全混合.

2.2 网格独立性验证

为了保证计算结果的准确性,微混合器模型采用六面体结构化网格进行划分. 对 $Re = 40$ =300 m, $M$ =200 m的微混合器模型分别取75万、150万、215万的计算网格验证网格独立性. 以300万网格为基准,取 $w_{2} = 150$ 下的微混合器出口截面的 $h = 150$ 作衡量标准,所得仿真结果的最大误差分别3.1%, 0.78%, 0.06%. 因此,本文微混合器模型选取的网格量为215万个. 当各个变量间相对残差值小于10 ^-8时,数值计算结果即被判定为收敛.

2.3 实验验证

本文借助高速摄影实验设备(如图2所示)对微混合器内流体的混合效果与流动状态进行研究分析,并与数值仿真结果进行对比. 高速摄影实验设备的型号为基恩士VW-9000. 微混合器芯片采用聚二甲基硅氧烷材质加工而成,去离子水和黑色墨水溶液(用去离子水稀释10倍)作为两种流体工质.

显示原图| 下载原图ZIP| 生成PPT

图2 高速摄影系统示意图

Fig. 2 Sketch of high-speed photography system

利用微注射泵将稀释后的黑色墨水溶液和去离子水等速注入微混合器入口1和入口2中,同时使用高速摄像镜头拍摄流体混合过程的光学图像.

对 $s = 500$ m, $Re$ m和 $~$ m的微混合器进行实验和数值模拟,并计算得到 $M$ 在0.1 $M$ 60下的 $M$ 值,如图3所示.

显示原图| 下载原图ZIP| 生成PPT

图3 实验和数值模拟 $Re$ 值对比

Fig. 3 Comparison of mixed index values in experimental and numerical simulation

从图3中可以看出,实验和数值模拟的 $~$ 趋势基本一致. 当 $M$ 在0.1 $Re$ 1区间时, $Re$ 随 $~$ 的增大呈下降趋势. 当 $M$ 在1 $Re$ 60区间时, $Re$ 随 $x - y$ 的增大迅速提高,两者吻合较好.

$Re = 0.1$ 分别为0.1, 1, 20和40时,微混合器 $M$ 截面上的实验和数值模拟混合溶液质量分数分布图,如图4所示,流体流动的状态、混合程度均吻合良好. $Re = 1$ 时,入口流速较小,流体在通道中混合时间相对较长,不同流体工质分子间的扩散时间延长,导致分子扩散较为充分,因此 $Re = 0.1$ 较高. $M$ 时,入口流速虽然比 $Re$ 时大,但混合时间变短,不同组分流体间的接触面变得清晰, $M$ 降低. $Re$ 继续变大,分子扩散逐渐被对流扩散所取代成为影响 $M$ 的主要因素,随着 $Re$ 的增大,通道内开始出现涡流,有助于提升不同流体工质间的混合,进而 $x$ 逐步提高.

显示原图| 下载原图ZIP| 生成PPT

图4 不同 $y$ 下微混合器 $x$ - $y$ 截面上混合溶液质量分数分布图

Fig. 4 Distribution of mass fraction of mixed solution in $Re$ - $w 2$ cross section of different $M$ under micro-mixer

3 结果与讨论

3.1 窄缝宽度对混合性能的影响

不同 $P$ 下, $Re$ 和 $~$ 随 $h = 150$ 在0.1 $s = 500$ 60范围时的变化( $w_{2}$ m, $M$ m)分别如图5和图6所示. 由图5可知,随着 $w_{2}$ 的减小, $w_{2}$ 增强. 主要有以下两个原因：一方面,随着 $M$ 的减小,窄缝对流体的挤压效果更加显著,导致分子之间的扩散距离和扩散时间缩短,分子扩散更为充分;另一方面, $w_{2}$ 减小,使射流效应增强,更有益于混沌对流,对流的产生更有助于打破层流,导致不同流体工质之间的接触面积增大,促进了流体混合,进而 $Re$ 增大. 当 $~$ 为50 m时, $M$ 在0.1和在10 $Re$ 60区间下, $w_{2}$ 均大于90%.

显示原图| 下载原图ZIP| 生成PPT

图5 不同 $M$ 下 $w 2$ 对 $M$ 的影响

Fig. 5 Effect of $Re$ on $Re$ under different $w 2$

由图6可知,随着 $P$ 的增大,不同 $w_{2}$ 下的 $P$ 均呈上升趋势. 随着 $w_{2}$ 的减小,内部涡流导致能量损失增加, $P$ 增幅剧烈. 当 $w_{2}$ 为150 m, 200 m, 250 m和300 m时,微混合器通道内的 $Re = 60$ 没有明显变化. 当 $P$ 为100 m且 $w_{2}$ 时, $M$ 的最大值为32 kPa,对压力泵和微混合器强度的要求并不大;同时相对于前四种 $w_{2} = 100$ ,在 $Re$ 时的 $w_{2}$ 均大于90%. 因此,将 $P$ m确定为最佳值.

显示原图| 下载原图ZIP| 生成PPT

图6 不同 $w 2$ 下 $P$ 对 $Re$ 的影响

Fig. 6 Effect of $s$ on $M$ under different $P$

3.2 窄缝长度对混合性能的影响

不同 $Re$ 下, $~$ 和 $h = 150$ 随 $w_{2} = 150$ 在0.1 $s$ 60时的变化规律( $M$ m, $s = 600$ m)分别如图7和图8所示. 图7表明 $Re$ 越长, $M$ 越大. 当 $Re = 60$ m时,在不同 $s$ 下,微混合器的 $M$ 均最大,其中在 $P$ 达到最大值,可达99.2%. 这是由于在此 $s$ 下流体混合时间相对延长,更有助于分子扩散, $s$ 最高. 由图8可知, $Re$ 与 $s = 600$ 成正比, $P = 17$ 越长,窄缝对流体的挤压时间相对增长,阻力增加. 其中 $s = 100$ =60时, $s = 600$ m的 $P$ kPa,相比于 $s$ m时增加了48%. 虽然 $s =$ m时的 $Re$ 略微高于其他 $s$ 下的情况,但并不会影响到实际应用,因此选择最优的 $M$ 600 m.

显示原图| 下载原图ZIP| 生成PPT

图7 不同 $s$ 下 $M$ 对 $Re$ 的影响

Fig. 7 Effect of $Re$ on $s$ under different $P$

显示原图| 下载原图ZIP| 生成PPT

图8 不同 $s$ 下 $P$ 对 $Re$ 的影响

Fig. 8 Effect of $h$ on $Re$ under different $~$

3.3 扩展腔高度对混合性能的影响

取 $w_{2} = 150$ 为50, 100, 150, 200, 250和300 m,在 $s = 500$ 为0.1 $M$ 60下进行数值模拟( $P$ m, $Re$ m),得到 $~$ 和 $h$ 的结果分别如图9和图10所示.

图9表明 $M$ 在0.1 $h$ 1之间时,随着 $M$ 的增大, $Re = 0.1$ 不断下降,这是由于随着 $M$ 的增大,不同组分流体间的接触面积减小,导致分子扩散不充分, $h$ 降低. 当 $h = 300$ , $Re$ 在 $~$ 为50 m时达到最大,可达77%,在 $h = 150$ m时最小,为70%. $M$ 在5 $Re = 60$ 60时, $M$ m时的 $h$ 最大,其中 $Re$ 时, $h$ 可高达99%. 这是由于在此 $M$ 下涡流引起的对流扩散效应最为显著.

显示原图| 下载原图ZIP| 生成PPT

图9 不同 $h$ 下 $M$ 对 $Re$ 的影响

Fig. 9 Effect of $h$ on $P$ under different $Re = 40$

由图10可知,随着 $h = 300$ 的增大,微通道对流体的阻碍作用减弱, $P$ 逐渐降低. 其中 $h = 50$ 时, $P$ m下的 $Re$ 值与 $h = 150$ m下的 $Re$ 值相比,减小了5.2%. 考虑到实际工况条件, $h$ 一般不会小于5,所以选取最优的 $P$ m.

显示原图| 下载原图ZIP| 生成PPT

图10 不同 $h$ 下 $P$ 对 $Re$ 的影响

Fig. 10 Effect of $w 2 = 100$ on $s = 600$ under different $h = 150$

综上所述,带扩展腔的方波型微混合器的最佳结构参数为 $h = 150$ m, $w_{2} = 100$ m及 $s = 600$ m.

3.4 两种不同结构微混合器混合性能比较

对 $Re$ m, $~$ m及 $M$ m下的带扩展腔的方波型微混合器与相同尺寸下的方波型微混合器进行数值模拟,得到 $P$ 在0.1 $M$ 60下的 $Re$ 与 $Re = 1$ 的关系曲线图,分别如图11和图12所示.

由图11可知,两种微混合器的 $M$ 均随着 $Re$ 的增大先减小后增大,且均在 $M$ 时, $Re$ 达到最小. 随着 $~$ 继续增大,对流效应逐渐取代分子扩散成为混合的主导作用,通道内产生涡流, $M$ 逐步提高. $Re = 0.5$ 在0.1 $M$ 1时,方波型微混合器的 $Re$ 略高于带扩展腔的方波型微混合器,其中在 $M$ 时两者 $Re$ 最大

显示原图| 下载原图ZIP| 生成PPT

图11 $M$ 对 $M$ 的影响

Fig. 11 Effect of $Re$ on $~$

差值约为2.5%. 这是由于在此雷诺数阶段,混合主要基于分子扩散,而扩展腔的存在导致流体间的接触面积减小,致使带扩展腔的方波型微混合器的 $M$ 略微下降. $Re = 20$ 在5 $M$ 60时,带扩展腔的方波型微混合器的 $P$ 要明显高于方波型微混合器,其中在 $Re$ 时两者 $P$ 相差最多,可达12%. 这是由于随着入口速度的增大,扩展腔内逐渐出现涡流,对流效应明显增强,因此带扩展腔的方波型微混合器的混合效果更好.

图12表明两种微混合器模型的 $Re = 20$ 均随着 $P$ 的增大而增大,同时由于方波型微混合器的沿程阻力更大,导致方波型微混合器的 $Re$ 要高于带扩展腔的方波型微混合器. 其中,在 $P$ 时,方波型微混合器的 $Re$ 值要比带扩展腔的方波型微混合器高3.6%左右.

显示原图| 下载原图ZIP| 生成PPT

图12 $P$ 对 $M$ 的影响

Fig. 12 Effect of $P$ on $Re$

带扩展腔的方波型微混合器相对于方波型微混合器的 $w_{2} = 150$ 更高, $h = 150$ 更小. 对所测的生化反应结果更加准确,而且在复杂的微流控系统中更易集成,混合过程更为安全.

3.5 带扩展腔的方波型微混合器内部流场分析

采用如图2所示的高速摄影系统对微混合器的内部流场进行拍摄,得到两种流体浓度分布的实验结果.

当入口 $s = 500$ 分别取1, 10, 20和60时,微混合器前两个混合单元内两种流体浓度分布的数值模拟与实验分析之间的比较,如图13所示,此时 $Re = 1$ m, $Re$ m和 $h = 150$ m.

从图13中可以看出, $w_{2} = 100$ 时,数值模拟与实验中的流线均为层流状态,没有产生涡流,随着 $s = 600$ 的增大,微混合器内部开始产生涡流并且涡流逐渐增强的趋势在二者之间也基本吻合,特别是涡旋形状和涡心位置也保持一致.

对 $Re$ m, $~$ m及 $x$ m下的带扩展腔的方波型微混合器进行内部流场数值仿真,得到 $y$ 为0.1 $Re = 1$ 60之间微混合器前两个混合单元 $x$ - $y$ 截面上的模拟流线图,如图14所示.

显示原图| 下载原图ZIP| 生成PPT

图13 $x$ , 10, 20, 60时微混合器 $y$ - $Re = 1$ 截面上的流线分布

Fig. 13 Streamline distributions in $Re$ - $~$ plane of micro-mixer at $M$ , 10, 20, 60

从图14中可以看出, $Re = 5$ 在0.1 $Re = 1$ 1时,入口流体速度较小,图14(a)和图14(b)中通道内流体的流动状态均为层流,通道内没有涡流产生,混合主要基于分子扩散, $M$ 呈下降趋势. $Re = 10$ 时,入口流体速度相比于 $Re = 5$ 时增大,图14(c)中一半数目的扩展腔内开始出现涡流,对流效应开始逐渐替代分子扩散成为影响混合性能的主要因素, $M$ 开始上升. $Re$ 时,图14(d)除了在原先扩展腔内出现涡流外,主通道内开始出现较小的涡,扩展腔内的涡流与 $~$ 时相比,涡流强度增大且涡心逐渐向主通道方向移动,增大了流体混合的接触面积,对流明显增强, $~$ 继续提高. $Re = 10$ 在20 $Re$ 60下,图14(e) $~$ 图14(g)中所有的扩展腔内部和主通道内均出现了涡流,涡的数目比 $x - y$ 时增多. 扩展腔内部的涡随着入口流体速度的增大,涡心越来越偏离扩展腔,向主通道方向移动的趋势越来越明显,且主通道内涡流也随着入口速度的增大,涡心逐渐向主通道中间位置移动,导致涡的范围越来越大. 正是因为这些涡的移动,打破了通道内的层流状态,增大了不同流体工质间的接触面积,促进了物质交换,混合性能增强.

显示原图| 下载原图ZIP| 生成PPT

图14 $x - y$ 在0.1 $Re = 0.1 ~ 60$ 60下 $Re$ 截面上的流线图

Fig. 14 The streamline distributions in $~$ plane at $2$

4 结论

(1)带扩展腔结构的方波型微混合器能在流体层流的基础上引入涡流,有效地促进了不同流体之间的混合,有助于提高监测生化反应结果的安全性及准确性. 相比于方波型微混合器,带扩展腔的方波型微混合器在雷诺数 $_{2}$ 时具有更好的混合强度. 虽然两种微混合器在相同 $μ$ 范围0.1{Invalid MML}60区间内表现出类似的压降趋势,但是带扩展腔的方波型微混合器的压降略低于方波型微混合器.

(2)微混合器混合性能的变化规律为设计微混合器提供了重要的理论参考. 优化后的微混合器的体积缩小,流体间的扰动和对流效应增强,不同流体之间的接触面积增大,分子扩散更加充分,混合效果得到进一步提高,更有助于集成和生化反应的顺利进行.

The authors have declared that no competing interests exist.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]	Lee C, Wang W, Liu C, et al. Passive mixers in microfluidic systems: A review .Chemical Engineering Journal, 2016, 288: 146-160 [本文引用: 1]
[2]	李健, 夏国栋. 布置成涡结构微混合器内的流动与混合特性 . 化工学报, 2013, 64(7): 2328-2335 [本文引用: 2] (Li Jian, Xia Guodong. Fluid flow and mixing characteristics in micromixer with vortex-generated structures .CIESC Journal, 2013, 64(7): 2328-2335 (in Chinese)) [本文引用: 2]
[3]	李艺凡, 夏国栋, 王军. 结构参数对布置窄缝和挡板的微混合器内流体流动和混合的影响 . 化工学报, 2015, 66(10): 3857-3865 [本文引用: 5] (Li Yifan, Xia Guodong, Wang Jun. Effect of structural parameters on fluid flow and mixing characteristics in micromixer with gaps and baffles .CIESC Journal, 2015, 66(10): 3857-3865 (in Chinese)) [本文引用: 5]
[4]	任郑玲, 卢晨阳, 王安杰等. T型微混合器合成Cu{Invalid MML}O纳米颗粒 . 化工学报, 2017, 68(6): 2611-2617 [本文引用: 2] (Ren Zhengling, Lu Chenyang, Wang Anjie, et al. Synthesis of Cu{Invalid MML}O nanoparticles in T-shaped micro-mixer .CIESC Journal, 2017, 68(6): 2611-2617 (in Chinese)) [本文引用: 2]
[5]	刘赵淼, 杨洋. 几何构型对流动聚焦生成微液滴的影响 . 力学学报, 2016, 48(4): 867-876 [本文引用: 1] (Liu Zhaomiao, Yang Yang. Influence of geometry configurations on the microdroplets in flow focusing microfluidics .Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(4): 867-876 (in Chinese)) [本文引用: 1]
[6]	罗振兵,孙明波,张攀峰. 第九届全国流体力学青年研讨会报告综述 .力学学报, 2016, 48(2): 511-517 [本文引用: 1] (Luo Zhenbing, Sun Mingbo, Zhang Panfeng. Review of the Ninth National Symposium on Fluid Mechanics for Young Scholars .Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(2): 511-517 (in Chinese)) [本文引用: 1]
[7]	Li Y, Xu F, Liu C, et al. A novel microfluidic mixer based on dual-hydrodynamic focusing for interrogating the kinetics of DNA-protein interaction .The Analyst, 2013, 138(16): 4475 [本文引用: 1]
[8]	姜枫, 刘国君, 杨志刚等. 微混合器的研究现状 . 微纳电子技术, 2016, 53(3): 166-176 [本文引用: 2] (Jiang Feng, Liu Guojun, Yang Zhigang, et al. Research status in micro-mixers .Micronanoelectronic Technology, 2016, 53(3): 166-176 (in Chinese)) [本文引用: 2]
[9]	李战华, 吴健康, 胡国庆等. 微流控芯片中的流体流动. 北京: 科学出版社, 2012 [本文引用: 2] (Li Zhanhua, Wu Jiankang, Hu Guoqing, et al.Fluid Flow in Microfluidic Chips. Beijing: Science Press, 2012 (in Chinese)) [本文引用: 2]
[10]	Hardt S, Drese KS, Hessel V, et al. Passive micromixers for applications in the microreactor and {Invalid MML}TAS fields .Microfluidics & Nanofluidics, 2005, 1(2): 108-118 [本文引用: 1]
[11]	苏祥龙, 许文祥, 陈文. 基于分形导数对非牛顿流体层流的数值研究 . 力学学报, 2017, 49(5): 1020-1028 [本文引用: 1] (Su Xianglong, Xu Wenxiang, Chen Wen. Numerical study for laminar flow of non-Newtonian fluid based on fractal derivative .Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(5): 1020-1028 [本文引用: 1]
[12]	Xia G, Li J, Tian X, et al. Analysis of flow and mixing characteristics of planar asymmetric split-and-recombine (P-SAR) micromixers with fan-shaped cavities .Industrial & Engineering Chemistry Research, 2012, 51(22): 7816-7827 [本文引用: 1]
[13]	Wu C, Tsai R. Fluid mixing via multidirectional vortices in converging-diverging meandering microchannels with semi-elliptical side walls .Chemical Engineering Journal, 2013, 217: 320-328
[14]	Li J, Xia G, Li Y. Numerical and experimental analyses of planar asymmetric split-and-recombine micromixer with dislocation sub-channels .Journal of Chemical Technology & Biotechnology, 2013, 88(9): 1757-1765
[15]	Sabotin I, Tristo G, Junkar M, et al. Two-step design protocol for patterned groove micromixers .Chemical Engineering Research & Design, 2013, 91(5): 778-788 [本文引用: 1]
[16]	肖水云, 李鸣, 杨大勇. T型微混合器内混合强化的数值模拟 . 应用数学和力学, 2016, 37(3): 301-310 [本文引用: 1] (Xiao Shuiyun, Li Ming, Yang Dayong. Numerical Simulation of Mixing Enhancement in T-shaped Micromixers .Applied Mathematics and Mechanics, 2016, 37(3): 301-310 (in Chinese)) [本文引用: 1]
[17]	付强, 鲁聪达. 典型被动微混合器对雷诺数变化的敏感度研究 . 化学工程, 2016 (6): 18-21 (Fu Qiang, Lu Congda. Influence of typical passive micromixer on sensitivity of changing Re .Chemical Engineering, 2016 (6): 18-21 (in Chinese))
[18]	Santana HS, Tortola DS, Jr. Silva JL, et al. Biodiesel synthesis in micromixer with static elements .Energy Conversion & Management, 2017, 141(SI): 28-39
[19]	Wong SH, Ward M, Wharton CW. Micro T-mixer as a rapid mixing micromixer .Sensors & Actuators B Chemical, 2004, 100(3): 359-379 [本文引用: 2]
[20]	毛罕平, 施杰, 杨宁等. 纸质被动式农药微混合器混合特性试验 . 农业机械学报, 2016, 47(7): 73-79 [本文引用: 3] (Mao Hanping, Shi Jie, Yang Ning, et al. Mixing characteristics experiment of paper-based passive micro-mixer for mixing pesticide .Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(7): 73-79 (in Chinese)) [本文引用: 3]
[21]	Hsieh S, Lin J, Chen J. Mixing efficiency of Y-type micromixers with different angles .International Journal of Heat & Fluid Flow, 2013, 44: 130-139 [本文引用: 1]
[22]	Chen X, Li T. A novel passive micromixer designed by applying an optimization algorithm to the zigzag microchannel .Chemical Engineering Journal, 2017, 313: 1406-1414 [本文引用: 1]
[23]	Mengeaud V, Josserand J, Girault HH. Mixing processes in a zigzag microchannel: Finite element simulations and optical study .Analytical Chemistry, 2002, 74(16): 4279-4286 [本文引用: 1]
[24]	Hossain S, Ansari MA, Kim K. Evaluation of the mixing performance of three passive micromixers .Chemical Engineering Journal, 2009, 150(2-3): 492-501 [本文引用: 2]
[25]	Chen X, Li T, Zeng H, et al. Numerical and experimental investigation on micromixers with serpentine microchannels .International Journal of Heat & Mass Transfer, 2016, 98: 131-140 [本文引用: 3]
[26]	Chen X, Li T, Hu Z. A novel research on serpentine microchannels of passive micromixers .Microsyst Technol, 2017, 23: 2649-2656 [本文引用: 2]
[27]	刘贵亚, 张宇, 李民权等. 不同拓扑衍生结构微混合器的仿真与分析 . 微纳电子技术, 2015, 52(10): 633-638 (Liu Guiya, Zhang Yu, Li Minquan, et al. Simulation and analysis of the microfluidics mixers with different topology structures .Micronanoelectronic Technology, 2015, 52(10): 633-638 (in Chinese))
[28]	Pradeep A, Raveendran J, Ramachandran T, et al. Computational simulation and fabrication of smooth edged passive micromixers with alternately varying diameter for efficient mixing .Microelectronic Engineering, 2016, 165: 32-40 [本文引用: 1]
[29]	Johnson TJ, Ross D, Locascio LE. Rapid Microfluidic Mixing.Analytical Chemistry, 2002, 74(1): 45-51 [本文引用: 1]
[30]	Alam A, Kim K. Analysis of mixing in a curved microchannel with rectangular grooves. Chemical Engineering Journal, 2012, 181- 182: 708-716 [本文引用: 1]
[31]	陈方璐, 常洪龙, 洪水金等. 一种基于多重涡流新型微混合器的设计与实验 . 传感器与微系统, 2013, 32(9): 118-121 [本文引用: 1] (Chen Fanglu, Chang Honglong, Hong Shuijin, et al. Design and experiment for a novel micromixer based on multiple vortex .Transducer and Microsystem Technologies, 2013, 32(9): 118-121 (in Chinese)) [本文引用: 1]
[32]	Shih TR, Chung CK. A high-efficiency planar micromixer with convection and diffusion mixing over a wide Reynolds number range .Microfluidics & Nanofluidics, 2008, 5(2): 175-183 [本文引用: 1]
[33]	Hossain S, Husain A, Kim K. Shape optimization of a micromixer with staggered-herringbone grooves patterned on opposite walls .Chemical Engineering Journal, 2010, 162(2): 730-737 [本文引用: 1]
[34]	Tsai R, Wu C. An efficient micromixer based on multidirectional vortices due to baffles and channel curvature .Biomicrofluidics, 2011, 5(1): 14103
[35]	Afzal A, Kim K. Passive split and recombination micromixer with convergent-divergent walls .Chemical Engineering Journal, 2012, 203: 182-192
[36]	Wang R, Lijin B, Shi D, et al. Investigation on the splitting-merging passive micromixer based on Baker’s transformation .Sensors & Actuators B Chemical, 2017, 249: 395-404 [本文引用: 1]

Passive mixers in microfluidic systems: A review

2016

... 微流控芯片是微型化、集成化的检测与分析器件,凭借其比表面积大、流体表面张力效果明显、节约试剂和安全性高等优点,广泛应用于生物分析、化工合成和临床测试等领域^{[1,2,3,4,5,6]}. 生化反应顺利进行的前提是不同试剂之间的充分混合,设计具有较高混合强度的微混合器具有重要意义^[3,7-8]. 微尺度下流体流动的典型特征是入口雷诺数

P

非常低,流动状态通常为层流,分层流动造成混合困难且混合强度低,通过优化微混合器的结构,进而改变流体的流动状态成为实现微混合器高效混合的研究重点^[8,9,10,11]. ...

布置成涡结构微混合器内的流动与混合特性

2013

P

... 采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下 ...

布置成涡结构微混合器内的流动与混合特性

2013

P

... 采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下 ...

结构参数对布置窄缝和挡板的微混合器内流体流动和混合的影响

2015

P

... [3,7-8]. 微尺度下流体流动的典型特征是入口雷诺数

P

... 被动式微混合器通过优化设计微通道结构诱导产生对流,以打破流体层流状态,使流体间的接触面积增大,缩短分子扩散的距离,进而有助于提高混合强度. 同时,其具有结构简单、运行稳定、易于集成等优势^{[3, 9, 12-15]}. 目前,对于被动式微混合器的研究多集中在T型 ^{[4, 16-19]}、Y型^[20,21]、锯齿型^{[20, 22-25]}以及方波型微混合器^{[20, 24, 26-28]}中. ...

... 式中,

p

表示速度矢量,

μ

为流体工质的密度,

C

为压力,

D

为动力黏度,

\nabla V

为组分浓度,

D = 3.23 \times 1 0^{- 10}

为组分扩散系数,

ρ

为速度梯度. 两个入口分别采用不同流体组分的去离子水作为流体工质,其扩散系数

μ

m²/s ^{[3, 32]},其中流体工质密度

\cdot

和动力黏度

Re

分别为998 kg/m³和0.000 97 kg/(m

Re = ρ Ul μ

s),两种不同流体组分的去离子水之间不发生化学反应,忽略溶解热效应. ...

... 采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下 ...

结构参数对布置窄缝和挡板的微混合器内流体流动和混合的影响

2015

P

... [3,7-8]. 微尺度下流体流动的典型特征是入口雷诺数

P

... 式中,

p

表示速度矢量,

μ

为流体工质的密度,

C

为压力,

D

为动力黏度,

\nabla V

为组分浓度,

D = 3.23 \times 1 0^{- 10}

为组分扩散系数,

ρ

为速度梯度. 两个入口分别采用不同流体组分的去离子水作为流体工质,其扩散系数

μ

m²/s ^{[3, 32]},其中流体工质密度

\cdot

和动力黏度

Re

分别为998 kg/m³和0.000 97 kg/(m

Re = ρ Ul μ

s),两种不同流体组分的去离子水之间不发生化学反应,忽略溶解热效应. ...

... 采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下 ...

T型微混合器合成Cu{Invalid MML}O纳米颗粒

2017

P

T型微混合器合成Cu{Invalid MML}O纳米颗粒

2017

P

几何构型对流动聚焦生成微液滴的影响

2016

P

几何构型对流动聚焦生成微液滴的影响

2016

P

第九届全国流体力学青年研讨会报告综述

2016

P

第九届全国流体力学青年研讨会报告综述

2016

P

A novel microfluidic mixer based on dual-hydrodynamic focusing for interrogating the kinetics of DNA-protein interaction

2013

P

微混合器的研究现状

2016

P

... [8,9,10,11]. ...

微混合器的研究现状

2016

P

... [8,9,10,11]. ...

2012

P

2012

P

Passive micromixers for applications in the microreactor and {Invalid MML}TAS fields

2005

P

基于分形导数对非牛顿流体层流的数值研究

2017

P

基于分形导数对非牛顿流体层流的数值研究

2017

P

Analysis of flow and mixing characteristics of planar asymmetric split-and-recombine (P-SAR) micromixers with fan-shaped cavities

2012

Fluid mixing via multidirectional vortices in converging-diverging meandering microchannels with semi-elliptical side walls

2013

Numerical and experimental analyses of planar asymmetric split-and-recombine micromixer with dislocation sub-channels

2013

Two-step design protocol for patterned groove micromixers

2013

T型微混合器内混合强化的数值模拟

2016

T型微混合器内混合强化的数值模拟

2016

典型被动微混合器对雷诺数变化的敏感度研究

2016

典型被动微混合器对雷诺数变化的敏感度研究

2016

Biodiesel synthesis in micromixer with static elements

2017

Micro T-mixer as a rapid mixing micromixer

2004

... T型微混合器是结构最简单的被动式微混合器,其只有在高流速、高压下才能实现在短时间内流体的充分混合,满足不了绝大多数生化反应的要求^[19]. Mengeaud等^[23]通过实验对一种Y型入口的锯齿型微混合器进行了研究,在相同雷诺数下,混合效果显著高于T型微混合器. Hossain等^[24]评估了3种不同形状(方波型、多波型和锯齿型)微通道的混合效果,发现在所研究的雷诺数范围内方波型微混合器的混合强度相对较好. Chen等^[25]通过实验和数值模拟的方法研究了6种不同通道形状(方波型、多波型、锯齿型、T型、口型和圆环型)的微混合器对混合性能的影响,发现方波型微混合器的混合强度最大. Chen等^[26]通过数值模拟分析了方波型微混合器通道比和方波单元数对其混合性能的影响,发现通道比和方波单元数的增加有助于混合强度的提高. ...

纸质被动式农药微混合器混合特性试验

2016

... [20, 22-25]以及方波型微混合器^{[20, 24, 26-28]}中. ...

... [20, 24, 26-28]中. ...

纸质被动式农药微混合器混合特性试验

2016

... [20, 22-25]以及方波型微混合器^{[20, 24, 26-28]}中. ...

... [20, 24, 26-28]中. ...

Mixing efficiency of Y-type micromixers with different angles

2013

A novel passive micromixer designed by applying an optimization algorithm to the zigzag microchannel

2017

Mixing processes in a zigzag microchannel: Finite element simulations and optical study

2002

Evaluation of the mixing performance of three passive micromixers

2009

Numerical and experimental investigation on micromixers with serpentine microchannels

2016

... 本文参考并借鉴Chen等^[25]数值模拟及实验中的方波型微混合器,设计了一种带扩展腔的新型方波型微混合器,如图1所示. 微混合器的入口为相互垂直的Y型入口,不同组分的流体分别从入口1和入口2中流入. 为保证不同组分的流体等量流入微通道中,入口1和入口2的通道宽度相同. 具体通道尺寸如下：入口通道宽度

w_{1} = 300

m,出口通道宽度

H = 1 350

m,混合主通道宽度

L_{1} = 1000

m,混合通道高度

T = 200

m,入口通道长度

w_{2}

m,通道深度

h

m和微混合器总长度10.5 mm. 窄缝宽度

s

、扩展腔高度

s

和窄缝长度

L_{1}

作为结构优化的3个参数. 为确保微混合器总长度不会变化,当

L_{2}

改变时,

L_{3}

保持不变,混合通道长度

s

和通道出口长度

\nabla \cdot V = 0 (1)

随

ρ V \cdot \nabla V = - \nabla p + μ \nabla^{2} V (2)

的变化而变化. ...

A novel research on serpentine microchannels of passive micromixers

2017

不同拓扑衍生结构微混合器的仿真与分析

2015

不同拓扑衍生结构微混合器的仿真与分析

2015

Computational simulation and fabrication of smooth edged passive micromixers with alternately varying diameter for efficient mixing

2016

2002

... 为了进一步提高微混合器的混合强度,保证混合过程的安全性及生化反应结果的准确性,可以在微通道中加入不同类型的扩展腔. Johnson等^[29]在T型微通道的底部布置了一系列平行矩形凹槽,发现在

P

的条件下,不同组分的流体完全混合时间小于1 s. Alam等^[30]设计了一种带有矩形腔的多波型微混合器,矩形腔的引入有效地促进了流体之间快速混合. 同样地,陈方璐等^[31]在Y型微通道上设计了多个呈反向分布规则排列的扩展腔,仿真和实验结果表明：该混合器在低流速(0.01 L/s)和高流速(7 L/s)下均可实现流体之间高效快速的混合,同时,在相同入口流速下混合强度远远高于同等尺寸下的Y型微混合器. ...

181-

2012

P

一种基于多重涡流新型微混合器的设计与实验

2013

P

一种基于多重涡流新型微混合器的设计与实验

2013

P

A high-efficiency planar micromixer with convection and diffusion mixing over a wide Reynolds number range

2008

... 式中,

p

表示速度矢量,

μ

为流体工质的密度,

C

为压力,

D

为动力黏度,

\nabla V

为组分浓度,

D = 3.23 \times 1 0^{- 10}

为组分扩散系数,

ρ

为速度梯度. 两个入口分别采用不同流体组分的去离子水作为流体工质,其扩散系数

μ

m²/s ^{[3, 32]},其中流体工质密度

\cdot

和动力黏度

Re

分别为998 kg/m³和0.000 97 kg/(m

Re = ρ Ul μ

s),两种不同流体组分的去离子水之间不发生化学反应,忽略溶解热效应. ...

Shape optimization of a micromixer with staggered-herringbone grooves patterned on opposite walls

2010

... 采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下 ...

An efficient micromixer based on multidirectional vortices due to baffles and channel curvature

2011

Passive split and recombination micromixer with convergent-divergent walls

2012

Investigation on the splitting-merging passive micromixer based on Baker’s transformation

2017

... 采用混合强度评价微混合器混合性能^{[2, 3, 33-36]}. 混合强度的计算公式如下 ...

〈

〉

扩展腔对方波型微混合器混合性能的影响研究

INFLUENCE OF EXPANDED CAVITY ON MIXING PERFORMANCE OF SQUARE-WAVE MICRO-MIXER

引言

1 几何模型

2 研究方法及可行性验证

2.1 数值模拟方法

2.2 网格独立性验证

2.3 实验验证

3 结果与讨论

3.1 窄缝宽度对混合性能的影响

3.2 窄缝长度对混合性能的影响

3.3 扩展腔高度对混合性能的影响

3.4 两种不同结构微混合器混合性能比较

3.5 带扩展腔的方波型微混合器内部流场分析

4 结论

参考文献 文献选项 原文顺序 文献年度倒序 文中引用次数倒序 被引期刊影响因子

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子