Droplet Impact Research Articles

Predicting Leidenfrost Temperature and Effects of Impact Conditions on Droplet Size Distribution in Film-Boiling Regime

Abstract The interaction of a droplet with a solid wall is relevant in various engineering applications. The properties of the resulting secondary droplets are determined by the wall temperature, ambient pressure, impact momentum, and impact angle. This paper presents a comprehensive characterization of drop-wall interactions and the subsequent atomization as a function of the combined effects of such parameters. A drop-wall interaction model is derived for F-24 droplets using smoothed particle hydrodynamics (SPH). The model can predict different impact outcome regimes (deposition, rebound, contact-splash, and film-splash) for different ambient pressures, wall temperatures, and impact parameters. The model also addresses the effect of ambient pressure on the Leidenfrost behavior. Size distributions of secondary droplets are compared for vertical and nonvertical impacts of F-24 droplets on superheated surfaces in the film-boiling regime. The non-dimensional Sauter mean diameter (SMD) of the secondary droplets varies based on the position in the impact plane for all the nonvertical impacts but remains almost unchanged for vertical impacts. The leading arc for nonvertical impact consists of larger secondary droplets, and the size decreases toward the trailing arc. An empirical relation is proposed to represent this trend. This research sheds light on successive droplet impacts by studying the effects of impact frequency on SMD evolution. The results are compared to single droplet impact cases for different fuels and Weber numbers. The size of secondary droplets for successive impacts is observed to be nearly indistinguishable from that of single-droplet vertical impacts.

Investigation on dynamic characteristics of droplet impact on surfaces with different rough microstructures based on lattice Boltzmann method

Constructing microstructures on smooth surfaces is significant for 3D printing technology and surface self-cleaning effects in engineering applications. This paper uses the lattice Boltzmann method (LBM) to study the droplet impact behavior on different microstructured surfaces: right-triangle, semicircular, and rectangular microstructures. The formulas for calculating the surface wetting characteristics of different microstructures are derived theoretically, the difference between simulation results and those from theoretical analyses is less than 5 %. The increase of the microstructure column height h may eventually lead to the disappearance of the Wenzel mode. The droplet impact on the microstructured hydrophilic surface can accelerate the droplet to the stable stage and reduce the excessive impact effect. The droplet impacts on the superhydrophobic microstructured surface can accelerate the rebound separation and shorten the time to reach the stable stage. When droplets impact the superhydrophobic surface of the right-triangle microstructure, a secondary separation phenomenon will occur, resulting in the maximum rebound height and the maximum reduction of the near-wall density. The superhydrophobic right-triangle microstructured surface can more effectively promote self-cleaning.

System Dynamics of the Water Droplet Impacting Flexible Cantilever Beams.

In this work, a water droplet impacting superhydrophobic flexible cantilever beams is systematically studied via experimental methods, aimed at recognizing the significance of the system dynamics that arises from the interplay between substrate oscillation and droplet impact. Influences of the substrate stiffness and the impact Weber number on the substrate oscillation and droplet impact dynamic are the focus particularly. For substrate oscillations, the beam deflection increases with the Weber number but decreases with the beam stiffness, while the oscillation period of the beam is not affected by the impact dynamic. For the droplet impact dynamic, the spreading dynamic is independent of beam oscillation, while the retraction dynamic is closely related to the surface elasticity. The effect of the cantilever beams on the droplet (i.e., promoting or inhibiting the rebound behavior) is dependent on the coupling movement of the water drop and the cantilever beam, which is varied by changing the stiffness of the cantilever beam. The findings of this work will provide a theoretical reference for the application of flexible substrates in the fields of anti-icing and self-cleaning.

Dynamic characteristics of droplet impact on a cold cylindrical surface

The freezing phenomenon often leads to economic and safety losses in various industries. Studying the processes of droplets impact on cold surfaces is beneficial in preventing potential harm. In this study, the processes of droplet impact on a cold cylindrical surface are investigated by high-speed photography. How the dynamic characteristics are affected by Weber number and surface temperature are studied. The results indicate that when the surface temperature is between −25 to −5 °C, an increase in Weber number results in higher initial kinetic energy and greater viscous dissipation. During the spreading process, the initial kinetic energy predominates, leading to an increase in the maximum spreading ratio. However, the rise in viscous dissipation contributes to an increase in the stable ratio. As Weber number increases, both the maximum spreading time and stable time decrease. Over the same range of surface temperatures, a decrease in surface temperature results in increased viscous dissipation, which leads to a reduction in the maximum spreading ratio. The maximum spreading time and stable time, however, are not significantly affected by surface temperature. Moreover, the relationships of the maximum spreading ratios and stable ratios between axial and circumferential directions are constructed respectively. These relationships exhibit similar linear characteristics, which are not affected by the Weber number and surface temperature.

