Tank-treading, swinging, and tumbling of liquid-filled elastic capsules in shear flow

In silico modeling of patient-specific blood rheology in type 2 diabetes mellitus

Increased blood viscosity in type 2 diabetes mellitus (T2DM) is a risk factor for the development of insulin resistance and diabetes-related vascular complications; however, individuals with T2DM exhibit heterogeneous hemorheological properties, including cell deformation and aggregation. Using a multiscale red blood cell (RBC) model with key parameters derived from patient-specific data, we present a computational study of the rheological properties of blood from individual patients with T2DM. Specifically, one key model parameter, which determines the shear stiffness of the RBC membrane (μ) is informed by the high-shear-rate blood viscosity of patients with T2DM. At the same time, the other, which contributes to the strength of the RBC aggregation interaction (D0), is derived from the low-shear-rate blood viscosity of patients with T2DM. The T2DM RBC suspensions are simulated at different shear rates, and the predicted blood viscosity is compared with clinical laboratory-measured data. The results show that the blood viscosity obtained from clinical laboratories and computational simulations are in agreement at both low and high shear rates. These quantitative simulation results demonstrate that the patient-specific model has truly learned the rheological behavior of T2DM blood by unifying the mechanical and aggregation factors of the RBCs, which provides an effective way to extract quantitative predictions of the rheological properties of the blood of individual patients with T2DM.

Rheology and structure of elastic capsule suspensions within rectangular channels

Three-dimensional simulations of the pressure-driven flow dynamics of elastic capsule suspensions within both slit and rectangular cross-section channels are presented. The simulations utilize the Immersed Boundary Method and the Lattice-Boltzmann Method models. The capsule volume fraction is fixed at 0.1 (i.e., a semi-dilute suspension), while the channel Reynolds number (Re), the capillary number (Ca), and the cross-sectional channel dimensions are systematically varied. Comparing results for slit and rectangular channels, it is found that multi-directional confinement hinders inertial focusing due to the capsule-free layers that develop in the two transverse directions. Furthermore, the thicknesses of the capsule-free layers in the two transverse directions differ when the height and width of the channel are not equal. Both the size and aspect ratio of the channel impact the apparent viscosity. It is found that square channels exhibit maximal viscosity and that holding one dimension fixed while increasing or decreasing the other results in a decrease in viscosity. The results therefore represent an expansion of the Fahraeus-Lindqvist effect from 1D cylindrical channels to 2D rectangular channels.

αI-spectrin represents evolutionary optimization of spectrin for red blood cell deformability

Spectrin tetramers of the membranes of enucleated mammalian erythrocytes play a critical role in red blood cell survival in circulation. One of the spectrins, αI, emerged in mammals with enucleated red cells after duplication of the ancestral α-spectrin gene common to all animals. The neofunctionalized αI-spectrin has moderate affinity for βI-spectrin, whereas αII-spectrin, expressed in nonerythroid cells, retains ancestral characteristics and has a 10-fold higher affinity for βI-spectrin. It has been hypothesized that this adaptation allows for rapid make and break of tetramers to accommodate membrane deformation. We have tested this hypothesis by generating mice with high-affinity spectrin tetramers formed by exchanging the site of tetramer formation in αI-spectrin (segments R0 and R1) for that of αII-spectrin. Erythrocytes with αIIβI presented normal hematologic parameters yet showed increased thermostability, and their membranes were significantly less deformable; under low shear forces, they displayed tumbling behavior rather than tank treading. The membrane skeleton is more stable with αIIβI and shows significantly less remodeling under deformation than red cell membranes of wild-type mice. These data demonstrate that spectrin tetramers undergo remodeling in intact erythrocytes and that this is required for the normal deformability of the erythrocyte membrane. We conclude that αI-spectrin represents evolutionary optimization of tetramer formation: neither higher-affinity tetramers (as shown here) nor lower affinity (as seen in hemolytic disease) can support the membrane properties required for effective tissue oxygenation in circulation.

A neural network-based algorithm for high-throughput characterisation of viscoelastic properties of flowing microcapsules

Microcapsules, consisting of a liquid droplet enclosed by a viscoelastic membrane, have a wide range of biomedical and pharmaceutical applications and also serve as a popular mechanical model for biological cells. In this study, we develop a novel high throughput approach, by combining a machine learning method with a high-fidelity mechanistic capsule model, to accurately predict the membrane elasticity and viscosity of microcapsules from their dynamic deformation when flowing in a branched microchannel. The machine learning method consists of a deep convolutional neural network (DCNN) connected by a long short-term memory (LSTM) network. We demonstrate that with a superior prediction accuracy the present hybrid DCNN-LSTM network can still be faster than a conventional inverse method by five orders of magnitude, and can process thousands of capsules per second. We also show that the hybrid network has fewer restrictions compared with a simple DCNN.

Multiplicity of stable orbits for deformable prolate capsules in shear flow

