A review of numerical methods for red blood cell flow simulation

Colloid Transport in Bicontinuous Nanoporous Media

Colloid transport and retention in porous media are critical processes influencing various Earth science applications, from groundwater remediation to enhanced oil recovery. These phenomena become particularly complex in the confined spaces of nanoporous media, where strong boundary layer effects and nanoconfinement significantly alter colloid behavior. In this work, we use particle dynamics models to simulate colloid transport and retention processes in bicontinuous nanoporous (BNP) media under pressure gradients. By utilizing particle-based models, we track the movement of each colloid and elucidate the underlying colloid retention mechanisms. Under unfavorable attachment conditions, the results reveal two colloid retention mechanisms: physical straining and trapping in low-flow zone. Furthermore, we investigate the effects of critical factors including colloid volume fraction, d, pressure difference, ΔP, interaction between colloids and BNP media, Ec-p, and among colloids, Ec-c, on colloid transport. Analysis of breakthrough curves and colloid displacements demonstrates that higher values of d, lower values of ΔP, and strong Ec-p attractions significantly increase colloid retention, which further lead to colloid clogging and jamming. In contrast, Ec-c has minimal impact on colloid transport due to the limited colloid-colloid interaction in nanoporous channels. This work provides critical insights into the fundamental factors governing colloid transport and retention within stochastic nanoporous materials.

Accelerated simulation methodologies for computational vascular flow modelling

Vascular flow modelling can improve our understanding of vascular pathologies and aid in developing safe and effective medical devices. Vascular flow models typically involve solving the nonlinear Navier-Stokes equations in complex anatomies and using physiological boundary conditions, often presenting a multi-physics and multi-scale computational problem to be solved. This leads to highly complex and expensive models that require excessive computational time. This review explores accelerated simulation methodologies, specifically focusing on computational vascular flow modelling. We review reduced order modelling (ROM) techniques like zero-/one-dimensional and modal decomposition-based ROMs and machine learning (ML) methods including ML-augmented ROMs, ML-based ROMs and physics-informed ML models. We discuss the applicability of each method to vascular flow acceleration and the effectiveness of the method in addressing domain-specific challenges. When available, we provide statistics on accuracy and speed-up factors for various applications related to vascular flow simulation acceleration. Our findings indicate that each type of model has strengths and limitations depending on the context. To accelerate real-world vascular flow problems, we propose future research on developing multi-scale acceleration methods capable of handling the significant geometric variability inherent to such problems.

Diseased Erythrocyte Enrichment Based on I-Shaped Pillar DLD Arrays

Enrichment of erythrocytes is a necessary step in the diagnosis of blood diseases. Due to the high deformability and viscosity of erythrocytes, they cannot be regarded as stable point-like solids, so the influence of their deformability on fluid dynamics must be considered. Therefore, by using the special effect of an I-shaped pillar (I-pillar) on erythrocytes, erythrocytes with different deformability can be made to produce different provisional distances in the chip, so as to achieve the separation of the two kinds of erythrocytes. In this study, a microfluidic chip was designed to conduct a control test between erythrocytes stored for a long time and fresh erythrocytes. At a specific flow rate, the different deformable erythrocytes in the chip move in different paths. Then, the influence of erythrocyte deformability on its movement trajectory was analyzed by two-dimensional finite element flow simulation. DLD sorting technology provides a new method for the sorting and enrichment of diseased erythrocytes.

Numerical simulation of flow characteristics in a permeable liver sinusoid with leukocytes

