液滴撞击超疏水壁面反弹及破碎行为研究

液滴撞击壁面时,壁面亲水性对液滴撞击壁面后的变化历程具有重要的影响。利用相界面追踪的复合Level Set-VOF方法对液滴撞击超疏水壁面的运动进行了研究。研究结果表明,撞击速度较小时,液滴撞壁后发生反弹;撞击速度较大时,液滴撞壁后会发生破碎现象;初始粒径的增大和表面张力的减小,有利于液滴撞壁后产生铺展破碎现象;撞击角度对撞壁后的液滴行为具有较大的影响。通过数值模拟,给出了一定条件下液滴垂直及倾斜撞击超疏水壁面反弹及破碎的临界条件。     

含气泡油滴撞击油膜壁面时气泡的变形与破裂

油--气润滑过程中润滑油液滴受高速气流扰动易形成含气泡油滴,微气泡将对油滴撞击壁面时的运动过程以及壁面油膜 层的形成质量产生重要影响. 基于耦合的水平集--体积分数 方法,对含气泡油滴撞击油膜壁面行为进行数值模拟研究, 考察含气泡油滴撞击油膜壁面时气泡的变形运动过程,探讨气泡破裂的动力学机制,分析气泡大小、碰撞速度和液体黏度等因素对含气 泡油滴撞壁过程中气泡变形特征参数的影响规律. 研究表明:含气泡油滴撞击油膜壁面后气泡会发生变形,并破裂形成膜液滴;气泡随同 液滴运动过程中,气泡内外压力和速度梯度变化是使气泡发生破裂的主要诱因. 气泡大小对气泡破裂方式影响较大,气泡较小时发生单 点破裂,而气泡较大时更容易发生多处破裂. 不同大小气泡受力差异较大,气泡大小与破裂发生时刻没有明显相关性. 碰撞速度和液体 黏度对气泡的变形、破裂和破裂发生时刻都具有一定的影响. 碰撞速度越大,油滴动能越大,更容易产生气泡变形和破裂现象. 液体黏 度增大,在油滴撞壁运动前期促进气泡变形,而在运动后期可以阻延气泡破裂行为发生.      

Fluctuating flow in a liquid layer and secondary spray created by an impacting spray

A spray impacting onto a wall produces a flow of secondary droplets. For relatively sparse spray these secondary droplets are produced by the splashing of the impacting drops and their interactions. For dense sprays, like Diesel injection sprays, these secondary droplets are created by the fluctuating liquid film created on the wall. In the present paper hydrodynamic models are presented for these two extreme cases. The velocities of the secondary droplets produced by the crown splash in a sparse spray are described theoretically. Next, the fluctuations in the motion of the liquid film created by a dense impacting spray are analyzed statistically. This motion yields the formation of finger-like jets, as observed in experiments of a Diesel spray impacting onto a rigid wall. The characteristic size and velocity of the film fluctuations are estimated. These two theoretical models are validated by comparison with the experimental data.     

强磁场影响下镓液滴撞击固壁的实验研究

在材料的电磁冶金过程及磁约束核聚变装置中, 金属液滴在磁场和壁面温度影响下的撞击过程表现出复杂的动力学特性. 本文对水平磁场作用下液态镓(Ga)液滴撞击等温和过冷壁面的铺展和回弹特性进行了实验研究. 采用高速相机拍摄液滴撞击过程中轮廓的变化, 通过图像处理获得不同磁场强度、不同撞击速度和不同底板温度下的最大铺展因子、回弹过程中的最大高度以及产生的二次液滴的半径和速度. 碰撞速度由0.45 ~ 1.8 m/s, 磁场强度从0 ~ 1.6 T, 底板温度为30 °C, ?20 °C和?10 °C. 基于实验结果分析了磁场和壁面温度对液滴铺展和回弹的影响规律. 实验结果表明, 液滴撞击等温壁面和过冷壁面的最大铺展因子随We的变化均与理论预测关系式一致. 液滴撞击等温壁面的情况下, 不同的We下, 出现不同的回弹现象. 磁场抑制了平行于磁场方向的液滴铺展和回弹过程中二次液滴的产生, 而对回弹过程中的液滴在平行磁场方向上有拉伸作用. 液滴撞击过冷壁面时, 在一定的We值范围内, 同样会出现二次液滴分离现象, 此时产生的二次液滴的速度较小. 磁场的增强和We的增大都会导致液滴在高度方向的振荡减弱, 加速凝固过程.      

A two-dimensional splashing model for investigating impingement characteristics of supercooled large droplets

