A multiscale red blood cell model with accurate mechanics, rheology, and dynamics

Prediction of mechanical hemolysis in medical devices via a Lagrangian strain-based multiscale model

This work introduces a new Lagrangian strain-based model to predict the shear-induced hemolysis in biomedical devices. Current computational models for device-induced hemolysis usually utilize empirical fitting of the released free hemoglobin (Hb) in plasma from damaged red blood cells (RBCs). These empirical correlations contain parameters that depend on specific device and operating conditions, thus cannot be used to predict hemolysis in a general device. The proposed algorithm does not have any empirical parameters, thus can presumably be used for hemolysis prediction in various blood-wetting medical devices. In contrast to empirical correlations in which the Hb release is related to the shear stress and exposure time without considering the physical processes, the proposed model links flow-induced deformation of the RBC membrane to membrane permeabilization and Hb release. In this approach, once the steady-state numerical solution of blood flow in the device is obtained under a prescribed operating condition, sample path lines are traced from the inlet of the device to the outlet to calculate the history of the shear stress tensor. In solving the fluid flow, it is assumed that RBCs do not have any influence on the flow pattern. Along each path line, shear stress tensor will be input into a coarse-grained (CG) RBC model to calculate the RBC deformation. Then the correlations obtained from molecular dynamics (MD) simulations are applied to relate the local areal RBC deformation to the perforated area on the RBC membrane. Finally, Hb released out of transient pores is calculated over each path line via a diffusion equation considering the effects of the steric hindrance and increased hydrodynamic drag due to the size of the Hb molecule. The total index of hemolysis (IH) is calculated by integration of released Hb over all the path lines in the computational domain. Hemolysis generated in the Food and Drug Administration (FDA) nozzle and two blood pumps, that is, a CentriMag blood pump (a centrifugal pump) and HeartMate II (an axial pump), for different flow regimes including the laminar and turbulent flows are calculated via the proposed algorithm. In all the simulations, the numerical predicted IH is close to the range of experimental data. The results promisingly indicate that this multiscale approach can be used as a tool for predicting hemolysis and optimizing the hematologic design of other types of blood-wetting devices.

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.

Multiplicity of stable orbits for deformable prolate capsules in shear flow

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

Recent advances in computational modeling of fibrin clot formation: A review

The field of thrombosis and hemostasis is crucial for understanding and developing new therapies for pathologies such as deep vein thrombosis, diabetes related strokes, pulmonary embolisms, and hemorrhaging related diseases. In the last two decades, an exponential growth in studies related to fibrin clot formation using computational tools has been observed. Despite this growth, the complete mechanism behind thrombus formation and hemostasis has been long and rife with obstacles; however, significant progress has been made in the present century. The computational models and methods used in this context are diversified into different spatiotemporal scales, yet there is no single model which can predict both physiological and mechanical properties of fibrin clots. In this review, we list the major strategies employed by researchers in modeling fibrin clot formation using recent and existing computational techniques. This review organizes the computational strategies into continuum level, system level, discrete particle (DPD), and multi-scale methods. We also discuss strengths and weaknesses of various methods and future directions in which computational modeling of fibrin clots can advance.

The key events of thrombus formation: platelet adhesion and aggregation

Thrombus formation is a complex, dynamic and multistep process, involving biochemical reactions, mechanical stimulation, hemodynamics, and so on. In this study, we concentrate on its two crucial steps: (i) platelets adhered to a vessel wall, or simply platelet adhesion, and (ii) platelets clumping and arrested to the adherent platelets, named platelet aggregation. We report the first direct simulation of three modes of platelet adhesion, detachment, rolling adhesion and firm adhesion, as well as the formation, disintegration, arrestment and consolidation of platelet plugs. The results show that the bond dissociation in the detachment mode is mainly attributed to a high probability of rupturing bonds, such that any existing bond can be quickly ruptured and all bonds would be completely broken. In the rolling adhesion, however, it is mainly attributed to the strong traction from the shear flow or erythrocytes, causing that the bonds are ruptured at the trailing edge of the platelet. The erythrocytes play an important role in platelet activities, such as the formation, disintegration, arrestment and consolidation of platelet plugs. They exert an aggregate force on platelets, a repulsion at a near distance but an attraction at a far distance to the platelets. This aggregate force can promote platelets to form a plug and/or bring along a part of a platelet plug causing its disintegration. It also greatly influences the arrestment and consolidation of platelet plugs, together with the adhesive force from the thrombus.

Simulation of circulating tumor cell transport and adhesion in cell suspensions in microfluidic devices

Understanding cell transport and adhesion dynamics under flow is important for many biotransport problems. We investigated the influence of cell size, ligand coating density, micropost size, and intercellular collisions on circulating tumor cell adhesion and transport in microfluidic devices. The cells were modeled as coarse-grained cell membranes and the adhesion was modeled as pairwise interacting potentials, while the fluid was solved using the lattice Boltzmann method. The coupling between the cell and the fluid was achieved through the immersed boundary method. The cell showed transient rolling adhesion in high shear regions and firm adhesion in low shear regions. The adhesive force for rolling cells on a micropost was increasing before the cell reached the crest of the post and then decreasing afterward. The adhesive strength for cells increases with ligand coating density. Cell trajectories in a microfluidic device with a shifted post design were studied as well. At low concentrations, the majority of the cells follow streamlines closely. However, the intercellular collision and collision from red blood cells impacted the cell trajectories. An L 2 norm of | e | was defined to characterize the difference between the cell trajectories and the associated streamlines. It was shown that | e | L 2 increases with micropost sizes and cell concentrations.

Importance of Erythrocyte Deformability for the Alignment of Malaria Parasite upon Invasion

