制备复合液滴的微尺度流动方法

微尺度流动能够一步到位地制备不同结构和功能、尺寸在微米量级的复合液滴.文章回顾了几种常见的基于复合液滴的微尺度流动方法,包括同轴电雾化、复合流动聚焦、微流控芯片、玻璃微毛细管等,并对各种技术的原理和进展进行了简要概括和分析.在这类流动中,不同种类的流体在一定的几何结构通道或外力场作用下平稳地拉伸成微细射流并最终破碎成复合液滴.在同轴电雾化和复合流动聚焦技术中,从毛细管流出的流体能够形成稳定的锥-射流结构,当外力作用改变时能够形成不同的流动模式.在微流控芯片和玻璃微毛细管技术中,流体被约束在固定管道内,不同管道构型下能够形成不同的流动形态.这些方法都采用纯物理机理,过程稳定、易于操作,制备的复合液滴粒径可控,单分散性好,微观结构可设计,在科学研究和工程实际中具有重要的应用价值.      

液体横向射流在气膜作用下的破碎过程

为了研究液体横向射流在气膜作用下的破碎过程,采用背景光成像技术及VOF TO DPM方法进行了实验研究和仿真研究,模拟介质为水和空气.研究结果表明,液体射流在气膜作用下主要存在两种破碎过程:柱状破碎和表面破碎.Rayleigh-Taylor(R-T)不稳定性产生的表面波是液体射流发生柱状破碎的主要原因,气流穿透表面波的波谷导致射流柱破碎,破碎后的液丝沿流向逐渐发展呈带状分布.Kelvin-Helmholtz(K-H)不稳定性产生的表面波是液体射流发生表面破碎的主要原因,液丝和液滴从射流表面剥离.局部动量比对液体横向射流的破碎过程具有重要影响,当局部动量比较低时,液体射流的破碎由K-H不稳定性主导;随着局部动量比的增大液体射流的破碎逐渐由R-T不稳定性主导.液体射流的破碎长度及穿透深度均随局部动量比的增大而增大.     

基于格子Boltzmann方法的液液不混溶两相流动数值模拟

具有介观特性的格子Boltzmann方法能够准确方便地捕捉相界面细节,在两相流动领域具有广泛的应用前景.本文在简要介绍格子Boltzmann方法的基础上,利用格子Boltzmann方法的颜色模型对四类经典两相流动问题进行了模拟,界面张力的Laplace定律验证、单液滴松弛过程、两个液滴融合过程、水平通道内不混溶液液两相流动.结果表明,液滴界面张力符合Laplace定律;黏性越大,液滴松弛过程越稳定;界面张力越大,液滴融合速度越快;格子Boltzmann方法能够有效描述液液两相流动的界而信息。研究工作为应用格子Boltzmann方法分析两相流动问题奠定了理论基础.     

壁面吸附复合液滴的变形与运动特性

流动作用下壁面吸附复合液滴可呈现不同的动力学状态,了解其动力学状态转变及参数影响规律是操控液滴运移的关键.本文基于Front-tracking Method结合广义Navier边界条件建立了能够模拟移动接触线问题的数值方法。主要考察了层液滴毛细数、内层液滴大小、内层液滴毛细数对接触角滞后壁面吸附复合液滴动力学特性的影响。研究发现:外层液滴毛细数越大、内层液滴毛细数越小且内层液滴体积越大,复合液滴更容易从壁面脱附或被流体剪切发生破碎.     

气体中心式离心喷嘴喷雾实验与三维仿真

在实验的基础上, 基于RNG k-ε模型对常压下气体中心式同轴离心(gas-centered swirl coaxial,GCSC)喷嘴喷雾形态和破碎模式进行了三维仿真研究。采用网格自适应加密(adaptive mesh refinement,AMR)技术、耦合水平集和流体体积(coupled level-set and volume of fluid, CLSVOF)方法对气液界面进行捕捉。结果表明, 液体质量流率($\dot{m}_{\mathrm{l}}$)不变, 随着气体质量流率($\dot{m}_{\mathrm{g}}$)的增加, 中心气流的引射作用增强, 液膜内外压差增大, 雾化锥角减小, 并对其流动特性进行了分析; 而$\dot{m}_{\mathrm{g}}$不变时, 液膜在喷嘴出口的径向速度与切向速度随$\dot{m}_{\mathrm{l}}$的增大而增大, 导致雾化锥角增大。同时根据气液质量流率比(gas-liquid mass flow rate,GLR), 将喷雾的破碎模式分为穿孔破碎、气泡破碎和气动破碎。      