In this article, a two-dimensional (2D) splashing model is proposed to investigate the dynamics when supercooled large droplets (SLD) impinging on a wall surface in the aircraft-icing field. Energy conservation for droplet motion and impingement is used to capture the properties of the splashed droplets. A new statistical treatment of the droplet impinging energy and angle during the droplet-wall interaction is introduced in order to calculate the average dynamics of the SLD within a micro-control volume on wall surface. Based on the LEWICE predictions of droplet collection efficiencies and the available experimental ones, a new criterion for droplet splashing/deposition as well as a new formulation for the splashed mass is suggested. Lagrangian approach is adopted to describe the movement and impingement of SLD. The proposed model together with the previously developed droplet tracking method (DTM) for calculating droplet collection efficiency with the effect of droplet reimpingement constitute a relatively complete predicting approach of SLD impingement characteristics. Comparisons between the current predictions and the experimental observations, including SLD impingement over clean and contaminated airfoil surfaces as well as shapes of ice accretion in typical icing conditions, are carried out. Further, results obtained with the LEWICE splashing model are also plotted on the same graphs in order to assess the accuracy of the current splashing model in predicting SLD impingement. Results show that good agreement is achieved between the current predictions, including SLD impingement and ice accretion shapes, and the experimental ones. The predictions of the impingement distribution over contaminated surfaces obtained with the current splashing model show a much closer agreement with the experimental results than the ones obtained with LEWICE splashing model. For further investigation of SLD impingement, the properties of the droplet splashing and reimpingement during the ice accretion process are also addressed.     

Three-Dimensional droplet splashing dynamics measurement with a stereoscopic shadowgraph system

In the current study, a novel three-dimensional (3D) measurement technique is established using a high-speed stereoscopic shadowgraph system, which is applied to investigate the 3D splashing dynamics of silicone oil dropping on liquid films quantitatively. Both crown-type splashing and crown wall bottom breakdown splashing morphologies are involved for measurement and comparison. Based on the shadowgraph images, the secondary droplets diameter, number and mass fraction are determined. The three-dimensional (3D) coordinates of the secondary droplets during the splashing processes are reconstructed and tracked, based on which the 3D trajectories, velocities, ejecting angles and kinetic energy are calculated and analysed. It is found that the secondary droplets in crown bottom breakdown cases have a larger mass fraction and kinetic energy than that in crown-type splashing cases. The measurement indicates that the radial velocity increases with the decreasing of film thickness, while the vertical velocity does not vary too much. Significant disparities between two-dimensional (2D) and three-dimensional (3D) velocity magnitudes as well as kinetic energy have been identified, which indicate that accurate time-resolved 3D measurements are of great importance for quantitative investigation of splashing phenomena and the high-speed stereoscopic shadowgraph system has been proved to be able to play a role.     

The combined influence of a rough surface and thin fluid film upon the splashing threshold and splash dynamics of a droplet impacting onto them

Surface roughness can have a critical effect upon the splashing threshold and dynamics of a drop impacting on either a dry or rough solid surface or one coated by a thin fluid film. As most coating applications and spray systems quickly evolve to a state where the droplets impinge upon fluid deposited by preceding droplets, the combined contributions of surface roughness and a pre-deposited thin liquid film of comparable thickness upon droplet impingement dynamics are examined. For comparison, we include results for droplets impacting on a smooth, dry surface and a smooth surface wetted by a thin fluid film. The inclusion of surface roughness considerably lowers the splashing threshold and alters the splashing dynamics such that differences in fluid surface tensions between 20.1 and 72.8 dynes/cm or viscosities between 0.4 and 3.3 cP have little effect.     

椭球液滴撞击超疏水表面反弹过程数值分析

液滴在自由落体或受外力作用时常发生椭球形变,对其撞击超疏水壁面的运动形态及形成二次液滴有较大影响.本文假定具有不同轴向半径比值(AR)的椭球形液滴,采用耦合水平集-体积分数(CLSVOF)方法对椭球形液滴撞击超疏水平壁面进行数值模拟研究,对椭球形液滴撞击超疏水平壁面反弹过程运动形态和AR对二次液滴形成的影响进行分析.研究表明,不同AR液滴撞击超疏水壁面后反弹过程具有一定相似性,同时存在明显差异.液滴反弹过程中会拉伸形成长液柱,其未扰动界面直径随AR的增大而减小.在Plateau-Rayleigh不稳定性影响下,长液柱会发生破碎形成二次液滴,但较低的AR对二次液滴的形成有较明显的抑制作用,当AR小于临界值0.6时,反弹过程中液滴内部压力均匀稳定,最终不产生二次液滴.     

A STUDY ON DROPLET VELOCITY OF ETHANOL DURING ELECTROSPRAYING PROCESS AT CONE-JET MODE

