On the effects of membrane viscosity on transient red blood cell dynamics

Droplet dynamics in homogeneous isotropic turbulence with the immersed boundary-lattice Boltzmann method

We develop a numerical method for simulating the dynamics of a droplet immersed in a generic time-dependent velocity gradient field. This approach is grounded on the hybrid coupling between the lattice Boltzmann (LB) method, employed for the flow simulation, and the immersed boundary (IB) method, utilized to couple the droplet with the surrounding fluid. We show how to enrich the numerical scheme with a mesh regularization technique, allowing droplets to sustain large deformations. The resulting methodology is adapted to simulate the dynamics of droplets in homogeneous and isotropic turbulence, with the characteristic size of the droplet being smaller than the characteristic Kolmogorov scale of the outer turbulent flow. We report statistical results for droplet deformation and orientation collected from an ensemble of turbulent trajectories, as well as comparisons with theoretical models in the limit of small deformation.

Morphology, repulsion, and ordering of red blood cells in viscoelastic flows under confinement

Red blood cells (RBC), the primary carriers of oxygen in the body, play a crucial role across several biomedical applications, while also being an essential model system of a deformable object in the microfluidics and soft matter fields. However, RBC behavior in viscoelastic liquids, which holds promise in enhancing microfluidic diagnostic applications, remains poorly studied. We here show that using viscoelastic polymer solutions as a suspending carrier causes changes in the clustering and shape of flowing RBC in microfluidic flows when compared to a standard Newtonian suspending liquid. Additionally, when the local RBC concentration increases to a point where hydrodynamic interactions take place, we observe the formation of equally-spaced RBC structures, resembling the viscoelasticity-driven ordered particles observed previously in the literature, thus providing the first experimental evidence of viscoelasticity-driven cell ordering. The observed RBC ordering, unaffected by polymer molecular architecture, persists as long as the surrounding medium exhibits shear-thinning, viscoelastic properties. Complementary numerical simulations reveal that viscoelasticity-induced repulsion between RBCs leads to equidistant structures, with shear-thinning modulating this effect. Our results open the way for the development of new biomedical technologies based on the use of viscoelastic liquids while also clarifying fundamental aspects related to multibody hydrodynamic interactions in viscoelastic microfluidic flows.

A computational study of cell membrane damage and intracellular delivery in a cross-slot microchannel

We propose a three-dimensional computational framework to simulate the flow-induced cell membrane damage and the resulting enhanced intracellular mass transport in a cross-slot microchannel. We model the cell as a liquid droplet enclosed by a viscoelastic membrane and solve the cell deformation using a well-tested immersed-boundary lattice-Boltzmann method. The cell membrane damage, which is directly related to the membrane permeability, is considered using continuum damage mechanics. The transport of the diffusive solute into the cell is solved by a lattice-Boltzmann model. After validating the computational framework against several benchmark cases, we consider a cell flowing through a cross-slot microchannel, focusing on the effects of the flow strength, channel fluid viscosity and cell membrane viscosity on the membrane damage and enhanced intracellular transport. Interestingly, we find that under a comparable pressure drop across the device, for cells with low membrane viscosity, the inertial flow regime, which can be achieved by driving a low-viscosity liquid at a high speed, often leads to much larger membrane damage, compared with the high-viscosity low-speed viscous flow regime. However, the enhancement can be significantly reduced or even reversed by an increase of the cell membrane viscosity, which limits cell deformation, particularly in the inertial flow regime. Our computational framework and simulation results may guide the design and optimisation of microfluidic devices, which use cross-slot geometry to disrupt cell membranes to enhance intracellular delivery of solutes.

Lattice Boltzmann simulation of deformable fluid-filled bodies: progress and perspectives

With the rapid development of studies involving droplet microfluidics, drug delivery, cell detection, and microparticle synthesis, among others, many scientists have invested significant efforts to model the flow of these fluid-filled bodies. Motivated by the intricate coupling between hydrodynamics and the interactions of fluid-filled bodies, several methods have been developed. The objective of this review is to present a compact foundation of the methods used in the literature in the context of lattice Boltzmann methods. For hydrodynamics, we focus on the lattice Boltzmann method due to its specific ability to treat time- and spatial-dependent boundary conditions and to incorporate new physical models in a computationally efficient way. We split the existing methods into two groups with regard to the interfacial boundary: fluid-structure and fluid-fluid methods. The fluid-structure methods are characterised by the coupling between fluid dynamics and mechanics of the flowing body, often used in applications involving membranes and similar flexible solid boundaries. We further divide fluid-structure-based methods into two subcategories, those which treat the fluid-structure boundary as a continuum medium and those that treat it as a discrete collection of individual springs and particles. Next, we discuss the fluid-fluid methods, particularly useful for the simulations of fluid-fluid interfaces. We focus on models for immiscible droplets and their interaction in a suspending fluid and describe benchmark tests to validate the models for fluid-filled bodies.

