在两相流模式和在含水甲醇与CO2混合物在绝热的矩形微通道的过渡标准

1、在两相流模式和在含水甲醇与CO2混合物在绝热的矩形微通道的过渡标准C. Wei nmueller, N. Hotz, A. Mueller, D. P oulikakos *摘要：本文介绍了在对绝热双相流的流程作业图调查，200微米的硅制深的微通道，采用激光感受荧光检测。利用两个的管嘴宽度(管口分别为30微米和50微米)和以甲醇水和 CO2溶剂为气相，来研究表面张力和管嘴几何学对流动模式转变的影响。结果发现和量化的以几何尺寸更小的管口和更小的液体表面张力，可以在更低的表面气体和液体流速促进毛细管气泡传播。在表面气体的测量域(0.01 -0.625 m/s)和液体速率的测量域(0.0005 -0

2、.5000 m/s)内，我们观察到分散的气泡，有规律排列的气泡，楔入，蛞喻和环形流动，从而扩大了较小 表面液体速率近两个数量级的实验知识基础。在这里提出的流程图的帮助下，我们可以表征观察到的规律地有序气泡流动为气泡的分散流和楔入流的传质状态。本调查的结果是小规模甲醇燃料电池的运作的直接关系。1.简介:如今多相流是一个一直发展的课题，与微型流动的应用相比，它可以与归类于大规模的 应用。经典两相流的发现包括能源和工艺技术广泛领域的大型应用，比如热量传递(Kattan etal,1998; Kew and Cornwell,1997;Thonon et al, 1995;Webb,2004),化学和

3、生物反应器(Darmana et al,2005;Dudukovic et al,1999;Millies and Mewes,1999;Va nhoute n et al ， 1994)。由于电子工 业的小型化和分散的能源转换系统应用的趋势已经开始应用于主机的创新，比如应用于芯片(Romera-Guereca et al，2008)，微芯片冷却(Amon et al,2001;Escher et al,2009;Rosengarten et al,2006;Senn and Poulikakos,2004)或者是微型燃料电池 (Choban et al,2004; Hotz et al,20

4、06a,b;Senn and Poulikakos, 2005)。为了设计高效率的设备，我们需要对引入微管路得预期双相流 有一个好的认识。很多应用和科学技术的出现(Gu nther and Jen se n,2006;Joa nicot and Ajdari,2005; Kobayashiet al,2004)和引入的微观物理现象的复杂性结合在一起，组成了一个极具挑战的研究课题。即使管路中的大规模两相流很好的了解了(Baker,1954;Barnea and Taitel,1993;Lockhart andMarti nelli,1949; Ma ndha ne et al,1974;Tait

5、el and Dukler,1976)。然而宏观两相流受到重力，惯性(剪切)力和表面张力的相互作用很大的影响，微观或者毛细管的两相流偏离对于重力影响的逐渐减小和其他提到的力量的重要性的转变的描述。Kreutzer et al.( 2005)提出了微观和毛细管两相流的分类，该分类依赖于Bretherton (1961)分析地导出的临界邦德数：这个标准为了长气泡的垂直向上移动和在浮力推动的低雷诺数得圆管流而导出的。在低的邦德数和毛细管数下(Ca=n u/ q)泡沫运动可用润滑理论来描述。但是，邦德数低于阈值(关 键)时，泡沫运动终止。尽管在这个分类的近似性质中，我们不考虑去湿现象， 粘度或者不同的