Invasion of erythrocytes by merozoites is an essential step for the survival and progression of malaria parasites. To invade red blood cells (RBCs), apicomplexan parasites have to adhere with their apex to the RBC membrane. This necessary apex-membrane contact (or alignment) is not immediately established because the orientation of a free merozoite with respect to the RBC membrane is random when an adhesion contact first occurs. Therefore, it has been suggested that after the initial adhesion, merozoites facilitate their proper alignment by inducing considerable membrane deformations, frequently observed before the invasion process. This proposition is based on a positive correlation between RBC membrane deformation and successful invasion; however, the role of RBC mechanics and its deformation in the alignment process remains elusive. Using a mechanically realistic model of a deformable RBC, we investigate numerically the importance of RBC deformability for merozoite alignment. Adhesion between the parasite and RBC membrane is modeled by an attractive potential that might be inhomogeneous, mimicking possible adhesion gradients at the surface of a parasite. Our results show that RBC membrane deformations are crucial for successful merozoite alignment and require interaction strengths comparable to adhesion forces measured experimentally. Adhesion gradients along the parasite body further improve its alignment. Finally, an increased membrane rigidity is found to result in poor merozoite alignment, which can be a possible reason for a reduction in the invasion susceptibility of RBCs in several blood diseases associated with membrane stiffening.

Deformation of human red blood cells in extensional flow through a hyperbolic contraction

Flow-induced damage to red blood cells has been an issue of considerable importance since the introduction of the first cardiovascular devices. Early blood damage prediction models were based on measurements of damage by shear stress only. Subsequently, these models were extrapolated to include other components of the fluid stress tensor. However, the expanded models were not validated by measurements of damage in response to the added types of stress. Recent investigations have proposed that extensional stress might be more damaging to red cells than shear stress. In this study, experiments were conducted to compare human red cell deformation under laminar extensional stress versus laminar shear stress. It was found that the deformation caused by shear stress is matched by that produced by an extensional stress that is approximately 34 times smaller. Assuming that blood damage scales directly with cell deformation, this result indicates that mechanistic blood damage prediction models should weigh extensional stress more than shear stress.

Near-Newtonian Blood Behavior - Is It Good to Be a Camel?

From a certain level of exercise-intensity onward, hematocrit increases in horses, which brings more oxygen carriers into the bloodstream. Camels, however, when used in competitive racing could be even in need of iron supplementation and blood transfusions due to a severe reduction of their available hematocrit compared to their resting hematocrit. Since the extrinsic and intrinsic mechanical properties of camel erythrocytes (RBC) are so different compared to RBCs of other mammals, the question arises whether this observation might be a response to endurance exercise aiming at keeping the RBC count low. Rheometry indicated dromedary camel blood to behave almost Newtonian, which is unique amongst mammals. Shear thinning did increase with the hematocrit, but remained marginal compared to horses. As a result, camel whole blood viscosity (WBV) exceeded horse WBV at high shear rates, an effect, which was significantly augmented when the packed cell volume (PCV) was increased. Therefore, in camels any infusion of RBCs into the bloodstream can increase the cardiac work and the energy input into the endothelium more effectively, which should generate vascular remodeling in the long term. Yielding, however, was completely absent in camel blood, confirming low cohesion between its components at quasi-static flow. Camel blood remained a viscous liquid without a threshold even at unphysiologically high PCVs. This can help to washout lactate when camels start to dehydrate and might contribute to the sustained working ability of these animals. The subtle pseudoplastic behavior and the high viscosity contrast across the RBC membrane point to weak coupling between blood flow and red cell behavior. We predict that RBCs flow as separate entities and can show various types of motion, which can lead to friction instead of being collectively aligned to the flow direction. In comparison to horses, this behavior will become relevant at higher RBC counts in front of flow obstacles and possibly cause vascular remodeling if the PCV rises during strenuous exercise, a matter that should be avoided.

State diagram for wall adhesion of red blood cells in shear flow: from crawling to flipping

Red blood cells in shear flow show a variety of different shapes due to the complex interplay between hydrodynamics and membrane elasticity. Malaria-infected red blood cells become generally adhesive and less deformable. Adhesion to a substrate leads to a reduction in shape variability and to a flipping motion of the non-spherical shapes during the mid-stage of infection. Here, we present a complete state diagram for wall adhesion of red blood cells in shear flow obtained by simulations, using a particle-based mesoscale hydrodynamics approach, multiparticle collision dynamics. We find that cell flipping at a substrate is replaced by crawling beyond a critical shear rate, which increases with both membrane stiffness and viscosity contrast between the cytosol and suspending medium. This change in cell dynamics resembles the transition between tumbling and tank-treading for red blood cells in free shear flow. In the context of malaria infections, the flipping-crawling transition would strongly increase the adhesive interactions with the vascular endothelium, but might be suppressed by the combined effect of increased elasticity and viscosity contrast.

Deformation and dynamics of erythrocytes govern their traversal through microfluidic devices with a deterministic lateral displacement architecture

Deterministic lateral displacement (DLD) microfluidic devices promise versatile and precise processing of biological samples. However, this prospect has been realized so far only for rigid spherical particles and remains limited for biological cells due to the complexity of cell dynamics and deformation in microfluidic flow. We employ mesoscopic hydrodynamics simulations of red blood cells (RBCs) in DLD devices with circular posts to better understand the interplay between cell behavior in complex microfluidic flow and sorting capabilities of such devices. We construct a mode diagram of RBC behavior (e.g., displacement, zig-zagging, and intermediate modes) and identify several regimes of RBC dynamics (e.g., tumbling, tank-treading, and trilobe motion). Furthermore, we link the complex interaction dynamics of RBCs with the post to their effective cell size and discuss relevant physical mechanisms governing the dynamic cell states. In conclusion, sorting of RBCs in DLD devices based on their shear elasticity is, in general, possible but requires fine-tuning of flow conditions to targeted mechanical properties of the RBCs.

Ocular anterior chamber blood cell population differentiation using spectroscopic optical coherence tomography

There is potential clinical significance in identifying cellular responses in the anterior chamber (AC) of the eye, which can indicate hyphema (an accumulation of red blood cells [RBCs]) or aberrant intraocular inflammation (an accumulation of white blood cells [WBCs]). In this work, we developed a spectroscopic OCT analysis method to differentiate between populations of RBCs and subtypes of WBCs, including granulocytes, lymphocytes and monocytes, both in vitro and in ACs of porcine eyes. We developed an algorithm to track single cells within OCT data sets, and extracted the backscatter reflectance spectrum of each single cell from the detected interferograms using the short-time Fourier transform (STFT). A look-up table of Mie back-scattering spectra was generated and used to correlate the backscatter spectral features of single cells to their characteristic sizes. The extracted size distributions based on the best Mie spectra fit were significantly different between each cell type. We also studied theoretical backscattering models of single RBCs to further validate our experimental results. The described work is a promising step towards clinically differentiating and quantifying AC blood cell types.

