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.

A Decoupled Finite Element Scheme for Simulating the Dynamics of Red Blood Cells in an L-shaped Cavity.

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

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

A Numerical Scheme for Coupled CFD and CSD Simulation of 3-D Fluid-Membrane Interactions

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

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

A fictitious domain method with a hybrid cell model for simulating motion of cells in fluid flow

Numerical Simulation of Hemodynamic Interactions of Red Blood Cells in Microcapillary Flow

An Energy Stable $C^0$ Finite Element Scheme for A Phase-Field Model of Vesicle Motion and Deformation

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

A Simulation of Blood Cells in Branching Capillaries

Unraveling the motion and deformation characteristics of red blood cells in a deterministic lateral displacement device

Numerical approximations of the Navier–Stokes equation coupled with volume-conserved multi-phase-field vesicles system: Fully-decoupled, linear, unconditionally energy stable and second-order time-accurate numerical scheme

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

A Fully Decoupled Numerical Method for Cahn–Hilliard–Navier–Stokes–Darcy Equations Based on Auxiliary Variable Approaches

Unresolved RBCs: An upscaling strategy for the CFD-DEM simulation of blood flow with deformable cells

Numerical Simulation of the Transient Shape of the Red Blood Cell in Microcapillary Flow

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

Numerical investigation on red blood cell flow based on unstructured grid.

Developing Hybrid cell-edge and cell-node Dissipative Compact Scheme for Complex Geometry Flows

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