Prediction of anomalous blood viscosity in confined shear flow

Reply to the 'Comment on "Dynamics and rheology of vesicles under confined Poiseuille flow"' by G. Coupier and T. Podgorski, Soft Matter, 2024, , DOI: 10.1039/D3SM01679J

In this answer, we provide our arguments in support of the possibility to observe the single file-organization of red blood cells in microvessels and the resulting unexpectedly weak increase of blood viscosity with increasing hematocrit, the physiological relevance of which was questioned in the comment. The key element is that the equivalent diameter in 3D for the maximal hematocrit corresponding to a single file of red blood cells is about 10 µm and not 20 µm, as in 2D. In addition, the viscosity contrast (ratio between the cell internal and external viscosities) value must be chosen in our 2D simulation in a such a way that the effective viscosity (a linear combination of the internal, external and membrane viscosities) be close to that of a real RBC. Taking these two facts into account, we find a reasonable agreement between our 2D viscosity simulations data and experimental data, despite the crude 2D assumption.

Comment on "Dynamics and rheology of vesicles under confined Poiseuille flow" by Z. Gou, H. Zhang, A. Nait-Ouhra, M. Abbasi, A. Farutin and C. Misbah, Soft Matter, 2023, , 9101

In a recent paper, [Gou et al., Soft Matter, 2023, 19, 9101-9114] studied numerically the viscosity of a confined suspension of vesicles flowing in a channel as a function of vesicle concentration. In order to discuss the genericity of the observed behaviour, namely a nearly constant effective viscosity at low concentrations, we complement their study by a comparison with the few existing ones in the literature. In particular, we highlight that they fail to reproduce well established results for blood viscosity in microcirculation, thereby suggesting that the conclusions regarding the optimization of cell transport and oxygenation may not apply. We conclude with a quick discussion on potential improvements regarding numerical modeling, as long as physiological relevance is sought.

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.

Flow-induced "waltzing" red blood cells: Microstructural reorganization and the corresponding rheological response

We investigate flow-induced structural organization in a dilute suspension of tumbling red blood cells (RBCs) under confined shear flow. For small Reynolds (Re = 0.1) and capillary numbers (Ca), with fully coupled hydrodynamic interaction (HI) and without interparticle adhesion, we find that HI between the biconcave discoid particles prompts the formation of layered RBC chains and synchronized rotating RBC pairs, referred here as "waltzing doublets." As the volume fraction ϕ increases, more waltzing doublets appear in RBC files. Stronger shear stress disrupts structural arrangements at higher Ca. We find that the flow-induced organization of waltzing doublets changes how the suspension viscosity varies with ϕ qualitatively. The intrinsic viscosity is particularly sensitive to microstructural rearrangement, increasing (decreasing) with ϕ at low (high) Ca that correlates with the change in the fraction of doublets. We verified flow-induced collective motion with comparison to two-cell simulations in which the cell volume fraction is controlled by varying the domain volume.

Upstream wall vortices in viscoelastic flow past a cylinder

We report a novel inertia-less, elastic flow instability for a viscoelastic, shear-thinning wormlike micellar solution flowing past a microcylinder in a channel with blockage ratio B_R = 2R/W = 0.5 and aspect ratio α = H/W ≈ 5, where R ≈ 100 μm is the cylinder radius, W is the channel width, and H is the channel height. The instability manifests upstream of the cylinder and changes form with increasing Weissenberg number over the range 0.5 ≲ Wi = Uλ/R ≲ 900, where U is the average flow velocity and λ is the terminal relaxation time of the fluid. Beyond a first critical Wi, the instability begins as a bending of the streamlines near the upstream pole of the cylinder that breaks the symmetry of the flow. Beyond a second critical Wi, small, time-steady, and approximately symmetric wall-attached vortices form upstream of the cylinder. Beyond a third critical Wi, the flow becomes time dependent and pulses with a characteristic frequency commensurate with the breakage timescale of the wormlike micelles. This is accompanied by a breaking of the symmetry of the wall-attached vortices, where one vortex becomes considerably larger than the other. Finally, beyond a fourth critical Wi, a vortex forms attached to the upstream pole of the cylinder whose length fluctuates in time. The flow is highly time dependent, and the cylinder-attached vortex and wall-attached vortices compete dynamically for space and time in the channel. Our results add to the rapidly growing understanding of viscoelastic flow instabilities in microfluidic geometries.

