A comparative study of droplet impact dynamics on a dual-scaled superhydrophobic surface and lotus leaf

The impact dynamics of water droplets on an artificial dual-scaled superhydrophobic surface was studied and compared with that of a lotus leaf with impact velocity V up to 3 m/s. The lower critical impact velocity for the bouncing of droplets was about 0.08 m/s on both surfaces. At relatively low impact velocities, regular rebound of droplets and air bubble trapping and flow jetting on both surfaces were observed as V was increased. For intermediate V, partial pinning and rebound of droplets were found on the artificial dual-scaled surface due to the penetration of the droplets into the micro- and nano-scale roughness. On the lotus leaf, however, the droplets bounced off with intensive vibrations instead of being partially pinned on the surface because of the irregular distribution of microbumps on the leaf. As the impact velocity was sufficiently high, droplet splashing occurred on both surfaces. The contact time and restitution coefficient of the impinging droplets were also measured and discussed.     

液滴撞击加热壁面传热实验研究

本文采用高速摄像仪对水滴和乙醇液滴撞击加热壁面后的蒸发过程进行了实验观测, 分析了液滴撞击加热壁面后的蒸发特性参数. 实验中, 两种液体初始温度均为20 ℃, 不锈钢壁面初始温度范围为68-126℃. 水滴初始直径为2.07 mm, 撞击壁面时Weber 数为2-44; 乙醇液滴初始直径为1.64 mm, Weber数为3-88. 结果表明, 液滴受到重力、表面张力及流动性的影响, 在蒸发过程的大部分时间内, 水滴高度持续降低而接触直径几乎不变; 蒸发后期, 液滴发生回缩, 水滴的接触直径、高度和接触角出现振荡现象. 乙醇液滴的接触角随时间的增加呈现先减小随后保持不变的趋势, 而接触直径和高度则持续减小, 直到液滴完全蒸发. 液滴蒸发总时长与液体物性和壁面温度有关, 随壁面温度的升高而减小, 与液滴撞击壁面时的Weber 数无关. 同时, 随着壁面温度的升高, 液滴显热部分占总换热量的比重增大, 显热部分能量不可忽略, 本文实验条件下得到水滴的平均热流密度为0.014-0.110 W·mm^-2.     

液滴撞击固体表面铺展特性的实验研究

利用高速摄影仪记录了不同液滴撞击固体表面的形态变化过程, 并探讨了液滴撞击参数对撞击过程液滴形态的影响. 结果表明: 液滴黏度对液滴铺展过程起着决定性作用, 液滴表面张力对液滴铺展后的回缩起到主导作用, 两者的共同作用决定着液滴的震荡特性; 撞击速度的提高会增大液滴的最大铺展因子, 但达到最大铺展的时间因表面张力的不同呈现不同变化规律.     

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

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

实验观测液滴撞击倾斜表面液膜的特殊现象

采用高速摄像仪以10000帧/s 的拍摄速度对液滴撞击倾斜表面液膜的过程进行了实验观测, 分析了液滴撞击倾斜表面液膜后的铺展、水花形成以及飞溅等现象, 考察了撞击角对液滴震荡变形过程的影响; 在此基础上, 定量讨论了液滴铺展速度随时间的变化规律, 揭示了液滴撞击速度和撞击角对前、后铺展因子及初始铺展速度的影响.观测发现, 在撞击角为28.0°–74.7°范围内, 随着撞击角的减小, 液滴在液膜表面的震荡变形程度增大; 前铺展因子随撞击速度的增大而增大, 随撞击角的减小而增大; 后铺展因子随撞击速度的增大几乎不发生变化, 但是随撞击角的增大而增大; 液滴初始铺展速度随撞击速度和撞击角的升高而增大. 关键词： 液滴撞击 倾斜液膜 铺展因子 铺展速度     

Numerical and analytic investigation of the dynamics of hollow droplet impact onto substrate

A numerical and analytical model is proposed for the impact of a hollow melt droplet onto the surface of a solid polished substrate. The model is based on integral laws of the mass and energy conservation of the colliding droplet, it accounts for capillary and adhesion properties of the melt. The main parameters of the high-velocity deformation of a hollow particle have been computed: the variation of its height, shell thickness, and the contact spot diameter up to the moment of the spreading droplet solidification as well as the pressure variation inside the droplet until the moment of the shell rupture. The critical pressure value at which the rupture occurs is estimated by a formula characterizing the spherical shell strength. Quite a fair agreement of the computed values of the final diameter of the splat of a spread and solidified droplet with the data of physical experiment is shown. The work was financially supported by the Russian Foundation for Basic Research (Grant No. 06-01-00080) as well as within the framework of the Interdisciplinary Integration Project of SB RAS No. 90 for the years 2006–2008 “Scientific Fundamentals of the Creation of Multi-Layer Nano-Structural Coatings and Interfaces in Heterogeneous Materials for Operation in the Fields of Extreme External Effects”.     

