No Arabic abstract
The behaviour of microscopic swimmers has previously been explored near large scale confining geometries and in the presence of very small-scale surface roughness. Here we consider an intermediate case of how a simple microswimmer, the tangential spherical squirmer, behaves adjacent to singly- and doubly-periodic sinusoidal surface topographies that spatially oscillate with an amplitude that is an order of magnitude less than the swimmer size and wavelengths within an order of magnitude of this scale. The nearest neighbour regularised Stokeslet method is used for numerical explorations after validating its accuracy for a spherical tangential squirmer that swims stably near a flat surface. The same squirmer is then introduced to the different surface topographies. The key governing factor in the resulting swimming behaviour is the size of the squirmer relative to the surface topography wavelength. When the squirmer is much larger, or much smaller, than the surface topography wavelength then directional guidance is not observed. Once the squirmer size is on the scale of the topography wavelength limited guidance is possible, often with local capture in the topography troughs. However, complex dynamics can also emerge, especially when the initial configuration is not close to alignment along topography troughs or above topography crests. In contrast to sensitivity in alignment and the topography wavelength, reductions in the amplitude of the surface topography or variations in the shape of the periodic surface topography do not have extensive impacts on the squirmer behaviour. Our findings more generally highlight that detailed numerical explorations are necessary to determine how surface topographies may interact with a given swimmer that swims stably near a flat surface and whether surface topographies may be designed to provide such swimmers with directional guidance.
The emerging field of self-driven active particles in fluid environments has recently created significant interest in the biophysics and bioengineering communities owing to their promising future biomedical and technological applications. These microswimmers move autonomously through aqueous media where under realistic situations they encounter a plethora of external stimuli and confining surfaces with peculiar elastic properties. Based on a far-field hydrodynamic model, we present an analytical theory to describe the physical interaction and hydrodynamic couplings between a self-propelled active microswimmer and an elastic interface that features resistance toward shear and bending. We model the active agent as a superposition of higher-order Stokes singularities and elucidate the associated translational and rotational velocities induced by the nearby elastic boundary. Our results show that the velocities can be decomposed in shear and bending related contributions which approach the velocities of active agents close to a no-slip rigid wall in the steady limit. The transient dynamics predict that contributions to the velocities of the microswimmer due to bending resistance are generally more pronounced than to shear resistance. Our results provide insight into the control and guidance of artificial and synthetic self-propelling active microswimmers near elastic confinements.
The dynamics of an adhesive two-dimensional vesicle doublet under various flow conditions is investigated numerically using a high-order, adaptive-in-time boundary integral method. In a quiescent flow, two nearby vesicles move slowly towards each other under the adhesive potential, pushing out fluid between them to form a vesicle doublet at equilibrium. A lubrication analysis on such draining of a thin film gives the dependencies of draining time on adhesion strength and separation distance that are in good agreement with numerical results. In a planar extensional flow we find a stable vesicle doublet forms only when two vesicles collide head-on around the stagnation point. In a microfluid trap where the stagnation of an extensional flow is dynamically placed in the middle of a vesicle doublet through an active control loop, novel dynamics of a vesicle doublet are observed. Numerical simulations show that there exists a critical extensional flow rate above which adhesive interaction is overcome by the diverging stream, thus providing a simple method to measure the adhesion strength between two vesicle membranes. In a planar shear flow, numerical simulations reveal that a vesicle doublet may form provided that the adhesion strength is sufficiently large at a given vesicle reduced area. Once a doublet is formed, its oscillatory dynamics is found to depend on the adhesion strength and their reduced area. Furthermore the effective shear viscosity of a dilute suspension of vesicle doublets is found to be a function of the reduced area. Results from these numerical studies and analysis shed light on the hydrodynamic and rheological consequences of adhesive interactions between vesicles in a viscous fluid.
We investigate regular configurations of a small number of particles settling under gravity in a viscous fluid. The particles do not touch each other and can move relative to each other. The dynamics is analyzed in the point-particle approximation. A family of regular configurations is found with periodic oscillations of all the settling particles. The oscillations are shown to be robust under some out-of-phase rearrangements of the particles. In the presence of an additional particle above such a regular configuration, the particle periodic trajectories are horizontally repelled from the symmetry axis, and flattened vertically. The results are used to propose a mechanism how a spherical cloud, made of a large number of particles distributed at random, evolves and destabilizes.
A numerical study is presented to analyze the thermal mechanisms of unsteady, supersonic granular flow, by means of hydrodynamic simulations of the Navier-Stokes granular equations. For this purpose a paradigmatic problem in granular dynamics such as the Faraday instability is selected. Two different approaches for the Navier-Stokes transport coefficients for granular materials are considered, namely the traditional Jenkins-Richman theory for moderately dense quasi-elastic grains, and the improved Garzo-Dufty-Lutsko theory for arbitrary inelasticity, which we also present here. Both solutions are compared with event-driven simulations of the same system under the same conditions, by analyzing the density, the temperature and the velocity field. Important differences are found between the two approaches leading to interesting implications. In particular, the heat transfer mechanism coupled to the density gradient which is a distinctive feature of inelastic granular gases, is responsible for a major discrepancy in the temperature field and hence in the diffusion mechanisms.
The present article experimentally and theoretically probes the evaporation kinetics of sessile saline droplets. Observations reveal that presence of solvated ions leads to modulated evaporation kinetics, which is further a function of surface wettability. On hydrophilic surfaces, increasing salt concentration leads to enhanced evaporation rates, whereas on superhydrophobic surfaces, it first enhances and reduces with concentration. Also, the nature and extents of the evaporation regimes constant contact angle or constant contact radius are dependent on the salt concentration. The reduced evaporation on superhydrophobic surfaces has been explained based on observed via microscopy crystal nucleation behaviour within the droplet. Purely diffusion driven evaporation models are noted to be unable to predict the modulated evaporation rates. Further, the changes in the surface tension and static contact angles due to solvated salts also cannot explain the improved evaporation behaviour. Internal advection is observed using PIV to be generated within the droplet and is dependent on the salt concentration. The advection dynamics has been used to explain and quantify the improved evaporation behaviour by appealing to the concept of interfacial shear modified Stefan flows around the evaporating droplet. The analysis leads to accurate predictions of the evaporation rates. Further, another scaling analysis has been proposed to show that the thermal and solutal Marangoni advection within the system leads to the advection behaviour. The analysis also shows that the dominant mode is the solutal advection and the theory predicts the internal circulation velocities with good accuracy. The findings may be of importance to microfluidic thermal and species transport systems.