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.

Augmenting the Immersed Boundary Method with Radial Basis Functions (RBFs) for the Modeling of Platelets in Hemodynamic Flows

A Radial Basis Function (RBF)-Finite Difference Method for the Simulation of Reaction-Diffusion Equations on Stationary Platelets within the Augmented Forcing Method

Immersed boundary simulations of cell-cell interactions in whole blood

An Improved Meshfree Scheme Based on Radial Basis Functions for Solving Incompressible Navier–Stokes Equations

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

Benchmarking the Immersed Boundary Method for Viscoelastic Flows

Hyperviscosity-Based Stabilization for Radial Basis Function-Finite Difference (RBF-FD) Discretizations of Advection-Diffusion Equations

Immersed boundary finite element method for blood flow simulation

A radial basis function partition of unity method for steady flow simulations

Numerical Simulation of Hydrodynamic Focusing of Particles in Straight Channel Flows with the Immersed Boundary-Lattice Boltzmann Method

A Stable Radial Basis Function Partition of Unity Method with d-Rectangular Patches for Modelling Water Flow in Porous Media

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

Lagrangian radial basis function‐based particle hydrodynamics method and its application for viscous flows

An RBF-FD polynomial method based on polyharmonic splines for the Navier-Stokes equations: Comparisons on different node layouts

A voxelized immersed boundary (VIB) finite element method for accurate and efficient blood flow simulation

Local divergence-free velocity interpolation for the immersed boundary method using composite B-splines

Space-time reduced basis methods for parametrized unsteady Stokes equations

Robust radial basis function interpolation based on geodesic distance for the numerical coupling of multiphysics problems

A novel Mortar Method Integration using Radial Basis Functions

A modified lattice Boltzmann approach based on radial basis function approximation for the non‐uniform rectangular mesh