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

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

非牛顿流体液滴袋状破碎的数值模拟研究

本文利用简化的水平集–流体体积耦合算法(S-CLSVOF)及自适应网格加密技术,对非牛顿流体液滴袋状破碎过程进行三维数值模拟。基于Herschel-Bulkley型非牛顿流体模型,计算分析了在不同韦伯数下液滴袋状破碎的形变程度、气液界面发展过程及曳力系数变化趋势。结果表明:在低奥内佐格数下的液滴袋状破碎过程中,非牛顿流体相较于牛顿流体在破碎后出现了较为明显的液丝。而随着韦伯数增大,非牛顿流体液滴在流向与径向上的形变程度加剧,曳力系数增大,气液界面变化显著。     

超声速流体Kelvin-Helmholtz不稳定性速度梯度效应研究

利用加权本质上无振荡（WENO）方法模拟超声速流体Kelvin-Helmholtz（KH）不稳定性,研究速度梯度对KH不稳定性线性增长率和后期非线性演化的影响.模拟发现超声速流体中的速度梯度对KH不稳定性具有较强的致稳作用,给出了包含速度梯度致稳的线性增长率经验公式.数值模拟和经验公式符合得很好.模拟给出了清晰的流场密度等值线,这说明WENO方法模拟超声速流体KH不稳定性具有较好的界面变形捕捉能力.模拟结果表明速度梯度影响KH涡的演化,在给定密度梯度的情况下速度梯度越大KH涡的横向尺度越小. 关键词： Kelvin-Helmholtz不稳定性 超声速流体 速度梯度     

简单剪切流条件下液滴破碎和碰撞模拟

利用耗散粒子动力学(DPD)对简单剪切流条件下液滴的破碎和碰撞分离过程进行了模拟.液滴由多颗简单的DPD粒子构成,液滴与周围流体的不互溶性由提高它们之间的保守力系数获得.通道的上下壁面沿着相反方向运动从而使流体产生简单剪切流动.计算结果表明,在简单剪切流条件下,当界面张力数超过临界界面张力数后,液滴不稳定;液滴将被拉长并破碎成数个较小的稳定的液滴.当两个液滴相互碰撞并分离的现象发生时,两个液滴的质心在速度梯度方向上的距离与整个碰撞过程相关,碰撞之后的δz值比碰撞前的δz大.     

不同界面组分分布对Richtmyer-Meshkov不稳定性的影响

基于二维非定常Euler方程,对平面激波与不同界面组分分布下氦气气柱作用过程所引起的Richtmyer-Meshkov不稳定性现象进行了数值模拟,探讨了激波冲击轻质气柱后气柱界面形态的演变及流场波系结构,定量分析了气柱特征尺度(气柱长度、高度和中轴宽度)和气柱体积压缩率随时间变化.此外,结合流场压强、速度、环量和气体混合率,多角度分析了激波驱动界面气体混合的流动机制,获得了不同界面组分分布对界面不稳定性的影响.结果表明,随着气柱界面从完全扩散界面向间断界面的过渡,界面两侧的声反射系数随之增大,使入射激波与气柱界面的作用由常规透射转变为非常规透射,反射激波逐渐加强,透射激波逐渐减弱,使得Richtmyer-Meshkov不稳定性随之增强;同时,界面两侧阿特伍德数的增大,加强了Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性的发展.此外,界面不稳定性的加强使得流场环量增大,导致气体混合率的增长速率随之升高.     

同轴流动聚焦中射流不稳定性的理论研究

