纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.      

纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.     

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

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

液滴撞击具有浸润性分布的倾斜固壁铺展行为的格子Boltzmann研究

采用基于单组分多相伪势模型的格子Boltzmann方法,模拟了三维液滴撞击左右两侧浸润性不同的倾斜固壁的铺展过程,获得了液滴在壁面两侧的铺展因子、相对铺展宽度、相对高度和液滴运动速度随时间的变化情况,研究了壁面浸润性分布和壁面倾斜角度对液滴铺展过程的影响.结果表明,液滴在倾斜壁面的铺展过程受到重力和表面力的综合作用,重力影响液滴的铺展和沿壁面向下的滑动,壁面浸润性分布影响液滴向壁面亲水侧横向移动.     

平壁液滴静态铺展影响因素的研究

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

STATIC SPREADING OF DROPLET IMPACT ON SOLID SURFACE：INFLUENCE FACTOR

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

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

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

中性润湿平板上液膜的惯性收缩

液滴撞击平板的动力学机理研究具有重要的理论与工程价值, 对该过程中液滴的形貌变化及主要影响因素的研究是科学技术界关注的重点之一. 液滴在高速撞击平板达到最大铺展半径以后发生毛细-惯性收缩, 收缩速率满足类Taylor-Culick公式. 结合实验与有限元方法, 对平板上铺展液滴的收缩过程进行了研究. 结果表明, 在中性润湿(接触角约为90°)平板上液滴/液膜的收缩在经历上述收缩过程以后, 会有一个慢匀速收缩过程, 速度约为第一阶段收缩速度的1/10. 对后一阶段的撞击参数影响测试显示, 该收缩过程主要与液体密度、液膜初始形状有关; 而与液体黏性、壁面润湿条件等无关——其仍然是一种毛细-惯性机制主导的液面演化行为, 类似于液体射流的Rayleigh-Plateau失稳. 尽管黏性效应对于液滴撞击的铺展行为有明显影响, 但上述结论在10倍黏性的液体测试中仍然成立. 本研究可以为液滴反弹机理研究和相关工艺控制提供参考.      

用格子Boltzmann方法模拟液滴撞击固壁动力学行为

首次用格子Boltzmann方法中的伪势模型对液滴撞击固壁的动力学行为进行了数值模拟.详细研究了液滴在壁面上的流动状态以及各种因素对撞击过程的影响.通过数值模拟得到:壁面的可润湿性越小,液滴越容易发生反弹,液滴的回缩速度越快;液滴的撞击速度越大,所得到的相对直径越大,回缩速度越快;液滴的粘性越小,所得到的相对直径越大;液滴的表面张力越大,液滴越容易发生反弹现象.另外,液滴的最大相对直径与We数满足一定的线性关系,这些结果与前人的理论预测和实验结果完全吻合.     

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

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

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

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

Simultaneous droplet impingement dynamics and heat transfer on nano-structured surfaces

This study examines the hydrodynamics and temperature characteristics of distilled deionized water droplets impinging on smooth and nano-structured surfaces using high speed (HS) and infrared (IR) imaging at We = 23.6 and Re = 1593, both based on initial drop impingement parameters. Results for a smooth and nano-structured surface for a range of surface temperatures are compared. Droplet impact velocity, transient spreading diameter and dynamic contact angle are measured. The near surface average droplet fluid temperatures are evaluated for conditions of evaporative cooling and boiling. Also included are surface temperature results using a gold layered IR opaque surface on silicon. Four stages of the impingement process are identified: impact, boiling, near constant surface diameter evaporation, and final dry-out. For the boiling conditions there is initial nucleation followed by severe boiling, then near constant diameter evaporation resulting in shrinking of the droplet height. When a critical contact angle is reached during evaporation the droplet rapidly retracts to a smaller diameter reducing the contact area with the surface. This continues as a sequence of retractions until final dry out. The basic trends are the same for all surfaces, but the nano-structured surface has a lower dissipated energy during impact and enhances the heat transfer for evaporative cooling with a 20% shorter time to achieve final dry out.     

Spreading dynamics of droplet on an inclined surface

A three-dimensional unsteady theoretical model of droplet spreading process on an inclined surface is developed and numerically analyzed to investigate the droplet spreading dynamics via the lattice Boltzmann simulation. The contact line motion and morphology evolution for the droplet spreading on an inclined surface, which are, respectively, represented by the advancing/receding spreading factor and droplet wetted length, are evaluated and analyzed. The effects of surface wettability and inclination on the droplet spreading behaviors are examined. The results indicate that, dominated by gravity and capillarity, the droplet experiences a complex asymmetric deformation and sliding motion after the droplet comes into contact with the inclined surfaces. The droplet firstly deforms near the solid surface and mainly exhibits a radial expansion flow in the start-up stage. An evident sliding-down motion along the inclination is observed in the middle stage. And the surface-tension-driven retraction occurs during the retract stage. Increases in inclination angle and equilibrium contact angle lead to a faster droplet motion and a smaller wetted area. In addition, increases in equilibrium contact angle lead to a shorter duration time of the middle stage and an earlier entry into the retract stage.     

