Search

A point force acting on a Brinkman fluid in confinement is always counterbalanced by the force on the porous medium, the force on the walls and the stress at open boundaries. We discuss the distribution of those forces in different geometries: a long pipe, a medium with a single no-slip planar boundary, a porous sphere with an open boundary and a porous sphere with a no-slip wall. We determine the forces using the Lorentz reciprocal theorem and additionally validate the results with explicit analytical flow solutions. We discuss the relevance of our findings for cellular processes such as cytoplasmic streaming and centrosome positioning., Comment: 16 pages, 5 figures

We present a two-phase field model and a hybrid particle-phase field model to simulate dilute colloidal sedimentation and flotation near a liquid-gas interface (or fluid-fluid interface in general). Both models are coupled to the incompressible Stokes equation, which is solved numerically using a combination of sine and regular Fourier transforms to account for the no-slip boundary conditions at the boundaries. The continuum two-phase field model allows us to analytically solve the equilibrium interfacial profile using a perturbative approach, demonstrating excellent agreement with numerical simulations. Notably, we show that strong coupling to particle dynamics can significantly alter the liquid-gas interface, thereby modifying the liquid-gas interfacial tension. In particular, we show that the renormalized surface tension is monotonically decreasing with increasing colloidal particle concentration and decreasing buoyant mass.

Dynamics at low Reynolds numbers experiences recent revival in the fields of biophysics and active matter. While in bulk isotropic fluids it is exhaustively studied, this is less so in anisotropic fluids and in confined situations. Here, we combine the latter two by studying the rotation of a disk-like inclusion in a uniaxially anisotropic, globally oriented, incompressible two-dimensional fluid film. In terms of a perturbative expansion in parameters that quantify anisotropies in viscosity and in additional linear friction with a supporting substrate or other type of confinement, we derive analytical expressions for the resulting hydrodynamic flow and pressure fields as well as for the resistance and mobility coefficients of the rotating disk. It turns out that, in contrast to translational motion, the solutions remain well-behaved also in the absence of the additional linear friction. Comparison with results from finite-element simulations show very good agreement with those from our analytical calculations. Besides applications to describe technological systems, for instance, in the area of microfluidics and thin cells of aligned nematic liquid crystals, our solutions are important for quantitative theoretical approaches to fluid membranes and thin films in general featuring a preferred direction., Comment: 17 pages, 5 figures

Functionalized thin elastic films and membranes frequently feature internal sites of net forces or stresses. These are, for instance, active sites of actuation, or rigid inclusions in a strained membrane that induce counterstress upon externally imposed deformations. We theoretically analyze the geometry of isotropic, flat, thin, linearly elastic films or membranes of finite thickness, laterally extended to infinity. At the mathematical core of such characterizations are the fundamental solutions for localized force and stress singularities associated with corresponding Green's functions. We derive such solutions in three dimensions and place them into the context of previous two-dimensional calculations. To this end, we consider both no-slip and stress-free conditions at the top and/or bottom surfaces. We provide an understanding for why the fully free-standing thin elastic membrane leads to diverging solutions in most geometries and compare these situations to the truly two-dimensional case. A no-slip support of at least one of the surfaces stabilizes the solution, which illustrates that the divergences in the fully free-standing case are connected to global deformations. Within the mentioned framework, our results are important for associated theoretical characterizations of thin elastic films, whether supported or free-standing, and of membranes subject to internal or external forces or stresses., Comment: 13 pages

To make progress towards the development of a theory on the motion of inclusions in thin structured films and membranes, we here consider as an initial step a circular disk in a two-dimensional, uniaxially anisotropic fluid layer. We assume overdamped dynamics, incompressibility of the fluid, and global alignment of the axis of anisotropy. Motion within this layer is affected by additional linear friction with the environment, for instance, a supporting substrate. We investigate the induced flows in the fluid when the disk is translated parallel or perpendicular to the direction of anisotropy. Moreover, expressions for corresponding mobilities and resistance coefficients of the disk are derived. Our results are obtained within the framework of a perturbative expansion in the parameters that quantify the anisotropy of the fluid. Good agreement is found for moderate anisotropy when compared to associated results from finite-element simulations. At pronounced anisotropy, the induced flow fields are still predicted qualitatively correctly by the perturbative theory, although quantitative deviations arise. We hope to stimulate with our investigations corresponding experimental analyses, for example, concerning fluid flows in anisotropic thin films on uniaxially rubbed supporting substrates., Comment: 14 pages, 5 figures

Microorganisms and synthetic microswimmers often encounter complex environments consisting of networks of obstacles embedded into viscous fluids. Such settings include biological media, such as mucus with filamentous networks, as well as environmental scenarios, including wet soil and aquifers. A fundamental question in studying their locomotion is how the impermeability of these porous media impact their propulsion performance compared with the case that in a purely viscous fluid. Previous studies showed that the additional resistance due to the embedded obstacles leads to an enhanced propulsion of different types of swimmers, including undulatory swimmers, helical swimmers, and squirmers. In this work we employ a canonical three-sphere swimmer model to probe the impact of propulsion in porous media. The Brinkman equation is utilized to model a sparse network of stationary obstacles embedded into an incompressible Newtonian liquid. We present both a far-field theory and numerical simulations to characterize the propulsion performance of the swimmer in such porous media. In contrast to enhanced propulsion observed in other swimmer models, our results reveal that both the propulsion speed and efficiency of the three-sphere swimmer are largely reduced by the impermeability of the porous medium. We attribute the substantial reduction in propulsion performance to the screened hydrodynamic interactions among the spheres due to the more rapid spatial decays of flows in Brinkman media. These results highlight how enhanced or hindered propulsion in porous media is largely dependent on individual propulsion mechanisms. The specific example and physical insights provided here may guide the design of synthetic microswimmers for effective locomotion in porous media in their potential biological and environmental applications., Comment: 23 pages, 4 figures