This work investigates the orbital dynamics of a fluid-filled deformable prolate capsule in unbounded simple shear flow at zero Reynolds number using direct simulations. The motion of the capsule is simulated using a model that incorporates shear elasticity, area dilatation, and bending resistance. Here the deformability of the capsule is characterized by the nondimensional capillary number Ca, which represents the ratio of viscous stresses to elastic restoring stresses on the capsule. For a capsule with small bending stiffness, at a given Ca, the orientation converges over time towards a unique stable orbit independent of the initial orientation. With increasing Ca, four dynamical modes are found for the stable orbit, namely, rolling, wobbling, oscillating-swinging, and swinging. On the other hand, for a capsule with large bending stiffness, multiplicity in the orbit dynamics is observed. When the viscosity ratio λ ≲ 1, the long-axis of the capsule always tends towards a stable orbit in the flow-gradient plane, either tumbling or swinging, depending on Ca. When λ ≳ 1, the stable orbit of the capsule is a tumbling motion at low Ca, irrespective of the initial orientation. Upon increasing Ca, there is a symmetry-breaking bifurcation away from the tumbling orbit, and the capsule is observed to adopt multiple stable orbital modes including nonsymmetric precessing and rolling, depending on the initial orientation. As Ca further increases, the nonsymmetric stable orbit loses existence at a saddle-node bifurcation, and rolling becomes the only attractor at high Ca, whereas the rolling state coexists with the nonsymmetric state at intermediate values of Ca. A symmetry-breaking bifurcation away from the rolling orbit is also found upon decreasing Ca. The regime with multiple attractors becomes broader as the aspect ratio of the capsule increases, while narrowing as viscosity ratio increases. We also report the particle contribution to the stress, which also displays multiplicity.

Numerical simulation of deformable particles in a Coulter counter

In Coulter counters, cells counting and volumetry are achieved by monitoring their electrical print when they flow through a sensing zone. However, the volume measurement may be impaired by the cell dynamics, which may be difficult to control. In this paper, numerical simulations of the dynamics and electrical signature of red blood cells in a Coulter counter are presented, accounting for the deformability of the cells. In particular, a specific numerical pipeline is developed to overcome the challenge of the multi-scale nature of the problem. It consists in segmenting the whole computation of the cell dynamics and electrical response in a series of dedicated computations, with a saving of one order of magnitude in computational time. This numerical pipeline is used with rigid spheres and deformable red blood cells in an industrial Coulter counter geometry, and compared with experimental measurements. The simulations not only reproduce electrical signatures typical of those measured experimentally, but also allow an analysis of the electrical signature in terms of the heterogeneity of the electrical field and dynamics of the particles in the measurement zone. This study provides a methodology for computing the sizing of rigid or deformable particles by Coulter counters, opening the way to a better understanding of cells signatures in such devices.

Dynamic characteristics of a deformable capsule in a simple shear flow

The dynamic characteristics of a two-dimensional deformable capsule in a simple shear flow are studied with an immersed boundary-lattice Boltzmann method. Simulations are conducted by varying the Reynolds number (Re) from 0.0125 to 2000 and the dimensionless shear rate (G) from 0.001 to 0.5. The G-Re plane can be divided into four regions according to the deformation dependence on the parameters considered: viscous dominant, inertia dominant, transitional, and anomalous regions. There are four typical dynamic behaviors over the G-Re plane: steady deformation, prerupture state, quasisteady deformation, and continuous elongation. Analysis indicates that the pressure distribution and its variations due to the interplay of the fluid inertia force, the viscous shear stress, and the membrane elastic force determines the complex behaviors of the capsule. The effects of the bending rigidity and the internal-to-external viscosity ratio on the dynamics of the capsule are further studied. It is found that the capsule experiences smaller deformation when the higher bending rigidity is included, and the low bending rigidity does not have a remarkable influence on the capsule deformation. The capsule normally experiences smaller deformation due to the increase of the internal-to-external viscosity ratio.

Elastic particle deformation in rectangular channel flow as a measure of particle stiffness

In this study, we experimentally observed and characterized soft elastic particle deformation in confined flow in a microchannel with a rectangular cross-section. Hydrogel microparticles of pNIPAM were produced using two different concentrations of crosslinker. This resulted in particles with two different shear moduli of 13.3 ± 5.5 Pa and 32.5 ± 15.7 Pa and compressive moduli of 66 ± 10 Pa and 79 ± 15 Pa, respectively, as measured by capillary micromechanics. Under flow, the particle shapes transitioned from circular to egg, triangular, arrowhead, and ultimately parachute shaped with increasing shear rate. The shape changes were reversible, and deformed particles relaxed back to circular/spherical in the absence of flow. The thresholds for each shape transition were quantified using a non-dimensional radius of curvature at the tip, particle deformation, circularity, and the depth of the concave dimple at the trailing edge. Several of the observed shapes were distinct from those previously reported in the literature for vesicles and capsules; the elastic particles had a narrower leading tip and a lower circularity. Due to variations in the shear moduli between particles within a batch of particles, each flow rate corresponded to a small but finite range of capillary number (Ca) and resulted in a series of shapes. By arranging the images on a plot of Ca versus circularity, a direct correlation was developed between shape and Ca and thus between particle deformation and shear modulus. As the shape was very sensitive to differences in shear modulus, particle deformation in confined flow may allow for better differentiation of microparticle shear modulus than other methods.

Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

We review recent advances in multiscale modeling of the biomechanical characteristics of red blood cells (RBCs) in hematological diseases, and their relevance to the structure and dynamics of defective RBCs. We highlight examples of successful simulations of blood disorders including malaria and other hereditary disorders, such as sickle-cell anemia, spherocytosis, and elliptocytosis.

