Hydrodynamic interaction of microswimmers near a wall

Do hydrodynamically assisted binary collisions lead to orientational ordering of microswimmers?

We have investigated the onset of collective motion in systems of model microswimmers, by performing a comprehensive analysis of the binary collision dynamics using three-dimensional direct numerical simulations (DNS) with hydrodynamic interactions. From this data, we have constructed a simplified binary collision model (BCM) which accurately reproduces the collective behavior obtained from the DNS for most cases. Thus, we show that global alignment can mostly arise solely from binary collisions. Although the agreement between both models (DNS and BCM) is not perfect, the parameter range in which notable differences appear is also that for which strong density fluctuations are present in the system (where pseudo-sound mound can be observed (N. Oyama et al., Phys. Rev. E 93, 043114 (2016))).

Guidance of microswimmers by wall and flow: Thigmotaxis and rheotaxis of unsteady squirmers in two and three dimensions

The motions of an unsteady circular-disk squirmer and a spherical squirmer have been investigated in the presence of a no-slip infinite wall and a background shear flow in order to clarify the similarities and differences between two- and three-dimensional motions. Despite the similar bifurcation structure of the dynamical system, the stability of the fixed points differs due to the Hamiltonian structure of the disk squirmer. Once the unsteady oscillating surface velocity profile is considered, the disk squirmer can behave in a chaotic manner and cease to be confined in a near-wall region. In contrast, in an unsteady spherical squirmer, the dynamics is well attracted by a stable fixed point. Additional wall contact interactions lead to stable fixed points for the disk squirmer, and, in turn, the surface entrapment of the disk squirmer can be stabilized, regardless of the existence of the background flow. Finally, we consider spherical motion under a background flow. The separated time scales of the surface entrapment (thigmotaxis) and the turning toward the flow direction (rheotaxis) enable us to reduce the dynamics to two-dimensional phase space, and simple weather-vane mechanics can predict squirmer rheotaxis. The analogous structure of the phase plane with the wall contact in two and three dimensions implies that the two-dimensional disk swimmer successfully captures the nonlinear interactions, and thus two-dimensional approximation could be useful in designing microfluidic devices for the guidance of microswimmers and for clarifying the locomotions in a complex geometry.

Modeling of active swimmer suspensions and their interactions with the environment

In this article, we review mathematical models used to study the behaviour of suspensions of micro-swimmers and the accompanying biophysical phenomena, with specific focus on stimulus response. The methods discussed encompass a range of interactions exhibited by the micro-swimmers; including passive hydrodynamic (gyrotaxis) and gravitational (gravitaxis) effects, and active responses to chemical cues (chemotaxis) and light intensities (phototaxis). We introduce the simplest models first, and then build towards more sophisticated recent developments, in the process, identifying the limitations of the former and the new results obtained by the latter. We comment on the accuracy/validity of the models adopted, based on the agreement between theoretical results and experimental observations. We conclude by identifying some of the open problems and associated challenges faced by researchers in the realm of active suspensions.

Hydrodynamics Defines the Stable Swimming Direction of Spherical Squirmers in a Nematic Liquid Crystal

We present a study of the hydrodynamics of an active particle-a model squirmer-in an environment with a broken rotational symmetry: a nematic liquid crystal. By combining simulations with analytic calculations, we show that the hydrodynamic coupling between the squirmer flow field and liquid crystalline director can lead to reorientation of the swimmers. The preferred orientation depends on the exact details of the squirmer flow field. In a steady state, pushers are shown to swim parallel with the nematic director while pullers swim perpendicular to the nematic director. This behavior arises solely from hydrodynamic coupling between the squirmer flow field and anisotropic viscosities of the host fluid. Our results suggest that an anisotropic swimming medium can be used to characterize and guide spherical microswimmers in the bulk.

Scattering of biflagellate microswimmers from surfaces

