EFFECT OF TWO-DIMENSIONAL MICROPILLAR ARRAYED TOPOGRAPHY ON SPREADING OF INSOLUBLE SURFACTANT-LADEN DROPLET

针对二维微柱阵列壁面上含不溶性活性剂液滴的铺展过程,采用润滑理论建立了液膜厚度和浓度演化模型,采用数值计算方法得到了液滴的铺展特征及相关参数的影响. 研究表明：活性剂液滴在微柱阵列壁面上铺展时,在壁面凸起处衍生出隆起结构,壁面凹槽处衍生出凹陷结构,随时间持续,隆起和凹陷均向两侧移动,且数量不断增加. 活性剂液膜流经凸起时,隆起高度呈驼峰形变化. 增大预置液膜厚度或活性剂初始浓度,铺展区域隆起和凹陷数量增多,液滴铺展速度加快. 增加凹槽深度或减小斜度会使毛细力作用增强,液膜破断可能性加大;增大凹槽宽度可加速活性剂液滴的铺展,加剧液膜表面波动幅度.     

Impact of liquid drops on a rough surface comprising microgrooves

Impact of water drops on a stainless steel surface comprising rectangular shaped parallel grooves is studied experimentally. Geometric parameters of the surface groove structure such as groove depth, groove width and solid pillar width separating any two successive grooves were kept at 7.5, 136 and 66 μm, respectively. The study was confined to the impact of drops in inertia dominated flow regime with Weber number in the range 15–257. Experimental results of drop impact process obtained for the grooved surface were compared with those obtained for a smooth surface to elucidate the influence of surface grooves on the impact process. The grooves definitely influence both spreading and receding processes of impacting liquid drops. A more striking observation from this study is that the receding process of impacting liquid drops is dramatically changed by the groove structure for all droplet Weber number.     

镓铟锡液滴撞击泡沫金属表面的运动学特性研究

常温下为液态的镓铟锡合金以其优异的导热性能在具有特殊要求的传热领域有着重要的应用价值,与传统流动介质相比较大的表面张力使得其产生的流动现象必有所区别.本文研究镓铟锡所形成的液滴撞击泡沫金属表面后所产生的铺展、回缩及回弹现象.采用高速相机拍摄液滴投影轮廓随液滴运动的变化过程,并通过图像处理获得不同撞击速度、底板表面孔径下的液滴铺展系数、中心位置轮廓高度以及液滴回弹后在空中的振动特性.研究结果表明:具有较高表面张力的镓铟锡液滴的铺展系数随无量纲时间的变化在铺展初始阶段仍满足常规流体的1/2次幂关系,只在铺展后期与底板的无量纲孔径有关系;液滴的最大铺展系数在较小无量纲孔径底板大于在光滑镍板,且随底板无量纲孔径增大而逐渐减小;在回弹过程,由于底板孔隙结构的存在使得液滴回弹后在空中的振动呈现3种形态:规则的横向和纵向振动、带旋转的横向和纵向振动以及旋转振动;最后,通过对振动频率的拟合和分析,进一步拓展了传统振动频率理论公式在非规则振动过程预测中的应用.      

Effect of the contact-line dynamics on the oscillations of a compressed droplet

This paper studies the natural and forced oscillations of a deformed droplet of an inviscid liquid surrounded by a different liquid and bounded in the axial direction by solid planes. In equilibrium, the droplet is a figure of revolution and the ratio of its radius to height is significant. The equilibrium contact angle between the side surface of the droplet and the solid surface is different from a right angle. The motion of the contact line is taken into account by setting an effective boundary condition. It is shown that three characteristic ranges of natural frequencies exist.     

EXPERIMENTAL AND THEORETICAL INVESTIGATIONS ON THE VELOCITY OSCILLATIONS OF DYNAMIC CRACK PROPAGATING IN BRITTLE MATERIAL TENSION

