Mesoscopic simulation of a thinning liquid bridge using the dissipative particle dynamics method

Rupture process of liquid bridges: The effects of thermal fluctuations

Rupture of a liquid bridge is a complex dynamic process, which has attracted much attention over several decades. We numerically investigated the effects of the thermal fluctuations on the rupture process of liquid bridges by using a particle-based method know as many-body dissipative particle dynamics. After providing a comparison of growth rate with the classical linear stability theory, the complete process of thinning liquid bridges is captured. The transitions among the inertial regime (I), the viscous regime (V), and the viscous-inertial regime (VI) with different liquid properties are found in agreement with previous work. A detailed description of the thermal fluctuation regime (TF) and another regime, named the breakup regime, are proposed in the present study. The full trajectories of thinning liquid bridges are summarized as I→V→VI→TF→ breakup for low-Oh liquids and V→I→ Intermediate →V→VI→TF→ breakup for high-Oh liquids, respectively. Moreover, the effects of the thermal fluctuations on the formation of satellite drops are also investigated. The distance between the peaks of axial velocity is believed to play an important role in forming satellite drops. The strong thermal fluctuations smooth the distribution of axial velocity and change the liquid bridge shape into a double cone without generating satellite drops for low-Oh liquids, while for high-Oh liquids, this distance is extended and a large satellite drop is formed after the breakup of the liquid filament occurs on both ends, which might be due to strong thermal fluctuations. This work can provide insights on the rupture mechanism of liquid bridges and be helpful for designing superfine nanoprinting.

Application of the dissipative particle dynamics method to the instability problem of a liquid thread

We investigate the application of the dissipative particle dynamics method to the instability problem of a long liquid thread surrounded by another fluid. The dispersion curves obtained from simulations are compared with classic theoretical predictions. The results from standard dissipative particle dynamics (DPD) simulations at first have a tendency of gradually approaching to Tomotika's Stokes flow prediction when the Reynolds number is decreased. But they then abnormally deviate again when the viscosity is very large. The same phenomenon is also confirmed in droplet retraction simulations when also compared with theoretical Stokes flow results. On the other hand, when a hard-core DPD model is used, with the decrease of the Reynolds number the simulation results did finally approach Tomotika's predictions when Re≈0.1. A combined presentation of the hard-core DPD results and the standard DPD results, excluding the abnormal ones, demonstrates that they are approximately on a continuum when labeled with Reynolds number. These results suggest that the standard DPD method is a suitable method for investigation of the instability problem of immersed liquid thread in the inertioviscous regime (0.1<Re<10), which is relevant for microfluidics applications but there is currently no theory. It cannot reach the Re≈0.1 regime because of some irregular properties of highly viscous DPD fluid, while the hard-core DPD is suitable to overcome this inferiority with standard DPD.

Effect of Soluble Surfactants on the Kinetics of Thinning of Liquid Bridges during Drops Formation and on Size of Satellite Droplets

The results of an experimental study on thinning and breakage of liquid bridges during detachment of a drop from the tip of a capillary are presented for a series of surfactant solutions (including cationic, anionic, and nonionic surfactants) over a broad range of molecular masses, values of critical micelle concentration, and concentrations. The used experimental protocol revealed that the kinetics of the bridge thinning depends much more on the dynamics of adsorption at the surface of the drop before it destabilizes, rather than on the depletion of surfactant from the surface of the thinning bridge due to its stretching as the instability develops. The kinetics of the bridge thinning and the size of satellite droplets formed after the bridge breakage depend considerably on the surfactant concentration and the value of critical micelle concentration. It is proposed that the dynamic surface tension on the time scale of the drop formation can be used as an effective surface tension for the description of the bridge kinetics over the broad range of experimental conditions used.