Y型微通道两相流内部流动特性

利用显微粒子图像测速技术、高速度数码显微系统及数值模拟方法研究了Y 型微通道内液滴的形成. 主要考虑了Y 型角度（45°,90°,135°,180°）、两相流量大小等因素的影响. 发现在挤压机制中,Y 型微通道内分散相液滴的形成主要受到来自连续相的剪切作用,Y 型角度越小,分散相所受到的剪切作用越大. 在液滴生成过程中,连续相速度剖面呈非对称抛物线型分布. 当Y 型角度小于180°时,角度的变化对液滴直径大小影响较小,但角度的减小会加快液滴的生成时间. 当Y 型角度为180°时,生成的液滴体积最大且生成时间最长. 毛细数对液滴直径和生成时间的变化同时产生影响,连续相毛细数的增大使得连续相在两相交汇位置处对分散相的作用力更集中,导致分散相更易破裂.      

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

流动聚焦型微流控装置能够方便、高效地生成均一度好且大小精确可调的微液滴(气泡),故被广泛应用于颗粒材料合成、药物封装、细胞培养等诸多领域. 进一步优化通道结构有助于实现对合成微粒粒径、均一度和尺寸范围的精确调控. 本文数值研究了通道深度、缩颈段长度以及两相夹角等几何构型因素对流动聚焦生成微液滴直径及其生成周期各个阶段的影响. 控制液滴生成方式为滴流式,发现液滴直径随通道深度d 的增加近似呈线性增大,且当通道深度小于30 μm 时,随着通道深度的下降,微液滴生成周期在毛细力的强烈作用下出现骤升,通道深度超过80 μm 时,微液滴的生成周期基本接近恒定. 连续相和离散相的夹角θ接近90°时,液滴直径及其生成周期最短,夹角太大或太小均不利于生成均一度好且粒径微小可控的液滴. 调整缩颈段长度l引起液滴直径及其生成周期的变化幅度仅为其平均值的3%~5% 左右. 此外,缩颈段宽度也是影响流动聚焦生成微液滴直径及其生成周期的重要因素,在通道深度固定时,缩颈段越宽,微液滴直径及其生成周期越大.      

组装式微流控系统制备双重包裹微液滴的方法研究

介绍了一种组装式微流控系统制备单、双重包裹微液滴的方法。微系统中用三通接头构成T型微流体通道,使得分散相在连续相强烈的剪切力和压力差作用下断裂形成单个微米级液滴。在制备单个微液滴基础上,用毛细管将两个三通接头串联,通过调控三相流量,可产生双重包裹液滴。结合实验结果,分析了流体粘度比对液滴大小的影响,并得出液滴的尺寸与流量比之间的关系式,为制备不同尺寸的液滴提供了参考依据。对制备的样品进行统计分析,结果显示,液滴的多分散性指数均小于3.2%,表明微液滴的高度均匀性。此外,通过调节三相液体的流量不仅可以控制内外层液滴的大小,还可以调节内层包裹液滴的个数。本文提出的制备方法,设备组装拆卸简便,不需表面亲疏水性处理,装置利用率高,产生的单、双重包裹微液滴可满足高通量的测量分析要求。     

T型微通道反应器内气液两相流动机制及影响因素

基于液滴或气泡的多相微流控是近年来微流控技术中快速发展的重要分支之一.本文利用高速显微摄影技术和数字图像处理技术对T型微通道反应器内气液两相流动机制及影响因素进行实验研究.实验采用添加表面活性剂的海藻酸钠水溶液作为液相,空气作为气相.研究T型微通道反应器内气液两相流型的转变过程,并根据微通道内气泡的生成频率和生成气泡的长径比对气泡流进行分类.研究发现当前的进料方式下,可以观测到气泡流和分层流2种流型,且依据气泡生成频率和微通道内气泡的长径比可将气泡流划分为分散气泡流、短弹状气泡流和长弹状气泡流3种类型,并基于受力分析确定3种气泡流的形成机制分别为剪切机制、剪切–挤压机制和挤压机制.考察不同液相黏度和表面张力系数对不同类型气泡流范围的影响规律.结果表明:液相黏度相较于表面张力系数而言,对气泡流生成范围影响更大.给出不同类型气泡流流型转变条件的无量纲关系式,实现微通道生成微气泡过程的可控操作.     