A Marangoni surfer is an object embedded in a gas-liquid interface, propelled by gradients in surface tension. We derive an analytical theorem for the lower bound on the viscous dissipation by a Marangoni surfer in the limit of small Reynolds and capillary numbers. The minimum dissipation can be expressed with the reciprocal difference between drag coefficients of two passive bodies of the same shape as the Marangoni surfer, one in a force-free interface and the other in an interface with surface incompressibility. The distribution of surface tension that gives the optimal propulsion is given by the surface tension of the solution for the incompressible surface and the flow is a superposition of both solutions. For a surfer taking the form of a thin circular disk, the minimum dissipation is $16\mu a V^2$, giving a Lighthill efficiency of $1/3$. This places the Marangoni surfers among the hydrodynamically most efficient microswimmers., Comment: 17 pages, 5 figures

Swimming at the microscale has recently garnered substantial attention due to the fundamental biological significance of swimming microorganisms and the wide range of biomedical applications for artificial microswimmers. These microswimmers invariably find themselves surrounded by different confining boundaries, which can impact their locomotion in significant and diverse ways. In this work, we employ a widely used three-sphere swimmer model to investigate the effect of confinement on swimming at low Reynolds numbers. We conduct theoretical analysis via the point-particle approximation and numerical simulations based on the finite element method to examine the motion of the swimmer along the centerline in a capillary tube. The axisymmetric configuration reduces the motion to one-dimensional movement, which allows us to quantify how the degree of confinement affects the propulsion speed in a simple manner. Our results show that the confinement does not significantly affect the propulsion speed until the ratio of the radius of the tube to the radius of the sphere is in the range of $\mathcal{O}(1)-\mathcal{O}(10)$, where the swimmer undergoes substantial reduction in its propulsion speed as the radius of the tube decreases. We provide some physical insights into how reduced hydrodynamic interactions between moving spheres under confinement may hinder the propulsion of the three-sphere swimmer. We also remark that the reduced propulsion performance stands in stark contrast to the enhanced helical propulsion observed in a capillary tube, highlighting how the manifestation of confinement effects can vary qualitatively depending on the propulsion mechanisms employed by the swimmers., Comment: 10 pages, 7 figures, revised manuscript resubmitted to Phys Fluids

Green's function plays an important role in many areas of physical sciences and is a prime tool for solving diverse hydrodynamic equations in the linear regime. In the present contribution, the axisymmetric low-Reynolds-number Brinkman flow induced by monopole and dipole singularities in proximity of a stationary plate of circular shape is theoretically investigated. The flow singularities are directed along the central axis of the plate. No-slip boundary conditions are assumed to hold at the surface of the plate. The Green's functions are determined to a large extent analytically, reducing the solution of the linear hydrodynamic equations to well-behaved one-dimensional integrals amenable to numerical computation. In our approach, the Brinkman flow problem is formulated as a mixed boundary value problem that is subsequently mapped in the form of dual integral equations on the domain boundaries. Thereupon, the solution of the equations of fluid motion is eventually reduced to the solution of two independent Fredholm integral equations of the first kind. The overall flow structure and emerging eddy patterns are found to strongly depend on the magnitude of the relevant geometrical and physical parameters of the system. Moreover, the effect of the confining plate on the dynamics of externally driven or force-free particles is assessed through the calculation of the relevant hydrodynamic reaction functions. The effect of the plate on the locomotory behavior of a self-propelling active dipole swimmer is shown to be maximum when the radius of the plate is comparable to the distance separating the swimmer from the plate. Our results may prove useful for characterizing transport processes in microfluidic devices and may pave the way toward understanding and controlling of small-scale flows in porous media., Comment: 17 pages, 7 figures, resubmitted to Phys Rev Res

We theoretically and computationally study the low-Reynolds-number hydrodynamics of a linear active microswimmer surfing on a compressible thin fluid layer characterized by an odd viscosity. Since the underlying three-dimensional fluid is assumed to be very thin compared to any lateral size of the fluid layer, the model is effectively two-dimensional. In the limit of small odd viscosity compared to the even viscosities of the fluid layer, we obtain analytical expressions for the self-induced flow field, which includes non-reciprocal components due to the odd viscosity. On this basis, we fully analyze the behavior of a single linear swimmer, finding that it follows a circular path, the radius of which is, to leading order, inversely proportional to the magnitude of the odd viscosity. In addition, we show that a pair of swimmers exhibits a wealth of two-body dynamics that depends on the initial relative orientation angles as well as on the propulsion mechanism adopted by each swimmer. In particular, the pusher-pusher and pusher-puller-type swimmer pairs exhibit a generic spiral motion, while the puller-puller pair is found to either co-rotate in the steady state along a circular trajectory or exhibit a more complex chaotic behavior resulting from the interplay between hydrodynamic and steric interactions. Our theoretical predictions may pave the way toward a better understanding of active transport in active chiral fluids with odd viscosity, and may find potential applications in the quantitative microrheological characterization of odd-viscous fluids., Comment: 21 pages, 8 figures

The energy dissipation and entropy production by self-propelled microswimmers differ profoundly from passive particles pulled by external forces. The difference extends both to the shape of the flow around the swimmer, as well as to the internal dissipation of the propulsion mechanism. Here we derive a general theorem that provides an exact lower bound on the total, external and internal, dissipation by a microswimmer. The problems that can be solved include an active surface-propelled droplet, swimmers with an extended propulsive layer and swimmers with an effective internal dissipation. We apply the theorem to determine the swimmer shapes that minimize the total dissipation while keeping the volume constant. Our results show that the entropy production by active microswimmers is subject to different fundamental limits than the entropy production by externally driven particles., Comment: 14+5 pages, 6+1 figures