设计了一种精确测量裂纹在预拉伸带板中的传播速度的实验装置, 系统研究了裂纹在受不同预加载荷作用的不同尺寸的有机玻璃 板中的传播行为.实验结果表明:在一个细长、上下表面固定的带板中, 裂纹速度分为加速和稳定2个阶段.在稳定速度阶段, 裂纹平均速度v₀ 是预加载荷在板中储存的弹性能Wh的增函数.由于此时裂纹处于自相似传播阶段, 可以认为材料的动态断裂能G_c是裂纹传播速度v₀的增函数, 即材料存在"速度增韧性".当裂纹传播速度达到一定阈值, 裂纹的传播速度出现明显的振荡现象, 裂纹的振荡周期与断裂面出现的周期性沟槽的尺度一致.在更高的预加能量下, 裂纹传播路径出现弯曲、微分岔以及宏观分岔现象.建立了一个描述裂纹在带板中的直线传播行为的动力学模型, 并用这个模型模拟裂纹在传播过程中的速度振荡现象.     

STATIC SPREADING OF DROPLET IMPACT ON SOLID SURFACE：INFLUENCE FACTOR

通过建立液滴撞击固体平壁的静态铺展力学平衡的数学模型,从理论上得到了静态铺展半径与液滴物性参数、以及液滴与固体壁面接触角之间关系的数学表达式,将理论结果与数值模拟的结果进行了比较,两者吻合较好.比较了不同条件下液滴的静态铺展半径的变化规律,分别得到了液滴密度、体积、表面张力和接触角等因素对液滴静态铺展半径的影响规律.     

On dual-grid level-set method for contact line modeling during impact of a droplet on hydrophobic and superhydrophobic surfaces

In this paper, a numerical methodology for modeling contact line motion in a dual-grid level-set method (DGLSM) – solved on a uniform grid for interface which is twice that for the flow equations – is presented. A quasi-dynamic contact angle model – based on experimental inputs – is implemented to model the dynamic wetting of a droplet, impacting on a hydrophobic or a superhydrophobic surface. High-speed visualization experiments are also presented for the impact of a water droplet on hydrophobic surfaces, with non-bouncing at smaller and bouncing at larger impact velocity. The experimental results for temporal variation of the droplet shapes, wetted-diameter and maximum height of the droplet match very well with the DGLSM based numerical results. The validation of the numerical results is also presented with already published experimental results, for the non-bouncing on a hydrophobic and bouncing on a superhydrophobic surface, at a constant impact velocity. Finally, a qualitative as well as quantitative performance of the DGLSM as compared to the traditional level set method (LSM) is presented by considering our experimental results. The accuracy of the partially refined DGLSM is close to that of the fine-grid based LSM, at a computation cost which is close to that of the coarse-grid based LSM. The DGLSM is demonstrated as an improved LSM for the computational multi-fluid dynamics (CMFD) simulations involving contact line motion.     

Droplet impact dynamics for two liquids impinging on anisotropic superhydrophobic surfaces

Droplet impingement experiments were performed on grooved hydrophobic surfaces with cavity fractions of 0, 80, and 93?% using droplets of water and a 50?%/50?% water/glycerol mixture. The influence of liquid viscosity, cavity fraction, and spreading direction, relative to the surface grooves, is explored qualitatively and quantitatively. The maximum droplet spread diameter, velocity of the rebounding jet, and the time delay between droplet impact and jet emission were characterized for Weber numbers, We, based on droplet impact speed and diameter, up to 500. The unequal shear stresses and contact angles influence the maximum spread diameters in the two primary spread directions. At We?>?100, the ratio of the spread diameter along the direction of the grooves to the spread diameter perpendicular to the grooves increases above unity with increasing We. The maximum droplet spread diameter is compared to recent predictive models, and the data reveal differing behavior for the two fluids considered. The results also reveal the existence of very high relative jet velocities in the range 5????We????15 for water droplets, while such jets were not observed for the more viscous mixture. Further, in the range 115????We????265, the water/glycerol jet formation dynamics are radically different from the water behavior. Most evident is the existence of two-pronged jets, which arise from the anisotropy of the surface and the unequal shear stresses and contact angles that prevail on the surfaces. It is these influences that give rise to differences in the maximum spread diameters in the two primary spread directions. Similar two-pronged jet emission was observed for water over the very narrow range of We from 91 to 96. The issuing jet velocities were also observed to increase with increasing cavity fraction for both fluids and over the entire range of We explored. Lastly, the elapsed time between droplet impact and jet emission decreased with increasing cavity fraction.     