Red blood cell distribution in a microvascular network with successive bifurcations

Nonproportional RBC distribution is an important characteristic in microvascular networks, which can result in heterogeneity of oxygen supply that may cause ischemic death in severe cases. In this paper, we perform three-dimensional numerical simulations of a large number of RBCs in a microvascular network, by using a hybrid method of smoothed dissipative particle dynamic and immersed boundary method. The distribution of multiple RBCs in a T-bifurcation is first simulated as a validation study, and a reasonable agreement is observed both qualitatively and quantitatively on the RBC flux between our results and the previously published numerical and empirical results. Next, the distribution of a large number of RBCs in a microvascular network is investigated, including the effects of cell deformability, aggregation and tube hematocrit. The simulation results indicate that decreased deformability and increased aggregation strength have a similar effect on the RBC distribution: the large RBC flux becomes larger, but the small becomes smaller. A high hematocrit also causes a similar phenomenon that the RBCs are more apt to flow into a high RBC-flux branch, because they are arranged compactly into a rouleaux and difficultly broken up at a high hematocrit. These results imply that lower cell deformability, stronger aggregation or higher tube hematocrit would be conducive to the phase separation of hematocrit and plasma skimming processes in microcirculation.

Lattice Boltzmann simulations of nonequilibrium fluctuations in a nonideal binary mixture

In recent years the lattice Boltzmann (LB) methodology has been fruitfully extended to include the effects of thermal fluctuations. So far, all studied cases pertain to equilibrium fluctuations, i.e., fluctuations with respect to an equilibrium background state. In this paper we take a step further and present results of fluctuating LB simulations of a binary mixture confined between two parallel walls in the presence of a constant concentration gradient in the wall-to-wall direction. This is a paradigmatic setup for the study of nonequilibrium (NE) fluctuations, i.e., fluctuations with respect to a nonequilibrium state. We analyze the dependence of the structure factors for the hydrodynamical fields on the wave vector q in both the directions parallel and perpendicular to the walls, highlighting the long-range (∼|q|^{-4}) nature of correlations in the NE framework. Results at the small scales (high wave numbers) quantitatively agree with the predictions of fluctuating hydrodynamics without fitting parameters. At larger scales (low wave numbers), however, results show finite-size effects induced by confinement and call for further studies aimed at controlling boundary conditions in the fluctuating LB framework as well as compressibility effects. Moreover, in the presence of a nonideal equation of state of the mixture, we also observe that the (spatially homogeneous) average pressure changes, due to a genuinely new contribution triggered by NE fluctuations. These NE pressure effects are studied at changing the system size and the concentration gradient. Taken all together, we argue that the results of this article are useful and instrumental to boost the applicability of the fluctuating LB methodology in the framework of NE fluctuations, possibly in conjunction with experiments.

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.

High-Throughput Microfluidic Characterization of Erythrocyte Shapes and Mechanical Variability

The motion of red blood cells (RBCs) in microchannels is important for microvascular blood flow and biomedical applications such as blood analysis in microfluidics. The current understanding of the complexity of RBC shapes and dynamics in microchannels is mainly based on several simulation studies, but there are a few systematic experimental investigations. Here, we present a combined study that systematically characterizes RBC behavior for a wide range of flow rates and channel sizes. Even though simulations and experiments generally show good agreement, experimental observations demonstrate that there is no single well-defined RBC state for fixed flow conditions but rather a broad distribution of states. This result can be attributed to the inherent variability in RBC mechanical properties, which is confirmed by a model that takes the variation in RBC shear elasticity into account. This represents a significant step toward a quantitative connection between RBC behavior in microfluidic devices and their mechanical properties, which is essential for a high-throughput characterization of diseased cells.

Numerical simulation of optical control for a soft particle in a microchannel

Technologies that use optical force to actively control particles in microchannels are a significant area of research interest in various fields. An optical force is generated by the momentum change caused by the refraction and reflection of light, which changes the particle surface as a function of the angle of incidence of light and which in turn feeds back and modifies the force on the particle. Simulating this phenomenon is a complex task. The deformation of a particle, the interaction between the surrounding fluid and the particle, and the reflection and refraction of light should be analyzed simultaneously. Herein, a deformable particle in a microchannel subjected to optical interactions is simulated using the three-dimensional lattice Boltzmann immersed-boundary method. The laser from the optical source is analyzed by dividing it into individual rays. To calculate the optical forces exerted on the particle, the intensity, momentum, and ray direction are calculated. The optical-separator problem with one optical source is analyzed by measuring the distance traveled because of the optical force. The optical-stretcher problem with two optical sources is then studied by analyzing the relation between the intensity of the optical source and particle deformation. This simulation will help the design of sorting and measuring by optical force.

Dynamic characteristics of a deformable capsule in a simple shear flow

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

Dynamics of deformable straight and curved prolate capsules in simple shear flow

This work investigates the motion of neutrally-buoyant, slightly deformable straight and curved prolate fluid-filled capsules in unbounded simple shear flow at zero Reynolds number using direct simulations. The curved capsules serve as a model for the typical crescent-shaped sickle red blood cells in sickle cell disease (SCD). The effects of deformability and curvature on the dynamics are revealed. We show that with low deformability, straight prolate spheroidal capsules exhibit tumbling in the shear plane as their unique asymptotically stable orbit. This result contrasts with that for rigid spheroids, where infinitely many neutrally stable Jeffery orbits exist. The dynamics of curved prolate capsules are more complicated due to a combined effect of deformability and curvature. At short times, depending on the initial orientation, slightly deformable curved prolate capsules exhibit either a Jeffery-like motion such as tumbling or kayaking, or a non-Jeffery-like behavior in which the director (end-to-end vector) of the capsule crosses the shear-gradient plane back and forth. At long times, however, a Jeffery-like quasiperiodic orbit is taken regardless of the initial orientation. We further show that the average of the long-time trajectory can be well approximated using the analytical solution for Jeffery orbits with an effective orbit constant C _eff and aspect ratio ℓ _eff. These parameters are useful for characterizing the dynamics of curved capsules as a function of given deformability and curvature. As the capsule becomes more deformable or curved, C _eff decreases, indicating a shift of the orbit towards log-rolling motion, while ℓ _eff increases weakly as the degree of curvature increases but shows negligible dependency on deformability. These features are not changed substantially as the viscosity ratio between the inner and outer fluids is changed from 1 to 5. As cell deformability, cell shape, and cell-cell interactions are all pathologically altered in blood disorders such as SCD, these results will have clear implications on improving our understanding of the pathophysiology of hematologic disease.