Hybrid smoothed dissipative particle dynamics and immersed boundary method for simulation of red blood cells in flows

In biofluid flow systems, often the flow problems of fluids of complex structures, such as the flow of red blood cells (RBCs) through complex capillary vessels, need to be considered. The smoothed dissipative particle dynamics (SDPD), a particle-based method, is one of the easy and flexible methods to model such complex structure fluids. It couples the best features of the smoothed particle hydrodynamics (SPH) and dissipative particle dynamics (DPD), with parameters having specific physical meaning (coming from SPH discretization of the Navier-Stokes equations), combined with thermal fluctuations in a mesoscale simulation, in a similar manner to the DPD. On the other hand, the immersed boundary method (IBM), a preferred method for handling fluid-structure interaction problems, has also been widely used to handle the fluid-RBC interaction in RBC simulations. In this paper, we aim to couple SDPD and IBM together to carry out the simulations of RBCs in complex flow problems. First, we develop the SDPD-IBM model in details, including the SDPD model for the evolving fluid flow, the RBC model for calculating RBC deformation force, the IBM for treating fluid-RBC interaction, and the solid boundary treatment model as well. We then conduct the verification and validation of the combined SDPD-IBM method. Finally, we demonstrate the capability of the SDPD-IBM method by simulating the flows of RBCs in rectangular, cylinder, curved, bifurcated, and constricted tubes, respectively.

Hydrodynamic mobility of a solid particle near a spherical elastic membrane. II. Asymmetric motion

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 toward shearing and bending. Together with the results obtained in the first part [Daddi-Moussa-Ider and Gekle, Phys. Rev. E 95, 013108 (2017)2470-004510.1103/PhysRevE.95.013108], 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 nonvanishing resistance toward shearing. We further compute the 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.

Lattice Boltzmann Simulation of Healthy and Defective Red Blood Cell Settling in Blood Plasma

In this paper, an attempt has been made to study sedimentation of a red blood cell (RBC) in a plasma-filled tube numerically. Such behaviors are studied for a healthy and a defective cell which might be created due to human diseases, such as diabetes, sickle-cell anemia, and hereditary spherocytosis. Flow-induced deformation of RBC is obtained using finite-element method (FEM), while flow and fluid-membrane interaction are handled using lattice Boltzmann (LB) and immersed boundary methods (IBMs), respectively. The effects of RBC properties as well as its geometry and orientation on its sedimentation rate are investigated and discussed. The results show that decreasing frontal area of an RBC and/or increasing tube diameter results in a faster settling. Comparison of healthy and diabetic cells reveals that less cell deformability leads to slower settling. The simulation results show that the sicklelike and spherelike RBCs have lower settling velocity as compared with a biconcave discoid cell.

Deformation of a Capsule in a Power-Law Shear Flow

An immersed boundary-lattice Boltzmann method is developed for fluid-structure interactions involving non-Newtonian fluids (e.g., power-law fluid). In this method, the flexible structure (e.g., capsule) dynamics and the fluid dynamics are coupled by using the immersed boundary method. The incompressible viscous power-law fluid motion is obtained by solving the lattice Boltzmann equation. The non-Newtonian rheology is achieved by using a shear rate-dependant relaxation time in the lattice Boltzmann method. The non-Newtonian flow solver is then validated by considering a power-law flow in a straight channel which is one of the benchmark problems to validate an in-house solver. The numerical results present a good agreement with the analytical solutions for various values of power-law index. Finally, we apply this method to study the deformation of a capsule in a power-law shear flow by varying the Reynolds number from 0.025 to 0.1, dimensionless shear rate from 0.004 to 0.1, and power-law index from 0.2 to 1.8. It is found that the deformation of the capsule increases with the power-law index for different Reynolds numbers and nondimensional shear rates. In addition, the Reynolds number does not have significant effect on the capsule deformation in the flow regime considered. Moreover, the power-law index effect is stronger for larger dimensionless shear rate compared to smaller values.

Dynamic and rheological properties of soft biological cell suspensions

Quantifying dynamic and rheological properties of suspensions of soft biological particles such as vesicles, capsules, and red blood cells (RBCs) is fundamentally important in computational biology and biomedical engineering. In this review, recent studies on dynamic and rheological behavior of soft biological cell suspensions by computer simulations are presented, considering both unbounded and confined shear flow. Furthermore, the hemodynamic and hemorheological characteristics of RBCs in diseases such as malaria and sickle cell anemia are highlighted.

Computational fluid dynamics in the microcirculation and microfluidics: what role can the lattice Boltzmann method play?

Patient-specific simulations, efficient parametric analyses, and the study of complex processes that are otherwise experimentally intractable are facilitated through the use of Computational Fluid Dynamics (CFD) to study biological flows. This review discusses various CFD methodologies that have been applied across different biological scales, from cell to organ level. Through this discussion the lattice Boltzmann method (LBM) is highlighted as an emerging technique capable of efficiently simulating fluid problems across the midrange of scales; providing a practical analytical tool compared to methods more attuned to the extremities of scale. Furthermore, the merits of the LBM are highlighted through examples of previous applications and suggestions for future research are made. The review focusses on applications in the midrange bracket, such as cell-cell interactions, the microcirculation, and microfluidic devices; wherein the inherent mesoscale nature of the LBM renders it well suited to the incorporation of fluid-structure interaction effects, molecular/particle interactions and interfacial dynamics. The review demonstrates that the LBM has the potential to become a valuable tool across a range of emerging areas in bio-CFD, such as understanding and predicting disease, designing lab-on-a-chip devices, and elucidating complex biological processes.

Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method