We use a three-bead-spring model to investigate the dynamics of biflagellate microswimmers near a surface. While the primary dynamics and scattering are governed by geometric-dependent direct contact, the fluid flows generated by the swimmer locomotion are important in orienting it toward or away from the surface. Flagellar noise and in particular cell spinning about the main axis help a surface-trapped swimmer escape, whereas the time a swimmer spends at the surface depends on the incident angle. The dynamics results from a nuanced interplay of direct collisions, hydrodynamics, noise, and the swimmer geometry. We show that to correctly capture the dynamics of a biflagellate swimmer, minimal models need to resolve the shape asymmetry.

Hydrodynamic interactions in dense active suspensions: From polar order to dynamical clusters

We study the role of hydrodynamic interactions in the collective behavior of collections of microscopic active particles suspended in a fluid. We introduce a calculational framework that allows us to separate the different contributions to their collective dynamics from hydrodynamic interactions on different length scales. Hence we are able to systematically show that lubrication forces when the particles are very close to each other play as important a role as long-range hydrodynamic interactions in determining their many-body behavior. We find that motility-induced phase separation is suppressed by near-field interactions, leading to open gel-like clusters rather than dense clusters. Interestingly, we find a globally polar ordered phase appears for neutral swimmers with no force dipole that is enhanced by near-field lubrication forces in which the collision process rather than long-range interaction dominates the alignment mechanism.

Population splitting of rodlike swimmers in Couette flow

We present a quantitative analysis on the response of a dilute active suspension of self-propelled rods (swimmers) in a planar channel subjected to an imposed shear flow. To best capture the salient features of the shear-induced effects, we consider the case of an imposed Couette flow, providing a constant shear rate across the channel. We argue that the steady-state behavior of swimmers can be understood in the light of a population splitting phenomenon, occurring as the shear rate exceeds a certain threshold, initiating the reversal of the swimming direction for a finite fraction of swimmers from down- to upstream or vice versa, depending on the swimmer position within the channel. Swimmers thus split into two distinct, statistically significant and oppositely swimming majority and minority populations. The onset of population splitting translates into a transition from a self-propulsion-dominated regime to a shear-dominated regime, corresponding to a unimodal-to-bimodal change in the probability distribution function of the swimmer orientation. We present a phase diagram in terms of the swim and flow Péclet numbers showing the separation of these two regimes by a discontinuous transition line. Our results shed further light on the behavior of swimmers in a shear flow and provide an explanation for the previously reported non-monotonic behavior of the mean, near-wall, parallel-to-flow orientation of swimmers with increasing shear strength.

Traction reveals mechanisms of wall effects for microswimmers near boundaries

The influence of a plane boundary on low-Reynolds-number swimmers has frequently been studied using image systems for flow singularities. However, the boundary effect can also be expressed using a boundary integral representation over the traction on the boundary. We show that examining the traction pattern on the boundary caused by a swimmer can yield physical insights into determining when far-field multipole models are accurate. We investigate the swimming velocities and the traction of a three-sphere swimmer initially placed parallel to an infinite planar wall. In the far field, the instantaneous effect of the wall on the swimmer is well approximated by that of a multipole expansion consisting of a force dipole and a force quadrupole. On the other hand, the swimmer close to the wall must be described by a system of singularities reflecting its internal structure. We show that these limits and the transition between them can be independently identified by examining the traction pattern on the wall, either using a quantitative correlation coefficient or by visual inspection. Last, we find that for nonconstant propulsion, correlations between swimming stroke motions and internal positions are important and not captured by time-averaged traction on the wall, indicating that care must be taken when applying multipole expansions to study boundary effects in cases of nonconstant propulsion.

Lattice-Boltzmann simulations of microswimmer-tracer interactions

