Lailai Zhu

Self-propulsion in viscoelastic fluids: Pushers vs. pullers

We use numerical simulations to address locomotion at zero Reynolds number in viscoelastic (Giesekus) fluids. The swimmers are assumed to be spherical, to self-propel using tangential surface deformation, and the computations are implemented using a finite element method. The emphasis of the study is on the change of the swimming kinematics, energetics, and flow disturbance from Newtonian to viscoelastic, and on the distinction between pusher and puller swimmers. In all cases, the viscoelastic swimming speed is below the Newtonian one, with a minimum obtained for intermediate values of the Weissenberg number, We. An analysis of the flow field places the origin of this swimming degradation in non-Newtonian elongational stresses. The power required for swimming is also systematically below the Newtonian power, and always a decreasing function of We. A detail energetic balance of the swimming problem points at the polymeric part of the stress as the primary We-decreasing energetic contribution, while the contributions of the work done by the swimmer from the solvent remain essentially We-independent. In addition, we observe negative values of the polymeric power density in some flow regions, indicating positive elastic work by the polymers on the fluid. The hydrodynamic efficiency, defined as the ratio of the useful to total rate of work, is always above the Newtonian case, with a maximum relative value obtained at intermediate Weissenberg numbers. Finally, the presence of polymeric stresses leads to an increase of the rate of decay of the flow velocity in the fluid, and a decrease of the magnitude of the stresslet governing the magnitude of the effective bulk stress in the fluid.

Low-Reynolds-number swimming in a capillary tube

In this thesis, simulations are performed to study the motion ofindividual cells in flow, focusing on the hydrodynamics of actively swimming cells likethe self-propelling microorganisms, and of passively advected objects like the red bloodcells. In particular, we develop numerical tools to address the locomotion ofmicroswimmers in viscoelastic fluids and complex geometries, as well as the motion ofdeformable capsules in micro-fluidic flows.For the active movement, the squirmer is used as our model microswimmer. The finiteelement method is employed to study the influence of the viscoelasticity of fluid on theperformance of locomotion. A boundary element method is implemented to study swimmingcells inside a tube. For the passive counterpart, the deformable capsule is chosen as the modelcell. An accelerated boundary integral method code is developed to solve thefluid-structure interaction, and a global spectral method is incorporated to handle theevolving cell surface and its corresponding membrane dynamics.We study the locomotion of a neutral squirmer with anemphasis on the change of swimming kinematics, energetics, and flowdisturbance from Newtonian to viscoelastic fluid. We also examine the dynamics of differentswimming gaits resulting in different patterns of polymer deformation, as well as theirinfluence on the swimming performance. We correlate the change of swimming speed withthe extensional viscosity and that of power consumption with the phase delay of viscoelasticfluids. Moreover, we utilise the boundary element method to simulate the swimming cells in astraight and torus-like bent tube, where the tube radius is a few times the cell radius. Weinvestigate the effect of tube confinement to the swimming speed and power consumption. Weanalyse the motions of squirmers with different gaits, which significantly affect thestability of the motion. Helical trajectories are produced for a neutralsquirmer swimming, in qualitative agreement with experimental observations, which can beexplained by hydrodynamic interactions alone.We perform simulations of a deformable capsule in micro-fluidic flows. We look atthe trajectory and deformation of a capsule through a channel/duct with a corner. Thevelocity of capsule displays an overshoot as passing around the corner, indicating apparentviscoelasticity induced by the interaction between the deformable membrane and viscousflow. A curved corner is found to deform the capsule less than the straight one. In addition, we propose a new cell sorting device based on the deformability of cells. Weintroduce carefully-designed geometric features into the flow to excite thehydrodynamic interactions between the cell and device. This interaction varies andclosely depends on the cell deformability, the resultant difference scatters the cellsonto different trajectories. Our high-fidelity computations show that the new strategy achievesa clear and robust separation of cells. We finally investigate the motion of capsule in awall-bounded oscillating shear flow, to understand the effect of physiological pulsation to thedeformation and lateral migration of cells. We observe the lateral migration velocity of a cellvaries non-monotonically with its deformability.

Locomotion by tangential deformation in a polymeric fluid

In several biologically relevant situations, cell locomotion occurs in polymeric fluids with Weissenberg number larger than 1. 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 (the 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 1. 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.

Micropropulsion and microrheology in complex fluids via symmetry breaking

Related Articles Characterization and simulation of the first extensional mode of rectangular micro-plates in liquid media Appl. Phys. Lett. 101, 151904 (2012) A microfluidic platform for measuring electrical activity across cells Biomicrofluidics 6, 034121 (2012) Wall-less liquid pathways formed with three-dimensional microring arrays Appl. Phys. Lett. 101, 114108 (2012) Acoustophoretic microfluidic chip for sequential elution of surface bound molecules from beads or cells Biomicrofluidics 6, 034115 (2012) Sound wave resonances in micro-electro-mechanical systems devices vibrating at high frequencies according to the kinetic theory of gases Phys. Fluids 24, 092001 (2012)

Shape optimization of the diffuser blade of an axial blood pump by computational fluid dynamics.

Computational fluid dynamics (CFD) has been a viable and effective way to predict hydraulic performance, flow field, and shear stress distribution within a blood pump. We developed an axial blood pump with CFD and carried out a CFD-based shape optimization of the diffuser blade to enhance pressure output and diminish backflow in the impeller-diffuser connecting region at a fixed design point. Our optimization combined a computer-aided design package, a mesh generator, and a CFD solver in an automation environment with process integration and optimization software. A genetic optimization algorithm was employed to find the pareto-optimal designs from which we could make trade-off decisions. Finally, a set of representative designs was analyzed and compared on the basis of the energy equation. The role of the inlet angle of the diffuser blade was analyzed, accompanied by its relationship with pressure output and backflow in the impeller-diffuser connecting region.

Squirming through shear-thinning fluids

Many micro-organisms 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 non-trivial 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.

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.

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 best reproduces the behavior of the three-dimensional surface area. This is maximum at the centre of the constriction and shows a second peak after it, whose location depends on the Ca number. 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.Graphical abstract

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.

HEMOLYSIS ANALYSIS OF AXIAL BLOOD PUMPS WITH VARIOUS STRUCTURE IMPELLERS

Low hemolysis is an important factor for axial blood pumps that has been used in patients with heart failure. The structure of impellers plays a key role in the hemolytic properties of axial blood pumps. Axial blood pumps with various structure impellers exhibit different hemolytic characteristic. In the present study, we aimed to investigate the type of impellers structures in axial blood pumps that contain the best low hemolytic properties. Also, it is expensive and time-consuming to validate the axial blood pumps hemolytic property by in vivo experiments. Therefore, in the present study, the numerical method was applied to analyze the hemolytic property in a blood pump. Specifically, the hemolysis of the pump was calculated by using a forward Euler approach based on the changes in shear stress and related exposure times along the particle trace lines. The different vane structures and rotational speed that affect hemolysis were analyzed and compared. The results showed that long–short alternant vanes exhibited the best hemolytic property which could be utilized in the optimization design of axial blood pumps.