Nanodrug-carrier delivery in the blood stream is strongly influenced by nanoparticle (NP) dispersion. This paper presents a numerical study on NP transport and dispersion in red blood cell (RBC) suspensions under shear and channel flow conditions, utilizing an immersed boundary fluid-structure interaction model with a lattice Boltzmann fluid solver, an elastic cell membrane model and a particle motion model driven by both hydrodynamic loading and Brownian dynamics. The model can capture the multiphase features of the blood flow. Simulations were performed to obtain an empirical formula to predict NP dispersion rate for a range of shear rates and cell concentrations. NP dispersion rate predictions from the formula were then compared to observations from previous experimental and numerical studies. The proposed formula is shown to accurately predict the NP dispersion rate. The simulation results also confirm previous findings that the NP dispersion rate is strongly influenced by local disturbances in the flow due to RBC motion and deformation. The proposed formula provides an efficient method for estimating the NP dispersion rate in modeling NP transport in large-scale vascular networks without explicit RBC and NP models.

Dynamics of a single red blood cell in simple shear flow

This work describes simulations of a red blood cell (RBC) in simple shear flow, focusing on the dependence of the cell dynamics on the spontaneous curvature of the membrane. The results show that an oblate spheroidal spontaneous curvature maintains the dimple of the RBC during tank-treading dynamics as well as exhibits off-shear-plane tumbling consistent with the experimental observations of Dupire et al. [J. Dupire, M. Socol, and A. Viallat, Proc. Natl. Acad. Sci. USA 109, 20808 (2012)] and their hypothesis of an inhomogeneous spontaneous shape. As the flow strength (capillary number Ca) is increased at a particular viscosity ratio between inner and outer fluid, the dynamics undergo transitions in the following sequence: tumbling, kayaking or rolling, tilted tank-treading, oscillating-swinging, swinging, and tank-treading. The tilted tank-treading (or spinning frisbee) regime has been previously observed in experiments but not in simulations. Two distinct classes of regime are identified: a membrane reorientation regime, where the part of membrane that is at the dimple at rest moves to the rim and vice versa, is observed in motions at high Ca such as tilted tank-treading, oscillating-swinging, swinging, and tank-treading, and a nonreorientation regime, where the part of the membrane starting from the dimple stays at the dimple, is observed in motions at low Ca such as rolling, tumbling, kayaking, and flip-flopping.

Three-dimensional dynamics of oblate and prolate capsules in shear flow

We study computationally the dynamics of oblate and prolate spheroidal capsules in simple shear flow with small inertia for a range of dimensionless shear rates. The capsule is modelled as a liquid droplet enclosed by a hyperelastic membrane, and its equatorial plane is initially tilted out of the plane of shear. We find, at low shear rates, the well-accepted tumbling motion is not always stable for both oblate and prolate capsules. For an oblate capsule, the dominant stable modes for increasing dimensionless shear rate are as follows: rolling with the equatorial plane staying in the plane of shear, precessing following Jeffery's orbit [Proc. R. Soc. London A 102, 161 (1922)], and tumbling. Interestingly, the order of modes is reversed for a prolate capsule: tumbling, precessing, and rolling with increasing dimensionless shear rate. At transitional regimes, we find the stable motion of a capsule can depend on its initial titled angle, even at the same shear rate. At high dimensionless shear rates, a spheroidal capsule undergoes a complicated oscillating-swinging motion: Its major axis oscillates about the plane of shear in addition to the swinging about a mean angle with flow direction found previously, and the amplitudes of both oscillations decrease when increasing the dimensionless shear rate towards a steady tank treading motion asymptotically. We summarize the results in phase diagrams and discuss the reorientation of both oblate and prolate capsules in a wide range of dimensionless shear rates.

Elastic behavior of a red blood cell with the membrane's nonuniform natural state: equilibrium shape, motion transition under shear flow, and elongation during tank-treading motion

Direct numerical simulations of the mechanics of a single red blood cell (RBC) were performed by considering the nonuniform natural state of the elastic membrane. A RBC was modeled as an incompressible viscous fluid encapsulated by an elastic membrane. The in-plane shear and area dilatation deformations of the membrane were modeled by Skalak constitutive equation, while out-of-plane bending deformation was formulated by the spring model. The natural state of the membrane with respect to in-plane shear deformation was modeled as a sphere ([Formula: see text]), biconcave disk shape ([Formula: see text]) and their intermediate shapes ([Formula: see text]) with the nonuniformity parameter [Formula: see text], while the natural state with respect to out-of-plane bending deformation was modeled as a flat plane. According to the numerical simulations, at an experimentally measured in-plane shear modulus of [Formula: see text] and an out-of-plane bending rigidity of [Formula: see text] of the cell membrane, the following results were obtained. (i) The RBC shape at equilibrium was biconcave discoid for [Formula: see text] and cupped otherwise; (ii) the experimentally measured fluid shear stress at the transition between tumbling and tank-treading motions under shear flow was reproduced for [Formula: see text]; (iii) the elongation deformation of the RBC during tank-treading motion from the simulation was consistent with that from in vitro experiments, irrespective of the [Formula: see text] value. Based on our RBC modeling, the three phenomena (i), (ii), and (iii) were mechanically consistent for [Formula: see text]. The condition [Formula: see text] precludes a biconcave discoid shape at equilibrium (i); however, it gives appropriate fluid shear stress at the motion transition under shear flow (ii), suggesting that a combined effect of [Formula: see text] and the natural state with respect to out-of-plane bending deformation is necessary for understanding details of the RBC mechanics at equilibrium. Our numerical results demonstrate that moderate nonuniformity in a membrane's natural state with respect to in-plane shear deformation plays a key role in RBC mechanics.