A non-linear fluid suspension model for blood flow

Motivated by the complex rheological behaviors observed in small/micro scale blood vessels, such as the Fahraeus effect, plasma-skimming, shear-thinning, etc., we develop a non-linear suspension model for blood. The viscosity is assumed to depend on the volume fraction (hematocrit) and the shear rate. The migration of the red blood cells (RBCs) is studied using a concentration flux equation. A parametric study with two representative problems, namely simple shear flow and a pressure driven flow demonstrate the ability of this reduced-order model to reproduce several key features of the two-fluid model (mixture theory approach), with much lower computational cost.

Effects of membrane reference state on shape memory of a red blood cell

By using a three-dimensional continuum model, we simulate the shape memory of a red blood cell after the remove of external forces. The purpose of this study is to illustrate the effect of membrane reference state on cell behavior during the recovery process. The reference state of an elastic element is the geometry with zero stress. Since the cell membrane is composed of cytoskeleton and lipid bilayer, both the reference states of cytoskeleton (RSC) and lipid bilayer (RSL) are considered. Results show that a non-spherical RSC can result in shape memory. The energy barrier due to non-spherical RSC is determined by the ratio of the equator length to the meridian length of the RSC. Thus different RSCs can have similar energy barrier and leading to identical recovery response. A series of simulations of more intermediate RSCs show that the recovery time scale is inversely proportional to the energy barrier. Comparing to spherical RSL, a spheroid RSL contributes to the energy barrier and recovery time. Furthermore, we observe a folding recovery due to the biconcave RSL which is different from the tank treading recovery. These results may motivate novel numerical and experimental studies to determine the exact RSC and RSL.

Membrane Wrapping Efficiency of Elastic Nanoparticles during Endocytosis: Size and Shape Matter

Using coarse-grained molecular dynamics simulations, we systematically investigate the receptor-mediated endocytosis of elastic nanoparticles (NPs) with different sizes, ranging from 25 to 100 nm, and shapes, including sphere-like, oblate-like, and prolate-like. Simulation results provide clear evidence that the membrane wrapping efficiency of NPs during endocytosis is a result of competition between receptor diffusion kinetics and thermodynamic driving force. The receptor diffusion kinetics refer to the kinetics of receptor recruitment that are affected by the contact edge length between the NP and membrane. The thermodynamic driving force represents the amount of required free energy to drive NPs into a cell. Under the volume constraint of elastic NPs, the soft spherical NPs are found to have similar contact edge lengths to rigid ones and to less efficiently be fully wrapped due to their elastic deformation. Moreover, the difference in wrapping efficiency between soft and rigid spherical NPs increases with their sizes, due to the increment of their elastic energy change. Furthermore, because of its prominent large contact edge length, the oblate ellipsoid is found to be the least sensitive geometry to the variation in NP's elasticity among the spherical, prolate, and oblate shapes during the membrane wrapping. In addition, simulation results indicate that conflicting experimental observations on the efficiency of cellular uptake of elastic NPs could be caused by their different mechanical properties. Our simulations provide a detailed mechanistic understanding about the influence of NPs' size, shape, and elasticity on their membrane wrapping efficiency, which serves as a rational guidance for the design of NP-based drug carriers.

Blood flow and transport in the human placenta

The placenta is a multi-functional organ that exchanges blood gases and nutrients between a mother and her developing fetus. In humans, fetal blood flows through intricate networks of vessels confined within villous trees, the branches of which are bathed in pools of maternal blood. Fluid mechanics and transport processes play a central role in understanding how these elaborate structures contribute to the function of the placenta, and how their disorganization may lead to disease. Recent advances in imaging and computation have spurred significant advances in simulations of fetal and maternal flows within the placenta, across a range of lengthscales. Models describe jets of maternal blood emerging from spiral arteries into a disordered and deformable porous medium, and solute uptake by fetal blood flowing through elaborate three-dimensional capillary networks. We survey recent developments and emerging challenges in modeling flow and transport in this complex organ.

Numerical design of a microfluidic chip for probing mechanical properties of cells

Microfluidic chips have been widely used to probe the mechanical properties of cells, which are recognized as a promising label-free biomarker for some diseases. In our previous work (Ye et al., 2018), we have studied the relationships between the transit time and the mechanical properties of a cell flowing through a microchannel with a single constriction, which potentially forms a basis for a microfluidic chip to measure cell's mechanical properties. Here, we investigate this microfluidic chip design and examine its potential in performances. We first develop the simultaneous dependence of the transit time on both the shear and bending moduli of a cell, and then examine the chip sensitivity with respect to the cell mechanical properties while serializing a single constriction along the flow direction. After that, we study the effect of the flow velocity on the transit time, and also test the chip's ability to identify heterogeneous cells with different mechanical properties. The results show that the microfluidic chip designed is capable of identifying heterogeneous cells, even when only one unhealthy cell is included. The serialization of chip can greatly increase the chip sensitivity with respect to the mechanical properties of cells. The flow with a higher velocity helps in not only promoting the chip throughput, but also in providing more accurate transit time measurements, because the cell prefers a symmetric deformation under a high velocity.

Quantifying Shear-Induced Deformation and Detachment of Individual Adherent Sickle Red Blood Cells