Hydrodynamic interactions in systems composed of self-propelled particles, such as swimming microorganisms and passive tracers, have a significant impact on the tracer dynamics compared to the equivalent "dry" sample. However, such interactions are often difficult to take into account in simulations due to their computational cost. Here, we perform a systematic investigation of swimmer-tracer interaction using an efficient force-counterforce-based lattice-Boltzmann (LB) algorithm [De Graaf et al., J. Chem. Phys. 144, 134106 (2016)JCPSA60021-960610.1063/1.4944962] in order to validate its ability to capture the relevant low-Reynolds-number physics. We show that the LB algorithm reproduces far-field theoretical results well, both in a system with periodic boundary conditions and in a spherical cavity with no-slip walls, for which we derive expressions here. The force-lattice coupling of the LB algorithm leads to a "smearing out" of the flow field, which strongly perturbs the tracer trajectories at close swimmer-tracer separations, and we analyze how this effect can be accurately captured using a simple renormalized hydrodynamic theory. Finally, we show that care must be taken when using LB algorithms to simulate systems of self-propelled particles, since its finite momentum transport time can lead to significant deviations from theoretical predictions based on Stokes flow. These insights should prove relevant to the future study of large-scale microswimmer suspensions using these methods.

Dynamics of a microorganism in a sheared viscoelastic liquid

In this paper, we investigate the dynamics of a model spherical microorganism, called squirmer, suspended in a viscoelastic fluid undergoing unconfined shear flow. The effect of the interplay of shear flow, fluid viscoelasticity, and self-propulsion on the orientational dynamics is addressed. In the limit of weak viscoelasticity, quantified by the Deborah number, an analytical expression for the squirmer angular velocity is derived by means of the generalized reciprocity theorem. Direct finite element simulations are carried out to study the squirmer dynamics at larger Deborah numbers. Our results show that the orientational dynamics of active microorganisms in a sheared viscoelastic fluid greatly differs from that observed in Newtonian suspensions. Fluid viscoelasticity leads to a drift of the particle orientation vector towards the vorticity axis or the flow-gradient plane depending on the Deborah number, the relative weight between the self-propulsion velocity and the flow characteristic velocity, and the type of swimming. Generally, pullers and pushers show an opposite equilibrium orientation. The results reported in the present paper could be helpful in designing devices where separation of microorganisms, based on their self-propulsion mechanism, is obtained.

Phase separation and coexistence of hydrodynamically interacting microswimmers

A striking feature of the collective behavior of spherical microswimmers is that for sufficiently strong self-propulsion they phase-separate into a dense cluster coexisting with a low-density disordered surrounding. Extending our previous work, we use the squirmer as a model swimmer and the particle-based simulation method of multi-particle collision dynamics to explore the influence of hydrodynamics on their phase behavior in a quasi-two-dimensional geometry. The coarsening dynamics towards the phase-separated state is diffusive in an intermediate time regime followed by a final ballistic compactification of the dense cluster. We determine the binodal lines in a phase diagram of Péclet number versus density. Interestingly, the gas binodals are shifted to smaller densities for increasing mean density or dense-cluster size, which we explain using a recently introduced pressure balance [S. C. Takatori, et al., Phys. Rev. Lett. 2014, 113, 028103] extended by a hydrodynamic contribution. Furthermore, we find that for pushers and pullers the binodal line is shifted to larger Péclet numbers compared to neutral squirmers. Finally, when lowering the Péclet number, the dense phase transforms from a hexagonal "solid" to a disordered "fluid" state.

Hydrodynamic interaction of swimming organisms in an inertial regime

We numerically investigate the hydrodynamic interaction of swimming organisms at small to intermediate Reynolds number regimes, i.e., Re∼O(0.1-100), where inertial effects are important. The hydrodynamic interaction of swimming organisms in this regime is significantly different from the Stokes regime for microorganisms, as well as the high Reynolds number flows for fish and birds, which involves strong flow separation and detached vortex structures. Using an archetypal swimmer model, called a "squirmer," we find that the inertial effects change the contact time and dispersion dynamics of a pair of pusher swimmers, and trigger hydrodynamic attraction for two pullers. These results are potentially important in investigating predator-prey interactions, sexual reproduction, and the encounter rate of marine organisms such as copepods, ctenophora, and larvae.

