The influence of Marangoni convection on fluid dynamics of oscillating single rising droplets

Numerical Simulation on Insoluble Surfactant Mass Transfer on Deformable Bubble Interface in a Couette Flow by Phase-Field Lattice Boltzmann Method-Finite-Difference Method Hybrid Approach

Elaborate management on bubble shape and transportations depends on the balance between multiple physical behaviors for two-phase flow with Marangoni stress and the interface mass transfer. In this paper, a new model combining PFLBM (phase-field lattice Boltzmann method) and FDM ( finite-difference method) coupling with the ghost-cell method was built. The PFLBM-FDM was validated for the high accuracy, less computational cost, and low mass loss compared to other methods. Based on the PFLBM-FDM, a surfactant-laden bubble deformed and transported in a laminar Couette flow was investigated. The deformation ratio and transportation velocity were explored with different density ratios, surface tensions, shear velocities, and diffusion coefficients. The numerical results showed that the equilibrium state of the bubble deformation was decided only by the dimensionless numbers when the Sh number was higher than 100. Moreover, the transportation velocity of the bubble can be controlled by the balance between the Marangoni stress and shear velocity. When the Sh is lower than 100, the Marangoni stress from the surfactant is not a long-range force, which only works at the early flow. Otherwise, the Marangoni stress will be a long-range force that provides a persistent force to accelerate the bubble by ∼10%. Increasing ReH will further intensify the effect. Based on all the data, a correlation of the bubble deformation including with the densities of two fluids was obtained and the error range is less 5%.

Single Contaminated Drops Falling through Stagnant Liquid at Low Reynolds Numbers

Numerical simulations of contaminated spherical drops falling through a stagnant liquid at low Reynolds numbers are carried out using the finite difference method. The numerical results are used to describe the behavior of the surfactant concentrations and to understand the surfactant effects on the fluid motions in detail. The predicted interfacial surfactant concentration, Γ, is almost zero for angles, θ, below a certain value (the stagnant-cap angle, θcap), whereas it steeply increases and reaches a large value for θ>θcap (the stagnant-cap region). The increase in the initial surfactant concentration, C0, in the drop enhances the adsorption from the drop to the interface, which results in the increase in Γ and the decrease in θcap. Peaks appear in the predicted Marangoni stresses around θcap, which causes similar peaks in the pressure distribution. The high-pressure spots prevent the fluid motion along the interface, which results in the formation of the stagnant-cap region and the attenuation of the tangential velocity in the continuous phase. The surfactant flux from the bulk to the interface decreases C in the vicinity of the interface for θ<θcap and the weak diffusion cannot compensate for the reduction in C by adsorption, which results in C at the interface smaller than C0. The pattern of the low C region is determined by the advection and does not smear out because of a small diffusive flux.