Molten droplet solidification and substrate remelting in microcasting Part I: numerical modeling and experimental verification

This paper presents a numerical model of a molten metal droplet impinging, solidifying and bonding to a solid substrate. The physical and numerical model includes dissimilar materials, multi-dimensional axisymmetric heat transfer, tracking of solid/liquid interfaces during remelting and solidification, and coupled treatment of the continuous droplet/substrate region. The numerical model solves for the evolution of the temperature distribution in the droplet and substrate, predicts the position of the remelting and solidification fronts, and accounts for convective motion. The effect of the convection induced by the droplet spreading is modeled through a time-dependent effective thermal conductivity. High-speed filming of the molten droplet impinging and spreading on the substrate is performed to obtain the required parameters to determine this time dependent effective conductivity. The accuracy of the model is investigated with experimental techniques. This research is directly related to the development of microcasting Shape Deposition Manufacturing (SDM) which is a process for automatically fabricating complex multi-material objects by sequentially depositing material layers. Microcasting is a molten metal droplet deposition process in SDM, which is able to create fully dense metal layers with controlled microstructure. Important issues in the production of high quality objects manufactured with microcasting SDM are: attainment of interlayer metallurgical bonding through substrate remelting, control of both substrate and droplet cooling rates, and minimization of residual thermal stresses. To validate experimentally the numerical modeling approach, predicted cooling rates are compared with thermocouple measurements and substrate remelting depths are verified through optical metallographic techniques. Received on 10 June 1998     

A numerical study of droplet motion/departure on condensation of mixture vapor using lattice Boltzmann method

Droplet motion/departure, which is governed by external force acceleration coefficient, droplet radius and surface wettability on solid surfaces under external forces such as gravitational force, play a significant role in characterizing condensation heat transfer, especially when high fractional non-condensable gases (NCG) present. However, due to the challenge in visualizing the vapor/steam velocity field imposed by droplet motion/departure, the detailed mechanism of droplet motion/departure on condensing surfaces has not been completely investigated experimentally. In this study, droplet motion/departures on solid surfaces under external forces and their interactions with steam flow are simulated using two dimensional (2D) multiphase lattice Boltzmann method (LBM). Large external force acceleration coefficient, droplet radius and contact angle, lead to large droplet deformation and high motion/departure velocity, which significantly shortens the droplet residual time on the solid surface. Our simulation shows that steam vortices (lateral velocity) induced by droplet motion/departure can greatly disturb the vapor flow and would be intensified by increasing external force acceleration coefficient, droplet radius, and contact angle. In addition, the location of vortex center shifts in the ascending direction with increase of these factors. The average lateral velocities induced by droplet motion/departure at various conditions are obtained. The mass transfer resistance is substantially reduced owing to the droplet motion/departure, leading to an enhanced heat flux. The experimental results are compared to validate the influence of droplet motion/departure on condensation heat transfer performance, especially for steam–air mixture with the presence of high fractional NCG.     

Lattice Boltzmann simulation of multiple droplets impingement and coalescence in an inkjet-printed line patterning process

Three-dimensional computations on the basis of the index-function lattice Boltzmann method are performed to simulate the process of multiple droplets impinging and coalescing into a line pattern on a solid substrate. The employed calculation model is validated by theoretical calculated values and experimental data from the literature. The influences of the equilibrium contact angle, droplet spacing and impinging velocity on the droplets impingement and coalescence behaviours are investigated. Numerical results demonstrate the width of the formed line depends significantly on the equilibrium contact angle and droplet spacing. The droplet spacing plays a significant role in controlling the coalescence moment of multiple droplets. The resolution of the printed pattern can be slightly increased with increase in impinging velocity.     

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.     

An experiment on the breakup of impinging droplets on a hot surface

Breakup characteristics of liquid droplets impinging on a hot surface are investigated experimentally with the wall temperatures in the Leidenfrost temperature range of 220–330°C for n-decane fuel. Factors influencing droplet breakup are wall temperature, impinging velocity, droplet diameter and impinging angle. The 50% breakup probability shows that the impinging velocity decreases linearly with the droplet diameter increase and there exists an optimum impinging angle near 80° having the minimum value in the impinging velocity for given wall temperature and droplet size. Near the wall temperature of 250°C corresponding to the Leidenfrost temperature, a peculiar nonlinear behavior in the breakup probability is observed.This work was supported by the Turbo and Power Machinery Research Center, Seoul National University.