Hydrodynamic oscillations and variable swimming speed in squirmers close to repulsive walls

We present a lattice Boltzmann study of the hydrodynamics of a fully resolved squirmer, confined in a slab of fluid between two no-slip walls. We show that the coupling between hydrodynamics and short-range repulsive interactions between the swimmer and the surface can lead to hydrodynamic trapping of both pushers and pullers at the wall, and to hydrodynamic oscillations in the case of a pusher. We further show that a pusher moves significantly faster when close to a surface than in the bulk, whereas a puller undergoes a transition between fast motion and a dynamical standstill according to the range of the repulsive interaction. Our results critically require near-field hydrodynamics and demonstrate that far-field hydrodynamics is insufficient to give even a qualitatively correct account of swimmer behaviour near walls. Finally our simulations suggest that it should be possible to control the density and speed of squirmers at a surface by tuning the range of steric and electrostatic swimmer-wall interactions.

Modeling a spheroidal microswimmer and cooperative swimming in a narrow slit

We propose a hydrodynamic model for a spheroidal microswimmer with two tangential surface velocity modes. This model is analytically solvable and reduces to Lighthill's and Blake's spherical squirmer model in the limit of equal major and minor semi-axes. Furthermore, we present an implementation of such a spheroidal squirmer by means of particle-based mesoscale hydrodynamics simulations using the multiparticle collision dynamics approach. We investigate its properties as well as the scattering of two spheroidal squirmers in a slit geometry. Thereby we find a stable fixed point, where two pullers swim cooperatively forming a wedge-like conformation with a small constant angle.

Taylor line swimming in microchannels and cubic lattices of obstacles

Microorganisms naturally move in microstructured fluids. Using the simulation method of multi-particle collision dynamics, we study in two dimensions an undulatory Taylor line swimming in a microchannel and in a cubic lattice of obstacles, which represent simple forms of a microstructured environment. In the microchannel the Taylor line swims at an acute angle along a channel wall with a clearly enhanced swimming speed due to hydrodynamic interactions with the bounding wall. While in a dilute obstacle lattice swimming speed is also enhanced, a dense obstacle lattice gives rise to geometric swimming. This new type of swimming is characterized by a drastically increased swimming speed. Since the Taylor line has to fit into the free space of the obstacle lattice, the swimming speed is close to the phase velocity of the bending wave traveling along the Taylor line. While adjusting its swimming motion within the lattice, the Taylor line chooses a specific swimming direction, which we classify by a lattice vector. When plotting the swimming velocity versus the magnitude of the lattice vector, all our data collapse on a single master curve. Finally, we also report more complex trajectories within the obstacle lattice.

Dimensionality matters in the collective behaviour of active emulsions

The behaviour of artificial microswimmers consisting of droplets of a mesogenic oil immersed in an aqueous surfactant solution depends qualitatively on the conditions of dimensional confinement; ranging from only transient aggregates in Hele-Shaw geometries to hexagonally packed, convection-driven clusters when sedimenting in an unconfined reservoir. We study the effects of varying the swimmer velocity, the height of the reservoir, and the buoyancy of the droplet swimmers. Two simple adjustments of the experimental setting lead to a suppression of clustering: either a decrease of the reservoir height below a certain value, or a match of the densities of droplets and surrounding phase, showing that the convection is the key mechanism for the clustering behaviour.

Understanding the onset of oscillatory swimming in microchannels

Self-propelled colloids (swimmers) in confining geometries follow trajectories determined by hydrodynamic interactions with the bounding surfaces. However, typically these interactions are ignored or truncated to the lowest order. We demonstrate that higher-order hydrodynamic moments cause rod-like swimmers to follow oscillatory trajectories in quiescent fluid between two parallel plates, using a combination of lattice-Boltzmann simulations and far-field calculations. This behavior occurs even far from the confining walls and does not require lubrication results. We show that a swimmer's hydrodynamic quadrupole moment is crucial to the onset of the oscillatory trajectories. This insight allows us to develop a simple model for the dynamics near the channel center based on these higher hydrodynamic moments, and suggests opportunities for trajectory-based experimental characterization of swimmers' hydrodynamic properties.

