矩形微通道中气–非牛顿流体两相流实验研究

本文通过实验研究了空气与非牛顿流体在小尺寸矩形微通道中的两相流动特性,给出了在不同气液两相流速下的流型图。同时比较了气液两相流速,液相黏度,表面张力,微通道尺寸对各流型分布区域以及Taylor气泡/液柱长度的影响。发现气液两相流流速变化对流型区域的影响极为明显,而黏性力和表面张力只是在局部范围改变了流型分布。Taylor气泡长度随气液两相流速比增大而增大,随液相黏度及表面张力增大而减小。液柱长度随液相黏度增大而增大,随两相流速比及表面张力增大而减小。最后,我们给出了基于无量纲数JG/JL,Re和Ca的Taylor气泡/液柱长度预测公式,公式预测结果与实验结果基本一致。     

集输管路上升管系统严重段塞流实验和理论模拟

严重段塞流的实验研究表明,在气泡进入上升管底部到运动至出口的过程中,上升管中气泡头部以下流型为弹状流型;当气泡头部流出上升管后,上升管中的流型可看作块状流型。根据实验结果,本文提出了采用漂移流模型简化计算上升管中两相流动、上游管道中气体膨胀满足质量守恒,同时考虑上升管内液体动量守恒的严重段塞流计算模型。计算值与测量值比较表明,模型可以正确预测出气体膨胀流动过程,气体流动时间不受入口气液流量的影响。模型可以准确计算出严重段塞流周期、液塞长度和倾斜管中液柱最大长度等参数。     

微推进器两相流型特征及产生机理

针对以低浓度过氧化氢作为液体燃料的硅基微推进器,采用高速可视化手段考察过氧化氢分解形成的气液两相流在微喷嘴扩张段的流动特征.发现流型总体特点是泡状流型阶段和环状流型喷气阶段交替出现.在环状流喷气阶段,微喷嘴扩张段内的气液界面呈现特殊的瓶颈形结构,形成瓶颈形流型结构的微液滴在微喷嘴出口处往复运动,形成单向止回阀效应,使得环状流喷气过程呈现周期性间断的特点.     

T型微通道气液两相流流动与传热的数值研究

为进行微通道气液两相流动及传热研究,建立了T型微通道模型,利用VOF方法进行数值模拟,对弹状流流型进行了数值分析,考察了弹状流的形成过程、速度、表面张力和粘度对弹状流流型的影响,计算了气泡速度和液膜厚度,并与经验关联式进行比较。结果表明:液体粘度对弹状流的影响较大;而表面张力对弹状流几乎没有影响,当气速一定时,液速越大,气泡的长度越小而液柱的长度越大。气泡速度和液膜厚度与经验关联式的计算结果误差分别为2.57%和4.87%。     

对流梯形微通道内气液两相流动的数值研究

采用VOF方法,对梯形微通道内不可压缩气液两相流动进行了数值模拟研究,详细分析了气泡形成过程,以及当量直径、截面形状、液体表面张力和粘度等对气泡液柱形成过程和长度的影响,拟合出微通道气泡液柱长度计算公式。结果表明:气泡液柱的长度受表观气速和表观液速的影响较大;表面张力对气泡尺寸的影响较小,当液体粘度增加为水粘度的10倍时,形成的气泡形状不规则。增大表面张力,形成气泡的时间增加;增大粘度,形成气泡的时间减小。     

矩形微通道内两相流流阻压降特性的可视化研究

以去离子水和质量分数为0.3%的水基纳米流体为工质,在当量直径为1.241mm的矩形微通道内进行了两相流流动沸腾的实验研究,并借助高速摄像仪对矩形微通道内流型随着质量流量及热流密度的变化进行了观察分析。实验结果表明:单位长度上的两相摩擦压降会随着质量流速的提升而提高;单位长度上的两相摩擦压降会随着热流密度的增大而升高;减小质量流速和提高其热流密度均会加快气泡的产生并提高气泡的脱离直径,当热流密度增大到一定程度时,通道内几乎为环状流与液态单相流交替出现,且环状流占周期中的较长时间。     

复杂微通道内气泡在浮力作用下上升行为的格子Boltzmann方法模拟