Microfluidics Approach to the Mechanical Properties of Red Blood Cell Membrane and Their Effect on Blood Rheology

In this article, we describe the general features of red blood cell membranes and their effect on blood flow and blood rheology. We first present a basic description of membranes and move forward to red blood cell membranes' characteristics and modeling. We later review the specific properties of red blood cells, presenting recent numerical and experimental microfluidics studies that elucidate the effect of the elastic properties of the red blood cell membrane on blood flow and hemorheology. Finally, we describe specific hemorheological pathologies directly related to the mechanical properties of red blood cells and their effect on microcirculation, reviewing microfluidic applications for the diagnosis and treatment of these diseases.

Tristability in Viscoelastic Flow Past Side-by-Side Microcylinders

Viscoelastic flows through microscale porous arrays exhibit complex path selection and switching phenomena. However, understanding this process is limited by a lack of studies linking between a single object and large arrays. Here, we report experiments on viscoelastic flow past side-by-side microcylinders with variable intercylinder gap. With increasing flow rate, a sequence of two imperfect symmetry-breaking bifurcations forces selection of either one or two of the three possible flow paths around the cylinders. Tuning the gap length through the value where the first bifurcation becomes perfect reveals regions of bistability and tristability in a dimensionless flow rate-gap length phase diagram.

Disorder Suppresses Chaos in Viscoelastic Flows

Viscoelastic flows through microstructured geometries transition from steady to time dependent and chaotic dynamics under critical flow conditions. However, the implications of geometric disorder for flow stability are unknown. We measure the onset of spatiotemporal velocity fluctuations for a viscoelastic flow through microfluidic pillar arrays, having controlled variations of geometric disorder. Introducing a small perturbation into the pillar array (∼10% of the lattice constant) delays the onset of the instability to higher flow speed, and yet larger disorders (≥25%) suppress the transition to chaos. We show that disorder introduces preferential flow paths that promote shear over extensional deformation and enhance flow stability by locally reducing polymer stretching.

Shear-induced non-monotonic viscosity dependence for model red blood cell suspensions in microvessels

The cell-free layer thickness of an aggregating red blood cell (RBC) suspension in a rectangular microchannel is investigated by hybrid fluid-particle numerical modeling. Several factors affect the suspension viscosity, cell-free layer thickness, and the cell aggregate distribution. These include the hematocrit, vessel size, red cell stiffness, aggregation interaction, and shear rate. In particular, the effect of the shear rate on the cell-free layer thickness is controversial. We found that the suspension viscosity increases along with a decrease in the cell-free layer thickness as the shear rate increases for aggregating model RBCs at low shear rates. At moderate to high shear rates, the cell-free layer thickness increases with the increasing shear rate from medium to strong shear flow for both 10% and 20% red blood cell suspensions.

ATP Release by Red Blood Cells under Flow: Model and Simulations

ATP is a major player as a signaling molecule in blood microcirculation. It is released by red blood cells (RBCs) when they are subjected to shear stresses large enough to induce a sufficient shape deformation. This prominent feature of chemical response to shear stress and RBC deformation constitutes an important link between vessel geometry, flow conditions, and the mechanical properties of RBCs, which are all contributing factors affecting the chemical signals in the process of vasomotor modulation of the precapillary vessel networks. Several in vitro experiments have reported on ATP release by RBCs due to mechanical stress. These studies have considered both intact RBCs as well as cells within which suspected pathways of ATP release have been inhibited. This has provided profound insights to help elucidate the basic governing key elements, yet how the ATP release process takes place in the (intermediate) microcirculation zone is not well understood. We propose here an analytical model of ATP release. The ATP concentration is coupled in a consistent way to RBC dynamics. The release of ATP, or the lack thereof, is assumed to depend on both the local shear stress and the shape change of the membrane. The full chemo-mechanical coupling problem is written in a lattice-Boltzmann formulation and solved numerically in different geometries (straight channels and bifurcations mimicking vessel networks) and under two kinds of imposed flows (shear and Poiseuille flows). Our model remarkably reproduces existing experimental results. It also pinpoints the major contribution of ATP release when cells traverse network bifurcations. This study may aid in further identifying the interplay between mechanical properties and chemical signaling processes involved in blood microcirculation.