Purely hydrodynamic origin for swarming of swimming particles

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the collective motion of model swimmers in bulk and confinement. Calculating the dynamic structure factor, we clarified that the swarming in bulk systems can be understood as a pseudoacoustic mode. Under confinement between flat parallel walls, this pseudoacoustic mode leads to a traveling wavelike motion. This swarming behavior is due purely to the hydrodynamic interactions between the swimmers and depends strongly on the type and strength of swimming (i.e., pusher or puller).

Getting drowned in a swirl: Deformable bead-spring model microswimmers in external flow fields

Deformability is a central feature of many types of microswimmers, e.g., for artificially generated self-propelled droplets. Here, we analyze deformable bead-spring microswimmers in an externally imposed solvent flow field as simple theoretical model systems. We focus on their behavior in a circular swirl flow in two spatial dimensions. Linear (straight) two-bead swimmers are found to circle around the swirl with a slight drift to the outside with increasing activity. In contrast to that, we observe for triangular three-bead or squarelike four-bead swimmers a tendency of being drawn into the swirl and finally getting drowned, although a radial inward component is absent in the flow field. During one cycle around the swirl, the self-propulsion direction of an active triangular or squarelike swimmer remains almost constant, while their orbits become deformed exhibiting an "egglike" shape. Over time, the swirl flow induces slight net rotations of these swimmer types, which leads to net rotations of the egg-shaped orbits. Interestingly, in certain cases, the orbital rotation changes sense when the swimmer approaches the flow singularity. Our predictions can be verified in real-space experiments on artificial microswimmers.

Swimming in a crystal

We study catalytic Janus particles and Escherichia coli bacteria swimming in a two-dimensional colloidal crystal. The Janus particles orbit individual colloids and hop between colloids stochastically, with a hopping rate that varies inversely with fuel (hydrogen peroxide) concentration. At high fuel concentration, these orbits are stable for 100s of revolutions, and the orbital speed oscillates periodically as a result of hydrodynamic, and possibly also phoretic, interactions between the swimmer and the six neighbouring colloids. Motile E. coli bacteria behave very differently in the same colloidal crystal: their circular orbits on plain glass are rectified into long, straight runs, because the bacteria are unable to turn corners inside the crystal.

Microalgae Scatter off Solid Surfaces by Hydrodynamic and Contact Forces

Interactions between microorganisms and solid boundaries play an important role in biological processes, such as egg fertilization, biofilm formation, and soil colonization, where microswimmers move within a structured environment. Despite recent efforts to understand their origin, it is not clear whether these interactions can be understood as being fundamentally of hydrodynamic origin or hinging on the swimmer's direct contact with the obstacle. Using a combination of experiments and simulations, here we study in detail the interaction of the biflagellate green alga Chlamydomonas reinhardtii, widely used as a model puller microorganism, with convex obstacles, a geometry ideally suited to highlight the different roles of steric and hydrodynamic effects. Our results reveal that both kinds of forces are crucial for the correct description of the interaction of this class of flagellated microorganisms with boundaries.

Biogenic mixing induced by intermediate Reynolds number swimming in stratified fluids