Double-layered channels of sinusoid lumen and Disse space separated by fenestrated liver sinusoidal endothelial cells (LSECs) endow the unique mechanical environment of the liver sinusoid network, which further guarantees its biological function. It is also known that this mechanical environment changes dramatically under liver fibrosis and cirrhosis, including the reduced plasma penetration and metabolite exchange between the two flow channels and the reduced Disse space deformability. The squeezing of leukocytes through narrow sinusoid lumen also affects the mechanical environment of liver sinusoid. To date, the detailed flow-field profile of liver sinusoid is still far from clear due to experimental limitations. It also remains elusive whether and how the varied physical properties of the pathological liver sinusoid regulate the fluid flow characteristics. Here a numerical model based on the immersed boundary method was established, and the effects of Disse space and leukocyte elasticities, endothelium permeability, and sinusoidal stenosis degree on fluid flow as well as leukocyte trafficking were specified upon a mimic liver sinusoid structure. Results showed that endothelium permeability dominantly controlled the plasma penetration velocity across the endothelium, whereas leukocyte squeezing promoted local penetration and significantly regulated wall shear stress on hepatocytes, which was strongly related to the Disse space and leukocyte deformability. Permeability and elasticity cooperatively regulated the process of leukocytes trafficking through the liver sinusoid, especially for stiffer leukocytes. This study will offer new insights into deeper understanding of the elaborate mechanical features of liver sinusoid and corresponding biological function.

Blood flow and diameter effect in the navigation process of magnetic nanocarriers inside the carotid artery

BACKGROUND AND OBJECTIVE

Serious side effects are occurred during the cancer therapy. Magnetic driving of nanoparticles is a novel method for the elimination of these effects by supplying with anticancer drug or increase the temperature of the infected area. For this reason, a numerical model for optimal guidance of nanoparticles, through the gradient magnetic field, inside the human artery system is presented in this study.

METHODS

The present method couples Computational Fluid Dynamics (CFD) and Discrete Element Method (DEM) techniques. In addition, the optimum magnetic intensity each time is evaluated by using the covariance matrix adaptation evolution strategy (CMA-ES). Under five feature blood flow velocities in cardiac cycle, the developed method evaluate and select the optimum gradient magnetic field in order to eliminate the deviation of the guided nanoparticles from a pre-described trajectory.

RESULTS

Results of the simulations indicate both the influence of the blood flow and the volume of nanocarriers in the magnetic driving process in real conditions. Specifically, the blood flow and the volume of particles are inversely proportional parameters in the magnetic navigation process. As the blood flow is decreased, the deviation of nanoparticles compared to the desired path is minimized. On the contrary, the decrease of the volume of nanocarriers increase the distance of particles from the described trajectory. However, greater magnetic gradient values are needed as the blood flow is increased. Furthermore, the imposed gradient magnetic values are strongly connected with the position of the nanoparticles and the blood blow velocity.

CONCLUSIONS

Based on the results of the present study, the most important parameter in the navigation process is the magnetic volume of particles. Under real conditions, the effect of the blood flow is insignificant compared to the volume of particles in the navigation process. In addition, great differences in the optimized magnetic sequence are presented both among the different blood flows and the volume of particles.

A new membrane formulation for modelling the flow of stomatocyte, discocyte, and echinocyte red blood cells

In this work, a numerical model that enables simulation of the deformation and flow behaviour of differently aged Red Blood Cells (RBCs) is developed. Such cells change shape and decrease in deformability as they age, thus impacting their ability to pass through the narrow capillaries in the body. While the body filters unviable cells from the blood naturally, cell aging poses key challenges for blood stored for transfusions. Therefore, understanding the influence RBC morphology and deformability have on their flow is vital. While several existing models represent young Discocyte RBC shapes well, a limited number of numerical models are developed to model aged RBC morphologies like Stomatocytes and Echinocytes. The existing models are also limited to shear and stretching simulations. Flow characteristics of these morphologies are yet to be investigated. This paper aims to develop a new membrane formulation for the numerical modelling of Stomatocyte, Discocytes and Echinocyte RBC morphologies to investigate their deformation and flow behaviour. The model used represents blood plasma using the Lattice Boltzmann Method (LBM) and the RBC membrane using the discrete element method (DEM). The membrane and the plasma are coupled by the Immersed Boundary Method (IBM). Previous LBM-IBM-DEM formulations represent RBC membrane response based on forces generated from changes in the local area, local length, local bending, and cell volume. In this new model, two new force terms are added: the local area difference force and the local curvature force, which are specially incorporated to model the flow and deformation behaviour of Stomatocytes and Echinocytes. To verify the developed model, the deformation behaviour of the three types of RBC morphologies are compared to well-characterised stretching and shear experiments. The flow modelling capabilities of the method are then demonstrated by modelling the flow of each cell through a narrow capillary. The developed model is found to be as accurate as benchmark Smoothed Particle Hydrodynamics (SPH) approaches while being significantly more computationally efficient.

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.

