Source author record

Lailai Zhu

Lailai Zhu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Optimising low-Reynolds-number predation via optimal control and reinforcement learning

We seek the best stroke sequences of a finite-size swimming predator chasing a non-motile point or finite--size prey at low Reynolds number. We use optimal control to seek the globally-optimal solutions for the former and RL for general situations. The predator is represented by a squirmer model that can translate forward and laterally, rotate and generate a stresslet flow. We identify the predator's best squirming sequences to achieve the time-optimal (TO) and efficiency-optimal (EO) predation. For a point prey, the TO squirmer executing translational motions favours a two-fold L-shaped trajectory that enables it to exploit the disturbance flow for accelerated predation; using a stresslet mode significantly expedites the EO predation, allowing the predator to catch the prey faster yet with lower energy consumption and higher predatory efficiency; the predator can harness its stresslet disturbance flow to suck the prey towards itself; compared to a translating predator, its compeer combining translation and rotation is less time--efficient, and the latter occasionally achieves the TO predation via retreating in order to advance. We also adopt RL to reproduce the globally-optimal predatory strategy of chasing a point prey, qualitatively capturing the crucial two--fold attribute of TO path. Using a numerically emulated RL environment, we explore the dependence of the optimal predatory path on the size of prey. Our results might provide useful information that help design synthetic microswimmers such as \textit{in vivo} medical micro-robots capable of capturing and approaching objects in viscous flows.

preprint2022arXiv

Viscoelastic levitation

The effects of viscoelasticity have been shown to manifest themselves via symmetry breaking. In this investigation, we show a novel phenomenon that arises from this idea. We observe that when a dense sphere is rotated near a wall (the rotation being aligned with the wall-normal direction and gravity), it levitates to a fixed distance away from the wall. Since the shear is larger in the gap (between the sphere and the wall) than in the open side of the sphere, the shear-induced elastic stresses are thus asymmetric, resulting in a net elastic vertical force that balances the weight of the sphere. We conduct experiments, theoretic models, and numerical simulations for rotating spheres of various sizes and densities in a Boger-type fluid. In the small Deborah number range, the results are collapsed into a universal trend by considering a dimensionless group of the ratio of elastic to gravitational forces.

preprint2020arXiv

Harnessing elasticity to generate self-oscillation via an electrohydrodynamic instability

Under a steady DC electric field of sufficient strength, a weakly conducting dielectric sphere in a dielectric solvent with higher conductivity can undergo spontaneous spinning (Quincke rotation) through a pitchfork bifurcation. We design an object composed of a dielectric sphere and an elastic filament. By solving an elasto-electro-hydrodynamic (EEH) problem numerically, we uncover an EEH instability exhibiting diverse dynamic responses. Varying the bending stiffness of the filament, the composite object displays three behaviours: a stationary state, undulatory swimming and steady spinning, where the swimming results from a self-oscillatory instability through a Hopf bifurcation. By conducting a linear stability analysis incorporating an elastohydrodynamic model, we theoretically predict the growth rates and critical conditions, which agree well with the numerical counterparts. We also propose a reduced model system consisting of a minimal elastic structure which reproduces the EEH instability. The elasto-viscous response of the composite structure is able to transform the pitchfork bifurcation into a Hopf bifurcation, leading to self-oscillation. Our results imply a new way of harnessing elastic media to engineer self-oscillations, and more generally, to manipulate and diversify the bifurcations and the corresponding instabilities. These ideas will be useful in designing soft, environmentally adaptive machines.

preprint2020arXiv

Particle motion nearby rough surfaces

We study the hydrodynamic coupling between particles and solid, rough boundaries characterized by random surface textures. Using the Lorentz reciprocal theorem, we derive analytical expressions for the grand mobility tensor of a spherical particle and find that roughness-induced velocities vary nonmonotonically with the characteristic wavelength of the surface. In contrast to sedimentation near a planar wall, our theory predicts continuous particle translation transverse and perpendicular to the applied force. Most prominently, this motion manifests itself in a variance of particle displacements that grows quadratically in time along the direction of the force. This increase is rationalized by surface roughness generating particle sedimentation closer to or farther from the surface, which entails a significant variability of settling velocities.