6、通道横截面，但是它还是被用于本文，因为它提供了将毛细血管/微观两相流分类得简单方法。临界邦德数在定义微流双水相现象的重要性时，其毋庸置疑的优点有表面张力的引入，在仅仅基于管路尺寸的再分类上的明显改进等。Tabk- 1PJrimcter SiMs.ClunricJ 曲Qd记try (pmjNozzle wkkh (pm)Liquid fluidSet 120050HzOSet?200501 M Me OHSCI 320050MeOHSet 4200301 M Me OH宏观实验得到的相关关系是否仍然成立，同一批实验发挥作用到何种程度， 这些问题都需要更深入的研究。为此，很多当代的论文(Colem

7、a n and Garimella, 1999; Cubaud and Ho, 2004;Hassan et al., 2005; Kreutzer et al., 2005; Trip lett et al., 1999; Waelchli and von Rohr, 2006; Zhao and Bi, 2001)报道了这一领域的进展，然而表面张力和喷嘴几何学性质对相转变的影响没有 考虑进去。本研究重点关注表面张力的影响，利用检测三种不同组成的甲醇-水液态混合溶液，CO2为气相，见表.1。除此之外，一组喷嘴集合体(表 .1)被研究来推断对流动过渡状 态的影响。工作液体的选择也要考虑它与小规

8、模直接甲醇燃料电池的相关性。在表面张力和喷嘴几何学的影响能够被论述清楚之前，一个参考情形(参数组1,见表1)在文献中被确立，来客观比较类似参数(流体性质，管道几何学和制造性质)。除此之外，对在低液体流动速率下的流动相图未来的需要(Gunther and Jensen, 2006; Hassan et al, 2005)，在该实验中有所阐述。两相流动转变的现有知识被扩展了近两个数量级在边长为200微米的方形微通道中。2.实验步骤2.1.精细加工一个基于光蚀刻影响标准的精细加工/微机电系统技术(Madou, 2002)被应用于生产微尺寸的管道网络，制造步骤见图.1。在第一步，一个500 25微米厚

9、，磷掺杂，双面抛光，4、阻抗达到30欧/厘米的硅圆片都被涂上了一层有着负性光致抗蚀剂组成的10微米厚圆片(AZR4562, MicroChem)。在随后的光刻过程中 (MA6, S ss Micro Tec)，展开高分辨率面膜极性的 掩饰极性模仿在UV照射下的光刻技术。一个干蚀刻程序用电感耦合等离子体系统(ICP备表面技术系统)来操作，最终形成了深和宽都是200微米的矩形管路截面。管道壁的表面粗糙度是由引入的“ Bosch process和掩饰分辨率决定的，并且不超过1微米。在硅圆片的背面，通过上面提到的步骤的循环，液体进口、出口和压力端口合并在一起。为了光通过，一个直径100mm , 700

10、 50厚的硼硅酸盐玻璃圆片被通过由衬底焊机(SB6，苏斯微技)被焊接到硅圆片上，在450摄氏度和1000V电压下进行。在将圆片切割成预先规划的为管路碎片后， 用环氧基树脂将聚醚醚酮 (Up church Scie ntific) 流体连接器与进口、出口和压力管密封。在芯片上的微流网络是按照典型的微型反应设计的，这是为了与微型燃料电池相适应(Lu etal, 2004)，并且相较于现有的微流调查结果(Cubaud and Ho, 2004; Waelchli and von Rohr,2006)。与Cubaud and Ho (2004)的设计一个十分独特的不同是本蛇纹石弯在这里的四舍五入 是被

11、刻意如此做的,目的是来自墙面动量变化降低的惯性作用,尤其是在急转弯处,来逐渐减小对两相流的影响。技术调查(Fries et al., 2008)提到了在拐弯处动量变化对在一个分段气液流的尾流的塞子处液体混合的重要影响，但是仅仅在自然地流动模式中得到印证。微管路得混合部分对两相流的形成具有特别的重要性(Cubaud et al., 2005)。液流被按照T行分路分成两个平行的通路(在长度和体积流动阻力相同下)。在混合区域，两相流将按直角加入的气流围住 (Cubaud et al., 2005; Dollet et al., 2008; Ganan-Calvo and Gordillo, 2001

12、)。如此设计的层流通过一个 200微米深和50微米宽的管口，来促进它的流动。在通过管口的突然扩大200*200微米，将导致产生气泡，具体的见(Dollet et al,2008;Ga nan-Calvo andGordillo,2001)。MaskK*卜 NtgeOvfl PruioreaiaiSilioqr WbterLiqiiid InJelGse JnletrXl imi 1000 VFiji# CcnnXtqrPrsssune Port 1WcaSurwnrnl DomainPn&Mure Pon 2-OuPelIn-*“kitFU. 1 Procetf Adw the btricjA

13、wn of the m ImotfuiifKl chip lllubtr int ptwdlitHff-ifrtiv. ptwbomm develDpincfiE. dry -cthinf b# ICE lucKsidc fMWL diwdK bcndi媲 iOTHl I he lEjcIuncnl ut ibc Muid- Jiid tup vkw J a lypiLdl niuj-ULhjnnd dhqiLSyringe PumpMeOH町0Gas CylinderPlI *IICO；Laser fnducedFluorescence Microscopy! CCD camera I ep

14、i fluorescentI filter cube、i n-iI r - n - i NW: YAGI包汐I Tiicroscope ler-sMicroChannelChip25 barMass Flow Meter DifferentialPressureSsnsofsIra2.5 IbarElectro nic Scafe02ibarData Acquisition SystemFig. 2, Sthtmdlit 血gwE oflHecxpcnnieiitil 沖阿:dtuh2.2实验装置实验室在准实验室条件下（22摄氏度和lOOKpa，实验装置如图.2）进行的。注入泵泵入 需要的水和

15、甲醇溶液 Vi进去，流量从0.0012到1.200mL/min。并对几个不同浓度的甲醇做了 研究。相关参数列于表1,相关流动性质列于表 2。Tdble 2Fluid properties.Fluid mixtureDensityViscosity mPd s)Surface tension (mN/m)H3O997JL001672.751 M MeOH980L09826610MeOH7SG.9.544022.95CO2LSI10.0150-储存于钢瓶(PanGas, purityP99.9 vol%)的CO2的体积流率Vg由质量流量控制器(Bronkhorst,EL-FLOW) 来调节，范围为