Impact of a compound droplet on a curved surface: Effects of Weber and Reynolds numbers

Droplet impact on curved surfaces is a common phenomenon in both industrial processes and the natural world. Prediction and design of advanced technologies, such as cell processing and biotechnology printing, require an in-depth understanding of the dynamics of compound droplet impacts. However, there is a significant scarcity of research concerning the dynamics of compound droplets impacting on curved surfaces. This study employs numerical simulations to investigate the impact dynamics of a compound droplet on a curved surface over a range of Weber number (We) and Reynolds number (Re). The numerical approach employed here is an immersed boundary lattice Boltzmann method. The numerical simulation results reveal that both the We and Re have a pronounced effect on the spreading, contraction, and rupturing, of compound droplets on a curved surface. An analysis, with We and Re as key factors, elucidates their influences on the impact behavior of compound droplets. Furthermore, based on the phase diagram of impact results on the We and Re coordinate axes, two split lines are introduced for the first time, delineated by inertia, viscosity, and surface tension forces. These split lines provide clear guidance on controlling the droplet impact pattern by changing either We or Re, thereby facilitating accurate predictions of compound droplet impact results.

Impact and freezing characteristics of deionized water droplets on cold curved surfaces

Processes involving droplet impact and subsequent freezing occur widely in practical engineering applications. In the present study, a visualization experimental setup is utilized to investigate the effects of the impact of single millimeter-scale droplets on curved surfaces at room and low temperatures. The influences of the Weber number We, wall temperature, and wall wettability on the dynamics of droplet impact and the characteristics of ice formation are examined. The morphological evolution of droplet impact and the variations of the dimensionless spreading coefficient are analyzed. The results indicate that at high We (We = 277), droplets reach their maximum spread on cold walls in a shorter time than on room-temperature walls, and their peak spreading coefficient is smaller. Upon impact with a cold wall, droplets exhibit a spread–splatter behavior. Low temperatures suppress the oscillatory behavior of droplets on a curved wall. In the case of a hydrophilic wall surface, as the impact We increases from 42 to 277, the impact mode gradually transitions from spread–retract–freeze to spread–splatter–freeze. The maximum spreading coefficient first increases and then decreases with increasing impact We. At high We (We = 277), the wall wettability has a minimal effect on the dynamics of droplet impact and freezing, with a spread–splatter–freeze mode being exhibited for both hydrophobic and hydrophilic walls, and the final freezing morphology is similar.

Damage evolution simulation and lifetime prediction for composite blades under continuous droplet impacts

Offshore wind power generation is a promising technology in renewable energy applications due to its high reserves of wind energy in sea areas. To improve the energy transformation efficiency, the blade length of the offshore wind turbines has become larger and larger, and it has made rain erosion be one of the most frequent failures during the turbine operation. As the natural rainfall is stochastic in spatial and time domains, it is difficult to depict the damage evolution process caused by rain impact exactly. Therefore, regular and continuous droplet impact simulation and experiments present an alternative methodology for this issue. With the composite structure of the blade and deformability of the liquid droplet in consideration, a fluid solid interaction model will be established to investigate the impact response and subsequent damage evolution. In which, the Smooth Particle Hydrodynamics (SPH) model is utilized to depict the constitutive relationship within the droplet, and Finite Element Method (FEM) is used to construct the Representative Volume Element (RVE) model of the blade leading edge. The impact process is simulated first to obtain the impact pressure distribution at the contact center and velocity field in the droplet. Furthermore, the stress wave propagation in the blade multilayer structure can be analyzed. Owing to the multiaxial fatigue feature of the continuous droplet impact, the continuum damage mechanics is integrated with the fatigue criterion and the Jump-in-Cycle procedure is used to simulate the high-cycle fatigue process. The damage factor distribution on the blade coating surface and its influence on mechanical properties are analyzed. Thereafter, the droplet impact fatigue life can be accumulated based on Miner’s linear rules. The theoretical achievements are validated by experimental data provided by Rain Erosion Testing (RET), which shows a good agreement between each other. As a result, V-N curves and D-N curves, i.e. quantitative relationship between droplet falling conditions and impact fatigue life, are established. The achievements in this study can provide an effective tool for rain erosion mechanism analysis and life prediction in industrial applications.