本文采用气-液两相流格子Boltzmann方法模拟了复杂微通道内气泡在浮力作用下的上升过程,主要研究障碍物表面润湿性、浮力大小、障碍物尺寸和气泡初始位置对气泡变形、分裂、合并的动力学行为以及对气泡上升速度、终端速度和气泡剩余质量的运动特性的影响.研究发现,障碍物表面接触角较小时气泡能够完整地通过障碍物通道,随着障碍物表面接触角增加,气泡通过障碍物通道时严重变形,并会发生分裂行为,使得部分气泡黏附在障碍物表面,从而导致气泡到达终端时质量减少.相应地,气泡上升速度以及终端速度也随着微通道表面接触角的增加而减小.另一方面,随着浮力的增加,气泡在上升过程中更容易发生分裂和合并现象,且气泡剩余质量和终端速度随着浮力的增加呈对数形式增加.此外,随着微通道障碍物半径增加,气泡剩余质量首先缓慢减小然后快速减小,而气泡终端速度近似呈线性减小.最后,数值结果还表明当气泡初始位置偏离管道中间时,其上升速度、气泡剩余质量以及气泡终端速度都与初始位置在管道中间时的变化趋势一致,然而对应的数值均减小,且气泡在上升过程中变形更严重.     

同向流动微气泡产生与分裂的实验研究

本文使用高速CCD可视化技术研究硅微通道内气液等温同向流动微气泡形成过程.作为离散相的空气和连续相的水,流经各自通道后汇聚在主通道前段;在一定的流量比值范围内,连续产生大小均匀的微气泡;发现存在bubbling和jetting两种微气泡形成模式.气泡产生频率随空气和水的流速增大而增加.最后研究了气泡在T型微通道出口处的行为,观察到分裂和不分裂两种形式并进行了分析.     

流动聚焦型微流控体系中的液滴的图案化排列和转变模式

流体在微流通道中形成剪切流场(低雷诺数)．不同于宏观体系,由于剪切力和表面张力的竞争作用,产生的液滴在微尺度下的微流通道中形成特殊的排列现象---周期性类似“晶格”排列现象．设计了新型流动聚焦型微流控芯片,分析研究在微流体系中液滴周期性图案化排列和转变机理性,液滴排列模式受两方面因素影响：水油两相的流速比值和微通道尺寸．当微通道宽度为250或300 μm时,液滴形成单层分散,双层和单层挤压排列．当微通道宽度为350 μm 时,液滴会形成单层分散到三层排列到双层挤压最后到单层挤压排列．当出口通道宽度增加到400 μm时,甚至出现了液滴四层排列的现象．同时研究了各个液滴排列模式的“转变点”.     

Y型微通道内双重乳液流动破裂机理

基于体积分数法建立了Y型微通道中双重乳液流动非稳态理论模型,数值模拟研究了Y型微通道内双重乳液破裂情况,详细分析了双重乳液流经Y型微通道时的流场信息以及双重乳液形变参数演化特性,定量地给出了双重乳液流动破裂的驱动以及阻碍作用,揭示了双重乳液破裂流型的内在机理.研究结果表明:流经Y型微通道时,双重乳液受上游压力驱动产生形变,形变过程中乳液两端界面张力差阻碍双重乳液形变破裂,两者正相关;隧道的出现将减缓双重乳液外液滴颈部收缩速率以及沿流向拉伸的速率,并减缓了内液滴沿流向拉伸的速率,其对于内液滴颈部收缩速率影响不大;隧道破裂和不破裂工况临界线可以采用幂律关系式l~*=βCa~b进行预测,隧道破裂和阻塞破裂工况临界线可以采用线性关系l~*=α描述;与单乳液运动相图相比,双重乳液运动相图各工况的分界线关系式系数α和β均相应增大.     

液滴碰撞Janus颗粒球表面的行为特征

为研究液滴碰撞Janus颗粒(双亲性)球表面的独特行为特征,以粒径为5.0 mm铜球为材料制备了Janus颗粒,用直径为2.0 mm的液滴,在韦伯数(We)为2.7,10,20,30的测试情况下对Janus颗粒球表面进行了碰撞实验.结果表明:液滴碰撞Janus颗粒球表面后的运动可分为铺展、回缩、振荡和回弹4个过程.在不同We下,液滴碰撞Janus颗粒后的运动状态主要与表面润湿性相关,在Janus颗粒亲水侧表现为铺展特性且铺展系数γ随着时间t的增大而逐渐增大并趋于稳定;但在疏水侧,表现为回弹现象,铺展系数γ会出现类似"抛物线"形状;当液滴碰撞Janus颗粒球表面亲-疏水分界线时,液滴铺展和回弹同时发生.基于能量平衡和受力分析发现,液滴动能和表面能的互相转化是液滴铺展的关键,液滴会在重力、惯性力、表面张力、黏性力、接触力等力的综合作用下展现其独特的行为特征并最终达到平衡状态.     