流体黏性及表面张力对气泡运动特性的影响

基于气泡边界层理论,引入黏性修正,采用边界积分法,考虑黏性效应和表面张力在单气泡以及双气泡耦合作用过程中的影响.首先将建立的数值模型与Rayleigh-Plesset的解析解进行对比,发现二者符合良好,验证了数值模型的有效性;在此基础上,建立考虑流体弱黏性效应的双气泡耦合模型,研究流体黏性和表面张力作用下,气泡表面变形、射流速度、流场能量转换等物理量的变化规律;最后研究雷诺数和韦伯数对于气泡脉动特性的影响规律.结果表明,流体黏性会抑制气泡脉动和气泡射流发展,降低气泡半径和射流速度;表面张力不改变气泡脉动幅值,但缩短了脉动周期,提升气泡势能.     

三相流体的轴对称格子Boltzmann模型及其在Rayleigh-Plateau不稳定性的应用

基于多组分相场理论提出了一类模拟三相流体流动的轴对称格子Boltzmann模型.该模型利用两个粒子分布函数来捕捉三种不同流体之间的相界面,另一个粒子分布函数来求解流体动力学方程以获得流场信息.为了刻画坐标变换引起的轴对称效应,巧妙地设计了演化方程中平衡态分布函数和外力项分布函数,从理论上保证本文模型可以正确恢复三相流体系统的宏观控制方程,并且轴对称效应产生的源项中不包含任何复杂的梯度项,从而比现有的轴对称格子Boltzmann模型更加简单高效.首先通过模拟一系列轴对称多相流的基准算例,包括静态的双液滴、液体透镜的扩展和二元流体Rayleigh-Plateau不稳定性,来验证本文模型的有效性与正确性.接下来,利用该模型研究了三相流体的Rayleigh-Plateau不稳定性的增长过程,定量分析了波数和液柱半径比对复合液体线程破裂过程中界面动力学行为、界面破裂时间以及生成子液滴尺寸的影响.可以发现复合的液体线程在波数较大时破裂生成一个复合主液滴和卫星液滴,而在波数较小时可以生成更多数量的卫星液滴,这导致复合主液滴和卫星液滴的尺寸随着波数的增加呈现先增大而后减少的趋势.另外,我们发现内部流体...     

粗糙纳通道内流体流动与传热的分子动力学模拟研究

为研究粗糙表面对纳尺度流体流动和传热及其流固界面速度滑移与温度阶跃的影响,本文建立了粗糙纳通道内流体流动和传热耦合过程的分子动力学模型,模拟研究了粗糙通道内流体的微观结构、速度和温度分布、速度滑移和温度阶跃并与光滑通道进行了比较,并分析了固液相互作用强度和壁面刚度对界面处速度滑移和温度阶跃的影响规律. 研究结果表明,在外力作用下,纳通道主流区域的速度分布呈抛物线分布,由于流体流动导致的黏性耗散使得纳通道内的温度分布呈四次方分布. 并且,在固体壁面处存在速度滑移与温度阶跃. 表面粗糙度的存在使得流体剪切流动产生了额外的黏性耗散,使得粗糙纳通道内的流体速度水平小于光滑通道,温度水平高于光滑通道,并且粗糙表面的速度滑移与温度阶跃均小于光滑通道. 另外,固液相互作用强度的增大和壁面刚度的减小均可导致界面处速度滑移和温度阶跃程度降低. 关键词： 速度滑移 温度阶跃 流固界面 粗糙度     

用改进的耦合型Level Set方法计算一维双介质可压缩流动

用带有虚拟流体(Ghost Fluid)修正的Level Set方法计算了一维可压缩双介质流动,把描述流动的Euler方程和描述流体界面运动的Level Set方程耦合起来,得到一个整体的守恒律系统,应用高分辨率差分格式求解;为了解决流体界面附近的数值跳动问题,在界面附近引入了虚拟流体方法的Isobaric修正,并给出了算例.     