Although the dynamics of colloids in the vicinity of a solid interface has been widely characterized in the past, experimental studies of Brownian diffusion close to an air-water interface are rare and limited to particle-interface gap distances larger than the particle size. At the still unexplored lower distances, the dynamics is expected to be extremely sensitive to boundary conditions at the air-water interface. There, ad hoc experiments would provide a quantitative validation of predictions. Using a specially designed dual wave interferometric setup, the 3D dynamics of 9 micrometers diameter particles at a few hundreds of nanometers from an air-water interface is here measured in thermal equilibrium. Intriguingly, while the measured dynamics parallel to the interface approaches expected predictions for slip boundary conditions, the Brownian motion normal to the interface is very close to the predictions for no-slip boundary conditions. These puzzling results are rationalized considering current models of incompressible interfacial flow and deepened developing an ad hoc model which considers the contribution of tiny concentrations of surface active particles at the interface. We argue that such condition governs the particle dynamics in a large spectrum of systems ranging from biofilm formation to flotation process., Comment: 33 pages, 5 figures. Submitted to Journal of Colloid and Interface Science

Casimir forces in quantum electrodynamics emerge between microscopic metallic objects because of the confinement of the vacuum electromagnetic fluctuations occurring even at zero temperature. Their generalization at finite temperature and in material media are referred to as Casimir--Lifshitz forces. These forces are typically attractive, leading to the widespread problem of stiction between the metallic parts of micro- and nanodevices. Recently, repulsive Casimir forces have been experimentally realized but their reliance on specialized materials prevents their dynamic control and thus limits their further applicability. Here, we experimentally demonstrate that repulsive critical Casimir forces, which emerge in a critical binary liquid mixture upon approaching the critical temperature, can be used to actively control microscopic and nanoscopic objects with nanometer precision. We demonstrate this by using critical Casimir forces to prevent the stiction caused by the Casimir--Lifshitz forces. We study a microscopic gold flake above a flat gold-coated substrate immersed in a critical mixture. Far from the critical temperature, stiction occurs because of dominant Casimir--Lifshitz forces. Upon approaching the critical temperature, however, we observe the emergence of repulsive critical Casimir forces that are sufficiently strong to counteract stiction. This experimental demonstration can accelerate the development of micro- and nanodevices by preventing stiction as well as providing active control and precise tunability of the forces acting between their constituent parts., Comment: 27 pages, 5 figures

Breaking spatial symmetry is an essential requirement for phoretic active particles to swim at low Reynolds number. This fundamental prerequisite for swimming at the micro-scale is fulfilled either by chemical patterning of the surface of active particles or alternatively by exploiting geometrical asymmetries to induce chemical gradients and achieve self-propulsion. In the present manuscript, a far-field analytical model is employed to quantify the leading-order contribution to the induced phoretic velocity of a chemically homogeneous isotropic active colloid near a finite-sized disk of circular shape resting on an interface separating two immiscible viscous incompressible Newtonian fluids. To this aim, the solution of the phoretic problem is formulated as a mixed-boundary-value problem which is subsequently transformed into a system of dual integral equations on the inner and outer domains. Depending on the ratio of different involved viscosities and solute solubilities, the sign of phoretic mobility and chemical activity, as well as the ratio of particle-interface distance to the radius of the disk, the isotropic active particle is found to be either repelled from the interface, attracted to it, or reach a stable hovering state and remains immobile near the interface. Our results may prove useful in controlling and guiding the motion of self-propelled phoretic active particles near aqueous interfaces.

Geometric confinements play an important role in many physical and biological processes and significantly affect the rheology and behavior of colloidal suspensions at low Reynolds numbers. On the basis of the linear Stokes equations, we investigate theoretically and computationally the viscous azimuthal flow induced by the slow rotation of a small spherical particle located in the vicinity of a rigid no-slip disk or inside a gap between two coaxially positioned rigid no-slip disks of the same radius. We formulate the solution of the hydrodynamic problem as a mixed-boundary-value problem in the whole fluid domain, which we subsequently transform into a system of dual integral equations. Near a stationary disk, we show that the resulting integral equation can be reduced into an elementary Abel integral equation that admits a unique analytical solution. Between two coaxially positioned stationary disks, we demonstrate that the flow problem can be transformed into a system of two Fredholm integral equations of the first kind. The latter are solved by means of numerical approaches. Using our solution, we further investigate the effect of the disks on the slow rotational motion of a colloidal particle and provide expressions of the hydrodynamic mobility as a function of the system geometry. We compare our results with corresponding finite-element simulations and observe very good agreement., Comment: 14 pages, 9 figures, revised manuscript resubmitted to Phys. Fluids (invited article)

The ability of microswimmers to deploy optimal propulsion strategies is of paramount importance for their locomotory performance and survival at low Reynolds numbers. Although for perfectly spherical swimmers minimum dissipation requires a neutral type swimming, any departure from the spherical shape may lead the swimmer to adopt a new propulsion strategy, namely those of puller- or pusher-type swimming. In this study, by using the minimum dissipation theorem for microswimmers, we determine the flow field of an optimal nearly spherical swimmer, and show that indeed depending on the shape profile, the optimal swimmer can be a puller, pusher, or neutral. Using an asymptotic approach, we find that amongst all the modes of the shape function, only the third mode determines to leading order the swimming type of the optimal swimmer., Comment: 9 pages, 3 figures

Contrasting the well explored problem on how to steer a macroscopic agent like an airplane or a moon lander to optimally reach a target, "optimal microswimming", i.e. the quest for the optimal navigation strategy for microswimmers, remains unsolved. Here, we systematically explore this problem and show that the characteristic flow field of microswimmers crucially influences the required navigation strategy to reach a target fastest. The resulting optimal trajectories can have remarkable and non-intuitive shapes, which qualitatively differ from those of dry active particles or motile macroagents. Our results provide generic insights into the role of hydrodynamics and fluctuations on optimal navigation at the microscale and suggest that microorganisms might have survival advantages when strategically controlling their distance to remote walls. In particular, when fluctuations are present, choosing the optimal strategy, which appropriately respects hydrdynamics, can halve the time to reach the target compared to cases microswimmers head straight towards it.