Confinement effect on the microcapillary flow and shape of red blood cells

The ability to change shape is essential for the proper functioning of red blood cells (RBCs) within the microvasculature. The shape of RBCs significantly influences blood flow and has been employed in microfluidic lab-on-a-chip devices, serving as a diagnostic biomarker for specific pathologies and enabling the assessment of RBC deformability. While external flow conditions, such as the vessel size and the flow velocity, are known to impact microscale RBC flow, our comprehensive understanding of how their shape-adapting ability is influenced by channel confinement in biomedical applications remains incomplete. This study explores the impact of various rectangular and square channels, each with different confinement and aspect ratios, on the in vitro RBC flow behavior and characteristic shapes. We demonstrate that rectangular microchannels, with a height similar to the RBC diameter in combination with a confinement ratio exceeding 0.9, are required to generate distinctive well-defined croissant and slipper-like RBC shapes. These shapes are characterized by their equilibrium positions in the channel cross section, and we observe a strong elongation of both stable shapes in response to the shear rate across the different channels. Less confined channel configurations lead to the emergence of unstable other shape types that display rich shape dynamics. Our work establishes an experimental framework to understand the influence of channel size on the single-cell flow behavior of RBCs, providing valuable insights for the design of biomicrofluidic single-cell analysis applications.

Dynamics and rheology of vesicles under confined Poiseuille flow

The rheological behavior and dynamics of a vesicle suspension, serving as a simplified model for red blood cells, are explored within a Poiseuille flow under the Stokes limit. Investigating vesicle response has led to the identification of novel solutions that complement previously documented forms like the parachute and slipper shapes. This study has brought to light the existence of alternative configurations, including a fully off-centered form and a multilobe structure. The study unveils the presence of two distinct branches associated with the slipper shape. One branch arises as a consequence of a supercritical bifurcation from the symmetric parachute shape, while the other emerges from a saddle-node bifurcation. Notably, the findings are represented through diagrams that display data collapsing harmoniously based on a combination of independent dimensionless parameters. Delving into the rheological implications, a remarkable observation emerges: the normalized viscosity (i.e. similar to intrinsic viscosity) exhibits a non-monotonic trend as a function of vesicle concentration. Initially, the normalized viscosity diminishes as the concentration increases, followed by a subsequent rise at higher concentrations. Noteworthy is the presence of a minimum value in the normalized viscosity at lower concentrations, aligning well with the concentrations observed in microcirculation scenarios. The intricate behavior of the normalized viscosity can be attributed to a delicate spatial arrangement within the suspension. Importantly, this trend echoes the observations made in a linear shear flow scenario, thereby underscoring the universality of the rheological behavior for confined suspensions.

Lift at low Reynolds number

Lift forces are widespread in hydrodynamics. These are typically observed for big and fast objects and are often associated with a combination of fluid inertia (i.e. large Reynolds numbers) and specific symmetry-breaking mechanisms. In contrast, the properties of viscosity-dominated (i.e. low Reynolds numbers) flows make it more difficult for such lift forces to emerge. However, the inclusion of boundary effects qualitatively changes this picture. Indeed, in the context of soft and biological matter, recent studies have revealed the emergence of novel lift forces generated by boundary softness, flow gradients and/or surface charges. The aim of the present review is to gather and analyse this corpus of literature, in order to identify and unify the questioning within the associated communities, and pave the way towards future research.

Numerical analysis of oxygen uptake processes by red blood cells in stopped-flow measurements: Effects of cell shape, membrane permeability and unstirred layer