研究液滴在静电喷雾下的速度特性是理解喷雾形态的形成及演化的关键.结合锥射流模式下乙醇静电喷雾实验数据,建立了静电喷雾二维轴对称模型.基于离散相液滴运动方程、连续相空气运动方程、电场方程以及用户自定义函数,进行了数值求解,获得了锥射流模式下的乙醇静电喷雾形态、空间电场分布以及液滴速度场分布.考虑了不同空气入口流速的影响,得到了乙醇/空气同轴射流静电喷雾形态的变化规律.结果表明,喷雾外围液滴与空气流场有较强的相互作用,导致喷雾中轴线附近的液滴速度分布变化较小,而在喷雾外围处的液滴速度分布沿径向剧烈变化;随着空气入口速度的增大,乙醇/空气同轴射流静电喷雾形态先趋于发散,当空气入口速度大于喷雾外围液滴轴向速度时,喷雾形态则趋于聚拢.因此,除改变施加电压、液体流量和电极结构外,通过控制空气入口速度来影响喷雾液滴速度场,也可获得不同的静电喷雾效果.     

INVESTIGATION ON THE ATTENUATION OF A ROTATING LIQUID MOTION IN A CYLINDER WITH BOUNDARY RESISTANCE

为了分析圆柱形容器内同步旋转液面在边界阻力作用下的衰减机理，采用层流摩阻与湍流摩阻的理论模型和实验测试的方法对同步旋转液面在边壁阻力作用下的衰减规律进行了研究. 通过分析壁面阻力的特性，采用当量湿周对圆柱形容器中旋转液体运动的能量损失进行模拟，建立了壁面糙率、旋转半径，静水深度和液体旋转速度之间的本构方程. 通过对圆柱容器中的同步旋转液面在边界阻力作用下的层流和湍流形态的衰减过程进行实验测试，验证了理论模型的正确性. 基于验证后的理论方法，获得了圆柱容器中旋转液体的角速度、高度等随时间的变化过程，分析了摩阻形态、糙率和直径对旋转液面衰减过程的影响. 结果表明，边壁糙率和容器半径是决定旋转液面的衰减过程主要因素，边壁糙率越大液体旋转速度衰减越快，圆柱形容器半径越大液体旋转的衰减越慢.     

CHARACTERISTIC OF DROPLET OSCILLATION ON THE SURFACE OF RECTANGULAR HYDROPHOBIC GROOVES

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

Atomized non-equilibrium two-phase flow in mesochannels: Momentum analysis

Two-phase pressure drop measurements are very difficult to make while the fluid is in non-equilibrium condition, i.e. while phase change is taking place. This is further complicated when an atomized liquid is introduced in the system at much higher velocity than other components such as liquid layer, vapor core, and entrained droplets. The purpose of this paper is to develop a model to predict the two-phase pressure characteristics in a mesochannel under various heat flux and liquid atomization conditions. This model includes the momentum effects of liquid droplets from entrainment and atomization. To verify the model, an in-house experimental setup consisting of a series of converging mesochannels, an atomization facility and a heat source was developed. The two-phase pressure of boiling PF5050 was measured along the wall of a mesochannel. The one-dimensional model shows good agreement with the experimental data. The effects of channel wall angle, droplet velocity and spray mass fraction on two-phase pressure characteristics are predicted. Numerical results show that an optimal spray cooling unit can be designed by optimizing channel wall angle and droplet velocity.     

Phase Doppler measurements of spray impact onto rigid walls

In this work, an experimental study of spray impact onto a horizontal flat and rigid surface is presented. The phase Doppler technique has been used to characterize both the impacting and the secondary spray in terms of mass and number flux, size distribution and velocities of the droplets above the target. A high-resolution CCD camera has been used to measure the average liquid film thickness formed due to spray impact, whereas a high-speed CMOS camera has been used to characterize the splashing droplets from the wall. This visualization of the splashing phenomenon and the knowledge about the liquid film thickness are used to formulate a new physical model of the crown evolution. Furthermore, information about the incident-to-ejected mass fraction and number fraction are novel contributions of this study. Considerable data are provided comparing the impact of single drops onto a liquid film to impact of drops in a spray, and the significance of the observed differences for modelling efforts is discussed. The measurements of this study are also shown to be rather sensitive to the placement of the phase Doppler measurement volume above the surface and to the operating parameters of the instrument. These effects have been documented and discussed for this particular measurement situation.     

Numerical study on dynamic characteristics of double droplets impacting a super-hydrophobic tube with different impact velocities

ABSTRACT