A deep understanding of the physical interactions between nanoparticles and target cell membranes is important in designing efficient nanocarrier systems for drug delivery applications. Here, we present a theoretical framework to describe the hydrodynamic flow field induced by a point-force singularity (Stokeslet) directed parallel to a finite-sized elastic membrane endowed with shear and bending rigidities. We formulate the elastohydrodynamic problem as a mixed-boundary-value problem, which we then reduce into a well-behaved system of integro-differential equations. It follows that shear and bending linearly decouple so that the solution of the overall flow problem can be obtained by linear superposition of the contributions arising from these modes of deformation. Additionally, we probe the effect of the membrane on the hydrodynamic drag acting on a nearby particle, finding that, in a certain range of parameters, translational motion near an elastic membrane with only energetic resistance toward shear can, surprisingly, be sped up compared to bulk fluid. Our results may find applications in microrheological characterizations of colloidal systems near elastic confinements., Comment: To appear in J. Phys. Soc. Jpn

We investigate theoretically on the basis of the steady Stokes equations for a viscous incompressible fluid the flow induced by a Stokeslet located on the centre axis of two coaxially positioned rigid disks. The Stokeslet is directed along the centre axis. No-slip boundary conditions are assumed to hold at the surfaces of the disks. We perform the calculation of the associated Green's function in large parts analytically, reducing the spatial evaluation of the flow field to one-dimensional integrations amenable to numerical treatment. To this end, we formulate the solution of the hydrodynamic problem for the viscous flow surrounding the two disks as a mixed-boundary-value problem, which we then reduce into a system of four dual integral equations. We show the existence of viscous toroidal eddies arising in the fluid domain bounded by the two disks, manifested in the plane containing the centre axis through adjacent counterrotating eddies. Additionally, we probe the effect of the confining disks on the slow dynamics of a point-like particle by evaluating the hydrodynamic mobility function associated with axial motion. Thereupon, we assess the appropriateness of the commonly-employed superposition approximation and discuss its validity and applicability as a function of the geometrical properties of the system. Additionally, we complement our semi-analytical approach by finite-element computer simulations, which reveals a good agreement. Our results may find applications in guiding the design of microparticle-based sensing devices and electrokinetic transport in small scale capacitors.

Geometric confinements are frequently encountered in the biological world and strongly affect the stability, topology, and transport properties of active suspensions in viscous flow. Based on a far-field analytical model, the low-Reynolds-number locomotion of a self-propelled microswimmer moving inside a clean viscous drop or a drop covered with a homogeneously distributed surfactant, is theoretically examined. The interfacial viscous stresses induced by the surfactant are described by the well-established Boussinesq-Scriven constitutive rheological model. Moreover, the active agent is represented by a force dipole and the resulting fluid-mediated hydrodynamic couplings between the swimmer and the confining drop are investigated. We find that the presence of the surfactant significantly alters the dynamics of the encapsulated swimmer by enhancing its reorientation. Exact solutions for the velocity images for the Stokeslet and dipolar flow singularities inside the drop are introduced and expressed in terms of infinite series of harmonic components. Our results offer useful insights into guiding principles for the control of confined active matter systems and support the objective of utilizing synthetic microswimmers to drive drops for targeted drug delivery applications., Comment: 19 pages, 7 figures. Regular article contributed to the Topical Issue of the European Physical Journal E entitled "Physics of Motile Active Matter" edited by Gerhard Gompper, Clemens Bechinger, Holger Stark, and Roland G. Winkler

The navigation of micro-robots in complex flow environments is controlled by rheotaxis, the reorientation with respect to flow gradients. Here we demonstrate how payloads can be exploited to enhance the motion against flows. Using fully resolved hydrodynamic simulations, the mechanisms are described that allow micro-robots of different shapes to reorient upstream. We find that cargo pullers are the fastest at most flow strengths, but pushers feature a non-trivial optimum as a function of the counter flow strength. Moreover, the rheotactic performance can be maximised by tuning the micro-robot shape or cargo size. These results may be used to control micro-swimmer navigation, but they also apply to rheotaxis in microbial ecology and the prevention of bacterial contamination dynamics., Comment: 10 pages, 4 figures

Guiding active microswimmers by external fields to requested target locations is a promising strategy to realize complex transport on the microscale. To this end, one possibility consists of attaching the microswimmers to orientable passive components. Accordingly, we analyze theoretically, using a minimal model, the dynamics of a microswimmer when rigidly attached tangentially to a (significantly larger) microplatelet, here represented by a thin circular disk. On this way, we determine the flow field in the whole space induced by a Stokeslet that is located above the center of a spatially fixed rigid disk of no-slip surface conditions. Finally, we determine and analyze possible trajectories of the overall composite. To this end, the platelet is additionally endowed with a permanent magnetic moment, which allows to steer the motion of the whole composite by a homogeneous external magnetic field. As previous experimental studies suggest, related setups may be helpful to guide sperm cells to requested targets or for the purpose of coordinated drug delivery., Comment: 18 pages, 4 figures. To appear in Phys. Fluids

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., Comment: 20 pages, 3 figures. To appear in PRE. Abstract shortened to comply with the arXiv limit of 1920 characters