The transport process of oxygen and other gas species across red blood cell (RBC) membrane is of great importance for better understanding the critical biological functions of RBCs, and the stopped-flow experiments have often been employed for such investigations. In previous stopped-flow analyses, the RBC had usually been represented by a spherical capsule based on the RBC volume, and an assumed unstirred layer (USL) thickness had been used to determine the membrane permeability. In this research, unlike these previous studies, we simulate the oxygen uptake process with different RBC shapes (shperical, ellipsoidal and biconcave) and examine the effects of USL thickness and membrane permeability over broad ranges based on literature values. Our results show that the excess membrane area can greatly improve the oxygen transport efficiency, and a same uptake half-time can be obtained using different combinations of USL thickness and membrane permeability. These findings raise concerns on the reliability and uncertainty for the results and conclusions in previous studies, and also call for more complete numerical models, for example, with the fluid flow and cell deformation considered, and more in-depth investigations on the oxygen transport processes.

Effect of Cell Age and Membrane Rigidity on Red Blood Cell Shape in Capillary Flow

Blood flow in the microcirculatory system is crucially affected by intrinsic red blood cell (RBC) properties, such as their deformability. In the smallest vessels of this network, RBCs adapt their shapes to the flow conditions. Although it is known that the age of RBCs modifies their physical properties, such as increased cytosol viscosity and altered viscoelastic membrane properties, the evolution of their shape-adapting abilities during senescence remains unclear. In this study, we investigated the effect of RBC properties on the microcapillary in vitro flow behavior and their characteristic shapes in microfluidic channels. For this, we fractioned RBCs from healthy donors according to their age. Moreover, the membranes of fresh RBCs were chemically rigidified using diamide to study the effect of isolated graded-membrane rigidity. Our results show that a fraction of stable, asymmetric, off-centered slipper-like cells at high velocities decreases with increasing age or diamide concentration. However, while old cells form an enhanced number of stable symmetric croissants at the channel centerline, this shape class is suppressed for purely rigidified cells with diamide. Our study provides further knowledge about the distinct effects of age-related changes of intrinsic cell properties on the single-cell flow behavior of RBCs in confined flows due to inter-cellular age-related cell heterogeneity.

Surface viscosity in simple liquids

The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dynamics of interfacial layers in molecular simulations of simple liquids. We estimate the surface viscosity to be 8-16 times smaller than that of the bulk fluid at the thermodynamic point considered. This result can have important implications for reactions at liquid surfaces in atmospheric chemistry and catalysis.

Effect of droplet deformability on shear thinning in a cylindrical channel

Droplets suspended in fluids flowing through microchannels are often encountered in different contexts and scales, from oil extraction down to microfluidics. They are usually flexible and deform as a product of the interplay between flexibility, hydrodynamics, and interaction with confining walls. Deformability adds distinct characteristics to the nature of the flow of these droplets. We simulate deformable droplets suspended in a fluid at a high volume fraction flowing through a cylindrical wetting channel. We find a discontinuous shear thinning transition, which depends on the droplet deformability. The capillary number is the main dimensionless parameter that controls the transition. Previous results have focused on two-dimensional configurations. Here we show that, in three dimensions, even the velocity profile is different. To perform this study, we improve and extend to three dimensions a multicomponent lattice Boltzmann method which prevents the coalescence between the droplets.

Effect of constitutive law on the erythrocyte membrane response to large strains

Three constitutive laws, that is the Skalak, neo-Hookean and Yeoh laws, commonly employed for describing the erythrocyte membrane mechanics are theoretically analyzed and numerically investigated to assess their accuracy for capturing erythrocyte deformation characteristics and morphology. Particular emphasis is given to the nonlinear deformation regime, where it is known that the discrepancies between constitutive laws are most prominent. Hence, the experiments of optical tweezers and micropipette aspiration are considered here, for which relationships between the individual shear elastic moduli of the constitutive laws can also be established through analysis of the tension-deformation relationship. All constitutive laws were found to adequately predict the axial and transverse deformations of a red blood cell subjected to stretching with optical tweezers for a constant shear elastic modulus value. As opposed to Skalak law, the neo-Hookean and Yeoh laws replicated the erythrocyte membrane folding, that has been experimentally observed, with the trade-off of sustaining significant area variations. For the micropipette aspiration, the suction pressure-aspiration length relationship could be excellently predicted for a fixed shear elastic modulus value only when Yeoh law was considered. Importantly, the neo-Hookean and Yeoh laws reproduced the membrane wrinkling at suction pressures close to those experimentally measured. None of the constitutive laws suffered from membrane area compressibility in the micropipette aspiration case.

Effects of membrane viscoelasticity on the red blood cell dynamics in a microcapillary

