Simplified particulate model for coarse-grained hemodynamics simulations

Hydrodynamic simulations of sedimenting dilute particle suspensions under repulsive DLVO interactions

We present guidelines to estimate the effect of electrostatic repulsion in sedimenting dilute particle suspensions. Our results are based on combined Langevin dynamics and lattice Boltzmann simulations for a range of particle radii, Debye lengths and particle concentrations. They show a simple relationship between the slope K of the concentration-dependent sedimentation velocity and the range χ of the electrostatic repulsion normalized by the average particle-particle distance. When χ → 0, the particles are too far away from each other to interact electrostatically and K = 6.55 as predicted by the theory of Batchelor. As χ increases, K likewise increases as if the particle radius increased in proportion to χ up to a maximum around χ = 0.4. Over the range χ = 0.4-1, K relaxes exponentially to a concentration-dependent constant consistent with known results for ordered particle distributions. Meanwhile the radial distribution function transitions from a disordered gas-like to a liquid-like form. Power law fits to the concentration-dependent sedimentation velocity similarly yield a simple master curve for the exponent as a function of χ, with a step-like transition from 1 to 1/3 centered around χ = 0.6.

Recent advances in blood rheology: a review

Due to the potential impact on the diagnosis and treatment of various cardiovascular diseases, work on the rheology of blood has significantly expanded in the last decade, both experimentally and theoretically. Experimentally, blood has been confirmed to demonstrate a variety of non-Newtonian rheological characteristics, including pseudoplasticity, viscoelasticity, and thixotropy. New rheological experiments and the development of more controlled experimental protocols on more extensive, broadly physiologically characterized, human blood samples demonstrate the sensitivity of aspects of hemorheology to several physiological factors. For example, at high shear rates the red blood cells elastically deform, imparting viscoelasticity, while at low shear rates, they form "rouleaux" structures that impart additional, thixotropic behavior. In addition to the advances in experimental methods and validated data sets, significant advances have also been made in both microscopic simulations and macroscopic, continuum, modeling, as well as novel, multiscale approaches. We outline and evaluate the most promising of these recent developments. Although we primarily focus on human blood rheology, we also discuss recent observations on variations observed across some animal species that provide some indication on evolutionary effects.

Numerical Simulations of Red-Blood Cells in Fluid Flow: A Discrete Multiphysics Study

In this paper, we present a methodological study of modelling red blood cells (RBCs) in shear-induced flows based on the discrete multiphysics (DMP) approach. The DMP is an alternative approach from traditional multiphysics based on meshless particle-based methods. The proposed technique has been successful in modelling multiphysics and multi-phase problems with large interfacial deformations such as those in biological systems. In this study, we present the proposed method and introduce an accurate geometrical representation of the RBC. The results were validated against available data in the literature. We further illustrate that the proposed method is capable of modelling the rupture of the RBC membrane with minimum computational difficulty.

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.

Computational study of radial particle migration and stresslet distributions in particle-laden turbulent pipe flow

Particle-laden turbulent flows occur in a variety of industrial applications as well as in naturally occurring flows. While the numerical simulation of such flows has seen significant advances in recent years, it still remains a challenging problem. Many studies investigated the rheology of dense suspensions in laminar flows as well as the dynamics of point-particles in turbulence. Here we employ a fully-resolved numerical simulation based on a lattice Boltzmann scheme, to investigate turbulent flow with large neutrally buoyant particles in a pipe flow at low Reynolds number and in dilute regimes. The energy input is kept fixed resulting in a Reynolds number based on the friction velocity around 250. Two different particle radii were used giving a particle-pipe diameter ratio of 0.05 and 0.075. The number of particles is kept constant resulting in a volume fraction of 0.54% and 1.83%, respectively. We investigated Eulerian and Lagrangian statistics along with the stresslet exerted by the fluid on the spherical particles. It was observed that the high particle-to-fluid slip velocity close to the wall corresponds locally to events of high energy dissipation, which are not present in the single-phase flow. The migration of particles from the inner to the outer region of the pipe, the dependence of the stresslet on the particle radial positions and a proxy for the fragmentation rate of the particles computed using the stresslet have been investigated.

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.

Inflow/Outflow Boundary Conditions for Particle-Based Blood Flow Simulations: Application to Arterial Bifurcations and Trees