Local hemodynamic analysis after coronary stent implantation based on Euler-Lagrange method

Coronary stents are deployed to treat the coronary artery disease (CAD) by reopening stenotic regions in arteries to restore blood flow, but the risk of the in-stent restenosis (ISR) is high after stent implantation. One of the reasons is that stent implantation induces changes in local hemodynamic environment, so it is of vital importance to study the blood flow in stented arteries. Based on regarding the red blood cell (RBC) as a rigid solid particle and regarding the blood (including RBCs and plasma) as particle suspensions, a non-Newtonian particle suspensions model is proposed to simulate the realistic blood flow in this work. It considers the blood's flow pattern and non-Newtonian characteristic, the blood cell-cell interactions, and the additional effects owing to the bi-concave shape and rotation of the RBC. Then, it is compared with other four common hemodynamic models (Newtonian single-phase flow model, Newtonian Eulerian two-phase flow model, non-Newtonian single-phase flow model, non-Newtonian Eulerian two-phase flow model), and the comparison results indicate that the models with the non-Newtonian characteristic are more suitable to describe the realistic blood flow. Afterwards, based on the non-Newtonian particle suspensions model, the local hemodynamic environment in stented arteries is investigated. The result shows that the stent strut protrusion into the flow stream would be likely to produce the flow stagnation zone. And the stent implantation can make the pressure gradient distribution uneven. Besides, the wall shear stress (WSS) of the region adjacent to every stent strut is lower than 0.5 Pa, and along the flow direction, the low-WSS zone near the strut behind is larger than that near the front strut. What's more, in the regions near the struts in the proximal of the stent, the RBC particle stagnation zone is easy to be formed, and the erosion and deposition of RBCs are prone to occur. These hemodynamic analyses illustrate that the risk of ISR is high in the regions adjacent to the struts in the proximal and the distal ends of the stent when compared with struts in other positions of the stent. So the research can provide a suggestion on the stent design, which indicates that the strut structure in these positions of a stent should be optimized further.

Effect of Micropolar Fluid Properties on the Blood Flow in a Human Carotid Model

Blood is a non-homogeneous fluid that flows inside the human artery system and provides the cells with nutrients. In this study the auto rotation effect of blood’s microstructure on its flow inside a human carotid model is studied by using a micropolar fluid model. The study aims to investigate the flow differences that occur due to its microstructure as compared to a Newtonian fluid. We focus on the vortex viscosity effect, i.e., the ratio of microrotation viscosity to the total one, because this is the only parameter that affects directly the fluid flow. Simulations in a range of vortex viscosities, are carried out in a 3D human carotid model that is computationally reconstructed. All of the simulations are conducted at the diastolic Reynolds number that occurs in the human carotid. Results indicate that micropolarity affects blood velocity in the range of parameters studied by 4%. As micropolarity is increased, higher velocities in the center of vessels and lower near the boundaries are found as compared to a Newtonian fluid consideration. This is an indication that the increase of the fluid’s micropolarity leads to an increase of the boundary layer thickness. More importantly, an increase in vortex viscosity and the resulting increase in microrotation result in decreased shear stress in the carotid’s walls; this finding can be significant in regards to the onset and the development of atherosclerosis. Finally, the flow distribution at the carotid seems to heavily be affected by the geometry and the micropolarity of the fluid.

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