T型微通道反应器内气液两相流动机制及影响因素

基于液滴或气泡的多相微流控是近年来微流控技术中快速发展的重要分支之一.本文利用高速显微摄影技术和数字图像处理技术对T型微通道反应器内气液两相流动机制及影响因素进行实验研究.实验采用添加表面活性剂的海藻酸钠水溶液作为液相,空气作为气相.研究T型微通道反应器内气液两相流型的转变过程,并根据微通道内气泡的生成频率和生成气泡的长径比对气泡流进行分类.研究发现当前的进料方式下,可以观测到气泡流和分层流2种流型,且依据气泡生成频率和微通道内气泡的长径比可将气泡流划分为分散气泡流、短弹状气泡流和长弹状气泡流3种类型,并基于受力分析确定3种气泡流的形成机制分别为剪切机制、剪切-挤压机制和挤压机制.考察不同液相黏度和表面张力系数对不同类型气泡流范围的影响规律.结果表明:液相黏度相较于表面张力系数而言,对气泡流生成范围影响更大.给出不同类型气泡流流型转变条件的无量纲关系式,实现微通道生成微气泡过程的可控操作.      

流动聚焦装置中液滴的形成:喉部尺寸的影响

在本文中,研究了液滴在一个流动聚焦微流体设备中的形成过程,分析了喉部长度和宽度以及连续相的流速和分散相的粘度对液滴尺寸的影响。在固定的分散相流速(Qd)下,连续相流速(Q_c)对于液滴尺寸有重要的影响。当Q_c0.7mL/h,液滴尺寸在喉部长度到达一个临界值之前先趋于增加,之后随着喉部尺寸的继续增加逐渐下降;当Q_c0.7mL/h,液滴尺寸随着喉部长度的增加而降低。而越大的喉部宽度会产生越大尺寸的液滴。在Q_c继续增长的过程中通常会出现从挤压模式到滴模式的转变,最终液滴尺寸呈现出随着Q_c的增加指数降低的特征。归因于流速控制破碎机制,低粘度分散相下,液滴尺寸随粘度的增加而增加。     

NUMERICAL STUDY ON THE SEDIMENTED MOTION CHARACTERISTICS OF THREE ALIGNED CIRCULAR PARTICLES IN THE INCLINED CHANNELS

主要应用浸没边界的格子玻尔兹曼方法（immersed boundary-lattice Boltzmann method, IB–LBM） 对处于不同倾斜角度通道内的三个刚体圆形颗粒在重力作用下下落的动力学特性进行了计算研究. 首先分析通道倾斜角度的影响, 结果显示当通道倾斜角处于59°90°的范围时会发生后一个颗粒超越前一个颗粒的现象. 其次, 研究了Re对颗粒沉降特性的影响, 结果表明Re 越大, 颗粒间发生聚集的时间越早. 研究还发现当3 个颗粒的直径大小不均匀时, 颗粒由大到小纵向依次排列, 或者出现中间小球直径较相邻两个小球直径大的排列情况, 均能促使颗粒加快聚集. 本文的研究结果可为环境工程及地质学中的颗粒沉降问题提供有价值的参考.     

压电式微滴按需喷射的过程控制和规律