preprint2016arXiv

A pancake droplet translating in a Hele-Shaw cell: lubrication film and flow field

We adopt a boundary integral method to study the dynamics of a translating droplet confined in a Hele-Shaw cell in the Stokes regime. The droplet is driven by the motion of the ambient fluid with the same viscosity. We characterize the three-dimensional (3D) nature of the droplet interface and of the flow field. The interface develops an arc-shaped ridge near the rear-half rim with a protrusion in the rear and a laterally symmetric pair of higher peaks; this pair of protrusions has been identified by recent experiments (Huerre et al., Phys. Rev. Lett., vol. 115 (6), 2015, 064501) and predicted asymptotically (Burgess and Foster, Phys. Fluids A, vol. 2 (7), 1990, pp. 1105-1117). The mean film thickness is well predicted by the extended Bretherton model (Klaseboer et al., Phys. Fluids, vol. 26 (3), 2014, 032107) with fitting parameters. The flow in the streamwise wall-normal middle plane is featured with recirculating zones, which are partitioned by stagnation points closely resembling those of a two-dimensional droplet in a channel. Recirculation is absent in the wall-parallel, unconfined planes, in sharp contrast to the interior flow inside a moving droplet in free space. The preferred orientation of the recirculation results from the anisotropic confinement of the Hele-Shaw cell. On these planes, we identify a dipolar disturbance flow field induced by the travelling droplet and its $1/r^2$ spatial decay is confirmed numerically. We pinpoint counter-rotating streamwise vortex structures near the lateral interface of the droplet, further highlighting the complex 3D flow pattern.

preprint2016arXiv

Motion of an elastic capsule in a constricted microchannel

We study the motion of an elastic capsule through a microchannel characterized by a localized constriction. We consider a capsule with a stress-free spherical shape and impose its steady state configuration in an infinitely long straight channel as the initial condition for our calculations. We report how the capsule deformation, velocity, retention time, and maximum stress of the membrane are affected by the capillary number, Ca, and the constriction shape. We estimate the deformation by measuring the variation of the three-dimensional surface area and a series of alternative quantities easier to extract from experiments. These are the Taylor parameter, the perimeter and the area of the capsule in the spanwise plane. We find that the perimeter is the quantity that reproduces the behavior of the three-dimensional surface area the best. We observe that, in general, area-deformation correlated quantities grow linearly with Ca, while velocity-correlated quantities saturate for large Ca but display a steeper increase for small Ca. The velocity of the capsule divided by the velocity of the flow displays, surprisingly, two different qualitative behaviors for small and large capillary numbers. Finally, we report that longer constrictions and spanwise wall bounded (versus spanwise periodic) domains cause larger deformations and velocities. If the deformation and velocity in the spanwise wall bounded domains are rescaled by the initial equilibrium deformation and velocity, their behavior is undistinguishable from that in a periodic domain. In contrast, a remarkably different behavior is reported in sinusoidally shaped and smoothed rectangular constrictions indicating that the capsule dynamics is particularly sensitive to abrupt changes in the cross section. In a smoothed rectangular constriction larger deformations and velocities occur over a larger distance.

preprint2015arXiv

Squirming through shear-thinning fluids

Many microorganisms find themselves immersed in fluids displaying non-Newtonian rheological properties such as viscoelasticity and shear-thinning viscosity. The effects of viscoelasticity on swimming at low Reynolds numbers have already received considerable attention, but much less is known about swimming in shear-thinning fluids. A general understanding of the fundamental question of how shear-thinning rheology influences swimming still remains elusive. To probe this question further, we study a spherical squirmer in a shear-thinning fluid using a combination of asymptotic analysis and numerical simulations. Shear-thinning rheology is found to affect a squirming swimmer in nontrivial and surprising ways; we predict and show instances of both faster and slower swimming depending on the surface actuation of the squirmer. We also illustrate that while a drag and thrust decomposition can provide insights into swimming in Newtonian fluids, extending this intuition to problems in complex media can prove problematic.