The coupled level set and volume of fluid method is applied to the numerical study on the successive impact of double droplets on a super-hydrophobic tube. The impact velocity varies from 0.25 to 2?m/s. These impact processes present spread, retract, rebound, breakup and splash. The out-of-phase impact takes place with the impact velocity from 0.25 to 1.25?m/s, while the in-phase impact takes place with the impact velocity from 1.44 to 2?m/s. With the impact velocity larger than 1.25?m/s, the liquid crown presents and deforms after the trailing droplet impact, then it would gather at the film edge, rebound or break up. When impact velocities range from 1.44 to 1.5?m/s, the finger liquid film presents before the liquid crown appearing. The finger head breaks with the impact velocity of 1.5?m/s during the leading droplet spreading. The zigzag liquid film becomes more obvious for larger velocities.     

基于SPH方法的湿润性固壁边界条件模拟研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

DROPLET INTERACTIONS IN THERMOCAPILLARY MIGRATION

通过密度匹配方法地面模拟微重力环境,以豆油和硅油作为实验系统母液和液滴,研究非均匀温度场中液滴热毛细迁移行为. 用马赫- 泽德干涉仪测量实验系统背景温度场. 用阴影法观测液滴热毛细迁移. 单液滴实验中,得到与YGB 理论预测趋势相一致的结果;双滴实验中,观察到小液滴超越大液滴现象;实验表明液滴迁移过程中液体的黏性系数起重要作用,小液滴超越大液滴过程中,两液滴水平位移都出现明显偏移,且大液滴横向位移更显著. 创新性的设计了三相多液滴热毛细迁移实验,为实际工业生产中多相分离问题提供科学研究基础.     

Maximum spreading of electrically charged droplets impacting on dielectric substrates

In this work, the electric charging effect on the spreading of droplet impacting on dielectric substrates has been investigated. The charged water droplets were directed on the paraffin wax and the Teflon-coated plates. The impinging behavior was visualized and recorded using a CCD camera to identify the maximum extent of the flattened droplets. Droplet diameter and velocity approaching the wall were measured as well. The diameter of the electrically charged droplet at the maximum spread turned out to be larger compared to that of neutral droplet (at the maximum spread), and the difference becomes larger with increasing of the electric charge ratio (defined as the ratio of the actual electric charge to the Rayleigh limit). This phenomenon is considered to be due to reduction of effective interfacial tensions between the liquid and the gas and between the liquid and the solid by electric charging. Finally, an improved model was proposed to predict the maximum spreading ratio for electrically charged droplets by introducing correlations on the liquid–gas and the liquid–solid interfacial tensions.     

Impact dynamics of drops on thin films of viscoelastic wormlike micelle solutions

The impact dynamics of water drops on thin films of viscoelastic wormlike micelle solutions is experimentally studied using a high-speed digital video camera at frame rates up to 4000 frame/s. The composition and thickness of the thin film is modified to investigate the effect of fluid rheology on the evolution of crown growth, the formation of satellite droplets and the formation of the Worthington jet. The experiments are performed using a series of wormlike micelle solutions composed of a surfactant, cetyltrimethylammonium bromide (CTAB), and a salt, sodium salicylate (NaSal), in deionized water. The linear viscoelastic shear rheology of the wormlike micelle solutions is well described by a Maxwell model with a single relaxation time while the steady shear rheology is found to shear thin quite heavily. In transient homogeneous uniaxial extension, the wormlike micelle solutions demonstrate significant strain hardening. The size and velocity of the impacting drop is varied to study the relative importance of Weber, Ohnesorge, and Deborah numbers on the impact dynamics. The addition of elasticity to the thin film fluid is found to suppress the crown growth and the formation of satellite drops with the largest effects observed at small film thicknesses. A new form of the splashing threshold is postulated which accounts for the effects of viscoelasticity and collapses the satellite droplet data onto a single master curve dependent only on dimensionless film thickness and the underlying surface roughness. Additionally, a plateau is observed in the growth of the maximum height of the Worthington jet height with increasing impact velocity. It is postulated that the complex behavior of the Worthington jet growth is the result of a dissipative mechanism stemming from the scission of wormlike micelles.     

A numerical investigation of droplet impact on a heated wall in the film boiling regime

In this research, a novel approach is proposed to simulate the impingement of a water droplet on a superheated wall by solving conservation equations of mass, momentum and energy. Due to high temperature of the surface, a vapor layer is formed between droplet and surface. The vapor layer is captured using a very fine mesh near the surface. The level set method is applied for interface tracking and appropriate jump conditions are imposed at the interface by the ghost fluid method. The proposed algorithm is validated by comparing numerical results to the available experimental and analytical solutions. The effect of impact velocity, surfactant and polymer additives on the heat removal is studied. Simulations show that an increase in impact velocity enhances the heat removal. On the other hand, polymer additive has no significant effect while surfactant enhances dissipated heat by increasing the contact time or even by sticking the droplet on the surface.     

基于LS-DYNA的热成型钢断裂失效预测研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。