壳核结构的微胶囊在医学药学材料食品农业等领域具有广泛的应用前景, 其制备方法一直是相关领域关注的焦点.同轴流动聚焦(co-flow focusing)是一种新型制备技术, 利用复合射流的破碎制备微胶囊具有包裹率高过程量化可控参数域广产率高等诸多优势.在实验中, 复合射流的破碎受到多个过程参数的影响, 并涉及了多层界面的耦合效应.利用简化的物理模型, 在时间和时空域中分析了三相水-油-水复合射流不稳定性的发展和演化.在黏性流体线性稳定性理论中, 同轴射流和驱动液体的基本速度型分别基于管流和误差函数构造, 并通过数值方法求解满足相应边界条件下的线化小扰动控制方程.结果表明:增加内外层界面的界面张力均有利于射流的破碎; 流体的黏性对同轴射流的稳定性均有着促进作用; 越大的黏性越小的内界面张力对应着越大的射流破碎波长; 内外界面的耦合作用以及复合液滴的包裹情况均与内外射流的半径比息息相关; 绝对-对流不稳定性转换的临界Weber数随Reynolds数内层界面张力的增大而增大, 随内层和驱动流体的黏性增大而减小.这些结果将有助于提高液体驱动下同轴流动聚焦技术的过程控制, 为实际应用提供理论指导.      

表面张力对高雷诺数Rayleigh-Taylor不稳定性后期增长的影响

采用多相流的相场格子Boltzmann方法数值研究了微通道内高雷诺数单模Rayleigh-Taylor (RT)不稳定性的后期演化规律,重点分析表面张力对相界面动力学行为以及气泡与尖钉增长的影响.数值实验表明,随着界面张力的增大,可以有效降低演化过程中相界面结构的复杂程度,并抑制不稳定性后期相界面破裂形成离散液滴.另外,增大表面张力可以先促进后抑制气泡振幅的增长,而当表面张力较小时,尖钉振幅增长曲线之间并无明显差别,当表面张力增大到一定值后,它对尖钉振幅的抑制效果可明显地被观察到.进一步,根据不稳定性速度增长曲线,将高雷诺数单模RT不稳定性的演化划分为线性增长、饱和速度增长、重加速、混沌混合四个发展阶段.数值计算获取气泡与尖钉的饱和速度符合包含界面张力效应的势流理论模型.另外还统计了不同表面张力和Atwood数下表征RT不稳定性后期演化的气泡与尖钉增长率,结果显示气泡与尖钉后期增长率随着表面张力的增大总体上呈现出先促进后抑制的规律.最后,从数值计算和理论分析两方面研究了不同Atwood数下RT不稳定性发生的临界表面张力,发现两者结果符合得很好,并且临界表面张力随着流体Atwood数的增大而增大.     

液滴变形-袋式-多模式破碎转换研究

液滴破碎是影响内燃机燃烧的核心因素,其直接决定喷雾特性、蒸发速率以及混合气形成。为了揭示宽参数空间下液滴破碎模式之间的转换机理,本文采用简化耦合的多相流模型S-CLSVOF以及自适应网格加密技术对液滴在气流中的破碎过程进行了数值模拟研究。同时进一步将模拟结果与基于Rayleigh-Taylor失稳机理的理论模型进行比较并对其参数进行了修正。结合理论分析和数值模拟探讨了液滴破碎模式之间的转换。     

微观结构尺度对液滴润湿状态的影响

在分析表面微观特征对液滴润湿状态的影响时,尺度效应是无法绕过的前提条件和分析基础。本文以矩形微结构表面液滴为研究对象,通过数值模拟,得到了给定微结构表面液滴在不同液滴半径与矩形微结构宽度比例下的润湿状态。此外,本文还模拟得到了液滴完全浸润矩形微结构的临界接触角及其相对于微结构尺寸的变化规律。结果表明,液滴半径与微结构宽度比越小,液滴下气液界面开始浸润微结构时的临界接触角越大,无重力条件下表面液滴下气液界面的稳定性越差。若重力不可忽略,则液滴半径与微结构宽度比越大,重力的影响越大,表面液滴下气液界面的稳定性越差。     

等直径微液滴碰撞过程的改进光滑粒子动力学模拟

运用一种改进光滑粒子动力学(SPH)方法模拟了相溶和不相溶两种情况下的等直径微液滴碰撞动力学过程. 为提高传统SPH方法的数值精度和稳定性, 采用一种不涉及核导数计算的核梯度改进形式; 为处理液滴界面张力采用修正的van der Waals表面张力模型. 通过模拟牛顿液滴碰撞聚并变形过程并与相关文献或试验结果进行对比, 验证了改进SPH 方法模拟微液滴碰撞过程的可靠性. 随后, 研究了基于van der Waals模型相溶聚合物微液滴碰撞聚并变形过程及不相溶微液滴碰撞后的反弹、分离过程, 讨论了碰撞过程中碰撞速度、碰撞角度、密度比等参数对碰撞变形过程的影响, 分析了流体桥、旋转角度等因素的变化情况. 关键词： 光滑粒子动力学 微液滴 聚合物液滴 碰撞     