基于范德瓦尔斯表面张力模式液滴撞击疏水壁面过程的研究

液滴撞击疏水壁面过程的研究在介观流体力学和微流体作用材料科学的研究中具有重要的理论意义和工程价值. 论文在SPH方法中引入范德瓦尔斯状态方程处理液滴表面张力, 考虑流体粒子之间远程吸引, 近程排斥的内部作用力, 提出了流体粒子与疏水壁面粒子间势能函数与表面张力相结合的作用模式. 通过模拟真空条件下两个静止的等体积液滴相互融合的过程, 验证了计算模式在模拟液滴的表面张力中的有效性. 采用该模式模拟的液滴撞击疏水壁面过程, 不仅能够有效地模拟液滴撞击壁面后的变形过程, 而且清晰地模拟出液滴的回弹、腾空以及二次撞壁现象的完整过程. 模拟结果与液滴撞击疏水壁面的实验结果以及VOF模拟结果符合较好, 表明本文所提出的表面张力和疏水壁面作用力处理模式对模拟液滴撞壁过程具有实际应用价值.     

基于VOSET方法模拟并排气泡的上升过程

采用VOSET方法和耦合表面张力模型的N-S方程,模拟竖直通道内并排气泡对的上升过程,模拟与实验结果吻合较好.重点研究表面张力系数对并排气泡上升轨迹和速度的影响.结果表明,随着表面张力系数的变化,并排气泡对的上升过程出现三种类型:两气泡融合,两气泡反复靠近、远离但未融合,两气泡碰撞反弹后逐渐远离.在未融合的情况下,并排气泡对的上升轨迹关于通道中心线对称,左右两个气泡的上升速度基本一致,水平速度大小相同,方向相反.     

Numerical simulation on dynamic behaviors of bubbles flowing through bifurcate T-junction in microfluidic device

Based on the volume of fluid(VOF) method, a numerical model of bubbles splitting in a microfluidic device with T-junction is developed and solved numerically. Various flow patterns are distinguished and the effects of bubble length,capillary number, and diameter ratio between the mother channel and branch are discussed. The break-up mechanism is explored in particular. The results indicate that the behaviors of the bubbles can be classified into two categories: break-up and non-break. Under the condition of slug flowing, the branches are obstructed by the bubbles that the pressure difference drives the bubbles into break-up state, while the bubbles that retain non-break state flow into an arbitrary branch under bubbling flow condition. The break-up of the short bubbles only occurs when the viscous force from the continuous phase overcomes the interfacial tension. The behavior of the bubbles transits from non-break to break-up with the increase of capillary number. In addition, the increasing of the diameter ratio is beneficial to the symmetrical break-up of the bubbles.     

Hydrodynamic binary coalescence of droplets under air flow in a hydrophobic microchannel

Based on the volume of fluid(VOF) method, we conduct a numerical simulation to study the hydrodynamic binary coalescence of droplets under air flow in a hydrophobic rectangular microchannel. Two distinct regimes, coalescence followed by sliding motion and that followed by detaching motion, are identified and discussed. Additionally, the detailed hydrodynamic information behind the binary coalescence is provided, based on which a dynamic mechanical analysis is conducted to reveal the hydrodynamic mechanisms underlying these two regimes. The simulation results indicate that the sliding motion of droplets is driven by the drag force and restrained by the adhesion force induced by the interfacial tension along the main flow direction. The detachment(i.e., upward motion) of the droplet is driven by the lift force associated with an aerodynamic lifting pressure difference imposed on the coalescent droplet, and also restrained by the adhesion force perpendicular to the main flow direction. Especially, the lift force is mainly induced by an aerodynamic lifting pressure difference imposed on the coalescent droplet. Two typical regimes can be quantitatively recognized by a regime diagram depending on Re and We. The higher Re and We respectively lead to relatively larger lift forces and smaller adhesion forces acting on the droplet, both of which are helpful to detachment of the coalesced droplet.     

Decompressing emulsion droplets favors coalescence