preprint2015arXiv

The dynamics of a capsule in a wall-bounded oscillating shear flow

The motion of an initially spherical capsule in a wall-bounded oscillating shear flow is investigated via an accelerated boundary integral implementation. The neo-Hookean model is used as the constitutive law of the capsule membrane. The maximum wall-normal migration is observed when the oscillation period of the imposed shear is of the order of the relaxation time of the elastic membrane; hence, the optimal capillary number scales with the inverse of the oscillation frequency and the ratio agrees well with the theoretical prediction in the limit of high-frequency oscillation. The migration velocity decreases monotonically with the frequency of the applied shear and the capsule-wall distance. We report a significant correlation between the capsule lateral migration and the normal stress difference induced in the flow. The periodic variation of the capsule deformation is roughly in phase with that of the migration velocity and normal stress difference, with twice the frequency of the imposed shear. The maximum deformation increases linearly with the membrane elasticity before reaching a plateau at higher capillary numbers when the deformation is limited by the time over which shear is applied in the same direction and not by the membrane deformability. The maximum membrane deformation scales as the distance to the wall to the power 1/3 as observed for capsules and droplets in near-wall steady shear flows.

preprint2015arXiv

The motion of a deforming capsule through a corner

A three-dimensional deformable capsule convected through a square duct with a corner is studied via numerical simulations. We develop an accelerated boundary integral implementation adapted to general geometries and boundary conditions. A global spectral method is adopted to resolve the dynamics of the capsule membrane developing elastic tension according to the neo-Hookean constitutive law and bending moments in an inertialess flow. The simulations show that the trajectory of the capsule closely follows the underlying streamlines independently of the capillary number. The membrane deformability, on the other hand, significantly influences the relative area variations, the advection velocity and the principal tensions observed during the capsule motion. The evolution of the capsule velocity displays a loss of the time-reversal symmetry of Stokes flow due to the elasticity of the membrane. The velocity decreases while the capsule is approaching the corner as the background flow does, reaches a minimum at the corner and displays an overshoot past the corner due to the streamwise elongation induced by the flow acceleration in the downstream branch. This velocity overshoot increases with confinement while the maxima of the major principal tension increase linearly with the inverse of the duct width. Finally, the deformation and tension of the capsule are shown to decrease in a curved corner.

preprint2015arXiv

The stability of a rising droplet: an inertialess nonmodal growth mechanism

Prior modal stability analysis (Kojima et al., Phys. Fluids, vol. 27, 1984) predicted that a rising or sedimenting droplet in a viscous fluid is stable in the presence of surface tension no matter how small, in contrast to experimental and numerical results. By performing a non-modal stability analysis, we demonstrate the potential for transient growth of the interfacial energy of a rising droplet in the limit of inertialess Stokes equations. The predicted critical capillary numbers for transient growth agree well with those for unstable shape evolution of droplets found in the direct numerical simulations of Koh & Leal (Phys. Fluids, vol. 1, 1989). Boundary integral simulations are used to delineate the critical amplitude of the most destabilizing perturbations. The critical amplitude is negatively correlated with the linear optimal energy growth, implying that the transient growth is responsible for reducing the necessary perturbation amplitude required to escape the basin of attraction of the spherical solution.

preprint2014arXiv

A microfluidic device to sort capsules by deformability: A numerical study