Modelling of mass-transfer induced instabilities at liquid–liquid interfaces based on molecular simulations

Interfaces and especially mass transfer across interfaces are of great importance in many fields of chemical engineering. Interfacial convection, which is generally called the Marangoni effect, may improve mass transfer across interfaces quite drastically and has not been investigated adequately in detail. In order to investigate the influence of mass transfer on a liquid–liquid interface molecular computer simulations have been performed. Since many molecules have to be considered for a significant modelling of the interface, cubic lattice systems have been chosen for the simulation which proceeds according to the Monte-Carlo scheme. The parameters that describe the thermodynamic and transport properties resemble those of realistic standard EFCE test systems for extraction. Results of various Monte-Carlo simulations show that under certain conditions mass transfer across interfaces induces the formation of nano droplets in the close vicinity of the interface. The different combinations of the nano droplet behaviour due to attractive or repulsive long-range forces together with the characteristics of coalescence may lead to different macroscopic interfacial instabilities such as spontaneous emulsification or eruptions. Based on diffusive and thermodynamic properties of the chosen lattice system a first stability criterion which allows the prediction of the onset of nano droplet formation is developed. The theoretical results compare well with experimental observations at a single drop and in a two-phase cell where the instabilities are investigated optically via Schlieren optics.     

Critical behavior of phase interfaces in porous media: Analysis of scaling properties with the use of noncoherent and coherent light

Image sequences of the surface of disordered layers of porous medium (paper) obtained under noncoherent and coherent illumination during capillary rise of a liquid are analyzed. As a result, principles that govern the critical behavior of the interface between liquid and gaseous phases during its pinning are established. By a cumulant analysis of speckle-modulated images of the surface and by the statistical analysis of binarized difference images of the surface under noncoherent illumination, it is shown that the macroscopic dynamics of the interface at the stage of pinning is mainly controlled by the power law dependence of the appearance rate of local instabilities (avalanches) of the interface on the critical parameter, whereas the growth dynamics of the local instabilities is controlled by the diffusion of a liquid in a layer and weakly depends on the critical parameter. A phenomenological model is proposed for the macroscopic dynamics of the phase interface for interpreting experimental data. The values of critical indices are determined that characterize the samples under test within this model. These values are compared with the results of numerical simulation for discrete models of directed percolation corresponding to the Kardar-Parisi-Zhang equation.     

倾斜壁面中Kelvin-Helmholtz不稳定性演化数值分析

基于FTM（Front Tracking Method）对倾斜壁面下的二维不混溶、不可压缩流体的Kelvin-Helmholtz（K-H）不稳定性进行数值模拟.研究壁面倾角,速度梯度层厚度以及理查德森数对K-H不稳定性发展的影响.研究表明,壁面倾角越大,K-H不稳定性发展越快,卷起的液体越多;倾斜壁面下速度梯度层厚度的增加对界面卷起表现出抑制作用.理查德森数重力项越大,界面卷起越缓慢,而理查德森数表面张力项对界面卷起高度的影响较小.     

液态金属自由表面膜流MHD效应的数值模拟

对液态金属自由表面膜流在强磁场下的磁流体力学效应进行了数值模拟研究,获得了液态金属自由表面的形状、截面流速分布及截面上的电动势分布,从而能对膜流的一些磁流体动力学行为作出解释。数值计算结果与理论分析和实验结果符合较好。由实验和数值模拟结果可以得出,液态金属膜流通过强磁场时,磁场会阻碍膜流的运动。     

基于不可压缩光滑粒子动力学的黏性液滴变形过程仿真