Vaso-occlusive crisis, a common painful complication of sickle cell disease, is a complex process triggered by intercellular adhesive interactions among blood cells and the endothelium in all human organs (e.g., the oxygen-rich lung as well as hypoxic systems such as liver and kidneys). We present a combined experimental-computational study to quantify the adhesive characteristics of sickle mature erythrocytes (SMEs) and irreversibly sickled cells (ISCs) under flow conditions mimicking those in postcapillary venules. We employed an in vitro microfluidic cell adherence assay, which is coated uniformly with fibronectin. We investigated the adhesion dynamics of SMEs and ISCs in pulsatile flow under well-controlled hypoxic conditions, inferring the cell adhesion strength by increasing the flow rate (or wall shear stress (WSS)) until the onset of cell detachment. In parallel, we performed simulations of individual SMEs and ISCs under shear. We introduced two metrics to quantify the adhesion process, the cell aspect ratio (AR) as a function of WSS and its rate of change (the dynamic deformability index). We found that the AR of SMEs decreases significantly with the increase of WSS, consistent between the experiments and simulations. In contrast, the AR of ISCs remains constant in time and independent of the flow rate. The critical WSS value for detaching a single SME in oxygenated state is in the range of 3.9-5.5 Pa depending on the number of adhesion sites; the critical WSS value for ISCs is lower than that of SMEs. Our simulations show that the critical WSS value for SMEs in deoxygenated state is above 6.2 Pa (multiple adhesion sites), which is greater than their oxygenated counterparts. We investigated the effect of cell shear modulus on the detachment process; we found that for the same cell adhesion spring constant, the higher shear modulus leads to an earlier cell detachment from the functionalized surface. These findings may aid in the understanding of individual roles of sickle cell types in sickle cell disease vaso-occlusion.

The sickle cell trait affects contact dynamics and endothelial cell activation in Plasmodium falciparum-infected erythrocytes

Sickle cell trait, a common hereditary blood disorder, protects carriers from severe disease in infections with the human malaria parasite Plasmodium falciparum. Protection is associated with a reduced capacity of parasitized erythrocytes to cytoadhere to the microvascular endothelium and cause vaso-occlusive events. However, the underpinning cellular and biomechanical processes are only partly understood and the impact on endothelial cell activation is unclear. Here, we show, by combining quantitative flow chamber experiments with multiscale computer simulations of deformable cells in hydrodynamic flow, that parasitized erythrocytes containing the sickle cell haemoglobin displayed altered adhesion dynamics, resulting in restricted contact footprints on the endothelium. Main determinants were cell shape, knob density and membrane bending. As a consequence, the extent of endothelial cell activation was decreased. Our findings provide a quantitative understanding of how the sickle cell trait affects the dynamic cytoadhesion behavior of parasitized erythrocytes and, in turn, endothelial cell activation.

Dissipative Coupling of Fluid and Immersed Objects for Modelling of Cells in Flow

Modelling of cell flow for biomedical applications relies in many cases on the correct description of fluid-structure interaction between the cell membrane and the surrounding fluid. We analyse the coupling of the lattice-Boltzmann method for the fluid and the spring network model for the cells. We investigate the bare friction parameter of fluid-structure interaction that is mediated via dissipative coupling. Such coupling mimics the no-slip boundary condition at the interface between the fluid and object. It is an alternative method to the immersed boundary method. Here, the fluid-structure coupling is provided by forces penalising local differences between velocities of the object's boundaries and the surrounding fluid. The method includes a phenomenological friction coefficient that determines the strength of the coupling. This work aims at determination of proper values of such friction coefficient. We derive an explicit formula for computation of this coefficient depending on the mesh density assuming a reference friction is known. We validate this formula on spherical and ellipsoidal objects. We also provide sensitivity analysis of the formula on all parameters entering the model. We conclude that such formula may be used also for objects with irregular shapes provided that the triangular mesh covering the object's surface is in some sense uniform. Our findings are justified by two computational experiments where we simulate motion of a red blood cell in a capillary and in a shear flow. Both experiments confirm our results presented in this work.

Shear rate dependent margination of sphere-like, oblate-like and prolate-like micro-particles within blood flow

This study investigates the shear rate dependent margination of micro-particles (MPs) with different shapes in blood flow through numerical simulations. We develop a multiscale computational model to handle the fluid-structure interactions involved in the blood flow simulations. The lattice Boltzmann method (LBM) is used to solve the plasma dynamics and a coarse-grained model is employed to capture the dynamics of red blood cells (RBCs) and MPs. These two solvers are coupled together by the immersed boundary method (IBM). The shear rate dependent margination of sphere MPs is firstly investigated. We find that margination of sphere MPs dramatically increases with the increment of wall shear rate [small gamma, Greek, dot above]ω under 800 s-1, induced by the breaking of rouleaux in blood flow. However, the margination probability only slowly grows when [small gamma, Greek, dot above]ω > 800 s-1. Furthermore, the shape effect of MPs is examined by comparing the margination behaviors of sphere-like, oblate-like and prolate-like MPs under different wall shear rates. We find that the margination of MPs is governed by the interplay of two factors: hydrodynamic collisions with RBCs including the collision frequency and collision displacement of MPs, and near wall dynamics. MPs that demonstrate poor performance in one process such as collision frequency may stand out in the other process like near wall dynamics. Specifically, the ellipsoidal MPs (oblate and prolate) with small aspect ratio (AR) outperform those with large AR regardless of the wall shear rate, due to their better performance in both the collision with RBCs and near wall dynamics. Additionally, we find there exists a transition shear rate region 700 s-1 < [small gamma, Greek, dot above]ω < 900 s-1 for all of these MPs: the margination probability dramatically increases with the shear rate below this region and slowly grows above this region, similar to sphere MPs. We further use the surface area to volume ratio (SVR) to distinguish different shaped MPs and illustrate their shear rate dependent margination in a contour in the shear rate-SVR plane. It is of significance that we can approximately predict the margination of MPs with a specific SVR. All these simulation results can be potentially applied to guide the design of micro-drug carriers for biomedical applications.

Flow-Induced Transitions of Red Blood Cell Shapes under Shear