The mechanical properties of red blood cells (RBCs) play key roles in their biological functions in microcirculation. In particular, RBCs must deform significantly to travel through microcapillaries with sizes comparable with or even smaller than their own. Although the dynamics of RBCs in microcapillaries have received considerable attention, the effect of membrane viscoelasticity has been largely overlooked. In this work, we present a computational study based on the boundary integral method and thin-shell mechanics to examine how membrane viscoelasticity influences the dynamics of RBCs flowing through straight and constricted microcapillaries. Our results reveal that the cell with a viscoelastic membrane undergoes substantially different motion and deformation compared with results based on a purely elastic membrane model. Comparisons with experimental data also suggest the importance of accounting for membrane viscoelasticity to properly capture the transient dynamics of an RBC flowing through a microcapillary. Taken together, these findings demonstrate the significant effects of membrane viscoelasticity on RBC dynamics in different microcapillary environments. The computational framework also lays the groundwork for more accurate quantitative modeling of the mechanical response of RBCs in their mechanotransduction process in subsequent investigations.

Energy Dissipation in the Human Red Cell Membrane

The membrane of the human red cell consists of a lipid bilayer and a so-called membrane skeleton attached on the cytoplasmic side of the bilayer. Upon the deformation of red cells, energy is dissipated in their cytoplasm and their membrane. As to the membrane, three contributions can be distinguished: (i) A two-dimensional shear deformation with the membrane viscosity as the frictional parameter; (ii) A motion of the membrane skeleton relative to the bilayer; (iii) A relative motion of the two monolayers of the bilayer. The frictional parameter in contributions (ii) and (iii) is a frictional coefficient specific for the respective contribution. This perspective describes the history up to recent advances in the knowledge of these contributions. It reviews the mechanisms of energy dissipation on a molecular scale and suggests new ones, particularly for the first contribution. It proposes a parametric fitting expected to shed light on the discrepant values found for the membrane viscosity by different experimental approaches. It proposes strategies that could allow the determination of the frictional coefficients pertaining to the second and the third contribution. It highlights the consequences characteristic times have on the state of the red cell membrane in circulation as well as on the adaptation of computer models to the red cell history in an in vitro experiment.

Tank-treading dynamics of red blood cells in shear flow: On the membrane viscosity rheology

In this article, extensive three-dimensional simulations are conducted for tank-treading (TT) red blood cells (RBCs) in shear flow with different cell viscous properties and flow conditions. Apart from recent numerical studies on TT RBCs, this research considers the uncertainty in cytoplasm viscosity, covers a more complete range of shear flow situations of available experiments, and examines the TT behaviors in more details. Key TT characteristics, including the rotation frequency, deformation index, and inclination angle, are compared with available experimental results of similar shear flow conditions. Fairly good simulation-experiment agreements for these parameters can be obtained by adjusting the membrane viscosity values; however, different rheological relationships between the membrane viscosity and the flow shear rate are noted for these comparisons: shear thinning from the TT frequency, Newtonian from the inclination angle, and shear thickening from the cell deformation. Previous studies claimed a shear-thinning membrane viscosity model based on the TT frequency results; however, such a conclusion seems premature from our results and more carefully designed and better controlled investigations are required for the RBC membrane rheology. In addition, our simulation results reveal complicate RBC TT features and such information could be helpful for a better understanding of in vivo and in vitro RBC dynamics.

Red blood cell shape transitions and dynamics in time-dependent capillary flows

The dynamics of single red blood cells (RBCs) determine microvascular blood flow by adapting their shape to the flow conditions in the narrow vessels. In this study, we explore the dynamics and shape transitions of RBCs on the cellular scale under confined and unsteady flow conditions using a combination of microfluidic experiments and numerical simulations. Tracking RBCs in a comoving frame in time-dependent flows reveals that the mean transition time from the symmetric croissant to the off-centered, nonsymmetric slipper shape is significantly faster than the opposite shape transition, which exhibits pronounced cell rotations. Complementary simulations indicate that these dynamics depend on the orientation of the RBC membrane in the channel during the time-dependent flow. Moreover, we show how the tank-treading movement of slipper-shaped RBCs in combination with the narrow channel leads to oscillations of the cell's center of mass. The frequency of these oscillations depends on the cell velocity, the viscosity of the surrounding fluid, and the cytosol viscosity. These results provide a potential framework to identify and study pathological changes in RBC properties.

Lattice Boltzmann simulations on the tumbling to tank-treading transition: effects of membrane viscosity