When blood flows through a bifurcation, red blood cells (RBCs) travel into side branches at different hematocrit levels, and it is even possible that all RBCs enter into one branch only, leading to a complete separation of plasma and RBCs. To quantify this phenomenon via particle-based mesoscopic simulations, we developed a general framework for open boundary conditions in multiphase flows that is effective even for high hematocrit levels. The inflow at the inlet is duplicated from a fully developed flow generated in a pilot simulation with periodic boundary conditions. The outflow is controlled by adaptive forces to maintain the flow rate and velocity gradient at fixed values, while the particles leaving the arteriole at the outlet are removed from the system. Upon validation of this approach, we performed systematic 3D simulations to study plasma skimming in arterioles of diameters 20 to 32 microns. For a flow rate ratio 6:1 at the branches, we observed the “all-or-nothing” phenomenon with plasma only entering the low flow rate branch. We then simulated blood-plasma separation in arteriolar bifurcations with different bifurcation angles and same diameter of the daughter branches. Our simulations predict a significant increase in RBC flux through the main daughter branch as the bifurcation angle is increased. Finally, we demonstrated the effectiveness of the new methodology in simulations of blood flow in vessels with multiple inlets and outlets, constructed using an angiogenesis model.

Blood tests, which provide a wealth of information on the state of human health, are often performed on cell-free samples. Therefore, blood-plasma separation needs to be achieved. A simple but effective solution for isolating plasma from blood utilizes capillary bifurcations. In a particle-based simulation study of plasma skimming in capillary bifurcations, the blood flow properties such as velocity and pressure fields differ drastically at the inlet and outlet regions. Therefore, a new open (non-periodic) boundary is required. In this paper, we have developed and validated a general parallel framework for open boundary conditions. This is a non-trivial enabling technology that could be used in all open boundary systems and all particle-based Lagrangian simulations. We performed systematic 3D simulations of blood flow in arteriolar bifurcations and elucidated the biophysical mechanism of blood-plasma separation as well as quantified the effects of branch size and bifurcation angle on cell separation efficiency, which have not been addressed before. We also demonstrated the applicability of the methodology in arterial trees with multiple inlets and outlets.

Interface deformations affect the orientation transition of magnetic ellipsoidal particles adsorbed at fluid-fluid interfaces

Manufacturing new soft materials with specific optical, mechanical and magnetic properties is a significant challenge. Assembling and manipulating colloidal particles at fluid interfaces is a promising way to make such materials. We use lattice-Boltzmann simulations to investigate the response of magnetic ellipsoidal particles adsorbed at liquid-liquid interfaces to external magnetic fields. We provide further evidence for the first-order orientation phase transition predicted by Bresme and Faraudo [Journal of Physics: Condensed Matter, 2007, 19, 375110]. We show that capillary interface deformations around the ellipsoidal particle significantly affect the tilt-angle of the particle for a given dipole-field strength, altering the properties of the orientation transition. We propose scaling laws governing this transition, and suggest how to use these deformations to facilitate particle assembly at fluid-fluid interfaces.

Timescales of emulsion formation caused by anisotropic particles

Particle stabilized emulsions have received much interest in the recent past, but our understanding of the dynamics of emulsion formation is still limited. For simple spherical particles, the time dependent growth of fluid domains is dominated by the formation of droplets, particle adsorption and coalescence of droplets (Ostwald ripening), which eventually can be almost fully blocked due to the presence of the particles. Ellipsoidal particles are known to be more efficient stabilizers of fluid interfaces than spherical particles and their anisotropic shape and the related additional rotational degrees of freedom have an impact on the dynamics of emulsion formation. In this paper, we investigate this point by means of simple model systems consisting of a single ellipsoidal particle or a particle ensemble at a flat interface as well as a particle ensemble at a spherical interface. By applying combined multicomponent lattice Boltzmann and molecular dynamics simulations we demonstrate that the anisotropic shape of ellipsoidal particles causes two additional timescales to be of relevance in the dynamics of emulsion formation: a relatively short timescale can be attributed to the adsorption of single particles and the involved rotation of particles towards the interface. As soon as the interface is jammed, however, capillary interactions between the particles cause a local reordering on very long timescales leading to a continuous change in the interface configuration and increase of the interfacial area. This effect can be utilized to counteract the thermodynamic instability of particle stabilized emulsions and thus offers the possibility to produce emulsions with exceptional stability.

Simulation of finite-size fibers in turbulent channel flows