A recent study of red blood cells (RBCs) in shear flow [Lanotte et al., Proc. Natl. Acad. Sci. U.S.A. 113, 13289 (2016)PNASA60027-842410.1073/pnas.1608074113] has demonstrated that RBCs first tumble, then roll, transit to a rolling and tumbling stomatocyte, and finally attain polylobed shapes with increasing shear rate, when the viscosity contrast between cytosol and blood plasma is large enough. Using two different simulation techniques, we construct a state diagram of RBC shapes and dynamics in shear flow as a function of shear rate and viscosity contrast, which is also supported by microfluidic experiments. Furthermore, we illustrate the importance of RBC shear elasticity for its dynamics in flow and show that two different kinds of membrane buckling trigger the transition between subsequent RBC states.

Numerical Investigation of the Effects of Red Blood Cell Cytoplasmic Viscosity Contrasts on Single Cell and Bulk Transport Behaviour

In-silico cellular models of blood are invaluable to gain understanding about the many interesting properties that blood exhibits. However, numerical investigations that focus on the effects of cytoplasmic viscosity in these models are not very prevalent. We present a parallelised method to implement cytoplasmic viscosity for HemoCell, an open-source cellular model based on immersed boundary lattice Boltzmann methods, using an efficient ray-casting algorithm. The effects of the implementation are investigated with single-cell simulations focusing on the deformation in shear flow, the migration due to wall induced lift forces, the characteristic response time in periodic stretching and pair collisions between red blood cells and platelets. Collective transport phenomena are also investigated in many-cell simulations in a pressure driven channel flow. The simulations indicate that the addition of a viscosity contrast between internal and external fluids significantly affects the deformability of a red blood cell, which is most pronounced during very short time-scale events. Therefore, modelling the cytoplasmic viscosity contrast is important in scenarios with high velocity deformation, typically high shear rate flows.

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.

Influence of particle size and shape on their margination and wall-adhesion: implications in drug delivery vehicle design across nano-to-micro scale

Intravascular drug delivery technologies majorly utilize spherical nanoparticles as carrier vehicles. Their targets are often at the blood vessel wall or in the tissue beyond the wall, such that vehicle localization towards the wall (margination) becomes a pre-requisite for their function. To this end, some studies have indicated that under flow environment, micro-particles have a higher propensity than nano-particles to marginate to the wall. Also, non-spherical particles theoretically have a higher area of surface-adhesive interactions than spherical particles. However, detailed systematic studies that integrate various particle size and shape parameters across nano-to-micro scale to explore their wall-localization behavior in RBC-rich blood flow, have not been reported. We address this gap by carrying out computational and experimental studies utilizing particles of four distinct shapes (spherical, oblate, prolate, rod) spanning nano- to-micro scale sizes. Computational studies were performed using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) package, with Dissipative Particle Dynamics (DPD). For experimental studies, model particles were made from neutrally buoyant fluorescent polystyrene spheres, that were thermo-stretched into non-spherical shapes and all particles were surface-coated with biotin. Using microfluidic setup, the biotin-coated particles were flowed over avidin-coated surfaces in absence versus presence of RBCs, and particle adhesion and retention at the surface was assessed by inverted fluorescence microscopy. Our computational and experimental studies provide a simultaneous analysis of different particle sizes and shapes for their retention in blood flow and indicate that in presence of RBCs, micro-scale non-spherical particles undergo enhanced 'margination + adhesion' compared to nano-scale spherical particles, resulting in their higher binding. These results provide important insight regarding improved design of vascularly targeted drug delivery systems.

Synergistic Integration of Laboratory and Numerical Approaches in Studies of the Biomechanics of Diseased Red Blood Cells

In red blood cell (RBC) disorders, such as sickle cell disease, hereditary spherocytosis, and diabetes, alterations to the size and shape of RBCs due to either mutations of RBC proteins or changes to the extracellular environment, lead to compromised cell deformability, impaired cell stability, and increased propensity to aggregate. Numerous laboratory approaches have been implemented to elucidate the pathogenesis of RBC disorders. Concurrently, computational RBC models have been developed to simulate the dynamics of RBCs under physiological and pathological conditions. In this work, we review recent laboratory and computational studies of disordered RBCs. Distinguished from previous reviews, we emphasize how experimental techniques and computational modeling can be synergically integrated to improve the understanding of the pathophysiology of hematological disorders.

Effect of spectrin network elasticity on the shapes of erythrocyte doublets

Red blood cell (RBC) aggregates play an important role in determining blood rheology. RBCs in plasma or polymer solution interact attractively to form various shapes of RBC doublets, where the attractive interactions can be varied by changing the solution conditions. A systematic numerical study on RBC doublet formation is performed, which takes into account the shear elasticity of the RBC membrane due to the spectrin cytoskeleton, in addition to the membrane bending rigidity. RBC membranes are modeled by two-dimensional triangular networks of linked vertices, which represent three-dimensional cell shapes. The phase space of RBC doublet shapes in a wide range of adhesion strengths, reduced volumes, and shear elasticities is obtained. The shear elasticity of the RBC membrane changes the doublet phases significantly. Experimental images of RBC doublets in different solutions show similar configurations. Furthermore, we show that rouleau formation is affected by the doublet structure.

Modeling thermal inkjet and cell printing process using modified pseudopotential and thermal lattice Boltzmann methods

Pseudopotential lattice Boltzmann methods (LBMs) can simulate a phase transition in high-density ratio multiphase flow systems. If coupled with thermal LBMs through equation of state, they can be used to study instantaneous phase transition phenomena with a high-temperature gradient where only one set of formulations in an LBM system can handle liquid, vapor, phase transition, and heat transport. However, at lower temperatures an unrealistic spurious current at the interface introduces instability and limits its application in real flow system. In this study, we proposed new modifications to the LBM system to minimize a spurious current which enables us to study nucleation dynamic at room temperature. To demonstrate the capabilities of this approach, the thermal ejection process is modeled as one example of a complex flow system. In an inkjet printer, a thermal pulse instantly heats up the liquid in a microfluidic chamber and nucleates bubble vapor providing the pressure pulse necessary to eject droplets at high speed. Our modified method can present a more realistic model of the explosive vaporization process since it can also capture a high-temperature/density gradient at nucleation region. Thermal inkjet technology has been successfully applied for printing cells, but cells are susceptible to mechanical damage or death as they squeeze out of the nozzle head. To study cell deformation, a spring network model, representing cells, is connected to the LBM through the immersed boundary method. Looking into strain and stress distribution of a cell membrane at its most deformed state, it is found that a high stretching rate effectively increases the rupture tension. In other words, membrane deformation energy is released through creation of multiple smaller nanopores rather than big pores. Overall, concurrently simulating multiphase flow, phase transition, heat transfer, and cell deformation in one unified LB platform, we are able to provide a better insight into the bubble dynamic and cell mechanical damage during the printing process.