Oscillations of the fluid flow and the free surface in a cavity with a submerged bifurcated nozzle

The free surface dynamics and sub-surface flow behavior in a thin (height and width much larger than thickness), liquid filled, rectangular cavity with a submerged bifurcated nozzle were investigated using free surface visualization and particle image velocimetry (PIV). Three regimes in the free surface behavior were identified, depending on nozzle depth and inlet velocity. For small nozzle depths, an irregular free surface is observed without clear periodicities. For intermediate nozzle depths and sufficiently high inlet velocities, natural mode oscillations consistent with gravity waves are present, while at large nozzle depths long term self-sustained asymmetric oscillations occur.For the latter case, time-resolved PIV measurements of the flow below the free surface indicated a strong oscillation of the direction with which each of the two jets issue from the nozzle. The frequency of the jet oscillation is identical to the free surface oscillation frequency. The two jets oscillate in anti-phase, causing the asymmetric free surface oscillation. The jets interact through a cross-flow in the gaps between the inlet channel and the front and back walls of the cavity.     

An icing physics study by using lifetime-based molecular tagging thermometry technique

An experimental investigation was conducted to quantify the unsteady heat transfer and phase changing process within small icing water droplets in order to elucidate underlying physics to improve our understanding of the important micro-physical process of icing phenomena. A novel, lifetime-based molecular tagging thermometry (MTT) technique was developed and implemented to achieve temporally-and-spatially resolved temperature distribution measurements to reveal the time evolution of the unsteady heat transfer and dynamic phase changing process within micro-sized water droplets in the course of icing process. It was found that, after a water droplet impinged onto a frozen cold surface, the liquid water at the bottom of the droplet would be frozen and turned to solid ice rapidly, while the upper portion of the droplet was still in liquid state. As the time goes by, the interface between the liquid phase water and solid phase ice was found to move upward continuously with more and more liquid water within the droplet turned to solid ice. Interestingly, the averaged temperature of the remaining liquid water within the small icing droplet was found to increase, rather than decrease, continuously in the course of icing process. The temperature increase of the remaining liquid water is believed to be due to the heat release of the latent heat during solidification process. The volume expansion of the water droplet during the icing process was found to be mainly upward to cause droplet height growth rather than radial to enlarge the contact area of the droplet on the test plate. As a result, the spherical-cap-shaped water droplet was found to turn to a prolate-spheroid-shaped ice crystal with cusp-like top at the end of the icing process. The required freezing time for the water droplets to turn to ice crystals completely was found to depend on the surface temperature of the test plate strongly, which would decrease exponentially as the surface temperature of the frozen cold test plate decreases.     

Droplet heat transfer on micro-post arrays: Effect of droplet size on droplet thermal characteristics

Heat transfer towards a water droplet from hydrophobic micro-post array surface is considered while mimicking the environmental temperatures. Micro-post arrays are created on a silicon wafer surface via lithography technique. The textured surfaces are replicated by polydimethylsiloxane (PDMS) to achieve an optical transmittance. The droplet adhesion on micro-post array surface is presented and the influence of droplet size on the heat transfer and droplet internal flow characteristics is examined. The flow predictions are validated via the particle image velocimetry data. It is found that adhesion force between the water droplet and the micro-post arrays surface depends on the geometric size and the orientation of the micro-post arrays on the surface. Temperature and flow fields are influenced by the droplet size. The Nusselt and the Bond numbers increase with the droplet volume; however, the Bond number remains less than unity indicating that the Marangoni current dominates over the buoyancy current in the droplet. The Nusselt number attains larger values for micro-post array surface than that of the plain surface. This is because of temperature and velocity oscillations along the contact lines at the droplet bottom due to the pitches of the micro-post arrays.     