Guided by extensive numerical simulations, we propose a microfluidic device that can sort elastic capsules by their deformability. The device consists of a duct embedded with a semi-cylindrical obstacle, and a diffuser which further enhances the sorting capability. We demonstrate that the device can operate reasonably well under changes in the initial position of the the capsule. The efficiency of the device remains essentially unaltered under small changes of the obstacle shape (from semi-circular to semi-elliptic cross-section). Confinement along the direction perpendicular to the plane of the device increases its efficiency. This work is the first numerical study of cell sorting by a realistic microfluidic device.

preprint2013arXiv

Low-Reynolds number swimming in a capillary tube

We use the boundary element method to study the low-Reynolds number locomotion of a spherical model microorganism in a circular tube. The swimmer propels itself by tangen- tial or normal surface motion in a tube whose radius is on the order of the swimmer size. Hydrodynamic interactions with the tube walls significantly affect the average swimming speed and power consumption of the model microorganism. In the case of swimming parallel to the tube axis, the locomotion speed is always reduced (resp. increased) for swimmers with tangential (resp. normal) deformation. In all cases, the rate of work nec- essary for swimming is increased by confinement. Swimmers with no force-dipoles in the far field generally follow helical trajectories, solely induced by hydrodynamic interactions with the tube walls, and in qualitative agreement with recent experimental observations for Paramecium. Swimmers of the puller type always display stable locomotion at a lo- cation which depends on the strength of their force dipoles: swimmers with weak dipoles (small α) swim in the centre of the tube while those with strong dipoles (large α) swim near the walls. In contrast, pusher swimmers and those employing normal deformation are unstable and end up crashing into the walls of the tube. Similar dynamics is observed for swimming into a curved tube. These results could be relevant for the future design of artificial microswimmers in confined geometries.

preprint2013arXiv

Micropropulsion and microrheology in complex fluids via symmetry breaking

Many biological fluids have polymeric microstructures and display non-Newtonian rheology. We take advantage of such nonlinear fluid behavior and combine it with geometrical symmetry-breaking to design a novel small-scale propeller able to move only in complex fluids. Its propulsion characteristics are explored numerically in an Oldroyd-B fluid for finite Deborah numbers while the small Deborah number limit is investigated analytically using a second-order fluid model. We then derive expressions relating the propulsion speed to the rheological properties of the complex fluid, allowing thus to infer the normal stress coefficients in the fluid from the locomotion of the propeller. Our simple mechanism can therefore be used either as a non-Newtonian micro-propeller or as a micro-rheometer.

preprint2012arXiv

Locomotion by tangential deformation in a polymeric fluid

In several biologically relevant situations, cell locomotion occurs in polymeric fluids with Weissenberg {number} larger than one. Here we present results of three-dimensional numerical simulations for the steady locomotion of a self-propelled body in a model polymeric (Giesekus) fluid at low Reynolds number. Locomotion is driven by steady tangential deformation at the surface of the body (so-called squirming motion). In the case of a spherical squirmer, we show that the swimming velocity is systematically less than that in a Newtonian fluid, with a minimum occurring for Weissenberg numbers of order one. The rate of work done by the swimmer always goes up compared to that occurring in the Newtonian solvent alone, but is always lower than the power necessary to swim in a Newtonian fluid with the same viscosity. The swimming efficiency, defined as the ratio between the rate of work necessary to pull the body at the swimming speed in the same fluid and the rate of work done by swimming, is found to always be increased in a polymeric fluid. Further analysis reveals that polymeric stresses break the Newtonian front-back symmetry in the flow profile around the body. In particular, a strong negative elastic wake is present behind the swimmer, which correlates with strong polymer stretching, and its intensity increases with Weissenberg number and viscosity contrasts. {The velocity induced by the squirmer is found to decay in space faster than} in a Newtonian flow, with a strong {dependence} on the polymer relaxation time and viscosity. Our computational results are also extended to prolate spheroidal swimmers and smaller polymer stretching are obtained for slender shapes compared to bluff swimmers. The swimmer with an aspect ratio of two is found to be the most hydrodynamically efficient.