微滴喷射增材制造技术中沉积微滴的大小与均匀性是影响成型件质量的关键因素.本文设计了一种用于生成均匀微滴的压电驱动式微滴喷射装置,通过压电材料带动柔性膜片振动,将液体从喷嘴中喷出生成微滴,采用数值模拟和实验相结合的方法,研究了不同控制参数下膜片振幅及其对生成微滴尺寸和均匀性的影响.研究结果表明:膜片振幅大小受到驱动电压和压电频率的共同影响,压电频率是导致膜片中心点振幅实验测量值小于理论计算值的主要原因,膜片振动会导致喷嘴内部压力发生变化从而影响微滴生成尺寸.在相同驱动电压条件下,压电频率为10 Hz时存在压电膜片振幅最大值.随着膜片振幅的增大,喷孔处液体速度和液柱长度增大到临界值时可以生成微滴,当喷孔处的液柱长度超过临界值时,会形成卫星液滴. 当膜片振幅区间在30 $\mu$m$\sim $42 $\mu $m可以稳定生成微滴,生成最小微滴尺寸为339.8$\mu$m,直径最大变化率为0.29%,相邻两微滴间距最大变化率为2.67%,生成微滴的尺寸及均匀性较好.研究结果有助于提高压电式微滴喷射装置的液滴生成质量.      

压电式微滴按需喷射的过程控制和规律

微滴喷射增材制造技术中沉积微滴的大小与均匀性是影响成型件质量的关键因素.本文设计了一种用于生成均匀微滴的压电驱动式微滴喷射装置,通过压电材料带动柔性膜片振动,将液体从喷嘴中喷出生成微滴,采用数值模拟和实验相结合的方法,研究了不同控制参数下膜片振幅及其对生成微滴尺寸和均匀性的影响.研究结果表明:膜片振幅大小受到驱动电压和压电频率的共同影响,压电频率是导致膜片中心点振幅实验测量值小于理论计算值的主要原因,膜片振动会导致喷嘴内部压力发生变化从而影响微滴生成尺寸.在相同驱动电压条件下,压电频率为10 Hz时存在压电膜片振幅最大值.随着膜片振幅的增大,喷孔处液体速度和液柱长度增大到临界值时可以生成微滴,当喷孔处的液柱长度超过临界值时,会形成卫星液滴.当膜片振幅区间在30μm～42μm可以稳定生成微滴,生成最小微滴尺寸为339.8μm,直径最大变化率为0.29%,相邻两微滴间距最大变化率为2.67%,生成微滴的尺寸及均匀性较好.研究结果有助于提高压电式微滴喷射装置的液滴生成质量.     

液滴冲击液面过程中变形特征的实验研究

用高速相机拍摄了液滴以不同速度冲击液面的运动过程,测量了冲击过程中典型时刻液坑、中心液柱、次生液滴和次生液柱的几何尺寸,对这些几何尺寸的最大值与韦伯数(We)的关系进行了回归分析,结果表明:空间上,液坑的最大垂直深度、液坑的最大水平长度、中心液柱的最大高度、次生液滴等效直径和次生液柱的最大高度随We数增加呈线性增加;次生液滴等效直径是初始液滴等效直径的1.2～2倍;当200We220时,次生液滴冲击液面没有次生液柱生成,当360We713时,有次生液柱形成,220≤We≤360为过渡区;时间上,We数越大,运动过程中液坑的大小和中心液柱的高度变化越快,液坑和中心液柱持续时间越长,中心液柱达到最大高度的时刻以及破碎生成次生液滴的时刻也越晚。     

Droplet coalescence in the shear flow of model emulsions