Analysis of membrane tank-tread of nonspherical capsules and red blood cells

We present an analysis of membrane motion of deformable capsules and red blood cells suspended in a linear shear flow and undergoing swinging and tumbling motions using three-dimensional numerical simulations. This study is motivated by the theory of the shape-preserving cells which predicts that the direction of the membrane rotation depends on the cell orientation and reverses at every 45° inclination angle of the cell major axis with respect to the external flow direction. By considering large deformation of capsules and red blood cells, here we investigate how the shape oscillation affects the time dependence and the direction reversal of the membrane rotation. We find that the membrane tank-tread is highly time-dependent in nature and synchronized with the time-dependent deformation. The maximum and minimum of the tank-tread velocity occur at and near the minimum and maximum deformation, respectively. For the swinging capsules and red blood cells, the direction of the membrane rotation is always along the direction of the external fluid rotation; however, a direction reversal occurs during the tumbling motion in which case the membrane rotates in the direction of the external fluid rotation when the major axis is mostly in the extensional quadrant of the shear flow, and in the opposite direction when it is mostly in the compressional quadrant. Unlike the theory which predicts the direction reversal at every 45° inclination angle irrespective of the control parameters, namely, the capillary number, viscosity ratio, and asphericity, we find that the angle at which the direction reversal occurs depends on these parameters. In particular, if the tumbling motion occurs by decreasing the capillary number, the membrane rotation is in the direction of the external flow rotation in the entire extensional quadrant, but in the opposite direction in the compressional quadrant, irrespective of the specific values of the capillary number. If the tumbling motion occurs by increasing the viscosity ratio and asphericity, the angle at which the direction reversal occurs depends on the specific values of these two parameters. The spatial variation of the tank-tread velocity also is analyzed and attributed to the straining motion of the external flow.

Effect of thermal noise on vesicles and capsules in shear flow

We add thermal noise consistently to reduced models of undeformable vesicles and capsules in shear flow and derive analytically the corresponding stochastic equations of motion. We calculate the steady-state probability distribution function and construct the corresponding phase diagrams for the different dynamical regimes. For fluid vesicles, we predict that at small shear rates thermal fluctuations induce a tumbling motion for any viscosity contrast. For elastic capsules, due to thermal mixing, an intermittent regime appears in regions where deterministic models predict only pure tank treading or tumbling.

Elastic capsule deformation in general irrotational linear flows

Knowledge of the response of elastic capsules to imposed fluid flow is necessary for predicting deformation and motion of biological cells and synthetic capsules in microfluidic devices and in the microcirculation. Capsules have been studied in shear, planar extensional, and axisymmetric extensional flows. Here, the flow gradient matrix of a general irrotational linear flow is characterized by two parameters, its strain rate, defined as the maximum of the principal strain rates, and by a new term, q, the difference in the two lesser principal strain rates, scaled by the maximum principal strain rate; this characterization is valid for ellipsoids in irrotational linear flow, and it gives good results for spheres in general linear flows at low capillary numbers. We demonstrate that deformable non-spherical particles align with the principal axes of an imposed irrotational flow. Thus, it is most practical to model deformation of non-spherical particles already aligned with the flow, rather than considering each arbitrary orientation. Capsule deformation was modeled for a sphere, a prolate spheroid, and an oblate spheroid, subjected to combinations of uniaxial, biaxial, and planar extensional flows; modeling was performed using the immersed boundary method. The time response of each capsule to each flow was found, as were the steady-state deformation factor, mean strain energy, and surface area. For a given capillary number, planar flows led to more deformation than uniaxial or biaxial extensional flows. Capsule behavior in all cases was bounded by the response of capsules to uniaxial, biaxial, and planar extensional flow.

Phase diagram and breathing dynamics of a single red blood cell and a biconcave capsule in dilute shear flow

We present phase diagrams of the single red blood cell and biconcave capsule dynamics in dilute suspension using three-dimensional numerical simulations. The computational geometry replicates an in vitro linear shear flow apparatus. Our model includes all essential properties of the cell membrane, namely, the resistance against shear deformation, area dilatation, and bending, as well as the viscosity difference between the cell interior and suspending fluids. By considering a wide range of shear rate and interior-to-exterior fluid viscosity ratio, it is shown that the cell dynamics is often more complex than the well-known tank-treading, tumbling, and swinging motion and is characterized by an extreme variation of the cell shape. As a result, it is often difficult to clearly establish whether the cell is swinging or tumbling. Identifying such complex shape dynamics, termed here as "breathing" dynamics, is the focus of this article. During the breathing motion at moderate bending rigidity, the cell either completely aligns with the flow direction and the membrane folds inward, forming two cusps, or it undergoes large swinging motion while deep, craterlike dimples periodically emerge and disappear. At lower bending rigidity, the breathing motion occurs over a wider range of shear rates, and is often characterized by the emergence of a quad-concave shape. The effect of the breathing dynamics on the tank-treading-to-tumbling transition is illustrated by detailed phase diagrams which appear to be more complex and richer than those of vesicles. In a remarkable departure from the vesicle dynamics, and from the classical theory of nondeformable cells, we find that there exists a critical viscosity ratio below which the transition is independent of the viscosity ratio, and dependent on shear rate only. Further, unlike the reduced-order models, the present simulations do not predict any intermittent dynamics of the red blood cells.

