1
|
Wu F, Lin J, Wang L, Lin S. Polymer Vesicles in a Nanochannel under Flow Fields: A DPD Simulation Study. MACROMOL THEOR SIMUL 2022. [DOI: 10.1002/mats.202200027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Fangsheng Wu
- Shanghai Key Laboratory of Advanced Polymeric Materials Key Laboratory for Ultrafine Materials of Ministry of Education Frontiers Science Center for Materiobiology and Dynamic Chemistry School of Materials Science and Engineering East China University of Science and Technology Shanghai 200237 China
| | - Jiaping Lin
- Shanghai Key Laboratory of Advanced Polymeric Materials Key Laboratory for Ultrafine Materials of Ministry of Education Frontiers Science Center for Materiobiology and Dynamic Chemistry School of Materials Science and Engineering East China University of Science and Technology Shanghai 200237 China
| | - Liquan Wang
- Shanghai Key Laboratory of Advanced Polymeric Materials Key Laboratory for Ultrafine Materials of Ministry of Education Frontiers Science Center for Materiobiology and Dynamic Chemistry School of Materials Science and Engineering East China University of Science and Technology Shanghai 200237 China
| | - Shaoliang Lin
- Shanghai Key Laboratory of Advanced Polymeric Materials Key Laboratory for Ultrafine Materials of Ministry of Education Frontiers Science Center for Materiobiology and Dynamic Chemistry School of Materials Science and Engineering East China University of Science and Technology Shanghai 200237 China
| |
Collapse
|
2
|
Numerical Simulations of Red-Blood Cells in Fluid Flow: A Discrete Multiphysics Study. CHEMENGINEERING 2021. [DOI: 10.3390/chemengineering5030033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
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.
Collapse
|
3
|
Livi C, Di Staso G, Clercx HJH, Toschi F. Influence of numerical resolution on the dynamics of finite-size particles with the lattice Boltzmann method. Phys Rev E 2021; 103:013303. [PMID: 33601495 DOI: 10.1103/physreve.103.013303] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Accepted: 12/17/2020] [Indexed: 11/07/2022]
Abstract
We investigate and compare the accuracy and efficiency of different numerical approaches to model the dynamics of finite-size particles using the lattice Boltzmann method (LBM). This includes the standard bounce-back (BB) and the equilibrium interpolation (EI) schemes. To accurately compare the different implementations, we first introduce a boundary condition to approximate the flow properties of an unbounded fluid in a finite simulation domain, taking into account the perturbation induced by a moving particle. We show that this boundary treatment is efficient in suppressing detrimental effects on the dynamics of spherical and ellipsoidal particles arising from the finite size of the simulation domain. We then investigate the performances of the BB and EI schemes in modeling the dynamics of a spherical particle settling under Stokes conditions, which can now be reproduced with great accuracy thanks to the treatment of the exterior boundary. We find that the EI scheme outperforms the BB scheme in providing a better accuracy scaling with respect to the resolution of the settling particle, while suppressing finite-size effects due to the particle discretization on the lattice grid. Additionally, in order to further increase the capability of the algorithm in modeling particles of sizes comparable to the lattice spacing, we propose an improvement to the EI scheme, the complete equilibrium interpolation (CEI). This approach allows us to accurately capture the boundaries of the particle also when located between two fluid nodes. We evaluate the CEI performance in solving the dynamics of an under-resolved particle under analogous Stokes conditions and also for the case of a rotating ellipsoid in a shear flow. Finally, we show that EI and CEI are able to recover the correct flow solutions also at small, but finite, Reynolds number. Adopting the CEI scheme it is not only possible to detect particles with zero lattice occupation, but also to increase up to one order of magnitude the accuracy of the dynamics of particles with a size comparable to the lattice spacing with respect to the BB and the EI schemes.
Collapse
Affiliation(s)
- C Livi
- Fluids and Flows Group and J.M. Burgers Centre for Fluid Dynamics, Department of Applied Physics, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - G Di Staso
- Fluids and Flows Group and J.M. Burgers Centre for Fluid Dynamics, Department of Applied Physics, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - H J H Clercx
- Fluids and Flows Group and J.M. Burgers Centre for Fluid Dynamics, Department of Applied Physics, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - F Toschi
- Fluids and Flows Group and J.M. Burgers Centre for Fluid Dynamics, Department of Applied Physics, Eindhoven University of Technology, Eindhoven, The Netherlands
| |
Collapse
|
4
|
Coclite A. Vascular journey and adhesion mechanics of micro-sized carriers in narrow capillaries. Microvasc Res 2020; 132:104069. [PMID: 32888941 DOI: 10.1016/j.mvr.2020.104069] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2020] [Revised: 08/26/2020] [Accepted: 08/28/2020] [Indexed: 10/23/2022]
Abstract
In this work a Lattice Boltzmann-Immersed Boundary method is used for predicting the dynamics of rigid and deformable adhesive micro-carriers (1 μm) navigating a capillary by the size of 10 μm with 20% hematocrit. Red cells and particles are modeled as a collection of mass-spring elements responding to a bending potential, an elastic potential and total enclosed area conservation constraint. Furthermore, particle surfaces are uniformly decorated with adhesive molecules (ligands) interacting with receptors disposed on the walls. Particle adhesion is modeled as a short-range ligad-receptor interaction and in term of formation and destruction probability functions that discriminate whether a chemical bond can be formed or destroyed. If a bond is established an attractive elastic force is activated. Particle transport and adhesion are characterized in terms of their ability to reach the capillary peripheries (margination rate) and firmly adhere the vasculature. This analysis is carried out systematically by varying particles' and cells' releasing positions and stiffness (Ca = 0 and 10-2). Moreover, three rigid and soft representative particles are transported on a finer mesh (Δx = 15 nm) and the chemical strength of their adhesive coating is varied (σ = 0.5, 1.0, and 2.0) to precisely analyze the resulting adhesion mechanics. Stiffness is found to weakly influence the margination rate while significantly affect the ability of such constructs to efficiently interact with the endothelium by forming stable chemical bonds.
Collapse
Affiliation(s)
- Alessandro Coclite
- Scuola di Ingegneria, Università degli Studi della Basilicata, Viale dell'Ateneo Lucano, 85100 Potenza, Italy.
| |
Collapse
|
5
|
Tiribocchi A, Bonaccorso F, Lauricella M, Melchionna S, Montessori A, Succi S. Curvature dynamics and long-range effects on fluid-fluid interfaces with colloids. SOFT MATTER 2019; 15:2848-2862. [PMID: 30816901 DOI: 10.1039/c8sm02396d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We investigate the dynamics of a phase-separating binary fluid, containing colloidal dumbbells anchored to the fluid-fluid interface. Extensive lattice Boltzmann-immersed boundary method simulations reveal that the presence of soft dumbbells can significantly affect the curvature dynamics of the interface between phase-separating fluids, even though the coarsening dynamics is left nearly unchanged. In addition, our results show that the curvature dynamics exhibits distinct non-local effects, which might be exploited for the design of new soft mesoscale materials. We point out that the inspection of the statistical dynamics of the curvature can disclose new insights into local inhomogeneities of the binary fluid configuration, as a function of the volume fraction and aspect ratio of the dumbbells.
Collapse
Affiliation(s)
- A Tiribocchi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy. and Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185, Rome, Italy.
| | - F Bonaccorso
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy.
| | - M Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185, Rome, Italy.
| | - S Melchionna
- ISC-CNR, Istituto Sistemi Complessi, Università Sapienza, P.le A. Moro 2, 00185 Rome, Italy.
| | - A Montessori
- Department of Engineering, University of Rome, "Roma Tre" Via Vito Volterra 62, 00146 Rome, Italy.
| | - S Succi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy. and Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185, Rome, Italy. and Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA.
| |
Collapse
|
6
|
Ye T, Phan-Thien N, Lim CT. Particle-based simulations of red blood cells—A review. J Biomech 2016; 49:2255-2266. [DOI: 10.1016/j.jbiomech.2015.11.050] [Citation(s) in RCA: 94] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2015] [Accepted: 11/07/2015] [Indexed: 11/29/2022]
|
7
|
A micro-scale simulation of red blood cell passage through symmetric and asymmetric bifurcated vessels. Sci Rep 2016; 6:20262. [PMID: 26830454 PMCID: PMC4735796 DOI: 10.1038/srep20262] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2015] [Accepted: 12/31/2015] [Indexed: 12/20/2022] Open
Abstract
Blood exhibits a heterogeneous nature of hematocrit, velocity, and effective viscosity in microcapillaries. Microvascular bifurcations have a significant influence on the distribution of the blood cells and blood flow behavior. This paper presents a simulation study performed on the two-dimensional motions and deformation of multiple red blood cells in microvessels with diverging and converging bifurcations. Fluid dynamics and membrane mechanics were incorporated. Effects of cell shape, hematocrit, and deformability of the cell membrane on rheological behavior of the red blood cells and the hemodynamics have been investigated. It was shown that the blood entering the daughter branch with a higher flow rate tended to receive disproportionally more cells. The results also demonstrate that red blood cells in microvessels experienced lateral migration in the parent channel and blunted velocity profiles in both straight section and daughter branches, and this effect was influenced by the shape and the initial position of the cells, the hematocrit, and the membrane deformability. In addition, a cell free region around the tip of the confluence was observed. The simulation results are qualitatively consistent with existing experimental findings. This study may provide fundamental knowledge for a better understanding of hemodynamic behavior of micro-scale blood flow.
Collapse
|
8
|
Bodnár T, Fasano A, Sequeira A. Mathematical Models for Blood Coagulation. FLUID-STRUCTURE INTERACTION AND BIOMEDICAL APPLICATIONS 2014. [DOI: 10.1007/978-3-0348-0822-4_7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
9
|
Fedosov DA, Dao M, Karniadakis GE, Suresh S. Computational biorheology of human blood flow in health and disease. Ann Biomed Eng 2013; 42:368-87. [PMID: 24419829 DOI: 10.1007/s10439-013-0922-3] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2013] [Accepted: 10/02/2013] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Dmitry A Fedosov
- Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425, Jülich, Germany
| | | | | | | |
Collapse
|
10
|
Melchionna S, Amati G, Bernaschi M, Bisson M, Succi S, Mitsouras D, Rybicki FJ. Risk assessment of atherosclerotic plaques based on global biomechanics. Med Eng Phys 2013; 35:1290-7; discussion 1290. [DOI: 10.1016/j.medengphy.2013.02.002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2012] [Revised: 01/31/2013] [Accepted: 02/05/2013] [Indexed: 01/14/2023]
|
11
|
Krüger T, Gross M, Raabe D, Varnik F. Crossover from tumbling to tank-treading-like motion in dense simulated suspensions of red blood cells. SOFT MATTER 2013; 9:9008-15. [PMID: 25353617 DOI: 10.1039/c3sm51645h] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
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.
Collapse
Affiliation(s)
- Timm Krüger
- Centre for Computational Science, University College London, 20 Gordon Street, London WC1H 0AJ, UK
| | | | | | | |
Collapse
|
12
|
Pontrelli G, Halliday I, Spencer TJ, König CS, Collins MW. Modelling the glycocalyx–endothelium–erythrocyte interaction in the microcirculation: a computational study. Comput Methods Biomech Biomed Engin 2013; 18:351-61. [DOI: 10.1080/10255842.2013.799146] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
13
|
Fedosov DA, Noguchi H, Gompper G. Multiscale modeling of blood flow: from single cells to blood rheology. Biomech Model Mechanobiol 2013; 13:239-58. [PMID: 23670555 DOI: 10.1007/s10237-013-0497-9] [Citation(s) in RCA: 120] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2012] [Accepted: 04/27/2013] [Indexed: 10/26/2022]
Abstract
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.
Collapse
Affiliation(s)
- Dmitry A Fedosov
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 , Jülich, Germany,
| | | | | |
Collapse
|
14
|
Ju M, Ye SS, Namgung B, Cho S, Low HT, Leo HL, Kim S. A review of numerical methods for red blood cell flow simulation. Comput Methods Biomech Biomed Engin 2013; 18:130-40. [PMID: 23582050 DOI: 10.1080/10255842.2013.783574] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
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.
Collapse
Affiliation(s)
- Meongkeun Ju
- a Department of Bioengineering , National University of Singapore , Singapore
| | | | | | | | | | | | | |
Collapse
|
15
|
Li, X, Vlahovska PM, Karniadakis GE. Continuum- and particle-based modeling of shapes and dynamics of red blood cells in health and disease. SOFT MATTER 2013; 9:28-37. [PMID: 23230450 PMCID: PMC3516861 DOI: 10.1039/c2sm26891d] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Xuejin Li,
- Division of Applied Mathematics, Brown University, Providence, RI 02912, USA
| | | | | |
Collapse
|
16
|
Cerdà JJ, Holm C, Kremer K. Novel Simulation Approaches for Polymeric and Soft Matter Systems. MACROMOL THEOR SIMUL 2011. [DOI: 10.1002/mats.201100072] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|