Shear thinning and shear thickening of a confined suspension of vesicles

Widely regarded as an interesting model system for studying flow properties of blood, vesicles are closed membranes of phospholipids that mimic the cytoplasmic membranes of red blood cells. In this study we analyze the rheology of a suspension of vesicles in a confined geometry: the suspension, bound by two planar rigid walls on each side, is subject to a shear flow. Flow properties are then analyzed as a function of shear rate γ[over ̇], the concentration of the suspension ϕ, and the viscosity contrast λ=η_{in}/η_{out}, where η_{in} and η_{out} are the fluid viscosities of the inner and outer fluids, respectively. We find that the apparent (or effective viscosity) of the suspension exhibits both shear thinning (decreasing viscosity with shear rate) or shear thickening (increasing viscosity with shear rate) in the same concentration range. The shear thinning or thickening behaviors appear as subtle phenomena, dependant on viscosity contrast λ. We provide physical arguments on the origins of these behaviors.

Particles on Demand for Kinetic Theory

A novel formulation of fluid dynamics as a kinetic theory with tailored, on-demand constructed particles removes restrictions on flow speed and temperature as compared to its predecessors, the lattice Boltzmann methods and their modifications. In the new kinetic theory, discrete particles are determined by a rigorous limit process which avoids ad hoc assumptions about their velocities. Classical benchmarks for incompressible and compressible flows demonstrate that the proposed discrete-particles kinetic theory opens up an unprecedented wide domain of applications for computational fluid dynamics.

Numerical-experimental observation of shape bistability of red blood cells flowing in a microchannel

Red blood cells flowing through capillaries assume a wide variety of different shapes owing to their high deformability. Predicting the realized shapes is a complex field as they are determined by the intricate interplay between the flow conditions and the membrane mechanics. In this work we construct the shape phase diagram of a single red blood cell with a physiological viscosity ratio flowing in a microchannel. We use both experimental in vitro measurements as well as 3D numerical simulations to complement the respective other one. Numerically, we have easy control over the initial starting configuration and natural access to the full 3D shape. With this information we obtain the phase diagram as a function of initial position, starting shape and cell velocity. Experimentally, we measure the occurrence frequency of the different shapes as a function of the cell velocity to construct the experimental diagram which is in good agreement with the numerical observations. Two different major shapes are found, namely croissants and slippers. Notably, both shapes show coexistence at low (<1 mm s^-1) and high velocities (>3 mm s^-1) while in-between only croissants are stable. This pronounced bistability indicates that RBC shapes are not only determined by system parameters such as flow velocity or channel size, but also strongly depend on the initial conditions.

Theory and algorithms to compute Helfrich bending forces: a review

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

Front microrheology of the non-Newtonian behaviour of blood: scaling theory of erythrocyte aggregation by aging

We introduce a new framework to study the non-Newtonian behaviour of fluids at the microscale based on the analysis of front advancement. We apply this methodology to study the non-linear rheology of blood in microchannels. We carry out experiments in which the non-linear viscosity of blood samples is quantified at different haematocrits and ages. Under these conditions, blood exhibits a power-law dependence on the shear rate. In order to analyse our experimental data, we put forward a scaling theory which allows us to define an adhesion scaling number. This theory yields a scaling behaviour of the viscosity expressed as a function of the adhesion capillary number. By applying this scaling theory to samples of different ages, we are able to quantify how the characteristic adhesion energy varies as time progresses. This connection between microscopic and mesoscopic properties allows us to estimate quantitatively the change in the cell-cell adhesion energies as the sample ages.