Axisymmetric oscillations of a cylindrical droplet with a moving contact line

Forced oscillations of a cylindrical droplet of an inviscid liquid surrounded by another liquid and bounded in the axial direction by rigid planes are investigated. The system is affected by vibrations whose force is directed parallel to the axis of symmetry of the droplet. The velocity of motion of the contact line is proportional to the deviation of the contact angle from the value at which the droplet is in equilibrium. Linear and nonlinear oscillations are considered. The conditions of the occurrence of resonance are determined.     

MOVEMENT AND SPLASHING OF A DROPLET IMPACTING ON A WET WALL

壁面液体层的存在对液滴撞击壁面的运动具有重要的影响。采用气液两相流动相界面追踪的水平集和流体体积复合方法和壁面润湿模型，实现了液滴撞击湿润壁面运动的数值求解；在此基础上，开展了液滴撞击湿润壁面运动的研究。研究结果表明：液滴以不同速度撞击湿润壁面时，会呈现出黏附铺展、波动运动、皇冠几何体运动以及飞溅运动等几种不同的运动形态，液滴撞击湿润壁面后的压力分布是不同运动形态形成的主要原因；飞溅运动是一定条件下皇冠几何体运动的一种特殊形态，液滴从皇冠几何体侧壁顶端的飞溅分离满足毛细破碎理论；撞击速度对分离液滴的运动方向影响较小，而对壁面液体层厚度的影响则较大；撞击速度和壁面液体层厚度对分离液滴形态、飞溅分离位置、飞溅速度以及飞溅发生时刻等都具有一定的影响。     

二维微柱阵列壁面对活性剂液滴铺展的影响

针对二维微柱阵列壁面上含不溶性活性剂液滴的铺展过程,采用润滑理论建立了液膜厚度和浓度演化模型,采用数值计算方法得到了液滴的铺展特征及相关参数的影响. 研究表明:活性剂液滴在微柱阵列壁面上铺展时,在壁面凸起处衍生出隆起结构,壁面凹槽处衍生出凹陷结构,随时间持续,隆起和凹陷均向两侧移动,且数量不断增加. 活性剂液膜流经凸起时,隆起高度呈驼峰形变化. 增大预置液膜厚度或活性剂初始浓度,铺展区域隆起和凹陷数量增多,液滴铺展速度加快. 增加凹槽深度或减小斜度会使毛细力作用增强,液膜破断可能性加大;增大凹槽宽度可加速活性剂液滴的铺展,加剧液膜表面波动幅度.      

Effect of the contact-line dynamics on the natural oscillations of a cylindrical droplet

The natural oscillations of a cylindrical droplet of an inviscid liquid surrounded by a different liquid and bounded in the axial direction by solid planes are studied. The motion of the contact line is described using an effective boundary condition. The dependence of the frequency and damping ratio on the capillary parameter is found. It is shown that the fundamental frequency of the translation mode vanishes beginning from a certain value of the capillary parameter. Depending on the ratio of the radial and axial dimensions of the droplet, the fundamental frequency of the axisymmetric mode and modes higher than the translation mode can vanish in a certain range of the capillary parameter. This dependence of the natural oscillation frequencies on the problem parameters allows one to determine the capillary parameter. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 78–86, September–October, 2007.     

SHOCK WAVE LOADING OF CONTACT SURFACE DETONATION MODE