Impact of a compound droplet on a solid surface: The effect of the shell on the core

The dynamic behavior of compound droplets impacting a solid surface was studied via experiments over κ (defined as the ratio of the compound droplet shell thickness h to the diameter D0 of compound droplet) ranging from 0 to 0.34, We ranging from 25 to 325 and Re ranging from 165.3 to 3405.2. The spreading diameter ratio, the maximum spreading dynamic contact angle and spreading speed of the core were investigated. Four modalities of the core of compound droplets were observed on the solid surface, including a) core rebound, b) no rebound, c) core splitting rebound, d) core splitting. The results revealed that the thickness of the shell, We, and the viscosity of the shell have a significant effect on the rebound and spreading processes of the core of the compound droplet. The high viscosity oil shell is conducive to its spreading. As the thickness of the oil shell increases, its cushioning effect on the water core also increases. In addition,κ=0.02254We0.503,κ=-0.336e-We121.056+0.319 and κ=0.06086e-We121.056+0.07716e-We121.70598+0.01595 were used to divide the modal boundary of the compound droplet core. Further analysis reveals the correlation between We, Re, βm and.κ,12+Weβm=8+31-cos22.78+84.57κβm3+0.955We1.05Reβm6.5.

Sound of natural and artificial rain impact on ETFE membranes

Thin membranes are known to produce a drum-like effect when subjected to rainfall. Rain noise affects acoustic comfort and might disrupt communication, speech intelligibility and cognitive task performance. Therefore, assessment of the characteristics of noise produced by impact of rain on membranes in different architectural settings is important. Given inherent disparities between natural and artificial precipitation-primarily in terms of droplet size distribution and terminal velocities-underscores the importance of conducting evaluations using natural rain noise. The focus of this study is on the comparative assessment of the sound generated by natural rainfall versus the sound generated by the impact of individual droplet in a semi-anechoic room on single Ethylene tetrafluoroethylene (ETFE) membranes. The study assesses the extent of spectral variations of the rain sound produced in both approaches and the corresponding absolute value of noise intensity.

Numerical simulations of freezing behaviors of water droplets impacting cold hydrophobic surfaces

Freezing, as a prominent phase transition phenomenon, occurs frequently in nature and industry. This study employs a combination of the Volume of Fluid(VOF) model, solidification model, and a dynamic contact angle model to simulate the dynamics and solidification process of droplets impacting supercooled surfaces. Two distinct freezing morphologies, central-concave and central-convex, are observed. These different morphologies are attributed to the competition between the time scales associated with droplet impact dynamics and solidification (retraction time tr and solidification time tf,p). Thus, we propose a scaling law to predict the occurrence of different morphologies: when tr < tf,p, the final freezing morphology exhibits a central-convex shape, while when tr > tf,p, it exhibits a central-concave shape. In order to explore effects of nucleation on solidification, the classical nucleation theory is adopted to determine the nucleation time, and thus the freezing behaviors can be well captured. Furthermore, we examine the effect of volume expansion during phase transition on the freezing behaviors, and find that the increase in freezing time due to volumetric expansion can be well predicted by the scaling law tf ∼ Δ(H2)/ΔT.

Analysis of droplet motion in cavity zone of rotating packed bed

Droplet motion in outer cavity zone of rotating packed bed is a complex phenomenon that requires analysis of liquid and gas flow to clarify influences of factors such as rotation speed, liquid flow rate, or gas flow rate. In this study, hydrodynamics is described using laser Doppler velocimetry and high-resolution visualization for droplet motion characterization. The results reveal that droplet motion initially depends on the peripheral velocity of the packing, while further from the packing, it loses kinetic energy mainly due to the drag force caused by the surrounding gas. Two-way coupling is investigated as both phases interact with each other. Rankine vortex theory is applied to simulate gas flow in the outer cavity zone. A model, based on Newton’s second law, is developed to predict droplet dynamics, including impingement velocity, trajectory length, and total residence time. Validation of the model shows that it can predict the liquid velocity in the outer cavity zone accurately. Upon impingement of droplets on the casing, part of the liquid rebounds as secondary droplets into the cavity, where they circulate with the gas flow. Overall, this study provides insights into droplet dynamics in the rotating packed bed and offers a basis for optimizing process efficiency and design in various applications, including chemical processing and environmental engineering.

Implementation of a high-frequency phosphor thermometry technique to study the heat transfer of a single droplet impingement