微通道内光热相变驱动流体运动特性研究

本文通过可视化实验,对微通道内光热效应致相变驱动流体运动特性进行了研究,通过红外聚焦激光跟随微通道内液柱气液界面进行加热持续产生的蒸发冷凝-聚合过程对液柱进行连续驱动。实验研究了激光功率、加热点位置对相变过程中的界面行为、冷凝液滴分布、驱动速率的影响规律。结果表明,激光功率越高,光斑距离界面越近,液柱蒸发速率越大,蒸汽浓度高,冷凝液滴分布越密集,驱动流体流动速度越快。     

Simple and double emulsions via coaxial jet electrosprays

We report for the first time the generation of electrified coaxial jets of micrometric diameter in liquid media. Scaling laws to predict the inner and outer diameter of the coaxial jet are given. We show some experiments illustrating the formation process of the coaxial jet, and demonstrate how this process can be used to yield either o/w (oil in water) or o/w/o (oil/water/oil) emulsions of micrometric size. Some interesting analogies with other hydrodynamic focusing processes are also pointed out.     

Absolute instability of a liquid jet in a coflowing stream

Cylindrical liquid jets are inherently unstable and eventually break into drops due to the Rayleigh-Plateau instability, characterized by the growth of disturbances that are either convective or absolute in nature. Convective instabilities grow in amplitude as they are swept along by the flow, while absolute instabilities are disturbances that grow at a fixed spatial location. Liquid jets are nearly always convectively unstable. Here we show that two-phase jets can breakup due to an absolute instability that depends on the capillary number of the outer liquid, provided the Weber number of the inner liquid is >O(1). We verify our experimental observations with a linear stability analysis.     

Dripping to jetting transitions in coflowing liquid streams

A liquid forced through an orifice into an immiscible fluid ultimately breaks into drops due to surface tension. Drop formation can occur right at the orifice in a dripping process. Alternatively, the inner fluid can form a jet, which breaks into drops further downstream. The transition from dripping to jetting is not understood for coflowing fluid streams, unlike the case of drop formation in air. We show that in a coflowing stream this transition can be characterized by a state diagram that depends on the capillary number of the outer fluid and the Weber number of the inner fluid.     

Interaction of inner and outer layers in plane and radial wall jets

Large eddy simulations of turbulent radial and plane wall jets were performed at different Reynolds numbers using the Lagrangian dynamic eddy viscosity subgrid-scale model. The results were validated with experimental data available in the literature. Compared to the plane ones, the radial wall jets have an extra direction for expansion, which causes faster decay rates. Thus, the resulting pressure gradient distributions are different. However, the comparison of the results with the turbulent boundary layers under adverse and favourable pressure gradients reveals that these pressure gradients are not strong enough to cause any fundamental physical difference between plane and radial wall jets. In both cases, the local Reynolds number is an important determining factor in characterisation of the flow. The joint probability density function analysis shows that the local Reynolds number determines the level of intrusion of the outer layer into the inner layer: the lower the local Reynolds number, the stronger is the interaction of the inner and outer layers. These results can be used to clarify the scatter of the reported log-law constants in the literature.     

Viscoelastic theory for nematic interfaces

A complete macroscopic theory for compressible nematic-viscous fluid interfaces is developed and used to characterize the interfacial elastic, viscous, and viscoelastic material properties. The derived expression for the interfacial stress tensor includes elastic and viscous components. Surface gradients of the interfacial elastic stress tensor generates tangential Marangoni forces as well as normal forces. The latter may be present even in planar surfaces, implying that in principle static planar interfaces may accommodate pressure jumps. The asymmetric interfacial viscous stress tensor takes into account the surface nematic ordering and is given in terms of the interfacial rate of deformation and interfacial Jaumann derivative. The material function that describes the anisotropic viscoelasticity is the dynamic interfacial tension, which includes the interfacial tension and dilational viscosities. Viscous dissipation due to interfacial compressibility is described by the anisotropic dilational viscosity, and it is shown to describe the Boussinesq surface fluid appropriate for Newtonian interfaces when the director is homeotropic. Three characteristic interfacial shear viscosities are defined according to whether the surface orientation is along the velocity direction, the velocity gradient, or the unit normal. In the last case the expression reduces to the interfacial shear viscosity of the Boussinesq surface fluid. The theory provides a theoretical framework to study interfacial stability, thin liquid film stability and hydrodynamics, and any other interfacial rheology phenomena.     