Computational Fluid Dynamics (CFD) is currently used to design and improve the hydraulic properties of biomedical devices, wherein the large scale blood circulation needs to be simulated by accounting for the mechanical response of red blood cells (RBCs) at the mesoscale. In many practical instances, biomedical devices work on time-scales comparable to the intrinsic relaxation time of RBCs: thus, a systematic understanding of the time-dependent response of erythrocyte membranes is crucial for the effective design of such devices. So far, this information has been deduced from experimental data, which do not necessarily adapt to the broad variety of fluid dynamic conditions that can be encountered in practice. This work explores the novel possibility of studying the time-dependent response of an erythrocyte membrane to external mechanical loads via mesoscale numerical simulations, with a primary focus on the detailed characterisation of the RBC relaxation time tc following the arrest of the external mechanical load. The adopted mesoscale model exploits a hybrid Immersed Boundary-Lattice Boltzmann Method (IB-LBM), coupled with the Standard Linear Solid (SLS) model to account for the RBC membrane viscosity. We underscore the key importance of the 2D membrane viscosity μm to correctly reproduce the relaxation time of the RBC membrane. A detailed assessment of the dependencies on the typology and strength of the applied mechanical loads is also provided. Overall, our findings open interesting future perspectives for the study of the non-linear response of RBCs immersed in time-dependent strain fields.

Modelling of Red Blood Cell Morphological and Deformability Changes during In-Vitro Storage

Storage lesion is a critical issue facing transfusion treatments, and it adversely affects the quality and viability of stored red blood cells (RBCs). RBC deformability is a key indicator of cell health. Deformability measurements of each RBC unit are a key challenge in transfusion medicine research and clinical haematology. In this paper, a numerical study, inspired from the previous research for RBC deformability and morphology predictions, is conducted for the first time, to investigate the deformability and morphology characteristics of RBCs undergoing storage lesion. This study investigates the evolution of the cell shape factor, elongation index and membrane spicule details, where applicable, of discocyte, echinocyte I, echinocyte II, echinocyte III and sphero-echinocyte morphologies during 42 days of in-vitro storage at 4 °C in saline-adenine-glucose-mannitol (SAGM). Computer simulations were performed to investigate the influence of storage lesion-induced membrane structural defects on cell deformability and its recoverability during optical tweezers stretching deformations. The predicted morphology and deformability indicate decreasing quality and viability of stored RBCs undergoing storage lesion. The loss of membrane structural integrity due to the storage lesion further degrades the cell deformability and recoverability during mechanical deformations. This numerical approach provides a potential framework to study the RBC deformation characteristics under varying pathophysiological conditions for better diagnostics and treatments.

Deformation behaviour of stomatocyte, discocyte and echinocyte red blood cell morphologies during optical tweezers stretching

The red blood cell (RBC) deformability is a critical aspect, and assessing the cell deformation characteristics is essential for better diagnostics of healthy and deteriorating RBCs. There is a need to explore the connection between the cell deformation characteristics, cell morphology, disease states, storage lesion and cell shape-transformation conditions for better diagnostics and treatments. A numerical approach inspired from the previous research for RBC morphology predictions and for analysis of RBC deformations is proposed for the first time, to investigate the deformation characteristics of different RBC morphologies. The present study investigates the deformability characteristics of stomatocyte, discocyte and echinocyte morphologies during optical tweezers stretching and provides the opportunity to study the combined contribution of cytoskeletal spectrin network and the lipid-bilayer during RBC deformation. The proposed numerical approach predicts agreeable deformation characteristics of the healthy discocyte with the analogous experimental observations and is extended to further investigate the deformation characteristics of stomatocyte and echinocyte morphologies. In particular, the computer simulations are performed to investigate the influence of direct stretching forces on different equilibrium cell morphologies on cell spectrin link extensions and cell elongation index, along with a parametric analysis on membrane shear modulus, spectrin link extensibility, bending modulus and RBC membrane-bead contact diameter. The results agree with the experimentally observed stiffer nature of stomatocyte and echinocyte with respect to a healthy discocyte at experimentally determined membrane characteristics and suggest the preservation of relevant morphological characteristics, changes in spectrin link densities and the primary contribution of cytoskeletal spectrin network on deformation behaviour of stomatocyte, discocyte and echinocyte morphologies during optical tweezers stretching deformation. The numerical approach presented here forms the foundation for investigations into deformation characteristics and recoverability of RBCs undergoing storage lesion.