Contributing to a better understanding of spray cooling systems, the heat transfer process underlying the event of a droplet impinging onto a uniformly heated stainless steel surface (SS304) was experimentally investigated. Since the heat transfer process is linked to the droplet’s hydrodynamics, high-speed videos were recorded to measure the deformation of the droplet. A series of isothermal and non-isothermal impacts were performed for Weber numbers (We) within the range 17.7≤We≤58.2. A strong relationship between the maximum spreading ratio reached by the droplet and its initial kinetic energy was found. The surface temperature directly affects the droplet hydrodynamic during the impact by promoting an oscillatory behavior of the droplet after the maximum spreading is reached. Given the spatial–temporal resolution of the heat transfer process, a high-frequency phosphor thermometry technique was implemented, finding that the temperature drop upon droplet impact was independent of impact velocity. The sharp temperature drop results in an intense thermal interaction that occurred during the first 10 ms of the impact. The maximum average heat flux registered was 98.56 W/cm2 with a cooling effectiveness of 3.5%.

Investigation into the mechanism of ignition of magnesium alloy plates by high-temperature molten droplet

The ignition and spread of magnesium alloy fires are mainly attributed to the ignition of high-temperature molten metal droplets, but the underlying ignition mechanism remains unclear. This study addresses this knowledge gap by employing metallic aluminum and copper particles heated to temperatures over 1200 °C. Systematically, these heated particles were released at a fixed rate onto a magnesium alloy plate located inside a controlled heating furnace that maintained a constant temperature of 500 °C. Experimental results show that whether there are irregular pores and cracks in the resulting molten product can be used as a criterion for judging whether the magnesium alloy sample ignites. Furthermore, it is worth noting that molten copper droplets of the same size exhibit a higher heat-carrying capacity than their aluminum counterparts when both are heated to the same temperature. In addition, this study expanded the research scope through numerical simulation and analysis, using the PATO (Porous-Material Analysis Toolbox) solver to clarify the transient heat conduction process between the high-temperature droplet and the magnesium alloy plate. Simulations show that heat transfer upon droplet impact is primarily longitudinal, directed toward the base of the material. This results in a significantly higher temperature in the center of the affected magnesium alloy plate than at its periphery. In summary, this study helps to improve the understanding of the dangers of high-temperature molten droplets, and it is of great significance to investigate fire accidents caused by high-temperature molten droplets.

Numerical and theoretical modeling of water droplet impact on hydrophilic and superhydrophobic cones

The phenomenon of droplet impingement on solid surfaces is prevalent in various natural and industrial contexts. Research on impact dynamics on conical surfaces keeps emerging, with superhydrophobic cones receiving more attention than hydrophilic ones. This study systematically investigates water droplet impact dynamics on both hydrophilic and superhydrophobic cones using a two-phase numerical solver under different Weber numbers (We) and cone angles (φ). Three distinct phases are identified in the We–φ map to describe the different outcomes on each type of cones. Generally, deposition occurs ultimately on hydrophilic cones, whereas rebounding is observed on superhydrophobic ones. The maximum spreading area βAmax on hydrophilic cones depends only slightly on φ but consistently increases with We, following a scaling law of We0.5 at higher We. In contrast, on superhydrophobic cones, βAmax increases significantly with both We and φ, and the exponent in the scaling laws of βAmax with respect to We increases notably as φ increases. Three characteristic times are defined to describe important motion states on both types of cones. Corresponding scaling laws for each time with We are established. Two theoretical models are developed to predict the maximum spreading position for droplets on hydrophilic cones and the rebound position on superhydrophobic cones, respectively. Gravitational potential energy is included in the energy budget for both models, and an auxiliary viscous dissipation due to spontaneous spreading is accounted for the hydrophilic case. Satisfactory agreement between the theoretical and numerical results is achieved.

Investigation of droplet splashing behavior during oblique jet impact onto a wall