With the rapid advent of biomedical and biotechnological innovations, a deep understanding of the nature of interaction between nanomaterials and cell membranes, tissues, and organs, has become increasingly important. Active penetration of nanoparticles through cell membranes is a fascinating phenomenon that may have important implications in various biomedical and clinical applications. Using a fully analytical theory supplemented by particle-based computer simulations, the penetration process of an active particle through a planar two-dimensional elastic membrane is studied. The membrane is modeled as a self-assembled sheet of particles, uniformly arranged on a square lattice. A coarse-grained model is introduced to describe the mutual interactions between the membrane particles. The active penetrating particle is assumed to interact sterically with the membrane particles. State diagrams are presented to fully characterize the system behavior as functions of the relevant control parameters governing the transition between different dynamical states. Three distinct scenarios are identified. These compromise trapping of the active particle, penetration through the membrane with subsequent self-healing, in addition to penetration with permanent disruption of the membrane. The latter scenario is accompanied by a partial fragmentation of the membrane and creation of a hole of a size exceeding the interaction range of the membrane components. Our analytical theory is based on a combination of a perturbative expansion technique and a discrete-to-continuum formulation. Our approach might be helpful for the prediction of the transition threshold between the trapping and penetration in real-space experiments involving motile swimming bacteria or artificial active particles., Comment: 13 pages, 7 figures. Abstract shortened to comply with the arXiv limit of 1920 characters