Feedback control system simulator for the control of biological cells in microfluidic cross slots and integrated microfluidic systems

Control systems for lab on chip devices require careful characterisation and design for optimal performance. Traditionally, this involves either extremely computationally expensive simulations or lengthy iteration of laboratory experiments, prototype design, and manufacture. In this paper, an efficient control simulation technique, valid for typical microchannels, Computed Interpolated Flow Hydrodynamics (CIFH), is described that is over 500 times faster than conventional time integration techniques. CIFH is a hybrid approach, utilising a combination of pre-computed flows and hydrodynamic equations and allows the efficient simulation of dynamic control systems for the transport of cells through micro-fluidic devices. The speed-ups achieved by using pre-computed CFD solutions mapped to an n-dimensional control parameter space, significantly accelerate the evaluation and improvement of control strategies and chip design. Here, control strategies for a naturally unstable device geometry, the microfluidic cross-slot, have been simulated and optimal parameters have been found for proposed devices capable of trapping and sorting cells.

Oscillatory tank-treading motion of erythrocytes in shear flows

In this paper, we investigate the oscillatory dynamics of the tank-treading motion of healthy human erythrocytes in shear flows with capillary number Ca = O(1) and small to moderate viscosity ratios 0.01 ≤ λ ≤ 1.5. These conditions correspond to a wide range of surrounding medium viscosities (4-600 m Pa s) and shear flow rates (2-560 s(-1)), and match those used in ektacytometry systems. For a given viscosity ratio, as the flow rate increases, the steady-state erythrocyte length L (in the shear plane) increases logarithmically while its depth W (normal to the shear plane) decreases logarithmically. In addition, the flow rate increase dampens the oscillatory erythrocyte inclination but not its length oscillations (which show relative variations of about 5-8%). For a given flow rate, as the viscosity ratio increases, the erythrocyte length L contracts while its depth W increases (i.e., the cell becomes less deformed) with a small decrease in the length variations. The average orientation angle of the erythrocyte shows a significant decrease with the viscosity ratio as does the angle oscillation while the oscillation period increases. These trends continue in higher viscosity ratios resulting eventually in the transition from a (weakly oscillatory) tank-treading motion to a tumbling motion. Our computations show that the erythrocyte width S, which exists in the shear plane, is practically invariant in time, capillary number, and viscosity ratio, and corresponds to a real cell thickness of about 2.5 μm. Comparison of our computational results with the predictions of (low degree-of-freedom) theoretical models and experimental findings, suggests that the energy dissipation due to the shape-memory effects is more significant than the energy dissipation due to the membrane viscosity. Our work shows that the oscillatory tank-treading motion can account for more than 50% of the variations found in ektacytometry systems; thus, researchers who wish to study inherent differences between erythrocytes within a population must devise a way of monitoring individual cells over time so that they can remove the oscillation effects.

Micro-capsules in shear flow

This paper deals with flow-induced shape changes of elastic capsules. The state of the art concerning both theory and experiments is briefly reviewed starting with dynamically induced small deformation of initially spherical capsules and the formation of wrinkles on polymerized membranes. Initially non-spherical capsules show tumbling and tank-treading motion in shear flow. Theoretical descriptions of the transition between these two types of motion assuming a fixed shape are at odds with the full capsule dynamics obtained numerically. To resolve the discrepancy, we expand the exact equations of motion for small deformations and find that shape changes play a dominant role. We classify the dynamical phase transitions and obtain numerical and analytical results for the phase boundaries as a function of viscosity contrast, shear and elongational flow rate. We conclude with perspectives on time-dependent flow, on shear-induced unbinding from surfaces, on the role of thermal fluctuations and on applying the concepts of stochastic thermodynamics to these systems.

Tank-treading and tumbling frequencies of capsules and red blood cells