During the flow of an emulsion, droplets of the dispersed phase can deform, break up, coalesce or migrate to other regions within the flow field. Understanding these different processes is relevant to morphology development in immiscible polymer blends. Here, emulsions of castor oil in silicone oil were employed to study shear-induced coalescence alone; the conditions chosen were such that drop breakup and drop migration did not occur. A cone-and-plate device and tubes of varying length were used to examine the influence of the average shear rate, the time of shearing, concentration of the dispersed phase, and temperature on the average droplet size. It was found that the extent of “demixing” was not influenced by the spatially non-homogeneous nature of flow in a tube; results correlated very well with the average shear rate. On the other hand, coalescence was significant even when the concentration of the dispersed phase was as low as 0.5%, and it became more important as the concentration was increased. Other results were that the extent of coalescence could be promoted by lowering the shear rate. In quantitative terms, it was found that available coalescence theory gave the correct order of magnitude for the average steady-state droplet size as a function of the imposed shear rate, but the actual variation of drop size with shear rate was gentler than that predicted by theory. An unusual observation was that, under some circumstances, the droplets did not coalesce but simply stuck to each other and maintained their separate identity. Received: 25 March 1999/Accepted: 22 July 1999     

水击驻波场中分散相颗粒的运动分析

根据水击在管段内形成的驻波场现象,分析了流体内分散相颗粒受到的驻波作用力;运用李雅普诺夫稳定判据研究了颗粒积聚与分离的机理;考虑到颗粒运动方程的严重刚性而很难进行数值求解,采用相空间和非对称分析方法获得了分散相颗粒的运动轨迹近似解,并进行了实验验证。结果表明:水击驻波场中分散相颗粒的受力方程中惯性项对颗粒初始运动速率的影响不可忽略;在水击驻波波节的±λ/4范围内,分散相颗粒经过一定的时间会发生积聚,其运动速度呈对称分布,最大速度出现在3λ/8位置处;随着分散相颗粒粒径和密度等物性参数以及水击驻波的频率和连续相初始速度的增大,颗粒达到平衡位置的时间呈减小趋势,且连续相的初始速度对颗粒到达波节时间的影响显著。     

Study of two-phase flows in reduced gravity using ground based experiments

Experimental studies have been carried out to support the development of a framework of the two-fluid model along with an interfacial area transport equation applicable to reduced gravity two-phase flows. The experimental study simulates the reduced gravity condition in ground based facilities by using two immiscible liquids of similar density namely, water as the continuous phase and Therminol 59^® as the dispersed phase. We have acquired a total of eleven data sets in the bubbly flow and bubbly to slug flow transition regimes. These flow conditions have area-averaged void (volume) fractions ranging from 3 to 30% and channel Reynolds number for the continuous phase between 2,900 and 8,800. Flow visualization has been performed and a flow regime map developed which is compared with relevant bubbly to slug flow regime transition criteria. The comparison shows that the transition boundary is well predicted by the criterion based on critical void fraction. The value of the critical void fraction at transition was experimentally determined to be approximately 25%. In addition, important two-phase flow local parameters, including the void fraction, interfacial area concentration, droplet number frequency and droplet velocity, have been acquired at two axial locations using state-of-the-art multi-sensor conductivity probe. The radial profiles and axial development of the two-phase flow parameters show that the coalescence mechanism is enhanced by either increasing the continuous or dispersed phase Reynolds number. Evidence of turbulence induced particle interaction mechanism is highlighted. The data presented in this paper clearly show the marked differences in terms of bubble (droplet) size, phase distribution and phase interaction in two-phase flow between normal and reduced gravity conditions.     

CHARACTERISTIC OF DROPLET OSCILLATION ON THE SURFACE OF RECTANGULAR HYDROPHOBIC GROOVES

液滴振荡行为是液滴运动中的重要伴随现象,具有重要科研价值.由于液滴撞击疏水沟槽板时运动行为与光滑表面明显不同,可以推测疏水沟槽表面液滴振荡特性也将会呈现与众不同的行为特点.采用高速摄像技术,研究了矩形疏水沟槽表面上水滴高度和接触线振荡行为随沟槽尺寸和撞击速度的变化规律.结果发现,矩形疏水沟槽造成的各向润湿异性使得振荡过程中水滴在平行沟槽方向上的接触线长度大于垂直方向,但并不影响水滴高度方向上衰减振荡的周期,即水滴振荡周期与沟槽间距无关;同时由于疏水沟槽表面上存在能垒束缚效应,致使水滴振荡过程中接触线的铺展和回缩运动不服从典型阻尼振荡规律,而呈现振荡数次后直接趋稳的特点.如水滴以0.61 m/s撞击时,接触线经历2次振荡后即维持稳定,但此时水滴仍在持续振荡中.另外,还初步分析了水滴振荡周期与沟槽间距无关的原因.     