We study fully resolved motion of interacting swimmers in density stratified fluids using an archetypal swimming model called “squirmer”. The intermediate Reynolds number regime is particularly important, because the vast majority of organisms in the aphotic ocean (i.e. regions that are 200 m beneath the sea surface) are small (mm-cm) and their motion is governed by the balance of inertial and viscous forces. Our study shows that the mixing efficiency and the diapycnal eddy diffusivity, a measure of vertical mass flux, within a suspension of squirmers increases with Reynolds number. The mixing efficiency is in the range of O(0.0001–0.04) when the swimming Reynolds number is in the range of O(0.1–100). The values of diapycnal eddy diffusivity and Cox number are two orders of magnitude larger for vertically swimming cells compared to horizontally swimming cells. For a suspension of squirmers in a decaying isotropic turbulence, we find that the diapycnal eddy diffusivity enhances due to the strong viscous dissipation generated by squirmers as well as the interaction of squirmers with the background turbulence.

Flagellar swimmers oscillate between pusher- and puller-type swimming

Self-propulsion of cellular microswimmers generates flow signatures, commonly classified as pusher and puller type, which characterize hydrodynamic interactions with other cells or boundaries. Using experimentally measured beat patterns, we compute that the flagellated green alga Chlamydomonas oscillates between pusher and puller, rendering it an approximately neutral swimmer, when averaging over its full beat cycle. Beyond a typical distance of 100μm from the cell, inertia attenuates oscillatory microflows. We show that hydrodynamic interactions between cells oscillate in time and are of similar magnitude as stochastic swimming fluctuations. From our analysis, we also find that the rate of hydrodynamic dissipation varies in time, which implies that flagellar beat patterns are not optimized with respect to this measure.

Locomotion of a microorganism in weakly viscoelastic liquids

In the present work we study the motion of microorganisms swimming by an axisymmetric distribution of surface tangential velocity in a weakly viscoelastic fluid. The second-order fluid constitutive equation is used to model the suspending fluid, while the well-known "squirmer model" [M. J. Lighthill, Comm. Pure Appl. Math. 5, 109 (1952); J. R. Blake, J. Fluid Mech. 46, 199 (1971)] is employed to describe the organism propulsion mechanism. A regular perturbation expansion up to first order in the Deborah number is performed, and the generalized reciprocity theorem from Stokes flow theory is then used, to derive analytical formulas for the squirmer velocity. Results show that "neutral" squirmers are unaffected by viscoelasticity, whereas "pullers" and "pushers" are slowed down and hastened, respectively. The power dissipated by the swimming microorganism and the "swimming efficiency" are also analytically quantified.

Rheotaxis of spherical active particles near a planar wall

For active particles the interplay between the self-generated hydrodynamic flow and an external shear flow, especially near bounding surfaces, can result in a rich behavior of the particles not easily foreseen from the consideration of the active and external driving mechanisms in isolation. For instance, under certain conditions, the particles exhibit "rheotaxis", i.e., they align their direction of motion with the plane of shear spanned by the direction of the flow and the normal of the bounding surface and move with or against the flow. To date, studies of rheotaxis have focused on elongated particles (e.g., spermatozoa), for which rheotaxis can be understood intuitively in terms of a "weather vane" mechanism. Here we investigate the possibility that spherical active particles, for which the "weather vane" mechanism is excluded due to the symmetry of the shape, may nevertheless exhibit rheotaxis. Combining analytical and numerical calculations, we show that, for a broad class of spherical active particles, rheotactic behavior may emerge via a mechanism which involves "self-trapping" near a hard wall owing to the active propulsion of the particles, combined with their rotation, alignment, and "locking" of the direction of motion into the shear plane. In this state, the particles move solely up- or downstream at a steady height and orientation.

Detention Times of Microswimmers Close to Surfaces: Influence of Hydrodynamic Interactions and Noise

After colliding with a surface, microswimmers reside there during the detention time. They accumulate and may form complex structures such as biofilms. We introduce a general framework to calculate the distribution of detention times using the method of first-passage times and study how rotational noise and hydrodynamic interactions influence the escape from a surface. We compare generic swimmer models to the simple active Brownian particle. While the respective detention times of source dipoles are smaller, the ones of pullers are larger by up to several orders of magnitude, and pushers show both trends. We apply our results to the more realistic squirmer model, for which we use lubrication theory, and validate them by simulations with multiparticle collision dynamics.