The destabilization process of an emulsion under flow is investigated in a microfluidic device. The experimental approach enables us to generate a periodic train of droplet pairs, and thus to isolate and analyze the basic step of the destabilization, namely, the coalescence of two droplets which collide. We demonstrate a counterintuitive phenomenon: coalescence occurs during the separation phase and not during the impact. Separation induces the formation of two facing nipples in the contact area that hastens the connection of the interfaces prior to fusion. Moreover, droplet pairs initially stabilized by surfactants can be destabilized by forcing the separation. Finally, we note that the fusion mechanism is responsible for a cascade of coalescence events in a compact system of droplets where the separation is driven by surface tension.     

Multi-bubble motion behavior of uniform magnetic field based on phase field model

Aiming at the interaction and coalescence of bubbles in gas–liquid two-phase flow, a multi-field coupling model was established to simulate deformation and dynamics of multi-bubble in gas–liquid two-phase flow by coupling magnetic field, phase field, continuity equation, and momentum equation. Using the phase field method to capture the interface of two phases, the geometric deformation and dynamics of a pair of coaxial vertical rising bubbles under the applied uniform magnetic field in the vertical direction were investigated. The correctness of results is verified by mass conservation method and the comparison of the existing results. The results show that the applied uniform magnetic field can effectively shorten the distance between the leading bubble and the trailing bubble, the time of bubbles coalescence, and increase the velocity of bubbles coalescence. Within a certain range, as the intensity of the applied uniform magnetic field increases, the velocity of bubbles coalescence is proportional to the intensity of the magnetic field, and the time of bubbles coalescence is inversely proportional to the intensity of the magnetic field.     

Interfacial Area Transport Equation for Bubble Coalescence and Breakup: Developments and Comparisons

Bubble coalescence and breakup play important roles in physical-chemical processes and bubbles are treated in two groups in the interfacial area transport equation (IATE). This paper presents a review of IATE for bubble coalescence and breakup to model five bubble interaction mechanisms: bubble coalescence due to random collision, bubble coalescence due to wake entrainment, bubble breakup due to turbulent impact, bubble breakup due to shearing-off, and bubble breakup due to surface instability. In bubble coalescence, bubble size, velocity and collision frequency are dominant. In bubble breakup, the influence of viscous shear, shearing-off, and surface instability are neglected, and their corresponding theory and modelling are rare in the literature. Furthermore, combining turbulent kinetic energy and inertial force together is the best choice for the bubble breakup criterion. The reviewed one-group constitutive models include the one developed by Wu et al., Ishii and Kim, Hibiki and Ishii, Yao and Morel, and Nguyen et al. To extend the IATE prediction capability beyond bubbly flow, two-group IATE is needed and its performance is strongly dependent on the channel size and geometry. Therefore, constitutive models for two-group IATE in a three-type channel (i.e., narrow confined channel, round pipe and relatively larger pipe) are summarized. Although great progress in extending the IATE beyond churn-turbulent flow to churn-annual flow was made, there are still some issues in their modelling and experiments due to the highly distorted interface measurement. Regarded as the challenges to be addressed in the further study, some limitations of IATE general applicability and the directions for future development are highlighted.     

枝晶生长和气泡形成的数值模拟

本文建立了二维的格子玻尔兹曼方法-元胞自动机(lattice Boltzmann method-cellular automaton, LBM-CA)耦合模型, 对凝固过程中枝晶生长和气泡形成进行模拟研究. 本模型采用CA方法模拟枝晶的生长, 根据界面溶质平衡法计算枝晶生长的驱动力. 采用基于Shan-Chen多相流的LBM模拟气泡在液相中的生长和运动. 在LBM-CA的耦合模型中包含了固-液-气三相之间的相互作用. 应用Laplace定理和模拟气-液-固三相之间的润湿现象对模型进行了验证. 应用所建立的LBM-CA耦合模型模拟研究了气-液相互作用系数对单气泡生长的影响. 发现单气泡的生长速度和平衡半径随气-液相互作用系数的增大而增大. 定向凝固过程中枝晶和气泡生长的模拟结果再现了枝晶的择优生长、 气泡的优先形核位置、气泡的长大、合并、在枝晶间受挤变形以及在液相通道中的运动等物理现象, 与实验结果符合良好. 此外, 初始气体含量越高, 凝固结束时气泡的体积分数也相对较高. 本模型的模拟结果可以揭示在凝固过程中气泡形核、 生长和运动演化以及与枝晶生长相互作用的物理机理.