Shallow surface etching of organic and inorganic compounds by electrospray droplet impact

The electrospray droplet impact (EDI) was applied to bradykinin, polyethylene terephthalate (PET), SiO₂/Si, and indium phosphide (InP). It was found that bradykinin deposited on the stainless steel substrate was ionized/desorbed without the accumulation of radiation products. The film thickness desorbed by a single collisional event was found to be less than 10 monolayers. In the EDI mass spectra for PET, several fragment ions were observed but the XPS spectra did not change with prolonged cluster irradiation. The etching rate for SiO₂ by EDI was measured to be ∼0.2 nm/min. The surface roughness of InP etched by EDI was found to be one order of magnitude smaller than that etched by 3 keV Ar⁺ for about the same etching depths. EDI is capable of shallow surface etching with little damage left on the etched surface.     

Experimental investigation of spray characteristics of a liquid jet in a turbulent subsonic gaseous crossflow

Two perforated plates with different solidity ratios, S=50% and 67%, were used to investigate the effect of the velocity fluctuations of a subsonic gaseous crossflow on the spray characteristics of a liquid jet including droplet size and velocity distributions. The experiments were conducted over a range of jet-to-crossflow momentum flux ratio of q=16.5-172, and two gas Weber numbers of ?We_g=2.7 and 5.9, corresponding to the enhanced capillary breakup and bag breakup regimes, respectively. The experimental results of this study revealed that the distribution of droplets size associated with a turbulent and a uniform crossflow for each specific breakup regime were approximately identical. The bimodal and single peak distributions of droplets size, respectively, associated with enhanced capillary and bag breakup regimes were generally consistent with the literature reports. However, the transition of the liquid primary breakup regime from enhanced capillary to bag breakup mode was delayed in a turbulent crossflow compared to its uniform counterpart. The general behavior of droplets size-velocity profiles were also consistent with the literature reports. Nonetheless, complex variations in the distribution of droplets velocity when changing the crossflow turbulence intensity were observed and linked with the presence of instabilities on the liquid jet's surface. Finally, the present experiments allowed shedding more light on the reason why the breakup mechanisms of a liquid jet in a conventional uniform crossflow should not be generalized to predict the distinct breakup process of a liquid jet in a turbulent crossflow.     

Experimental analysis of shape deformation of evaporating droplet using Legendre polynomials

Experiments involving heating of liquid droplets which are acoustically levitated, reveal specific modes of oscillations. For a given radiation flux, certain fluid droplets undergo distortion leading to catastrophic bag type breakup. The voltage of the acoustic levitator has been kept constant to operate at a nominal acoustic pressure intensity, throughout the experiments. Thus the droplet shape instabilities are primarily a consequence of droplet heating through vapor pressure, surface tension and viscosity. A novel approach is used by employing Legendre polynomials for the mode shape approximation to describe the thermally induced instabilities. The two dominant Legendre modes essentially reflect (a) the droplet size reduction due to evaporation, and (b) the deformation around the equilibrium shape. Dissipation and inter-coupling of modal energy lead to stable droplet shape while accumulation of the same ultimately results in droplet breakup.     

Simulation of impact of a hollow droplet on a flat surface

Despite many theoretical and experimental works dealing with the impact of dense continuous liquid droplets on a flat surface, the dynamics of the impact of hollow liquid droplets is not well addressed. In an effort to understand dynamics of the hollow droplet impingement, a numerical study for the impact of a hollow droplet on a flat surface is presented. The impingement model considers the transient flow dynamics during impact and spreading of the droplet using the volume of fluid surface tracking method (VOF) coupled with the momentum transport model within a one-domain continuum formulation. The model is used to simulate the hydrodynamic behaviour of the impact of glycerin hollow droplet. It is found that the impact and spreading of the hollow droplet on a flat surface is distinctly different from the conventional dense droplet and has some new hydrodynamic features. A phenomenon of formation of a central counter jet of the liquid is predicted. With the help of simulations the cause of this phenomenon is discussed. Comparison of the predicted length of the central counter jet and the velocity of the counter jet front shows good agreements with the experimental data. The influence of the droplet initial impact velocity and the hollow droplet shell thickness on the impact behaviour is highlighted.     

液滴碰撞液膜润湿壁面空气夹带数值分析