A finite difference method with subsampling for immersed boundary simulations of the capsule dynamics with viscoelastic membranes

The membrane or interfacial viscosity is an important property in many multiphase and biofluidic situations, such as the red blood cell dynamics and emulsion stability. The immersed boundary method (IBM), which incorporates the dynamic flow-membrane interaction via force distribution and velocity interpolation, has been extensively employed in simulations of such systems. Unfortunately, direct implementation of membrane viscosity in IBM suffers severe numerical instability, which causes an IBM calculation to break down before generating any useful results. Few attempts have been recently reported; however, several concerns exist in these attempts, such as the inconsistency to the classical definition of membrane viscosity, the inability to model the shear and dilatational viscosities separately, the unjustified mathematical formulations, and the complicated algorithms and computation. To overcome these concerns, in this paper, we propose a finite difference approach for implementing membrane viscosity in immersed boundary simulations. The viscous stress is obtained via finite difference approximations to the differential strain-stress relationship, with the help of a subsampling scheme to reduce the numerical noise in the calculated strain rates. This simple method has also avoided the complicated matrix calculations in previous attempts, and hence, a better computational efficiency is expected. Detailed mathematical description of the method and key steps for its implementation in immersed boundary programs are provided. Validation and illustration calculations are performed, and our results are compared with analytical solutions and previous publications with satisfactory agreement. The influences of membrane mesh resolution and simulation time step are also examined; and the results show no indication that our finite difference method has downgraded the general IBM accuracy. Based on these simulations and analysis, we believe that our method would be a better choice for future IBM simulations of capsule dynamics with viscoelastic membranes.

A coarse-grained red blood cell membrane model to study stomatocyte-discocyte-echinocyte morphologies

An improved red blood cell (RBC) membrane model is developed based on the bilayer coupling model (BCM) to accurately predict the complete sequence of stomatocyte-discocyte-echinocyte (SDE) transformation of a RBC. The coarse-grained (CG)-RBC membrane model is proposed to predict the minimum energy configuration of the RBC from the competition between lipid-bilayer bending resistance and cytoskeletal shear resistance under given reference constraints. In addition to the conventional membrane surface area, cell volume and bilayer-leaflet-area-difference constraints, a new constraint: total-membrane-curvature is proposed in the model to better predict RBC shapes in agreement with experimental observations. A quantitative evaluation of several cellular measurements including length, thickness and shape factor, is performed for the first time, between CG-RBC model predicted and three-dimensional (3D) confocal microscopy imaging generated RBC shapes at equivalent reference constraints. The validated CG-RBC membrane model is then employed to investigate the effect of reduced cell volume and elastic length scale on SDE transformation, to evaluate the RBC deformability during SDE transformation, and to identify the most probable RBC cytoskeletal reference state. The CG-RBC membrane model can predict the SDE shape behaviour under diverse shape-transforming scenarios, in-vitro RBC storage, microvascular circulation and flow through microfluidic devices.

A monolithic fluid-structure interaction framework applied to red blood cells

A parallel fully coupled (monolithic) fluid-structure interaction (FSI) algorithm has been applied to the deformation of red blood cells (RBCs) in capillaries, where cell deformability has significant effects on blood rheology. In the present FSI algorithm, fluid domain is discretized using the side-centered unstructured finite volume method based on the Arbitrary Lagrangian-Eulerian (ALE) formulation; meanwhile, solid domain is discretized with the classical Galerkin finite element formulation for the Saint Venant-Kirchhoff material in a Lagrangian frame. In addition, the compatible kinematic boundary condition is enforced at the fluid-solid interface in order to conserve the mass of cytoplasmic fluid within the red cell at machine precision. In order to solve the resulting large-scale algebraic linear systems in a fully coupled manner, a new matrix factorization is introduced similar to that of the projection method, and the parallel algebraic multigrid solver BoomerAMG is used for the scaled discrete Laplacian provided by the HYPRE library, which we access through the PETSc library. Three important physical parameters for the blood flow are simulated and analyzed: (1) the effect of capillary diameter, (2) the effect of red cell membrane thickness, and (3) the effect of red cell spacing (hematocrit). The numerical calculations initially indicate a shape deformation in which biconcave discoid shape changes to a parachute-like shape. Furthermore, the parachute-like cell shape in small capillaries undergoes a cupcake-shaped buckling instability, which has not been observed in the literature. The instability forms thin riblike features, and the red cell deformation is not axisymmetric but three-dimensional.