Effects of non-uniform interfacial tension in small Reynolds number flow past a spherical liquid drop

A singular perturbation solution is given for small Reynolds number flow past a spherical liquid drop. The interfacial tension required to maintain the drop in a spherical shape is calculated. When the interfacial tension gradient exceeds a critical value, a region of reversed flow occurs on the interface at the rear and the interior flow splits into two parts with reversed circulation at the rear. The magnitude of the interior fluid velocity is small, of order the Reynolds number. A thin transition layer attached to the drop at the rear occurs in the exterior flow. The effects could model the stagnant cap which forms as surfactant is added but the results apply however the variability in the interfacial tension might have been induced.     

Comparative visualized investigation of impact-driven high-speed liquid jets injected in submerged water and in ambient air

This paper is a comparative study on the characteristics of high-speed liquid jets injected in surrounding water and air using shadowgraph technique. One of the main objectives is to investigate the effects of liquid’s physical properties, used to generate the high-speed liquid jets, on jet generation’s characteristics. Moreover, comparative investigations on effects of those liquid jets after injected in water and air are reported. The high-speed liquid jets were generated by the impact of a projectile launched by a horizontal single-stage power gun. The impact-driven high-speed liquid jets were visualized by shadowgraph technique and images were recorded by a high-speed digital video camera. The process of impact-driven high-speed liquid jet injection in air and water, oblique shock waves, jet-induced shock waves, shock waves propagation, the bubble behavior, bubble collapse-induced rebound shock waves and bubble cloud regeneration were clearly observed. It was found that different properties of liquid (surface tension and kinematic viscosity) affect the jet maximum velocity and shape of the jet. Bubble behaviors were only found for the jet injected in water. From the shadowgraph images, it is found that the maximum average jet velocity, expansion and contraction velocities of bubble in axial direction increase when the value of the multiplied result of surface tension by kinematic viscosity increases. Therefore, surface tension and kinematic viscosity are the significant physical properties that affect characteristics of high-speed liquid jets.     

Numerical visualization of vortex flow behavior in square jets in cross-flow

Direct numerical simulations have been performed in this study to visualize the flow behavior of single and multiple square jets issuing normally into a cross-flow. Three configurations are considered, a single jet located in the centre of the domain, twin jets in side-by-side (SBS) arrangement in the spanwise direction and triple jets in tandem arrangement with twin jets at the front and a third jet in downstream along the centre line. Simulation uses a jet to cross-flow velocity ratio of 2.5 and the Reynolds number 225, based on the free-stream quantities and the jet width. While the vortical structures predicted from single jet case were in good qualitatively agreement with the findings of other researchers, our results show that the process of merging between two counter-rotating vortex pairs (CRVP) in twin jets configurations is strongly dependent on the jet-to-jet edge distance. Further downstream in the far-field, results from the SBS twin jets show a most dominating larger CRVP accompanied with a smaller inner vortex pair. The observations are in good qualitative agreement with the experimental findings in the literature. The resulting flow structures of triple jets in tandem configuration have revealed, for the first time, more complicated flow interactions between individual jets and cross-flow, providing further insights of complex flow physics and its potential engineering applications.     

Modified infiltration of solvated ions and ionic liquid in a nanoporous carbon

Infiltration of ions in a nanoporous carbon is responsive to the external electric field. If the liquid phase is an aqueous solution of electrolyte, the effective solid-liquid interfacial tension decreases as the voltage rises, similar to the electrowetting phenomenon at a large graphite surface. If the liquid phase is an ionic liquid, however, the effective interfacial tension increases with the voltage. The accessible nanopore volume is not dependent on the electric field. The unique phenomena should be related to the confinement effects of the nanopore inner surfaces.