采用复合水平集-流体体积法并综合考虑传热及接触热阻的作用, 对液滴碰撞液膜润湿壁面空气夹带现象进行了数值分析. 揭示了夹带空气形成机理, 探索了夹带空气特性参数随碰撞速度和液膜厚度的变化规律, 获得了夹带空气作用下液滴碰撞润湿壁面的传热机理. 研究结果表明: 撞壁前气液两相压力差是引起气液相界面拓扑结构变化以及夹带空气形成的主要原因; 液滴碰撞速度与压缩空气层内压力以及相界面形变高度密切相关; 液滴接触液膜时, 碰撞轴上液滴底部和液膜表面速度相等, 大约是碰撞速度的1/2; 碰撞速度对夹带空气层底部到破碎点的无量纲弧长和最大无量纲夹带空气直径均存在较大的影响; 液滴和液膜的无量纲形变高度与斯托克斯数密切相关; 液膜初始厚度对液滴和液膜的无量纲形变高度和最大无量纲夹带空气直径影响较大; 撞壁初始阶段, 碰撞中心区域夹带空气对壁面热流密度分布存在较大的影响.     

Estimating critical surface tension from droplet spreading area

Critical surface tension (CST) is a measure of solid surface tension and is mainly determined by measuring the contact angle of a droplet on a target solid surface. The concept of CST makes it possible to determine solid surface tension without any unprovable assumptions such as the Fowkes hypothesis. However, it requires somewhat special devices and skills for measuring the contact angle. In this work, we propose a simple method to determine the CST of a solid by measuring the droplet spreading area. This method is developed by combining the conventional CST with a simple analytical droplet model. The difference in estimated CSTs between our method and the conventional one is within 3.0%. Our method enables a quick and simple evaluation of the solid surface tension without special devices for measuring the contact angle.     

Experimental study of the effect of helium/nitrogen concentration and initial droplet diameter on nonane droplet combustion with minimal convection

In this paper, we study the influence of inert concentration and initial droplet diameter on nonane (C₉H₂₀) droplet combustion in an environment that promotes spherical droplet flames. The oxygen concentration is fixed while the inert is varied between nitrogen and helium. A range of initial droplet diameters (D_o) are examined in each ambient gas: 0.4 mm < D_o < 0.8 mm; and an oxidizing ambiance consisting of 30% oxygen (fixed) and 70% inert (fixed), with the inert in turn composed of mixtures of nitrogen and helium in concentrations of 0, 25, 50, 75, and 100% N₂. The experiments are carried out at normal atmospheric pressure in a cold ambiance (room temperature) under low gravity to minimize the influence of convection and promote spherical droplet flames. For burning within a helium inert (0% N₂), the droplet flames are entirely blue and there is no influence of initial droplet diameter on the local burning rate (K). With increasing dilution by nitrogen, droplet flames show significant yellow luminosity indicating the presence of soot and the individual burning histories show K reducing with increasing D_o. The evolution of droplet diameter D(t) is nonlinear for a given D_o in the presence of either helium or nitrogen inerts indicating that soot formation has little to do with nonlinear burning. A correlation is presented of the data in the form where the effective burning rate, K′, and ε are concentration-dependent. Correlations for these parameters are presented in the paper.     

Numerical simulation on partial coalescence of a droplet with different impact velocities

Partial coalescence is a complicated flow phenomenon. In the present study, the coalescence process is simulated with the volume of fluid(VOF) method. The numerical results reveal that a downward high-velocity region plays a significant role in partial coalescence. The high-velocity region pulls the droplet downward continuously which is an important factor for the droplet turning into a prolate shape and the final pinch-off. The shift from partial coalescence to full coalescence is explained based on the droplet shape before the pinch-off. With the droplet impact velocity increasing, the droplet shape will get close to a sphere before the pinch-off. When the shape gets close enough to a sphere, the partial coalescence shifts to full coalescence. The effect of film thickness on the coalescence process is also investigated. With large film thickness,partial coalescence happens, while with small film thickness, full coalescence happens. In addition, the results indicate that the critical droplet impact velocity increases with the increase of surface tension coefficient but decreases with the increase of viscosity and initial droplet diameter. And there is a maximum critical Weber number with the increase of surface tension coefficient and initial droplet diameter.     

钛表面液滴冻结的实验研究

基于可视化实验系统,研究了钛表面静止液滴的冻结行为,分析了底板倾斜角度对液滴的相变时间和接触直径的影响,以及液滴体积大小对相变时间的影响.结果 表明:底板倾斜角度对液滴相变时间的影响并不是单调的,随着底板倾斜角度增加相变时间呈现先减小再增加最后又减小的趋势,底板倾斜角度为0°以及35°时,相变时间最长;液滴的接触直径随...     

液滴撞击圆柱内表面的数值研究