基于叠加原理,利用炸药接触面起爆方式,实现对可压缩固体材料的矩形波加载. 通过联立爆炸产物和可压缩固体的（u,p）曲线以及综合考虑炸药与材料试件各自由边所受稀疏波干扰情况,从理论上给出冲击波平台压力和冲击波加载平台宽度的计算方法,为将来利用接触爆炸实现对材料试件的近似矩形波加载提供理论基础.     

粗糙表面上微液滴的数值模拟

本文使用VOF方法将微液滴在粗糙壁面上的接触现象转化为不可压缩两相流动问题,并对其进行三维数值模拟.选择具有柱形突起和槽道两种微结构的壁面进行模拟.计算了不同粗糙系数时液滴在突起结构表面的静止形态和接触角,计算结果和实验数据吻合得较好.和理论模型进行比较,分析了经典模型的适用范围.对于微槽道结构的壁面,计算给出不同方向测量得到的液滴接触角.实现了液滴在倾斜壁面上滑落过程的模拟.液滴沿斜面下滑时,前进角和后退角的变化存在周期性,这一周期性变化和表面粗糙结构密切相关.     

Spark-generated bubble collapse near or inside a circular aperture and the ensuing vortex ring and droplet formation

This paper aims to study the oscillation of a sparkgenerated submerged bubble located near or inside a circular aperture made in a flat plate using high-speed visualization technique. In the case of a bubble oscillating near an aperture the initial free surface of the water was set at the bottom surface of the plate. The effects of aperture size and bubblefree surface distance on the bubble behavior as well as on the ensuing droplet dynamics are investigated. It was found that the direction of the bubble reentrant jet was towards the aperture or away from it respectively when the normalized aperture size was smaller or greater than a certain critical value. In addition, a toroidal vortex ring was observed to form, which rotated inwards as it moved away from the aperture. It was also found that if the bubble was incepted at a distance sufficiently away from a supercritical size aperture a single droplet could be produced. In the case of a bub- ble initiated in the middle of a circular aperture submerged just beneath the water free surface, the bubble was found to take the shape of an ellipsoid during its expansion. Then a reentrant jet was initiated and pierced the bubble from its top side.     

THE OSCILLATION MODE OF INCOMPRESSIBLE CAVITY FLOW

采用二维大涡模拟方法进行了空腔水流场的数值计算, 考察空腔前缘动量损失厚度及来流速度等因素如何影响空腔流的振荡, 同时考察了空腔长深比与空腔流振荡模式的关系. 用空腔水流场的粒子图像测速测量结果验证了数值计算的可信性.结果表明, 空腔水流是否发生振荡取决于壁面摩擦速度.瞬时涡结构和空腔阻力系数2个方面的特征显示空腔水流场有2种典型的振荡模式, 剪切层模式与尾流模式, 确定振荡模式的关键因素是空腔长深比.     

Stretchable hydrophobic surfaces and droplet heating

Silicon elastomer surface is treated towards achieving the hydrophobic state. Functionalized nano-silica units are coated onto elastomer surface and resulting texture characteristics are examined prior to stretching, stretched and after stretching. The droplet heating of the hydrophobic elastomer surface is carried out when the surface is subjected to unstretching, stretching and stretch releasing conditions. The thermal-flow field in the liquid is simulated and validated incorporating high speed recording system. Nano-size silica units coated elastomer surface demonstrates the hydrophobic wetting state. The hydrophobic wetting state changes slightly for stretched and stretched released surface. The contact angle is about 154° ± 2° for unstretched surface while it is 152° ± 2° for the stretched released surface; hence, stretch relaxing provides reversible change of the surface wetting state of the elastomer surface. The contact angle reduces to 142° ± 2° when surface is under stretched, which is related to increased pillar spacing on the surface. The droplet heating results in development of Marangoni current in the fluid, which significantly affects the flow and temperature fields and it becomes more apparent for the large size droplets. The maximum flow velocity increases almost 9% in 45 µL as the surface is stretched. The Nusselt number increases with droplet size and the Bond number has the values less than unity; hence, stretching increases the Nusselt number by 60% for droplet of 45 µL.