The dynamical behavior of almost neutrally buoyant finite-size rigid fibers or rods in turbulent channel flow is studied by direct numerical simulations. The time evolution of the fiber orientation and translational and rotational motions in a statistically steady channel flow is obtained for three different fiber lengths. The turbulent flow is modeled by an entropy lattice Boltzmann method, and the interaction between fibers and carrier fluid is modeled through an external boundary force method. Direct contact and lubrication force models for fiber-fiber interactions and fiber-wall interaction are taken into account to allow for a full four-way interaction. The density ratio is chosen to mimic cellulose fibers in water. It is shown that the finite size leads to fiber-turbulence interactions that are significantly different from earlier reported results for pointlike particles (e.g., elongated ellipsoids smaller than the Kolmogorov scale). An effect that becomes increasingly accentuated with fiber length is an accumulation in high-speed regions near the wall, resulting in a mean fiber velocity that is higher than the mean fluid velocity. The simulation results indicate that the finite-size fibers tend to stay in the high-speed streaks due to collisions with the wall. In the central region of the channel, long fibers tend to align in the spanwise direction. Closer to the wall the long fibers instead tend to toward to a rotation in the shear plane, while very close to the wall they become predominantly aligned in the streamwise direction.

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.

Crossover from tumbling to tank-treading-like motion in dense simulated suspensions of red blood cells

Via computer simulations, we provide evidence that the shear rate induced red blood cell tumbling-to-tank-treading transition also occurs at quite high volume fractions, where collective effects are important. The transition takes place as the ratio of effective suspension stress to the characteristic cell membrane stress exceeds a certain value and does not explicitly depend on volume fraction or cell deformability. This value coincides with that for a transition from an orientationally less ordered to a highly ordered phase. The average cell deformation does not show any signature of transition, but rather follows a simple scaling law independent of volume fraction.

Multiscale modeling of blood flow: from single cells to blood rheology

Mesoscale simulations of blood flow, where the red blood cells are described as deformable closed shells with a membrane characterized by bending rigidity and stretching elasticity, have made much progress in recent years to predict the flow behavior of blood cells and other components in various flows. To numerically investigate blood flow and blood-related processes in complex geometries, a highly efficient simulation technique for the plasma and solutes is essential. In this review, we focus on the behavior of single and several cells in shear and microcapillary flows, the shear-thinning behavior of blood and its relation to the blood cell structure and interactions, margination of white blood cells and platelets, and modeling hematologic diseases and disorders. Comparisons of the simulation predictions with existing experimental results are made whenever possible, and generally very satisfactory agreement is obtained.

A review of numerical methods for red blood cell flow simulation

In this review, we provide an overview of the simulation techniques employed for modelling the flow of red blood cells (RBCs) in blood plasma. The scope of this review omits the fluid modelling aspect while focusing on other key components in the RBC-plasma model such as (1) describing the RBC deformation with shell-based and spring-based RBC models, (2) constitutive models for RBC aggregation based on bridging theory and depletion theory and (3) additional strategies required for completing the RBC-plasma flow model. These include topics such as modelling fluid-structure interaction with the immersed boundary method and boundary integral method, and updating the variations in multiphase fluid property through the employment of index field methods. Lastly, we summarily discuss the current state and aims of RBC modelling and suggest some research directions for the further development of this field of modelling.

Continuum- and particle-based modeling of shapes and dynamics of red blood cells in health and disease

We review recent advances in multiscale modeling of the mechanics of healthy and diseased red blood cells (RBCs), and blood flow in the microcirculation. We cover the traditional continuum-based methods but also particle-based methods used to model both the RBCs and the blood plasma. We highlight examples of successful simulations of blood flow including malaria and sickle cell anemia.

Rotational behaviour of red blood cells in suspension: a mesoscale simulation study

The nature of blood as a suspension of red blood cells makes computational haemodynamics a demanding task. Our coarse-grained blood model, which builds on a lattice Boltzmann method for soft particle suspensions, enables the study of the collective behaviour of the order of 10(6) cells in suspension. After demonstrating the viscosity measurement in Kolmogorov flow, we focus on the statistical analysis of the cell orientation and rotation in Couette flow. We quantify the average inclination with respect to the flow and the nematic order as a function of shear rate and haematocrit. We further record the distribution of rotation periods around the vorticity direction and find a pronounced peak in the vicinity of the theoretical value for free model cells, even though cell-cell interactions manifest themselves in a substantial width of the distribution.