For the issue of jet impingement on the wall in industrial cooling processes, an experimental setup based on high-speed photography for oblique jet impingement onto the wall was constructed. The experimental focus was on the study of liquid droplet splashing behavior after oblique jet impingement on the wall, discussing the liquid droplet splashing behavior under three jet impingement modes: Rayleigh regime, first wind-induced regime, and secondary wind-induced regime. By employing methods such as trajectory imaging and particle image velocity for liquid droplet parameter image measurement, obtaining the particle size, velocity, and distribution of splashed droplets after oblique jet impact on walls under different working conditions. The impact of jet impingement velocity and angle on droplet splashing parameters was analyzed. The results showed that when the impingement point is before the breakup length, with increasing flow velocity, the surface wave of the liquid column, and the spreading liquid film became more pronounced, but the loss of liquid-phase components due to splashing was relatively small. When the impingement point is after the breakup length, the secondary breakup resulting in a “crown”-shaped liquid film after droplet impingement leads to a significant loss of liquid-phase components through splashing. As the inlet velocity of the jet increases, there is a decreasing trend in droplet size and an increasing trend in droplet velocity. With an increase in jet angle, there is a decreasing trend in droplet size and velocity. Based on the concentration, size, and velocity distribution characteristics of splashing droplets, the area after oblique jet impingement on the wall can be divided into the impingement zone, low-concentration low-velocity zone, high-concentration high-velocity zone, and lateral splashing zone. This has significant implications for understanding the splashing mechanism after oblique jet impingement on the wall and optimizing operating conditions.

Numerical simulation of a droplet impact onto a curved liquid film

Numerical simulation of a droplet impact onto a curved liquid film

Direct numerical simulation of droplet impact onto dry stationary and moving walls at low to high Weber numbers

The accurate prediction of the drop-wall impact dynamics and regime are of interest in many engineering applications. Considerable progress has been made in different elements of Direct Numerical Simulation (DNS) approaches. However, identifying the optimal combination of improved submodels that perform well under a wide range of operating conditions, covering low to high Weber numbers with varying surface wettabilities and wall motions, still requires further investigation. Here, we conduct a comparative study on the performance of key DNS elements, including the Interface Tracking (IT) algorithm for the widely-used volume-of-fluid approach, Contact Angle (CA), Contact-Line Velocity (CLV), Instability Actuation (IA), and numerical discretization. The study spans a broad range of Weber numbers, encompassing various impact phenomena such as drop oscillation, partial rebound, fingering, and splash. It is concluded that the Kistler CA model provides the most accurate predictions among the models considered here. In terms of IT, the Multidimensional Universal Limiter with Explicit Solution (MULES) algorithm is the most efficient one, especially for moderate to high Weber numbers. For high Weber numbers, involving the fingering instability and splash events, while applying the IA mechanism slightly improves the results, using a high-quality fine-enough grid and appropriate numerical discretization scheme to control the dispersion and dissipation errors are more important ingredients. It is shown that the recommended model combination is also able to predict the important features of drop impact on moving walls with reasonable accuracy.

Experimental study on the influence of surface properties on droplet collision dynamics, from adhesion to rebound to breakup

Liquid droplet impact on dry surfaces often results in bouncing or breakup beyond a certain threshold. Surface contact angles, especially dynamic ones present during impact, significantly affect this process. Our experimental study underscores that advancing and receding contact angles influence droplet behaviors like rebounding and different types of breakup. This discovery provides new insights and criteria for understanding liquid droplet impact on surfaces. Special characteristics were found in the breakup on microstructured surfaces: the size of fractured droplets notably decreases, and the spreading–breakup occurs more easily and earlier. Additionally, microstructured surfaces reduce contact time to some extent. Furthermore, the uniqueness of oblique impacts is mainly reflected in how they lower the threshold of the receding contact angle for rebound. Studying the correlations and differences in droplet rebound and breakup related to these surface characteristics will contribute to improving research on liquid–solid interactions and the design of hydrophobic surfaces, including microstructured surfaces.

Study of the Impact of Droplets on Superhydrophobic Surfaces

The potential applications of superhydrophobicity for self-cleaning, microfluidic systems, and so on, inspired by the “lotus leaf effect”, have recently attracted great attention from academia and industry. To date, however, neither experimental nor theoretical studies have reached the scalable application of superhydrophobic surfaces, and no simple methods have been developed to produce durable superhydrophobic surfaces. In this project, conical, cylindrical, and cube microstructures with different shapes and dimensions were prepared by 3D printing. Among these structures, the conical microstructures without further modification exhibited the best hydrophobic performance with a contact angle of 120°. However, the hydrophobic performance of the cube microstructures was significantly improved with the addition of an FAS coating to the surface, reaching 131°. This indicates that, with the help of the coating, the 3D printed materials with special features can be transformed into hydrophobic surfaces. The outcome of this project will be useful in promoting the efficient and low-cost preparation of superhydrophobic functional surfaces and their engineering applications.

Dynamics of impact of polymer droplets on viscoelastic surfaces

Dynamics of impact of polymer droplets on viscoelastic surfaces