针对液滴撞击圆柱内表面的过程,利用基于相场的格子Boltzmann方法模拟液滴以不同初速度、从不同初始高度、撞击不同大小的圆柱内表面时液滴的形态变化,分析了液滴自身物性(如密度和黏性等)和圆柱内表面润湿性等因素对撞击现象的具体影响.研究发现:撞击韦伯数、密度比及动力黏性比、圆柱半径等对液滴撞击后沿圆柱内表面的铺展均有一定影响,较高的韦伯数下液滴可能会发生分裂;液滴初始高度对大密度比和动力黏性比的撞击影响较小;液滴反弹现象可能出现在接触角较大时;重力作用会抑制撞击后液滴的振荡.     

Simulation solution for micro droplet impingementon a flat dry surface

This paper presents a computational fluid dynamics approach for micro droplet impacting on a flat dry surface. A two-phase flow approach is employed using FLUENT VOF multiphase model to calculate the flow distributions upon impact. The contact line velocity is tracked to calculate the dynamic contact angle through user defined function program. The study showed that the treatment of contact line velocity is crucial for the accurate prediction of droplet impacting on poor wettability surfaces. On the other hand, it has much less influence on the simulation of droplet impacting on good wettability surfaces. Good fit between simulation results and experimental data is obtained using this model.     

Water droplet impact on superhydrophobic surfaces with microstructures and hierarchical roughness

Quantitative correlation between the critical impact velocity of droplet and geometry of superhydrophobic surfaces with microstructures is systematically studied.Experimental data shows that the critical impact velocity induced wetting transition of droplet on the superhydrophobic surfaces is strongly determined by the perimeter of single micropillar,the space between the repeat pillars and the advancing contact angle of the sidewall of the micropillars.The proposed model agrees well with the experimental results,and clarifies that the underlying mechanism which is responsible for the superhydrophobic surface with hierarchical roughness could sustain a higher liquid pressure than the surfaces with microstructures.     

Microgravity experiments on droplet motion during flame spreading along three-fuel-droplet array

Flame spreading along a fuel droplet array at microgravity has been studied as a simple model of spray combustion. A three droplet array with a pendulum suspender was employed to investigate interactions between flame spreading and droplet motion in the array direction. Initial droplet diameter was 0.8 mm and fuel was n-heptane. A silicon carbide pendulum suspender of 15 μm in diameter and 30 mm in length was used for the third droplet. The first fixed droplet was ignited by electric spark. Behavior of the flame and the third droplet was observed using a high-speed video camera. Dimensionless span, which is the averaged droplet span divided by the averaged initial diameter of the three droplets, was varied from 2.7 to 10. Large displacement of the movable droplet was observed after group flame grew around the movable droplet. As the initial dimensionless span increased, the averaged droplet speed after the occurrence of flame spreading to the movable droplet increased steeply, taking the maximum value around 5 in initial dimensionless span, and then decreased gradually. The movable droplet advanced toward the second droplet in small spans and moved away from the second droplet in large spans. The direction of the motion changed around 4.6 in initial dimensionless span. Flame spread induction time from the second to the third droplet increased exponentially as the initial dimensionless span was increased. The induction time of flame spreading to a movable droplet was longer than that of flame spreading to a fixed droplet. From calculations of flame spreading along a 20-droplet array, it was predicted that the droplet speed nearly converged after flame spread to the sixteenth droplet. The maximum speed of the nineteenth droplet appeared around 7.5 in the initial dimensionless span.     

Binary collision of a burning droplet and a non-burning droplet of xylene: Outcome regimes and flame dynamics

The binary collisions of a burning droplet and a non-burning droplet of xylene are experimentally investigated. The experimental parameters span an extensive range of Weber number and impact parameter, covering the collision outcome regimes of coalescence, reflexive separation, and stretching separation. A high-speed camera captures the temporal details of the collision process, involving flame spread, visible radiation, and flame distributions around droplets. For reflexive separation and stretching separation, the flame from the droplet spreads to the ligament, surrounding it during the interaction process, and then spreads around separated droplets and satellite droplets. Highly-interactive flames are formed in-between the droplets, with very sooty flames generated for most collisions. For the coalescence case, a swirling flame forms around the rotating coalesced droplet. For similar Weber numbers, visible flame radiation is compared for different collision regimes. The visible flame radiation changes more significantly for the reflexive and stretching separation cases than it does for the coalescence case. The change of the averaged visible flame radiation for reflexive separation and stretching separation is more than two times higher than that for coalescence. The map of three different collision regimes is plotted in the Weber number versus impact parameter domain and compared with available theoretical model predictions. Although the different outcomes of collision with the presence of flame can be well predicted by the model, using fluid properties determined by the averaged properties of the two droplets, the dynamics of the detailed processes involved in the collisions are very interesting and have strong implications on overall combustion behavior that go well beyond the mapped regimes.