Clusters of red blood cells in microcapillary flow: hydrodynamic versus macromolecule induced interaction

We present experiments on RBCs that flow through micro-capillaries under physiological conditions. The strong flow-shape coupling of these deformable objects leads to a rich variety of cluster formation. We show that the RBC clusters form as a subtle imbrication between hydrodynamic interactions and adhesion forces because of plasma proteins, mimicked by the polymer dextran. Clusters form along the capillaries and macromolecule-induced adhesion contributes to their stability. However, at high yet physiological flow velocities, shear stresses overcome part of the adhesion forces, and cluster stabilization due to hydrodynamics becomes stronger. For the case of pure hydrodynamic interaction, cell-to-cell distances have a pronounced bimodal distribution. Our 2D-numerical simulations on vesicles capture the transition between adhesive and non-adhesive clusters at different flow velocities.

Lattice Kinetic Theory in a Comoving Galilean Reference Frame

We prove that the fully discrete lattice Boltzmann method is invariant with respect to Galilean transformation. Based on this finding, a novel class of shifted lattices is proposed which dramatically increases the operating range of lattice Boltzmann simulations, in particular, for gas dynamics applications. A simulation of vortex-shock interaction is used to demonstrate the accuracy and efficiency of the proposed lattices. With one single algorithm it is now possible to simulate a broad range of applications, from low Mach number flows to transonic and supersonic flow regimes.

Rheology of Confined Non-Brownian Suspensions

We study the rheology of confined suspensions of neutrally buoyant rigid monodisperse spheres in plane-Couette flow using direct numerical simulations. We find that if the width of the channel is a (small) integer multiple of the sphere diameter, the spheres self-organize into two-dimensional layers that slide on each other and the effective viscosity of the suspension is significantly reduced. Each two-dimensional layer is found to be structurally liquidlike but its dynamics is frozen in time.

Viscometry of single nanoliter-volume droplets using dynamic force spectroscopy

We present an atomic force microscope-based platform for viscometry of ‘nanoliter' volume fluids.

Viscoelastic transient of confined red blood cells

The unique ability of a red blood cell to flow through extremely small microcapillaries depends on the viscoelastic properties of its membrane. Here, we study in vitro the response time upon flow startup exhibited by red blood cells confined into microchannels. We show that the characteristic transient time depends on the imposed flow strength, and that such a dependence gives access to both the effective viscosity and the elastic modulus controlling the temporal response of red cells. A simple theoretical analysis of our experimental data, validated by numerical simulations, further allows us to compute an estimate for the two-dimensional membrane viscosity of red blood cells, η(mem)(2D) ∼ 10(-7) N ⋅ s ⋅ m(-1). By comparing our results with those from previous studies, we discuss and clarify the origin of the discrepancies found in the literature regarding the determination of η(mem)(2D), and reconcile seemingly conflicting conclusions from previous works.

Vesicle dynamics in a confined Poiseuille flow: from steady state to chaos

Red blood cells (RBCs) are the major component of blood, and the flow of blood is dictated by that of RBCs. We employ vesicles, which consist of closed bilayer membranes enclosing a fluid, as a model system to study the behavior of RBCs under a confined Poiseuille flow. We extensively explore two main parameters: (i) the degree of confinement of vesicles within the channel and (ii) the flow strength. Rich and complex dynamics for vesicles are revealed, ranging from steady-state shapes (in the form of parachute and slipper shapes) to chaotic dynamics of shape. Chaos occurs through a cascade of multiple periodic oscillations of the vesicle shape. We summarize our results in a phase diagram in the parameter plane (degree of confinement and flow strength). This finding highlights the level of complexity of a flowing vesicle in the small Reynolds number where the flow is laminar in the absence of vesicles and can be rendered turbulent due to elasticity of vesicles.