This study is motivated in part by the discrepancy that exists in the literature with regard to the dependence of the tank-treading frequency of red blood cells on the shear rate and suspending medium viscosity. Here we consider three-dimensional numerical simulations of deformable capsules of initially spherical and oblate spheroidal shapes and biconcave discoid representing the red blood cell resting shape. By considering a much broader range of the viscosity ratio (ratio of capsule or cell interior to suspending fluid viscosity), shear rate, and aspect ratio (ratio of minor to major axes) than that considered in the previous experiments, we find several new characteristics of the tank-treading and tumbling frequencies that have not been reported earlier. These new characteristics are the result of the large shape deformation and the coupling between shape and angular oscillations of the capsules or cells. For the spherical and oblate spheroidal capsules, the tank-treading frequency shows a nonmonotonic trend that is characterized by an initial decrease leading to a minimum followed by an increase with increasing viscosity ratio. For red blood cells, we find two regimes of the viscosity dependence of the tank-treading frequency: an exponential regime in which the tank-treading frequency decreases at a slower rate with increasing viscosity ratio, and a logarithmic range in which it decreases at a much faster rate. While this trend agrees well with different theoretical models of shape-preserving capsules, it was not evident in previous experimental results. When the shear rate dependence is considered, the tank-treading frequency of red blood cells and capsules of highly elongated initial shapes exhibits a nonmonotonic trend that is characterized by an initial increase leading to a maximum followed by a sharp decrease with decreasing shear rate. This anomalous behavior of the tank-treading frequency is shown to be due to a breathing-like dynamics of the capsule or cell that is characterized by a repeated emergence and absence of deep, crater-like dimples, and a large swinging motion. We further observe that the tumbling frequency exhibits a decreasing trend with increasing viscosity ratio that is in contrast to the theoretical result for the shape-preserving capsules and is due to the periodic deformation and preferential alignment of the capsules in the extensional quadrant of the flow.

Effect of bending stiffness on the deformation of liquid capsules enclosed by thin shells in shear flow

Shear-induced deformation of liquid capsule enclosed by thin shell causes the development of in-plane tensions and bending moments due to the shell thickness or to a preferred three-dimensional unstressed configuration. This paper considers the effect of bending stiffness due to a preferred three-dimensional structure on the deformation and motion of the liquid capsule. To perform the numerical simulations, an improved formulation for computing the forces generated on the capsule surface during deformation is proposed. This formulation takes full account of large deformation kinematics and the development of in-plane tensions and bending moments. The deformation and orientation dynamics of capsules with different reference shapes are studied under various shear rates, viscosity ratios, and bending modulus. The numerical results show that the bending stiffness not only restricts the deformation but also affects the motion mode of the capsules. In addition, raising bending stiffness amplifies the shape deformation oscillations in tank-treading mode but reduces the oscillations in tumbling mode.

Dynamic modes of red blood cells in oscillatory shear flow

The dynamics of red blood cells (RBCs) in oscillatory shear flow was studied using differential equations of three variables: a shape parameter, the inclination angle θ, and phase angle ϕ of the membrane rotation. In steady shear flow, three types of dynamics occur depending on the shear rate and viscosity ratio. (i) tank-treading (TT): ϕ rotates while the shape and θ oscillate. (ii) tumbling (TB): θ rotates while the shape and ϕ oscillate. (iii) intermediate motion: both ϕ and θ rotate synchronously or intermittently. In oscillatory shear flow, RBCs show various dynamics based on these three motions. For a low shear frequency with zero mean shear rate, a limit-cycle oscillation occurs, based on the TT or TB rotation at a high or low shear amplitude, respectively. This TT-based oscillation well explains recent experiments. In the middle shear amplitude, RBCs show an intermittent or synchronized oscillation. As shear frequency increases, the vesicle oscillation becomes delayed with respect to the shear oscillation. At a high frequency, multiple limit-cycle oscillations coexist. The thermal fluctuations can induce transitions between two orbits at very low shear amplitudes. For a high mean shear rate with small shear oscillation, the shape and θ oscillate in the TT motion but only one attractor exists even at high shear frequencies. The measurement of these oscillatory modes is a promising tool for quantifying the viscoelasticity of RBCs, synthetic capsules, and lipid vesicles.

Dynamic modes of microcapsules in steady shear flow: effects of bending and shear elasticities

The dynamics of microcapsules in steady shear flow were studied using a theoretical approach based on three variables: the Taylor deformation parameter αD , the inclination angle θ , and the phase angle ϕ of the membrane rotation. It is found that the dynamic phase diagram shows a remarkable change with an increase in the ratio of the membrane shear and bending elasticities. A fluid vesicle (no shear elasticity) exhibits three dynamic modes: (i) tank treading at low viscosity ηin of internal fluid (αD and θ relaxes to constant values), (ii) tumbling (TB) at high ηin (θ rotates), and (iii) swinging (SW) at middle ηin and high shear rates γ (θ oscillates). All of three modes are accompanied by a membrane (ϕ) rotation. For microcapsules with low shear elasticity, the TB phase with no ϕ rotation and the coexistence phase of SW and TB motions are induced by the energy barrier of ϕ rotation. Synchronization of ϕ rotation with TB rotation or SW oscillation occurs with integer ratios of rotational frequencies. At high shear elasticity, where a saddle point in the energy potential disappears, intermediate phases vanish and either ϕ or θ rotation occurs. This phase behavior agrees with recent simulation results of microcapsules with low bending elasticity.

Rheology of a dilute suspension of liquid-filled elastic capsules

