2013年第32卷第4期 传感器与微系统(Transducer and Microsystem Technologies)83 SPR传感器阵列中微流通池的流场仿真 张晓娟,曾祥华,刘琳琳,胡 宁,杨 军 (重庆大学生物工程学院生物流变科学与技术教育部重点实验室 视觉损伤与再生修复重庆市重点实验室,重庆400030) 摘要:在自制表面等离子体激元共振(SPR)传感器阵列检测装置的微流通池的研制中,为了获得更好 的检测精度和重复性,需要流经各检测点的流量尽可能一致。设计了几种不同的微流通池结构,并利用 Comsol Muhiphysics软件对其内部的流场进行仿真分析。首先,分析了常用梭状流通池的内部流场分布, 并研究了不同长度过渡段对流场分布的影响。在此基础上,讨论了不同流通池形状,特别是流通池外形轮 廓,如椭圆状流通池和边界流线型更好的流通池对流场分布的影响。此后,还研究了过渡区中微结构阵列 设计对检测区流场分布的影响。仿真研究结果表明:适当延长过渡段长度,有利于改善流体流动的稳定性 和均匀性,但这种改善将随长度的进一步增加逐渐减弱。而不同的边界外形轮廓对检测区域流场的影响 不大;当采用一定的微结构阵列扰流设计后,流场均匀性一定程度上得到改善。 关键词:流通池;仿真;表面等离子体激元共振;流场 中图分类号:Q813.2 文献标识码:A 文章编号:1000-9787(2013)04-0083-05 Simulation 0f flOW field in micro flOW pool in SPR :}= SenS0r arraV ZHANG Xiao.juan。ZENG Xiang.hua,LIU Lin—lin,HU Ning,YANG Jun (Key Laboratory of Biorheological Science and Technology,Ministry of Education, Key Laboratory of Vision Loss,Regeneration and Restoration of Chongqing, School of Bioengineering,Chongqing University,Chongqing 400030,China) Abstract:In research and fabrication of micro flow pool for a self-made SPR sensor array device,in order to achieve high precision and repeatability,flow through each detecting point is required as coincident as possible. Several different flow pool stmetures are designed and simulated by using Comsol Muhiphysics software.At first, inner flow field distribution of common shuttle—shape flow pool is analyzed,and impact of different transition section is analyzed.Based on this simulation,impact of different shapes of flow pools,especially their outlines,on the flow field distribution are also discussed.And then,effect of some microstmctures array design on flow field distribution of detecting area.Simulation results show that appropriately extending the length of the transition section will help to improve the stability and uniformity of the flow of fluid.However,the improvement become smaller when the length exceeded a certain value.While effect of different shapes,on flow field of detecting area is not obvious.Microstructure array design can bring some improvement on the homogeneity of the flow field. Key words:flow pool;simulation;surface plasmon resonance(SPR);flow ifeld 0引 言 等 ]。其中,微流通池容纳样本溶液并使其中的待测物能 有足够的停留时间与敏感膜上的靶标分子发生反应。在 SPR传感器阵列装置的设计中,反应物需要与分布于流通 表面等离子体激元共振(surface plasmon resonance, SPR)传感技术是一项应用最为广泛的光学传感技术 ,具 有无需标记、实时和灵敏等优点 。近年来,随着SPR传 感技术的不断发展,它已经广泛应用于生物医学、食品安全 以及环境监测等领域 。在常见的基于棱镜的SPR生物 池中不同位置的多种靶标分子发生反应。为了克服样本液 在不同位置停留时间(待测分子与靶标分子反应时间)差 异所带来的误差,需要合理的流通池设计,使流经各靶标点 位置的溶液流速尽可能一致。因此,针对一定的分子靶标 传感器中,主要结构单元包含棱镜、传感芯片和微流通池 收稿日期:2012—11 2 基金项目:国家自然科学基金资助项目(81071278,31070882,81101168);教育部“新世纪优秀人才支持计划”资助项目(NCET-09-0842);重 庆大学研究生创新团队建设项目(200909A1002) 84 传感器与微系统 第32卷 阵列分布,研究合适的流通池设计很有必要。 由于传感器设计中高精度、低消耗的要求,越来越多地 采用微型的流通池结构。针对这类微型流通池中的流动情 况的研究也越发引起广泛关注 J。要实测微流通池中流 体速度的分布,如采用常用的粒子图像测速(particle image velocimetry,PIV)法 仍然十分困难。而计算流体力学 (computational fluid dynamics,CFD)仿真则在这一方面具有 传统实验方法难以比拟的优势 。很多研究表明,CFD仿 真所得到的结果能很好地模拟微流场中的实际流动情 况 -l 0_,对相关研究,特别是微流动结构的设计具有重要指 导作用。 本文针对自制SPR传感器阵列检测装置,利用CFD仿 真分析方法来探究合适的微流通池结构。针对SPR传感器 阵列检测的要求,设计了几种不同的微流通池结构,利用 Comsol Muhiphysics软件对其内部的流场进行仿真分析。 根据仿真分析结果可以得到适于SPR传感器阵列检测的微 流通池设计。 1微流通池的仿真模型 实验室自制SPR传感器阵列所采用敏感芯片的金膜尺 寸为12 mii1×8 mill,其中靶点阵列(检测区)范围为2 mm x 2 am,位于流通池中心。在流通池设计时,其尺寸必须覆盖 检测区,但又不能超过金膜的范围。本文设计了5种不同 形状的流通池模型用于分析测试,包括SPR传感器检测中 常用的正方形(图1(a))、圆形(图1(b))、梭形(图1(e))、 椭圆形(图1(d)),以及一种边界流线型更好的结构 (图1(e))。5种流通池进口、出口宽度均为1 mm,进出口 段长度也为1 mm。流通池宽度均为6 mm,前面2种的长度 也为6mm,后面3种为对称结构,长度相等,并可以做一定 的调整,以分析长度变化是否会对流场产生影响。 (a)正方形 (b)圆形 (a)square fb)circle l— I (c)梭形 (d)椭圆形 (c)fusiform fd1 ellipse \ / (e)流线型更好结构 (e)good streamline 图1 5种不同形状的微流通池模型 Fig 1 Five micro flow pool model with diferent geometries 2微流通池的仿真分析 2.1 仿真的控制方程和边界条件 与传感检测实验条件一致,仿真分析中流体介质为水, 流通池内流体运动可视为不可压缩流体的层流,整个过程 为稳态。对于液体的微尺度流动,其分子问距离相比气体 小得多,典型微流控器件尺寸(10—100 m)都远远大于其 分子平均自由程,因此,无滑移条件仍可使用。仿真中,流 通池入口条件均为速度人口,流速设为100 L/min;流通 池出口条件为压力P 0;壁面边界采用无滑移条件。基本 控制方程由Navier—Stokes连续性方程(1)和动量方程(2) 构成 Vu=0, (1) P[詈+“ ]=一 +pg+ “+F, (2) 式中P为流体密度,U为流体速度矢量,P为流体 力, 卵为有效粘性系数,F为附加动力源 、水的流体密度为 1 000 kg/m。,20℃水的有效'枯性系数为l×10一Pa・ 酸 力可忽略不计,仿真过程中无附加动力源, 此,等式(2) 进一步简化为 =叼Wu. (3) 2.2 CFD仿真分析 根据上述控制方程和边界条件,利川Comsol Muhiphy— sies软件对图l中的流通池模型进行二维仿 这类微 流控结构中,二维仿真的结果已足以准确描述其中的流场 分布 。其中,图1中后面3种模型的K度选为宽度的 2倍,5种不同结构的仿真结果如图2所爪. 仿真结果丧 明:溶液在从进口通道进入流通池的过程巾,由 】:流路突然 扩大,相应的流速会逐渐减小。因为静I卜壁 和粘性的影 响,越靠近壁面的地方速度越慢。这样,在微流通池【{ll不同 位置就产生了不同的流速。为_『研究小同流通?也形状对流 场分布的影响,选取通过不同流通池中心( 2(a))的轴 向(图3(a))和横向(图3(b))上的流场进行比较. .仿真结果表明:梭状、椭圆形、流线 更好的流通池,轴 向流速在检测区稳定性要比正方形和圆肜流通池好僻多, 几乎没有太大变化(『皋J 3(a))。在横 速度分布上,这3种 结构的变化率也要小很多(图3(b))。从图2可以看m:[}{ 于流通池四角方向滞留区(流速接近_r 0)的存 ,形状的 影响相对较小,因此,后面3种结构流场分布的改善很 】‘能 是由于过渡区(人口到检测区,如图1(e) )的延K所rjl 起的。因此,下面将讨论这一改变的影响。 3流通池过渡段的影响 以梭形流通池为例来研究过渡段长度 改变对检测区 流场分布的影响。假设,J长度分别为1 lllln(图4【1J曲 线1),2Intn(图4中曲线2),3 mm(图4中曲线3),4 nIn 传感器与微系统 第32卷 ∞.口 0_【\嘲醣 T∞.蛊 0一\ 煺 (a)r=0.15 mm micropillar array (b)Y-方向 (b)Yl direction (d) 方向 (d) 1 direction y/mm (d) 方向 (d) 2 direction 图5不同微柱尺寸对流场的影响 Fig 5 Impact of size of micropillar on flow ifeld 4.2微柱个数对检测区流场分布的影响 前面的研究表明,半径0.15 mm微柱阵列对流场均匀 性的效果较好,因此,在后面的研究中采用该尺寸的微柱。 为了比较微柱个数的影响,选用了如图5(a)的2×3(2列 3行,图6中曲线1)、2×4(图6中曲线2)、3×3(图6中曲 线3)、3×4(图6中曲线4)的微阵列进行仿真,图6中曲 线5为未加微柱结构。相比于未添加微阵列时,在轴向方 向前2种情况的速度分布更均匀(图6(a)和(b)),横向方 向后2种更均匀(图6(C)和(d))。这表明增加阵列列数有 利于横向流速均匀性,但同时也降低了轴向流速的稳定性。 增加阵列的行数影响相对较小,3行的阵列效果稍好。 4.3微柱分布对检测区流场分布的影响 采用微柱半径为0.15mm的2×3阵列,研究微柱位置 s.茸 0一\ 斌 T∞.盘 0一~ 避 呈 媛 (b)y 方向 (b) direction 呈 避 呈 y/mm (d】 。方向 (d) 2 direction 图6不同微柱个数对流场的影响 Fig 6 Impact of number of micropillar on flow field 分布改变对流场的影响。选用了5种不同的位置分布,行、 列间距分别为0.5,0.5mm(图7中曲线1),0.5与0.75mE, 第二列与曲线1一致(图7中曲线3),0.5,0.75mm,第一列 与曲线l一致(图7中曲线3),0.75,0.5mill,中间行保持在 中轴线上(图7中曲线4),0.75,0.75 mm,行分布与曲线4 一致,列分布与曲线2一致(图7叶1曲线4),未加微柱结构 (图7中曲线6)。仿真结果表明:微柱阵列有利于稳定横 向流场,增大微阵列中排与排的距离对轴向流场有一定的 改善作用,而对横向流场则不利,增大列与列之间的距离对 流场的影响不明显。 5结论 本文利用仿真软件分析了5种不同边界的流通池的内 部流场分布,并研究了不同长度的过渡段对梭状流通池流 场分布的影响。在此基础上讨论了通池流中增加微结构阵 列产生扰流后对流场分布的影响。研究结果表明:对于不 同的边界形状,检测区域流场变化不明显;而进口到出口的 长度不同,检测区流速的稳定性差异较大;如果在微流通池 第4期 T∞.口J卜0一\煳蛙 24 2 2 2 O 1.8 张晓娟,等:SPR传感器阵列中微流通池的流场仿真 T∞.鼻-(】【/ 87 参考文献: [1]Shagnfa H K,Kriszta F,Raj K,et a1.A versatile method to men— sure the binding to basic proteins by surface plasmon resonan— ee[J].Analytical Biochemistry,2012,421(2):385--390. 1.6 1 5 —1.0 —0.5 0 0.5 1.0 1.5 [2]“Y,“u X,Lin Z.Rencent developments and applications of sur. face plasmon resonance biosensors for the detection of mycotoxins (a)y1方向 (a)Y1 direction 2.4 2 2 in foodstuffs[J].Food Chemisty,r2012,132(3):1549--1554. [3]0u H C,Luo Z F,Jiang H,et a1.Indirect inhibitive immunoassay or estfradiol using suI ̄'aee plasmon resonance coupled to online in. 2 2 0 tube SPME[J].Analytical Letters,2009,42(17):2758—2773. 1.6 —[4]杨一I.0—0.5 0 0 5 l 0 1.5 军,黄小玲,张丽果,等.一种小型SPR生物传感检测装 i.5置的研制[J].仪器仪表学报,2011,32(7):1578--1584. [5] Gaspar A,Gomez F A.Development of an ultra-low volume flow cell or suIfface plasmon resonance detection in a miniaturized ca. (b)y2方向 (b)y2 direction 2 5 pillary electrophorosis system[J].Electrophoresis,2012,33(12): 2.0 1723--1728. 1.5 1.0 0 5 1.5 —1 0 [6]Sheppard G,Oseki T,Baba A,et a1.Thiolene—based microfluidic lfow cells for surface plasmon resonance imaging[J].Biomicro— —0 5 0 0.5 1 0 1 5 lfuidics,2011,5(2):026501. (C) 。方向 fc1 1 direction 2・5 — [7]Hirono T,Arimoto H,Okawa S.Microfluidic image cytometry for measuring number and sizes of biological cells flowing through a microchannel using the micro-PIV technique[J].Measurement Science&Technology,2008,19(2):025401. 娄 2。0 呈1.5 o 0 5 [8] Glatzel T,Litterst C,Cupelli C.Computational fluid dynamic. s(CFD)software tools for microfluidic applieations--A case stu— 一l 0 -0.5 0 0 5 1.0 1.5 1. dy[J].Computers&Fluids,2008,37(3):218--235. [9] Yang J,IJi C W,Yang M S.Hydrodynamic simulation of cell doc— king in microfluidic channels with different dam structures[J]. Lab Chip,2004,4(1):53-59. y/mm (d) :方向 fd1 X2 direction 图7不同微柱分布对流场的影响 Fig 7 Impact of layout of micropillar on flow field [10]rang M S,Yang J,Li C W,et a1.Generation of concentration gra— dient by controlled flow distribution and diffusive mixing in a mi— 中设立适当长度的过渡段,有利于改善流体流场的均匀分 crofluidic chip[J].Lab Chip,2002,2(3):158--163. 布;在梭状流通池中添加微结构阵列扰流后,流场均匀性一 定程度上得到改善。这一结果对SPR传感器阵列检测中流 通池阵列的设计有重要指导作用。 作者简介: 张晓娟(1987一),女,四川I洪雅人,硕士研究生,主要研究方向 微流控方法及生物传感器。 \ ; t声 ) (上接第82页) 时间误差传递模型。本文以5O s的同步精度需求设计为 例,采用理论分析和实验验证相结合的方式对时间管理模 型进行了验证,结果表明:本文采用的时间管理方法和误差 出版社,2009. [3] GJB 289A~1997数字式时分制指令/响应型多路传输数据总 线[s]. [4] GJB 2994-1997航天器数据系统时间码格式[S]. [5]谭维炽,顾莹琦.空间数据系统[M].北京:中国科学技术出 版社,2004. 模型能够满足遥感卫星高精度时间同步的设计和分析需 求,为遥感卫星星上时间管理的方法选取和方案设计提供 了便利的条件。 参考文献: [1] caHam W G,Goodman R S,Majewski R M.Spotlight synthetic aperture radar signal processing algorithms[M].Boston,London: Artech House,1995:357--496. [6]李新贝,谭超,高山,等.星上RS122接口电路的建模与 仿真[J].航天器工程,2011,20(1):109--113. [7]Gardner F M.Phaseloek techniques[M],New York:Wiley& Sons,2005. 作者简介: 田贺祥(1981一)男,河北保定人,博士研究生,主要研究方向 为遥感卫星数据系统与综合电子技术。 [2] 谭维炽,胡金刚.航天器系统工程[M].北京:中国科学技术