16、0.02到1.5mL/min。两个液流都是在 PEEK毛细管中流动的，并 且通过相同的入口注射到微管路芯片中。在混合区域，依据毛细管流动原理的喷嘴主要应用这个项技术(Cubaud et al., 2005; Dollet et al., 2008; Ganan-Calvo andGordillo, 2001),来形成两相流模式。在随后的微管路部分，最终的流动模式是通过一个蛇形结构的管路来研究的。正如前面提到的，弯处被设计的光滑和园面的，是为了降低壁面对动量变化的影响。在喷嘴后的5mm(25*dhyd)处，在水力的入口长度以外，第一个压力点与在研究的最后位置 的第二个压力点是合并的。在这两个点相

17、距精确地70mm。3个不同他的压力传感器(HubaControl 692;0.2 bar, 2.5 bar, 25 bar)被用来记录在两个点之间的压力降与流率的关系。在流过 观察的区域，喷嘴后70mm的下游，两相流流过PEEK的真空管并且收集于一个电子称(MettlerToledo XS204 DeltaRange)上的玻璃烧杯里，为了液体流率的重复试验。2.3荧光镜检法现在有的几个非浸入式的用来描述再微管路中流动模式的方法，比如Brightfield显微镜法，荧光镜检法磁共振成像法等。与平常应用的Brightfield显微镜法相区别，此处的实验应用了一个荧光镜检的技术(LaVision F

18、low Master MITAS),见图2.液体混合物掺杂了 0.01wt%Rhodamine 610 Chloride (Exciton),氟染色也是 Rhodamine B 的 情况下进行的，以得到一个荧光液体和飞荧光气相可分别的标志。用于降低表面张力的浓度一般认为为5%。一个双腔 ND：脉冲波长为10nm的掺钕钇铝石榴石激光器(New Wave,SoloII-15,532 nm)作为光源用于液体荧光检测。Rhodamine 610 Chloride在合适的溶剂里好izai波长为594nm的时候发出荧光。一个向上的荧光滤波器(dichromatic prism)可以预防激光脉冲造成危害；只