应用基于投影算法的不可压缩光滑粒子动力学(incompressible smoothed particle hydrodynamics, ISPH)法对黏性液滴变形过程进行了数值仿真. 对于张力失稳导致的粒子非物理簇集问题, 采用粒子移位技术加以解决. 为了验证本文ISPH 算法的精度和稳定性, 分别模拟了圆形黏性液滴的拉伸变形过程以及方形液滴的旋转变形过程, 得到了不同时刻液滴内部的压力变化特征, 准确地捕捉了液滴自由面演化过程, 并将数值计算结果与文献中的解析解进行了比较.分析结果表明, 基于投影算法的不可压缩光滑粒子动力学方法结合粒子移位技术, 能够有效地模拟黏性液滴变形过程, 可以得到精确和稳定的结果. 关键词： 不可压缩光滑粒子动力学 黏性液滴 自由面流动 数值仿真     

Fragmentation instabilities of a drop as it falls in a miscible fluid

We review here a series of experiments on the fragmentation instabilities that a liquid drop undergoes as it falls inside a fluid with which it is miscible, so called the solvent. Motivated by the original experiments initiated by Thomson and Newall in 1885, we started to investigate this subject more than one decade ago, encountering up to date a number of challenging problems in hydrodynamical instabilities, complicated by the presence of transient interfaces between the drop and the solvent. In particular, we have shown that when a drop of liquid is deposited over the surface of the same liquid, it falls down inside the solvent because the energy associated to its surface tension against air is instantaneously converted into kinetic energy. As a consequence, a very fast fluid injection takes place as the drop touches the free surface of the solvent and the drop enters inside the solvent. Due to the hydrodynamical instabilities related to the large velocity gradients, it develops into a ring that expands radially. The ring continues to go downwards inside the solvent until it stops at a certain height due to viscous dissipation. In the first stages of the ring expansion, a fluid membrane remains attached to the ring, so-called “turban" for its shape, which is concave or convex depending on the sign of the density difference between the drop and the solvent. When a small density difference is introduced between the drop and the solvent, the ring becomes unstable because of density gradients and it fragments into smaller droplets. On their turn, the secondary droplets may undergo the same instability and may fragment again, so that a cascade of fragmentation takes place. If the density difference is positive, that is, the drop is heavier than the solvent, the secondary droplets continue to go down deeper inside the solvent, until the whole process is washed out by the slow diffusion of the concentration gradients. When the density difference between the drop and the solvent is negative, then the secondary droplets rise up to the free surface of the solvent, where they are distorted by the equivalent of an impact with a rigid wall. Universal scaling laws for the cascade of fragmentation and for the dynamical behavior of the drop have been derived and accompany the experimental observations.     

On a moving liquid film and its instability on textured surfaces

Recently, two regimes of viscous friction on textured surfaces were proposed in the context of penetration of liquid film into the texture (EPL 79, 56005 (2007)): the Poiseuille and Stokes regimes. With this idea on viscous friction, we theoretically discuss instabilities on a liquid film on textured surfaces when the film is forced to move with external forces. When a film recedes due to a pressure drop, we find scaling laws for instabilities to be checked in future experiments. When a circular film expands due to centrifugal force we find that the expanding film is stable against rim fluctuations (within the linear stability analysis) with its radius determined by a simple equation. Our discussion sheds light on the curvature of the front of the moving liquid film on textured surfaces and how the film thickness is kept fixed to the texture height on textured surfaces, aspects which have not been discussed in previous studies.     

蒸发两层流系统的对流不稳定性分析

关于蒸发液层的Rayleigh-Marangoni-Benard不稳定性的研究中,早期文献中普遍采用的是单层流模型。近年来,一些学者采用两层流模型对蒸发稳定性进行了理论分析,有的文献中没有考虑蒸发率与饱和蒸汽压的耦合关系,所以得到的结果不能完全反应蒸发对系统稳定性的影响。本文建立了一种新的两层流模型,考虑了界面变形对系统稳定性的影响。采用线性稳定性方法对带有蒸发界面的两层流的Rayleigh-Marangoni-Benard对流不稳定性进行了分析,得到了临界 Marangoni数与波数的关系,重点讨论了蒸发系数以及重力对汽液两层流系统的不稳定性的影响。