Morphological and mechanical determinants of cellular uptake of deformable nanoparticles

Understanding the interactions of nanoparticles (NPs) with cell membranes and regulating their cellular uptake processes are of fundamental importance to the design of drug delivery systems with minimum toxicity, high efficiency and long circulation time. Employing the procedure of coarse-graining, we built an elastically deformable NP model with tunable morphological and mechanical properties. We found that the cellular uptake of deformable NPs depends on their shape: an increase in the particle elasticity significantly slows the uptake rate of spherical NPs, slightly retards that of prolate NPs, and promotes the uptake of oblate NPs. The intrinsic mechanisms have been carefully investigated through analysis of the endocytic mechanisms and free energy calculations. These findings provide unique insights into how deformable NPs penetrate across cell membranes and offer novel possibilities for designing effective NP-based carriers for drug delivery.

Understanding receptor-mediated endocytosis of elastic nanoparticles through coarse grained molecular dynamic simulation

For nanoparticle (NP)-based drug delivery platforms, the elasticity of the NPs has a significant influence on their blood circulation time and cellular uptake efficiency. However, due to the complexity of the endocytosis process and the inconsistency in the definition of elasticity for NPs in experiments, the understanding about the receptor-mediated endocytosis process of elastic NPs is still limited. In this work, we developed a coarse-grained molecular dynamics (CGMD) model for elastic NPs. The energy change of the elastic NPs can be precisely controlled by the bond, area, volume and bending potentials of this CGMD model. To represent liposomes with different elasticities, we systematically varied the bending rigidity of elastic NPs in CGMD simulations. Additionally, we changed the radius of the elastic NPs to explore the potential size effect. Through virtual nano-indentation tests, we found that the effective stiffness of elastic NPs was determined by their bending rigidity and size. Afterwards, we investigated the receptor-mediated endocytosis process of elastic NPs with different sizes and bending rigidities. We found that the membrane wrapping of soft NPs was faster than that of the stiff ones at the early stage, due to the NP deformation induced large contact area between the NPs and the membrane. However, because of the large energy penalties induced by the NP deformation, the membrane wrapping speed of soft NPs slows down during the late stage. Eventually, the soft NPs are wrapped less efficiently than the stiff ones during the membrane wrapping process. Through systematic CGMD simulations, we found a scaling law between the cellular uptake efficiency and the phenomenal bending rigidity of elastic NPs, which agrees reasonably well with experimental observations. Furthermore, we observed that the membrane wrapping efficiencies of soft and stiff NPs with large sizes were close to each other, due to the stronger ligand-receptor binding force and smaller difference in the stiffness of elastic NPs. Our computational model provides an effective tool to investigate the receptor-mediated endocytosis of elastic NPs with well controlled mechanical properties. This study can also be applied to guide the design of NP-based drug carriers with high efficacy, by utilizing their elastic properties.

Label-free sorting of soft microparticles using a bioinspired synthetic cilia array

Isolating cells of interest from a heterogeneous population has been of critical importance in biological studies and clinical applications. In this study, a novel approach is proposed for utilizing an active ciliary system in microfluidic devices to separate particles based on their physical properties. In this approach, the bottom of the microchannel is covered with an equally spaced cilia array of various patterns which is actuated by an external stimuli. 3D simulations are carried out to study cilia-particle interaction and isolation dynamic in a microfluidic channel. It is observed that these elastic hair-like filaments can influence particle's trajectories differently depending on their biophysical properties. This modeling study utilizes immersed boundary method coupled with the lattice Boltzmann method. Soft particles and cilia are implemented through the spring connected network model and point-particle scheme, respectively. It is shown that cilia array with proper stimulation is able to continuously and non-destructively separate cells into subpopulations based on their size, shape, and stiffness. At the end, a design map for fabrication of a programmable microfluidic device capable of isolating various subpopulations of cells is developed. This biocompatible, label-free design can separate cells/soft microparticles with high throughput which can greatly complement existing separation technologies.

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.

Adaptive resolution simulations coupling atomistic water to dissipative particle dynamics

Multiscale methods are the most efficient way to address the interlinked spatiotemporal scales encountered in soft matter and molecular liquids. In the literature reported hybrid approaches span from quantum to atomistic, coarse-grained, and continuum length scales. In this article, we present the hybrid coupling of the molecular dynamics (MD) and dissipative particle dynamics (DPD) methods, bridging the micro- and mesoscopic descriptions. The interfacing is performed within the adaptive resolution scheme (AdResS), which is a linear momentum conserving coupling technique. Our methodology is hence suitable to simulate fluids on the micro/mesoscopic scale, where hydrodynamics plays an important role. The presented approach is showcased for water at ambient conditions. The supramolecular coupling is enabled by a recently developed clustering algorithm SWINGER that assembles, disassembles, and reassembles clusters as needed during the course of the simulation. This allows for a seamless coupling between standard atomistic MD and DPD models. The developed framework can be readily applied to various applications in the fields of materials and life sciences, e.g., simulations of phospholipids and polymer melts, or to study the red blood cells behavior in normal and disease states.

Extended photoacoustic transport model for characterization of red blood cell morphology in microchannel flow

The dynamic response behavior of red blood cells holds the key to understanding red blood cell related diseases. In this regard, an understanding of the physiological functions of erythrocytes is significant before focusing on red blood cell aggregation in the microcirculatory system. In this work, we present a theoretical model for a photoacoustic signal that occurs when deformed red blood cells pass through a microfluidic channel. Using a Green's function approach, the photoacoustic pressure wave is obtained analytically by solving a combined Navier-Stokes and photoacoustic equation system. The photoacoustic wave expression includes determinant parameters for the cell deformability such as plasma viscosity, density, and red blood cell aggregation, as well as involving laser parameters such as beamwidth, pulse duration, and repetition rate. The effects of aggregation on blood rheology are also investigated. The results presented by this study show good agreements with the experimental ones in the literature. The comprehensive analytical solution of the extended photoacoustic transport model including a modified Morse type potential function sheds light on the dynamics of aggregate formation and demonstrates that the profile of a photoacoustic pressure wave has the potential for detecting and characterizing red blood cell aggregation.