Modeling Cell Adhesion and Extravasation in Microvascular System

The blood flow behaviors in the microvessels determine the transport modes and further affect the metastasis of circulating tumor cells (CTCs). Much biochemical and biological efforts have been made on CTC metastasis; however, precise experimental measurement and accurate theoretical prediction on its mechanical mechanism are limited. To complement these, numerical modeling of a CTC extravasation from the blood circulation, including the steps of adhesion and transmigration, is discussed in this chapter. The results demonstrate that CTCs prefer to adhere at the positive curvature of curved microvessels, which is attributed to the positive wall shear stress/gradient. Then, the effects of particulate nature of blood on CTC adhesion are investigated and are found to be significant in the microvessels. Furthermore, the presence of red blood cell (RBC) aggregates is also found to promote the CTC adhesion by providing an additional wall-directed force. Finally, a single cell passing through a narrow slit, mimicking CTC transmigration, was examined under the effects of cell deformability. It showed that the cell shape and surface area increase play a more important role than the cell elasticity in cell transit across the narrow slit.

The dynamics of a healthy and infected red blood cell in flow through constricted channels: A DPD simulation

Understanding the dynamics of red blood cell (RBC) motion under in silico conditions is central to the development of cost-effective diagnostic tools. Specifically, unraveling the relationship between the rheological properties and the nature of shape change in the RBC (healthy or infected) can be extremely useful. In case of malarial infection, RBC progressively loses its deformability and tends to occlude the microvessel. In the present study, detailed mesoscopic simulations are performed to investigate the deformation dynamics of an RBC in flow through a constricted channel. Specifically, the manifestation of viscous forces (through flow rates) on the passage and blockage characteristics of a healthy red blood cell (hRBC) vis-á-vis an infected red blood cell (iRBC) are investigated. A finite-sized dissipative particle dynamics framework is used to model plasma in conjunction with a discrete model for the RBC. Instantaneous wall boundary method was used to model no-slip wall boundary conditions with a good control on the near-wall density fluctuations and compressibility effects. To investigate the microvascular occlusion, the RBC motion through 2 types of constricted channels, viz, (1) a tapered microchannel and (2) a stenosed-type microchannel, were simulated. It was observed that the deformation of an infected cell was much less compared with a healthy cell, with an attendant increase in the passage time. Apart from the qualitative features, deformation indices were obtained. The deformation of hRBC was sudden, while the iRBC deformed slowly as it traversed through the constriction. For higher flow rates, both hRBC and iRBC were found to undergo severe deformation. Even under low flow rates, hRBC could easily traverse past the constricted channel. However, for sufficiently slow flow rates (eg, capillary flows), the microchannel was found to be completely blocked by the iRBC.

Structural modelling of the cardiovascular system

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

Investigation of red blood cell mechanical properties using AFM indentation and coarse-grained particle method

Background

Red blood cells (RBCs) deform significantly and repeatedly when passing through narrow capillaries and delivering dioxygen throughout the body. Deformability of RBCs is a key characteristic, largely governed by the mechanical properties of the cell membrane. This study investigated RBC mechanical properties using atomic force microscopy (AFM) with the aim to develop a coarse-grained particle method model to study for the first time RBC indentation in both 2D and 3D. This new model has the potential to be applied to further investigate the local deformability of RBCs, with accurate control over adhesion, probe geometry and position of applied force.

Results

