Abstract:We present a computational framework for the simulation of blood flow with fully resolved red blood cells (RBCs) using a modular approach that consists of a lattice Boltzmann solver for the blood plasma, a novel finite element based solver for the deformable bodies and an immersed boundary method for the fluid-solid interaction. For the RBCs, we propose a nodal projective FEM (npFEM) solver which has theoretical advantages over the more commonly used mass-spring systems (mesoscopic modeling), such as an unconditional stability, versatile material expressivity, and one set of parameters to fully describe the behavior of the body at any mesh resolution. At the same time, the method is substantially faster than other FEM solvers proposed in this field, and has an efficiency that is comparable to the one of mesoscopic models. At its core, the solver uses specially defined potential energies, and builds upon them a fast iterative procedure based on quasi-Newton techniques. For a known material, our solver has only one free parameter that demands tuning, related to the body viscoelasticity. In contrast, state-of-the-art solvers for deformable bodies have more free parameters, and the calibration of the models demands special assumptions regarding the mesh topology, which restrict their generality and mesh independence. We propose as well a modification to the potential energy proposed by Skalak et al. 1973 for the red blood cell membrane, which enhances the strain hardening behavior at higher deformations. Our viscoelastic model for the red blood cell, while simple enough and applicable to any kind of solver as a post-convergence step, can capture accurately the characteristic recovery time and tank-treading frequencies. The framework is validated using experimental data, and it proves to be scalable for multiple deformable bodies.

Coarse-grained Area-Difference-elasticity Membrane Model Coupled with IB–LB Method for Simulation of Red Blood Cell Morphology

An Efficient Red Blood Cell Model in the Frame of Ib-Lbm and Its Application

A Lattice Boltzmann Fictitious Domain Method For Modeling Red Blood Cell Deformation And Multiple-Cell Hydrodynamic Interactions In Flow

Modeling of Light Scattering by Biconcave and Deformed Red Blood Cells with the Invariant Imbedding T-matrix Method

Effects of Membrane Reference State on Shape Memory of a Red Blood Cell

A Multiscale Framework for Large Deformation Modeling of RBC Membranes

Numerical Simulation of Cell Migration and Deformation with Flow Using a Coupled Lattice Boltzmann and Phase-Field Method

Coarse-grained Simulation of Red Blood Cells Stretching Deformation

A Lattice Boltzmann Method for Simulating the Separation of Red Blood Cells at Microvascular Bifurcations

Three-dimensional Numerical Simulation of the Deformation and the Aggregation of Human Red Blood Cells

General coarse-grained red blood cell models: I. Mechanics

Shape changes and motion of a vesicle in a fluid using a lattice Boltzmann model

Deformation and Motion of A Red Blood Cell in A Shear Flow Simulated by A Lattice Boltzmann Model

Modeling Biological Membrane and Red Blood Cells by Coarse- Grained Particle Method

Continuum- and Particle-Based Modeling of Human Red Blood Cells

Coarse-Grained Molecular Dynamics Simulation of A Red Blood Cell

Using 3D Fluid–structure Interaction Model to Analyse the Biomechanical Properties of Erythrocyte

Numerical Computation Of The Elastic And Mechanical Properties Of Red Blood Cell Membrane Using The Higher-Order Cauchy-Born Rule

Numerical modeling of biomechanical responses of healthy red blood cell membrane under various loading conditions using a 3D quasicontinuum approach

Bridging the computational gap between mesoscopic and continuum modeling of red blood cells for fully resolved blood flow

Numerical Simulation of the Red Blood Cell Aggregation and Deformation Behaviors in Ultrasonic Field