Rheology of a dilute suspension of liquid-filled elastic capsules in linear shear flow is studied by three-dimensional numerical simulations using a front-tracking method. This study is motivated by a recent discovery that a suspension of viscous vesicles exhibits a shear viscosity minimum when the vesicles undergo an unsteady vacillating-breathing dynamics at the threshold of a transition between the tank-treading and tumbling motions. Here we consider capsules of spherical resting shape for which only a steady tank-treading motion is observed. A comprehensive analysis of the suspension rheology is presented over a broad range of viscosity ratio (ratio of internal-to-external fluid viscosity), shear rate (or, capillary number), and capsule surface-area dilatation. We find a result that the capsule suspension exhibits a shear viscosity minimum at moderate values of the viscosity ratio, and high capillary numbers, even when the capsules are in a steady tank-treading motion. It is further observed that the shear viscosity minimum exists for capsules with area-dilating membranes but not for those with nearly incompressible membranes. Nontrivial results are also observed for the normal stress differences which are shown to decrease with increasing capillary number at high viscosity ratios. Such nontrivial results neither can be predicted by the small-deformation theory nor can be explained by the capsule geometry alone. Physical mechanisms underlying these results are studied by decomposing the particle stress tensor into a contribution due to the elastic stresses in the capsule membrane and a contribution due to the viscosity differences between the internal and suspending fluids. It is shown that the elastic contribution is shear-thinning, but the viscous contribution is shear thickening. The coupling between the capsule geometry and the elastic and viscous contributions is analyzed to explain the observed trends in the bulk rheology.

Effect of the natural state of an elastic cellular membrane on tank-treading and tumbling motions of a single red blood cell

A two-dimensional computer simulation model was proposed for tank-treading and tumbling motions of an elastic biconcave red blood cell (RBC) under steady shear flow. The RBC model consisted of an outer cellular membrane and an inner fluid; the membrane's elastic properties were modeled by springs for stretch/compression and bending to consider the membrane's natural state in a practical manner. Membrane deformation was coupled with incompressible viscous flow of the inner and outer fluids of the RBC using a particle method. The proposed simulation model was capable of reproducing tank-treading and tumbling motions of an RBC along with rotational oscillation, which is the transition between the two motions. In simulations using the same initial RBC shape with different natural states of the RBC membrane, only tank-treading motion was exhibited in the case of a uniform natural state of the membrane, and a nonuniform natural state was necessary to generate the rotational oscillation and tumbling motion. Simulation results corresponded to published data from experimental and computational studies. In the range of simulation parameters considered, the relative membrane elastic force versus fluid viscous force was approximately 1 at the transition when the natural state nonuniformity was taken into account in estimating the membrane elastic force. A combination of natural state nonuniformity and elastic spring constant determined that change in the RBC deformation at the transition is that from a large compressive deformation to no deformation, such as rigid body.

Swinging and synchronized rotations of red blood cells in simple shear flow

The dynamics of red blood cells (RBCs) in simple shear flow was studied using a theoretical approach based on three variables: a shape parameter, the inclination angle theta, and phase angle phi of the membrane rotation. At high shear rate and low viscosity contrast of internal fluid, RBCs exhibit tank-treading motion, where phi rotates with swinging oscillation of shape and theta . At low shear rate, tumbling motion occurs and theta rotates. In the middle region between these two phases, it is found that synchronized rotation of phi and theta with integer ratios of the frequencies occurs in addition to intermittent rotation. These dynamics are robust to the modification of the potential of the RBC shape and membrane rotation. Our results agree well with recent experiments.

Elastic capsules in shear flow: analytical solutions for constant and time-dependent shear rates

We investigate the dynamics of microcapsules in linear shear flow within a reduced model with two degrees of freedom. In previous work for steady shear flow, the dynamic phases of this model, i.e. swinging, tumbling and intermittent behaviour, have been identified using numerical methods. In this paper, we integrate the equations of motion in the quasi-spherical limit analytically for time-constant and time-dependent shear flow using matched asymptotic expansions. Using this method, we find analytical expressions for the mean tumbling rate in general time-dependent shear flow. The capsule dynamics is studied in more detail when the inverse shear rate is harmonically modulated around a constant mean value for which a dynamic phase diagram is constructed. By a judicious choice of both modulation frequency and phase, tumbling motion can be induced even if the mean shear rate corresponds to the swinging regime. We derive expressions for the amplitude and width of the resonance peaks as a function of the modulation frequency.

Dynamics of nonspherical capsules in shear flow

Three-dimensional numerical simulations using a front-tracking method are presented on the dynamics of oblate shape capsules in linear shear flow by considering a broad range of viscosity contrast (ratio of internal-to-external fluid viscosity), shear rate (or capillary number), and aspect ratio. We focus specifically on the coupling between the shape deformation and orientation dynamics of capsules, and show how this coupling influences the transition from the tank-treading to tumbling motion. At low capillary numbers, three distinct modes of motion are identified: a swinging or oscillatory (OS) mode at a low viscosity contrast in which the inclination angle theta(t) oscillates but always remains positive; a vacillating-breathing (VB) mode at a moderate viscosity contrast in which theta(t) periodically becomes positive and negative, but a full tumbling does not occur; and a pure tumbling mode (TU) at a higher viscosity contrast. At higher capillary numbers, three types of transient motions occur, in addition to the OS and TU modes, during which the capsule switches from one mode to the other as (i) VB to OS, (ii) TU to VB to OS, and (iii) TU to VB. Phase diagrams showing various regimes of capsule dynamics are presented. For all modes of motion (OS, VB, and TU), a large-amplitude oscillation in capsule shape and a strong coupling between the shape deformation and orientation dynamics are observed. It is shown that the coupling between the shape deformation and orientation is the strongest in the VB mode, and hence at a moderate viscosity contrast, for which the amplitude of shape deformation reaches its maximum. The numerical results are compared with the theories of Keller and Skalak, and Skotheim and Secomb. Significant departures from the two theories are discussed and related to the strong coupling between the shape deformation, inclination, and transition dynamics.