The model considers the linear stretch capacity of the cytoskeleton, bending resistance and areal incompressibility of the bilayer, and volumetric incompressibility of the internal fluid. The model’s performance was validated against force–deformation experiments performed on RBCs under spherical AFM indentation. The model was then used to investigate the mechanisms which absorbed energy through the indentation stroke, and the impact of varying stiffness coefficients on the measured deformability. This study found the membrane’s bending stiffness was most influential in controlling RBC physical behaviour for indentations of up to 200 nm.

Conclusions

As the bilayer provides bending resistance, this infers that structural changes within the bilayer are responsible for the deformability changes experienced by deteriorating RBCs. The numerical model presented here established a foundation for future investigations into changes within the membrane that cause differences in stiffness between healthy and deteriorating RBCs, which have already been measured experimentally with AFM.

Electronic supplementary material

The online version of this article (10.1186/s12938-017-0429-5) contains supplementary material, which is available to authorized users.

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.

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.

Modeling microcirculatory blood flow: current state and future perspectives

Microvascular blood flow determines a number of important physiological processes of an organism in health and disease. Therefore, a detailed understanding of microvascular blood flow would significantly advance biophysical and biomedical research and its applications. Current developments in modeling of microcirculatory blood flow already allow to go beyond available experimental measurements and have a large potential to elucidate blood flow behavior in normal and diseased microvascular networks. There exist detailed models of blood flow on a single cell level as well as simplified models of the flow through microcirculatory networks, which are reviewed and discussed here. The combination of these models provides promising prospects for better understanding of blood flow behavior and transport properties locally as well as globally within large microvascular networks.

Understanding particle margination in blood flow - A step toward optimized drug delivery systems

Targeted delivery of drugs and imaging agents is very promising to develop new strategies for the treatment of various diseases such as cancer. For an efficient targeted adhesion, the particles have to migrate toward the walls in blood flow - a process referred to as margination. Due to a huge diversity of available carriers, a good understanding of their margination properties in blood flow depending on various flow conditions and particle properties is required. We employ a particle-based mesoscopic hydrodynamic simulation approach to investigate the margination of different carriers for a wide range of hematocrits (volume fraction of red blood cells) and flow rates. Our results show that margination strongly depends on the thickness of the available free space close to the wall, the so-called red blood cell-free layer (RBC-FL), in comparison to the carrier size. The carriers with a few micrometers in size are comparable with the RBC-FL thickness and marginate better than their sub-micrometer counterparts. Deformable carriers, in general, show worse margination properties than rigid particles. Particle margination is also found to be most pronounced in small channels with a characteristic size comparable to blood capillaries. Finally, different margination mechanisms are discussed.

Computational biorheology of human blood flow in health and disease

Hematologic disorders arising from infectious diseases, hereditary factors and environmental influences can lead to, and can be influenced by, significant changes in the shape, mechanical and physical properties of red blood cells (RBCs), and the biorheology of blood flow. Hence, modeling of hematologic disorders should take into account the multiphase nature of blood flow, especially in arterioles and capillaries. We present here an overview of a general computational framework based on dissipative particle dynamics (DPD) which has broad applicability in cell biophysics with implications for diagnostics, therapeutics and drug efficacy assessments for a wide variety of human diseases. This computational approach, validated by independent experimental results, is capable of modeling the biorheology of whole blood and its individual components during blood flow so as to investigate cell mechanistic processes in health and disease. DPD is a Lagrangian method that can be derived from systematic coarse-graining of molecular dynamics but can scale efficiently up to arterioles and can also be used to model RBCs down to the spectrin level. We start from experimental measurements of a single RBC to extract the relevant biophysical parameters, using single-cell measurements involving such methods as optical tweezers, atomic force microscopy and micropipette aspiration, and cell-population experiments involving microfluidic devices. We then use these validated RBC models to predict the biorheological behavior of whole blood in healthy or pathological states, and compare the simulations with experimental results involving apparent viscosity and other relevant parameters. While the approach discussed here is sufficiently general to address a broad spectrum of hematologic disorders including certain types of cancer, this paper specifically deals with results obtained using this computational framework for blood flow in malaria and sickle cell anemia.