The tumbling to tank-treading (TB-TT) transition for red blood cells (RBCs) has been widely investigated, with a main focus on the effects of the viscosity ratio [Formula: see text] (i.e., the ratio between the viscosities of the fluids inside and outside the membrane) and the shear rate [Formula: see text] applied to the RBC. However, the membrane viscosity [Formula: see text] plays a major role in a realistic description of RBC dynamics, and only a few works have systematically focused on its effects on the TB-TT transition. In this work, we provide a parametric investigation on the effect of membrane viscosity [Formula: see text] on the TB-TT transition for a single RBC. It is found that, at fixed viscosity ratios [Formula: see text], larger values of [Formula: see text] lead to an increased range of values of capillary number at which the TB-TT transition occurs; moreover, we found that increasing [Formula: see text] or increasing [Formula: see text] results in a qualitatively but not quantitatively similar behaviour. All results are obtained by means of mesoscale numerical simulations based on the lattice Boltzmann models. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Loading and relaxation dynamics of a red blood cell

We use mesoscale numerical simulations to investigate the unsteady dynamics of a single red blood cell (RBC) subjected to an external mechanical load. We carry out a detailed comparison between the loading (L) dynamics, following the imposition of the mechanical load on the RBC at rest, and the relaxation (R) dynamics, allowing the RBC to relax to its original shape after the sudden arrest of the mechanical load. Such a comparison is carried out by analyzing the characteristic times of the two corresponding dynamics, i.e., t_L and t_R. When the intensity of the mechanical load is small enough, the two kinds of dynamics are symmetrical (t_L≈t_R) and independent of the typology of mechanical load (intrinsic dynamics); otherwise, in marked contrast, an asymmetry is found, wherein the loading dynamics is typically faster than the relaxation one. This asymmetry manifests itself with non-universal characteristics, e.g., dependency on the applied load and/or on the viscoelastic properties of the RBC membrane. To deepen such a non-universal behaviour, we consider the viscosity of the erythrocyte membrane as a variable parameter and focus on three different typologies of mechanical load (mechanical stretching, shear flow, elongational flow): this allows to clarify how non-universality builds up in terms of the deformation and rotational contributions induced by the mechanical load on the membrane. Finally, we also investigate the effect of the elastic shear modulus on the characteristic times t_L and t_R. Our results provide crucial and quantitative information on the unsteady dynamics of RBC and its membrane response to the imposition/cessation of external mechanical loads.

Similar but Distinct Roles of Membrane and Interior Fluid Viscosities in Capsule Dynamics in Shear Flows

PURPOSE

The dynamics of biological capsules and red blood cells in shear flows has been studied extensively with experimental, analytical, and numerical methods. In particular, the effects of various parameters, including the shear rate or shear stress, membrane elasticity, capsule shape, and interior fluid viscosity, have been investigated carefully. The role of the membrane viscosity for capsule deformation dynamics has not been examined adequately. In previous studies, the so-called energy dissipation ratio has been used to account for the membrane viscosity effect by increasing the interior viscosity; however, the applicability and accuracy of this treatment have not been evaluated carefully.

METHODS

In this study, using the recently developed finite-difference scheme for immersed boundary simulations of viscoelastic membranes, we conduct comprehensive numerical simulations of the deformation processes of an originally spherical capsule in shear flows with various combinations of membrane and interior fluid viscosities.

RESULTS

Our results show that the membrane and interior fluid viscosity have similar however different effects on the capsule deformation dynamics. While the capsule deformation decreases with both membrane and interior fluid viscosities, a typical decrease-then-increase variation is observed for the inclination angle as the membrane viscosity increases, instead of the monotonic decrease in the inclination angle with the interior fluid viscosity increase. Also, although both large membrane and interior fluid viscosity values can introduce oscillations in the capsule deformation and inclination, larger aptitudes and slow decay processes are noticed at larger membrane viscosities. The variations of other dynamic parameters of the capsule, including the circumference, average membrane velocity, and rotation frequency, are also analyzed, and an intuitive mechanism is proposed to relate the membrane velocity and rotation frequency to the capsule deformation and inclination angle. The simple mechanism is then applied to explain the spoon-like variation patterns for membrane velocity and rotation frequency observed in our results. Furthermore, we examine the validity of the energy dissipation ratio approach based on the mathematical functional dependence.

CONCLUSIONS

Our results and analysis show that the dissipation ratio is a system and process dependent variable and it cannot be treated as a constant even for the same capsule. This research is valuable for a better understanding of the complex capsule dynamics in flows and also suggests that the membrane viscosity needs to be considered explicitly for accurate and reliable results in future studies.