A CLOSED SYSTEM OF EQUATIONS FOR DENSE TWO-PHASE FLOW AND EXPRESSIONS OF SHEARING STRESS OF DISPERSED PHA’E AT A WALL

In this paper, each of the two phases in dense two-phase flow is considered as continuous medium and the fundamental equations for two-phase flow are described in Eulerian form. The generalized constitutive relation of the Bingham fluid is applied to the dispersed phase with the analysis of physical mechanism of dense two-phase flow. The shearing stress of dispersed phase at a wall is used to give a boundary condition. Then a mathematical model for dense two-phase flow is obtained. In addition, the expressions of shearing stress of dispersed phase at a wall is derived according to the fundamental model of the factional collision between dispersed-phase particles and the wall.     

两串列扑翼的相位差对平均推力影响机理的实验研究

在一个低雷诺数的循环水洞中,实验研究了前后翅翼之间的相位差对两串列扑翼平均推力的影响.利用一个三分量的Kistler 压力传感器来测量扑翼的瞬时力;利用一个数字粒子测速仪系统(TSI DPIV) 来测量扑翼的前缘涡以及其周围的流场. 当相位差从0° 增加到360°,前翅的平均推力随着相位差正弦变化;前翅平均推力的增加是由于后翅的前缘涡和滞止区域增加了前翅的有效攻角. 后翅平均推力曲线有一个明显的V 字形低谷.低谷处较小的平均推力是由于前翅的脱落涡抑制了后翅前缘涡的形成并且减小了其有效攻角.当间距为0.5倍弦长相位差约为290°时,前后翅翼平均推力系数的合值能达到最大值0.667,明显大于两倍的单翼平均推力系数(2×0.255).      

Flame Structure and Propagation in Turbulent Flame-Droplet Interaction: A Direct Numerical Simulation Analysis

Three-dimensional Direct Numerical Simulations (DNS) in canonical configuration have been employed to study the combustion of mono-disperse droplet-mist under turbulent flow conditions. A parametric study has been performed for a range of values of droplet equivalence ratio ?_d, droplet diameter a_d and root-mean-square value of turbulent velocity u^′. The fuel is supplied entirely in liquid phase such that the evaporation of the droplets gives rise to gaseous fuel which then facilitates flame propagation into the droplet-mist. The combustion process in gaseous phase takes place predominantly in fuel-lean mode even for ?_d>1. The probability of finding fuel-lean mixture increases with increasing initial droplet diameter because of slower evaporation of larger droplets. The chemical reaction is found to take place under both premixed and non-premixed modes of combustion: the premixed mode ocurring mainly under fuel-lean conditions and the non-premixed mode under stoichiometric or fuel-rich conditions. The prevalence of premixed combustion was seen to decrease with increasing droplet size. Furthermore, droplet-fuelled turbulent flames have been found to be thicker than the corresponding turbulent stoichiometric premixed flames and this thickening increases with increasing droplet diameter. The flame thickening in droplet cases has been explained in terms of normal strain rate induced by fluid motion and due to flame normal propagation arising from different components of displacement speed. The statistical behaviours of the effective normal strain rate and flame stretching have been analysed in detail and detailed physical explanations have been provided for the observed behaviour. It has been found that the droplet cases show higher probability of finding positive effective normal strain rate (i.e. combined contribution of fluid motion and flame propagation), and negative values of stretch rate than in the stoichiometric premixed flame under similar flow conditions, which are responsible for higher flame thickness and smaller flame area generation in droplet cases.