Numerical study of the steady core-annular flow in a focusing geometry in the presence of soluble surfactants

Abstract

The steady state core annular flow of two immiscible fluids in the presence of soluble surfactants is investigated numerically in an axisymmetric tubular or flow focusing geometry. The two immiscible fluids form an interface onto which surfactants can be adsorbed/desorbed from the annular fluid where they are dissolved. The finite element method is employed to discretize the governing equations using the spine method to create the computational mesh, follow the shape of the interface and the axisymmetric wall and capture the flow domain. An unknown pressure distribution is imposed along the interface that quantifies the pressure discontinuity, while allowing the bulk pressure to remain continuous in each phase. GMRES is employed to accommodate the nearly singular nature of the jacobian matrix. The bulk surfactant concentration varies only inside a thin boundary layer of thickness δ∼Pe−1, which justifies omission of the advection-diffusion equation in the annulus bulk. Simulations reveal that the onset of Marangoni stresses tends to decelerate the velocity at the interface which is seen to translate towards the wall/axis of symmetry when the more viscous fluid is the inner/outer fluid, respectively. Finally, simulations in flow-focusing nozzles capture the development of spatially growing instabilities in the downstream region of the interface, with wavelength on the order of the nozzle radius, especially when the outer fluid is more viscous or when the surfactant concentration increases. Raising the volumetric flow rate in the annular region suppresses growth of waves due to interfacial friction. Spatial stability analysis verifies downstream growth of waves and recovers the relevant wavelength. Possible repercussions on the onset of jetting instabilities in devices with liquid focusing are discussed.