19、有发射光到达同步的电耦合摄像机芯片(LaVision Imager ProX 4 M)，一系列500ns的曝光时间形成独特的两相流图像，见图3。(8)(C)Fluoresang LiqukS FilmGas Coref* ManIlltttftfVUIIl仏WWW*“* (dlMcrocTiarviel Wale3流动模式3.1传统型分类两种不同浓度且互不相溶的液体注入一个绝热的微观引线，进过一个特定的进入长度，对于一个依赖于流体性质的特殊流动体系，管路几何学和两种液体的体积流率。对于一个单一气液混合物在特定的管路中的情形，可以得到一个基于气液体积流率流体相图。Table 3Flow patt

20、errtFlow p a ItemMajor flow regimeDisp ersedliiLermiltent (Taylor)Ann LI LitStiaLificcIDispersed bubblyRegularly ordered bubblyWedging (elongated bubbly, plug) Slug flowSmooth jninilarWavy annularSmooth siratifiedWavy stmtified不同的体积流率会产生一系列可区别的状态。试图在以发表的过量的信息中分别对各个状态和它们的描述找到一个共同的标准，本研究依赖于Kreutzer et

21、 al. (2005), Coleman andGarimella (1999) and Cubaud and Ho (2004)的在分类方面的结合的小修改，见表3.这个选择被选来作为一个比较的基础，因为在管路相似的几何特性，这包括矩形断裂面相同的数量级和高宽比。在另外一些文献中的报道(Coleman and Garimella, 1999),一个流动状态的典型不同可以通过圆形管路和矩形断裂面得到。这个额外的复杂性和流动状态的潜在偏差可通过一个尺寸 与另一个尺寸之比的作用描述。两相流的研究已经从第一次在微间距矩形截面切片通道中， 发展到在窄宽度或裂缝形态中研究,这些被导出的额流动相图对矩形断裂

22、面管路的使用性还 有待考证。流动状态的修改时因为除存在分散流之外，还存在有规律的气泡流模式。在不同的环境里，水动力学产生的气泡例如重点的应用于单次小批量反应 谓的有序泡状流被称为气泡悬浮液(Joanicot and Ajdari, 2005)，所(Gordillo et al., 2004)单分散阵列泡沫 (Cubaud et al., 2005)，或者格泡沫(Ganan-Calvo and Gordillo, 2001),但是还并没有合并到任何一个相图分类 中，这取决于做着的知识。既然描述的气泡排列(见下一章节)包括球形滴(见表3b),几何尺寸相同，泡沫一词，描述和关于多面体，改变气体体积，

23、在这里都是不适用的。文献中研究两相流的典型的工作流是水-空气或者水-氮混合物，因此规律地气泡在流动相图的环境中观察不到了。 在下一章节中我们将讨论和描述流动状态。3.2观察到的流动状态的描述 3.2.1分散流状态分散流可以被分为两个部分。在一个低的气体流率开3( a)。气相可以被描述为直径小于管路高度的离散气体流率的提高会使得气泡尺寸的增大，直到直径超正如前面提到的和表 3中列举的， 始，将会产生一个分散的气泡流，见图 球体，而在连续的液相中却是均匀的。 过管路高度。从那一刻起，既然气泡被管路截面所限制，流动状态被认为是间歇的。然而，在特定的条件下，气体流率的增加不会导致气泡尺寸的增大，但是会

24、导致气体数量到达一定的水平，单一气界会按照管路几何尺寸密集的设置于系统规定的参数范围内，见表3(b)。一些作者报道的气泡的合并(Ostwald Ripening)现象，还不能在特殊的流率下得到。本文章不讨论合并的标准。本文提到的是有规律地气泡流流动模式。3.2.2间歇流动状态已知的间歇流动状态是在宽松的不同表现形式，拉长气泡，楔入，塞，塞或分段流下测 得的，然而经常归咎于泰勒和他的工作(Taylor,1960)。这个流动模式的特性是它的气泡直径比管路横截面还大。因为表面张力的主导作用，气泡沿着管路轴线延展而不是分散或者破裂 成球型小液滴。两个可区别的流动模式可以再间歇流动，楔入和段塞流动状态下

25、被检测到，见表3( C)和(d)。虽然两相流的几何外形很相似都可以用一个被细液体膜围绕的轴对称气体堆芯来描 述，并且被与液体活塞分成几段，形成拉长的气泡，楔形流在特定的情况下会出现(Cubaudand Ho, 2004)沿着相邻管路壁的失润和干透情形下。相对的，活塞流模式以一个要更大，因此，加强拉长的气泡为特征，这个气泡没有显现 出任何的干燥区域。这些气体中心将会被高速流动的液体塞段一直润滑着。正如Cubaud and Ho (2004)所描绘的楔形流状态可以依据不同的气泡流速分为两种可以 互相转变的不同润湿状态。如果气泡的速率低于去湿速率，那么细长气泡和管壁之间的液膜会变干。如果气泡速率大于

26、临界速度，液膜将不会破裂，并且气泡将会被液膜包围。3.2.3环形流在低液体流速下，跨过管路横截面的液体活塞将会随着气体速度形成环或波状和最终变 平而消失；一个持续的轴对称的气体中心形成，见表3( e)。液体流仅仅在被管壁覆盖的细液体膜中流动。3.2.4层状流这里并没有观测层状流，但是很快就详尽的阐述了完整性的利处，并且它在一些流动相 图中的出现将会在下面谈给予讨论。层流状态可以完全分开液体和气体两相。在大规模两相流中，这种不同用浮力来感应， 反之对于微观管路规模， 在表面与表面张力的合并中可湿性 将会一直层状流。3.3结论 3.3.1参数集1-4的流动相图我们实验的结果将在这一节中予以讨论，然

27、后与以前报道的研究结果做一个彻底的比较。在图4( a) -( d)显示了四个参数集的结果。参数集1被解释为一个与相关文献中液气混合物比较的案例(Colemanand Garimella, 1999; Cubaud and Ho, 2004; Triplett et al.,1999;Zhao and Bi, 2001)。在参数集1-3中，只有甲醇-水溶液的浓度有改变，见于表2,参数集1 一蒸馏水开始，随后的是参数集2是1M的甲醇溶液，参数集3是纯的甲醇溶液。在参数集 4中，液体浓度与参数集2中一样，但是喷嘴的几何尺寸改变为30微米宽，目的是鉴定喷嘴的尺寸对结果的影响。何甲EM网ADniutar

28、 rOlAlft dfOMDl实验是在一个设定的注射泵控制的特定液体流率并且逐步增加液体流率的条件下进行的。液体流率可以达到稳定的条件，但是在喷嘴出口下游70mm处，还有水压入口长度和压力孔之见的地方，流动模式被记录和确定前至少有300s可以通过。记录的流动模式以表观气体流率 jg为横坐标和表观液体流率ji为纵坐标被列于表中， 见图 4（ a） -（ d）。表观气体流率jg定义为气体体积流率除以管路横截面积 A，表观液体流率ji的 定义是类似的：我们在参数集1-4中可以发现一个有趣的现象，关于起泡，楔入和活塞流之见的转移线。在参数集中表观气体流率 jg=0.08 m/s（the rms of

29、jg=0.06 -.0 m/s）时，转移线T-1和T-2在假象的 流动相图中时以相同德倾角和位置平行的，这个假设是在文中研究的液体性质的改变和喷嘴喷口轻微的改变对两个转移线在起泡和楔形/活塞流转变过程中整体趋势没影响。转移线把转移线T-1:j0.7397* j：.8333楔形流状态和气泡流分开来，精确地表示如下：与活塞流（T-2）平行的楔形流转移线得到一个相似的，关于转移线T-1的指数函数，两者有着相同的指数，但是不同的偏差：转移线T-2:jl =1.086* j：8333在图4 （a）中，在低于临界值jg = 0.08 m/s （只有参数集1和2在jg = 0.08 m/s ）时，气泡流和

30、楔形流之间的传递过程在气体流率中独立出来。这个过渡仅仅依赖于液体体积流率。这个过渡发生在:jl =0.0736m/s与参数集3中相对应，过渡更加依赖于两液体流率，因为甲醇作为一个液体其表面张力远低 于水，在同样的喷嘴尺寸下也就可以形成更小的气泡。这个可以用流动重点技术的工作原理来解释。气体中心的解离形成气泡(或者称泰勒气泡)的情况依赖于惯性力和表面张力的相互作用。表面张力越小，越早发生气泡中心的解离，形成的气泡液越小。这些也可以更方便的用小店的喷嘴来得到，正如参数集4中显示的结果，列于图4 (d)。总的来说，可以这么说两个设计的效果在低于临界值 jg = 0.08 m/g时产生稳定的较小气泡；

31、这些可以要么通过降低 喷嘴的尺寸使其低于50微米或者降低液体的表面张力。对于从楔形流到活塞流(T-2)过渡的一个类似的理由的转移线可以在如下条件进行：ji =0.0159m/s稍微降低1M的甲醇溶液对表面张力的影响可以再参数集2中看出，图4 (b)，但是并没有在参数集1中表现出明显的影响，图4 (a)。然而有趣的是，有规律的气泡状态可以在楔形流状 态中检测到，在转移线 T-1和T-2之间。这些表明甲醇量的很小变化都是使得在喷嘴处气体核 心破碎产生的气泡更小更稳定。同时，气体表面张力又高到足以防止个别气泡的合并，形成有规律的状态。有规律地气泡状态在参数集 4中更深入的检测到， 见图4( d)，这

32、个可以与1M 甲醇的溶液相似的解释。对于参数集3和4过渡装态在高于jg = 0.04 m/s时是十分有趣的。在参数集3的清况中，在溶液中甲醇的混合在从气泡流到楔形流的微管产生一个相变。这个国度完全发生在喷嘴下游70mm的范围内，并且可以用下面的转移线来表征：转移线 T-1a j =27.388j1.9503对于参数集4，一个有规律的气泡流状态作为一个T-1下的转移状态，遵守如下转移线T-1bjl =0.1228j；2778仅仅出现的环形流也在参数集 3和4中检测到，见图4( C)。楔形流和环形流过渡的回归线基 于三个具有很大不确定性的环形流数据点，但是为了完整期间：转移线T-3ji =0.01

33、2j2.78243.4与已有的流动相图的比较在后面的章节，最相关的，最近的和最有影响力的微型两相路研究将与我们的实验数据 做个比较。与实验明显最相似的实践计划由Cubaud and Ho (2004)提出。cTransition Une T_i ； Transition Lino T_2： Transition Line H_l：Transition Line H_2：Transition Line H_3：Iranion Line H_4：bubbfy/wedan。 wedginskig bubbly/wedging we旳ng相u slugtannular annuiarryirregul

34、ar bubbly regularly orcienKl bubbly wedgingslug annutar(CLibaud, and Ho, 2004) (Ctibaud and Ho 20G4i (Cubaud and Ho 2004) (Cubaud and Ho 2004)魅 S_ Row m* cxnvfUTiKn ten QihHid Ita (战Jtrfji3.4.1流动相图与 Cubaud and Ho (2004)的比较Cubaud and Ho (2004).研究了在矩形管路彻底的两相流，管路的横截面长宽分别为 500微米和200微米。微型管路是用反粒子蚀刻技术连接到覆盖

35、着阳性保护玻璃允许光通过的硅上。使用的工作液体是流动聚焦喷嘴混合的去离子水和压缩空气。这个喷嘴被证明在生产大范围尺寸气泡的有效措施。它们遵守喷嘴部分下游的气泡，楔形，活塞，环形，干流动列于典型 的两相流相图，见图5.Cubaud and Ho (2004)基于200微米的管路的转移线流动相图，是画来和我们实验数据做 比较的。在Cubaud and Ho (2004)中，做了这样的假想，在它们衡量领域内的流动相图是独 立于管路横截面尺寸的。在(Cubaud and Ho, 2004)中，做了这样的假设，在它们的衡量领域边界之内的流动相图是与管路横截面的尺寸是相对独立的。在转移线的对500微米管路

36、进行插补法之后，一般性的功能性的依赖于表观速率：其中对于气泡/楔形流转移线(T-1) |c =O.75|，对于楔形/活塞流转移线(T-2) C = 0,20 对于活塞/环状流转移线(T-4)|g = 0.04|，对于环状/干流转移线怎=0,005。与Cubaud and Ho (2004)在200微米横截面的微管得到的数据与我们的实验数据进行比 较，在气泡/楔形流(T-1)和楔形/活塞流(T-2)之间得到很好的一致性，见图5.转移线T-1的偏差也可以观察到。然而，对Cubaud and Ho (2004)实验中的转移线H-1进行深入回归检查后，可以看出这个转移线并不能单独精确获得对于200微米

37、横截面的微观的相转移。他们气泡流数据集中的几个点也并不能很好的符合于提出的气泡流状态。如果对他们数据集提出的相同数据(仅仅200微米时的值)进行分别解释，那会是抓一线会更加的渐渐倾斜，实际上与我 们T-1的结果相同。转移线T-1和T-2的平稳区在我们研究中低于 jg=0.08m/s时出现；这在Cubaud and Ho (2004) 试验中是检测不到得。 原因是我们更大的喷嘴设计， 它后面的孔口也是相同的深度， 可以是 更大的拉长的气泡在更小的气体体积流率下更稳定。在Cubaud and Ho (2004)中描述的活塞/环形流和环形流/干流之间的转化不在我们的测量 范围内，所以，不能得到公平的

38、比较。如果假设一个转移线的指数(线性对数刻度)插补法 到我们的测量范围， 可以估算误差。 一个可能的解释可以从我们活塞流的分类中推断出， 包括一个与几乎不出现流过管路的活塞流环形 -相似的流动状态。(Coleman and与 Cubaud and Ho (2004) 的工作不同，干流装态并没有在其他很多研究 Garimella, 1999; Triplett et al.,1999) 中的20微米管路中观测到。总的来说， 从我们与 Cubaud and Ho (2004) 研究的比较中得到的结论， 可以这么说， 管路 尺寸也对转移线有着一定的影响。在 Cubaud and Ho (2004)

39、中对两个不同管路尺寸的插补法 并没有显示出我们研究发现的平稳状态来支撑下面这个结论，微型规模的两相流不仅仅由管路几何性质 (Coleman and Garimella 1999) 决定，而且还由管路尺寸 (Triplett et al., 1999; Zhao and Bi, 2001)决定。总之，可以这么说 Cubaud and Ho (2004)出的流动相图印证了我们相变的 现象。还需要说明的是我们可以测量两个 Cubaud and Ho (2004)低得多的数据， 拓展了已知的微 型规模两相流的流动状态相图，可以观察到在临界值jl=0.0159m/s以下的活塞流主导现象。3.4.2与Tr

40、iplett et al. (1999) and Zhao and Bi(2001)流动相图的比较与经常归类和引用的文献的比较式必要的，目的是理解和提取对转移线有影响的参数。在下文中， Triplett et al. (1999) and Zhao and Bi(2001) 的工作被用作一个更深入的方法来达到 这个目的，并且与我们的数据进行比较，见图 6.Triplett et al. (1999)在他直径分别为1.1mm和1.45mm横截面为圆形和三角形的高硅管路， 四个流动状态：气泡，活塞，搅拌，环形波状流。设置的方向是水平的，液体，去离子水和 空气在一个十字形的搅拌器中混合，一般伴同通过

41、上面提到的几何尺寸200mm长的管路。Zhao and Bi (2001) 测验了一个由透明合成树脂制造的等边三角形形状的管路。与 Triplett et al. (1999) and Cubaud and Ho (2004) 相比，通过水和空气同时垂直向上的两相流注射通过多 孔介质得到。 流动状态结果是： 毛细管气泡， 活塞， 搅拌和环形流。从直接比较我们的测量数据和 Triplett et al. (1999) and Zhao and Bi (2001) 提出的相变线， 可以明显的观测数据是测量域不同，见图6.为此目的，除了观测领域的转移 (Triplett et al.,1999: j

42、g = 0.02 -80 m/s and jl = 0.02 -8 m/s; Zhao and Bi, 2001: jg = 0.1-00 m/s and jl = 0.08 -m/s)外，还要记住不同的管路几何参数，当从比较中得到结论时。十分有趣的是我们全部的测量范围都在 Triplett et al. (1999)预期的循环槽钢截面的活塞 流流动相图和 Zhao and Bi (2001)的三角微管流动相图内。 除了不同管路几何参数的偏差的说 明外，喷嘴还需要特别的注意。另外在 (Cubaud et al., 2005; Ganan-Calvo and Gordillo, 2001) 中也

43、得到了喷嘴尤其对微管 (Bo w nup CDmiursort to Sim wd &(711)1 i dTriplett r jILT 1999k这个解释支持我们对在比 Triplett et al. (1999)报道要低的j|下的气泡和楔形流。换句话说，考 虑到在前面叙说的偏差系数，在 Zhao and Bi(2001)中检测到的气泡状态更加符合我们的结 果。不过，偶然重合的可能性也不能完全排除。有趣的是，在活塞流和环状流之间相转换时进行转移线的插补，在Tripl ettet al. (1999)和Zhao and Bi (2001)得文献中，对于甲醇溶剂的相转移竟然极其的相似，见图4(

44、c)。与Coleman and Garimella (1999)的比较并没有得到任何的信息；我们的流动相图完全与 其楔形/活塞流流动状态一致，Zhao and Bi (2001)和Triplett et al. (1999)也对其应用做了相似的讨论。3.4.3与Taitel -Dukler模型的比较Taitel and Dukler (1976)在一篇有影响力的文献中提出了一个来描述在水平圆管中的宏 观两相流中的相转移物理模型。这个模型基于一个对切向力指数定律摩擦系数的假设 (Agrawal et al., 1973)两相层流中每一相(分别是气相和液相)的动量平衡。一个动量平衡无因次的表征为X

45、的函数(Lockhart and Martinelli, 1949),无因次参数Y反应的是倾斜的影响。倾角为0时，正如在本文中研究和比较的，Lockhart -Martinelli-Coeficient X仅仅与液相因次的平衡高度hl =h / D有关。Taitel and Dukler (1976)通过一个特殊无因次数发展出了每一个相转移的标准。间歇层流的相转移是基于Kelvin -Helmholtz为了发展在液/气界面的极小波的稳定性判其中C2是无穷小波的一个干扰因数;其中F是修正的弗劳德数:ug,h ,A,厲 分别是无因此气体速度，液体装填高度,液体和气体中心横截面。所有的变量都会被管路

46、几何约束规范化。间歇到环状流状态过渡的准则是基于Butterworth (1972)的结果，该结果显示一波足够高度掠过管壁形成在初次蓄水高度低于hl =0-5时的环状流装态。结果显示，对于水平管，X为常数值X=1.6。Taitel and Dukler (1976)提出的对从波动过渡到层状光滑流的准则是局限的或者说与小 规模两相流没有十分密切的关系。换句话说，从气泡流到间歇自喷流状态的转变时相关的， 并且是基于他们在服从于下面无因次准则气泡的流动和浮力稳流脉动的平衡中工作得到的 结论：U| , Ag , D| , S列阴)J其中分别为无因次液体流率，横截面积，水力直径和气液接触面积。所有的变量

47、都 由液体装填高度和管路直径规范化。即使这个模型被证明可以很适当地推测在宏观管路中的相转移，就如Taitel andDukler(1976)为Mandhane et al. (1974)研究的一个直径25mm的圆管做的流动相图中显示的， Taitel -Dukler模型没能成功推测在更小直径管路中的流动相转移，比如Damiandides andWestwater (1988)和 Coleman and Garimella (1999)比较的直径小于 5mm时的情况。Transition Lino 匸匕 Transition Lirl_: Transition Line F:Transition

48、 Line K: TransitMjfi LineTbubblyvedging dging/ug stralified/i ntermidfirt inlonTTtljenvannular bubblylntBrTnitterit0 irwgularbubb regularly oncSerdd bub My 口 wodging aluQ* annutar(Taiisl and Dukler 1076)laiiel and Dukler 1976) (Taial and DuWef 1376)1. 7. Fbw nup ampuufliDu Uh rdelTaitel -Dukler模型（Ta

49、itel and Dukler,1976）与我们实验的参数集 1的比，见图7,转移线F（层状流和间歇流之间的相变，基于无量纲准则F）的产量不一致。同时，需要指出的是参数集1没有描述任何的层流；也没有光滑流或者波状流。拓展了的Kelvin -Helmholtz理论的应用不能够圆满的接受再微流体中波动的增强。换句换说，在以Taitel -Dukler为基础的转移线吧间歇式流动从环状流（转移线X，基于X=1.6 ）中完全的排除，显示出了与我们实验中得到的从楔形流向活塞流过渡相似的坡度。表面上的一致能显示出关于模型有效性的一点有趣的信息：即使管路几何特性可能导致一些误差，但是基本的动量平衡状态还是有效

50、地。间歇流到分散气泡流的转移线（转移线T-1，基于无量纲准则F ）仅仅与我们在将较小的表观气体速率与转移线平衡匹配时才与我们的数据一致。但是，这个一致被认为是纯粹的巧合，因为既没有基于重力的浮力也没有不稳定的起伏现象在完全忽视便面张力的重要效果的 微流体中起重要的作用。关于原Taitel -Dukler模型的结束语市这样的，在模型中缺少对表面张力效应的研究，导 致它缺少了在低于临界邦德数Bo=3.368是的微流体流动的应用。3.4.4与修改了的Taitel Dukler模型的比较Barnea et al. （1983）做了一个工作，他将 Taitel -Dukler模型对证实了的直径 4mm或者更大 的较小直径管路的分类。表面张力作为一个相关作用力被考虑了进来，在层状流到间歇自喷流中遵循下面准则：图8显示显示了这个准则与参数集3的比较。可以发现