Interplay of physical mechanisms and biofilm processes: review of microfluidic methods

Bacteria in natural and artificial environments often reside in self-organized, integrated communities known as biofilms. Biofilms are highly structured entities consisting of bacterial cells embedded in a matrix of self-produced extracellular polymeric substances (EPS). The EPS matrix acts like a biological 'glue' enabling microbes to adhere to and colonize a wide range of surfaces. Once integrated into biofilms, bacterial cells can withstand various forms of stress such as antibiotics, hydrodynamic shear and other environmental challenges. Because of this, biofilms of pathogenic bacteria can be a significant health hazard often leading to recurrent infections. Biofilms can also lead to clogging and material degradation; on the other hand they are an integral part of various environmental processes such as carbon sequestration and nitrogen cycles. There are several determinants of biofilm morphology and dynamics, including the genotypic and phenotypic states of constituent cells and various environmental conditions. Here, we present an overview of the role of relevant physical processes in biofilm formation, including propulsion mechanisms, hydrodynamic effects, and transport of quorum sensing signals. We also provide a survey of microfluidic techniques utilized to unravel the associated physical mechanisms. Further, we discuss the future research areas for exploring new ways to extend the scope of the microfluidic approach in biofilm studies.

Effect of solid boundaries on swimming dynamics of microorganisms in a viscoelastic fluid

We numerically study the effect of solid boundaries on the swimming behavior of a motile microorganism in viscoelastic media. Understanding the swimmer-wall hydrodynamic interactions is crucial to elucidate the adhesion of bacterial cells to nearby substrates which is precursor to the formation of the microbial biofilms. The microorganism is simulated using a squirmer model that captures the major swimming mechanisms of potential, extensile, and contractile types of swimmers, while neglecting the biological complexities. A Giesekus constitutive equation is utilized to describe both viscoelasticity and shear-thinning behavior of the background fluid. We found that the viscoelasticity strongly affects the near-wall motion of a squirmer by generating an opposing polymeric torque which impedes the rotation of the swimmer away from the wall. In particular, the time a neutral squirmer spends at the close proximity of the wall is shown to increase with polymer relaxation time and reaches a maximum at Weissenberg number of unity. The shear-thinning effect is found to weaken the solvent stress and therefore, increases the swimmer-wall contact time. For a puller swimmer, the polymer stretching mainly occurs around its lateral sides, leading to reduced elastic resistance against its locomotion. The neutral and puller swimmers eventually escape the wall attraction effect due to a releasing force generated by the Newtonian viscous stress. In contrast, the pusher is found to be perpetually trapped near the wall as a result of the formation of a highly stretched region behind its body. It is shown that the shear-thinning property of the fluid weakens the wall-trapping effect for the pusher squirmer.

Locomotion of microorganisms near a no-slip boundary in a viscoelastic fluid

Locomotion of microorganisms plays a vital role in most of their biological processes. In many of these processes, microorganisms are exposed to complex fluids while swimming in confined domains, such as spermatozoa in mucus of mammalian reproduction tracts or bacteria in extracellular polymeric matrices during biofilm formation. Thus, it is important to understand the kinematics of propulsion in a viscoelastic fluid near a no-slip boundary. We use a squirmer model with a time-reversible body motion to analytically investigate the swimming kinematics in an Oldroyd-B fluid near a wall. Analysis of the time-averaged motion of the swimmer shows that both pullers and pushers in a viscoelastic fluid swim towards the no-slip boundary if they are initially located within a small domain of "attraction" in the vicinity of the wall. In contrast, neutral swimmers always move towards the wall regardless of their initial distance from the wall. Outside the domain of attraction, pullers and pushers are both repelled from the no-slip boundary. Time-averaged locomotion is most pronounced at a Deborah number of unity. We examine the swimming trajectories of different types of swimmers as a function of their initial orientation and distance from the no-slip boundary.