An analytical method is proposed for computing the low-Reynolds-number hydrodynamic mobility function of a small colloidal particle asymmetrically moving inside a large spherical elastic cavity, the membrane of which is endowed with resistance toward shear and bending. In conjunction with the results obtained in the first part [Daddi-Moussa-Ider, L\"{o}wen, and Gekle, Eur. Phys. J. E 41, 104 (2018)], in which the axisymmetric motion normal to the surface of an elastic cavity is investigated, the general motion for an arbitrary force direction can be addressed. The elastohydrodynamic problem is formulated and solved using the classic method of images through expressing the hydrodynamic flow fields as a multipole expansion involving higher-order derivatives of the free-space Green's function. In the quasi-steady limit, we demonstrate that the particle self-mobility function of a particle moving tangent to the surface of the cavity is larger than that predicted inside a rigid stationary cavity of equal size. This difference is justified by the fact that a stationary rigid cavity introduces additional hindrance to the translational motion of the encapsulated particle, resulting in a reduction of its hydrodynamic mobility. Furthermore, the motion of the cavity is investigated, revealing that the translational pair (composite) mobility, which linearly couples the velocity of the elastic cavity to the force exerted on the solid particle, is solely determined by membrane shear properties. Our analytical predictions are favorably compared with fully-resolved computer simulations based on a completed-double-layer boundary integral method., Comment: 14 pages, 4 figures. This is a pre-print of an article published in The European Physical Journal E. The final authenticated version is available online at: https://doi.org/10.1140/epje/i2019-11853-4

Elastic confinements play an important role in many soft matter systems and affect the transport properties of suspended particles in viscous flow. On the basis of low-Reynolds-number hydrodynamics, we present an analytical theory of the axisymmetric flow induced by a point-force singularity (Stokeslet) directed along the symmetry axis of a finite-sized circular elastic membrane endowed with resistance toward shear and bending. The solution for the viscous incompressible flow surrounding the membrane is formulated as a mixed boundary value problem, which is then reduced into a system of dual integral equations on the inner and outer sides of the domain boundary. We show that the solution of the elastohydrodynamic problem can conveniently be expressed in terms of a set of inhomogeneous Fredholm integral equations of the second kind with logarithmic kernel. Basing on the hydrodynamic flow field, we obtain semi-analytical expressions of the hydrodynamic mobility function for the translational motion perpendicular to a circular membrane. The results are valid to leading-order in the ratio of particle radius to the distance separating the particle from the membrane. In the quasi-steady limit, we find that the particle mobility near a finite-sized membrane is always larger than that predicted near a no-slip disk of the same size. We further show that the bending-related contribution to the hydrodynamic mobility increases monotonically upon decreasing the membrane size, whereas the shear-related contribution displays a minimum value when the particle-membrane distance is equal to the membrane radius. Accordingly, the system behavior may be shear or bending dominated, depending on the geometric and elastic properties of the system. Our results may find applications in the field of nanoparticle-based sensing and drug delivery systems near elastic cell membranes.

The interaction between nano- or micro-sized particles and cell membranes is of crucial importance in many biological and biomedical applications such as drug and gene delivery to cells and tissues. During their cellular uptake, the particles can pass through cell membranes via passive endocytosis or by active penetration to reach a target cellular compartment or organelle. In this manuscript, we develop a simple model to describe the interaction of a self-driven spherical particle (moving through an effective constant active force) with a minimal membrane system, allowing for both penetration and trapping. We numerically calculate the state diagram of this system, the membrane shape, and its dynamics. In this context, we show that the active particle may either get trapped near the membrane or penetrates through it, where the membrane can either be permanently destroyed or recover its initial shape by self-healing. Additionally, we systematically derive a continuum description allowing to accurately predict most of our results analytically. This analytical theory helps identifying the generic aspects of our model, suggesting that most of its ingredients should apply to a broad range of membranes, from simple model systems composed of magnetic microparticles to lipid bilayers. Our results might be useful to predict mechanical properties of synthetic minimal membranes., Comment: 16 pages, 6 figures. Revised manuscript resubmitted to J. Chem. Phys

Several recent experiments investigate the orientational and transport behavior of self-driven bacteria and colloidal particles in nematic liquid crystals. Correspondingly, we study theoretically the dynamics of a minimal model microswimmer in a uniaxially anisotropic fluid. As a first step, the hydrodynamic Green's function providing the resulting fluid flow in response to a localized force acting on the anisotropic fluid is derived analytically. On this basis, the behavior of both puller- and pusher-type microswimmers in the anisotropic fluid is analyzed. Depending on the propulsion mechanism and the relative magnitude of different involved viscosities, we find alignment of the swimmers parallel or perpendicular to the anisotropy axis. Particularly, also an oblique alignment is identified under certain circumstances. The observed swimmer reorientation results from the hydrodynamic coupling between the self-induced fluid flow and the anisotropy of the surrounding fluid, which distorts the self-generated flow field. We support parts of our results by a simplified linear stability analysis. Our theoretical predictions are in qualitative agreement with recent experimental observations on swimming bacteria in nematic liquid crystals. They support the objective of utilizing the, possibly switchable, anisotropy of a host fluid to guide individual microswimmers and active particles along a requested path, enabling controlled active transport., Comment: 26 pages, 7 figures. Revised manuscript resubmitted to Phys. Rev. Fluids

We study theoretically and numerically the coupling and rotational hydrodynamic interactions between spherical particles near a planar elastic membrane that exhibits resistance towards shear and bending. Using a combination of the multipole expansion and Faxen's theorems, we express the frequency-dependent hydrodynamic mobility functions as a power series of the ratio of the particle radius to the distance from the membrane for the self mobilities, and as a power series of the ratio of the radius to the interparticle distance for the pair mobilities. In the quasi-steady limit of zero frequency, we find that the shear- and bending-related contributions to the particle mobilities may have additive or suppressive effects depending on the membrane properties in addition to the geometric configuration of the interacting particles relative to the confining membrane. To elucidate the effect and role of the change of sign observed in the particle self and pair mobilities, we consider an example involving a torque-free doublet of counterrotating particles near an elastic membrane. We find that the induced rotation rate of the doublet around its center of mass may differ in magnitude and direction depending on the membrane shear and bending properties. Near a membrane of only energetic resistance toward shear deformation, such as that of a certain type of elastic capsules, the doublet undergoes rotation of the same sense as observed near a no-slip wall. Near a membrane of only energetic resistance toward bending, such as that of a fluid vesicle, we find a reversed sense of rotation. Our analytical predictions are supplemented and compared with fully resolved boundary integral simulations where a very good agreement is obtained over the whole range of applied frequencies., Comment: 14 pages, 7 figures. Revised manuscript resubmitted to J. Chem. Phys

When an elastic object is dragged through a viscous fluid tangent to a rigid boundary, it experiences a lift force perpendicular to its direction of motion. An analogous lift mechanism occurs when a rigid symmetric object translates parallel to an elastic interface or a soft substrate. The induced lift force is attributed to an elastohydrodynamic coupling that arises from the breaking of the flow reversal symmetry induced by the elastic deformation of the translating object or the interface. Here we derive explicit analytical expressions for the quasi-steady state lift force exerted on a rigid spherical particle translating parallel to a finite-sized membrane exhibiting a resistance toward both shear and bending. Our analytical approach proceeds through the application of the Lorentz reciprocal theorem so as to obtain the solution of the flow problem using a perturbation technique for small deformations of the membrane. We find that the shear-related contribution to the normal force leads to an attractive interaction between the particle and the membrane. This emerging attractive force decreases quadratically with the system size to eventually vanish in the limit of an infinitely-extended membrane. In contrast, membrane bending leads to a repulsive interaction whose effect becomes more pronounced upon increasing the system size, where the lift force is found to diverge logarithmically for an infinitely-large membrane. The unphysical divergence of the bending-induced lift force can be rendered finite by regularizing the solution with a cut-off length beyond which the bending forces become subdominant to an external body force., Comment: 15 pages, 4 figures, 80 references. Under review

Geometric confinements are frequently encountered in soft matter systems and in particular significantly alter the dynamics of swimming microorganisms in viscous media. Surface-related effects on the motility of microswimmers can lead to important consequences in a large number of biological systems, such as biofilm formation, bacterial adhesion and microbial activity. On the basis of low-Reynolds-number hydrodynamics, we explore the state diagram of a three-sphere microswimmer under channel confinement in a slit geometry and fully characterize the swimming behavior and trajectories for neutral swimmers, puller- and pusher-type swimmers. While pushers always end up trapped at the channel walls, neutral swimmers and pullers may further perform a gliding motion and maintain a stable navigation along the channel. We find that the resulting dynamical system exhibits a supercritical pitchfork bifurcation in which swimming in the mid-plane becomes unstable beyond a transition channel height while two new stable limit cycles or fixed points that are symmetrically disposed with respect to the channel mid-height emerge. Additionally, we show that an accurate description of the averaged swimming velocity and rotation rate in a channel can be captured analytically using the method of hydrodynamic images, provided that the swimmer size is much smaller than the channel height., Comment: 19 pages, 9 figures. Article contributed to the Topical Issue of the Journal of Physics: Condensed Matter, entitled "Transport in Narrow Channels" edited by Paolo Malgaretti, Gleb Oshanin, and Julian Talbot

On the basis of the linear hydrodynamic equations, we present an analytical theory for the low-Reynolds-number motion of a solid particle moving inside a larger spherical elastic cavity which can be seen as a model system for a fluid vesicle. In the particular situation where the particle is concentric with the cavity, we use the stream function technique to find exact analytical solutions of the fluid motion equations on both sides of the elastic cavity. In this particular situation, we find that the solution of the hydrodynamic equations is solely determined by membrane shear properties and that bending does not play a role. For an arbitrary position of the solid particle within the spherical cavity, we employ the image solution technique to compute the axisymmetric flow field induced by a point force (Stokeslet). We then obtain analytical expressions of the leading order mobility function describing the fluid-mediated hydrodynamic interactions between the particle and confining elastic cavity. In the quasi-steady limit of vanishing frequency, we find that the particle self-mobility function is higher than that predicted inside a rigid no-slip cavity. Considering the cavity motion, we find that the pair-mobility function is determined only by membrane shear properties. Our analytical predictions are supplemented and validated by fully-resolved boundary integral simulations where a very good agreement is obtained over the whole range of applied forcing frequencies., Comment: 15 pages, 5 figures, 90 references. To appear in Eur. Phys. J. E

The hydrodynamic flow field generated by self-propelled active particles and swimming microorganisms is strongly altered by the presence of nearby boundaries in a viscous flow. Using a simple model three-linked sphere swimmer, we show that the swimming trajectories near a no-slip wall reveal various scenarios of motion depending on the initial orientation and the distance separating the swimmer from the wall. We find that the swimmer can either be trapped by the wall, completely escape, or perform an oscillatory gliding motion at a constant mean height above the wall. Using a far-field approximation, we find that, at leading order, the wall-induced correction has a source-dipolar or quadrupolar flow structure where the translational and angular velocities of the swimmer decay as inverse third and fourth powers with distance from the wall, respectively. The resulting equations of motion for the trajectories and the relevant order parameters fully characterize the transition between the states and allow for an accurate description of the swimming behavior near a wall. We demonstrate that the transition between the trapping and oscillatory gliding states is first order discontinuous, whereas the transition between the trapping and escaping states is continuous, characterized by non-trivial scaling exponents of the order parameters. In order to model the circular motion of flagellated bacteria near solid interfaces, we further assume that the spheres can undergo rotational motion around the swimming axis. We show that the general three-dimensional motion can be mapped onto a quasi-two-dimensional representational model by an appropriate redefinition of the order parameters governing the transition between the swimming states., Comment: 17 pages, 9 figures. Manuscript accepted for publication in J. Chem. Phys

Elastic confinements are an important component of many biological systems and dictate the transport properties of suspended particles under flow. In this chapter, we review the Brownian motion of a particle moving in the vicinity of a living cell whose membrane is endowed with a resistance towards shear and bending. The analytical calculations proceed through the computation of the frequency-dependent mobility functions and the application of the fluctuation-dissipation theorem. Elastic interfaces endow the system with memory effects that lead to a long-lived anomalous subdiffusive regime of nearby particles. In the steady limit, the diffusional behavior approaches that near a no-slip hard wall. The analytical predictions are validated and supplemented with boundary-integral simulations., Comment: 16 pages, 7 figures and 161 references. Contributed chapter to the flowing matter book

Elastic channels are an important component of many soft matter systems, in which hydrodynamic interactions with confining membranes determine the behavior of particles in flow. In this work, we derive analytical expressions for the Green's functions associated to a point-force (Stokeslet) directed parallel or perpendicular to the axis of an elastic cylindrical channel exhibiting resistance against shear and bending. We then compute the leading order self- and pair mobility functions of particles on the cylinder axis, finding that the mobilities are primarily determined by membrane shear and that bending does not play a significant role. In the quasi-steady limit of vanishing frequency, the particle self- and pair mobilities near a no-slip hard cylinder are recovered only if the membrane possesses a non-vanishing shear rigidity. We further compute the membrane deformation, finding that deformation is generally more pronounced in the axial (radial) directions, for the motion along (perpendicular) to the cylinder centerline, respectively. Our analytical calculations for the Green's functions in an elastic cylinder can serve as a fundamental building block for future studies and are verified by fully resolved boundary integral simulations where a very good agreement is obtained., Comment: 28 pages, 7 figures and 129 references. Manuscript under review

In this paper, we present an analytical calculation of the rotational mobility functions of a particle rotating on the centerline of an elastic cylindrical tube whose membrane exhibits resistance towards shearing and bending. We find that the correction to the particle rotational mobility about the cylinder axis depends solely on membrane shearing properties while both shearing and bending manifest themselves for the rotational mobility about an axis perpendicular to the cylinder axis. In the quasi-steady limit of vanishing frequency, the particle rotational mobility nearby a no-slip rigid cylinder is recovered only if the membrane possesses a non-vanishing resistance towards shearing. We further show that for the asymmetric rotation along the cylinder radial axis, a coupling between shearing and bending exists. Our analytical predictions are compared and validated with corresponding boundary integral simulations where a very good agreement is obtained., Comment: 23 pages, 7 figures and 107 references. Revised manuscript resubmitted to Acta Mech

In this paper, we derive analytical expressions for the leading-order hydrodynamic mobility of a small solid particle undergoing motion tangential to a nearby large spherical capsule whose membrane possesses resistance towards shearing and bending. Together with the results obtained in the first part (Daddi-Moussa-Ider and Gekle, Phys. Rev. E {\bfseries 95}, 013108 (2017)) where the axisymmetric motion perpendicular to the capsule membrane is considered, the solution of the general mobility problem is thus determined. We find that shearing resistance induces a low-frequency peak in the particle self-mobility, resulting from the membrane normal displacement in the same way, although less pronounced, to what has been observed for the axisymmetric motion. In the zero frequency limit, the self-mobility correction near a hard sphere is recovered only if the membrane has a non-vanishing resistance towards shearing. We further compute the particle in-plane mean-square displacement of a nearby diffusing particle, finding that the membrane induces a long-lasting subdiffusive regime. Considering capsule motion, we find that the correction to the pair-mobility function is solely determined by membrane shearing properties. Our analytical calculations are compared and validated with fully resolved boundary integral simulations where a very good agreement is obtained., Comment: 17 pages, 9 figures and 64 references. Manuscript accepted for publication in Phys. Rev. E

We use the image solution technique to compute the leading order frequency-dependent self-mobility function of a small solid particle moving perpendicular to the surface of a spherical capsule whose membrane possesses shearing and bending rigidities. Comparing our results with those obtained earlier for an infinitely extended planar elastic membrane, we find that membrane curvature leads to the appearance of a prominent additional peak in the mobility. This peak is attributed to the fact that the shear resistance of the curved membrane involves a contribution from surface-normal displacements which is not the case for planar membranes. In the vanishing frequency limit, the particle self-mobility near a no-slip hard sphere is recovered only when the membrane possesses a non-vanishing resistance towards shearing. We further investigate capsule motion, finding that the pair-mobility function is solely determined by membrane shearing properties. Our analytical predictions are validated by fully resolved boundary integral simulations where a very good agreement is obtained., Comment: 16 pages, 8 figures and 93 references. Manuscript accepted for publication in Phys. Rev. E

Using a fully analytical theory, we compute the leading order corrections to the translational, rotational and translation-rotation coupling mobilities of an arbitrary axisymmetric particle immersed in a Newtonian fluid moving near an elastic cell membrane that exhibits resistance towards stretching and bending. The frequency-dependent mobility corrections are expressed as general relations involving separately the particle's shape-dependent bulk mobility and the shape-independent parameters such as the membrane-particle distance, the particle orientation and the characteristic frequencies associated with shearing and bending of the membrane. This makes the equations applicable to an arbitrary-shaped axisymmetric particle provided that its bulk mobilities are known, either analytically or numerically. For a spheroidal particle, these general relations reduce to simple expressions in terms of the particle's eccentricity. We find that the corrections to the translation-rotation coupling mobility are primarily determined by bending, whereas shearing manifests itself in a more pronounced way in the rotational mobility. We demonstrate the validity of the analytical approximations by a detailed comparison with boundary integral simulations of a truly extended spheroidal particle. They are found to be in a good agreement over the whole range of applied frequencies., Comment: Published as: J. Fluid Mech. 811, 210-233 (2017)

We present an analytical calculation of the hydrodynamic interaction between two spherical particles near an elastic interface such as a cell membrane. The theory predicts the frequency dependent self- and pair-mobilities accounting for the finite particle size up to the 5th order in the ratio between particle diameter and wall distance as well as between diameter and interparticle distance. We find that particle motion towards a membrane with pure bending resistance always leads to mutual repulsion similar as in the well-known case of a hard-wall. In the vicinity of a membrane with shearing resistance, however, we observe an attractive interaction in a certain parameter range which is in contrast to the behavior near a hard wall. This attraction might facilitate surface chemical reactions. Furthermore, we show that there exists a frequency range in which the pair-mobility for perpendicular motion exceeds its bulk value, leading to short-lived superdiffusive behavior. Using the analytical particle mobilities we compute collective and relative diffusion coefficients. The appropriateness of the approximations in our analytical results is demonstrated by corresponding boundary integral simulations which are in excellent agreement with the theoretical predictions., Comment: 16 pages, 7 figures and 109 references. Manuscript accepted for publication in J. Chem. Phys

We study the motion of a solid particle immersed in a Newtonian fluid and confined between two parallel elastic membranes possessing shear and bending rigidity. The hydrodynamic mobility depends on the frequency of the particle motion due to the elastic energy stored in the membrane. Unlike the single-membrane case, a coupling between shearing and bending exists. The commonly used approximation of superposing two single-membrane contributions is found to give reasonable results only for motions in the parallel, but not in the perpendicular direction. We also compute analytically the membrane deformation resulting from the motion of the particle, showing that the presence of the second membrane reduces deformation. Using the fluctuation-dissipation theorem we compute the Brownian motion of the particle, finding a long-lasting subdiffusive regime at intermediate time scales. We finally assess the accuracy of the employed point-particle approximation via boundary-integral simulations for a truly extended particle. They are found to be in excellent agreement with the analytical predictions., Comment: 14 pages, 8 figures and 96 references. Revised version resubmitted to Phys. Fluids

The physical approach of a small particle (virus, medical drug) to the cell membrane represents the crucial first step before active internalization and is governed by thermal diffusion. Using a fully analytical theory we show that the stretching and bending of the elastic membrane by the approaching particle induces a memory in the system which leads to anomalous diffusion, even though the particle is immersed in a purely Newtonian liquid. For typical cell membranes the transient subdiffusive regime extends beyond 10 ms and can enhance residence times and possibly binding rates up to 50\%. Our analytical predictions are validated by numerical simulations., Comment: 13 pages and 5 figures. The Supporting Information is included. Manuscript accepted for publication in Phys. Rev. E

Understanding the dynamics of particles in turbulent flow is important in many environmental and industrial applications. In this paper, the statistics of particle pair orientation is numerically studied in homogeneous isotropic turbulent flow, with Taylor microscale Reynolds number of 300. It is shown that the Kolmogorov scaling fails to predict the observed probability density functions (PDFs) of the pair rotation rate and the higher order moments accurately. Therefore, a multifractal formalism is derived in order to include the intermittent behavior that is neglected in the Kolmogorov picture. The PDFs of finding the pairs at a given angular velocity for small relative separations reveals extreme events with stretched tails and high kurtosis values. Additionally, The PDFs are found to be less intermittent and follow a complementary error function distribution for larger separations., Comment: 16 pages, 3 figures, accepted for publication in European Journal of Mechanics / B Fluids

The environmental impact of pollutants and effluents discharged into the atmosphere or the oceans has led researchers to conduct studies related to this issue. Several works have been carried out in this context in order to reduce the effect on the local environment. These types of ejections in nature are modeled as jets in the presence of a density gradient. In this study we treated the problem of inclined round turbulent buoyant jets and plumes ejected in a homogeneous or stratified fluid, at rest or in motion. The prediction of the flow behavior, i.e. the evolution of its variables, is first treated theoretically from a model whose formalism is valid in both the Boussinesq case as well as in the non-Boussinesq general case. Solving the equations governing the plumes is performed numerically using a Runge-Kutta 4th order. To validate the model, laboratory experiments are performed with round jets of air and helium for a wide range of densities. The confrontation theory-experience aims here to fix the limits of validity of the theoretical model. Keywords: Buoyant jet, plume, non-Boussinesq, laser tomography, image processing, perturbation method., Comment: 6 pages, in French, 5 figures