Broken detailed balance and non-equilibrium dynamics in living systems: a review

Living systems operate far from thermodynamic equilibrium. Enzymatic activity can induce broken detailed balance at the molecular scale. This molecular scale breaking of detailed balance is crucial to achieve biological functions such as high-fidelity transcription and translation, sensing, adaptation, biochemical patterning, and force generation. While biological systems such as motor enzymes violate detailed balance at the molecular scale, it remains unclear how non-equilibrium dynamics manifests at the mesoscale in systems that are driven through the collective activity of many motors. Indeed, in several cellular systems the presence of non-equilibrium dynamics is not always evident at large scales. For example, in the cytoskeleton or in chromosomes one can observe stationary stochastic processes that appear at first glance thermally driven. This raises the question how non-equilibrium fluctuations can be discerned from thermal noise. We discuss approaches that have recently been developed to address this question, including methods based on measuring the extent to which the system violates the fluctuation-dissipation theorem. We also review applications of this approach to reconstituted cytoskeletal networks, the cytoplasm of living cells, and cell membranes. Furthermore, we discuss a more recent approach to detect actively driven dynamics, which is based on inferring broken detailed balance. This constitutes a non-invasive method that uses time-lapse microscopy data, and can be applied to a broad range of systems in cells and tissue. We discuss the ideas underlying this method and its application to several examples including flagella, primary cilia, and cytoskeletal networks. Finally, we briefly discuss recent developments in stochastic thermodynamics and non-equilibrium statistical mechanics, which offer new perspectives to understand the physics of living systems.

Cell Damage Index as Computational Indicator for Blood Cell Activation and Damage

Shear‐induced hemolysis is a major concern in the design and optimization of blood‐contacting devices. Even with a small amount of mechanical stress, inflammatory reactions can be triggered in the cells. Blood damage is typically estimated using continuum fluid dynamics simulations. In this study, we report a novel cell damage index (CDI) obtained by simulations on the single‐cell level in a lattice Boltzmann fluid flow. The change of the cell surface area gives important information on mechanical stress of individual cells as well as for whole blood. We are using predefined basic channel designs to analyze and compare the newly developed CDI to the conventional blood damage calculations in very weak shear stress scenarios. The CDI can incorporate both volume fraction and channel geometry information into a single quantitative value for the characterization of flow in artificial chambers.

Method for the simulation of blood platelet shape and its evolution during activation

We present a simple physically based quantitative model of blood platelet shape and its evolution during agonist-induced activation. The model is based on the consideration of two major cytoskeletal elements: the marginal band of microtubules and the submembrane cortex. Mathematically, we consider the problem of minimization of surface area constrained to confine the marginal band and a certain cellular volume. For resting platelets, the marginal band appears as a peripheral ring, allowing for the analytical solution of the minimization problem. Upon activation, the marginal band coils out of plane and forms 3D convoluted structure. We show that its shape is well approximated by an overcurved circle, a mathematical concept of closed curve with constant excessive curvature. Possible mechanisms leading to such marginal band coiling are discussed, resulting in simple parametric expression for the marginal band shape during platelet activation. The excessive curvature of marginal band is a convenient state variable which tracks the progress of activation. The cell surface is determined using numerical optimization. The shapes are strictly mathematically defined by only three parameters and show good agreement with literature data. They can be utilized in simulation of platelets interaction with different physical fields, e.g. for the description of hydrodynamic and mechanical properties of platelets, leading to better understanding of platelets margination and adhesion and thrombus formation in blood flow. It would also facilitate precise characterization of platelets in clinical diagnosis, where a novel optical model is needed for the correct solution of inverse light-scattering problem.

A parallel fluid-solid coupling model using LAMMPS and Palabos based on the immersed boundary method

The study of viscous fluid flow coupled with rigid or deformable solids has many applications in biological and engineering problems, e.g., blood cell transport, drug delivery, and particulate flow. We developed a partitioned approach to solve this coupled Multiphysics problem. The fluid motion was solved by Palabos (Parallel Lattice Boltzmann Solver), while the solid displacement and deformation was simulated by LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). The coupling was achieved through the immersed boundary method (IBM). The code modeled both rigid and deformable solids exposed to flow. The code was validated with the Jeffery orbits of an ellipsoid particle in shear flow, red blood cell stretching test, and effective blood viscosity flowing in tubes. It demonstrated essentially linear scaling from 512 to 8192 cores for both strong and weak scaling cases. The computing time for the coupling increased with the solid fraction. An example of the fluid-solid coupling was given for flexible filaments (drug carriers) transport in a flowing blood cell suspensions, highlighting the advantages and capabilities of the developed code.

Relationship between transit time and mechanical properties of a cell through a stenosed microchannel

The changes in the mechanical properties of a cell are not only the cause of some diseases, but can also be a biomarker for some disease states. In recent times, microfluidic devices with built-in constrictions have been widely used to measure these changes. The transit time in such devices, defined as the time that a cell takes to pass through a constriction, has been found to be a crucial factor associated with the cell mechanical properties. Here, we use smoothed dissipative particle dynamics (SDPD), a particle-based numerical method, to explore the relationship between the transit time and mechanical properties of a cell. Three expressions of the transit time are developed from our simulation data, with respect to the stenosed size of constrictions, the shear modulus and bending modulus of cells, respectively. We show that a convergent constriction (the inlet is wider than the outlet), and a sharp-corner constriction (the constriction outlet is narrow) are better in identifying the differences in the transit time of cells. Moreover, the transit time increases and gradually approaches a constant as the shear modulus of cells increases, but increases first and then decreases as the bending modulus increases. These results suggest that the mechanical properties of cells can indeed be measured by analyzing their transit time, based on the recommended microfluidic device.